Gerhard P. Aigner
Veronika Peer
Martina Höckner- Department of Zoology, Center for Molecular Biosciences, University of Innsbruck, Innsbruck, Austria
Earthworms are affected by physical stress, like injury, and by exposure to xenobiotics, such as the toxic metal cadmium (Cd), which enters the environment mainly through industry and agriculture. The stress response to the single and the combination of both stressors was examined in regenerative and unharmed tissue of Lumbricus terrestris to reveal if the stress response to a natural insult like injury (amputation) interferes with Cd detoxification mechanisms. We characterized the roles of metallothionein 1 (MT1) and MT2 isoforms, heat shock protein 70 as well as immune biomarkers such as the toll-like receptors (TLR) single cysteine cluster TLR and multiple cysteine cluster TLR. The role of the activated transcription factors (ATFs) ATF2, ATF7, and the cAMP responsive element binding protein as putative regulatory intersection as well as a stress-dependent change of the essential trace elements zinc and calcium was analyzed. Phosphorylated AMP activated protein kinase, the cellular energy sensor, was measured to explore the energy demand, while the energy status was determined by detecting carbohydrate and protein levels. Taken together, we were able to show that injury rather than Cd is the driving force that separates the four treatment groups – Control, Cd exposure, Injury, Cd exposure and injury. Interestingly, we found that gene expression differed regarding the tissue section that was analyzed and we hypothesize that this is due to the migration of coelomocytes, earthworm immune cells, that take over a key role in protecting the organism from a variety of environmental challenges. Surprisingly, we discovered a role for MT1 in the response to multiple stressors and an isoform-specific function for the two newly characterized TLRs. In conclusion, we gathered novel information on the relation of innate immunity, wound healing, and Cd detoxification mechanisms in earthworms.
1 Introduction
Earthworms inhabit most of the Earth’s terrestrial ecosystems (1), and their soil-dwelling lifestyle makes them vulnerable to both environmental and anthropogenic sources of stress. They have therefore developed a strong stress response and a potent innate immune system. The latter is composed of a humoral and a cellular part, so-called coelomocytes, which consist of different cell types located in the coelomic cavity. Upon external stimulation coelomocytes are extruded through dorsal pores consequently covering and protecting the earthworm body from external challenges that can range from soil pathogens to xenobiotics (2, 3).
One of the major xenobiotics, which is brought into the environment is the toxic metal cadmium (Cd) with an estimated amount of >5.5 g Cd ha−1 yr−1 in Europe (4). Cd concentrations in European soils vary depending on the history of land use. In France, for example, Cd concentrations range from a median of 0.13 mg/kg dry soil for forest areas to a median of 1.3 mg/kg for soils derived from urban, mining, industrial or military areas, where locally extreme values of up to 162 mg/kg dry soil were detected (5). Cd is a proven carcinogen and is able to substitute zinc (Zn) and calcium (Ca) due to the same charge and ionic radius (ionic mimicry), which allows Cd to hijack transport mechanisms across the membrane (6). Subsequently, Cd exhibits both immune- and cytotoxic effects, interfering with a broad range of processes from epigenetic mechanisms to energy metabolism (7–9). However, organisms have developed protection strategies to counteract adverse effects of cellular stress, such as modulating the expression of heat shock proteins (HSPs) (10, 11) and more specifically in response to Cd, the production of so-called metallothioneins (MTs) (12). Three isoforms (MT1-3) have been identified in Lumbricus rubellus. The most studied isoform is the Cd-inducible MT2, which is responsible for Cd detoxification. On the contrary, MT1 has been suggested to function in metal homeostasis (13) and a third isoform, MT3, has been found mainly in embryonic tissues (14). Recently, also non-MT proteins have been discovered that function in Cd detoxification (15).
Coelomocytes take over key roles after physical stress (like injury), in regeneration, and wound healing. The regenerative capacity differs largely in the annelid phylum (16). It seems that also gene expression during regeneration is species-specific as shown in Lampito mauritii, Drawida calebispecies (17), and Eisenia fetida (18). Furthermore, earthworm transcriptome studies reveal transcripts that cannot be assigned to any homologs from other species, although there have been overlaps with genes expressed in vertebrate model organisms like zebrafish and mouse (17, 18). Therefore, we see a high potential in studying earthworms in more detail to unearth novel regenerative mechanisms. Regeneration can be defined as periodic replacement of cells and tissue, the regeneration of lost body parts or tissue is called post-traumatic regeneration (19). According to Kostyuchenko and Kozin, the latter includes processes that can partially overlap, like closure and wound healing; immune and/or nonspecific defense response; recruitment of cells that are to form a regenerate; growth of the regenerate; patterning and differentiation. Annelid regeneration generally involves cell proliferation and the formation of a blastema, which is a mass of undifferentiated cells. There are two ways to close a wound described in E. fetida, namely alimentary tract eversion and formation of a coelomic plug (20). A study on wounded Dendrobaena veneta revealed that irradiation led to wound closure by muscle contraction (21).
It has already been shown that coelomocytes accumulate in the regenerating blastema in Eisenia andrei. The latter study also revealed that specific subsets of coelomocytes form clusters in the blastema, while others were not only present in the blastema but also in the body wall, alimentary canal and in the muscle, probably to remove damaged cells (22). A previously published study suggested that wound healing and the detoxification of Cd are cellular processes that interact with each other, since in amputated earthworms the expression of MT2 was drastically decreased (23). Few information is available on Lumbricus terrestris (24) wound healing or regeneration specifically related to Cd exposure. By introducing a three-sided cut into the integument posterior to the clitellum, it has been suggested that Cd impairs enzymes that are responsible for wound healing in L. terrestris (25). There are several studies on the effects of polychlorinated biphenyls (PCBs), which showed delayed healing effects and accelerated allograft rejection (26, 27).
Pattern recognition receptors (PRRs) activate the immune system as a first response to and protection against pathogens (28). The best known PRRs are the Toll-like receptors (TLRs), which consist of three domains – a leucine-rich repeat extracellular domain, a transmembrane domain, and an intracellular Toll/IL-1 receptor (TIR) domain (29). TLR-mediated downstream signaling cascades are known to activate transcription factors of the cAMP responsive element binding protein (CREB)/activating transcription factor (ATF) family, subsequently leading to the induction of genes related to innate immunity and cell survival (30, 31). In earthworms, the exact number of CREB/ATF proteins and TLRs is unknown. Two TLRs have been characterized in E. andrei, a single cysteine cluster (scc) TLR and a multiple cysteine cluster (mcc) TLR, which are mainly expressed in reproductive organs and could be linked to gregarine infection (32, 33). Gram-positive bacteria induced the expression of sccTLR in E. andrei coelomocytes (32).
The effect of Cd on energy allocation has been evaluated in Enchytraeus albidus, showing that after four to eight days, carbohydrate reserves were significantly reduced, however, protein levels showed no change (34). AMP activated protein kinase (AMPK), the metabolic checkpoint in the cell, could provide valuable clues about the metabolic state of a tissue. In the activated form (phosphorylated) it is involved in regulating cell growth, apoptosis, glucose and lipid metabolism (35). Upon energy stress AMPK inhibits anabolic or stimulates catabolic pathways and it is also known to increase the number of mitochondria (36, 37). Little information is available on the energy demands in earthworms during wound healing or regenerative processes.
The Cd stress response, PRRs, regulatory factors, essential trace elements as well as the energy metabolism were evaluated to reveal factors that might be related to both stressors, revealing a putative link between wound healing and detoxification mechanisms. Furthermore, as mentioned above, we analyzed the suppression of Cd-induced MT2 expression in amputated earthworms as previously suggested (23).
2 Materials and methods
2.1 Experimental setup and sampling
The earthworm species L. terrestris and soil (mixture of humus and peat) were ordered from Proinsects GmbH (Minden, Germany). For the experimental setup, amputated and unharmed earthworms were kept in Cd-contaminated and non-contaminated soil, leading to a total of four treatment groups. Earthworms were amputated by excising ~1 cm of their posterior end with a sharp scalpel directly before the start of the experiment. The group Cd_Cut consisted of amputated earthworms that were kept in soil spiked with 50 mg CdCl2/kg dry soil, whereas unharmed individuals in Cd-contaminated soil were defined as the Cd group. Amputated individuals in non-contaminated soil were considered the Cut group and unharmed earthworms in non-contaminated soil served as control. All groups were fed once a week with 0.9 g of horse manure and were kept at 75% soil humidity, 15°C room temperature and a day/night cycle of 12/12 h. Before the start of the experiment, all earthworms were stored in control soil for at least two weeks for acclimation.
Tissue samples were taken once a week over a total of three weeks. In amputated earthworms the regenerative tissue (Tissue 1) and the anterior section (Tissue 2) were sampled. In unharmed earthworms the tissue sections corresponding to the Cut groups were sampled. The number of biological replicates is shown in Supplementary Table S1. Tissue samples for RNA extraction were stored in 1 ml absolute Ethanol (Merck, Darmstadt, Germany) at -20°C. For Zn quantification and protein extraction, tissue samples were shock frozen in liquid nitrogen and stored at -80°C until further use. For Ca quantification, tissue was processed immediately after sampling.
2.2 RNA extraction and cDNA synthesis
For RNA extraction, each tissue sample was placed in 1.5 ml screw-cap tubes containing 1 ml TRI REAGENT® (Sigma Aldrich, St. Louis, MO, USA) and RNase-free glass beads. Tissue was disrupted using a FastPrep-24™ 5G benchtop homogenizer (MP Biomedicals, Santa Ana, CA, USA). The following settings were applied for homogenization: three cycles of 40 s at 6 m/s with a pause of 120 s between each cycle. RNA extraction and DNA digestion were performed as described in the user manual. DNA was digested using DNAse I (Thermo Fisher Scientific, Waltham, MA, USA), including the RNAse inhibitor RiboLock (Thermo Fisher Scientific, Waltham, MA, USA). The quality of the extracted RNA was checked by gel electrophoresis to assess the integrity of the rRNA bands. RNA concentration was measured in technical triplicates using the Quant-iT™ RiboGreen™ RNA Assay Kit (Thermo Fisher Scientific, Waltham, MA, USA) on a Victor X4 plate reader (Perkin Elmer, Waltham, MA, USA). 450 ng of total RNA was reverse transcribed using RevertAid™ H Minus Reverse Transcriptase (Thermo Fisher Scientific, Waltham, MA, USA), Random Hexamer Primers (Thermo Fisher Scientific, Waltham, MA, USA), and RiboLock (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s instructions.
2.3 Sequence generation
For sequence generation of L. terrestris sccTLR, a previously published CD842 primer pair from D. veneta was used (38). PCRs were performed on a T 100 thermal cycler (Bio-Rad, USA) using Titanium Taq DNA polymerase (Takara, Japan). After gel electrophoresis, bands were excised and were purified using the QIAquick Gel Extraction Kit (Qiagen GmbH, Netherlands) according to the manufacturer’s instructions. TA cloning was conducted with the purified PCR products using the pGEM®-T vector systems (Promega, USA) as described in the manual. Colony PCR was performed using SP6 and T7 primer sequences. PCR products were purified (QIAquick PCR Purification Kit, Qiagen, Netherlands) and were then sent for sequencing (Microsynth, Switzerland).
We used homolog sequences to search the NCBI Blast+ assembled RNA-Seq database of L. terrestris to find sequences for mccTLR, MT1, CREB, ATF2, and ATF7 and to design primers (Supplementary Table S2). Amino acid sequences were analyzed using the SMART™ online protein domain characterization tool. Sequences (Supplementary Table S3) were published in Genbank under the following Accession Numbers: sccTLR (OR596221), mccTLR (OR596222), MT1 (OR596220), ATF2 (OR610157), ATF7 (OR610158), CREB (OR610159).
2.4 Quantitative real-time PCR
A primer matrix was applied to determine optimal primer concentrations and a standard curve was generated from a series of dilutions of known template concentrations to quantify absolute RNA copy numbers. Copy numbers were measured using a QuantStudio™ 3 Real-Time PCR system (Thermo Fisher Scientific, Waltham, MA, USA). The list of primer sequences used for qPCR is given in Supplementary Table S2. 5 µl of Power SYBR™ Green PCR Master Mix (Thermo Fisher Scientific, Waltham, MA USA) was used per sample. CREB, ATF2, ATF7, MT2, sccTLR and mccTLR were amplified as follows: 1 µl 10 × BSA, 1 µl forward primer (9 µM), 1 µl reverse primer (9 µM) and 2 µl cDNA. MT1 was amplified as described above, but using a reverse primer concentration of 3 µM. qPCRs were conducted in technical triplicates using the following settings: 50°C for 2 min; 95°C for 10 min; 40 repeats of 95°C for 15 s and 60°C for 1 min. Absolute copy numbers per 100 ng RNA were determined from the standard curve. Values were normalized to the threshold. Values below the detection range of our standard curve were set to zero.
2.5 Zinc quantification
Total Zn levels were quantified using the Zinc Quantification Kit (Fluorometric) (ab176725, Abcam®, United Kingdom). In detail, tissue samples were weighed and placed in microcentrifuge screw-cap tubes containing 250 µl EDTA-free lysis buffer [50 mM Tris, 100 mM NaCl, 1 mM sodium vanadate, 0.5 mM PMSF, 1 µg/ml aprotinin, 1 µg/ml leupeptin, 1 µg/ml pepstatin, 10 µM MG132, 0.5 mM DTT] and glass beads. Tissue was disrupted as described for the RNA extraction method. After homogenization, a centrifugation step was performed at 4°C for 1 min at the maximum speed to eliminate air bubbles. 250 µl of 7% TCA solution was added to the sample, vortexed and centrifuged at 4°C for 5 min at maximum speed. 400 µl of the supernatant were transferred to a new vial. For deprotonation, 10 µl of 100% TCA containing 1 M Na2CO3 was added to the sample, vortexed and incubated on ice for 5 min. Afterwards, pH neutrality was determined using pH test strips. Zn standard dilutions were measured at room temperature at Ex/Em = 485/525 nm using the EnSpire™ Multimode Plate Reader (PerkinElmer Inc, USA) to establish a standard curve. 50 µl of sample was applied for quantification. Data from different plates were normalized according to the Zn standards.
2.6 Calcium quantification
Free Ca levels were quantified using the Calcium Assay Kit (Colorimetric) (ab102505, Abcam®, United Kingdom) following the instructions for tissue samples with slight modifications: Tissue samples were weighed and were suspended in 500 µl Calcium Assay Buffer (provided in the kit) in 1.5 ml screw-cap tubes containing glass beads. Tissue was disrupted as described for the RNA extraction method. The optical density (OD) of the Ca standard dilutions was determined at room temperature at 575 nm using the EnSpire™ Multimode Plate Reader (Perkin Elmer Inc, USA) to establish a standard curve. 50 µl of sample was applied for quantification. Data from different plates were normalized according to the Ca standards.
2.7 Protein extraction and western blots
For total protein extraction, tissue samples were homogenized in 2 ml Eppendorf tubes containing 500 µl lysis buffer [25% glycerol, 1.5 mM MgCl2, 0.2 mM EDTA Titriplex III, 40 mM Hepes, 430 mM NaCl2] including protease inhibitors [1 mM sodium vanadate, 0.5 mM PMSF, 1 µg/ml aprotinin, 1 µg/ml leupeptin, 1 µg/ml pepstatin, 0.5 mM DTT] and glass beads using the Precellys® Evolution tissue homogenizer (Bertin Instruments, France) as follows: 5000 rpm, 2 x15 s with a 15 s pause in between, which was repeated after a 3 min pause on ice. The homogenates were centrifuged at 2348 rcf for 5 min at 4°C, the supernatant was transferred to a fresh tube, which was repeated once. Protein concentration was determined using a NanoDrop™ 2000 spectrophotometer (ThermoFisher Scientific, USA). 25 μg total protein for P-AMPK and 40 μg total protein for HSP70 were mixed 1:1 with 2x Laemmli buffer (Bio-Rad Laboratories, Inc., USA) and were heated at 95°C for 5 min. Samples were separated on a 12% Criterion TGX Stain-free precast gel (Bio-Rad Laboratories, Inc., USA) in a Criterion vertical electrophoresis cell (Bio-Rad Laboratories, Inc., USA) containing 1x Tris-glycin buffer for 90 min at 125 V. The PageRuler™ Prestained Protein Ladder (Thermo Fisher Scientific, USA) was used as a size standard. The gel was then UV activated for 1-2.5 min using a Chemidoc (Bio-Rad Laboratories Inc, USA) for total protein visualization. A PVDF blotting membrane (Bio-Rad Laboratories Inc, USA) was activated with methanol for 5 min and filter paper was soaked in 1x semidry buffer containing 20% Methanol. Then, the gel was blotted to the PVDF membrane for 30 min at 25 V using a Trans-Blot® Turbo Transfer System (Bio-Rad Laboratories Inc, USA). After transfer, the membrane was washed with 1x TBS-T (Tris buffered saline with Tween) and was then added to 1x TBS-T containing 5% milk powder for 1 h at room temperature on a shaker to block the remaining binding sites. The membrane was incubated overnight at 4°C on a shaker in 1x TBS-T containing 5% milk powder and the primary antibody (rabbit anti-phospho-AMPKα [Thr172] from Merck KGaA (Product Number 07-626), Germany, diluted 1:1000; mouse monoclonal anti-heat shock protein 70 from SIGMA Aldrich (Product Number H 5147), USA, diluted 1:2500). The membrane was rinsed twice with 1x TBS-T, followed by three 10 min washes with 1x TBS-T. The membrane was then incubated for 1 h on a shaker at room temperature in 1x TBS-T containing 5% milk powder with the secondary antibody (goat anti-rabbit IgG H&L (HRP), Abcam, UK; anti-mouse IgG (whole molecule) – peroxidase antibody produced in rabbit, Sigma-Aldrich, USA, both diluted 1:10000). Detection was performed using Amersham™ ECL Select™ Western Blotting Detection Reagent (Cytiva, USA) in a Chemidoc (Bio-Rad Laboratories, Inc., USA). To test specificity of the antibodies, mouse protein from NIH 3T3 cells was used as a positive control to make sure that the correct band was chosen for analysis in earthworm samples (Supplementary Figure S1). Total protein was quantified from the stain-free precast gels and used as loading control for normalization. The positive controls and an example for a loading control are shown in Supplementary Figure S1. All membranes that were used for analysis are shown in the Supplementary Material (Supplementary Figure S2). All samples from one treatment group and time point were loaded onto the same gel together with the respective control samples.
2.8 Protein & carbohydrate quantification
Tissue samples were taken from individuals of each treatment group, stored at -80°C and homogenized in Millipore grade water using the Precellys® Evolution Tissue Homogenizer (settings as in 2.7). 200 µl of homogenate was incubated with 30 µl 100% TCA for 10 min at -20°C. After centrifugation (3000 g, 10 min, 4°C) the supernatant was transferred to a new tube. The pellet was dissolved in 100 µl 5% TCA and centrifuged again under the same conditions. The supernatant was added to the first supernatant. The pellet was resuspended in 100 µl 1 M NaOH and the protein content was measured using a NanoDrop™ 2000 spectrophotometer (ThermoFisher Scientific, USA). Carbohydrates were determined according to a previous study (39) with the following specifications. D(+)-Glucose (Roth, Germany) was used as the standard. 150 µl H2SO4 and 30 µl 5% phenol were added directly to 50 µl supernatant and incubated at 90°C for 5 min. Measurements were accomplished at 490 nm using a Victor X4 plate reader (Perkin Elmer, Waltham, MA, USA), since most of the sugars can be detected at this wavelength.
2.9 Statistical analysis
Statistical analyses were performed using GraphPad Prism 9.5.1 software (GraphPad Software, San Diego, CA, USA). Technical outliers were identified and removed using the Grubbs test (p < 0.05). Normality was tested using the Kolmogorov-Smirnov test. For normally distributed data, ANOVA (including Bonferroni correction) or student’s t-test was performed. If the variance was significantly different (p < 0.05) Welch’s correction was performed. Kruskal-Wallace test (including Dunn`s test) or Mann-Whitney U test was calculated for non-normally distributed data. For comparing Tissue 1 and Tissue 2 samples mixed-effects analysis was performed in GraphPad Prism 9.5.1 software.
Principal component analysis (PCA) was performed in RStudio 2022.07.1 using the PCA script from the FactoMineR package. For visualization we used fviz_pca from the factoextra package. An ellipse with a confidence level of 0.95 was added. Correlation was accomplished using the cor script from the stats version 4.2.1. package. To overcome the issue of missing values for principal component and correlation analysis we used the mean values of each analysis, time point and treatment group.
3 Results and discussion
Many research disciplines consider earthworms as important model organisms, such as ecotoxicology and regenerative biology (40, 41). In the present study, we aimed to bridge both disciplines and investigated the effects of cadmium exposure, wound healing, and the combination of both stressors in the earthworm species L. terrestris. Moreover, we resolved the question whether detoxification and wound healing are cellular processes that interfere with each other, which has previously been proposed in a study showing that injury inhibits Cd-induced MT2 gene expression in L. terrestris and L. rubellus (23). However, the latter study did not distinguish between harmed and unharmed tissue, in contrast to the present study. Herein, we observed a similar effect in the anterior sections (Tissue 2), namely that Cd-induced MT2 expression was suppressed, but only in amputated earthworms.
PCA revealed a clear separation of the treatment groups along Dim 1 (47.8%) by the health status of the earthworms (amputated or unharmed) rather than by Cd exposure (Figure 1). Except for a significantly increased MT2 gene expression level, which is known to be the major detoxification mechanism of Cd (14) at week one, significantly decreased HSP70 and P-AMPK levels, and a significant increase in P-AMPK after two weeks, no other exclusively Cd-related effects were detected. PCA showed that P-AMPK, HSP70, MT1, and MT2 were the main factors to separate the groups along Dim 2 (Supplementary Table S4), explaining 19.2% of the variation. The factors metal and time did not contribute to the group’s separation (Supplementary Figure S3). The results of the correlation analysis are shown in Supplementary Figure S4.
Figure 1 Principal component analysis (PCA) based on the mean values of each analysis from amputated and unharmed L. terrestris. The percentage of variation between the treatment groups for dimension 1 (Dim1) and dimension 2 (Dim2) are shown. The ellipse represents a confidence level of 0.95.
Cd is generally known to be an immunotoxic agent (8). In a previous study, we were able to show that earthworm immune cells – coelomocytes - are affected regarding their metabolism (42), but also chromosomal aberrations and single-strand breaks have been observed upon Cd exposure (43). In earthworms, a dose-dependent change in the expression of antimicrobial peptides by Cd (44) has been shown, as well as changes in lysenin and coelomic cytolytic factor 1 (CCF-1) levels, two factors which are important in immune defense (45). A study on the blue mussel showed that Cd led to the upregulation of several TLRs (46). We identified two partial L. terrestris TLR sequences (Supplementary Table S3). Both, sccTLR and mccTLR, did not reveal Cd-dependent changes. SMART™ protein domain analysis revealed one sccTLR containing an intracellular Toll -interleukin 1-resistance (TIR) domain (141 aa), a transmembrane region (22 aa), a single leucine rich repeat C-terminal (LRR_CT) domain (54 aa) and six internal extracellular leucine-rich repeats (LRRs) (23-24 aa) (Supplementary Figure S5A). The second TLR was identified as a mccTLR, containing an intracellular TIR domain (144 aa), a transmembrane region (22 aa), two LRR_CT domains (48 aa and 60 aa) with two LRR-C-terminal domains (22 and 23 aa) and one LRR-N-terminal domain between them (38 aa), followed by 12 internal extracellular LRRs (19 - 23 aa) (Supplementary Figure S5C). Supplementary Figures S5B, S5D shows the comparison to the E. andrei TLRs. We observed a significant induction of sccTLR in amputated- in contrast to Cd-exposed earthworms (Figure 2A). Interestingly, we could hardly observe any gene expression of sccTLR in Tissue 2 (Figure 3A). On the other hand, mccTLR revealed constitutive gene expression in both tissue sections (Figure 3B) and a significant difference between the control and Cd_Cut group at week 1 (Figure 2B). Therefore, the time and dose of Cd exposure might play an important role regarding the change of immunity-related genes.
Figure 2 Absolute mRNA copy numbers of (A) single cysteine cluster toll-like receptor (sccTLR) and (B) multiple cysteine cluster toll-like receptor (mccTLR) in L. terrestris tissue from Control and exposed (50 mg CdCl2/kg dry soil (Cd)) individuals; regenerative tissue of amputated (Cut) and exposed amputated [50 mg CdCl2/kg dry soil (Cd_Cut)] individuals at week one (w1), week two (w2), and week three (w3). Asterisks indicate significant differences between treatments and the control of the respective week [p < 0.01 (**); p < 0.001 (***)]. Mean values ± SEM are presented.
Figure 3 Comparison of absolute mRNA copy numbers in Tissue 1 and Tissue 2 of (A) single cysteine cluster toll-like receptor (sccTLR); (B) multiple cysteine cluster toll-like receptor (mccTLR); (C) metallothionein 1 (MT1) and (D) metallothionein 2 (MT2) in L. terrestris tissue from Control; exposed [50 mg CdCl2/kg dry soil (Cd)]; amputated (Cut) and exposed amputated [50 mg CdCl2/kg dry soil (Cd_Cut)] individuals at week one (w1), week two (w2), and week three (w3). Grey bars represent gene expression in Tissue 1 (including regenerative tissue) and colored bars represent gene expression in Tissue 2 (anterior section, no regenerative tissue). Asterisks indicate significant differences between Tissue 1 and Tissue 2 using mixed-effects analysis [p < 0.05 (*)]. Mean values ± SEM are presented.
P-AMPK decreased after one week of Cd exposure, indicating that sufficient energy is on hand. This was followed by an increase in all treatment groups after two weeks, P-AMPK peaked after three weeks in the combined stress group, which hints towards low ATP levels and an elevated energy need (Figure 4A). It seemed that Cd had no effect on the energy metabolism. Interestingly, carbohydrate and protein levels showed a quite strong negative correlation (Supplementary Figure S4). Only the combined treatment group showed significant changes in carbohydrate and protein levels. In the Cd_Cut group, after one and two weeks, carbohydrate levels were low and protein levels were high, whereas after three weeks it was the opposite – we observed lower protein and higher carbohydrate levels (Figures 4C, D).
Figure 4 Comparison of fold change normalized to control of (A) phosphorylated AMP-activated protein kinase (P-AMPK); (B) heat shock protein-70 (HSP70); (C) total carbohydrates; (D) total protein in L. terrestris tissue from Control and exposed [50 mg CdCl2/kg dry soil (Cd)] individuals; regenerative tissue of amputated (Cut) and exposed amputated [50 mg CdCl2/kg dry soil (Cd_Cut)] individuals at week one (w1), week two (w2), and week three (w3). Asterisks indicate significant differences between treatments and the control of the respective week [p < 0.05 (*); p < 0.01 (**); p < 0.001 (***); p < 0.0001 (****)]. Mean values ± SEM error are presented.
Cd exposure indeed triggered a carbohydrate depletion after 8 days in study on E. albidus (34). In the present study, only the combined stress group revealed an increased energy demand and a depletion of carbohydrates at first. Interestingly, the respective protein levels showed an opposite effect. An important role of lipids cannot be excluded, which we, however did not determine herein. The situation changed after three weeks – carbohydrate levels seemed to recover while protein levels decreased. P-AMPK levels remained elevated after three weeks of combined exposure suggesting that only amputated earthworms exposed to Cd still had an increased energy need.
A previous study revealed that AMPK negatively regulates HSP70 expression in HepG2 cells and that Cd induces HSP70 expression (47). After one week, the most obvious result herein was that Cd caused a downregulation of HSP70. After two weeks and most obviously after three weeks, HSP70 levels were significantly increased in the combined stress group (Figure 4B).
The heat shock response is a protective mechanism under environmental stress that is involved in protein folding, degradation, etc. We observed a similar expression pattern between P-AMPK and HSP70 levels – if these two proteins are linked at a regulatory level remains to be elucidated. Recently, a review article revealed that Cd leads to an upregulation of HSP70 (48), however, there is also evidence that Cd decreases HSP70 levels like for example in Cyprinidae fish (49) and zebrafish (50). The latter study revealed that the Cd-induced decrease in HSP70 could be reversed by pre-exposure of animals to Zn. We were not able to link Zn levels to the expression of HSP70 in the tissues of the earthworm L. terrestris. Earthworms exposed to Cd for three days showed increased HSP70 expression in coelomocytes (51) as well as sea star coelomocytes six hours post-traumatic stress (52). Extracellular HSP70 has been earlier shown to regulate innate immune cell function (53). Arsenite exposure induced the expression of HSP70, which subsequently reduced the anti-inflammatory response to lipopolysaccharides in canine 030D cells (54). It would be interesting to evaluate a putative link between HSP70 and the anti-inflammatory response in earthworms in future studies.
We identified partial sequences of ATF2, ATF7 and CREB (Supplementary Table S3). These factors contributed the most to the separation of amputated earthworms (Cut and Cd_Cut) along Dim1 (Supplementary Table S4). ATF2 copy numbers revealed a significant increase in amputated earthworms after 1 and 2 weeks (Figure 5A). ATF7 gene expression increased significantly after two weeks in amputated earthworms (Control and Cd-exposed), and showed an additional induction in the Cd_Cut group after one week. CREB showed the same pattern (Figures 5B, C). A comparison of ATF2, ATF7 and CREB between Tissue 1 and Tissue 2 is presented in Supplementary Figure S6. It has recently been described in mice that ATF2 and ATF7 protect the intestinal epithelium from apoptosis during regeneration (55), whereas CREB is known to regulate cell proliferation and differentiation (30, 56, 57). Both factors revealed a strong positive correlation to sccTLR. Scc- and mccTLRs have been previously described for the earthworm species E. andrei (33), the sccTLR (EaTLR) was elevated in response to bacterial stress (32), while another study revealed that EaTLR showed a decreased expression anterior- but no changes in posterior regeneration (22). Moreover, mccTLRs were shown to be expressed in sea urchin coelomocytes after bacterial challenges (58). It is known from the literature that TLRs directly interact with ATF family members, such as ATF4 and TLR4 leading to cytokine production (59). In addition, several TLRs activate ATF2 leading to the stimulation of genes related to cellular motility (31) and that ATF2 protein is highly expressed in infiltrating murine macrophages in adipose tissue (60). Since we could show that sccTLR is increased in regenerative tissue, we suggest that sccTLR signaling recruits ATF7 and CREB in earthworms. mccTLR, in contrast to sccTLR, was hardly induced in regenerative tissue and did not show the same strong positive correlation to ATF7 and CREB, suggesting that mccTLR has a different function, probably not directly related to wound healing.
Figure 5 Absolute mRNA copy numbers of (A) activating transcription factor 2 (ATF2), (B) activating transcription factor 7 (ATF7) and (C) cAMP response element-binding protein (CREB) in L. terrestris tissue from Control and exposed [50 mg CdCl2/kg dry soil (Cd)] individuals; regenerative tissue of amputated (Cut) and exposed amputated [50 mg CdCl2/kg dry soil (Cd_Cut)] individuals at week one (w1), week two (w2), and week three (w3). Asterisks indicate significant differences between treatments and the control of the respective week [p < 0.05 (*); p < 0.01 (**); p < 0.001 (***); p < 0.0001 (****)]. Mean values ± SEM are presented.
Previously, ATF family members were hypothesized to play a role in Cd-induced MT2 gene expression (23). Since we did not observe a correlation between ATF2, ATF7 or CREB and MT2 gene activity, we suggest that these factors are not involved in MT2 gene regulation. A study on Cd-exposed E. andrei proposed for MT2 not only a role in Cd detoxification, but also a protective role during regeneration since in amputated and Cd-exposed earthworms MT2 expression was increased and regeneration was even faster in the presence of Cd (61). If wound healing in L. terrestris is accelerated in the presence of Cd, has to be addressed in future studies. MT2 expression increased in earthworms upon exposure to Cd (Figure 6B). However, a Cd-induced induction of the MT2 gene failed in Tissue 2 of amputated earthworms as shown in the Cd_Cut group (Figure 3D).
Figure 6 Absolute mRNA copy numbers of (A) metallothionein 1 (MT1) and (B) metallothionein 2 (MT2) in L. terrestris tissue from Control and exposed [50 mg CdCl2/kg dry soil (Cd)] individuals; regenerative tissue of amputated (Cut) and exposed amputated [50 mg CdCl2/kg dry soil (Cd_Cut)] individuals at week one (w1), week two (w2), and week three (w3). Asterisks indicate significant differences between treatments and the control of the respective week [p < 0.05 (*); p < 0.01 (**); p < 0.0001 (****)]. Mean values ± SEM are presented.
We identified a partial L. terrestris MT1 sequence (Supplementary Table S3). The highest induction of MT1 was observed in the Cd_Cut group after three weeks (p = 0.053) but only in Tissue 1 (Figure 3C). The Cd and Cut groups did not show significant differences from the respective control (Figure 6A). Numerous immunomodulatory functions have been suggested for MTs in humans, such as the role of MT1 as an enhancer in cells of innate immunity (62). Few functional studies on earthworm MT1 have been accomplished so far. A role for mt1 in Dendrobaena octaedra has been suggested in copper resistance but only in adapted populations (63). A combined treatment of copper exposure and freezing also lead to an upregulation of mt1 in D. octaedra (64). The involvement of L. rubellus MT1 in metal homeostasis, e.g. as a donor for essential metals like Zn, was proposed several years ago (13). Zn is required at all stages of wound healing (65) and it is essential for several functions in innate immunity (66). Zn is also known to be relevant for TLR/nuclear factor k-light-chain-enhancer of activated B cells (NF-kB) mediated inflammatory signaling and cytokine production (65). Total Zn and free Ca have been measured. Both, Zn and Ca levels increased significantly in amputated earthworms (Control and Cd-exposed), whereas Cd alone did not lead to significant changes in either Zn or Ca levels (Figures 7A, B). Interestingly, Zn levels were the highest in the combined stress group after two weeks, but dropped to control levels after three weeks. Interestingly, at the same time, MT1 levels were increased. Since we measured the total amount of Zn, we can exclude a role for MT1 in Zn homeostasis. Taken together, we propose a function of MT1 in the response to multiple stressors and a putative contribution to stress adaptation.
Figure 7 Levels of (A) total zinc and (B) free calcium in L. terrestris tissue from Control and exposed [50 mg CdCl2/kg dry soil (Cd)] individuals; regenerative tissue of amputated (Cut) and exposed amputated [50 mg CdCl2/kg dry soil (Cd_Cut)] individuals at week one (w1), week two (w2), and week three (w3). Asterisks indicate significant differences between treatments and the control of the respective week [p < 0.05 (*); p < 0.01 (**); p < 0.001 (***); p < 0.0001(****)]. Mean values ± SEM are presented.
A previous study showed that short-term exposure to Cd increased Ca levels in L. terrestris coelomocytes (42), which is known to be necessary for coelomocyte activation (67). Injuries in Drosophila embryos and zebrafish larvae showed that the first signal in the wound inflammatory response is Ca flashes leading to the release of hydrogen peroxide, which acts as an attractant for immune cells that are relocated to the site of injury within minutes. Moreover, Ca is known to regulate cell differentiation, contraction to reduce the wound size, and angiogenesis, and remains elevated during the inflammatory and proliferative phases while Ca levels decrease during the remodeling phase (68). Herein, we measured free Ca in tissue samples revealing that Ca levels peaked after two weeks and declined in week three, although levels were still significantly increased in amputated earthworms compared to controls. Considering human models, we conclude that earthworms transition from the proliferative into the remodeling phase between week two and three.
4 Conclusion
Earthworms are well prepared for environmental stress due to their soil-dwelling lifestyle, which entails different challenges ranging from physical stress, soil pathogens to xenobiotics. Mainly, these challenges do not occur as single events, but as a combination of stressors. We were able to show that the combination of stressors had the most severe effect. We also revealed spatial expression of MT2 and sccTLR along the posterior part of the earthworm. We think that the difference in gene expression between Tissue 1 and Tissue 2 is due to the distribution of coelomocytes during wound healing. In a previous study, we showed that MT2 expression is higher in coelomocytes than in tissue samples (9) and it is also known that TLRs are expressed in coelomocytes (58). Coelomocytes infiltrate regenerating tissue and therefore relocate to the wound site, which could explain the differential gene expression levels of MT2 and sccTLR in Tissue 1 compared to Tissue 2. However, the reason why we still don’t observe Cd-induced MT2 induction in Tissue 2 after three weeks, when the numbers of coelomocyte should have recovered, still needs to be clarified. We can only speculate that a combination of stressors slows down proliferation of coelomocytes, which might be reinforced by a challenged energy metabolism.
Data availability statement
The original contributions presented in the study are publicly available. This data can be found here: https://www.ncbi.nlm.nih.gov/nuccore; OR596220-OR596222, OR610157-OR610159.
Ethics statement
Ethical approval was not required for the study involving animals in accordance with the local legislation and institutional requirements because there is no ethical approval required for invertebrate species.
Author contributions
MH: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. GA: Data curation, Formal Analysis, Investigation, Methodology, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. VP: Data curation, Investigation, Methodology, Visualization, Writing – review & editing. BF: Data curation, Methodology, Writing – original draft. CP: Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The work was supported by project P33835-B from the Austrian Science Foundation (FWF).
Acknowledgments
We would like to thank Tobias Sapper for his contribution to generate gene expression data.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1272191/full#supplementary-material
References
1. Phillips HRP, Guerra CA, Bartz MLC, Briones MJI, Brown G, Crowther TW, et al. Global distribution of earthworm diversity. Science (2019) 366(6464):480–5. doi: 10.1126/SCIENCE.AAX4851/SUPPL_FILE/AAX4851-PHILLIPS-SM.PDF
2. Dales RP, Kalaç Y. Phagocytic defence by the earthworm Eisenia foetida against certain pathogenic bacteria. Comp Biochem Physiol Part A: Physiol (1992) 101(3):487–90. doi: 10.1016/0300-9629(92)90499-G
3. Kokhanyuk B, Bodó K, Sétáló G, Németh P, Engelmann P. Bacterial engulfment mechanism is strongly conserved in evolution between earthworm and human immune cells. Front Immunol (2021) 12:733541/BIBTEX. doi: 10.3389/FIMMU.2021.733541/BIBTEX
4. Six L, Smolders E. Future trends in soil cadmium concentration under current cadmium fluxes to European agricultural soils. Sci Total Environ (2014) 485–486(1):319–28. doi: 10.1016/J.SCITOTENV.2014.03.109
5. Joimel S, Cortet J, Jolivet CC, Saby NPA, Chenot ED, Branchu P, et al. Physico-chemical characteristics of topsoil for contrasted forest, agricultural, urban and industrial land uses in France. Sci Total Environ (2016) 545–546:40–7. doi: 10.1016/j.scitotenv.2015.12.035
6. Thévenod F, Fels J, Lee WK, Zarbock R. Channels, transporters and receptors for cadmium and cadmium complexes in eukaryotic cells: myths and facts. BioMetals (2019) 32(3):469–89. doi: 10.1007/S10534-019-00176-6
7. Šrut M, Drechsel V, Höckner M. Low levels of Cd induce persisting epigenetic modifications and acclimation mechanisms in the earthworm Lumbricus terrestris. PloS One (2017) 12(4):e0176047. doi: 10.1371/journal.pone.0176047
8. Mirkov I, Popov Aleksandrov A, Ninkov M, Tucovic D, Kulas J, Zeljkovic M, et al. Immunotoxicology of cadmium: Cells of the immune system as targets and effectors of cadmium toxicity. Food Chem Toxicol (2021) 149:112026. doi: 10.1016/j.fct.2021.112026
9. Aigner GP, Nenning P, Fiechtner B, Šrut M, Höckner M. DNA methylation and detoxification in the earthworm lumbricus terrestris exposed to cadmium and the DNA demethylation agent 5-aza-2′-deoxycytidine. Toxics (2022) 10(2):100. doi: 10.3390/toxics10020100
10. Spurgeon DJ, Weeks JM, Van Gestel CAM. A summary of eleven years progress in earthworm ecotoxicology: The 7th international symposium on earthworm ecology · Cardiff · Wales · 2002. Pedobiologia (2003) 47(5–6):588–606. doi: 10.1078/0031-4056-00234
11. Homa J, Stürzenbaum SR, Kolaczkowska E. Metallothionein 2 and heat shock protein 72 protect allolobophora chlorotica from cadmium but not nickel or copper exposure: body malformation and coelomocyte functioning. Arch Environ Contamination Toxicol (2016) 71(2):267–77. doi: 10.1007/S00244-016-0276-6/FIGURES/5
12. Sandbichler AM, Höckner M. Cadmium protection strategies–A hidden trade-off? Int J Mol Sci (2016) 17(1). doi: 10.3390/IJMS17010139
13. Stürzenbaum SR, Winters C, Galay M, Morgan AJ, Kille P. Metal ion trafficking in earthworms. J Biol Chem (2001) 276(36):34013–8. doi: 10.1074/jbc.M103605200
14. Stürzenbaum SR, Georgiev O, Morgan AJ, Kille P. Cadmium detoxification in earthworms: from genes to cells. Environ Sci Technol (2004) 38(23):6283–9. doi: 10.1021/ES049822C
15. Hussain N, Chatterjee SK, Maiti TK, Goswami L, Das S, Deb U, et al. Metal induced non-metallothionein protein in earthworm: A new pathway for cadmium detoxification in chloragogenous tissue. J Hazardous Mater (2021) 401:123357. doi: 10.1016/J.JHAZMAT.2020.123357
16. Özpolat BD, Bely AE. Developmental and molecular biology of annelid regeneration: a comparative review of recent studies. Curr Opin Genet Dev (2016) 40:144–53. doi: 10.1016/J.GDE.2016.07.010
17. Sahu R, Sahu SK, Nishank SS. de novo transcriptome profile of two earthworms Lampito mauritii and Drawida calebi during regeneration. Biochem Biophys Rep (2021) 27. doi: 10.1016/J.BBREP.2021.101092
18. Bhambri A, Dhaunta N, Singh Patel S, Hardikar M, Bhatt A, Srikakulam N, et al. Large scale changes in the transcriptome of Eisenia fetida during regeneration. PloS One (2018) 13(9):e0204234. doi: 10.1371/journal.pone.0204234
19. Kostyuchenko RP, Kozin VV. Comparative aspects of annelid regeneration: towards understanding the mechanisms of regeneration. Genes (2021) 12(8):1148. doi: 10.3390/genes12081148
20. Burke JM. Wound healing in Eisenia foetida (oligochaeta). 1. Histology and 3H-thymidine radioautography of the epidermis. J Exp Zool (1974) 188(1):49–63. doi: 10.1002/JEZ.1401880106
21. Amaroli A, Ferrando S, Pozzolini M, Gallus L, Parker S, Benedicenti S. The earthworm Dendrobaena veneta (Annelida): A new experimental-organism for photobiomodulation and wound healing. Eur J Histochem (2018) 62(1):51–61. doi: 10.4081/ejh.2018.2867
22. Bodó K, Kellermayer Z, László Z, Boros Á, Kokhanyuk B, Németh P, et al. Injury-induced innate immune response during segment regeneration of the earthworm, eisenia andrei. Int J Mol Sci (2021) 22(5):1–23. doi: 10.3390/IJMS22052363
23. Drechsel V, Schauer K, Šrut M, Höckner M. Regulatory plasticity of earthworm wMT-2 gene expression. Int J Mol Sci (2017) 18(6). doi: 10.3390/ijms18061113
24. Linnaeus C. (1758) Systema naturae per regna tria naturae :secundum classes, ordines, genera, species, cum characteribus, differentiis, synonymis, locis. Editio Decima, Reformata [10th Revised Edition],. Vol. 1: 824 Pp. Laurentius Salvius: Holmiae. https://doi.org/10.5962/bhl.title.542
25. Cikutovic MA, Fitzpatrick LC, Goven AJ, Venables BJ, Giggleman MA, Cooper EL. Wound healing in earthworms lumbricus terrestris: A cellular-based biomarker for assessing sublethal chemical toxicity. Bull Environ Contamination Toxicol (1999) 62:508–14. doi: 10.1007/s001289900905
26. Cooper EL, Roch P. The capacities of earthworms to heal wounds and to destroy allografts are modified by polychlorinated biphenyls (PCB). J Invertebrate Pathol (1992) 60(1):59–63. doi: 10.1016/0022-2011(92)90154-V
27. Ville P, Roch P, Cooper EL, Masson P, Narbonne JF. PCBs increase molecular-related activities (lysozyme, antibacterial, hemolysis, proteases) but inhibit macrophage-related functions (phagocytosis, wound healing) in earthworms. J Invertebrate Pathol (1995) 65(3):217–24. doi: 10.1006/JIPA.1995.1033
28. Janeway CA, Medzhitov R. Innate immune recognition. Annu Rev Immunol (2002) 20:197–216. doi: 10.1146/ANNUREV.IMMUNOL.20.083001.084359
29. Akira S, Takeda K. Toll-like receptor signalling. Nat Rev Immunol (2004) 4(7):499–511. doi: 10.1038/nri1391
30. Wen AY, Sakamoto KM, Miller LS. The role of the transcription factor CREB in immune function. J Immunol (2010) 185(11):6413–9. doi: 10.4049/jimmunol.1001829
31. Yu T, Li YJ, Bian AH, Zuo H, Zhu TW, Ji SX, et al. The regulatory role of activating transcription factor 2 in inflammation. Mediators Inflammation (2014) 2014. doi: 10.1155/2014/950472
32. Škanta F, Roubalová R, Dvořák J, Procházková P, Bilej M. Molecular cloning and expression of TLR in the Eisenia andrei earthworm. Dev Comp Immunol (2013) 41(4):694–702. doi: 10.1016/j.dci.2013.08.009
33. Prochazkova P, Roubalova R, Skanta F, Dvorak J, Pacheco NIN, Kolarik M, et al. Developmental and immune role of a novel multiple cysteine cluster TLR from eisenia andrei earthworms. Front Immunol (2019) 10:1277(JUN). doi: 10.3389/fimmu.2019.01277
34. Novais SC, Soares AMVM, De Coen W, Amorim MJB. Exposure of Enchytraeus albidus to Cd and Zn - Changes in cellular energy allocation (CEA) and linkage to transcriptional, enzymatic and reproductive effects. Chemosphere (2013) 90(3):1305–9. doi: 10.1016/j.chemosphere.2012.09.030
35. Mihaylova MM, Shaw RJ. The AMPK signalling pathway coordinates cell growth, autophagy and metabolism. Nat Cell Biol (2011) 13(9):1016–23. doi: 10.1038/ncb2329
36. Hardie DG, Ross FA, Hawley SA. AMPK: a nutrient and energy sensor that maintains energy homeostasis. Nat Rev Mol Cell Biol (2012) 13(4):251–62. doi: 10.1038/nrm3311
37. Herzig S, Shaw RJ. AMPK: guardian of metabolism and mitochondrial homeostasis. Nat Rev Mol Cell Biol (2018) 19(2):121–35. doi: 10.1038/NRM.2017.95
38. Fjøsne T. F., Stenseth E. B., Myromslien F., Rudi K. (2015) Gene expression of TLR homologues identified by genome-wide screening of the earthworm Dendrobaena veneta. Innate Immunity Vol. 21(2): 161–166. https://doi.org/10.1177/1753425914523056
39. Masuko T, Minami A, Iwasaki N, Majima T, Nishimura SI, Lee YC. Carbohydrate analysis by a phenol-sulfuric acid method in microplate format. Analytical Biochem (2005) 339(1):69–72. doi: 10.1016/J.AB.2004.12.001
40. Stürzenbaum SR, Andre J, Kille P, Morgan AJ. Earthworm genomes, genes and proteins: the (re)discovery of Darwin’s worms. Proc Biol Sci (2009) 276(1658):789–97. doi: 10.1098/RSPB.2008.1510
41. Shao Y, Wang XB, Zhang JJ, Li ML, Wu SS, Ma XY, et al. Genome and single-cell RNA-sequencing of the earthworm Eisenia andrei identifies cellular mechanisms underlying regeneration. Nat Commun (2020) 11(1):1–15. doi: 10.1038/s41467-020-16454-8
42. Höckner M, Tomanek L, Piechnik CA, Fiechtner B, Weinberger B. Cadmium-related effects on cellular immunity comprises altered metabolism in earthworm coelomocytes. Int J Mol Sci (2020) 21(2):599. doi: 10.3390/ijms21020599
43. Muangphra P, Gooneratne R. Comparative genotoxicity of cadmium and lead in earthworm coelomocytes. Appl Environ Soil Sci (2011) 2011:1–7. doi: 10.1155/2011/218929
44. Mo X, Qiao Y, Sun Z, Sun X, Li Y. Molecular toxicity of earthworms induced by cadmium contaminated soil and biomarkers screening. J Environ Sci (China) (2012) 24(8):1504–10. doi: 10.1016/S1001-0742(11)60957-1
45. Wang X, Chang L, Sun Z, Zhang Y, Yao L. Analysis of earthworm Eisenia fetida proteomes during cadmium exposure: An ecotoxicoproteomics approach. PROTEOMICS (2010) 10(24):4476–90. doi: 10.1002/PMIC.201000209
46. Granger Joly de Boissel P, Fournier M, Rodriguez-Lecompte JC, McKenna P, Kibenge F, Siah A. Functional and molecular responses of the blue mussel Mytilus edulis’ hemocytes exposed to cadmium - An in vitro model and transcriptomic approach. Fish Shellfish Immunol (2017) 67:575–85. doi: 10.1016/J.FSI.2017.06.001
47. Wang T, Yu Q, Chen J, Deng B, Qian L, Le Y. PP2A mediated AMPK inhibition promotes HSP70 expression in heat shock response. PloS One (2010) 5(10). doi: 10.1371/JOURNAL.PONE.0013096
48. Ehsan M, Yazdi T, Sadegh Amiri M, Nourbakhsh F, Rahnama M, Forouzanfar F, et al. Bio-indicators in cadmium toxicity: Role of HSP27 and HSP70. Environ Sci pollut Res (2021) 26:26359–79. doi: 10.1007/s11356-021-13687-y/Published
49. Jing J, Liu H, Chen H, Hu S, Xiao K, Ma X. Acute effect of copper and cadmium exposure on the expression of heat shock protein 70 in the Cyprinidae fish Tanichthys albonubes. Chemosphere (2013) 91(8):1113–22. doi: 10.1016/J.CHEMOSPHERE.2013.01.014
50. Wang CC, Si LF, Guo SN, Zheng JL. Negative effects of acute cadmium on stress defense, immunity, and metal homeostasis in liver of zebrafish: The protective role of environmental zinc dpre-exposure. Chemosphere (2019c) 222:91–7. doi: 10.1016/J.CHEMOSPHERE.2019.01.111
51. Homa J, Klimek M, Kruk J, Cocquerelle C, Vandenbulcke F, Plytycz B. Metal-specific effects on metallothionein gene induction and riboflavin content in coelomocytes of Allolobophora chlorotica. Ecotoxicol Environ Saf (2010) 73(8):1937–43. doi: 10.1016/j.ecoenv.2010.07.001
52. Pinsino A, Thorndyke MC, Matranga V. Coelomocytes and post-traumatic response in the common sea star Asterias rubens. Cell Stress Chaperones (2007) 12(4):331. doi: 10.1379/CSC-288.1
53. Borges TJ, Wieten L, Van Herwijnen MJC, Broere F, Van derZee R, Bonorino C, et al. The anti-inflammatory mechanisms of Hsp70. Front Immunol (2012) 3:95(MAY). doi: 10.3389/FIMMU.2012.00095
54. Lyu Q, Wawrzyniuk M, Rutten VPMG, van Eden W, Sijts AJAM, Broere F. Hsp70 and NF-kB mediated control of innate inflammatory responses in a canine macrophage cell line. Int J Mol Sci (2020) 21(18):6464. doi: 10.3390/ijms21186464
55. Meijer BJ, Giugliano FP, Baan B, van der Meer JHM, Meisner S, van Roest M, et al. ATF2 and ATF7 are critical mediators of intestinal epithelial repair. Cell Mol Gastroenterol Hepatol (2020) 10(1):23–42. doi: 10.1016/j.jcmgh.2020.01.005
56. Montminy M. Transcriptional regulation by cyclic AMP. Annu Rev Biochem (1997) 66(1):807–22. doi: 10.1146/annurev.biochem.66.1.807
57. Wang XS, Zhang S, Xu Z, Zheng SQ, Long J, Wang DS. Genome-wide identification, evolution of ATF/CREB family and their expression in Nile tilapia. Comp Biochem Physiol Part - B: Biochem Mol Biol (2019) 237:110324. doi: 10.1016/j.cbpb.2019.110324
58. Buckley KM, Rast JP. Dynamic evolution of toll-like receptor multigene families in echinoderms. Front Immunol (2012) 3:136(JUN). doi: 10.3389/fimmu.2012.00136
59. Zhang C, Bai N, Chang A, Zhang Z, Yin J, Shen W, et al. ATF4 is directly recruited by TLR4 signaling and positively regulates TLR4-trigged cytokine production in human monocytes. Cell Mol Immunol (2013) 10(1):84–94. doi: 10.1038/CMI.2012.57
60. Miyata Y., Fukuhara A., Otsuki M., Shimomura I. (2013) Expression of activating transcription factor 2 in inflammatory macrophages in obese adipose tissue. Obesity 21(4): 731–736. https://doi.org/10.1002/oby.20274
61. Rorat A, Vandenbulcke F, Gałuszka A, Klimek B, Plytycz B. Protective role of metallothionein during regeneration in Eisenia andrei exposed to cadmium. Comp Biochem Physiol Part - C: Toxicol Pharmacol (2017) 203:39–50. doi: 10.1016/j.cbpc.2017.10.002
62. Dai H, Wang L, Li L, Huang Z, Ye L. Metallothionein 1: A new spotlight on inflammatory diseases. Front Immunol (2021) 12:739918. doi: 10.3389/fimmu.2021.739918
63. Fisker KV, Holmstrup M, Sørensen JG. Variation in metallothionein gene expression is associated with adaptation to copper in the earthworm Dendrobaena octaedra. Comp Biochem Physiol - C Toxicol Pharmacol (2013) 157(2):220–6. doi: 10.1016/j.cbpc.2012.11.007
64. Fisker KV, Holmstrup M, Sørensen JG. Freezing of body fluids induces metallothionein gene expression in earthworms (Dendrobaena octaedra). Comp Biochem Physiol Part - C: Toxicol Pharmacol (2016) 179:44–8. doi: 10.1016/j.cbpc.2015.08.008
65. Lin PH, Sermersheim M, Li H, Lee PHU, Steinberg SM, Ma J. Zinc in wound healing modulation. In: Nutrients, vol. 10. MDPI AG (2018). doi: 10.3390/nu10010016
66. Maares M, Haase H. Zinc and immunity: An essential interrelation. Arch Biochem Biophys (2016) 611:58–65. doi: 10.1016/j.abb.2016.03.022
67. Opper B, Németh P, Engelmann P. Calcium is required for coelomocyte activation in earthworms. Mol Immunol (2010) 47(11–12):2047–56. doi: 10.1016/J.MOLIMM.2010.04.008
Keywords: Annelida, stress, injury, amputation, innate immunity, immune cells, metallothionein, toll-like receptors
Citation: Aigner GP, Peer V, Fiechtner B, Piechnik CA and Höckner M (2023) Wound healing and Cadmium detoxification in the earthworm Lumbricus terrestris – a potential case for coelomocytes? Front. Immunol. 14:1272191. doi: 10.3389/fimmu.2023.1272191
Received: 03 August 2023; Accepted: 16 November 2023;
Published: 05 December 2023.
Edited by:
Satya Sundar Bhattacharya, Tezpur University, IndiaCopyright © 2023 Aigner, Peer, Fiechtner, Piechnik and Höckner. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Martina Höckner, martina.hoeckner@uibk.ac.at