Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Dietary homogenization and spatial distributions of carbon, nitrogen, and sulfur isotope ratios in human hair in South Korea

  • Han-Seul Lee ,

    Contributed equally to this work with: Han-Seul Lee, Ji-Yu Shim

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Software, Writing – original draft

    Affiliations Research Center for Geochronology and Isotope Analysis, Korea Basic Science Institute, Cheongju-si, Chungbuk, Republic of Korea, Graduate School of Analytical Science and Technology, Chungnam National University, Daejeon, Republic of Korea

  • Ji-Yu Shim ,

    Contributed equally to this work with: Han-Seul Lee, Ji-Yu Shim

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Validation

    Affiliations Research Center for Geochronology and Isotope Analysis, Korea Basic Science Institute, Cheongju-si, Chungbuk, Republic of Korea, Graduate School of Analytical Science and Technology, Chungnam National University, Daejeon, Republic of Korea

  • Woo-Jin Shin,

    Roles Data curation, Funding acquisition, Resources, Visualization, Writing – review & editing

    Affiliation Research Center for Geochronology and Isotope Analysis, Korea Basic Science Institute, Cheongju-si, Chungbuk, Republic of Korea

  • Seung-Hyun Choi,

    Roles Formal analysis, Investigation

    Affiliation Center for Research Equipment, Korea Basic Science Institute, Cheongju-si, Chungbuk, Republic of Korea

  • Yeon-Sik Bong,

    Roles Formal analysis

    Affiliation Center for Research Equipment, Korea Basic Science Institute, Cheongju-si, Chungbuk, Republic of Korea

  • Kwang-Sik Lee

    Roles Conceptualization, Funding acquisition, Methodology, Supervision, Writing – original draft, Writing – review & editing

    kslee@kbsi.re.kr

    Affiliations Research Center for Geochronology and Isotope Analysis, Korea Basic Science Institute, Cheongju-si, Chungbuk, Republic of Korea, Graduate School of Analytical Science and Technology, Chungnam National University, Daejeon, Republic of Korea

Abstract

Dietary homogenization has progressed worldwide due to westernization and the globalization of food production systems. We investigated dietary heterogeneity in South Korea by examining the spatial distribution of carbon (C), nitrogen (N), and sulfur (S) isotope ratios using 264 human hair samples. Overall, variation in isotope values was small, indicating low dietary heterogeneity. We detected differences in δ13C, δ15N, and δ34S values between administrative provinces and metropolitan cities; inter-regional differences were typically < 1 ‰. Values of δ34S were significantly lower in hair samples from inland regions relative to those from coastal locations, and a similar pattern was observed in δ15N values. Understanding geographic variation in δ34S and δ15N values in human hair is useful for provenancing humans in South Korea.

Introduction

Stable isotope analyses of bio-elements, including carbon (C), hydrogen (H), nitrogen (N), oxygen (O), and sulfur (S), are commonly used to investigate the dietary protein sources consumed by both contemporary and past human populations [14]. Dietary intake patterns are recorded in body tissues such as hair, nail, bones, and teeth, which allow for the reconstruction of lifestyle attributes related to geographic origin and dietary habits [57]. In particular, human hair is a continually growing, readily sampled tissue that can serve as a longitudinal record of movement history [811]. In contrast to other body tissues, hair has no metabolic activity after being produced, and diet characteristics, region-of-origin, and recent travel history are thus stored in this tissue [1214].

Over the past two decades, scientists have investigated the spatial distribution of bio-element stable isotopes in human tissues for forensic analyses [2, 10]. Among stable isotopes, H and O isotope ratios in human hair are closely related to drinking water source [5]. Regional-scale distributions of H and O isotope ratios in modern human hair were documented in the United States (USA) by Ehleringer et al. [5] and showed geographically predictable variations. Isoscapes, visualizations of spatial distributions of isotopic signatures, are a valuable tool for mapping heterogeneity against various environmental and geographic parameters and provide useful baseline information for provenance and forensic inquiries. Isoscape studies focused on H and O values in water and hair have been conducted worldwide (e.g., [5, 15]). Recently, Gautam et al. [16, 17] demonstrated that H and O isotope signatures in human hair can be used in provenance- and forensic-related activities, such as residential history and migration history in South Korea, based on a clear latitudinal gradient in isotope compositions. Their study demonstrated the predictive ability of H and O isoscapes to determine the regional origin of hair in South Korea.

Unlike H and O isotopes, the stable isotope ratios of C, N, and S are typically used for dietary reconstruction studies because their ratios in human tissues are entirely related to food consumption [2]. Therefore, stable isotopes of C, N, and S in human tissues can complement H and O isotope analyses [2]. Hülsemann et al. [18] provided an excellent global overview of C and N isotope patterns using approximately 4,000 individual human hair samples and reported significant differences in δ13C and δ15N values among countries or regions. Broadly, human dietary patterns have homogenized due to rapid globalization in food production and distribution processes, resulting in reduced differences in isotope values among regions [3]. However, dietary heterogeneity is still preserved among European countries [6], North American countries [2, 19], and some regions in Japan [3].

The global spatial distribution of δ13C values in human hair is strongly related to the amount of C4 plants in human diets [18, 20] and δ15N values are reflective of trophic level, and can thus indicate whether an individual is omnivorous or vegetarian if they consume only terrestrial foods [4]. Individuals who consume large amounts of fish and seafood are expected to have higher δ15N values than those who mainly eat terrestrial foods [21]. Generally, δ34S values are reflective of the consumption of foods from different geological regions. Specifically, marine-derived foods exhibit approximately +20 ‰ of δ34S values [2]. Thus, δ34S values can reflect both diet and geographic origin [19, 22, 23]. Factors potentially explaining differences in δ34S values include bedrock geochemistry [11, 22], atmospheric deposition of sea-spray salts on coastal soils and vegetation [19, 24], and the consumption of fish- or meat-derived proteins [2, 15]. Valenzuela et al. [2] demonstrated distinct geographic pattern in δ34S values in the USA using human hair samples, where inland values were much lower than those from coastal regions. These regional differences can then be applied in forensic provenancing studies [2, 18, 19].

Datasets of stable isotope values determined from human hair are highly biased to North America and Europe, with limited data available from other countries [18]. Minimal data are available from East Asia, with only limited information about δ13C and δ15N values in human hair available from South Korea. To our knowledge, no study has examined spatial variation and heterogeneity in C, N, and S isotopes in human hair at a national scale in South Korea.

The use of multi-isotopes from human tissues improves the predictive power of isoscapes [22]. Building from the recently published H and O isoscapes created using human hair samples from South Korea [16, 17], we investigated spatial distributions of δ13C, δ15N, and δ34S values in human hair samples to extend the applicability of isotopes for human provenancing using a multi-isotope approach. We further investigated dietary heterogeneity and characterized the dietary patterns of contemporary Koreans using C, N, and S isotopes. Finally, we compared isotope ratios from the Korean hair samples with those from other countries to characterize dietary protein sources in Korea in a global context.

Materials and methods

Ethical statement and data availability

This study was approved by the Bioethics Committee at the Korea Basic Science Institute (IRB number: 005). Hair samples were collected from trash bins with oral consent of the owners of barbershops. Only samples identified as male hair were collected from barbershops. Biological sex was the only donor information collected. Isotope data and sample locations are available in S1 Table and S2 Table.

Sample collection and preparation

We collected 264 human scalp hair samples from trash bins at barbershops for men in South Korea. Hair samples were placed in individual plastic bags at the time of collection. To ensure that the isotopic signature of hair samples reflected the geographic location of sampling regions, hair samples were obtained from barbershops located in residential areas. We assumed that the collected hair samples were representative of men residing in the area around the collection site.

Prior to analysis, to remove lipids and surface contaminants, approximately 3 g of hair from a bundle of multiple hairs from the same individual were placed in a 100-mL beaker, washed using 0.5% Triton-X (Sigma Aldrich, St Louis, MO, USA), Milli-Q water (> 18 MΩ cm), and acetone (Honeywell Burdick & Jackson, Muskegon, MI, USA), and then rinsed twice with Milli-Q water. The amount of solvent used for each wash was sufficient to completely submerge the hair samples, and each step was treated in an ultrasonic bath for 2 min. Washed samples were then fully dried in a drying oven at 60°C overnight. Each dried sample was pulverized into a fine powder using a ball mill (MM400; Retsch, Haan, Germany). Homogenized samples were stored in vials prior to analysis. For the isotope analysis, 34 food samples were collected from supermarkets in metropolitan cities, subsequently, lipids were removed and samples were freeze-dried at −80°C for 72 h and then ground using a ceramic ball mill.

Stable isotope analyses

C, N, and S isotope ratios were determined using an isotope ratio mass spectrometer (IRMS; Isoprime VisION, Isoprime Ltd., Manchester, UK) that was equipped with a vario PyroCube elemental analyzer (Elementar, Hesse, Germany) that included a combustion tube maintained at 1150°C and a U-shaped adsorption column to separate the generated CO2 gas; the analyzer was operated in continuous flow mode. For the C isotope ratio, approximately 0.1 mg of prepared hair sample was encapsulated in a tin container (3.5 × 1.5 mm; Elemental Microanalysis, Okehampton, UK). For the N and S isotope ratios, approximately 0.8 mg of prepared hair sample was encapsulated in tin capsules, respectively. The tin capsules were stored in a desiccator until analysis to avoid vapor absorption. The encapsulated sample was combusted at 1150°C in a combustion reactor and excess oxides were removed in a reduction reactor at 850°C. Thereafter, N2, CO2, and SO2 gases that had been generated were separated using chromatography and the isotope ratios were measured using the IRMS. Because of varying content of C, N, and S in hair samples, C, N, and S isotopes were analyzed separately.

Stable isotope ratios were reported using the standard delta notation (δ) relative to an international standard unit per mill (‰) as follows: δ (‰) = (Rsample/Rreference − 1) × 1000, where Rsample and Rreference are the molar ratios of the heavy to light isotopes of the sample and standard, respectively, representing 13C/12C, 15N/14N, or 34S/32S. The resulting δ13C, δ15N, and δ34S values were reported against Vienna Pee Dee Belemnite (VPDB), atmospheric nitrogen (Air), and Vienna Canyon Diablo Troilite (VCDT), respectively. The analytical precision was within ± 0.1 ‰ for C, ± 0.2 ‰ for N, and ± 0.2 ‰ for S.

The following standard reference materials were used for the calibration of C, N, and S isotope ratios: NBS-22 (oil, δ13CVPDB = −29.8 ‰), IAEA-CH-3 (cellulose, δ13CVPDB = −24.73 ‰), IAEA-CH-6 (sucrose, δ13CVPDB = −10.45 ‰), IAEA-600 (caffeine, δ 13CVPDB = −27.77 ‰, δ15NAir = +1.0 ‰), USGS-40 (l-glutamic acid, δ 13CVPDB = −26.39 ‰, δ15NAir = −4.52 ‰), USGS-42 (Tibetan human hair, δ13CVPDB = −21.09 ‰, δ15NAir = +8.05 ‰, δ34SVCDT = +7.84 ‰), USGS-43 (Indian human hair, δ13CVPDB = −21.28 ‰, δ15NAir = +8.44 ‰, δ34SVCDT = +10.46 ‰), IAEA-N-2 (ammonium sulfate, δ15NAir = +20.3 ‰), IAEA-SO-5 (barium sulfate, δ34SVCDT = +0.5 ‰), and NBS-127 (barium sulfate, δ34SVCDT = +21.1 ‰).

Geostatistical analyses

Normality was tested using the Shapiro-Wilk test. P-values were considered significant at the α = 0.05 level. We used Student’s t-tests to compare mean isotope values between two groups and used analysis of variance (ANOVA) to compare isotope values among provinces. Statistical significance was interpreted based on p < 0.05. These analyses were conducted using XLSTAT Pro software (Addinsoft, New York, NY, USA).

The spatial distributions of the δ13C, δ15N, and δ34S values obtained from the hair samples were visualized across South Korea (including the province of Jeju Island) using interpolation mapping. Maps were produced using ordinary kriging (OK), a common geostatistical interpolation method, and we did not consider the application of other interpolation methods (e.g., empirical Bayesian kriging, simple kriging). The OK was applied using ArcGIS version 10.8 (ESRI, Redlands, CA, USA). In addition, spatial autocorrelation (Moran’s index) was performed using the software to examine the spatial distribution pattern of human hair (clustered, dispersed, or random).

Results and discussion

C, N, and S isotope ratios in human hair and food

The range and mean values of C, N, and S isotope ratios obtained from hair samples collected from nine administrative provinces and seven metropolitan cities in South Korea are summarized in Table 1.

thumbnail
Table 1. Carbon, nitrogen and sulfur isotope values of human hair samples collected from nine administrative provinces (GG through JJ) and major Metropolitan Cities (MC) in South Korea.

https://doi.org/10.1371/journal.pone.0256404.t001

Overall, the mean δ13C value among all hair samples was −19.2 ± 0.5 ‰ (range: −20.9 ‰ to −17.6 ‰). The mean δ15N value was +9.7 ± 0.5 ‰ (range: +8.4 to +11.7 ‰). The mean δ34S value was +7.3 ± 0.8 ‰ (range: +3.4 to +10.0 ‰). The δ13C, δ15N, and δ34S values obtained from food from local stores and supermarkets are shown in Table 2 [2532].

thumbnail
Table 2. Carbon, nitrogen and sulfur isotope values of food samples collected from South Korea.

https://doi.org/10.1371/journal.pone.0256404.t002

Frequency distributions of C, N, and S isotope ratios in both hair and food samples are shown in Fig 1. The δ13C values in food samples ranged from −29.1 ‰ to −9.9 ‰ (n = 545), indicative of both C3 and C4 plants. The majority of the food samples were reflective of C3 plants (Fig 1A). Generally, the results indicated that Koreans consume a wide variety of foods, inclusive of C3 and C4 plants [3337], because the δ13C values from hair samples corresponded to the mid-range value between C3 and C4 plants. The typical δ13C ranges of C3 plants and C4 plants are -33 ‰ to -23 ‰ and -16 ‰ to -9 ‰, respectively [36]. According to the Korea Rural Economic Institute [38], rice, wheat, and corn accounted for 55.4%, 22.7%, and 18.9% of all grains consumed in South Korea, respectively, in 2003–2018. Rice and wheat are C3 plants, while corn is a C4 plant; thus, the δ13C values obtained from hair samples are at least partially reflective of the relative abundance of these foods in an individual’s diet.

thumbnail
Fig 1.

Frequency distributions of (A) δ13C, (B) δ15N, and (C) δ34S values in human hair and food samples collected across South Korea.

https://doi.org/10.1371/journal.pone.0256404.g001

Unlike δ13C values, the distribution of δ15N values differed significantly among food and hair samples (Fig 1B), where those from food samples were significantly lower. The δ15N values from food samples ranged from −6.2 ‰ to +16.8 ‰ (mean +4.3 ± 4.0 ‰, n = 481), with the exception of one sample greatly enriched in 15N (+23.0 ‰), and approximately 60% of the food samples occurred within the range of +2 to +8 ‰. Thus, the average δ15N value in human hair samples was +5.4 ‰ higher than that from food samples. An elevated δ15N value in hair is strongly related to animal protein consumption and is reflective of individual dietary patterns. The supply of animal-based protein (approximately 56% of the total dietary protein) is significantly higher than plant-based protein (approximately 44% of total protein) in South Korea [39]. Thus, animal protein consumption is likely responsible for high δ15N values in Korean hair samples.

The δ34S values from both food and hair samples showed a similar pattern to that of δ15N values, and seafood samples with relatively higher δ34S values were more clearly separated from human hair and other food samples (Fig 1C). Generally, high consumption of marine-based and animal proteins is reflected by high δ15N and δ34S values in hair [2], and high δ34S values are associated with high intake of marine-based protein or the sea spray effect (Fig 1C; [19, 40]). We visualized the spatial distributions of δ15N and δ34S values of rice (unpublished data) and garlic [41] samples collected across South Korea, and little or no differences were observed in δ15N and δ34S values between coastal and inland crops. Such spatial distribution may indicate little or no sea-spray effects around coastal areas in the mainland. However, Jeju Island (JJ), located in the southern part of the Korean Peninsula is the administrative province most affected by sea-spray. Jeju Island samples were collected close to the coast, so local foods may have contained particularly high amounts of marine aerosols. Garlic samples collected from Jeju Island had slightly higher δ34S values than those collected from the mainland, indicating a weak sea-spray effect. Meanwhile, δ15N values of pork meat samples produced on Jeju Island did not differ from those produced on the mainland (p = 0.23, n = 16) [29]. No δ34S data for animal meats such as pork, beef and chicken are currently available for South Korea. Therefore, given the δ15N and δ 34S data of agro-products available for South Korea, high values of δ34S and δ15N in coastal human hairs from the mainland can primarily be interpreted as a result of fish and seafood consumption [19].

The results of the isotope analyses of hair and food samples are shown in Fig 2. Broadly, the δ13C and δ15N values from food samples were scattered across a wide range, whereas those from hair samples were confined to a narrow range. This may be because hair reflects the average isotope values of various food samples.

thumbnail
Fig 2.

Relationships between (A) δ13C and δ15N, and (B) δ15N and δ34S values obtained from human hair and food samples collected across South Korea.

https://doi.org/10.1371/journal.pone.0256404.g002

Samples of C3 plants (apples, onions, and garlic) and tofu made from soybean (C3 plant) had low δ13C values (Table 2); δ13C values in livestock products were higher and those from milk samples were similar to C3 plants. The δ13C values in marine fish and seaweeds ranged from −24.2 ‰ to −14.3 ‰ and −25.1 ‰ to −16.7 ‰, respectively. The δ15N values of marine fish were high relative to C3 plants, livestock products, and seaweeds, excluding mushroom samples, which have exceptionally high δ15N values due to the use of organic fertilizers such as animal manure [35]. In South Korea, A. bisporus mushrooms are usually cultivated with a compost medium in a greenhouse system, and not with soil-based cultivation [35]. The high δ15N value observed for A. bisporus may be associated with the δ15N values of the animal manure in the medium used for mushroom cultivation, which could reflect the N isotope fractionation that occurs between the medium and the mushroom fruit body during mushroom cultivation [35]. In addition, δ15N values for hair samples exhibited a similar range to those of marine fish, which were markedly higher than values in other samples. This finding indicates that humans are top-level predators in food chains, and marine fishes, which are predators in more complex food chains, can exhibit δ15N values similar to human hair.

Isotope ratios in livestock are reflective of their diets. Australian and New Zealand cattle are typically fed C3 plants, resulting in low δ13C values, whereas cattle in the USA and Mexico are typically fed corn, resulting in high δ13C values. In Korea, C3 and C4 plants are used to feed cattle; thus, δ13C values are intermediate to beef from the other continents [25, 26]. In our study, beef had a δ15N value that was approximately +3 ‰ higher than that of pork or chicken (Fig 2A) and the δ34S values in fish and seaweed were much higher than those in C3 plants and livestock products.

Regional homogeneity and localization of isotope ratios

We found no significant difference in the intake rate of animal and marine foods between metropolitan cities and administrative provinces in South Korea. On average, animal and marine food intake is reportedly only 1.2% higher in metropolitan cities in South Korea [37]. Data regarding the intake rates of animal and marine products were not available for all administrative provinces considered here. Although we expected minimal variation in δ13C, δ15N, and δ34S values in human hair from cities and provinces, we found significant differences among some provinces (Table 1) and small standard deviations within provinces (± 0.6 ‰ or less for both δ13C and δ15N). Comparisons of isotope values for all provinces and cities are provided in Fig 3.

thumbnail
Fig 3.

Boxplots showing (A) δ13C, (B) δ15N, and (C) δ34S values in human hair samples collected from nine administrative provinces (GG through JJ) and major metropolitan cities (MC) in Korea. Provinces are abbreviated as: Gyeonggi (GG), Gangwon (GW), Chungbuk (CB), Chungnam (CN), Jeonbuk (JB), Jeonnam (JN), Gyeongbuk (GB), Gyeongnam (GN), and Jeju (JJ). Seven metropolitan cities (hatched areas) are also shown. Dashed lines represent the mean isotope values of all samples. Colored bars represent the standard deviations of isotope values for all samples.

https://doi.org/10.1371/journal.pone.0256404.g003

Among all provinces, mean δ13C values were slightly higher in Jeju Island samples (−18.9 ± 0.4 ‰) and lowest in Jeonbuk (−19.6 ± 0.6 ‰). Mean δ15N was highest in Jeonnam (+10.0 ± 0.6 ‰) and lowest in Chungbuk (+9.3 ± 0.3 ‰), the only inland province that is not in contact with the sea. The mean δ34S value was highest in Jeju Island (+7.8 ± 0.7 ‰) and lowest in Chungbuk (+6.8 ± 0.3 ‰). Intra-region variation in δ34S was highest in Gangwon (± 1.2 ‰). Overall, the variation in δ34S values was low and did not deviate significantly from the mean.

Dependency on imported food, rather than food sovereignty, is increasing annually in South Korea [39], which would be expected to decrease regional differences in isotope ratios. Dietary patterns also tend to change and homogenize with industrialization and economic development. Consequently, food products in Korean grocery stores reflect global supermarkets due to changes in consumption patterns driven by the expansion of food markets and services. This trend toward globalization is likely driving the relatively homogenous pattern we observed in hair isotope values (Fig 3). According to Korea’s Food Balance Sheet of 2018 [39], Korea produced only 21.9% of the grain its population consumed, with the remaining 78.1% imported from the USA, Australia, and Brazil. Excluding rice, of which Korea produced 82.5% of its own amount consumed, the bulk of the major agricultural products on today’s shelves are imported. In addition, the self-sufficiency rate of animal food items, excluding eggs, has sharply declined.

Comparison of isotope values in hair among countries

To characterize the dietary patterns of contemporary Koreans on a global scale, we compared the isotope ratios obtained from the hair samples with those from other countries (Fig 4). We focused on comparisons with findings from East Asia (South Korea, Japan, and China) (Fig 4).

thumbnail
Fig 4. Plots of δ13C and δ15N values of human hair.

(A) δ13C and δ15N values found in human hair in South Korea, China, and Japan and (B) the same data including observations from Mongolia, Oceania, USA, Canada and Europe for comparison. In the figure on the left, symbols without error bars represent the average values of individual studies in each country (circle: South Korea, square: Japan, triangle: China). Data sources: South Korea [3, 42], Japan [3, 4244], China [18, 44, 45], Mongolia [3, 42], Oceania [18, 44, 4648], USA [2, 4, 6, 18, 42, 49, 50], Canada [18, 19, 44], Europe [6, 8, 18, 44, 51, 52].

https://doi.org/10.1371/journal.pone.0256404.g004

The δ13C and δ15N values in Korean hair obtained in this study are similar to those reported previously [3, 42]; thus, it is assumed that isotopic variation in hair samples is not large. In particular, δ13C and δ15N values in Korean hair were similar to those in Japanese hair [3, 4245]. Given the δ13C values, the consumption rates of C3 and C4 plants are similar between Japan and Korea (Fig 4A).

For comparison, the USA, which mainly consumes C4 plants reflected by high δ13C values [18], and the European Union (EU), which mainly consumes C3 plants reflected by low δ13C values [4], are also shown (Fig 4B). Canada [18, 19, 44] and Oceania [18, 44, 4648] are also illustrated (Fig 4B). The δ13C values of Korean hair were intermediate to those of the USA (+ Canada) [2, 4, 6, 18, 19, 42, 44, 49, 50] and EU (+ Oceania) [6, 8, 18, 44, 51, 52], suggesting proportional relationships with the consumption of C3 and C4 plants and seafood. Based on both δ13C and δ15N values, male Koreans eat both C3 and C4 plants and are omnivorous.

According to the Food and Agriculture Organization of the United Nations [53], Korean people consumed significantly less fish and animal protein than Japanese people on a gram per day basis decades ago. However, this gap has narrowed since 2004 and fish and animal protein intake between these two countries has been similar in recent years (S1 Fig; [53]). Thus, δ15N values obtained from contemporary Korean and Japanese hair samples should be similar, as was observed here.

The δ13C and δ15N values were higher in Korean hair relative to Chinese samples (Fig 4A). This may be because Chinese people are more dependent on plant protein intake and thus have lower δ15N values than Korean and Japanese people (S1 Fig; [53]). The low δ15N values of Chinese hair are also due to consumption of legumes or poultry with soy protein [22]. We found a positive correlation in δ13C and δ15N values among China, Japan, and South Korea (δ15N = 1.27*δ13C + 33.98, R2 = 0.996, p < 0.05). However, data from Mongolia deviated considerably from the correlation line (Fig 4B). The δ13C values of Mongolian hair were significantly lower than those in Japanese and Korean hair, whereas δ15N values were similar or slightly elevated [3, 45]. This may be reflective of a distinct lifestyle and dietary habit among Mongolians, who prefer a diet with a high intake of animal protein and little fish protein (S2 Fig; [5355]).

The range of variation of δ13C and δ15N values in Korean and Japanese hair samples were smaller than those of Chinese hair samples. This is likely because foods consumed by Korean and Japanese people are characteristic of a global supermarket, with more homogenous isotope compositions [3]. By contrast, subtle differences in food types in China and Mongolia can be identified by city [45], indicating that diets including local foods have been maintained.

Spatial distributions of C, N, and S isotopes

The spatial distributions of δ13C, δ15N, and δ34S values in human hair samples for South Korea are shown in Fig 5. A weak but interesting pattern was observed in the C isoscape (Fig 5A). In Seoul, the largest metropolitan city in South Korea, and its surrounding areas, δ13C values were slightly higher than in other regions. This was a consistent, albeit weak, pattern that was reflected in other metropolitan cities. Jeju Island, which sees high international traffic, also had higher δ13C values. Globalization has had the most rapid advancement in metropolitan cities and on Jeju Island within South Korea. Westernization and globalization tend to relate to increased consumption of foods containing corn [2, 4]. Corn, a C4 plant, relates to high δ13C values [56, 57]. Consequently, consuming corn-fed animals also increases the δ13C values in human hair. Excluding the areas mentioned above, we saw little effect of geographic and environmental factors on δ13C values (Fig 5A).

thumbnail
Fig 5.

Isoscapes of (A) δ13C, (B) δ15N, and (C) δ34S in South Korea, with sampling locations indicated by open circles and seven metropolitan cities shown by bold lines. Maps were produced using ordinary kriging (OK), a common geostatistical interpolation method using ArcGIS version 10.8 (ESRI, Redlands, CA, USA). The method used in this classification is natural break (Jenks). The administrative boundaries of nine administrative provinces are shown. Provinces are abbreviated as: Gyeonggi (GG), Gangwon (GW), Chungbuk (CB), Chungnam (CN), Jeonbuk (JB), Jeonnam (JN), Gyeongbuk (GB), Gyeongnam (GN), and Jeju (JJ).

https://doi.org/10.1371/journal.pone.0256404.g005

By contrast, there was a distinct geographic pattern in δ15N values. Specifically, N isotope values are markedly low in the northern to western inland and high along the southern coast (including the southeastern coast), that is, distinctively low in mountainous areas and high in agricultural areas (Fig 5B). Organic fertilizers such as animal manures (chicken and farmyard manures; δ15N = +14.7 ± 7.7 ‰, n = 30) are primarily used in agricultural lands in South Korea [36, 58]; thus, the spatial variation of δ15N values seems to be related to the consumption of food produced in these areas. Generally, δ15N values are also determined by the consumption of fish and meat [1, 56, 57, 59, 60] and the high values along the coastal regions are unsurprising, as marine-derived proteins have high δ15N values [61, 62].

Finally, δ34S values of human scalp hair in South Korea showed a distinct geographic pattern similar to that of δ15N values (Fig 5C), significantly lower in inland relative to coastal regions. Studies have shown that proximity to the sea can have a substantial influence on N and S isotope values, with higher δ15N and δ34S values in people who consume a large amount of marine products [22, 24]. A weak correlation was detected between δ13C and δ15N values (R2 = 0.089, p < 0.001, n = 264, Fig 6) and between δ13C and δ34S values (R2 = 0.128, p < 0.001, n = 264) from human hair samples. However, a stronger correlation was observed between δ15N and δ34S values (R2 = 0.443, p < 0.001, n = 264).

thumbnail
Fig 6.

Correlation plots showing the relationships between (A) δ13C and δ15N, (B) δ13C and δ34S, and (C) δ15N and δ34S values in human hair samples collected from barbershops across South Korea. A stronger correlation between δ15N and δ34S values was observed.

https://doi.org/10.1371/journal.pone.0256404.g006

Spatial autocorrelation was employed to evaluate the distributions of δ15N and δ34S values for human hair. For the δ15N and δ34S values, the Moran’s I index was 0.356 (z = 1.328, p = 0.184) and 0.141 (z = 0.538, p = 0.591), respectively. Data lacking significant autocorrelations suggest that human hairs in terms of δ15N and δ34S values are not spatially autocorrelated but are instead randomly dispersed. This observation can be explained by the possibility that the isotopic compositions of human hair could reflect both geographical and topographical features. Likewise, agricultural practices more extensively occurred in the lower elevation areas than in mountainous areas. The distribution of lowlands and highlands is not consistent with distance from marine environments. Thus, the isotopic value in each sampling region would be compensated or integrated, statistically resulting in the absence of spatial autocorrelation. From this perspective, the randomly dispersed- spatial distribution would also be attributed to westernization and globalization. Nonetheless, the interpolated map illustrates a distribution pattern that reflects land use and distance from marine environments.

Given the patterns we observed among all isotope values, we conclude that dietary patterns have become homogenized in Korea due to westernization and globalization, irrespective of region, but that dietary heterogeneity still exists in some provinces. Our results indicate that C, N, and S isotopic signatures in human hair are strongly correlated with the food consumed by individuals. These regional differences in mean δ13C, δ15N, and δ34S values are useful for provenancing studies and allow us a more holistic view of geographic patterns following the publication of H and O isotope data.

Conclusions

We collected 264 human scalp hair samples across South Korea to investigate dietary heterogeneity and spatial variation in C, N, and S isotopes at a national scale. To our knowledge, this is the first national-scale study assessing C, N, and S isotope ratios in human hair in South Korea. The mean values of δ13C, δ15N, and δ34S in hair samples were −19.2 ± 0.5 ‰, +9.7 ± 0.5 ‰, and +7.3 ± 0.8 ‰, respectively. The inter-regional difference in C isotope ratios was < 0.5 ‰, indicating dietary homogeneity due to rapid westernization and globalization. Nevertheless, regional dietary heterogeneity still exists in some areas in South Korea.

Spatial products like isoscapes are valuable tools because they allow us to visualize spatial heterogeneity in isotopes given environmental parameters, such as distance to sea. We found that δ34S values in human hair were greatly influenced by proximity to marine environments in South Korea, with higher values in coastal relative to inland regions. Interestingly, the N isoscape showed a similar pattern to that of the S isoscape. The consistency between the two illustrates the potential of isoscapes to act as a proxy for food source in South Korea, and can provide useful data for forensics and migration tracking. Further, understanding C, N, and S isoscapes in South Korea contributes to the growing literature on the identification of residential history using isotopic analyses of human hair.

Supporting information

S1 Table. δ13C, δ15N, and δ34S values of human hair samples (n = 264) collected from South Korea.

https://doi.org/10.1371/journal.pone.0256404.s001

(XLSX)

S2 Table. δ13C, δ15N, and δ34S values of food samples (n = 34) collected from metropolitan cities supermarkets in South Korea.

https://doi.org/10.1371/journal.pone.0256404.s002

(XLSX)

S1 Fig.

(A) Fish and (B) meat per capita supply for East Asian countries since 1961.

https://doi.org/10.1371/journal.pone.0256404.s003

(TIF)

S2 Fig.

Fish and animal contribution to protein supply for East Asian countries since 1961: (A) fish/total proteins, (B) animal/total proteins and (C) (fish+animal)/total proteins.

https://doi.org/10.1371/journal.pone.0256404.s004

(TIF)

References

  1. 1. O’Connell TC, Hedges RE. Investigations into the effect of diet on modern human hair isotopic values. Am J Phys Anthropol. 1999;108(4):409–25. https://doi.org/10.1002/(SICI)1096-8644(199904)108:4<409::AID-AJPA3>3.0.CO;2-E pmid:10229386
  2. 2. Valenzuela LO, Chesson LA, O’Gray SP, Cerling TE, Ehleringer JR. Spatial distributions of carbon, nitrogen and sulfur isotope ratios in human hair across the central United States. Rapid Commun Mass Spectrom. 2011;25(7):861–8. pmid:21416522
  3. 3. Kusaka S, Ishimara E, Hyodo F, Gakuhari T, Yonega M, Yumoto T, et al. Homogeneous diet of contemporary Japanese inferred from stable isotope ratios of hair. Sci Rep. 2016;6(1):33122. pmid:27616586
  4. 4. Ehleringer JR, Avalos SC, Tipple BJ, Valenzuela LO, Cerling TE. Stable isotopes in hair reveal dietary protein sources with links to socioeconomic status and health. Proc Natl Acad Sci U S A. 2020;117(33):20044–51. pmid:32747534
  5. 5. Ehleringer JR, Bowen GJ, Chesson LA, West AG, Podlesak DW, Cerling TE. Hydrogen and oxygen isotope ratios in human hair are related to geography. Proc Natl Acad Sci U S A. 2008;105(8):2788–93. pmid:18299562
  6. 6. Valenzuela LO, Chesson LA, Bowen G, Cerling TE, Ehleringer JR. Dietary heterogeneity among western industrialized countries reflected in the stable isotope ratios of human hair. PLoS One. 2012;7(3):e34234. pmid:22479574
  7. 7. Jackson GP, An Y, Konstantynova KI, Rashaid AH. Biometrics from the carbon isotope ratio analysis of amino acids in human hair. Sci Justice. 2015;55(1):43–50. pmid:25577006
  8. 8. Bol R, Marsh J, Heaton TH. Multiple stable isotope (18O, 13C, 15N and 34S) analysis of human hair to identify the recent migrants in a rural community in SW England. Rapid Commun Mass Spectrom. 2007;21(18):2951–4. pmid:17705338
  9. 9. Font L, Van Der Peijl G, Van Wetten I, Vroon P, Van Der Wagt B, Davies G. Strontium and lead isotope ratios in human hair: investigating a potential tool for determining recent human geographical movements. J Anal At Spectrom. 2012;27:719–32. https://doi.org/10.1039/c2ja10361c
  10. 10. Cerling TE, Barnette JE, Bowen GJ, Chesson LA, Ehleringer JR, Remien CH, et al. Forensic stable isotope biogeochemistry. Annu Rev Earth Planet Sci. 2016;44:175–206. https://doi.org/10.1146/annurev-earth-060115-012303
  11. 11. Kabalika Z, Morrison TA, McGill RAR, Munishi LK, Ekwem D, Mahene WL, et al. Tracking animal movements using biomarkers in tail hairs: a novel approach for animal geolocating from sulfur isoscapes. Mov Ecol. 2020;8:37. pmid:32968486
  12. 12. Petzke KJ, Boeing H, Metges CC. Choice of dietary protein of vegetarians and omnivores is reflected in their hair protein 13C and 15N abundance. Rapid Commun Mass Spectrom. 2005;19(11):1392–400. pmid:15880664
  13. 13. Hatch KA, Crawford MA, Kunz AW, Thomsen SR, Eggett DL, Nelson ST, et al. An objective means of diagnosing anorexia nervosa and bulimia nervosa using 15N/14N and 13C/12C ratios in hair. Rapid Commun Mass Spectrom. 2006;20(22):3367–73. pmid:17044129
  14. 14. Correia MA, Foley R, O’Connell TC, Ramírez‐Rozzi F, Lahr MM. Carbon and nitrogen isotopic signatures of hair, nail, and breath from tropical African human populations. Rapid Commun Mass Spectrom. 2019;33(22):1761–73. pmid:31287915
  15. 15. Valenzuela LO, Chesson LA, Bowen G, Cerling TE, Ehleringer JR. Spatial distribution of stable isotope values of human hair: Tools for region-of‐origin and travel history assignment. In: Forensic Science and Humanitarian Action: Interacting with the Dead and the Living. 1st ed. Parra RC, Zapico SC, Ubelaker DH, editors. John Wiley & Sons Ltd; 2020. pp. 385–410.
  16. 16. Gautam MK, Song BY, Shin WJ, Bong YS, Lee KS. Spatial variations in oxygen and hydrogen isotopes in waters and human hair across South Korea. Sci Total Environ. 2020;726:138365. pmid:32320869
  17. 17. Gautam MK, Song BY, Shin WJ, Bong YS, Lee KS. Datasets for spatial variation of O and H isotopes in waters and hair across South Korea. Data in Brief. 2020;30:105666. pmid:32426436
  18. 18. Hülsemann F, Lehn C, Schneiders S, Jackson G, Hill S, Rossmann A, et al. Global spatial distributions of nitrogen and carbon stable isotope ratios of modern human hair. Rapid Commun Mass Spectrom. 2015;29(22):2111–21. pmid:26467223
  19. 19. Bataille CP, Chartrand MM, Raposo F, St-Jean G. Assessing geographic controls of hair isotopic variability in human populations: A case-study in Canada. PLoS ONE. 2020;15(8):e0237105. pmid:32776947
  20. 20. Katzenberg MA, Krouse HR. Application of stable isotope variation in human tissues to problems in identification. Can Soc Forens Sci J, 1989;22(1):7–19. http://dx.doi.org/10.1080/00085030.1989.10757414
  21. 21. Mützel E, Lehn C, Peschel O, Hölzl S, Roßmann S. Assignment of unknown persons to their geographical origin by determination of stable isotopes in hair samples. Int J Legal Med. 2009;123(1):35–40. pmid:18813939
  22. 22. Ammer S, Bartelink E, Vollner J, Anderson B, Cunha E. Socioeconomic and geographic implications from carbon, nitrogen, and sulfur isotope ratios in human hair from Mexico. Forensic Sci Int. 2020;316:110455. pmid:32858379
  23. 23. Krouse HR, Levinson AA, Piggott D, Ueda A. Further stable isotope investigations of human urinary stones: comparison with other body components. Appl Geochem. 1987;2:205–11.
  24. 24. Zazzo A, Monahan FJ, Moloney AP, Green S, Schmidt O. Sulphur isotopes in animal hair track distance to sea. Rapid Commun Mass Spectrom. 2011;25(17):2371–8. pmid:21818798
  25. 25. Bong YS, Shin WJ, Lee AR, Kim YS, Kim KJ, Lee KS. Tracing the geographical origin of beefs being circulated in Korean markets based on stable isotopes. Rapid Commun Mass Spectrom. 2010;24(1):155–9. pmid:19960503
  26. 26. Bong YS, Gautam MK, Lee KS. Origin assessment of domestic and imported beef sold in the Korean markets using stable carbon and oxygen isotopes. Food Sci Biotechnol. 2012;21(1):233–7. https://doi.org/10.1007/s10068-012-0030-4
  27. 27. Kim KS, Kim JS, Hwang IM, Jeong IS, Khan N, Lee SI, et al. Application of stable isotope ratio analysis for origin authentication of pork. Korean J Food Sci Anim Resour. 2013;33(1):39–44. https://doi.org/10.5851/kosfa.2013.33.1.39
  28. 28. Park YM, Lee CM, Hong JH, Jamila N, Khan N, Jung JH, et al. Origin discrimination of defatted pork via trace elements profiling, stable isotope ratios analysis, and multivariate statistical techniques. Meat Sci. 2018;143:93–103. pmid:29715666
  29. 29. Shin WJ, Choi SH, Ryu JS, Song BY, Song JH, Park S, et al. Discrimination of the geographic origin of pork using multi‐isotopes and statistical analysis. Rapid Commun Mass Spectrom. 2018;32(21):1843–50. pmid:30074644
  30. 30. Chung IM, Park I, Yoon JY, Yang YS, Kim SH. Determination of organic milk authenticity using carbon and nitrogen natural isotopes. Food Chem. 2014;160:214–8. pmid:24799230
  31. 31. Chung IM, Kim JK, Lee KJ, Son NY, An MJ, Lee JH, et al. Discrimination of organic milk by stable isotope ratio, vitamin E, and fatty acid profiling combined with multivariate analysis: A case study of monthly and seasonal variation in Korea for 2016–2017. Food Chem. 2018;261:112–23. pmid:29739571
  32. 32. Chung IM, Kim JK, Yang YJ, An YJ, Kim SY, Kwon C, et al. A case study for geographical indication of organic milk in Korea using stable isotope ratios-based chemometric analysis. Food Control. 2020;107:106755. https://doi.org/10.1016/j.foodcont.2019.106755
  33. 33. Chung IM, Kim JK, Prabakaran M, Yang JH, Kim SH. Authenticity of rice (Oryza sativa L.) geographical origin based on analysis of C, N, O and S stable isotope ratios: a preliminary case report in Korea, China and Philippine. J Sci Food Agric. 2015;96(7):2433–9. pmid:26243037
  34. 34. Chung IM, Kim JK, Jin YI, Oh YT, Prabakaran M, Youn KJ, et al. Discriminative study of a potato (Solanum tuberosum L.) cultivation region by measuring the stable isotope ratios of bio-elements. Food Chem. 2016;212:48–57. pmid:27374505
  35. 35. Chung I M, Han JG, Kong WS, Kim JK, An MJ, Lee JH, et al. Regional discrimination of Agaricus bisporus mushroom using the natural stable isotope ratios. Food Chem. 2018;264:92–100. pmid:29853410
  36. 36. Park JH, Choi SH, Bong YS. Geographical origin authentication of onions using stable isotope ratio and compositions of C, H, O, N, and S. Food Control. 2019;101:121–5. https://doi.org/10.1016/j.foodcont.2019.02.023
  37. 37. Choi SH, Bong YS, Park JH, Lee KS. Geographical origin identification of garlic cultivated in Korea using isotopic and multi-elemental analyses. Food Control. 2020;111:107064. https://doi.org/10.1016/j.foodcont.2019.107064
  38. 38. Sharp Z. Principles of Stable Isotope Geochemistry, 2nd ed. Prentice Hall; 2017 https://doi.org/10.25844/h9q1-0p82
  39. 39. KREI. 2018 Food Balance Sheet. Korea Rural Economic Institute. Korean; 2019. pp. 317. Available from: https://www.krei.re.kr/krei/researchReportView.do?key=67&pageType=010101&biblioId=525942
  40. 40. Richards MP, Fuller BT, Hedges RE. Sulphur isotopic variation in ancient bone collagen from Europe: implications for human palaeodiet, residence mobility, and modern pollutant studies. Earth Planet Sci Lett. 2001;191(3–4):185–190. https://doi.org/10.1016/S0012-821X(01)00427-7
  41. 41. Choi S.H., Shin W.J., Bong Y.S., Lee K.S., 2021. Determination of the geographic origin of garlic using the bioelement content and isotope signatures. Food Control 130, 108339. https://doi.org/10.1016/j.foodcont.2021.108339
  42. 42. Endo T, Hayasaka M, Ogasawra H, Kimura O, Kotaki Y, Haraguchi K. Relationships among mercury concentration, and stable isotope ratios of carbon and nitrogen in the scalp hair of residents from seven countries: Effects of marine fish and C4 plants consumption. PLoS ONE. 2015;10(6):e0128149. pmid:26065892
  43. 43. Endo T, Hayasaka M, Hisamichi Y, Kimura O, Haraguchi K. Carbon and nitrogen stable isotope ratios and mercury concentration in the scalp hair of residents from Taiji a whaling town. Mar Pollut Bull. 2013;69(1–2):116–21. pmid:23453817
  44. 44. Lehn C, Rossmann A, Graw M. Provenancing of unidentified corpses by stable isotope techniques–presentation of case studies. Sci Justice. 2015;55(1):72–88. pmid:25577010
  45. 45. Thompson AH, Chesson LA, Podlesak DW, Bowen GJ, Cerling TE, Ehleringer JR. Stable isotope analysis of modern human hair collected from Asia (China, India, Mongolia, and Pakistan). Am J Phys Anthropol. 2010;141(3):440–51. pmid:19918992
  46. 46. Yoshinaga J, Minagawa M, Suzuki T, Ohtsuka R, Kawabe T, Inaoka T, et al. Stable carbon and nitrogen isotopic composition of diet and hair of Gidra‐speaking Papuans. Am J Phys Anthropol. 1996;100(1):23–34. https://doi.org/10.1002/(SICI)1096-8644(199605)100:1%3C23::AID-AJPA3%3E3.0.CO;2-E pmid:8859951
  47. 47. Hedges R, Rush E, Aalbersberg W. Correspondence between human diet, body composition and stable isotopic composition of hair and breath in Fijian villagers. Isotopes Environ Health Stud. 2009;45(1):1–17. pmid:19191122
  48. 48. Umezaki M, Naito YI, Tsutaya T, Baba J, Tadokoro K, Odani S, et al. Association between sex inequality in animal protein intake and economic development in the Papua New Guinea highlands: The carn and nitrogen isotopic composition of scalp hair and fingernail. Am J Phys Anthropol. 2016;159(1):164–73. pmid:26331634
  49. 49. Schoeller DA, Minagawa M, Slater R, Kaplan IR. Stable isotopes of carbon, nitrogen and hydrogen in the contemporary North American human food web. Ecol Food Nutr. 1986;18(3):159–170. https://doi.org/10.1080/03670244.1986.9990922
  50. 50. Reynard LM, Burt N, Kwon HEC, Tuross N. Limits and possibilities in the geolocation of humans using multiple isotope ratios (H, O, N, C) of hair from east coast cities of the USA. Isot Environ Health Stud. 2016;52(4–5):498–512. https://doi.org/10.1080/10256016.2016.1143821
  51. 51. Bol R, Pflieger C. Stable isotope (13C, 15N and 34S) analysis of the hair of modern humans and their domestic animals. Rapid Commun Mass Spectrom. 2002;16(23):2195–200. pmid:12442295
  52. 52. Lehn C, Kalbhenn EM, Rossmann A, Graw M. Revealing details of stays abroad by sequential stable isotope analyses along human hair strands. Int J Legal Med. 2018;133(3):935–47. pmid:29876636
  53. 53. FAO. FAO yearbook, Fishery and Aquaculture Statistics, Section 2 –Food balance sheets and fish contribution to protein supply–by country. Food and Agriculture Organization of the United Nations, New York; 2019. pp. 43–277.
  54. 54. Ritchie H, Roser M. Meat and dairy production; 2019. Available from: https://ourworldindata.org/meat-production#citationhttps://ourworldindata.org/meat-production#per-capita-meat-consumption
  55. 55. Dugee O, Khor GL, Lye MS, Luvsannyam L, Janchiv O, Jamyan B, et al. Association of major dietary patterns with obesity risk among Mongolian men and women. Asia Pac J Clin Nutr. 2009;18(3):433–40. pmid:19786392
  56. 56. Vogel JC. Fractionation of the carbon isotopes during photosynthesis. In: Fractionation of the Carbon Isotopes During Photosynthesis. Heidelberg: Springer; 1980. pp. 5–29. https://doi.org/10.1007/978-3-642-46428-7_1
  57. 57. Chesson LA, Meier-Augenstein W, Berg GE, Bataille CP, Bartelink EJ, Richards MP. Basic principles of stable isotope analysis in humanitarian forensic science. In: Forensic Science and Humanitarian Action: Interacting with the Dead and the Living. 1st ed. Parra RC, Zapico SC, Ubelaker DH, editors. John Wiley & Sons Ltd; 2020. pp. 285–310. https://doi.org/10.1002/9781119482062.ch20
  58. 58. Shin WJ, Kim I, Ryu JS, Lee KS. Dual isotopes of nitrate in Korean fertilizers and their application for identifying nitrate sources. Episodes. 2019;42(4):313–9. https://doi.org/10.18814/epiiugs/2019/019025
  59. 59. Focken U. Stable isotopes in animal ecology: the effect of ration size on the trophic shift of C and N isotopes between feed and carcass. Isot Environ Health Stud. 2001;37(3):199–211. pmid:11924851
  60. 60. Petzke KJ, Boeing H, Klaus S, Metges CC. Carbon and nitrogen stable isotopic composition of hair protein and amino acids can be used as biomarkers for animal-derived dietary protein intake in humans. Am Soc Nut Sci. 2005;135(6):1515–20. pmid:15930462
  61. 61. Richards MP, Pettitt PB, Stiner MC, Trinkaus E. Stable isotope evidence for increasing dietary breadth in the European mid-Upper Paleolithic. Proc Natl Acad Sci U S A. 2001;98(11):6528–32. pmid:11371652
  62. 62. Zanden MJV, Fetzer WW. Global patterns of aquatic food chain length. Oikos 2007;116(8):1378–88. https://doi.org/10.1111/j.0030-1299.2007.16036.x