Next Article in Journal
Presence of Arsenic in Potential Sources of Drinking Water Supply Located in a Mineralized and Mined Area of the Sierra Madre Oriental in Mexico
Previous Article in Journal
Variable Fitness Response of Two Rotifer Species Exposed to Microplastics Particles: The Role of Food Quantity and Quality
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxics
Volume 9
Issue 11
10.3390/toxics9110306
toxics-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Cadmium Sulfate on the Brown Garden Snail Cornu aspersum: Implications for DNA Methylation

by
Marius Georgescu
1,2,†,
George Andrei Drăghici
3,4,†,
Eliza-Florentina Oancea
5,
Cristina Adriana Dehelean
3,4,*,
Codruţa Şoica
3,4,
Nicolae-Valentin Vlăduţ
6 and
Dragoș Vasile Nica
3,6,*
1
Department of Functional Sciences, Physiology Discipline, “Victor Babeș” University of Medicine and Pharmacy Timișoara, Eftimie Murgu Sq. no. 2, 300041 Timișoara, Romania
2
Center of Immuno-Physiology and Biotechnologies (CIFBIOTEH), “Victor Babeș” University of Medicine and Pharmacy Timișoara, Eftimie Murgu Sq. no. 2, 300041 Timișoara, Romania
3
Faculty of Pharmacy, “Victor Babeș” University of Medicine and Pharmacy Timișoara, Eftimie Murgu Square No. 2, 300041 Timișoara, Romania
4
Research Center for Pharmaco-Toxicological Evaluations, Faculty of Pharmacy, “Victor Babes” University of Medicine and Pharmacy Timișoara, Eftimie Murgu Square no. 2, 300041 Timișoara, Romania
5
Faculty of Pharmacy, “Vasile Goldiș” Western University of Arad, B-dul Revoluției no. 94, 310025 Arad, Romania
6
INMA Bucharest, Bulevardul Ion Ionescu de la Brad 6, 077190 București, Romania
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Toxics 2021, 9(11), 306; https://doi.org/10.3390/toxics9110306
Submission received: 25 September 2021 / Revised: 1 November 2021 / Accepted: 8 November 2021 / Published: 15 November 2021
(This article belongs to the Section Ecotoxicology)
Browse Figures
Versions Notes

Abstract

:
An extensive literature exists regarding the cellular, physiological, and genetic effects of cadmium (Cd)—A highly toxic, but commonly used trace metal in modern industry. However, limited data are available on its epigenetic effects, especially for terrestrial sentinel invertebrates. We determined Cd retention, total DNA methylation, and the methylation status of 5′ end of the Cd-MT gene in the hepatopancreas of the brown garden snail, Cornu aspersum, fed Cd sulfate for four weeks. Bodyweight changes and survival were also measured. Hepatopancreas cadmium increased in a dose-dependent manner from the third-lowest dose onward, with very large amounts being found for the highest treatment group. However, no mortalities occurred, irrespective of dietary Cd dose. We identified significant genome-wide hypermethylation in specimens given the highest dose, which overlapped with a significant bodyweight decrease. The Cd-MT gene showed an unmethylated 5′ end of the Cd-MT gene and this status was not affected by cadmium exposure. Hepatopancreas DNA methylation is as sensitive as bodyweight to non-lethal concentrations of dietary Cd given as cadmium sulfate but less responsive than tissue accumulation. Such an exposure event, by contrast, does not affect the methylation status of the Cd-MT gene 5′ end.

1. Introduction

Cadmium (Cd) is a non-essential trace metal classified as a global priority pollutant due to its environmental persistence, bio-accumulative nature, and high toxicity at low ambient concentrations [1,2,3]. Decades of toxicological research have yielded a compelling body of evidence for the cellular, physiological, and genetic implications of exposure to Cd and its compounds [2,4,5]. Such events are also known to affect epigenetic mechanisms shaping gene expression, including DNA methylation [6,7]. 5-Methylcytosine (5mC) is the best-characterized and most common form of DNA methylation, being linked to stable, long-term transcriptional silencing [8]. This evolutionary conserved epigenetic mark [9] occurs via covalent addition of a methyl (CH3) group to the 5′-carbon position of cytosine pyrimidine ring primarily found in cytosine-guanine (CpG) dinucleotides [10]. Global and gene-specific changes in 5mC content associated with Cd exposure have been described in most vertebrate classes [11,12,13,14,15]. However, limited information exists for terrestrial invertebrates routinely used as sentinel organisms, such as annelids [16] or land snails [7]. The latter species are pivotal for terrestrial ecosystem functioning, serving as key players in calcium cycling and acting on trophic level, both as herbivorous and detritivore organisms [17,18].
CpG islands are often located at the 5′-end of genes, with promoters and transcription start sites embedded within these genomic regions of high CpG dinucleotide densities [19,20]. The methylation status of these gene-regulatory elements can play critical roles in modulating vertebrate gene expression, especially for inducible genes responsive to various biotic and abiotic environmental factors [20]. Cadmium exposure has been associated with altered patterns of DNA methylation at the 5′-end of vertebrate genes [21,22]. For most invertebrates, by contrast, there is only limited information regarding the occurrence of cytosine methylation at the 5′-region of genes and the effects of Cd on DNA methylation [20,23,24]. Of particular importance in the context of cadmium and terrestrial gastropods is to identify the methylation status in the 5′-region of the Cd-selective metallothionein (Cd-MT) gene and to determine whether such an exposure event can affect it. Unlike constitutively expressed genes with housekeeping functions, the Cd-MT gene is inducible and pivotal for cadmium detoxification in many species of land snails, including the brown garden snail, Cornu aspersum (Müller 1774)—an appropriate invertebrate model for investigating biological responses to cadmium [7,25,26,27].
To the best of our knowledge, only one study to date has investigated the effect of Cd on gastropod methylome. When given as cadmium chloride (CdCl2), dietary Cd affected the global 5mC content of genomic DNA in the hepatopancreas of mature C. aspersum snails [7]. Different anions/cations could exert different methylomic effects [28], but the impact of other cadmium compounds on epigenetics in this specific snail and other gastropod species is still not clear. These data are important for broadening our understanding of Cd ecological risk, especially for compounds commonly used in modern industry, but with less studied epigenotoxic profiles. Such an example is cadmium sulfate (CdSO4), a key compound in the manufacture of fluorescent screens, vacuum tubes, transistors, and photovoltaic and solar cells; as well as a pigment in paints and inks [3].
Experimental data also showed that the Cd-MT gene transcription is initiated after Cd exposure and shows a persistent induction rate [26]. However, the methylation status of the 5′ end of this gene or the effect of Cd on this gene segment has not been yet examined. At least two important questions derive from the foregoing considerations: (1) Whether and to which extent does cadmium sulfate exposure modify global 5mC levels of C. aspersus hepatopancreas cells? (2) What is the relevance of the methylation status of CG pairs at the 5′ region of the Cd-MT gene as an endpoint of Cd exposure. To test these hypotheses, mature specimens of C. aspersum were maintained under laboratory conditions for four weeks and given Cd-spiked diets at doses covering a wide concentration range, including the 0.2 milligrams per kilogram fresh weight benchmark level allowed in vegetal foods regularly eaten by terrestrial gastropods, such as fruits and leafy vegetables [29,30]; and those for Cd toxicity on key cellular endpoints in these invertebrates [7,31]. Genome-wide DNA methylation, methylation of the 5′-UTR of the Cd-MT gene, and Cd levels in the hepatopancreas of Cornu aspersum were determined. Bodyweight changes and snail death rates were used as a sublethal endpoint and a lethal endpoint, respectively. Our findings are important because they not only expand previous knowledge about the effect of Cd on an ecologically relevant study system but also fill a gap concerning stress-dependent changes in methylation patterns in gastropod species.

2. Materials and Methods

2.1. Landsnail Origin, Maintenance and Exposure

Cornu aspersum is one of the most widely used gastropod models in ecotoxicological studies because it can be easily maintained under laboratory conditions and is commercially available [7]. A total of 180 adult specimens (with a reinforced lip at shell aperture) were obtained from the “Mokry Dwór” snail farm (Krzymów, Wielkopolska, Poland). The snails were kept in groups of 10 in 20 L aerated polypropylene containers (50 × 20 × 20 cm) with perforated side walls at 20 °C in a 12/12 light/dark cycle for four weeks prior to the beginning of the experiments. The container bottom was covered with a layer of ash-free filter paper (50 × 20 cm), which was wetted twice daily with double distilled water. Exposure experiments were conducted under controlled laboratory conditions similar to those used in our previous investigations, using Cd-spiked, agar-based diets at concentrations between 0 and 100 mg per kg dry weight (mg/kg d. wt) [1,7]. Hepatopancreas samples were taken at the end of the exposure period to assess the effect of Cd on genome-wide 5mC levels and the methylation status of the Cd-MT gene. All experiments were performed at the Laboratory Animal Facility from the “Victor Babes” University of Medicine and Pharmacy from Timisoara.
At the start of the exposure experiments (September 2018), the most homogeneous 120 specimens were sorted based on their size, as determined via measuring the shell height (with a digital caliper) and the snail weight (with an analytical balance). The mean shell size of selected gastropods was 2.95 ± 0.34 cm in length, whereas the mean weight was 8.79 ± 1.85 g. The snails were given ad libidum Cd-spiked diets for four weeks. The eight nominal Cd treatments used to prepare the agar-based diets were 0, 0.05, 0.1, 0.2, 1, 10 and 100 mg/L (abbreviated as Cd0, Cd0.05, Cd0.1, Cd0.2, Cd1, Cd10, and Cd100). In contrast to our previous studies [1,4,7], cadmium sulfate (CdSO4 ≥ 99.99% trace metal basis(Sigma-Aldrich, St. Louis, MO, USA) was used here as a source of cadmium instead of cadmium chloride.
Exposures were run in two independent experiments per treatment group (Replica I and Replica II), with each replica consisting of 10 snails. Body weight was measured at the start and the end of the 4-week experimental period for four individuals per replicate jar. Four specimens for each experimental group (i.e., two specimens per replicate) were then dissected. This relatively modest sample size is not unusual for global DNA methylation studies due to cost constraints. In addition, our previous work showed that even smaller sample sizes, that is three specimens for each treatment group, can be used to detect Cd-related changes in total 5mC content in DNA of the hepatopancreas of adult C. aspersum [7]. The soft body was removed from the shell with a hemostat to collect hepatopancreas samples. These samples were used for performing chemical analyses and extracting DNA. Weekly food samples were randomly taken to check the feed cadmium concentration.

2.2. Snail Fitness Parameters

Snail fitness was assessed by determining the weight gain and survival rates. Weight gain was calculated as the difference between the snail weight at the start and the end of the exposure period. All snails were weighed with an analytical balance to the nearest 0.01 mg.

2.3. Chemical Analyses

Hepatopancreas cadmium concentrations were assessed as previously described [1,4,7]. In short, the frozen samples (n = 4 for each treatment) were thawed and oven-dried at 105 °C for 24 h. After being weighed and calcinated, these samples were dissolved by wet digestion with 65% nitric acid (Merck KGaAe) and heated on a hot plate to dryness. The ash was dissolved in 20 mL of 0.5 N HNO3, filtered through ash-free filter paper and the filtrate volume was brought to 30 mL with 10 mL HNO3 0.5 N. A similar approach was used for determining Cd levels in food samples, with triplicate samples per each treatment being analyzed at one, two, and three weeks.
Measurements were performed via flame (air-acetylene) atomic absorption spectrometry using a VARIAN AA240FS (Agilent, Santa Clara, CA, USA) fitted with a Cd-specific hollow-cathode lamp as a source of radiation. The results were expressed as milligram per kilogram dry weight (mg/kg d. wt). Standard solutions (1000 mg/L) of Cd-ICP multi-element standard solution IV CertiPUR were purchased from Merck. Only pure double distilled water was used to prepare the reagents and standard solutions. The cadmium quantification limit was 0.01 mg/kg d. wt. It was determined via the calibration curve used to validate the accuracy of the method by measuring the standard reference material obtained from the China National Analysis Center for Iron and Steel (NCS Certified Reference Material85105a). The average percent recovery was 98%, with the corresponding variation coefficients being less than 6%. The spectrometer was recalibrated every 10 samples, and blank samples were run every 15 samples. Blank samples were prepared in a similar manner to the aforementioned samples, using pure double distilled water instead of standard solutions.

2.4. DNA Extraction and 5mC Quantification Analysis

Genomic DNA was isolated from hepatopancreas samples using the GenElute™ Mammalian Genomic DNA Miniprep Kit (Merck KGaA, Darmstadt, Germany) according to the manufacturer’s directions. The DNA purity was checked by UV absorption (260/280 nm, DS-11 Spectrophotometer, DeNovix, Wilmington, NC, USA) [7]. Total 5mC DNA level was quantified using the MethylFlash Global DNA Methylation (5mC) ELISA Easy Kit (EpiGentek, Farmingdale, NY, USA). This kit provides a fast and precise method to determine global DNA methylation changes, with a detection limit of 5 nanograms (ng) of fully methylated DNA. In addition, the use of such ELISA-based assays for serial measurements of 5mC can serve as a cost-effective and reliable alternative to commonly used methods for whole-genome DNA methylation quantitation, such as chromatography, radioactive filter-binding, or bisulfite sequencing-based methods [32,33]. Optical density (OD) was measured at a wavelength of 450 nanometers (nm) with a Stat Fax 4200 microplate reader (Awareness Technology, Palm City, FL, USA). For absolute 5mC quantifications, a standard curve was created by graphing the corresponding ODs versus the log10-transformed concentrations of the positive controls. The coefficients of determination (R2) for ELISA calibration curves were above 0.97; an example is shown in Figure 1.

2.5. Methyl Specific Polymerase Chain Reaction (MSP)

Until now, only a 1128-base-pair sequence of the Cd-MT gene promoter of Cornu aspersum is known. This gene segment includes the metal responsive element binding sites (MREs), which are central to metal-dependent upregulation of MT gene expression in vertebrate species [34]. However, it is not known how far the promoter region may extend upstream beyond the known sequence. The cDNA sequence of the C. aspersum hepatopancreas encoding the entire Cd-MT protein (GenBank: EF152281.1) was hence used to manually design suitable primers for methylation analysis by bisulfite sequencing. The 5′-UTR sequence of the Cd-MT cDNA contains several CpG sites upstream of the coding ATG. Initially, one forward primer (CdMTbs-F1/ggatttatYGtaggatattaattaagg, promoter) and two reverse primers (CdMTbs-R1/ctaaaaataaaaccaaataccaatcctac, position 361; CdMTbs-R2/Ccttaccacacttacaaccatc, position 300) were designed to determine optimal bisulfite PCR conditions. Details regarding the CpG and GpC sites analyzed are given in Figure S1 (Supplementary Material). Bisulfite conversion of genomic DNA isolated from hepatopancreas was performed using the ‘high molarity’ protocol as previously described [35]. Typically, PCRs were run in a 40 µL reaction volume, with ~0.25–0.5 µg of bisulfite converted DNA using the HotStarTaq Master Mix reagent (Qiagen, Hilden, Germany), following the specifications of the supplier’s manual. A 10 µL aliquot of the PCR product was run on a 1.5% TAE (Tris-acetate-EDTA)-agarose gel to check for specific amplification. The amplicons were visualized via ethidium bromide staining following agarose gel electrophoresis (Merck, Darmstadt, Germany)

2.6. Statistical Analysis

Genome-wide DNA methylation in the hepatopancreas of brown garden snails, C. aspersum, were analyzed using a one-way ANOVA. Before conducting ANOVA, the metal and DNA methylation data sets were log10-transformed and then checked for normality and heteroscedasticity using Anderson-Darling tests and Bartlett’s tests, respectively. If the ANOVA results were significant, posthoc analysis was performed against controls using Dunnet’s tests. The same statistical approach was used for analyzing the effect of dietary cadmium on absolute values for weight gain/loss. This parameter was expressed as the difference between the measured values of the bodyweight at the end and the start of the experiment. Finally, a Chi2 test was conducted on a 7 × 2 contingency table to compare snail survival across different treatment groups, with tests based on 2 × 2 contingency tables being applied against controls in case of significant differences. All statistical analyses were performed using Statistica version 8 (StatSoft Inc., Tulsa, OK, USA). Statistical significance was defined at p less than 0.05.

3. Results

The mean levels of Cd in fodder samples were (i) for the Cd0 treatment, below the detection limit (0.01 mg/kg d. wt); (ii) for the Cd0.05 treatment, 0.05 ± 0.02 mg/kg d. wt; (iii) for the Cd0.1 treatment, 0.11 ± 0.04 mg/kg d. wt; (iv) for the Cd0.2 treatment, 0.19 ± 0.05 mg/kg d. wt; (v) for the Cd1 treatment, 1.07 ± 0.23 mg/kg d. wt; (vi) for the Cd10 treatment,10.79 ± 4.55 mg/kg d. wt; (vii) for the Cd100 treatment, 98.87 ± 8.99 mg/kg d. wt. The consumption rates were not determined, but 3 weeks after the start of the experiment most specimens from the highest treatment group tended to avoid food and enter into aestivation.

3.1. Hepatopancreas Cadmium

Average concentrations of hepatopancreas cadmium are shown in Figure 2a. The measured values varied between 5.31 ± 2.59 mg/kg d. wt in controls and 419.70 ± 84.87 mg/kg d. wt in the highest treatment group (Cd100) (Figure 2a). The log10-transformed data sets were normally distributed (Anderson-Darling tests, p ≥ 0.472) and showed similar variances (Bartlett’s test, p = 0.646). Dietary dose had a significant influence on Cd retention in the hepatopancreas (ANOVA, (F(6, 27) = 87.89, p < 0.001). Post hoc testing using the Dunnet’s approach revealed that the measured values increased significantly (Figure 2a starting from the third-lowest dose (Cd0.2) onward.

3.2. Fitness Parameters

3.2.1. Body Weight Gain

The mean values of body weight gain (loss) are given in Figure 2b. The average weight at the start of the experiment was between 8.20g in the snails exposed to the third-highest dose (Cd1) and 9.82g in the specimens from the highest treatment group (Cd100). At the end of the experiment, the weight gain/loss ranged from -0.90g in gastropods given the highest Cd dose (Cd100) up to 3.23g in controls. The data sets related to the later variable were normally distributed (Anderson-Darling test, p ≥ 0.209) and homoscedastic (Bartlett’s test, p = 0.793). We identified a significant effect of dietary Cd uptake on weight gain (ANOVA, (F(6, 55) = 12.87, p < 0.001). Post hoc analyses using the Dunnet’s procedure revealed a similar weight gain compared to controls for most treatment groups (Figure 2b). The only exception was noted for the gastropods fed the highest Cd dose (Cd100), which showed a significant 10% reduction in body weight.

3.2.2. Snail Survival

No mortalities were observed both for Cd-fed individuals and controls. Statistical analysis revealed no significant difference in survival among different groups (Chi2 test, p = 1). Therefore, the survival of mature snails, C. aspersum, appears not to be affected by continuous four-week exposure of up to 100 mg/kg d. wt Cd in the food.

3.3. Genome-Wide DNA Methylation in the Hepatopancreas

Mean 5mC levels are depicted in Figure 3. The measured values ranged from 0.29% in the snails exposed to the lowest Cd dose (Cd0.05) up to 0.99% in the specimens from the highest treatment group (Cd100), with an overall weighted mean for all groups of 0.40% (Figure 3). These data sets met the assumptions of normality (Anderson-Darling tests, p ≥ 0.189) and homogeneity of variance (Bartlett’s tests, p = 0.912). Cadmium exposure was associated with significant changes in DNA methylation levels (ANOVA, (F(6, 27) = 5.89, p = 0.001). Post hoc testing using the Dunnet’s approach revealed significant hypermethylation for gastropods consuming the highest Cd dose (Cd100), but no effect was identified for the other treatment groups (Figure 3).

3.4. Methylation Status of CG Pairs at the 5′ Region of the Cd-MT Gene

In the case of controls, all cytosines from the seven analyzed CpG sites were converted to thymines after bisulfite treatment and PCR amplification, irrespective of the primer pairs used. Identical results were obtained for specimens given different doses of cadmium in the diet for 28 days. These data demonstrate that the 7 CpG sites at the 5′-UTR of the Cd-MT mRNA are unmethylated. The 5′ region close to the presumed Cd-MT transcription start site is not constitutively methylated and the methylation status does not appear to be influenced by cadmium exposure.

4. Discussion

We analyzed the effect of environmentally relevant concentrations of dietary cadmium (given as cadmium sulfate) on terrestrial ecosystems using the brown garden snail, C. aspersum, as an invertebrate (eco) toxicological model. The accumulation of this non-essential trace metal into the hepatopancreas of mature specimens started to increase significantly from doses as low as 0.2 mg/kg d. wt. This trend is comparable with those reported in studies using specimens of similar age, but cadmium chloride as a source of dietary Cd [1,7,36]. Hence, it appears that chloride and sulfate ions may affect in a similar manner the dynamics of cadmium retention in the hepatopancreas of mature snails, C. aspersum.
Cadmium exerted a significant effect on body weight only in adult individuals from the highest treatment group. For comparable exposure doses and times, several studies with younger specimens revealed a significant dose-dependent decrease or higher reductions of body weight (>35%) [36,37,38]. These divergent results could be attributed to an age-dependent decrease in metabolic rate [39] leading to higher gastrointestinal adsorption of Cd, and subsequently its elevated retention in the hepatopancreas [2].
All snails survived throughout the four-week experimental period, irrespective of the Cd dose administered. These results are in agreement with data from previous investigations [1,40,41] and attest to the exceptional ability of mature gastropods, C. aspersum, to tolerate substantial amounts of dietary cadmium without notable lethal effects. Low death rates (0% to 6.7%) were also identified in younger specimens given Cd-enriched diets [36,37,38,42]. Overall, these data do not render mortality or bodyweight of adult gastropods, C. aspersum, as ecologically relevant endpoints for assessing the hazardous effects of environmental cadmium over a one-month period.
The range of genome-wide DNA methylation levels observed in the present study is close to those determined using a similar ELISA-based approach in mature brown garden snails exposed to dietary cadmium as cadmium chloride [7], or in adult specimens of acute bladder snail, Physella acuta (Draparnaud, 1805), and southern creeper, Zeacumantus subcarinatus (G. B. Sowerby II, 1855), that is between 0.2% and 0.9% [43,44]. Only the gastropods fed the greatest Cd dose as cadmium sulfate exhibited marked differences compared to controls, revealing significant hypermethylation of hepatopancreas DNA at 28 days (week 4). This is in contrast with methylation data available for cadmium chloride [7]. Thus, individuals reared under identical exposure conditions (snail age, exposure duration, rearing conditions) showed global cytosine hypermethylation at doses of 10 mg/kg d. wt Cd and higher [7]. Hence, the effect of Cd on 5mC levels in hepatopancreas DNA seems to be influenced by the type of inorganic anion bound to Cd. Indeed, there is evidence that different inorganic salts of the same metal can have different effects on genomic DNA methylation levels [28]. In the hepatopancreas DNA of C. aspersum, cytosine methylation seems to be constitutively absent at the 5′-end of the Cd-MT gene. The 5′-UTR of the Cd-MT mRNA is likely located close to the promoter and our data suggest that the Cd-MT gene in this snail species is not regulated by this type of epigenetic mechanism. We also found that cadmium uptake via food had no effect on the methylation status of the Cd-MT gene 5′-region. Consistent with our findings, Drechsel et al. observed no methylation in CG pairs of the Cd-inducible metallothionein gene wMT-2 promoter from the common earthworm, Lumbricus terrestris (Linnaeus, 1758), for both Cd-exposed and control individuals [45]. However, we cannot exclude a potential effect of Cd on the DNA methylation status of promoter or introns (non-coding regions) without examining in detail these regions of the Cd-MT gene. With this in mind, the methylation status of MREs and adjacent areas in the Cd-MT promoter region warrant future investigations, given the importance of these gene segments as binding sites for metal-response element-binding transcription factors (MTFs) [34]. This topic is of interest not only for C. aspersum, but also for other commonly used gastropod models, such as the Roman snail, Helix pomatia (Linnaeus, 1758). Both species share MRE binding sites, including CG pairs in their Cd-MT gene promoter region, but their Cd-MT genes show distinct transcription patterns. These variable expression patterns appear to be partly related to structural differences in their respective Cd-MT gene promoter regions, like the variable number of MRE binding sites [34]. Their MTF-2 nucleotide sequence contains CG pairs [34], but their methylation status is yet to be determined. Another topic to be addressed in future studies is the impact of cadmium on methylation of the gene body due to its important role in modulating stress and environmental responses in many invertebrates [19], including mollusks [46].
Based on the current results, the global 5mC content of genomic DNA in hepatopancreas tissue appears to be a less sensitive biomarker of cadmium exposure compared to hepatopancreas cadmium level. However, this endpoint was as sensitive as body weight to non-lethal concentrations of dietary cadmium given as cadmium sulfate. Moreover, significant hypermethylation of hepatopancreas DNA was associated with a significant decrease in body weight and snail tendency to enter into aestivation, with all these events occurring only for the highest treatment group. In the sea cucumber Apostichopus japonicus, the aestivation-induced intestinal hypometabolism was found to be driven by transcriptional suppression of metabolic pathways mediated via DNA hypermethylation [47]. In the near relative species Helix pomatia, the Cd-MT gene expression is not only upregulated following Cd exposure but can also vary owing to desiccation, where gastropods switch into aestivation [48,49]. The changes observed here in total 5mC levels in hepatopancreas tissue may hence be, at least partly, connected to Cd-induced lethargy, and not necessarily indicate a direct effect of cadmium on the hepatopancreas DNA methylation. More precisely, Cd exposure is perhaps not the main cause for hepatopancreas DNA hypermethylation, but rather aestivation due to Cd exposure.
The main limitation of this study is the fact that we did not investigate the effect of Cd on DNA methylation profiles and/or expression of genes relevant to such an exposure event, e.g., the key genes of the DNA methylation machinery (related to DNA methyltransferases, Ten Eleven Translocation enzymes, or Methyl-CpG-binding domain proteins). By using RT-PCR and primers specific to the mammalian DNMT genes, however, we have recently found pertinent evidence for the presence of the DNMT1 gene, but not the DNMT3 genes in the C. aspersum genome [50]. This gene is central for the maintenance of DNA methylation in eukaryotes [20]. The candidate gene identified by us appears to be functionally expressed and sensitive to Cd exposure. Unfortunately, the genome of C. aspersum is yet to be sequenced and only a part of its transcriptome was characterized [51]. In fact, such high-resolution analysis is currently prohibited not only for this species but also for other terrestrial mollusks due to the limited number of reference genomes available [51] and little information existing about their methylome [9]. Moreover, the snail hepatopancreas contains different cell types [31], which may exhibit different DNA methylation levels and profiles. Therefore, future studies should aim to generate genome-wide and gene-specific, high-resolution DNA methylation profiling and conduct functional analysis of 5mC and other DNA methylation variants on individuals and different species of land snails, with emphasis on cell type, environmental status, and ontogenetic stage. Such fundamental knowledge is a sine qua non prerequisite for fully determining the impact of Cd on the hepatopancreas DNA methylome and deciphering the mechanisms underlying its effect.

5. Conclusions

To summarize, the results show that genome-wide hepatopancreas DNA methylation is a less sensitive biomarker of Cd exposure compared to tissue accumulation, but as responsive as body weight for non-lethal concentrations of dietary cadmium given as cadmium sulfate. In addition, our findings provide an indication that cytosine methylation could be constitutively lacking at 5′-end of the Cd-MT gene, and this status is independent of the effect of Cd exposure. Overall, the data obtained fill a gap concerning stress-dependent changes in methylation patterns in gastropod species and improve the current state of knowledge linking environmental pollutants and gastropod methylome.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/toxics9110306/s1, Figure S1: CpG and GpC sites (double-stranded) in the the 5′-UTR of the Cd-MT mRNA in C. aspersum.

Author Contributions

Conceptualization, M.G. and D.V.N.; methodology, E.-F.O., G.A.D. and C.Ş.; software, C.Ş. and D.V.N.; validation, M.G., C.A.D. and N.-V.V.; formal analysis, E.-F.O. and G.A.D.; investigation, M.G., E.-F.O., G.A.D. and D.V.N.; resources, C.A.D., C.Ş., N.-V.V. and D.V.N.; data curation, M.G. and G.A.D.; writing—original draft preparation, E.-F.O. and G.A.D.; writing—review and editing, M.G., C.Ş. and D.V.N.; visualization, M.G. and G.A.D.; supervision, D.V.N.; project administration, C.A.D., N.-V.V. and D.V.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nica, D.; Popescu, C.; Draghici, G.; Privistirescu, I.; Suciu, M.; Stöger, R. Effect of cadmium on cytosine hydroxymethylation in gastropod hepatopancreas. Environ. Sci. Pollut. Res. 2017, 24, 15187–15195. [Google Scholar] [CrossRef] [PubMed]
  2. ATSDR. Toxicological Profile for Cadmium. Agency for Toxic Substances and Disease Registry. Division of Toxicology and Environmental Medicine/Applied Toxicology Branch, Atlanta, Georgia. 2012. Available online: http://www.atsdr.cdc.gov/toxprofiles/tp5.pdf (accessed on 27 July 2021).
  3. Mikhailenko, A.V.; Ruban, D.A.; Ermolaev, V.A.; van Loon, A.J. Cadmium pollution in the tourism environment: A literature review. Geosciences 2020, 10, 242. [Google Scholar] [CrossRef]
  4. Nica, D.V.; Draghici, G.A.; Andrica, F.M.; Popescu, S.; Coricovac, D.E.; Dehelean, C.A.; Gergen, I.I.; Kovatsi, L.; Coleman, M.D.; Tsatsakis, A. Short-term effects of very low dose cadmium feeding on copper, manganese and iron homeostasis: A gastropod perspective. Environ. Toxicol. Pharmacol. 2019, 65, 9–13. [Google Scholar] [CrossRef] [Green Version]
  5. Genchi, G.; Sinicropi, M.S.; Lauria, G.; Carocci, A.; Catalano, A. The effects of cadmium toxicity. Int. J. Environ. Res. Public Health 2020, 17, 3782. [Google Scholar] [CrossRef]
  6. Wang, B.; Li, Y.; Shao, C.; Tan, Y.; Cai, L. Cadmium and its epigenetic effects. Curr. Med. Chem. 2012, 19, 2611–2620. [Google Scholar] [CrossRef]
  7. Nica, D.V.; Popescu, C.; Draghici, G.A.; Andrica, F.M.; Privistirescu, I.A.; Gergen, I.I.; Stöger, R. High-level dietary cadmium exposure is associated with global DNA hypermethylation in the gastropod hepatopancreas. PLoS ONE 2017, 12, e0184221. [Google Scholar] [CrossRef] [Green Version]
  8. Jones, P.A. Functions of DNA methylation: Islands, start sites, gene bodies and beyond. Nat. Rev. Genet. 2012, 13, 484–492. [Google Scholar] [CrossRef]
  9. Fallet, M.; Luquet, E.; David, P.; Cosseau, C. Epigenetic inheritance and intergenerational effects in mollusks. Gene 2020, 729, 144–166. [Google Scholar] [CrossRef] [Green Version]
  10. Liu, J.; Hu, H.; Panserat, S.; Marandel, L. Evolutionary history of DNA methylation related genes in chordates: New insights from multiple whole genome duplications. Sci. Rep. 2020, 10, 1–14. [Google Scholar] [CrossRef] [Green Version]
  11. Takiguchi, M.; Achanzar, W.E.; Qu, W.; Li, G.; Waalkes, M.P. Effects of cadmium on DNA-(Cytosine-5) methyltransferase activity and DNA methylation status during cadmium-induced cellular transformation. Exp. Cell. Res. 2003, 286, 350–365. [Google Scholar] [CrossRef]
  12. Martinez-Zamudio, R.; Ha, H.C. Environmental epigenetics in metal exposure. Epigenetics 2011, 6, 820–827. [Google Scholar] [CrossRef] [Green Version]
  13. Pierron, F.; Baillon, L.; Sow, M.; Gotreau, S.; Gonzalez, P. Effect of low-dose cadmium exposure on DNA methylation in the endangered European eel. Environ. Sci. Technol. 2013, 48, 797–803. [Google Scholar] [CrossRef]
  14. Head, J.A.; Dolinoy, D.C.; Basu, N. Epigenetics for ecotoxicologists. Environ. Toxicol. Chem. 2012, 31, 221–227. [Google Scholar] [CrossRef]
  15. Simoniello, P.; Trinchella, F.; Filosa, S.; Scudiero, R.; Magnani, D.; Theil, T.; Motta, C.M. Cadmium contaminated soil affects retinogenesis in lizard embryos. J. Exp. Zool. A Ecol. Genet. Physiol. 2014, 321, 207–219. [Google Scholar] [CrossRef] [PubMed]
  16. Š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, e0176047. [Google Scholar] [CrossRef] [Green Version]
  17. Nica, D.V.; Bordean, D.M.; Borozan, A.B.; Gergen, I.; Bura, M.; Banatean-Dunea, I. Use of land snails (Pulmonata) for monitoring copper pollution in terrestrial ecosystems. Rev. Environ. Contam. Toxicol. 2013, 95–137. [Google Scholar]
  18. Kerney, M.P.; Cameron, R.A.D. A Field Guide to the Land Snails of Britain and Northwestern Europe; William Collins Sons & Co., Ltd.: London, UK, 1979; pp. 37–98. [Google Scholar]
  19. Keller, T.E.; Han, P.; Soojin, V.Y. Evolutionary transition of promoter and gene body DNA methylation across invertebrate-vertebrate boundary. Mol. Biol. Evol. 2015, 33, 1019–1028. [Google Scholar] [CrossRef]
  20. Zhang, S.C.; Wang, M.Y.; Feng, J.R.; Chang, Y.; Ji, S.R.; Wu, Y. Reversible promoter methylation determines fluctuating expression of acute phase proteins. eLife 2020, 9, e51317. [Google Scholar] [CrossRef]
  21. Arita, A.; Costa, M. Epigenetics in metal carcinogenesis: Nickel, arsenic, chromium and cadmium. Metallomics 2009, 1, 222–228. [Google Scholar] [CrossRef] [Green Version]
  22. Sanders, A.; Smeester, L.; Rojas, D.; DeBussycher, T.; Wu, M.; Wright, F.; Zhou, Y.H.; Laine, J.; Rager, J.; Swamy, G.; et al. Cadmium exposure and the epigenome: Exposure-associated patterns of DNA methylation in leukocytes from mother-baby pairs. Epigenetics 2014, 9, 212–221. [Google Scholar] [CrossRef]
  23. Riviere, G.; Wu, G.C.; Fellous, A.; Goux, D.; Sourdaine, P.; Favrel, P. DNA methylation is crucial for the early development in the oyster C. gigas. Mar. Biotechnol. 2013, 15, 739–753. [Google Scholar] [CrossRef]
  24. Rajasethupathy, P.; Antonov, I.; Sheridan, R.; Frey, S.; Sander, C.; Tuschl, T.; Kandel, E.R. A role for neuronal piRNAs in the epigenetic control of memory-related synaptic plasticity. Cell 2012, 149, 693–707. [Google Scholar] [CrossRef] [Green Version]
  25. Hispard, F.; Schuler, D.; de Vaufleury, A.; Scheifler, R.; Badot, P.M.; Dallinger, R. Metal distribution and metallothionein induction after cadmium exposure in the terrestrial snail Helix aspersa (Gastropoda, Pulmonata). Environ. Toxicol. Chem. Int. J. 2008, 27, 1533–1542. [Google Scholar] [CrossRef]
  26. Höckner, M.; Stefanon, K.; de Vaufleury, A.; Monteiro, F.; Pérez-Rafael, S.; Palacios, O.; Capdevila, M.; Atrian, S.; Dallinger, R. Physiological relevance and contribution to metal balance of specific and non-specific metallothionein isoforms in the garden snail, Cantareus aspersus. Biometals 2011, 24, 1079–1092. [Google Scholar] [CrossRef]
  27. Baurand, P.E.; Pedrini-Martha, V.; de Vaufleury, A.; Niederwanger, M.; Capelli, N.; Scheifler, R.; Dallinger, R. Differential expression of metallothionein isoforms in terrestrial snail embryos reflects early life stage adaptation to metal stress. PLoS ONE 2015, 10, e0116004. [Google Scholar] [CrossRef] [Green Version]
  28. Gao, X.; Cao, D.; Liu, J.; Wang, X.; Geng, S.; Liu, B.; Shi, D. Tissue-specific and cation/anion-specific DNA methylation variations occurred in C. virgata in response to salinity stress. PLoS ONE 2013, 8, e78426. [Google Scholar] [CrossRef]
  29. MWFEP (Romania Ministry of Waters, Forests and Environmental Protection). Ordinul nr. 756/1997 al Ministerului Apelor, Padurilor şi Protecţiei Mediului Pentru Aprobarea Reglementării Privind Evaluarea Poluarii Mediului Modificat de Ordinul nr. 1144/2002 al Ministerului Apelor si Protectiei Mediului. 2002. Available online: http://www.unimed.ro/Ordin%20nr.%201144-2002.pdf (accessed on 28 July 2021).
  30. Comission European. Setting maximum levels for certain contaminants in foodstuffs. Regulation 2006, 364, 5–24. [Google Scholar]
  31. Hödl, E.; Felder, E.; Chabicovsky, M.; Dallinger, R. Cadmium stress stimulates tissue turnover in Helix pomatia: Increasing cell proliferation from metal tolerance to exhaustion in molluscan midgut gland. Cell Tissue Res. 2010, 341, 159–171. [Google Scholar] [CrossRef]
  32. Jurkowska, R.Z.; Ceccaldi, A.; Zhang, Y.; Arimondo, P.B.; Jeltsch, A. DNA methyltransferase assays. Methods Mol. Biol. 2011, 791, 157–177. [Google Scholar]
  33. Kurdyukov, S.; Bullock, M. DNA methylation analysis: Choosing the right method. Biology 2016, 5, 3. [Google Scholar] [CrossRef]
  34. Höckner, M.; Stefanon, K.; Schuler, D.; Fantur, R.; De Vaufleury, A.; Dallinger, R. Coping with cadmium exposure in various ways: The two Helicid snails Helix pomatia and Cantareus aspersus share the metal transcription factor-2, but differ in promoter organization and transcription of their Cd-metallothionein genes. J. Exper. Zool. A Ecol. Genet. Physiol. 2009, 311, 776–787. [Google Scholar] [CrossRef]
  35. Genereux, D.P.; Johnson, W.C.; Burden, A.F.; Stöger, R.; Laird, C.D. Errors in the bisulfite conversion of DNA: Modulating inappropriate-and failed-conversion frequencies. Nucleic Acids Res. 2008, 36, e150. [Google Scholar] [CrossRef] [Green Version]
  36. Laskowski, R.; Hopkin, S.P. Effect of Zn, Cu, Pb, and Cd on fitness in snails (Helix aspersa). Ecotoxicol. Environ. Saf. 1996, 34, 59–69. [Google Scholar] [CrossRef] [Green Version]
  37. Russell, L.K.; DeHaven, J.I.; Botts, R.P. Toxic effects of cadmium on the Garden snail (Helix aspersa). Bull. Environ. Contam. Toxicol. 1981, 26, 634–640. [Google Scholar] [CrossRef]
  38. Gomot, A. Dose-dependent efects of cadmium on the growth of snails in toxicity bioassays. Arch. Environ. Contam. Toxicol. 1997, 33, 209–216. [Google Scholar] [CrossRef]
  39. Czarnołęski, M.; Kozłowski, J.; Dumiot, G.; Bonnet, J.C.; Mallard, J.; Dupont-Nivet, M. Scaling of metabolism in Helix aspersa snails: Changes through ontogeny and response to selection for increased size. J. Exp. Biol. 2008, 211, 391–400. [Google Scholar] [CrossRef] [Green Version]
  40. Draghici, G.A.; Dehelean, C.A.; Pinzaru, I.; Bordean, D.M.; Pop, G.; Nica, D. An 112-Days Experiment on Dietary Cadmium Retention in Hepatopancreas in Adult Cantareus aspersus Snails. Rev. Chim. 2019, 70, 2803–2804. [Google Scholar] [CrossRef]
  41. Gomot-de Vaufleury, A.; Kerhoas, I. Effects of cadmium on the reproductive system of the land snail Helix aspersa. Bull. Environ. Contam. Toxicol. 2000, 64, 434–442. [Google Scholar] [CrossRef]
  42. Scheifler, R.; Gomot-De Vaufleury, A.; Toussaint, M.L.; Badot, P.M. Transfer and effects of cadmium in an experimental food chain involving the snail Helix aspersa and the predatory carabid beetle Chrysocarabus splendens. Chemosphere 2002, 48, 571–579. [Google Scholar] [CrossRef]
  43. Joe, M. The Epigenetic Effect of Trematode Infection on the Snail Host Zeacumantus subcarinatus. Master’s Thesis, The University of Otago, Dunedin, New Zeeland, 2013. Available online: https://ourarchive.otago.ac.nz/handle/10523/4163 (accessed on 29 July 2021).
  44. Müller, R.; Charaf, S.; Scherer, C.; Oppold, A.; Oehlmann, J.; Wagner, M. Phenotypic and epigenetic effects of vinclozolin in the gastropod Physella acuta. J. Molluscan Stud. 2016, 82, 320–327. [Google Scholar] [CrossRef] [Green Version]
  45. Drechsel, V.; Schauer, K.; Šrut, M.; Höckner, M. Regulatory plasticity of earthworm wMT-2 gene expression. Int. J. Mol. Sci. 2017, 18, 1113. [Google Scholar] [CrossRef] [Green Version]
  46. Gavery, M.R.; Roberts, S.B. DNA methylation patterns provide insight into epigenetic regulation in the Pacific oyster (Crassostrea gigas). BMC Genom. 2010, 1, 483. [Google Scholar] [CrossRef] [Green Version]
  47. Li, Y.; Wang, R.; Xun, X.; Wang, J.; Bao, L.; Thimmappa, R.; Ding, J.; Jiang, J.; Zhang, L.; Li, T.; et al. Sea cucumber genome provides insights into saponin biosynthesis and aestivation regulation. Cell Discov. 2018, 4, 1–17. [Google Scholar] [CrossRef]
  48. Egg, M.; Hoeckner, M.; Brandstaetter, A.; Schuler, D.; Dallinger, R. Structural and bioinformatic analysis of the Roman snail Cd-Metallothionein gene uncovers molecular adaptation towards plasticity in coping with multifarious environmental stress. Mol. Ecol. 2009, 18, 2426–2443. [Google Scholar] [CrossRef]
  49. Pedrini-Martha, V.; Niederwanger, M.; Kopp, R.; Schnegg, R.; Dallinger, R. Physiological, diurnal and stress-related variability of cadmium-metallothionein gene expression in land snails. PLoS ONE 2016, 11, e0150442. [Google Scholar] [CrossRef] [Green Version]
  50. Draghici, G.A.; Nica, D.; Dehelean, C.; Susan, R.; Susan, M. DNMT gene expression in Cantareus aspersus snails exposed to low cadmium dose. Naunyn Schmiedebergs Arch. Pharmacol. 2020, 393, 53–54. [Google Scholar]
  51. Parmakelis, A.; Kotsakiozi, P.; Kontos, C.K.; Adamopoulos, P.G.; Scorilas, A. The transcriptome of a “sleeping” invader: De novo assembly and annotation of the transcriptome of aestivating Cornu aspersum. BMC Genom. 2017, 18, 491. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Toxics 09 00306 g001 550
Figure 1. ELISA standard curve. The upper right shows the standard curve equation and the corresponding coefficient of determination. The x axis illustrates the 5mC values of the negative control (0%) and the positive controls (0.1%, 0.2%, 0.5%, 1%, 2% and 5%) expressed as ln (5mC% + 1), whereas the y-axis shows the corresponding absorbances.
Figure 1. ELISA standard curve. The upper right shows the standard curve equation and the corresponding coefficient of determination. The x axis illustrates the 5mC values of the negative control (0%) and the positive controls (0.1%, 0.2%, 0.5%, 1%, 2% and 5%) expressed as ln (5mC% + 1), whereas the y-axis shows the corresponding absorbances.
Toxics 09 00306 g001
Toxics 09 00306 g002 550
Figure 2. Effect of dietary cadmium on (a) hepatopancreas cadmium levels (n = 4) and (b) body weight gain/loss (n = 8). Metal data are shown on a log10 scale, whereas weight data are given as absolute values. Data sets are depicted as mean (point) with one standard error (box) and one standard deviation (error bar) (Dunnet’s tests, ***—p < 0.001, **—p < 0.01, *—p < 0.05).
Figure 2. Effect of dietary cadmium on (a) hepatopancreas cadmium levels (n = 4) and (b) body weight gain/loss (n = 8). Metal data are shown on a log10 scale, whereas weight data are given as absolute values. Data sets are depicted as mean (point) with one standard error (box) and one standard deviation (error bar) (Dunnet’s tests, ***—p < 0.001, **—p < 0.01, *—p < 0.05).
Toxics 09 00306 g002
Toxics 09 00306 g003 550
Figure 3. Effect of dietary cadmium on genome-wide 5mC levels in DNA of the hepatopancreas of C. aspersum adults (n = 4). Data are shown on a log10 scale as mean (point) with one standard error (box) and one standard deviation (error bar) (Dunnet’s tests, ***—p < 0.001, **—p < 0.01, *—p < 0.05).
Figure 3. Effect of dietary cadmium on genome-wide 5mC levels in DNA of the hepatopancreas of C. aspersum adults (n = 4). Data are shown on a log10 scale as mean (point) with one standard error (box) and one standard deviation (error bar) (Dunnet’s tests, ***—p < 0.001, **—p < 0.01, *—p < 0.05).
Toxics 09 00306 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Georgescu, M.; Drăghici, G.A.; Oancea, E.-F.; Dehelean, C.A.; Şoica, C.; Vlăduţ, N.-V.; Nica, D.V. Effects of Cadmium Sulfate on the Brown Garden Snail Cornu aspersum: Implications for DNA Methylation. Toxics 2021, 9, 306. https://doi.org/10.3390/toxics9110306

AMA Style

Georgescu M, Drăghici GA, Oancea E-F, Dehelean CA, Şoica C, Vlăduţ N-V, Nica DV. Effects of Cadmium Sulfate on the Brown Garden Snail Cornu aspersum: Implications for DNA Methylation. Toxics. 2021; 9(11):306. https://doi.org/10.3390/toxics9110306

Chicago/Turabian Style

Georgescu, Marius, George Andrei Drăghici, Eliza-Florentina Oancea, Cristina Adriana Dehelean, Codruţa Şoica, Nicolae-Valentin Vlăduţ, and Dragoș Vasile Nica. 2021. "Effects of Cadmium Sulfate on the Brown Garden Snail Cornu aspersum: Implications for DNA Methylation" Toxics 9, no. 11: 306. https://doi.org/10.3390/toxics9110306

APA Style

Georgescu, M., Drăghici, G. A., Oancea, E. -F., Dehelean, C. A., Şoica, C., Vlăduţ, N. -V., & Nica, D. V. (2021). Effects of Cadmium Sulfate on the Brown Garden Snail Cornu aspersum: Implications for DNA Methylation. Toxics, 9(11), 306. https://doi.org/10.3390/toxics9110306

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop