Next Article in Journal
Burnout and Depression in Portuguese Healthcare Workers during the COVID-19 Pandemic—The Mediating Role of Psychological Resilience
Previous Article in Journal
Exploring the Role of Personal Demands in the Health-Impairment Process of the Job Demands-Resources Model: A Study among Master Students
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Correlation between Polycyclic Aromatic Hydrocarbons in Wharf Roach (Ligia spp.) and Environmental Components of the Intertidal and Supralittoral Zone along the Japanese Coast

1
Botanical Garden, Institute of Nature and Environmental Technology, Kanazawa University, Kakuma, Kanazawa, Ishikawa 920-1192, Japan
2
Faculty of Agriculture, Kyushu University, 744 Motooka, Nishi-ku, Fukuoka 819-0395, Japan
3
Department of Biological Sciences, Graduate School of Science and Technology, Kumamoto University, 2-39-1 Kurokami, Chuo-ku, Kumamoto 860-8555, Japan
4
Graduate School of Life and Environmental Sciences, Shimane University, 1060 Nishitsugawa-machi, Matsue, Shimane 690-8504, Japan
5
Education and Research Center for Marine Resources and Environment, Faculty of Fisheries, Kagoshima University, 4-50-20 Shimoarata, Kagoshima 890-0056, Japan
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(2), 630; https://doi.org/10.3390/ijerph18020630
Submission received: 14 December 2020 / Revised: 29 December 2020 / Accepted: 5 January 2021 / Published: 13 January 2021

Abstract

:
Polycyclic aromatic hydrocarbon (PAH) concentrations in wharf roach (Ligia spp.), as an environmental indicator, and in environmental components of the intertidal and supralittoral zones were determined, and the PAH exposure pathways in wharf roach were estimated. Wharf roaches, mussels, and environmental media (water, soil and sand, and drifting seaweed) were collected from 12 sites in Japan along coastal areas of the Sea of Japan. PAH concentrations in wharf roaches were higher than those in mussels (median total of 15 PAHs: 48.5 and 39.9 ng/g-dry weight (dw), respectively) except for samples from Ishikawa (wharf roach: 47.9 ng/g-dw; mussel: 132 ng/g-dw). The highest total PAH concentration in wharf roach was from Akita (96.0 ng/g-dw), followed by a sample from Niigata (85.2 ng/g-dw). Diagnostic ratio analysis showed that nearly all PAHs in soil and sand were of petrogenic origin. Based on a correlation analysis of PAH concentrations between wharf roach and the environmental components, wharf roach exposure to three- and four-ring PAHs was likely from food (drifting seaweed) and from soil and sand, whereas exposure to four- and five-ring PAHs was from several environmental components. These findings suggest that the wharf roach can be used to monitor PAH pollution in the supralittoral zone and in the intertidal zone.

1. Introduction

Polycyclic aromatic hydrocarbons (PAH) are composed of two or more fused aromatic rings without modification and include more than 100 chemical species. PAHs found in the environment are either of petrogenic origin, such as from natural oil seeps and accidental oil spills [1,2,3], or are of pyrogenic origin, which includes volcanic activity [4], the combustion of fossil fuels [5] and organic matter [6], and automobile exhaust [7]. PAHs are ubiquitous contaminants in the atmosphere [8,9], terrestrial soil [10,11], marine and freshwater sediment [12,13,14], seawater [15], and fresh water [16,17], and they can accumulate in a broad range of wildlife [18,19,20,21]. In addition, PAHs are a common contaminant found in the human body [22,23]. However, field data of PAH pollution in marine invertebrates are limited [24], except in bivalves [2,3,25,26]. Field surveys of PAH contamination in invertebrates in the supralittoral zone are particularly sparse [27].
PAH toxicity has been extensively studied. These compounds are well known for their carcinogenicity [28,29], mutagenicity [30], endocrine-disrupting activity [31], and developmental toxicity [32,33]. Owing to their toxicity and distribution, the United States Environmental Protection Agency has selected 16 representative PAHs as priority compounds for environmental monitoring. The International Agency for Research on Cancer has classified benzo[a]pyrene into group 1, indicating that it is “carcinogenic to humans”; dibenz[a,h]anthracene is in group 2A, indicating it is “probably carcinogenic to humans”.
Surveying the intertidal and supralittoral zones contributes to understanding the toxicity in nearshore organisms. Furthermore, many contaminants accumulate in this area. Therefore, a wide range of pollutants, such as heavy metals [34], pesticides [35], persistent organic pollutants [36,37], pharmaceuticals [38], antibiotics [39], and PAHs in seawater and sediment have been monitored. These contaminants are continuously deposited in coastal areas, with vertical distributions in the intertidal and supralittoral zones. While the intertidal zone is a common focus, the supralittoral zone is also a site of contamination and should also be monitored.
For this purpose, researchers often sample mussels [19,34,40] to conduct biomonitoring in coastal areas. However, mussels may not be appropriate or accessible for environmental biomonitoring in the supralittoral zone. To address this, our research group has identified the wharf roach (Ligia spp.) as a suitable environmental indicator for biomonitoring in the supralittoral zone. The wharf roach has a global distribution [41]. It has a relatively small habitat area [42] and is omnivorous [43]; therefore, it is exposed to environmental pollutants because it feeds on land above the tidal line [44]. Environmental pollutants such as radiocesium and radiosilver from the Fukushima Daiichi Nuclear Power Plant [45], tributyltin, dibutyltin [44], and PAHs and alkylated PAHs [46] have been assessed by sampling wharf roaches in the supralittoral zone. Although appropriate as a novel environmental indicator of PAHs in the supralittoral zone, wharf roaches live in the boundary region between the marine and terrestrial environments, and their exposure routes to PAH compounds are unclear.
In this study, the correlations between PAH concentrations in wharf roaches and in environmental components were investigated using field samples, and the exposure pathways by which wharf roaches are exposed to PAHs are estimated. Wharf roaches, mussels, seawater, soil and sand, and drifting seaweed were collected at 12 sampling sites in Japan along coastal areas around the Sea of Japan. These sites were chosen because transboundary PAH air pollution from the continent has been reported in this area [8].

2. Materials and Methods

2.1. Chemicals and Reagents

A PAH standard mixture containing the following 16 PAHs was purchased from Supelco (Bellefonte, PA, USA): naphthalene (Nap); acenaphthylene (Acy); acenaphthene (Ace); fluorene (Flu); phenanthrene (Phe); anthracene (Ant); fluoranthene (Flut); pyrene (Pyr); benzo[a]anthracene (BaA); chrysene (Chr); benzo[b]fluoranthene (BbF); benzo[k]fluoranthene (BkF); benzo[a]pyrene (BaP); benzo[g,h,i]perylene (BP); dibenzo[a,h]anthracene (DA); and indeno[1,2,3-cd]pyrene (IP). The following five deuterium isotopes of PAHs were purchased from Wako Pure Chemical Industries, Ltd. (Osaka, Japan): D8-Nap; D10-Ace; D10-Phe; D10-Pyr; and D12-BaP. All organic solvents were of pesticide residue and polychlorinated biphenyl test grade. Other chemicals and ultrapure water were of analytical grade.

2.2. Sampling

Wharf roaches (Ligia spp.) were collected in September 2018 from 12 sampling sites in Japan located along coastal areas of the Sea of Japan (Figure 1, Table 1). Wharf roaches were collected by hand and using tools.
Mussels (Mytilus spp.), seawater, soil and sand, and drifting seaweed were also collected. The sampling sites were located in the habitable zone for both wharf roaches and mussels and were suitable for sampling seawater, soil and sand, and drifting seaweed. Mussels were collected from the shore area by hand or using a metal spatula. We collected four liters of seawater from the water surface at each location using brown glass bottles. After collection, the seawater samples were spiked with PAH internal standard mixture and extracted on-site using an Empore C18 extraction disc (3M, Maplewood, MN, USA). Soil and sand and drifting seaweed samples were randomly collected from the supralittoral zone near the wharf roach sampling location using a metal shovel or tweezers. All of the samples were immediately stored at 4 °C and sent to the laboratory where they were kept at −30 °C. Debris was frequently observed at the sampling sites (supplementary materials: Figure S1).

2.3. PAH Analysis

The analytical methods were conducted as in Ito et al. [19] and Honda et al. [46], with slight modifications. Briefly, the samples (wharf roaches, mussels, drifting seaweed, and soil and sand) were placed in a glass container and freeze-dried for 48 h (FDU-1200; Eyela, Tokyo, Japan). After freeze-drying, 5–10 wharf roach samples were pooled and homogenized (whole body), 2–5 mussels were pooled and homogenized (whole soft tissue), and approximately 1 g of drifting seaweed was pooled and homogenized. Large-grain sand, pebbles, and other foreign matter were removed from the soil and sand samples to produce a roughly uniform grain size. Approximately 0.3 g of wharf roach, mussel, or drifting seaweed and 1 g of soil and sand was weighed and placed into a glass tube. After spiking with the internal standard mixture (20 ng of D8-Nap, 4 ng each of D10-Ace, D10-Phe, D10-Pyr, and D12-BaP), PAHs were extracted twice with 20 mL of dichloromethane/hexane (50%:50%, v/v). The extracts were saponified with 6 mL of 1 M potassium hydroxide in ethanol at 90 °C for 1 h. After saponification, the samples were extracted with hexane. The extracts were washed with ultrapure water and dehydrated with anhydrous sodium sulfate. The dehydrated samples were purified with silica gel (3% water content, v/v) and eluted with hexane and 3% acetone in hexane (v/v). The eluates were concentrated to 0.2 mL with N2 gas and transferred into gas chromatography (GC) vials for analysis. Extraction of PAHs from the C18 Empore discs was performed by cutting the disc into pieces and extracting it similarly to the other samples. The uncertainty of single measurements of individual samples was alleviated by sample pooling.
PAHs were quantified using a gas chromatograph (8890 GC system; Agilent, Santa Clara, CA, USA) equipped with a mass spectrometer (5977B GC/MSD; Agilent). Target analytes were separated using a DB-5ms column (0.25 mm i.d. × 30 m length, film thickness 0.25 μm; Agilent). The target analytes were quantified in selected ion monitoring mode. Details of the GC and mass spectrometry parameters are provided in Table S1. The monitored mass of target analytes, the internal standards, and the approximate retention times are provided in Table S2.

2.4. Quality Assurance and Quality Control

A method blank and a matrix-spiked sample were analyzed with each batch of 10 samples. Method blanks contained PAH concentrations that ranged from less than the limit of detection (LOD) to 2.79 ng/mL. Matrix effects were corrected using the isotope dilution method of quantification and the retention time. The recovery of PAHs in spiked samples (spiked content: 2–40 ng) ranged from 84.2 (Flu) to 118% (BaA). The instrumental LODs of each analyte ranged from 0.20 to 9.34 ng/mL. Because of the high background and relatively low recovery, Nap was excluded from the analysis.

2.5. Statistical Analysis

PAHs with concentrations below the LOD were substituted with the value of LOD/2 for the statistical analysis. R software (ver. 4.0.2, R Development Core Team) was used to conduct the statistical analysis, and the significance level was set at p ≤ 0.05. As a non-parametric statistical test, Spearman’s rank–order correlation coefficient was used to examine significant correlations of PAH concentrations among samples. Principal component analysis (PCA) was used to discern whether there were differences between the various sample types. The PCA was performed using SPSS statistics software (ver. 25; IBM, Armonk, NY, USA). Three factors retained based on the scree plot. PCA scores and loading plots were constructed using Microsoft Excel.

3. Results and Discussion

3.1. PAH Concentrations in Organisms

In wharf roach, the highest concentration of the total of 15 PAHs (ΣPAH) was detected in the sample from Akita (96.0 ng/g-dry weight (dw); Figure 2, Table 2), followed by samples from Niigata, Saga, Yamagata, Fukuoka, Aomori, Ishikawa, Shimane, Hyogo, Kyoto, Nagasaki, and Yamaguchi. The detected ΣPAH concentrations (median: 48.5 ng/g-dw, detected range: 26.9 to 96.0 ng/g-dw) were similar to previously reported concentrations, except for samples that were previously collected in highly polluted areas (median: 47.0 ng/g-dw, detected range: 28.6 to 72.2 ng/g-dw) [46]. The PAH composition was largely consistent among the 12 sampling sites (Figure 2). The predominant PAH species was Pyr (median: 24.9 ng/g-dw), followed by Phe, Flut, and Acy (median: 7.71, 6.44, and 5.40 ng/g-dw, respectively). In a previous study, the dominant PAH species were Phe, followed by Ace and Pyr [46]. Thus, the current results differ from previous observations. The composition differences reflect differences in pollution among the sampling sites and differences in contamination sources between sampling seasons (current study: summer, previous study: fall–winter; [47]). High molecular weight PAHs (HMW-PAHs) (i.e., five- and six-ring PAHs) were detected at lower concentrations than low molecular weight PAHs (LMW-PAHs) (i.e., three- and four-ring PAHs), which were the most common. Three-ring PAHs were contributed at 30% and four-ring PAHs at 68% of ΣPAH concentration. It is generally believed that HMW-PAHs do not highly accumulate in organisms because of a lower intake efficiency [48].
A consistent trend of high Pyr concentrations (median: 11.6 ng/g-dw) was observed in mussel tissue. The highest ΣPAH concentration was detected in a sample from Ishikawa (132 ng/g-dw; Figure 3), followed by samples from Aomori, Hyogo, Akita, Kyoto, Fukuoka, Nagasaki, Shimane, Niigata, Yamagata, Yamaguchi, and Saga. The detected range of ΣPAH among the 12 sampling sites was 16.6 to 132 ng/g-dw. This finding is similar to that of previous studies; 16 PAHs, 87.3–361 ng/g-dw [19]; 18 PAHs, 15.2 to 527 ng/g-dw (both values were converted from initial reported values of 2.6 to 90 ng/g-ww by using an 82.9% moisture content) [40]. Of the 15 PAHs, Pyr was the predominant species (median: 11.6 ng/g-dw), followed by Phe and Flut (7.10 and 5.83 ng/g-dw, respectively). Additionally, the concentrations of HMW-PAHs were greater in mussels than in wharf roaches. The relative standard deviation of ΣPAH among the 12 sampling sites was larger for mussels than for wharf roaches (0.63 and 0.39, respectively). These differences reflect differences in the PAH pollution in the supralittoral and intertidal zones and suggest that wharf roaches and mussels differ in their exposure pathways and bioaccumulation capability.

3.2. PAH Concentrations in Environmental Media

The ΣPAH composition was consistent across seawater samples (Figure 4). The detected range of ΣPAH among the 12 sampling sites was 7.79 to 19.1 ng/L (mean 13.4 ng/L). This finding is similar to that of previous studies; 13 PAHs, 6.83 to 13.81 ng/L (mean 9.4 ng/L) in the Japan Sea [49]; 13 PAHs, 10.9 to 29.7 ng/L (mean 19.6 ng/L) in the Sea of Japan and East Sea [50]. Phe was the predominant species (median: 5.67 ng/L). In general, there was a high contribution of three-ring PAHs to ΣPAH (87%). Comparatively, four-, five-, and six-ring PAHs had very low contributions to ΣPAH (11%, 1.5%, and 0.4%, respectively). The concentrations of LMW-PAHs were higher than those of HMW-PAHs owing to their increased water solubility [47,51].
The detected range of ΣPAH in soil and sand among the 12 sampling sites was 2.06 to 506 ng/g-dw (median 22.0 ng/g-dw). This finding is lower than that of previous studies; 18 PAHs in sea sediment, 6.40 to 7765 ng/g-dw in Osaka Bay, Japan [12]; 16 PAHs in wharf soil, 842 ± 203 ng/g-dw in Guangdong, China [52]. Phe was the predominant species (median: 5.67 ng/g-dw). At seven sampling sites, that had relatively small particles and a larger soil fraction, the soil and sand samples contained a larger amount of HMW-PAHs than the seawater samples (Aomori, Ishikawa, Kyoto, Hyogo, Fukuoka, Saga, and Nagasaki; Figure 5). The other five sites, which had relatively large sand particles, had PAH compositions similar to that of the seawater samples (Akita, Yamagata, Niigata, and Shimane). The PAH with the highest median concentration detected among the 12 sites was Pyr (4.91 ng/g-dw), followed by Flut and Chr (2.35 and 1.57 ng/g-dw, respectively). There were four-ring PAHs predominant in the ΣPAH at 61%, which was much higher than the percentage of three-ring PAHs (20%). Additionally, five- and six-ring PAHs comprised a much higher percentage of ΣPAH in soil and sand (12.1% and 6.9%, respectively) than in seawater (1.5% and 0.4%, respectively). PAH hydrophobicity is positively correlated with ring number and molecular weight. Higher molecular weight PAHs tend to partition into the particle phase. Therefore, HMW-PAHs were more likely to be found in the soil and sand samples than in the seawater samples. Additionally, the sampling sites of Ishikawa and Kyoto were located in a wharf area, while the other sites were located in beach or rocky areas. Soil samples were collected from a small ditch near a wharf in Ishikawa and Kyoto, and they contained higher concentrations of PAHs (506 and 138 ng/g-dw, respectively) compared with the other 10 sites (median: 15.7 ng/g-dw). This difference was not observed for the other analyzed samples.
Drifting seaweed is a major food source for wharf roaches. The detected range of ΣPAH in drifting seaweed among the 12 sampling sites was 11.8 to 109 ng/g-dw (mean 47.6 ng/g-dw). This finding is higher than that of previous studies; 16 PAHs in living seaweed, 1.0 to 56.4 ng/g-dw (mean 7.0 ng/g-dw) in Venice, Italy [53]. The PAH composition in drifting seaweed was roughly separated into two groups, with two exceptions (Figure 6) (notably, drifting seaweed samples could not be found at the Shimane site). One group was largely consistent with the PAH trends in seawater (Akita, Yamagata, Niigata, Kyoto, and Shimane). The other group was more similar to the composition found in soil and sand samples (Yamaguchi, Fukuoka, Saga, and Nagasaki). There were two exceptions: Aomori and Ishikawa. As with other environmental media, Pyr was the predominant species (median: 14.4 ng/g-dw). However, Aomori had a higher concentration and percentage of BP (6.95 ng/g-dw, 20%), and Ishikawa had a higher concentration and percentage of Flut (38%). Beyond Pyr, Flut and Phe also had substantial concentrations (6.85 and 4.11 ng/g-dw, respectively). Seaweed can accumulate PAHs from seawater during its lifecycle [53]. Therefore, it is likely that these PAHs were transferred from the seawater and nearby soil and sand into the seaweed. Additionally, aerial deposition on drifting seaweed is a possible exposure source [54].

3.3. PCA

The PCA confirmed some of the similarities found among the various sample types (Figure 7). In the first principal component (PC) (accounting for 73.4% of the variance), wharf roach and seawater were closely related, whereas soil and sand and mussels formed a related group. The loading plot indicated that the separations in the first PC were largely driven by differences between the three-ring PAHs (i.e., Acy and Ace) and the larger molecular weight PAHs (Figure S2). In general, concentrations of the larger five- and six-ring PAHs were higher in soil and sand, seaweed, and mussels than in wharf roach and seawater (Table S3). By contrast, seawater samples contained greater concentrations of the smaller three-ring PAHs. These differences may account for the separation between these groups in PC1. Compared with the other samples, the relative similarities between soil and sand and mussels were expected given the hydrophobic tendency of PAHs, whereby they partition onto particles or into the high lipid content of mussels. This partitioning also explains the results of the correlation analysis, where four- and five-ring PAHs were particularly well correlated between seawater and wharf roaches. The close association between wharf roaches and seawater supports the idea that the wharf roach can be a useful proxy for assessing surface-layer PAH contamination. However, there appear to be some differences, as seen in PC2, though this accounts for only 14.6% of the variance. Although the loading plot indicated that these separations are caused by differences between the four-ring PAHs and the five and six-ring PAHs, it is not clear why these differences exist.

3.4. Identification of Pollution Sources

PAHs are derived from several sources, such as petrogenic, pyrogenic, biogenic, and diagenetic sources [25,55]. PAH pollution in the marine environment is generally thought to be mainly derived from pyrolytic and petrogenic sources. Diagnostic ratios comparing PAH isomers are commonly used to identify pollution sources [2,56,57]. Based on these ratios, petrogenic sources were found to predominate in the soil and sand samples collected at nine of the 12 sampling sites (Figure 8). Pyrolytic sources were detected in samples from three sites (Ishikawa, Kyoto, and Nagasaki), though these associations were weak. Pollution sources for the seawater samples were different from those of the soil and sand samples. The seawater samples likely had a high contribution from pyrolytic sources. Intertidal and supralittoral zones collect floating oil derived from oil contamination from boats or domestic wastewater [5] in the marine environment and are easily polluted by combustion sources located in the surrounding terrestrial area. The ratios found in this study may reflect these phenomena.

3.5. Correlation between PAH Concentrations in Wharf Roaches and Environmental Components

Wharf roaches were presumed to be exposed to PAHs through the soil and sand, seawater, and dietary seaweed exposure pathways. Wharf roaches are terrestrial isopods that occupy the supralittoral zone, but they take in seawater through their legs [58]. They feed on drifting seaweed and biofilms [43]. The correlations between PAH ring numbers and organisms and substrates were analyzed. The three-ring PAHs were significantly correlated between wharf roach and drifting seaweed and soil and sand (Spearman’s rank correlation, p = 0.01, <0.01; rho = 0.34, 0.51, respectively); the four-ring PAHs were significantly correlated between wharf roach and drifting seaweed, seawater, and soil and sand (p = 0.01, <0.01, 0.01; rho = 0.62, 0.59, 0.36, respectively); the five-ring PAHs were significantly negatively correlated between wharf roach and seawater (p = 0.02; rho = −0.33). Detected concentrations of five-ring PAHs were relatively low in all analyzed samples. Therefore, it is likely that there was no significant variation in five-ring PAH concentrations to result in a significant correlation between wharf roach and seawater. Because several samples had concentrations of six-ring PAHs that were less than the LOD, correlations for these PAHs could not be analyzed (detection frequency: 80%).
For mussels, the major exposure pathway of LMW-PAHs is generally considered to be water, and that of HMW-PAHs is considered to be particles [59]. This is explained by the hydrophobic properties of PAHs, which make them more likely to partition onto nonpolar particulate surfaces. Mussels continuously filter seawater containing organic matter, which constitutes the majority of its food. For this reason, mussels are liable to accumulate PAHs that are present in the dissolved and particulate phases. The ΣPAH concentration in mussels from Ishikawa (132 ng/g-dw) was the highest among sampling sites, likely because Ishikawa also had the highest soil and sand ΣPAH concentration (506 ng/g-dw).
However, wharf roaches are omnivorous scavengers and are not directly exposed to marine sediment. Food is generally considered the major pathway by which organisms take in PAHs [60]. Surprisingly, the statistical analysis results indicated that LMW-PAHs (i.e., three- and four-ring PAHs) were significantly correlated between wharf roach and food source (i.e., drifting seaweed; rho = 0.34, 0.62). There was also a significant correlation between LMW-PAH concentrations in soil and sand and in wharf roach (rho = 0.51, 0.36). The LMW-PAH composition results indicated that the LMW compounds in wharf roaches had a pattern similar to that in soil and sand (especially in samples from Akita, Niigata, and Shimane; Figure 2 and Figure 5) and to that in drifting seaweed (especially in samples from Akita, Yamagata, Niigata, Kyoto, Hyogo, and Saga; Figure 2 and Figure 6). Some terrestrial isopods accumulate PAHs from contaminated soil [61,62]. In the present study, particles of soil and sand were frequently detected in the intestines of wharf roaches (data not shown). It is possible that wharf roaches directly ingest soil and sand particles to obtain the organic matter adhered to them. This may have been a possible exposure route in the current study. Additionally, van Brummelen et al. [61] reported a negative correlation between biota-to-soil accumulation factors and the PAH Kow. Therefore, soil and sand may have contributed less than other routes to the HMW-PAH exposure of wharf roaches. However, the actual contribution ratio of each exposure pathway and the metabolism of PAHs in isopods remain unclear. It is possible that efficient LMW-PAH metabolism lowered their concentrations in the samples, thereby masking a detectable correlation between wharf roach and seawater.
Compared with LMW-PAHs, HMW-PAHs are metabolized more slowly in the body [63,64]; however, HMW-PAHs (more than five rings) are highly hydrophobic and are found less frequently in aquatic samples. Therefore, it is difficult to evaluate the wharf roach exposure pathway to HMW-PAHs. Wharf roaches take in seawater from the ocean surface [58] where there is often floating oil and related hydrophobic substances. Furthermore, the water surface is relatively hydrophobic compared with the underlying water. Therefore, in addition to drifting seaweed and soil and sand, wharf roaches can possibly accumulate a broad range of PAHs via seawater. However, HMW-PAH concentrations in wharf roaches were negatively correlated with those in seawater. This phenomenon may not have indicated a negative association between these two groups, rather, it could have resulted from the relatively low concentrations of five-ring PAHs in the wharf roach and seawater samples. In general, HMW-PAHs are not efficiently taken up and accumulated by organisms because of their high hydrophobicity [63]. Although not significant, in the samples taken from Saga, the relatively high BP concentration in wharf roach may have been related to the high BP concentrations in the soil/sand and drifting seaweed samples. Therefore, it is possible that wharf roaches are exposed to HMW-PAHs through several environmental components; however, no single component has a predominant contribution because of the low intake efficiency.
Bioconcentration, biomagnification, and biota-sediment accumulation have different exposure pathways and contribution ratios that affect the total accumulation. In addition to drifting seaweed, wharf roaches feed on biofilms and other organic matter distributed throughout the supralittoral zone. Additionally, PAH metabolism varies considerably among species. Therefore, it is necessary to expose wharf roaches to several PAHs under laboratory conditions to assess the actual accumulation system and calculate the bioconcentration and biomagnification factors.

4. Conclusions

The PAH correlation analysis between wharf roaches and the environmental media (seawater, soil and sand, and drifting seaweed) indicated that exposure to LMW-PAHs occurred mainly via soil and sand and drifting seaweed. HMW-PAH exposure occurred via several environmental components. These findings suggest that wharf roaches can reflect PAH pollution in a broad range of coastal components. To our knowledge, this is the first study estimating PAH exposure pathways for an organism in the supralittoral zone.
There are several limitations to this study. Sorption can occur across many particulate types (e.g., microplastics) and requires further investigation. Additionally, exposure may vary by season and should be further clarified.

Supplementary Materials

The following are available online at https://www.mdpi.com/1660-4601/18/2/630/s1, Figure S1: Pollution in the intertidal and supralittoral zones in the Japan coastal area, Figure S2: Loading plot for the first two principal components (accounting for 88% of the variance), Table S1: Gas chromatography (GC) and mass spectrometry (MS) conditions, Table S2: Mass transition and approximate retention time of target analytes and internal standards, Table S3: Total concentration of PAHs grouped by ring number in samples from 12 sampling sites (ng/g-dw or ng/L).

Author Contributions

Conceptualization: M.H., S.U., and Y.O.; methodology: M.H. and K.M.; formal analysis: M.H. and E.N.; investigation: M.H. and K.M.; writing—original draft preparation: M.H. and K.M.; writing—review and editing: M.H., K.M., E.N., S.U., and Y.O.; supervision: Y.O.; project administration: M.H. All authors have read the manuscript and agree with the published version.

Funding

This work was partly supported by a Sasakawa Scientific Research Grant from The Japan Science Society (grant number: 2019-6037) and by the Cooperative Research Program of the Institute of Nature and Environmental Technology, Kanazawa University (acceptance number: 19050, 2019).

Institutional Review Board Statement

Experiments were conducted in accordance with the Japanese Ministry of Education, Culture, Sports, Science and Technology and Kanazawa University guidelines for the protection of animal welfare.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank Kazuichi Hayakawa (Kanazawa University) for advice regarding the study design. We also thank Katherine Thieltges, from Edanz Group (https://en-author-services.edanz.com/) for editing a draft of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Beyer, J.; Trannum, H.C.; Bakke, T.; Hodson, P.V.; Collier, T.K. Environmental effects of the Deepwater Horizon oil spill: A review. Mar. Pollut. Bull. 2016, 110, 28–51. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Koyama, J.; Uno, S.; Nagai, Y.; Anukorn, B. Early monitoring of spilled oil contamination in Rayong, Thailand. Jpn. J. Environ. Toxicol. 2016, 19, 25–33. [Google Scholar] [CrossRef]
  3. Uno, S.; Koyama, J.; Kokushi, E.; Monteclaro, H.; Santander, S.; Cheikyula, J.O.; Miki, S.; Añasco, N.; Pahila, I.G.; Taberna, H.S., Jr.; et al. Monitoring of PAHs and alkylated PAHs in aquatic organisms after 1 month from the Solar I oil spill off the coast of Guimaras Island, Philippines. Environ. Monit. Assess. 2010, 165, 501–515. [Google Scholar] [CrossRef] [PubMed]
  4. Abdel-Shafy, H.I.; Mansour, M.S. A review on polycyclic aromatic hydrocarbons: Source, environmental impact, effect on human health and remediation. Egypt. J. Pet. 2016, 25, 107–123. [Google Scholar] [CrossRef] [Green Version]
  5. Mojiri, A.; Zhou, J.L.; Ohashi, A.; Ozaki, N.; Kindaichi, T. Comprehensive review of polycyclic aromatic hydrocarbons in water sources, their effects and treatments. Sci. Total Environ. 2019, 696, 133971. [Google Scholar] [CrossRef]
  6. Kim, E.J.; Oh, J.E.; Chang, Y.S. Effects of forest fire on the level and distribution of PCDD/Fs and PAHs in soil. Sci. Total Environ. 2003, 311, 177–189. [Google Scholar] [CrossRef]
  7. Liu, Y.; Wang, S.; Lohmann, R.; Yu, N.; Zhang, C.; Gao, Y.; Zhao, J.; Ma, L. Source apportionment of gaseous and particulate PAHs from traffic emission using tunnel measurements in Shanghai, China. Atmos. Environ. 2015, 107, 129–136. [Google Scholar] [CrossRef] [Green Version]
  8. Hayakawa, K.; Tang, N.; Nagato, E.G.; Toriba, A.; Sakai, S.; Kano, F.; Goto, S.; Endo, O.; Arashidani, K.I.; Kakimoto, H. Long term trends in atmospheric concentrations of polycyclic aromatic hydrocarbons and nitropolycyclic aromatic hydrocarbons: A study of Japanese cities from 1997 to 2014. Environ. Pollut. 2018, 233, 474–482. [Google Scholar] [CrossRef]
  9. Slezakova, K.; Pires, J.C.M.; Castro, D.; Alvim-Ferraz, M.C.M.; Delerue-Matos, C.; Morais, S.; Pereira, M.C. PAH air pollution at a Portuguese urban area: Carcinogenic risks and sources identification. Environ. Sci. Pollut. Res. 2013, 20, 3932–3945. [Google Scholar] [CrossRef]
  10. Offiong, N.A.O.; Inam, E.J.; Etuk, H.S.; Essien, J.P. Current status and challenges of remediating petroleum-derived PAHs in soils: Nigeria as a case study for developing countries. Remediation 2019, 30, 65–75. [Google Scholar] [CrossRef]
  11. Titaley, I.A.; Simonich, S.L.M.; Larsson, M. Recent Advances in the Study of the Remediation of Polycyclic Aromatic Compound (PAC)-Contaminated Soils: Transformation Products, Toxicity, and Bioavailability Analyses. Environ. Sci. Technol. Lett. 2020, 7, 873–882. [Google Scholar] [CrossRef]
  12. Miki, S.; Uno, S.; Ito, K.; Koyama, J.; Tanaka, H. Distributions of polycyclic aromatic hydrocarbons and alkylated polycyclic aromatic hydrocarbons in Osaka Bay, Japan. Mar. Pollut. Bull. 2014, 85, 558–565. [Google Scholar] [CrossRef] [PubMed]
  13. Sampei, Y.; Uraoka, S.I.; Ono, T.; Dettman, D.L. Polycyclic aromatic hydrocarbons (PAHs) in sediment cores from lakes Shinji and Nakaumi, SW Japan: A proxy of recent fire events in the watershed. Estuar. Coast. Shelf Sci. 2019, 226, 106269. [Google Scholar] [CrossRef]
  14. Wang, Q.; Peng, F.; Chen, Y.; Jin, L.; Lin, J.; Zhao, X.; Yin, J.; Li, J.Y. Heavy metals and PAHs in an open fishing area of the East China Sea: Multimedia distribution, source diagnosis, and dietary risk assessment. Environ. Sci. Pollut. Res. 2019, 26, 21140–21150. [Google Scholar] [CrossRef]
  15. Sarria-Villa, R.; Ocampo-Duque, W.; Páez, M.; Schuhmacher, M. Presence of PAHs in water and sediments of the Colombian Cauca River during heavy rain episodes, and implications for risk assessment. Sci. Total Environ. 2016, 540, 455–465. [Google Scholar] [CrossRef]
  16. Chizhova, T.; Koudryashova, Y.; Prokuda, N.; Tishchenko, P.; Hayakawa, K. Polycyclic Aromatic Hydrocarbons in the Estuaries of Two Rivers of the Sea of Japan. Int. J. Environ. Res. Public Health 2020, 17, 6019. [Google Scholar] [CrossRef]
  17. Guo, W.; He, M.; Yang, Z.; Lin, C.; Quan, X.; Men, B. Distribution, partitioning and sources of polycyclic aromatic hydrocarbons in Daliao River water system in dry season, China. J. Hazard. Mater. 2009, 164, 1379–1385. [Google Scholar] [CrossRef]
  18. Hamilton, P.B.; Cowx, I.G.; Oleksiak, M.F.; Griffiths, A.M.; Grahn, M.; Stevens, J.R.; Carvalho, G.R.; Nicol, E.; Tyler, C.R. Population-level consequences for wild fish exposed to sublethal concentrations of chemicals–a critical review. Fish Fish. 2016, 17, 545–566. [Google Scholar] [CrossRef]
  19. Ito, K.; Uno, S.; Tanaka, H.; Miki, S.; Kokushi, E.; Yamamoto, M.; Koyama, J.; Añasco, N.C. Distribution of parent and alkylated PAHs in bivalves collected from Osaka bay, Japan. Jpn. J. Environ. Toxicol. 2015, 18, 11–24. [Google Scholar] [CrossRef]
  20. Kannan, K.; Perrotta, E. Polycyclic aromatic hydrocarbons (PAHs) in livers of California sea otters. Chemosphere 2008, 71, 649–655. [Google Scholar] [CrossRef]
  21. Nakata, H.; Uehara, K.; Goto, Y.; Fukumura, M.; Shimasaki, H.; Takikawa, K.; Miyawaki, T. Polycyclic aromatic hydrocarbons in oysters and sediments from the Yatsushiro Sea, Japan: Comparison of potential risks among PAHs, dioxins and dioxin-like compounds in benthic organisms. Ecotoxicol. Environ. Saf. 2014, 99, 61–68. [Google Scholar] [CrossRef] [PubMed]
  22. Guo, Y.; Senthilkumar, K.; Alomirah, H.; Moon, H.B.; Minh, T.B.; Mohd, M.A.; Nakata, H.; Kannan, K. Concentrations and profiles of urinary polycyclic aromatic hydrocarbon metabolites (OH-PAHs) in several Asian countries. Environ. Sci. Technol. 2013, 47, 2932–2938. [Google Scholar] [CrossRef] [PubMed]
  23. Moon, H.B.; Lee, D.H.; Lee, Y.S.; Kannan, K. Occurrence and accumulation patterns of polycyclic aromatic hydrocarbons and synthetic musk compounds in adipose tissues of Korean females. Chemosphere 2012, 86, 485–490. [Google Scholar] [CrossRef] [PubMed]
  24. Lawson, M.C.; Cullen, J.A.; Nunnally, C.C.; Rowe, G.T.; Hala, D.N. PAH and PCB body-burdens in epibenthic deep-sea invertebrates from the northern Gulf of Mexico. Mar. Pollut. Bull. 2020, 111825. [Google Scholar] [CrossRef]
  25. Hylland, K. Polycyclic aromatic hydrocarbon (PAH) ecotoxicology in marine ecosystems. J. Toxicol. Environ. Health Part A 2006, 69, 109–123. [Google Scholar] [CrossRef]
  26. Noh, J.; Kim, H.; Lee, C.; Yoon, S.J.; Chu, S.; Kwon, B.O.; Ryu, J.; Kim, J.J.; Lee, H.; Yim, U.H.; et al. Bioaccumulation of polycyclic aromatic hydrocarbons (PAHs) by the marine clam, Mactra veneriformis, chronically exposed to oil-suspended particulate matter aggregates. Environ. Sci. Technol. 2018, 52, 7910–7920. [Google Scholar] [CrossRef]
  27. Bejarano, A.C.; Michel, J. Oil spills and their impacts on sand beach invertebrate communities: A literature review. Environ. Pollut. 2016, 218, 709–722. [Google Scholar] [CrossRef]
  28. Balcıoğlu, E.B. Potential effects of polycyclic aromatic hydrocarbons (PAHs) in marine foods on human health: A critical review. Toxin Rev. 2016, 35, 98–105. [Google Scholar] [CrossRef]
  29. Moorthy, B.; Chu, C.; Carlin, D.J. Polycyclic aromatic hydrocarbons: From metabolism to lung cancer. Toxicol. Sci. 2015, 145, 5–15. [Google Scholar] [CrossRef] [Green Version]
  30. Lemieux, C.L.; Long, A.S.; Lambert, I.B.; Lundstedt, S.; Tysklind, M.; White, P.A. In vitro mammalian mutagenicity of complex polycyclic aromatic hydrocarbon mixtures in contaminated soils. Environ. Sci. Technol. 2015, 49, 1787–1796. [Google Scholar] [CrossRef]
  31. Yadetie, F.; Zhang, X.; Hanna, E.M.; Aranguren-Abadía, L.; Eide, M.; Blaser, N.; Brun, M.; Jonassen, I.; Goksøyr, A.; Karlsen, O.A. RNA-Seq analysis of transcriptome responses in Atlantic cod (Gadus morhua) precision-cut liver slices exposed to benzo[a]pyrene and 17-ethynylestradiol. Aquat. Toxicol. 2018, 201, 174–186. [Google Scholar] [CrossRef] [PubMed]
  32. Honda, M.; Suzuki, N. Toxicities of Polycyclic Aromatic Hydrocarbons for Aquatic Animals. Int. J. Environ. Res. Public Health 2020, 17, 1363. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Suzuki, N.; Ogiso, S.; Yachiguchi, K.; Kawabe, K.; Makino, F.; Toriba, A.; Kiyomoto, M.; Sekiguchi, T.; Tabuchi, Y.; Kondo, T.; et al. Monohydroxylated polycyclic aromatic hydrocarbons influence spicule formation in the early development of sea urchins (Hemicentrotus pulcherrimus). Comp. Biochem. Physiol. Part C 2015, 171, 55–60. [Google Scholar] [CrossRef] [PubMed]
  34. Azizi, G.; Akodad, M.; Baghour, M.; Layachi, M.; Moumen, A. The use of Mytilus spp. mussels as bioindicators of heavy metal pollution in the coastal environment. A review. J. Mater. Environ. Sci. 2018, 9, 1170–1181. [Google Scholar] [CrossRef]
  35. Hano, T.; Ito, K.; Ohkubo, N.; Sakaji, H.; Watanabe, A.; Takashima, K.; Sato, T.; Sugaya, T.; Matsuki, K.; Onduka, T.; et al. Occurrence of neonicotinoids and fipronil in estuaries and their potential risks to aquatic invertebrates. Environ. Pollut. 2019, 252, 205–215. [Google Scholar] [CrossRef]
  36. Loganathan, B.G.; Kumar, K.S.; Masunaga, S.; Sajwan, K.S. Polychlorinated dibenzo-p-dioxins, dibenzofurans, and dioxin-like polychlorinated biphenyls in sediment and mussel samples from Kentucky Lake, USA. Arch. Environ. Contam. Toxicol. 2008, 54, 20–30. [Google Scholar] [CrossRef]
  37. Meng, J.; Hong, S.; Wang, T.; Li, Q.; Yoon, S.J.; Lu, Y.; Giesy, J.P.; Khim, J.S. Traditional and new POPs in environments along the Bohai and Yellow Seas: An overview of China and South Korea. Chemosphere 2017, 169, 503–515. [Google Scholar] [CrossRef]
  38. Rodríguez-Navas, C.; Björklund, E.; Bak, S.A.; Hansen, M.; Krogh, K.A.; Maya, F.; Forteza, R.; Cerdà, V. Pollution pathways of pharmaceutical residues in the aquatic environment on the island of Mallorca, Spain. Arch. Environ. Contam. Toxicol. 2013, 65, 56–66. [Google Scholar] [CrossRef]
  39. Xue, X.; Xue, J.; Liu, W.; Adams, D.H.; Kannan, K. Trophic magnification of parabens and their metabolites in a subtropical marine food web. Environ. Sci. Technol. 2017, 51, 780–789. [Google Scholar] [CrossRef]
  40. Tanaka, H.; Onduka, T. Background levels of PAHs in the coastal waters of Japan based on residual concentrations of bivalves. J. Environ. Chem. 2010, 20, 137–148. [Google Scholar] [CrossRef] [Green Version]
  41. Tsai, M.L.; Dai, C.F. Life history plasticity and reproductive strategy enabling the invasion of Ligia exotica (Crustacea: Isopoda) from the littoral zone to an inland creek. Mar. Ecol. Prog. Ser. 2001, 210, 175–184. [Google Scholar] [CrossRef]
  42. Hurtado, L.A.; Mateos, M.; Santamaria, C.A. Phylogeography of supralittoral rocky intertidal Ligia isopods in the Pacific region from central California to central Mexico. PLoS ONE 2010, 5, e11633. [Google Scholar] [CrossRef] [PubMed]
  43. Kang, C.; Choy, E.J.; Son, Y.; Lee, J.; Kim, J.K.; Kim, Y.; Lee, K. Food web structure of a restored macroalgal bed in the eastern Korean peninsula determined by C and N stable isotope analyses. Mar. Biol. 2008, 153, 1181–1198. [Google Scholar] [CrossRef]
  44. Undap, S.L.; Matsunaga, S.; Honda, M.; Sekiguchi, T.; Suzuki, N.; Khalil, F.; Qiu, X.; Shimasaki, Y.; Ando, H.; Sato-Okoshi, W.; et al. Accumulation of organotins in wharf roach (Ligia exotica Roux) and its ability to serve as a biomonitoring species for coastal pollution. Ecotoxicol. Environ. Saf. 2013, 96, 75–79. [Google Scholar] [CrossRef] [PubMed]
  45. Qiu, X.; Undap, S.L.; Honda, M.; Sekiguchi, T.; Suzuki, N.; Shimasaki, Y.; Ando, H.; Sato-Okoshi, W.; Wada, T.; Sunobe, T.; et al. Pollution of radiocesium and radiosilver in wharf roach (Ligia sp.) by the Fukushima Dai-ichi Nuclear Power Plant accident. J. Radioanal. Nucl. Chem. 2016, 311, 121–126. [Google Scholar] [CrossRef]
  46. Honda, M.; Qiu, X.; Koyama, J.; Uno, S.; Undap, S.L.; Shimasaki, Y.; Oshima, Y. The wharf roach, Ligia sp., A novel indicator of polycyclic aromatic hydrocarbon contamination in coastal areas. Int. J. Environ. Res. 2018, 12, 1–11. [Google Scholar] [CrossRef]
  47. Niu, L.; Cai, H.; Van Gelder, P.H.; Luo, P.; Liu, F.; Yang, Q. Dynamics of polycyclic aromatic hydrocarbons (PAHs) in water column of Pearl River estuary (China): Seasonal pattern, environmental fate and source implication. Appl. Geochem. 2018, 90, 39–49. [Google Scholar] [CrossRef]
  48. Thomann, R.V.; Komlos, J. Model of biota-sediment accumulation factor for polycyclic aromatic hydrocarbons. Environ. Toxicol. Chem. Int. J. 1999, 18, 1060–1068. [Google Scholar] [CrossRef]
  49. Hayakawa, K.; Makino, F.; Yasuma, M.; Yoshida, S.; Chondo, Y.; Toriba, A.; Kameda, T.; Tang, N.; Kunugi, M.; Nakase, H.; et al. Polycyclic aromatic hydrocarbons in surface water of the Southeastern Japan Sea. Chem. Pharm. Bull. 2016, 64, 625–631. [Google Scholar] [CrossRef] [Green Version]
  50. Koudryashova, Y.; Chizhova, T.; Tishchenko, P.; Hayakawa, K. Seasonal Variability of Polycyclic Aromatic Hydrocarbons (PAHs) in a Coastal Marine Area in the Northwestern Region of the Sea of Japan/East Sea (Possiet Bay). Ocean Sci. J. 2019, 54, 635–655. [Google Scholar] [CrossRef]
  51. Zhang, H.; Sun, L.; Sun, T.; Li, H.; Luo, Q. Spatial distribution and seasonal variation of polycyclic aromatic hydrocarbons (PAHs) contaminations in surface water from the Hun River, Northeast China. Environ. Monit. Assess. 2013, 185, 1451–1462. [Google Scholar] [CrossRef]
  52. Xiao, R.; Bai, J.; Wang, J.; Lu, Q.; Zhao, Q.; Cui, B.; Liu, X. Polycyclic aromatic hydrocarbons (PAHs) in wetland soils under different land uses in a coastal estuary: Toxic levels, sources and relationships with soil organic matter and water-stable aggregates. Chemosphere 2014, 110, 8–16. [Google Scholar] [CrossRef] [PubMed]
  53. Pavoni, B.; Caliceti, M.; Sperni, L.; Sfriso, A. Organic micropollutants (PAHs, PCBs, pesticides) in seaweeds of the lagoon of Venice. Oceanol. Acta 2003, 26, 585–596. [Google Scholar] [CrossRef] [Green Version]
  54. Jiang, Y.; Lin, T.; Wu, Z.; Li, Y.; Li, Z.; Guo, Z.; Yao, X. Seasonal atmospheric deposition and air–sea gas exchange of polycyclic aromatic hydrocarbons over the Yangtze River Estuary, East China Sea: Implications for source–sink processes. Atmos. Environ. 2018, 178, 31–40. [Google Scholar] [CrossRef]
  55. Inam, E.; Offiong, N.A.; Essien, J.; Kang, S.; Kang, S.Y.; Antia, B. Polycyclic aromatic hydrocarbons loads and potential risks in freshwater ecosystem of the Ikpa River Basin, Niger Delta—Nigeria. Environ. Monit. Assess. 2016, 188, 49. [Google Scholar] [CrossRef] [PubMed]
  56. Francioni, E.; Wagener, A.L.R.; Scofield, A.L.; Depledge, M.H.; Cavalier, B. Evaluation of the mussel Perna perna as a biomonitor of polycyclic aromatic hydrocarbon (PAH) exposure and effects. Mar. Pollut. Bull. 2007, 54, 329–338. [Google Scholar] [CrossRef] [PubMed]
  57. Soclo, H.H.; Garrigues, P.; Ewald, M. Origin of polycyclic aromatic hydrocarbons (PAHs) in coastal marine sediments: Case studies in Cotonou (Benin) and Aquitaine (France) areas. Mar. Pollut. Bull. 2000, 40, 387–396. [Google Scholar] [CrossRef]
  58. Horiguchi, H.; Hironaka, M.; Meyer-Rochow, V.B.; Hariyama, T. Water uptake via two pairs of specialized legs in Ligia exotica (Crustacea, Isopoda). Biol. Bull. 2007, 213, 196–203. [Google Scholar] [CrossRef]
  59. Meador, J.P.; Casillas, E.; Sloan, C.A.; Varanasi, U. Comparative bioaccumulation of polycyclic aromatic hydrocarbons from sediment by two infaunal invertebrates. Mar. Ecol. Prog. Ser. 1995, 123, 107–124. [Google Scholar] [CrossRef] [Green Version]
  60. Hickey, C.W.; Roper, D.S.; Holland, P.T.; Trower, T.M. Accumulation of organic contaminants in two sediment-dwelling shellfish with contrasting feeding modes: Deposit-(Macomona liliana) and filter-feeding (Austrovenus stutchburyi). Arch. Environ. Contam. Toxicol. 1995, 29, 221–231. [Google Scholar] [CrossRef]
  61. Van Brummelen, T.C.; Verweij, R.A.; Wedzinga, S.A.; Van Gestel, C.A. Polycyclic aromatic hydrocarbons in earthworms and isopods from contaminated forest soils. Chemosphere 1996, 32, 315–341. [Google Scholar] [CrossRef]
  62. Van Gestel, C.A.; Loureiro, S. Terrestrial isopods as model organisms in soil ecotoxicology: A review. ZooKeys 2018, 801, 127–162. [Google Scholar] [CrossRef] [PubMed]
  63. Haritash, A.K.; Kaushik, C.P. Biodegradation aspects of polycyclic aromatic hydrocarbons (PAHs): A review. J. Hazard. Mater. 2009, 169, 1–5. [Google Scholar] [CrossRef] [PubMed]
  64. Liang, Y.; Tse, M.F.; Young, L.; Wong, M.H. Distribution patterns of polycyclic aromatic hydrocarbons (PAHs) in the sediments and fish at Mai Po Marshes Nature Reserve, Hong Kong. Water Res. 2007, 41, 1303–1311. [Google Scholar] [CrossRef]
Figure 1. Locations of the 12 sampling sites.
Figure 1. Locations of the 12 sampling sites.
Ijerph 18 00630 g001
Figure 2. The concentrations of 15 PAHs in wharf roaches (ng/g-dw) collected from 12 sampling sites.
Figure 2. The concentrations of 15 PAHs in wharf roaches (ng/g-dw) collected from 12 sampling sites.
Ijerph 18 00630 g002
Figure 3. Concentrations of 15 PAHs in mussels (ng/g-dw) collected from 12 sampling sites.
Figure 3. Concentrations of 15 PAHs in mussels (ng/g-dw) collected from 12 sampling sites.
Ijerph 18 00630 g003
Figure 4. Concentrations of 15 PAHs in seawater (ng/L) collected from 12 sampling sites.
Figure 4. Concentrations of 15 PAHs in seawater (ng/L) collected from 12 sampling sites.
Ijerph 18 00630 g004
Figure 5. Concentrations of 15 PAHs in soil and sand (ng/g-dw) collected from 12 sampling sites.
Figure 5. Concentrations of 15 PAHs in soil and sand (ng/g-dw) collected from 12 sampling sites.
Ijerph 18 00630 g005
Figure 6. Concentrations of 15 PAHs in drifting seaweed (ng/g-dw) collected from 11 sampling sites. Samples could not be found at the Shimane Prefecture site.
Figure 6. Concentrations of 15 PAHs in drifting seaweed (ng/g-dw) collected from 11 sampling sites. Samples could not be found at the Shimane Prefecture site.
Ijerph 18 00630 g006
Figure 7. PCA score plot in which the first two PCs account for 88% of the variance.
Figure 7. PCA score plot in which the first two PCs account for 88% of the variance.
Ijerph 18 00630 g007
Figure 8. Diagnostic ratios calculated for soil and sand samples and seawater samples.
Figure 8. Diagnostic ratios calculated for soil and sand samples and seawater samples.
Ijerph 18 00630 g008
Table 1. Site locations and sampling date.
Table 1. Site locations and sampling date.
Site No.City, PrefectureGPS CoordinateSampling Date
1Ajigasawa, Aomori40.779312, 140.21566012 September 2018
2Yurihonjo, Akita39.501981, 140.04514613 September 2018
3Tsuruoka, Yamagata38.722078, 139.68516313 September 2018
4Nagaoka, Niigata37.646270, 138.76476714 September 2018
5Suzu, Ishikawa37.307115, 137.23149228 September 2018
6Miyazu, Kyoto35.558870, 135.18658815 September 2018
7Kami, Hyogo35.652180, 134.60514516 September 2018
8Hamada, Shimane34.904503, 132.06074117 September 2018
9Hagi, Yamaguchi34.437040, 131.41682117 September 2018
10Fukuoka, Fukuoka33.603976, 130.27492318 September 2018
11Karatsu, Saga33.484003, 129.94434819 September 2018
12Hirado, Nagasaki33.361646, 129.62639519 September 2018
Table 2. Total concentration of 15 PAHs in samples from 12 sampling sites (ng/g-dw or ng/L).
Table 2. Total concentration of 15 PAHs in samples from 12 sampling sites (ng/g-dw or ng/L).
PrefectureWharf Roach
(ng/g-dw)
Mussel
(ng/g-dw)
Soil/Sand
(ng/g-dw)
Drifting Seaweed
(ng/g-dw)
Sea Water
(ng/L)
Aomori49.265.123.834.412.4
Akita96.044.59.9512.419.1
Yamagata53.030.137.631.07.79
Niigata85.233.420.211.89.90
Ishikawa47.913250682.412.8
Kyoto42.143.213889.515.6
Hyogo44.153.358.150.616.4
Shimane44.135.447.5na10.0
Yamaguchi26.929.22.0621.310.2
Fukuoka52.841.13.4562.917.8
Saga85.016.63.8518.711.9
Nagasaki37.738.611.110916.4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Honda, M.; Mukai, K.; Nagato, E.; Uno, S.; Oshima, Y. Correlation between Polycyclic Aromatic Hydrocarbons in Wharf Roach (Ligia spp.) and Environmental Components of the Intertidal and Supralittoral Zone along the Japanese Coast. Int. J. Environ. Res. Public Health 2021, 18, 630. https://doi.org/10.3390/ijerph18020630

AMA Style

Honda M, Mukai K, Nagato E, Uno S, Oshima Y. Correlation between Polycyclic Aromatic Hydrocarbons in Wharf Roach (Ligia spp.) and Environmental Components of the Intertidal and Supralittoral Zone along the Japanese Coast. International Journal of Environmental Research and Public Health. 2021; 18(2):630. https://doi.org/10.3390/ijerph18020630

Chicago/Turabian Style

Honda, Masato, Koki Mukai, Edward Nagato, Seiichi Uno, and Yuji Oshima. 2021. "Correlation between Polycyclic Aromatic Hydrocarbons in Wharf Roach (Ligia spp.) and Environmental Components of the Intertidal and Supralittoral Zone along the Japanese Coast" International Journal of Environmental Research and Public Health 18, no. 2: 630. https://doi.org/10.3390/ijerph18020630

APA Style

Honda, M., Mukai, K., Nagato, E., Uno, S., & Oshima, Y. (2021). Correlation between Polycyclic Aromatic Hydrocarbons in Wharf Roach (Ligia spp.) and Environmental Components of the Intertidal and Supralittoral Zone along the Japanese Coast. International Journal of Environmental Research and Public Health, 18(2), 630. https://doi.org/10.3390/ijerph18020630

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