Regional Variations of Insulin Secretion and Insulin Sensitivity in Japanese Participants With Normal Glucose Tolerance

Introduction

Diabetes mellitus is a metabolic disease characterized by hyperglycemia resulting from deficiency in insulin secretory capacity, insulin action or both (1). Kodama et al. compared the racial/ethnic differences for an optimal balance between insulin secretion and insulin sensitivity among Africans, Caucasians, and East Asians (2), and found that East Asians had better insulin sensitivity and lower insulin secretory capacity than Africans and Caucasians. However, a change in the living environment, such as through immigration, affects the balance between insulin secretion and insulin sensitivity even in individuals with the same racial/ethnic features. Fujimoto et al. reported that, after adjusting for the body mass index (BMI), the insulin response of Japanese-Americans was higher than that of native Japanese, suggesting that factors other than body size/adiposity play a distinct role in mediating insulin secretion (3).

Nutrients are plausible candidates as factors that mediate an altered insulin response (4). Immigration from Asia to Europe or the USA might change cultural dietary patterns, such as the intake of staple food, meat, and fish, with resultant changes in the balance among protein, fat, and carbohydrate intake (5). The elevation of plasma concentrations of free fatty acid (FFA) stimulated insulin secretion during the euglycemic or hyperglycemic clamp state (6). McGarry et al. confirmed that, in the insulin secretory pathway, a fatty acid-derived lipid moiety plays a pivotal role that is common to the action of a wide variety of secretagogues: highly saturated FFA (SFA), but not n-3 polyunsaturated fatty acids (PUFA), overstimulates the insulin response independent of the glucose levels (7). The Western diet is characterized by a high intake of meat and reduced intake of fish, which results in high SFA and low n-3 PUFA plasma concentrations (5). Taken together, the westernized dietary pattern could be closely related to an altered insulin response and insulin sensitivity through the hazardous effects of SFA, which is considered to be a form of “lipotoxicity” (8, 9). Regional differences in dietary patterns in Asian countries might affect the balance of insulin response and sensitivity. However, this notion is yet to be validated.

To clarify the regional differences in the insulin response and sensitivity and their relationship to nutrients, we compared the insulin secretory response during an oral glucose tolerance test in Japanese participants. We selected participants with normal glucose tolerance (NGT) to minimize the effects of personal glucose intolerance and chose four distinct areas of Japan with regard to the food environment: Fukushima, Nagano, Tokushima, and Okinawa (Supplementary Figure 1).

Materials and Methods

Study Design and Study Participants

This observational retrospective cohort study included an analysis in the Kokuho Database (KDB) which includes data form Japanese National Health Insurance and the data obtained during annual specific health checkups according to “The Specific Health Check and Guidance in Japan” (10, 11). Of the 47 regions (prefectures) in Japan, we selected four distinct regions based on the food environment (12) and structure of non-communicable diseases (13): Fukushima, Nagano, Tokushima, and Okinawa. From the overall list of 4,393 candidates who had undergone an oral 75-g glucose tolerance test (OGTT) between April 2005 and December 2016 in the KDB of four regions, we randomly selected study participants with the assumption that the participants were representative of candidates from each region (Supplementary Figure 2).

To obtain representative data for the daily intake of nutrients and food in four regions, we used data from the National Health and Nutrition Survey (NHNS), a series of nationally representative cross-sectional nutrition surveys conducted by local public health centers (12). The personalized data of the participants in the NHNS database were not available for their privacy policy, therefore, the researchers could only report aggregate data. However, since the nutritional samples in NHNS had been extracted and collected with a statistically rigorous manner to mimic real distributions in 47 regions (prefectures), we considered these as representative values for each region.

The study protocol was approved by the Fukushima Medical University Ethics Committee (approval number 30168). An opt-out notice for tacit consent was posted to publish the research information and ensure the opportunity for refusal. This study was conducted in accordance with the Ethical Guidelines for Medical and Health Research Involving Human Subjects enacted by the MHLW of Japan.

Anthropometric and Biochemical Measurements

The procedures for the anthropometric and biochemical measurements have been described previously (11). Briefly, trained staff measured height, body weight, waist circumference, and systolic and diastolic blood pressure (SBP and DBP) using a standard sphygmomanometer or an automated device on the right arm after resting for 5 min in a seated position. The participants reported their age, sex, smoking and alcohol intake, medications, and history of cardiovascular diseases and cerebrovascular disease in the questionnaire. Blood samples were collected after overnight fasting >12 h. Centrifuged samples were analyzed using an automatic clinical chemistry analyzer within 3 h of sample collection. All blood samples were analyzed at local laboratories by using the standard methods for laboratory tests recommended by the Japan Society of Clinical Chemistry. The analyzed factors were fasting plasma glucose (FPG), low density lipoprotein (LDL)-cholesterol, high density lipoprotein (HDL)-cholesterol, and triglyceride. Hypertension was defined as SBP ≥ 140 mmHg, DBP ≥ 90 mmHg, or use of antihypertensive medication. Dyslipidemia was defined as LDL-cholesterol 140 ≥ mg/dL, HDL- cholesterol ≤ 40 mg/dL, TG ≥ 150 mg/dL, or use of anti-dyslipidemic drugs. If necessary, the total cholesterol was calculated by the Friedewald's method (14), and participants whose TG was ≥400 mg/dL were excluded from the analysis.

Indices of Dietary Factors

Descriptive reports of the NHNS data, including health status, food and nutrient intakes, and lifestyles of Japanese civilians, are available (15). Household members were randomly selected from 475 out of ~1,00,000 nationwide census units in Japan. The average indices of dietary factors are available publicly in each region. Data on the dietary intake were collected using a 1-day semi-weighted household dietary record. Participants recorded food consumption at each household for a typical day. Trained fieldworkers, mainly registered dietitians, visited each household and checked the completeness of the data recording forms and, if necessary, confirmed portion sizes by using commercially available food models or food booklets, and corrected any missing and/or illogical information. In accordance with a survey manual of the NHNS, the trained fieldworkers converted these estimates of portion sizes or quantity of foods that were recorded into the weights of the foods and coded each food item according to the NHNS food number lists to calculate the energy and nutrient intake based on the record of household food consumption. The energy intake was calculated by multiplying the intake amounts of protein, lipids, and carbohydrates by applying the Atwater factors. Energy and nutrient intakes from food items with these values that are listed in the Standard Tables of Food Composition in Japan were calculated accordingly (16).

OGTT and Definition of Glucose Tolerance

After overnight fasting (>12 h), venous blood samples were collected before and after 30, 60, and 120 min of the ingestion of 75 g glucose (Trelan-G75; AY Pharmaceuticals Co., Ltd., Tokyo, Japan). Plasma samples for the assessment of insulin were collected in EDTA tubes on ice, separated after centrifugation (~3,000 Åg for 15 min at 4°C), and stored at −80°C until the assays were performed. Insulin was assayed routinely using a commercially available electrochemiluminescence immunoassay (CLEIA). The status of glucose tolerance was defined according to the WHO 2006 criteria (13): normal glucose tolerance (NGT; fasting plasma glucose <110 mg/dL and 2 h post-load plasma glucose <140 mg/dL), impaired glucose tolerance (IGT; fasting plasma glucose levels of 110–125 md/dL and/or 2 h post-load glucose levels of 140–199 mg/dL) and type 2 diabetes (fasting plasma glucose ≥ 126 mg/dL and/or 2 h post-load glucose ≥ 126 mg/dL).

Insulin Secretion and Insulin Sensitivity

The OGTT-based indices of insulin secretion and insulin sensitivity were extracted from the literature. Indices of insulin secretion included insulin₀, insulin_AUC, insulin_AUC/glucose_AUC, homeostatic model assessment estimates of β-cell function (HOMA-β), and the insulinogenic index. Insulin secretion included homeostatic model assessment estimates of insulin resistance (HOMA-IR), insulin sensitivity index (Matsuda index), quantitative insulin sensitivity check index (QUICKI), and the overall index included insulin₁₂₀ and disposition index (1, 17–19). Kahn et al. (1) first reported that the relationship between insulin sensitivity and β-cell function (acute insulin response, AIR_glucose) in human participants with NGT showed a hyperbolic curve in relation to the results of the intravenous glucose tolerance test (IVGTT) based on the minimum model of glucose kinetics developed by Bergman et al. (20). We used the indices based on OGTT to measure the hyperbolic relationship between the insulin sensitivity index (Matsuda index) (18) and β-cell function (insulinogenic index and insulin_AUC/glucose_AUC) (17, 21).

Statistical Analysis

Categorical variables are expressed as frequencies with percentages, and continuous variables are expressed as means with standard deviation (SD) or medians with interquartile ranges (IQRs). Intergroup comparisons of non-parametric variables were undertaken using the Kruskal–Wallis test, followed by Dunn's multiple comparisons test. Categorical variables are shown as percentages and were analyzed using Fisher's exact test. We used the values of participants from Tokushima as a control because the indices of insulin secretion and sensitivity were different from those of the other three regions. The sex-stratified Pearson product correlation coefficients were calculated for the correlation between indices of insulin response and sensitivity and that of nutrients for participants from the four regions. We performed a multiple regression analysis to estimate insulin_AUC/glucose_AUC. Although OGTT data sample size was 2,259, variables on the National nutritional survey were not available for each of 2,259 due to restrictions by a protection law on personal information. Instead, a total of 8 samples were obtained for the mean values of men and women in 4 prefectures. Given limited numbers in our samples (n = 8), overfitting multiple regression models could produce misleading coefficients, R-squared and p-values. We therefore adopted two independent variables in the model: one dietary composition and the other variable such as age, sex and BMI. Statistical analyses were performed using Prism version 8 or SPSS statistics version 27.0 for Windows (SPSS, IBM Corp, Armonk, NY, USA). All P-values were two-sided, and P < 0.1 was considered statistically significant.

Results

Of the 4,393 patients in the four regions who underwent an OGTT, we selected 2,259 (men, 1,190; women, 1,069) participants with NGT for inclusion in this study analysis (Supplementary Figure 2). The general characteristics of the participants are shown in Supplementary Figure 3. The median age, BMI, and waist circumference of men and women, respectively, were 56.0 and 58.0 (47.0–63.0 and 50.0–63.0) years, 25.5 and 25.9 (23.8–27.3 and 24.0–27.8) kg/m², and 89.2 and 90.0 (86.0–94.0 and 84.9–94.0) cm. Among the men, the age, BMI, waist circumference, SBP, DBP, and triglycerides, except LDL-C and HDL-C were significantly different among the four regions. For women, all items evaluated showed significant differences among the regions.

The glucose and insulin levels during the OGTT are shown in Supplementary Figure 4 and Figure 1A (men) and Figure 1B (women). The glucose levels were comparable among the four regions in the groups of both men and women. For the insulin responses in men, Fukushima ≈ Nagano ≈ Okinawa > Tokushima in all, 25 ≤ BMI <30 and BMI ≥ 30 and Nagano ≈ Okinawa > Fukushima ≈ Tokushima in BMI <25 (Supplementary Figure 4 and Figure 1A). For the insulin response in women, Nagano ≈ Okinawa > Fukushima ≈ Tokushima in all and BMI <25 and 25 ≤ BMI <30 (Supplementary Figure 4 and Figure 1B).

FIGURE 1

Figure 1. Changes in glucose and insulin levels in men (A) and women (B) during the 75-g oral glucose tolerance test. Glucose levels are shown in left panels and insulin levels are shown in right panels for all or stratified categories across three body mass index (BMI), BMI <25, 25 ≤ BMI <30, and BMI ≥ 30. The participants were randomly selected from the National Health Insurance database in Fukushima (°), Nagano (Δ), Tokushima (□) and Okinawa (•), Japan. Data are presented as median (25–75th percentile). To compare deference between the regions at each point, P-values were obtained by Kruskal–Wallis test, followed by Dunn's multiple comparisons test using Tokushima as a control. *P < 0.1, **P < 0.05, and ***P < 0.01 vs. Tokushima.

The indices of insulin secretory capacity and insulin sensitivity are shown in Table 1A (men) and Table 1B (women). Indices of the insulin secretory capacity, such as insulin₀, insulin₁₂₀, insulin_AUC, HOMA-β, and the insulinogenic index were larger overall in men and women in Fukushima, Nagano, and Okinawa than in those from Tokushima. The insulin_AUC/glucose_AUC ratio during the OGTT in men and women in the three BMI categories are shown in Supplementary Figure 5. Insulin sensitivity indices, such as HOMA-IR, Matsuda index, QUICKI, and disposition index, were better in participants from Tokushima than from those in the other regions (Tables 1A, 1B).

TABLE 1A

Table 1A. Indices of insulin secretion and insulin sensitivity in men.

TABLE 1B

Table 1B. Indices of insulin secretion and insulin sensitivity in women.

The region means of anthropometry, biochemical parameters and the daily intake of nutrients and food are shown in Table 2. Ages were different in men (Tokushima > Fukushima > Nagano ≈ Okinawa), but were comparable in women. Among anthropometry and biochemical parameters, BMI, waist circumference, systolic and diastolic blood pressure and triglycerides were higher in Okinawa than in Tokushima and Fukushima and systolic blood pressure (Fukushima > Tokushima) and dietary fiber (Fukushima and Nagano > Tokushima) were different. Characteristics for nutritional indices in four regions were summarized as follows: Fukushima (high intake: grain, fish, carbohydrate calorie, and dietary fiber; low intake: meat and fat calories), Nagano (high intake: grain, fish, and dietary fiber; low intake: meat), Tokushima (at the midpoint in food factors), and Okinawa (high intake: fat calories, low intake: grain, fish, carbohydrate and carbohydrate calories, and dietary fiber).

TABLE 2

Table 2. The mean values of anthropometry, biochemical parameters and the daily intake of nutrients and food in four regions obtained from the 2016 National Health and Nutrition Survey in Japan.

The correlation between the indices of nutrients and food intake (Table 3) and insulin secretion and sensitivity in four regions are shown in Supplementary Figure 6 and Figure 2. Overall, insulin₁₂₀ and insulin_AUC/glucose_AUC correlated negatively with fish intake, carbohydrate intake, and carbohydrate calories and positively correlated with fat calories. However, HOMA-β was negatively correlated with fish intake and carbohydrate calories and positively correlated with fat calories.

TABLE 3

Table 3. Simple and multiple regression analysis to estimate insulin_AUC/glucose_AUC.

FIGURE 2

Figure 2. The correlation between the indices of dietary factors and indices of insulin secretion and sensitivity in Fukushima men (°) and women (°), Nagano men (Δ) and women (Δ), Tokushima men (□) and women (□), and Okinawa men (•) and women (•). r: The Pearson product correlation coefficient and p: probability. Solid and dotted lines indicate r with 95% confidential intervals.

As shown in Figure 3, the plot between the insulin sensitivity index (Matsuda index) and insulin_AUC/glucose_AUC (left panel) or the insulinogenic index (right panel) for all participants showed hyperbolic relationships. In men or women, as shown in the figure insets, the points shifted from the lower right toward the upper left in the following order: Tokushima < Fukushima < Nagano < Okinawa.

FIGURE 3

Figure 3. The plot between the insulin sensitivity index (ISI; Matsuda index) and insulin_AUC/glucose_AUC (left panel) or the insulinogenic index (right panel) in all participants including those from four regions. The stabilization points of insulin sensitivity and the insulin response are depicted on plots in all participants and are shown in the inset as Fukushima men (°) and women (°), Nagano men (Δ) and women (Δ), Tokushima men (□) and women (□), and Okinawa men (•) and women (•).

Finally, we performed a multiple regression analysis to estimate insulin_AUC/glucose_AUC (Table 3). After corrected for age (Model 2) or sex (Model 3), energy intake, fish intake, protein intake, fat calorie and carbohydrate intake and calorie were significantly associated with insulin_AUC/glucose_AUC. However, after corrected for BMI (Model 4), only fat calorie showed a borderline significance.

Discussion

We compared the indices of insulin response and sensitivity of men and women with NGT in four different areas of Japan, and noted two major findings. First, there were regional variations in insulin secretion and sensitivity. This difference was observed even within the same BMI category. Second, the indices of insulin secretion correlated positively with fat intake and negatively with the intake of fish, carbohydrate calories, and dietary fiber. Although racial/ethnic differences in the optimal balance for insulin secretion and insulin sensitivity are well-known, this study is the first to report that the optimal balance differs significantly among participants from different regions in Japan and that nutritional factors could be linked to these differences, independently of body size/adiposity. Furthermore, this study suggests that the insulin response and insulin sensitivity can vary among adult individuals, even within the same race and the same country, and is therefore likely affected by environmental/lifestyle factors independently of genetic traits.

Regional Variations in Indices of Insulin Secretion and Insulin Sensitivity

The insulin secretory capacity and insulin sensitivity differ between races. Kodama et al. discovered that the stabilization points of insulin sensitivity and insulin response in a hyperbolic relationship model differ among African, Caucasian, and East Asian NGT cohorts (2). For healthy NGT subjects, differences in body size/compositions are well-known among Africans (less visceral fat and more skeletal muscle mass than Caucasians), Caucasians (more visceral fat and less skeletal muscle mass than Africans), and East Asians (less visceral fat and smaller height, weight, and BMI than the other two ethnic groups). Therefore, it is possible that the differences in insulin secretion and sensitivity may be affected by differences in body size/composition. Another explanation for this difference is the innate capacity for insulin secretion, namely genetic susceptibility (22, 23). A recent meta-analysis of genome-wide association studies in 77,418 East Asians with type 2 diabetes mellitus revealed that East Asians share common variants for type 2 diabetes mellitus with Europeans, but also exhibit unique undescribed associations of those variants (24). In this study, we selected only Japanese participants, and examined the relationship between dietary factors and indices of insulin secretion and sensitivity in participants from distinct regions. We found significant regional differences in insulin response and sensitivity, even within participants of the same race and country: the order of the insulin responses was Okinawa > Nagano > Fukushima > Tokushima (Figure 1A).

As obesity/adiposity might have an impact on insulin secretion and sensitivity, we subdivided the participants according to their BMI. In men, the insulin responses were Nagano ≈ Okinawa > Fukushima ≈ Tokushima in BMI <25 and Fukushima ≈ Nagano ≈ Okinawa > Tokushima in 25 ≤ BMI <30 and BMI ≥ 30 (Figure 1A). In women, the insulin responses were Nagano ≈ Okinawa > Fukushima ≈ Tokushima in BMI <25 and 25 ≤ BMI <30 (Figure 1B). Thus, regional differences in insulin response were assumed partially due to body size/adiposity, but also due to other unrecognized factors. The larger variability [median (25–75th percentile)] of insulin levels in men from Fukushima, Nagano and Okinawa may suggest an intra-regional heterogeneity. We assessed representative indices for insulin secretion and insulin sensitivity proposed (17–19, 21) and found interregional difference in some indices (Insulin_AUC, Insulin_AUC/Glucose_AUC, HOMA-β, insulinogenic index, HOMA-IR, Matsuda index, and QUICKI), but not in some (Disposition index). Although we cannot explain the discrepancy from our results, several possibilities need to be considered such as difference in diagnostic quality and/or mechanistic background which might be affectete by different foods or nutrients, i.e., fish, fat, carbohydrate in these indices.

Racial/ethnic differences in the optimal balance for insulin secretion and insulin sensitivity are dependent on body size/adiposity and/or genetic background (2). We found that the specified points of the insulin sensitivity index (Matsuda index) and insulin secretion (insulin_AUC/glucose_AUC) were clearly different for the four regions, and these variations cannot be explained only by the degree of obesity and visceral fat given the variations in the differences in the distribution of insulin response among the four regions even after correction for BMI.

Potential Mechanisms: Dietary Factors and Insulin Secretion

We evaluated the correlation between dietary factors and indices of insulin secretory capacity and sensitivity in men and women from the four Japanese regions. Insulin₁₂₀, insulin_AUC/glucose_AUC, and HOMA-β were negatively correlated with the intake of fish, carbohydrate calories, and dietary fiber and were positively correlated with fat calories (Figure 2 and Supplementary Figure 6). Three possible mechanisms might explain the regional difference in insulin response as differences in the intake of: fatty acids; carbohydrate, and fish.

Difference in SFA Intake

In this present study, the indices of insulin secretion did not correlate with the daily energy intake, but positively correlated with fat-derived energy and negatively with carbohydrate calories and dietary fiber. In the multiple regression analysis to estimate insulin_AUC/glucose_AUC, only fat calorie showed a borderline significance after corrected for BMI (Table 3). Marshall et al. reported that fasting insulin levels were positively correlated with fat intake and negatively associated with carbohydrate intake in Hispanic and Caucasian cohorts from Colorado, USA (25). A report in the Netherlands reported that fasting C-reactive peptide was associated positively with the total fat intake and negatively with carbohydrates (17). Maron et al. reported that total fat and SFA intakes positively and negatively correlated with fasting insulin secretion and carbohydrate intake, respectively, in non-diabetic individuals with coronary artery disease (CAD) (26). In contrast to the abovementioned studies (17, 25, 26) that evaluated such relationships in fixed regions, our study compared the average values in multiple regions, suggesting that the geographical differences in the insulin response may be linked to the degree of the intake of fat, carbohydrate, and dietary fiber.

Insulin secretion was augmented in Japanese immigrants in Seattle, USA (3, 22) and Hawaii, USA (27) as compared to that in native Japanese. In Pima Indians in the USA, the transition from an agrarian to a modern society was associated with the consumption of increasing amounts of dietary fat, decreasing amounts of dietary carbohydrate, and a decline in insulin sensitivity (28). In Japan, the era of westernization is fairly different among the four regions tested. In the earlier era (c.1945–1972) after World War II, the American food culture under the presence of US military bases in Okinawa caused a rapid shift in food culture from the Japanese (29) to a westernized dietary style (30, 31). In contrast, the other three regions more gradually shifted from the Japanese to a westernized dietary style for topographic, historical, and cultural reasons (12, 29, 32).

This study did not conclusively ascertain the reason why fat calories are related to insulin secretion, but one can presume two potential mechanisms. First, the SFA in the high-fat diet may cause insulin resistance. Second, the SFA may directly alter insulin secretion. Previously, we reported that SFA stimulates insulin secretion in normal animal models independently of glucose concentration (8, 9). McGarry et al. have shown that FFA directly stimulates insulin secretion in cultured pancreatic β-cells by glucose-insulin coupling, and that the response is affected by chain length and the FFA degree of saturation (7). Incretins are hormones secreted from the digestive tract on nutrient intake that act on pancreatic β-cells (33). Gastric inhibitory polypeptide (GIP), but not glucagon-like peptide-1 (GLP-1), can be stimulated by dietary fat in individuals with NGT (34). Currently, it is unknown whether postprandial incretin secretion is linked to the relationship between fat intake and insulin response. Taken above all, insulin sensitivity and insulin secretion are mutually and closely interlinked, and FFA effects might modulate the stabilization points of insulin sensitivity and insulin response (1, 21).

Difference in Carbohydrate Intake

The fasting insulin level correlated negatively with carbohydrate intake in healthy Hispanic and White cohorts (25), normal healthy Netherlanders (17), and Americans with coronary artery disease (26). Insulin levels positively correlated with sucrose intake and negatively correlated with complex carbohydrates and dietary fiber (25). Our results agree with previous studies (17, 25, 26). Low carbohydrate intake supposedly lowers insulin levels, increases circulating FFA levels, and enhances fatty acid oxidation (25, 26). There is a possibility that insulin sensitivity and insulin response may differ by the type of carbohydrates. Diets with high dietary fiber intake show a low glycemic index, wherein even the same amount of calories can reduce postprandial hyperinsulinemia, easily reduce body fat, and improve insulin sensitivity (35). The differences in carbohydrate intake in the four regions tested may reflect the differences in habits of consuming grains and fats (12, 29–32). Because there was no difference in protein intake among the four regions, there is a possibility that a difference in low carbohydrate calories and high fat calories may be linked to differences in insulin secretion. Low carbohydrate calories can reduce the circulating insulin level, which promotes higher circulating FFA levels for use in oxidation and production of ketone bodies (36). Therefore, it would be interesting to examine whether the difference in the glycemic index is linked to regional differences in insulin secretion.

Difference in Fish Intake

We found that fish intake was associated with decreased insulin secretion. The intake of meat (beef and chicken), when matched for macronutrient and energy content, is more insulinogenic than that of fish (37). Insulin secretion may be reduced because of improved sensitivity from a higher fish intake (38). Conversely, fish n-3 PUFA may play a role in transcriptional regulation for adipocytokines and protection of the insulin secretory capacity through membrane receptors (39). The four regions of Japan differed greatly in the intake ratio of meat and fish, which may cause changes in fatty acid composition (40), thereby affecting insulin secretion and sensitivity (37–39).

Other Potential Mechanisms: Discrepancies Between the Different Indexes and Interindividual Variability of the Indexes for Insulin Secretion and Sensitivity

As discussed above, we found interregional difference in some indices, but not in some, used to measure insulin secretion and sensitivity. We have no biological explanation for the discrepancies between the different indexes. The notion may be also considered that these discrepancies, together with the large interindividual variability in the results and in the cardiometabolic phenotypic characteristics of the participants could result in false interregional differences. True mechanisms for “interregional difference” need to be clarified in future studies which will investigate possible determinants for insulin secretion and sensitivity including more precise individual characteristics.

Strengths and Limitations

The strengths of this study are the use of nationally representative data on dietary intake in Japan that were obtained by a standardized and unified protocol for population-level surveys. Moreover, as we used the large database of national health insurance in Japan, we could compare the relationship between dietary factors and indices of insulin response and sensitivity in the general population. However, this study has several limitations. First, we compared the dietary factors and OGTT-based indices of insulin secretion and sensitivity from different datasets and combined the data into representative variables for the four regions. Therefore, the results obtained from the two datasets should be interpreted cautiously. The correlation between dietary factors and variables observed in this study need to be confirmed by future studies using the same dataset for the relationship between dietary intake and insulin (secretion/sensitivity) indices. Second, we selected four regions for comparison, but the selection may be biased and, perhaps, non-representative of the nationwide data. Third, the utility of the household-based dietary record in estimating energy and nutrient intake is limited and the household response rates for the 2016 NHNS were low (44.4% in 2016) (12), which possibly could cause bias. Fourth, measurements of OGTT were optional at routine annual specific health checkups, and this is another possible source of bias. Fifth, the intra-regional variability for potential confounding factors such as age, BMI and other metabolic alterations may be relevant for the differences in the insulin response and insulin sensitivity. However, we did not perform multivariate analysis including these potential confounders because limited numbers for nutritional indices (n = 8) may cause overfitting multiple regression models such as misleading coefficients. Therefore, we limited explanatory parameters to two in the models (Table 3).

Conclusion

This is the first study to evaluate regional differences in insulin response and insulin sensitivity in Japanese participants and to show significant regional differences in these indices and that nutritional factors may be linked to these differences independently of body size/adiposity. Furthermore, insulin response and insulin sensitivity can vary among adult individuals, even within the same race and the same country, and are therefore likely affected by environmental/lifestyle factors independently of genetic traits.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by the Fukushima Medical University Ethics Committee (approval number 30168). An opt-out notice for tacit consent was posted to publish the research information and ensure the opportunity for refusal. The ethics committee waived the requirement of written informed consent for participation.

Author Contributions

KW analyzed data and wrote first manuscript. MSh designed the study, analyzed data, and wrote the manuscript with input from all authors. MH, YH, AK, RA, BW, DW, MSa, and HM contributed to the discussion and edited the manuscript. MSh is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by Japan Society for the Promotion of Science (JPSP) (Grant Number JP16K01823 and JP17K00924 to AK and MSh).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors sincerely acknowledge Ms. Katsuko Kumagai and the colleagues in the study group for public health nurse activity for longstanding partnership and excellent assistance to this study.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2021.632422/full#supplementary-material

References

1. Kahn SE, Prigeon RL, McCulloch DK, Boyko EJ, Bergman RN, Schwartz MW, et al. Quantification of the relationship between insulin sensitivity and β-cell function in human subjects. Evidence for a hyperbolic function. Diabetes. (1993) 42:1663–72. doi: 10.2337/diabetes.42.11.1663

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Kodama K, Tojjar D, Yamada S, Toda K, Patel CJ, Butte AJ. Ethnic differences in the relationship between insulin sensitivity and insulin response: a systematic review and meta-analysis. Diabetes Care. (2013) 36:1789–96. doi: 10.2337/dc12-1235

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Fujimoto WY, Akanuma Y, Kanazawa Y, Mashiko S, Leonetti D, Wahl P. Plasma insulin levels in Japanese and Japanese-American men with type 2 diabetes may be related to the occurrence of cardiovascular disease. Diabetes Res Clin Pract. (1989) 6:121–7. doi: 10.1016/0168-8227(89)90116-2

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Prentki M, Matschinsky Franz M, Madiraju SRM. Metabolic signaling in fuel-induced insulin secretion. Cell Metabolism. (2013) 18:162–85. doi: 10.1016/j.cmet.2013.05.018

PubMed Abstract | CrossRef Full Text | Google Scholar

5. LeCroy MN, Stevens J. Dietary intake and habits of South Asian immigrants living in Western countries. Nutr Rev. (2017) 75:391–404. doi: 10.1093/nutrit/nux023

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Stein DT, Esser V, Stevenson BE, Lane KE, Whiteside JH, Daniels MB, et al. Essentiality of circulating fatty acids for glucose-stimulated insulin secretion in the fasted rat. J Clin Invest. (1996) 97:2728–35. doi: 10.1172/JCI118727

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Dobbins RL, Chester MW, Stevenson BE, Daniels MB, Stein DT, McGarry JD. A fatty acid- dependent step is critically important for both glucose- and non-glucose-stimulated insulin secretion. J Clin Invest. (1998) 101:2370–6. doi: 10.1172/JCI1813

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Unger RH. Lipotoxicity in the pathogenesis of obesity-dependent NIDDM. Genetic and clinical implications. Diabetes. (1995) 44:863–70. doi: 10.2337/diab.44.8.863

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Shimabukuro M, Zhou YT, Levi M, Unger RH. Fatty acid-induced beta cell apoptosis: a link between obesity and diabetes. Proc Natl Acad Sci USA. (1998) 95:2498–502. doi: 10.1073/pnas.95.5.2498

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Nagai K, Iseki C, Iseki K, Kondo M, Asahi K, Saito C, et al. Higher medical costs for CKD patients with a rapid decline in eGFR: a cohort study from the Japanese general population. PLoS One. (2019) 14:e0216432. doi: 10.1371/journal.pone.0216432

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Hirai H, Asahi K, Yamaguchi S, Mori H, Satoh H, Iseki K, et al. New risk prediction model of coronary heart disease in participants with and without diabetes: assessments of the Framingham risk and Suita scores in 3-year longitudinal database in a Japanese population. Sci Rep. (2019) 9:2813. doi: 10.1038/s41598-019-39049-w

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Saito A, Imai S, Htun NC, Okada E, Yoshita K, Yoshiike N, et al. The trends in total energy, macronutrients and sodium intake among Japanese: findings from the 1995-2016 National Health and Nutrition Survey. Br J Nutr. (2018) 120:424–34. doi: 10.1017/S0007114518001162

PubMed Abstract | CrossRef Full Text | Google Scholar

13. World Health Organization. Definition and Diagnosis of Diabetes Mellitus and Intermediate Hyperglycaemia: Report of a WHO/IDF Consultation. Geneva: WHO (2006).

Google Scholar

14. Friedewald WT, Levy RI, Fredrickson DS. Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clin Chem. (1972) 18:499–502. doi: 10.1093/clinchem/18.6.499

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Ikeda N, Takimoto H, Imai S, Miyachi M, Nishi N. Data resource profile: The Japan National Health and Nutrition Survey (NHNS). Int J Epidemiol. (2015) 44:1842–9. doi: 10.1093/ije/dyv152

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Watanabe T. Food composition tables of Japan and the nutrient table/database. J Nutr Sci Vitaminol. (2015) 61:S25–7. doi: 10.3177/jnsv.61.S25

PubMed Abstract | CrossRef Full Text | Google Scholar

17. den Biggelaar LJ, Sep SJ, Eussen SJ, Mari A, Ferrannini E, van Greevenbroek MM, et al. Discriminatory ability of simple OGTT-based beta cell function indices for prediction of prediabetes and type 2 diabetes: the CODAM study. Diabetologia. (2017) 60:432–41. doi: 10.1007/s00125-016-4165-3

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Matsuda M, DeFronzo RA. Insulin sensitivity indices obtained from oral glucose tolerance testing: comparison with the euglycemic insulin clamp. Diabetes Care. (1999) 22:1462–70. doi: 10.2337/diacare.22.9.1462

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Bergman RN, Ader M, Huecking K, Van Citters G. Accurate assessment of β-cell function. Diabetes. (2002) 51:S212. doi: 10.2337/diabetes.51.2007.S212

CrossRef Full Text | Google Scholar

20. Bergman RN, Ider YZ, Bowden CR, Cobelli C. Quantitative estimation of insulin sensitivity. Am J Physiol. (1979) 236:E667–77. doi: 10.1152/ajpendo.1979.236.6.E667

CrossRef Full Text | Google Scholar

21. Retnakaran R, Shen S, Hanley AJ, Vuksan V, Hamilton JK, Zinman B. Hyperbolic relationship between insulin secretion and sensitivity on oral glucose tolerance test. Obesity. (2008) 16:1901–7. doi: 10.1038/oby.2008.307

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Fujimoto WY, Boyko EJ, Hayashi T, Kahn SE, Leonetti DL, McNeely MJ, et al. Risk factors for type 2 diabetes: lessons learned from Japanese Americans in Seattle. J Diabetes Investig. (2012) 3:212–24. doi: 10.1111/j.2040-1124.2012.00195.x

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Hulman A, Simmons RK, Brunner EJ, Witte DR, Færch K, Vistisen D, et al. Trajectories of glycaemia, insulin sensitivity and insulin secretion in South Asian and white individuals before diagnosis of type 2 diabetes: a longitudinal analysis from the Whitehall II cohort study. Diabetologia. (2017) 60:1252–60. doi: 10.1007/s00125-017-4275-6

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Spracklen CN, Horikoshi M, Kim YJ, Lin K, Bragg F, Moon S, et al. Identification of type 2 diabetes loci in 433,540 East Asian individuals. Nature. (2020) 582:240–5. doi: 10.1038/s41586-020-2263-3

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Marshall JA, Bessesen DH, Hamman RF. High saturated fat and low starch and fibre are associated with hyperinsulinaemia in a non-diabetic population: the San Luis Valley Diabetes Study. Diabetologia. (1997) 40:430–8. doi: 10.1007/s001250050697

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Maron DJ, Fair JM, Haskell WL. Saturated fat intake and insulin resistance in men with coronary artery disease. The Stanford Coronary Risk Intervention Project Investigators and Staff. Circulation. (1991) 84:2020–7. doi: 10.1161/01.CIR.84.5.2020

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Hara H, Egusa G, Yamakido M, Kawate R. The high prevalence of diabetes mellitus and hyperinsulinemia among the Japanese-Americans living in Hawaii and Los Angeles. Diabetes Res Clin Pract. (1994) 24:S37–42. doi: 10.1016/0168-8227(94)90225-9

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Boyce VL, Swinburn BA. The traditional Pima Indian diet. Composition and adaptation for use in a dietary intervention study. Diabetes Care. (1993) 16:369–71. doi: 10.2337/diacare.16.1.369

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Murakami K, Shinozaki N, Fujiwara A, Yuan X, Hashimoto A, Fujihashi H, et al. A systematic review of principal component analysis-derived dietary patterns in Japanese adults: are major dietary patterns reproducible within a country? Adv Nutr. (2019) 10:237–49. doi: 10.1093/advances/nmy079

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Tanaka H, Shimabukuro T, Shimabukuro M. High prevalence of metabolic syndrome among men in Okinawa. J Atheroscler Thromb. (2005) 12:284–8. doi: 10.5551/jat.12.284

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Willcox DC, Willcox BJ, Todoriki H, Suzuki M. The Okinawan diet: health implications of a low-calorie, nutrient-dense, antioxidant-rich dietary pattern low in glycemic load. J Am Coll Nutr. (2009) 28:500s−16s. doi: 10.1080/07315724.2009.10718117

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Katanoda K, Matsumura Y. National Nutrition Survey in Japan–its methodological transition and current findings. J Nutr Sci Vitaminol. (2002) 48:423–32. doi: 10.3177/jnsv.48.423

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Seino Y, Kuwata H, Yabe D. Incretin-based drugs for type 2 diabetes: focus on East Asian perspectives. J Diabetes Investig. (2016) 7 (Suppl. 1):102–9. doi: 10.1111/jdi.12490

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Yamane S, Harada N, Hamasaki A, Muraoka A, Joo E, Suzuki K, et al. Effects of glucose and meal ingestion on incretin secretion in Japanese subjects with normal glucose tolerance. J Diabetes Investig. (2012) 3:80–5. doi: 10.1111/j.2040-1124.2011.00143.x

CrossRef Full Text | Google Scholar

35. Ludwig DS, Ebbeling CB. The carbohydrate-insulin model of obesity: beyond “calories in, calories out”. JAMA Intern Med. (2018) 178:1098–103. doi: 10.1001/jamainternmed.2018.2933

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Adam-Perrot A, Clifton P, Brouns F. Low-carbohydrate diets: nutritional and physiological aspects. Obes Rev. (2006) 7:49–58. doi: 10.1111/j.1467-789X.2006.00222.x

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Comerford KB, Pasin G. Emerging evidence for the importance of dietary protein source on glucoregulatory markers and type 2 diabetes: different effects of dairy, meat, fish, egg, and plant protein foods. Nutrients. (2016) 8:446. doi: 10.3390/nu8080446

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Ouellet V, Marois J, Weisnagel SJ, Jacques H. Dietary cod protein improves insulin sensitivity in insulin-resistant men and women. Diabetes Care. (2007) 30:2816. doi: 10.2337/dc07-0273

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Wang X, Chan CB. n-3 polyunsaturated fatty acids and insulin secretion. J Endocrinol. (2015) 224:R97–106. doi: 10.1530/JOE-14-0581

CrossRef Full Text | Google Scholar

40. Okami Y, Ueshima H, Nakamura Y, Kondo K, Kadota A, Okuda N, et al. Time-related changes in relationships between the keys score, dietary lipids, and serum total cholesterol in Japan - NIPPON DATA80/90/2010. Circ J. (2018) 83:147–55. doi: 10.1253/circj.CJ-18-0628

PubMed Abstract | CrossRef Full Text | Google Scholar