Blog

Mycelial nutrient transfer promotes bacterial co-metabolic organochlorine pesticide degradation in nutrient-deprived environments | The ISME Journal

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Carousel with three slides shown at a time. Use the Previous and Next buttons to navigate three slides at a time, or the slide dot buttons at the end to jump three slides at a time. Organic Garden Pesticides

Mycelial nutrient transfer promotes bacterial co-metabolic organochlorine pesticide degradation in nutrient-deprived environments | The ISME Journal

Edivaldo Wilson de Lima, Bruno Pinheiro Brunaldi, … Paulo Renato Matos Lopes

François Thomas, Erwan Corre & Aurélie Cebron

David Rex, Timothy J. Clough, … Gary J. Lanigan

Peng Li, Shuifeng Ye, … Yunfei Wu

Kirsten Lønne Enggrob, Thomas Larsen, ... Jim Rasmussen

Chuan Chen, Lingyan Li, … Fang-Jie Zhao

Marta Gil-Martínez, Carmen M. Navarro-Fernández, … Teodoro Marañón

Johannes Friedl, Clemens Scheer, … Katharina M. Keiblinger

Zhihuang Xie, Zhenhua Yu, … Jian Jin

The ISME Journal (2023 )Cite this article

Biotransformation of soil organochlorine pesticides (OCP) is often impeded by a lack of nutrients relevant for bacterial growth and/or co-metabolic OCP biotransformation. By providing space-filling mycelia, fungi promote contaminant biodegradation by facilitating bacterial dispersal and the mobilization and release of nutrients in the mycosphere. We here tested whether mycelial nutrient transfer from nutrient-rich to nutrient-deprived areas facilitates bacterial OCP degradation in a nutrient-deficient habitat. The legacy pesticide hexachlorocyclohexane (HCH), a non-HCH-degrading fungus (Fusarium equiseti K3), and a co-metabolically HCH-degrading bacterium (Sphingobium sp. S8) isolated from the same HCH-contaminated soil were used in spatially structured model ecosystems. Using 13C-labeled fungal biomass and protein-based stable isotope probing (protein-SIP), we traced the incorporation of 13C fungal metabolites into bacterial proteins while simultaneously determining the biotransformation of the HCH isomers. The relative isotope abundance (RIA, 7.1–14.2%), labeling ratio (LR, 0.13–0.35), and the shape of isotopic mass distribution profiles of bacterial peptides indicated the transfer of 13C-labeled fungal metabolites into bacterial proteins. Distinct 13C incorporation into the haloalkane dehalogenase (linB) and 2,5-dichloro-2,5-cyclohexadiene-1,4-diol dehydrogenase (LinC), as key enzymes in metabolic HCH degradation, underpin the role of mycelial nutrient transport and fungal-bacterial interactions for co-metabolic bacterial HCH degradation in heterogeneous habitats. Nutrient uptake from mycelia increased HCH removal by twofold as compared to bacterial monocultures. Fungal-bacterial interactions hence may play an important role in the co-metabolic biotransformation of OCP or recalcitrant micropollutants (MPs).

Despite its ban or restricted use in many parts of the world, the organochloride pesticide hexachlorocyclohexane (HCH) continues to pose a serious environmental risk [1, 2] due to its toxicity and environmental persistence [1, 3,4,5]. Commercial production of HCH typically results in a mixture of four major isomers (α-HCH (60–70%), β-HCH (5–12%), γ-HCH (10–12%), δ-HCH (6–10%)) [3, 6], of which only γ-HCH has insecticidal activity while the other isomers are discarded as HCH-muck, thereby generating stockpiles of persistent toxic waste [3, 6,7,8]. Microorganisms including fungi [5, 9, 10] and bacteria [7, 11, 12] are able to degrade HCH isomers. However, an elaborate degradation pathway and the respective catabolic enzymes for complete mineralization of HCH have only been reported for bacterial members of the family Sphingomonadaceae under ideal laboratory conditions [3, 7, 11, 12]. Microbial degradation in heterogeneous environments though emerges from manifold biotic and abiotic factors including redox conditions, moisture, and – in cases of co-metabolic degradation – the presence and spatial availability of auxiliary organic growth substrates to potential degraders [13,14,15]. Contrary to bacteria, mycelial fungi are highly adapted to heterogeneous habitats due to their space-filling mycelial structures, and their ability to translocate resources within their myceliafrom nutrient-rich to nutrient-deficient patches [16, 17].

Fungi and bacteria thereby often share the same microhabitat and contribute to nutrient cycling and other biogeochemical processes [17,18,19], such as the biotransformation of organic contaminants [14, 20, 21]. In the mycosphere (i.e. the microhabitat that surrounds fungal hyphae and mycelia), bacteria and fungi may compete for the same substrates and, depending on the organisms and habitat conditions, their interactions may be either ecologically neutral, competitive, or cooperative [22, 23] and range from apparently random physical to specific commensal or symbiotic associations [19]. For instance, mycelia have been shown to actively transport polycyclic aromatic hydrocarbons (PAHs) thus improving their bioavailability to bacterial degradation in unsaturated environments [24,25,26]. They may further serve as transport vectors for bacteria and phages [27], thus improving the access of bacteria to otherwise poorly available nutrients and contaminants [20, 28]. Fungi are also known to shape bacterial community structures surrounding the hyphosphere by secreting carbonaceous compounds that can be utilized by fungus-associated bacteria [17, 25, 29, 30]. They thereby exude carbon substrates such as phenols, citrate, or oxalate that may serve as auxiliary carbon sources for pollutant-degrading bacteria [21].

Co-metabolic degradation offers the benefit of contaminant removal even at trace quantities since the degrader is not reliant on the contaminant for carbon or energy source [15, 31]. However, the role of fungal resource transfer to bacteria for improved bacterial activity and co-metabolic contaminant degradation in contaminated oligotrophic microhabitats remains vastly underexplored. In microbial ecology, stable isotope probing of proteins (protein-SIP) has become a valuable technique that allows the identification of key metabolic actors and functional interactions in microbial communities [32,33,34]. Protein-SIP combines meta-proteomics based on mass spectrometry and stable isotope probing techniques to sensitively detect the assimilation of labeled target substrate (such as 13C) at peptide level [32, 33]. The information gleaned from mass spectrometry data includes relative isotope abundance (RIA), labeling ratio (LR), and the shape of the isotopic distribution [33]. The RIA describes the percentage of the heavy isotope (i.e., 13C) atoms in a peptide, which provides information on the incorporation of the labeled substrate into biomass. The LR describes the ratio of labeled to unlabeled peptides which relates to protein turnover following the addition of the labeled substrate [33, 35].

In the present study, we used spatially fragmented laboratory model ecosystems to test the hypothesis whether mycelial nutrient transfer from eutrophic habitats promotes co-metabolic degradation of HCH by bacteria in oligotrophic contaminated environments. A non-HCH-degrading soil fungus Fusarium equiseti K3 and a co-metabolically HCH degrading bacterium Sphingobium sp. S8 isolated from the same HCH contaminated soil were used. The strain S8 was placed onto a nutrient-free HCH-treated agar patch that was separated from the fungal inoculum by an air gap. 13C-labeled mycelia of strain K3 were allowed to overgrow the bacterial region. Using a combination of metaproteomics, stable isotope labeling, and chemical analyses [32], we traced the incorporation of 13C of fungal origin into bacterial proteins responsible for both housekeeping and co-metabolic HCH degradation, respectively and simultaneously followed by chemical analytics. Our results demonstrate that mycelial fungi facilitate bacterial degradation of contaminants in nutrient-deficient environments by providing C-substrates to bacteria, thus enhancing their biotransformation capacity.

Hexane (LiChrosolv grade), Acetone (99.8%), and γ-HCH (97%) were either purchased from Sigma (Munich, Germany) or Merck (Darmstadt, Germany) and were of analytical grade. 13C-glucose (D-glucose13C6, 99%) was obtained from Cambridge Isotope Laboratories, Inc, (USA) The analytical standard containing a mixture of the HCH isomers, α: β: γ: δ = 1:1:1:1 was obtained from Merck (Darmstadt, Germany).

The ascomycete Fusarium equiseti K3 was isolated from HCH-contaminated soil collected from a former obsolete pesticide store in Kitengela, Kenya (GPS: 01.49 S, 37.048E) and its draft genome sequence was described earlier [36]. The fungus was not found to degrade HCH (see below). Fungal pre-cultures were cultivated at 25 °C on potato dextrose agar (PDA; BD Difco, USA) prior to the transfer to 13C glucose/glucose-containing minimum salt medium (MSM) agar. The HCH degrading bacterium Sphingobium sp. S8 was isolated from the same soil as the fungus and its draft genome sequence was described by Khan et al [37]. The bacterium was cultivated on Luria-Bertani (LB) (DSMZ medium 381) agar medium treated with γ-HCH (97%; Sigma) to a final concentration of 100 µg mL−1. To prepare the inoculum of the S8 strain, a 72-h old culture (OD600 ~1.5) of strain S8 in LB medium was centrifuged at 12,000 × g at 15 °C for 10 min, the cell pellet washed twice with MSM medium, and resuspended in the same medium to a concentration of 7.2 × 108 cells mL−1.

The absence of mutual fungal-bacterial inhibition was assessed using a standard procedure as described before [38]. Briefly, a fungal mycelium was inoculated in the middle of a PDA plate, and four spots of bacterial suspension (10 µL) inoculated onto the agar surface around the fungal inoculum at a distance of 1 cm. The plates were later incubated at 25 °C and mutual growth or growth inhibition was analyzed by visual inspection and microscopy.

Laboratory microcosms consisted of spatially separated nutrient-rich (NR) and nutrient-deficient (ND) zones that were separated by 2 mm air gaps (Fig. 1A). To create NR zones, two patches (Ø 1 cm) were cut from an MSM [39] agar plate supplemented with 10% w/v 13C glucose/glucose as sole carbon source as described below. The two patches (patches 1 & 3; Fig. 1A) were placed in a Petri dish (Ø 8.5 cm) at a distance of 1.4 cm apart. To represent a ND zone an MSM agar patch (l x w x h: 1 × 1 × 0.5 cm; patch 2) was placed between them, leaving a 2 mm air gap on either side (Fig. 1A). Patch 2 served as a spot for HCH spiking (in the case of HCH degradation setups) and bacterial inoculation. A circular disc (Ø 0.5 cm) punched from a freshly grown PDA plate of the fungal strain K3 was inoculated on patch 1 and allowed to grow until its hyphae touched the edges of patch 2 (after ca. 3 d). Fifty µL (i.e., ≈3.6 × 107 cells) HCH degrading strain S8 in MSM medium was applied to patch 2 of the microcosm. The microcosms were then incubated at 25 °C in the dark until the mycelia had overgrown the patches (ca. 10 d) and subsequently subjected to protein and HCH extraction and analysis.

A Scheme of the microcosm consisting of two nutrient rich (NR) agar patches (patch1 and patch 3) and a nutrient-deficient agar patch 2 (ND zone) separated by air gaps to mimic unsaturated soil habitats. The fungus was inoculated to patch 1 and allowed to overgrow bacteria on patch 2 in presence of HCH. B Bird’s view of the setup placed in a Petri dish and micrographs of Sphingobium sp. S8 growing as biofilm in close proximity to growing hyphae of F. equiseti K3 on patch 2 viewed under the light microscope AZ100M (Nikon, Düsseldorf, Germany). C Glucose content measured on patch 2 after 10 days in the presence of F. equiseti K3 and glass fiber. Different lowercase letters indicate statistically significant difference between the treatment groups (p < 0.05, n = 3), while the whiskers represent data variability outside the upper and the lower quantiles.

Patch 1 (Fig. 1A) was spiked with 13C glucose by liquefying MSM medium agar at 95 °C in a water bath and subsequently transferring a five mL aliquot into a sterile tube containing 0.5 g of 13C glucose. The tube was then warmed at constant shaking at 750 rpm in a thermomixer (Eppendorf, Germany) at 95 °C until complete glucose dissolution. Two droplets of 200 µL each of the 13C-glucose-agar (i.e., 40 mg glucose) were then placed into the Petri dishes (Ø 8.5 cm) at a distance of 1.4 cm apart and allowed to solidify under a sterile laminar flow hood to form patches (1 & 3) (Fig. 1A) while patch 2 was used as an inoculation spot for the HCH degrading bacteria. Parallel control setups including glass fiber experiments to mimic mycelial structures and control experiments with bacteria directly growing on a 13C-glucose enriched (~20 mg) agar patch were prepared (Fig. S1.1). Following a 10-d incubation period, the microcosms were subjected to protein and glucose extraction from patch 2 and analyzed as described below. All experiments were performed in triplicate. Identical setups in the absence of bacteria served to quantify the glucose amounts transported by diffusion or capillary forces along glass fibers from patch 1 to patch 2 in the presence and absence of glass fibers after 10 d. To do so, patch 2 was weighed and transferred into a 15 mL Falcon tube (Corning, USA). Subsequently, 5 mL of hot water and 3 mL of 1 N HCl were added and the tube vigorously shaken for 5 minutes. The suspension was then filtered and the glucose concentration measured by pulsed amperometric detection (PAD).

After harvest, patch 2 (Fig. 1A) was transferred to sterile 10 mL glass tubes containing 2 mL of 1x PBS buffer pH 7.4, the tubes vortexed for 1 min and then sonicated (3 times for 30 s with 1 min pause each time) in a sonication bath (SONOREX SUPER RK 255 H, BANDELIN, Germany) at 35 kHz ultrasound wave. The cell suspension was then transferred into 2 mL tubes and centrifuged at 16,000 × g and 4 °C (Eppendorf, Germany) to pellet the cells. The cells were lysed in 50 µL lysis buffer (10 mM Tris; 5 mM EDTA; 0.29% NaCl; and 0.4% SDS) and processed as previously described [40]. A Bradford assay was used to estimate protein concentrations. For the SDS-PAGE, 75 µg of protein was precipitated with an equal volume of 20% trichloroacetic acid (TCA, Sigma) for 16 h at 4 °C. The proteins were then separated in a 12% acrylamide gel using the HoeferTM mini–Vertical Electrophoresis system (Thermo Fisher Scientific, USA) and the Laemmli-buffer system [40]. The protein bands were visualized by staining the gel in colloidal Coomassie Brilliant blue G-250 (Roth, Kassel, Germany). Protein bands were excised from the gel and subjected to in-gel tryptic digestion as described previously [40]. Tryptic peptides of all gel slices were then extracted for LC-MS/MS analysis.

Peptides were analyzed by LC-MS/MS as described before [33]: after tryptic digestion of the samples, the resulting peptides were subjected to a shotgun proteomics workflow (nanoLC-MS/MS). Briefly, peptide lysate was injected into a nanoLC system (UltiMate 3000, Dionex, Thermo Fisher Scientific). Peptides were first trapped for 3 min on a C18-reverse phase trapping column (Acclaim PepMap 100, 75 µm × 2 cm, particle size 3 µM, nanoViper, Thermo Fisher Scientific), followed by separation on a C18-reverse phase analytical column (Acclaim PepMap 100, 75 µm × 25 cm, particle size 3 µM, nanoViper, Thermo Fisher Scientific) with a solvent flow-rate of 300 nL/min and a column temperature of 35 °C. Eluting peptides were ionized by a nano ion source (Advion, TriVersa Nanomate, Ithaca, NY, USA) and analyzed at the Q Exactive HF mass spectrometer (Thermo Fisher Scientific) (cf. SI S1.2 for detailed description). Proteome Discoverer (v1.4.1.14, Thermo Scientific) was used for protein identification and the acquired MS/MS spectra were searched with the SequestHT algorithm against the UniProt reference proteome of Sphingobium sp. and Fusarium equiseti. Trypsin served as cleavage enzyme, allowing a maximum of two miss cleavages. The precursor mass tolerance (MS) was set to 10 ppm, and the fragment mass tolerance (MS/MS) was 0.05 Da. Carbamidomethylation of cysteine was considered as fixed and oxidation of methionine was set as dynamic modification. Peptide spectrum matches (PSM’s) were validated using percolator with a false discover rate (FDR) < 1% and quality filtered for XCorr ≥2.25 [for charge state +2] and ≥2.5 [for charge state +3]. Identified proteins were grouped by applying the strict parsimony principle. Identification of 13C-labeled peptides and quantification of 13C incorporation was done by comparing measured and expected isotopologue patterns, chromatographic retention time, and fragmentation patterns using MetaProSIP as previously described [33]. Proteome data from the genomes of Sphingobium sp. S8 and Fusarium equiseti K3 were used as a decoy input file in the MetaProSIP pipeline. Data was evaluated by calculating the relative isotope abundance (RIA) and the labeling ratio (LR) of the labeled peptides.

Similar setups as described above in presence of unlabeled glucose were used to quantify the effects of mycelial C-substrate transfer on bacterial degradation of α-, β-, γ-, and δ-HCH isomers, where 10 µL of a HCH mix dissolved in acetone containing 25 µg of each isomer and 50 µL (i.e., ≈3.6 × 107 cells) of bacterial cells suspended in MSM media were simultaneously placed to patch 2 (Fig. 1A). Four circular agar patches containing activated carbon (60 mg mL−1) were included in the microcosms to decrease the gaseous phase concentration of HCH. After a 10-d incubation period, the patches (patch 1–3 and the activated carbon agar patches) were separately transferred into clean 40 mL screw cap septum glass vials (Thermo Scientific, USA), crushed with 10–15 g of dried sodium sulfate, and then extracted with hexane/acetone (3:1, v/v) containing 0.31 µg L−1 dichlorodiphenyldichloroethylene (DDE) as an internal standard. One hundred µL of the extract was mixed with 10 µL of 4-chlorobenzotrichloride as injection standard and analyzed by GC-ECD (HP 7890 Series GC, Agilent, USA) with a 20 m 0.18 mm (18 µm film thickness) HP5MS capillary column (Agilent, USA). Helium was used as a carrier gas at a flow rate of 1 mL min−1. Quadruplicate setups including an abiotic setup, a fungus-only, and a bacteria-only setup, were also tested as controls (cf. Fig. S1.2).

To account for possible variations of the amounts of the HCH isomer added (e.g., due to volatilization, varying dissolution properties of the isomers), the amounts of the individual isomers recovered immediately after application onto patch 2 were taken as the initial amounts of HCH (At0). Subsequently, HCH amounts (Adeg, t10) as well asthe residual (Fres, t10) and degraded fractions (Fdeg, t10) of added HCH after 10 d were derived from Eqs. (1)–(4) (cf. SI S2.3 for details). At10 is the residual amounts of HCH, Aabiotic, t10 the residual HCH amounts in the abiotic control, AP1, t10 or P2, t10 or P3, t10 the residual HCH amounts on patches 1, 2 or 3, and AAC,t10 the residual HCH amount on activated charcoal.

Statistical testing of the experimental treatments (abiotic, fungal-monoculture, bacterial-monoculture, and fungal-bacterial co-culture; p < 0.05) was performed on residual HCH concentrations (At10) by analysis of variance (one-way ANOVA) and Tukey’s HSD (honestly significant difference) in R (v 4.2.0) using multcompView (v 0.1–8) [41], dblyr (v 1.0.10) [42], and ggplot2 (v 3.4.0) [43] packages. The significance of differences between individual means was computed using the student’s t test.

In a spatially heterogeneous laboratory ecosystem, we assessed whether mycelial nutrient transfer from nutrient-rich to nutrient-deprived areas facilitates bacterial co-metabolic HCH degradation in a nutrient-deficient habitat. Using 13C-labeled fungal biomass and protein-SIP, we observed higher protein contents and stable isotope enriched proteins in Sphingobium sp. S8 in presence of fungal hyphae (protein content: 0.45 mg) as compared to controls in presence of glass fiber (0.37 mg) and the absence of hyphal or glass fiber transport (0.29 mg; SI Table S2.1). The relative isotope abundance (RIA, Fig. 2), the labeling ratio (LR; Fig. 2), and the shape of isotopic spectral patterns (Fig. 3) evidenced 13C enrichment of the total bacterial proteome (RIA = 10.9–93.9%; LR = 0.15–0.57) and selected proteins (RIA = 10.1–92%; LR = 0.17–0.59) (SI Table S2.2). Constitutively expressed bacterial proteins from core cellular and metabolic processes were selected out of 87 labeled proteins for assessing 13C incorporation in our treatments (Fig. 2 and Table S2.2): the regulatory protein chaperone DnaK [44] (RIA = 9.4%; LR = 0.13), the ATP synthase subunit-β, involved in carbon assimilation [45] (RIA = 10.4%; LR = 0.24), and the peptidoglycan-associated lipoprotein, a structural protein [46] (RIA = 14.2%; LR = 0.14). To quantify fungal impacts on the expression of central enzymes of bacterial HCH degradation, the 13C incorporation into two enzymes responsible for the initial steps of HCH degradation was quantified (Figs. 2, 3); i.e. the haloalkane dehalogenase (LinB; RIA = 9.5%; LR = 0.35) and the 2,5-dichloro-2,5-cyclohexadiene-1,4-diol dehydrogenase (Lin C; RIA = 7.1%; LR = 0.34). LinB and LinC exhibited significant labeling and hence were built on 13C enriched substrates upon hyphal transfer. As expected, observed RIA (88% – 96%) and LR (0.38–0.82) values in hyphal nutrient transfer experiments were different (p < 9.8 × 10–5) in cells cultivated on 13C glucose agar. Control experiments applying glass fibers to link patches 1 and 2 instead of mycelia were further performed assuming that abiotic fibers may not actively transfer 13C glucose and, hence, would not lead to increased RIA and LR. Likely due to capillary effects, efficient glucose transport to patch 2 was observed leading to clear glucose enrichment (≈24.8 mg) as found in additional abiotic glucose transport experiments (Fig. 1C). By contrast, no hyphal glucose transfer to bacteria-free patch 2 (< 0.14 mg) was measured (Fig. 1C), hence suggesting that 13C enrichment of the proteins in the fungal transfer experiments was due to fungal metabolites. Such assumption is further backed by RIA and LR values of glass fiber experiments, where RIA (86–97%) was similar to bacterial cultivation experiments on 13C glucose agar. As RIA reflects the extent to which a labeled substrate is used for biomass production [33], high RIA values in glass fiber transfer and 13C glucose growth experiment reflect the direct assimilation of 13C glucose by the bacteria [47]. Lower RIA values, however, were observed in presence of mycelia as a likely result of the assimilation of partially labeled fungal metabolites [33, 47]. The LR in the mycelial transfer experiments were similar or lower than LR in the glass fiber transfer experiments (LR = 0.15–0.34; Fig. 2) yet clearly lower than on 13C glucose agar (LR = 0.38–0.82). Peptide mass spectra patterns also reflected 13C protein incorporation (Fig. 3) and allowed for the visual distinction between direct assimilation of the labeled substrate (Fig. 3B, E, H, K, N & C, F, I, L, O) or utilization from partially labeled metabolites (Fig. 3A, D, G, J, M) secreted by other organisms [33]. The skewed isotopic incorporation pattern (Fig. 3A, D, G, J, M) in presence of mycelia, therefore, suggests the assimilation of partially labeled substrates [47].

Box plots of relative isotope abundance (RIA, A, C, E, G, I) and labeling ratio (LR, B, D, F, H, J) of housekeeping and key HCH degrading proteins of Sphingobium sp. S8 in presence of mycelia, glass fibers, and control (bacteria grown on 10% w/v 13C glucose agar) experiments. A, B: Chaperone DnaK protein, (C, D): ATP synthase β-subunit, (E, F): peptidoglycan associated lipoprotein, (G, H): haloalkane dehalogenase (LinB), and (I, J): dichloro-2,5-cyclohexadiene-1,4-diol dehydrogenase (LinC) respectively. Different lowercase letters indicate statistically significant difference between the treatment groups (p < 0.05, n = 3), while the whiskers represent data variability outside the upper and the lower quantiles.

A–C peptide sequence TTPSIVAFTK of Chaperone DnaK protein; (D–F): peptide sequence LVLEVAQHLGENTVR of ATP synthase β subunit; (G–I): peptide sequence VTIEGHADER of peptidoglycan associated lipoprotein; (J, K, L): peptide sequence DYAGWLSESPIPK of haloalkane dehalogenase; (M, N, O) peptide sequence RGANVVVADIVEAAAEETVR of 2,5-dichloro-2,5-cyclohexadiene-1,4-diol dehydrogenase.

We further quantified the effects of mycelial C-substrate transfer on bacterial HCH degradation. Lower residual HCH fractions (Fres, t10; Eq. 4) of all HCH isomers were found in fungal-bacterial co-cultures (Fres, t10 = 19%–45%) than in bacterial (Fres, t10 = 50%–72%) or fungal (Fres, t10 = 92%–99%) monocultures after 10 days of incubation (SI Table S2.4). No difference in Fres, t10 (p < 0.05) was observed between the fungal monoculture and abiotic control setups (SI, Fig. S2.1) suggesting that increased HCH degradation in presence of fungal-bacterial co-cultures was due to increased bacterial activity. To distinguish between abiotic and biotic losses, the HCH fractions lost by microbial activity (Fdeg, t10) were calculated (Eq. 3) revealing that 71% (α-HCH), 57%, (β-HCH), 75% (γ-HCH), and 81% (δ-HCH) were degraded by fungal-bacterial activity (Fig. 4). For the α-, β-, γ-, and δ-HCH isomers this results in ca. 1.9-, 1.1-, 1.6- and 1.6-fold better degradation than by a bacterial monoculture (Fig. 4, Table S2.6) and reveals significant (p < 0.05) influence of C-substrate transfer on the degradation of the α-, γ-, and δ-HCH, yet not of β-HCH. Batch biodegradation experiments performed in parallel revealed that Sphingobium sp. S8 co-metabolically degraded β-HCH less efficiently (~59% of 34.4 µM β-HCH in 24 h; SI S1.4 & Fig S2.2) than the other isomers (~100 % removal; SI S1.4 & Fig S2.2). The low degradation efficiency of β-HCH that is independent of additional substrates may hence explain the poor effect of hyphal transfer on the β-HCH biodegradation.

The HCH fractions degraded in presence of, a fungal monoculture (A), bacterial monoculture (B), and fungal-bacterial co-culture (C). Different lowercase letters indicate statistically significant difference between the treatment groups (p < 0.05, n = 4), while the whiskers represent data variability outside the upper and the lower quantiles.

Using 13C-labeled fungal biomass and protein-SIP, we traced the incorporation of 13C fungal metabolites into proteins of HCH degrading Sphingobium sp. S8. Our findings suggest that mycelial metabolite transfer to nutrient-deficient habitats promoted bacterial growth and co-metabolic HCH in a spatially structured model ecosystem. The 13C-metabolite transfer was detected by relative isotope abundance, the labeling ratio, and the shape of isotopic mass distribution profiles of the bacterial proteins. Mycelial transport typically exhibited lower RIA than in fungus-free controls that provided 13C glucose at different bioavailability (Fig. 1), either by glass fiber transfer or upon growth on 13C glucose agar (Fig. 2A, C, E, G, I). The LR upon metabolite transport were similar or even higher (Chaperone DnaK (Fig. 2D) and LinC (Fig. 2J)) than in glass fiber experiments. As the LR reflects the extent of protein turnover and hence bacterial growth [33], suggests efficient activation of the bacterial metabolism and HCH biotransformation by fungal metabolites. 13C incorporation into key enzymes of co-metabolic HCH degradation (i.e., LinB and LinC) thereby provides a direct link between hyphal metabolite transport and co-metabolic bacterial HCH degradation in fragmented habitats (Fig. 4). Our study hence provides for the first time direct molecular evidence for the role of fungal-bacterial interactions for co-metabolic contaminant degradation in an otherwise nutrient-deprived environment. In natural environments fungal-bacterial associations [18, 22] can develop neutral, competitive (antagonistic) or cooperative [22, 48] relationships. In a synergistic relationship as in our experiment, fungal compounds have been found to stimulate bacterial growth in the hyphosphere [14, 17, 22, 49] by providing low molecular weight organic acids (LMWOAs) such as acetic acid and oxalic acid [50], carbohydrates, polyols, peptides, and amino acids [50, 51]. Although we didn’t characterize the metabolome of fungal strain K3, we speculate that similar metabolites may have led to the observed isotope HCH degradation effects. Observed active mycelial growth and the proximity of Sphingobium sp. S8 biofilms to fungal hyphae (Fig. 1B) suggests that C-substrate assimilation of strain S8 occurred by extracellular biotrophy [52] in the hyphosphere of Fusarium equiseti K3 as has been reported earlier in other systems [18, 53, 54]. Such an assumption is further supported by the fact that Sphingobium sp. S8 [37] and Fusarium equiseti K3 [36] had been isolated from the same organochloride pesticide-contaminated soil.

similar setups as used for protein-SIP experiments (Fig. S1.2), we determined mycelial transport effects on bacterial HCH degradation by quantifying residual concentrations of the α-, β-, γ-, and δ-isomers of HCH. After 10 days of incubation, we did not find HCH removal by Fusarium equiseti K3, yet clearly increased bacterial HCH degradation in presence of mycelia (Fig. 4 & Fig S2.1) led to 1.6–1.9-fold faster biodegradation of α-, γ-, and δ-HCH. The least benefit (1.1-fold increase) was observed for β-HCH. β-HCH is the most recalcitrant of the HCH isomers tested, due to its equatorially substituted chlorine atoms that form a barrier to the enzymatic dehydrochlorination preferring axial chlorine atoms [7, 8, 55, 56]. The Sphingobium sp. S8 could not effectively degrade HCH isomers in absence of an additional carbon source (Fig. 4). In the presence of 1 % glucose in liquid batch experiments, however, it showed 59–100 % removal of all HCH isomers (cf. SI Fig. S2.2) indicating its reliance on additional carbon sources for HCH degradation. Combining protein-SIP with the HCH degradation, our study, therefore, underpins the role of fungal exudates for efficient co-metabolic bacterial HCH degradation as also had been proposed for pesticide degradation in porous media [57, 58]. Co-metabolic degradation is thought to be a critical process for the removal of xenobiotic compounds at both trace and bulk concentrations [6, 59, 60]. Benimeli et al. (2007) for instance, demonstrated simultaneous consumption of glucose and γ-HCH by Streptomyces sp. M7 with increased γ-HCH removal at high glucose concentrations [61]. Similarly, Álvarez et al. (2012) demonstrated increased γ-HCH removal by Streptomyces strains in the presence of root exudates as an additional carbon source [62]. In addition, HCH degradation in Sphingomonas have mostly been conducted in the presence of varying percentages of Luria-Bertani (LB) broth or glucose as additional carbon sources [1, 63].

In our study, we observed commensal incorporation of fungus-derived 13C-labels into bacterial peptides. The presence of 13C labeled bacterial peptides indicated the availability of fungal exudates for bacterial uptake and hence nutrient transfer. The labeled bacterial peptides belonged to both housekeeping enzymes and enzymes relevant to HCH degradation. The labeling of bacterial haloalkane dehalogenase (LinB) and 2,5-dichloro-2,5-cyclohexadiene-1,4-diol dehydrogenase (LinC) involved in HCH degradation indicates the functional relevance of mycelial nutrient transfer to HCH degradation by bacteria. This was indicated by improved HCH degradation (p < 0.05) in the presence of fungus and bacteria compared to either bacteria (up to 2 times increase) or mycelia (up to 12 times increase) alone for all four HCH isomers. Our data thereby underpin the relevance of the mycosphere as a habitat for the efficient degradation of organic contaminants. Since co-metabolic contaminant degrading microorganisms in the mycosphere may not rely on the contaminants for growth, they also may degrade chemicals at minute concentrations (“micropollutants”) and down to proposed cleanup endpoints in the parts per trillion range [31].

The mass spectrometry proteomics data generated during the current study are available in the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD038571.

Mohn WW, Mertens B, Neufeld JD, Verstraete W, de Lorenzo V. Distribution and phylogeny of hexachlorocyclohexane-degrading bacteria in soils from Spain. Environ Microbiol. 2006;8:60–8.

Sharma P, Raina V, Kumari R, Malhotra S, Dogra C, Kumari H, et al. Haloalkane Dehalogenase LinB is responsible for β- and δ-hexachlorocyclohexane transformation in Sphingobium indicum B90A. Appl Environ Microbiol. 2006;72:5720–7.

Lal R, Dogra C, Malhotra S, Sharma P, Pal R. Diversity, distribution and divergence of lin genes in hexachlorocyclohexane-degrading sphingomonads. Trends Biotechnol. 2006;24:121–30.

Singh A, Lal R. Sphingobium ummariense sp. nov., a hexachlorocyclohexane (HCH)-degrading bacterium, isolated from HCH-contaminated soil. Int J Syst Evol Microbiol. 2009;59:162–6.

Singh BK, Kuhad RC. Biodegradation of lindane (γ-hexachlorocyclohexane) by the white-rot fungus Trametes hirsutus. Lett Appl Microbiol. 1999;28:238–41.

Kumar D, Pannu R. Perspectives of lindane (γ-hexachlorocyclohexane) biodegradation from the environment: a review. Bioresour Bioprocess. 2018;5:29.

Lal R, Pandey G, Sharma P, Kumari K, Malhotra S, Pandey R, et al. Biochemistry of microbial degradation of hexachlorocyclohexane and prospects for bioremediation. Microbiol Mol Biol Rev. 2010;74:58–80.

Nayyar N, Lal R. Hexachlorocyclohexane contamination and solutions: Brief history and beyond. J Bioremed Biodeg. 2016;07:338.

Bumpus JA, Tien M, Wright D, Aust SD. Oxidation of persistent environmental pollutants by a white rot fungus. Science. 1985;228:1434–6.

Phillips TM, Seech AG, Lee H, Trevors JT. Biodegradation of hexachlorocyclohexane (HCH) by microorganisms. Biodegr. 2005;16:363–92.

Dogra C, Raina V, Pal R, Suar M, Lal S, Gartemann KH, et al. Organization of lin genes and IS6100 among different strains of hexachlorocyclohexane-degrading Sphingomonas paucimobilis: Evidence for horizontal gene transfer. J Bacteriol. 2004;186:2225–35.

Pal R, Bala S, Dadhwal M, Kumar M, Dhingra G, Prakash O, et al. Hexachlorocyclohexane-degrading bacterial strains Sphingomonas paucimobilis B90A, UT26 and Sp+, having similar lin genes, represent three distinct species, Sphingobium indicum sp. nov., Sphingobium japonicum sp. nov. and Sphingobium francense sp. nov., and reclassification of [Sphingomonas] chungbukensis as Sphingobium chungbukense comb. nov. Int J Syst Evol Microbiol. 2005;55:1965–72.

Rijnaarts HHM, Bachmann A, Jumelet JC, Zehnder AJB. Effect of desorption and intraparticle mass transfer on the aerobic biomineralization of. alpha-hexachlorocyclohexane in a contaminated calcareous soil. Environ Sci Technol. 1990;24:1349–54.

Harms H, Schlosser D, Wick LY. Untapped potential: exploiting fungi in bioremediation of hazardous chemicals. Nat Rev Microbiol. 2011;9:177–92.

Wick LY Bioavailability as a microbial system property: Lessons learned from biodegradation in the mycosphere. In: Ortega-Calvo JJ, Parsons JR, (eds). Bioavailability of organic chemicals in soil and sediment. Cham: Springer International Publishing; 2020 p. 267–89.

Jennings DH. Translocation of Solutes in fungi. Biol Rev. 1987;62:215–43.

Nazir R, Warmink JA, Boersma H, van Elsas JD. Mechanisms that promote bacterial fitness in fungal-affected soil microhabitats. FEMS Microbiol Ecol. 2010;71:169–85.

Worrich A, Stryhanyuk H, Musat N, König S, Banitz T, Centler F, et al. Mycelium-mediated transfer of water and nutrients stimulates bacterial activity in dry and oligotrophic environments. Nat Commun. 2017;8:15472.

Deveau A, Bonito G, Uehling J, Paoletti M, Becker M, Bindschedler S, et al. Bacterial–fungal interactions: ecology, mechanisms and challenges. FEMS Microbiol Rev. 2018;42:335–52.

Kohlmeier S, Smits THM, Ford RM, Keel C, Harms H, Wick LY. Taking the fungal highway: Mobilization of pollutant-degrading bacteria by fungi. Environ Sci Technol. 2005;39:4640–6.

Wick LY, Furuno S, Harms H Fungi as transport vectors for contaminants and contaminant-degrading bacteria. In: Timmis KN, (eds). Handbook of hydrocarbon and lipid microbiology. Berlin, Heidelberg: Springer Berlin Heidelberg; 2010. p. 1555–61.

Espinosa-Ortiz EJ, Rene ER, Gerlach R. Potential use of fungal-bacterial co-cultures for the removal of organic pollutants. Crit Rev Biotechnol. 2021;0:1–23.

Haq IU, Zhang M, Yang P, van Elsas JD The interactions of bacteria with fungi in soil. In: Adv Appl Microbiol. Elsevier; 2014. p. 185–215.

Furuno S, Foss S, Wild E, Jones KC, Semple KT, Harms H, et al. Mycelia promote active transport and spatial dispersion of polycyclic aromatic hydrocarbons. Environ Sci Technol. 2012;46:5463–70.

Schamfuß S, Neu TR, van der Meer JR, Tecon R, Harms H, Wick LY. Impact of mycelia on the accessibility of fluorene to PAH-degrading bacteria. Environ Sci Technol. 2013;47:6908–15.

Schamfuß S, Neu TR, Harms H, Wick LY. A Whole Cell Bioreporter approach to assess transport and bioavailability of organic contaminants in water unsaturated systems. J Vis Exp. 2014;24:52334.

You X, Kallies R, Kühn I, Schmidt M, Harms H, Chatzinotas A, et al. Phage co-transport with hyphal-riding bacteria fuels bacterial invasion in a water-unsaturated microbial model system. ISME J. 2022;16:1275–83.

Wick LY, Remer R, Würz B, Reichenbach J, Braun S, Schäfer F, et al. Effect of fungal hyphae on the access of bacteria to phenanthrene in soil. Environ Sci Technol. 2007;41:500–5.

Boer W, de, Folman LB, Summerbell RC, Boddy L. Living in a fungal world: impact of fungi on soil bacterial niche development. FEMS Microbiol Rev. 2005;29:795–811.

Frey‐Klett P, Garbaye J, Tarkka M. The mycorrhiza helper bacteria revisited. N. Phytol. 2007;176:22–36.

Hazen TC Cometabolic Bioremediation. In: Timmis KN, editor. Handbook of hydrocarbon and lipid microbiology. Berlin, Heidelberg: Springer; 2010. p. 2505–14.

Jehmlich N, Vogt C, Lünsmann V, Richnow HH, von Bergen M. Protein-SIP in environmental studies.Curr Opin Biotechnol.2016;41:26-33.

Sachsenberg T, Herbst FA, Taubert M, Kermer R, Jehmlich N, von Bergen M, et al. MetaProSIP: automated inference of stable isotope incorporation rates in proteins for functional metaproteomics. J Proteome Res. 2015;14:619–27.

Taubert M, Jehmlich N, Vogt C, Richnow HH, Schmidt F, von Bergen M, et al. Time resolved protein-based stable isotope probing (Protein-SIP) analysis allows quantification of induced proteins in substrate shift experiments. Proteomics 2011;11:2265–74.

Lünsmann V, Kappelmeyer U, Benndorf R, Martinez-Lavanchy PM, Taubert A, Adrian L, et al. In situ protein-SIP highlights Burkholderiaceae as key players degrading toluene by para ring hydroxylation in a constructed wetland model: Protein-SIP in a toluene-degrading wetland. Environ Microbiol. 2016;18:1176–86.

Khan N, Toscan RB, Lunayo A, Wamalwa B, Muge E, Mulaa FJ, et al. Draft genome sequence of Fusarium equiseti K3, a fungal species isolated from hexachlorocyclohexane-contaminated soil. Microbiol Resour Announc. 2021;10:e00885–21.

Khan N, Brizola Toscan R, Lunayo A, Wamalwa B, Muge E, Mulaa FJ, et al. Draft genome sequences of two Sphingobium species associated with hexachlorocyclohexane (HCH) degradation isolated from an HCH-Contaminated Soil. Microbiol Resour Announc. 2022;11:e00886–21.

Yang CH, Menge JA, Cooksey DA. Mutations Affecting Hyphal Colonization and Pyoverdine Production in Pseudomonads Antagonistic toward Phytophthora parasitica. Appl Environ Microbiol. 1994;9:473–81.

Senoo K, Wada H. Isolation and identification of an aerobic γ-HCH-decomposing bacterium from soil. Soil Sci Plant Nutr. 1989;35:79–87.

Jehmlich N, Schmidt F, Taubert M, Seifert J, Bastida F, von Bergen M, et al.Protein-based stable isotope probing.Nat Protoc.2010;5:1957-66.

Graves S, Dorai-Raj HPP and LS with help from S. multcompView: Visualizations of paired comparisons. 2019. Available from: https://CRAN.R-project.org/package=multcompView

Wickham H, François R, Henry L, Müller K, RStudio. dplyr: A Grammar of data manipulation. 2022. Available from: https://CRAN.R-project.org/package=dplyr

ggplot2: Create elegant data visualisations using the grammar of graphics—ggplot2-package. Available from: https://ggplot2.tidyverse.org/reference/ggplot2-package.html

Schramm FD, Heinrich K, Thüring M, Bernhardt J, Jonas K. An essential regulatory function of the DnaK chaperone dictates the decision between proliferation and maintenance in Caulobacter crescentus. PLoS Genet. 2017;13:e1007148.

Li SX, Wu HT, Liu YT, Jiang YY, Zhang YS, Liu WD, et al. The F1Fo-ATP synthase β subunit is required for Candida albicans pathogenicity due to its role in carbon flexibility. Front Microbiol. 2018;9:1025.

Godlewska R, Wiśniewska K, Pietras Z, Jagusztyn-Krynicka EK. Peptidoglycan-associated lipoprotein (Pal) of Gram-negative bacteria: function, structure, role in pathogenesis and potential application in immunoprophylaxis. FEMS Microbiol Lett. 2009;298:1–11.

Taubert M, Vogt C, Wubet T, Kleinsteuber S, Tarkka MT, Harms H, et al. Protein-SIP enables time-resolved analysis of the carbon flux in a sulfate-reducing, benzene-degrading microbial consortium. ISME J. 2012;6:2291–301.

Little AEF, Robinson CJ, Peterson SB, Raffa KF, Handelsman J. Rules of engagement: Interspecies interactions that regulate microbial communities. Annu Rev Microbiol. 2008;62:375–401.

Morris BEL, Henneberger R, Huber H, Moissl-Eichinger C. Microbial syntrophy: interaction for the common good. FEMS Microbiol Rev. 2013;37:384–406.

Van Hees PAW, Rosling A, Essén S, Godbold DL, Jones DL, Finlay RD. Oxalate and ferricrocin exudation by the extramatrical mycelium of an ectomycorrhizal fungus in symbiosis with Pinus sylvestris. N Phytol. 2006;169:367–78.

Sun YP, Unestam T, Lucas SD, Johanson KJ, Kenne L, Finlay R. Exudation-reabsorption in a mycorrhizal fungus, the dynamic interface for interaction with soil and soil microorganisms. Mycorrhiza 1999;9:137–44.

Leveau JHJ, Preston GM. Bacterial mycophagy: definition and diagnosis of a unique bacterial-fungal interaction. N Phytol. 2008;177:859–76.

Bassler BL, Losick R. Bacterially speaking. Cell 2006;125:237–46.

Kim HJ, Boedicker JQ, Choi JW, Ismagilov RF. Defined spatial structure stabilizes a synthetic multispecies bacterial community. Proc Natl Acad Sci. 2008;105:18188–93.

Nayyar N, Sangwan N, Kohli P, Verma H, Kumar R, Negi V, et al. Hexachlorocyclohexane: persistence, toxicity and decontamination. Rev Environ Health. 2014;29:49–52.

Willett KL, Ulrich EM, Hites RA. Differential toxicity and environmental fates of hexachlorocyclohexane isomers. Environ Sci Technol. 1998;32:2197–207.

Ellegaard-Jensen L, Knudsen BE, Johansen A, Albers CN, Aamand J, Rosendahl S. Fungal-bacterial consortia increase diuron degradation in water-unsaturated systems.Sci Total Environ.2014;466–467:699–705.

Knudsen BE, Ellegaard-Jensen L, Albers CN, Rosendahl S, Aamand J. Fungal hyphae stimulate bacterial degradation of 2,6-dichlorobenzamide (BAM).Environment Pollut.2013;181:122–7.

Saez JM, Alvarez A, Fuentes MS, Amoroso MJ, Benimeli CS An Overview on Microbial Degradation of Lindane. In: Singh SN, (eds). Microbe-induced degradation of pesticides. Cham: Springer International Publishing; 2017. p. 191–212.

Nzila A. Update on the cometabolism of organic pollutants by bacteria. Environ Pollut. 2013;178:474–82.

Benimeli CS, González AJ, Chaile AP, Amoroso MJ. Temperature and pH effect on lindane removal by Streptomyces sp. M7 in soil extract. J Basic Microbiol. 2007;47:468–73.

Álvarez A, Yañez ML, Benimeli CS, Amoroso MJ. Maize plants (Zea mays) root exudates enhance lindane removal by native Streptomyces strains. Int Biodeterior Biodegrad. 2012;66:14–8.

Boltner D, Moreno-Morillas S, Ramos JL. 16S rDNA phylogeny and distribution of lin genes in novel hexachlorocyclohexane-degrading Sphingomonas strains. Environ Microbiol. 2005;7:1329–38.

We acknowledge the Department of Molecular Systems Biology (Proteomics Lab) at the Helmholtz Center for Environmental Research (UFZ) for providing their analytical facilities for this study. The work presented in this article is partly supported by Novo Nordisk Foundation grant NNF21OC0066551. Funding by Deutscher Akademischer Austausch Dienst (DAAD) and a grant from the International Foundation for Science (IFS grant W/5798-1) are also acknowledged. The authors thank Birgit Würz for her help with GC-MS analysis. We are also grateful to Jana Reichenbach and Rita Remer for their skilled experimental help. MvB and NJ are grateful for funding of the UFZ for the ProMetheus platform for proteomics and metabolomics.

Open Access funding enabled and organized by Projekt DEAL.

University of Nairobi, Department of Biochemistry, 00200-30197, Nairobi, Kenya

Nelson Khan, Edward Muge & Francis J. Mulaa

Helmholtz Centre for Environmental Research UFZ, Department of Environmental Microbiology, 04318, Leipzig, Germany

Nelson Khan & Luke Y. Wick

University of Nairobi, Department of Chemistry, 00200-30197, Nairobi, Kenya

Helmholtz Centre for Environmental Research UFZ, Department of Molecular Systems Biology, 04318, Leipzig, Germany

Martin von Bergen & Nico Jehmlich

German Centre for Integrative Biodiversity Research, (iDiv) Halle-Jena-Leipzig, Puschstraße 4, 04103, Leipzig, Germany

University of Leipzig, Faculty of Life Sciences, Institute of Biochemistry, Brüderstraße 34, 04103, Leipzig, Germany

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

NK and LYW designed the study. NK performed the experiments. NK and NJ performed the Protein SIP analysis. NJ and MvB interpreted the MetaProSIP data. NK, LYW and NJ wrote the manuscript. LYW, EM, BW, and FJM supervised the experiments. All authors read, revised, and approved the final manuscript.

Correspondence to Lukas Y. Wick.

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Khan, N., Muge, E., Mulaa, F.J. et al. Mycelial nutrient transfer promotes bacterial co-metabolic organochlorine pesticide degradation in nutrient-deprived environments. ISME J (2023). https://doi.org/10.1038/s41396-023-01371-7

DOI: https://doi.org/10.1038/s41396-023-01371-7

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Mycelial nutrient transfer promotes bacterial co-metabolic organochlorine pesticide degradation in nutrient-deprived environments | The ISME Journal

Deltamethrin Pesticide The ISME Journal (ISME J) ISSN 1751-7370 (online) ISSN 1751-7362 (print)