Socio-economic inequalities in health have been widely studied in industrialized countries. Many health indicators of morbidity and mortality are positively associated with low socio-economic status (SES)( Reference Marmot and Wilkinson 1 ). Furthermore, lower socio-economic groups have a higher incidence of non-communicable diseases than those in more socio-economically advantaged groups( Reference James, Nelson and Ralph 2 , Reference Smith and Brunner 3 ). Dietary factors may contribute to the impact of morbidity and mortality related to chronic diseases. Adults of high SES tend to consume a healthy diet, with greater consumption of fruit, vegetables and whole-grain foods, whereas those of low SES are more likely to consume fatty meats and refined cereals( Reference Darmon and Drewnowski 4 – Reference Hassen, Castetbon and Cardon 9 ).
Many studies have shown that children and adolescents from low-SES backgrounds in Western countries are more likely to be obese( Reference El-Sayed, Scarborough and Galea 10 , Reference Barrisuo, Miqueleiz and Albaladejo 11 ) and low SES in childhood is associated with obesity in adult women( Reference Senese, Almeida and Fath 12 ). It is suggested that children from low-SES households in Western countries have less healthy diets, such as lower consumption of vegetables, fruits, fish, fibre-rich bread and low-fat milk, and higher consumption of meat products, soft drinks and energy-dense foods( Reference Eloranta, Lindi and Schwab 13 – Reference Northstone, Smith and Cribb 17 ). In addition to cross-sectional studies, a prospective study found causal associations between lower SES and unhealthy dietary profiles among European children( Reference Fernandez-Alvia, Bornhorst and Bammann 18 ). These studies used FFQ to evaluate dietary quality between different SES. Previous studies using dietary recall/records suggest that children in households with lower incomes consume less Ca, vitamin D and dairy foods in Canada( Reference Mark, Lambert and O’Loughlin 19 ), and a Korean study showed that they obtain less energy from protein and higher energy from carbohydrates( Reference Kwon, Oh and Park 20 ). A systematic review revealed a positive association between SES and micronutrient intake in children; however, only a few studies have been published( Reference Navakovic, Cavelaars and Geelen 21 ).
Very few studies on SES and health/diet have been reported in Japan, because Japanese society is considered an egalitarian society( Reference Marmot and Smith 22 ). However, income inequality in Japan has risen since the 1980s( 23 ) and studies on SES and health, mostly in adults, have been conducted since the 1990s( Reference Kagamimori, Gaina and Nasermoaddeli 24 ). Some studies show an association between household SES and diet in adults( Reference Fukuda and Hiyoshi 25 , Reference Nishi, Horikawa and Murayama 26 ). The child poverty rate has increased in Japan( 23 ); thus, an act to accelerate measures for disadvantaged children was enacted in 2013. However, there is little evidence of an association between household income and the diets of children in Japan. Only one study of pre-school children has been published( Reference Okubo, Miyake and Sasaki 27 ) and the results suggested that household income is not associated with child diet after adjustment for maternal education. However, the study was limited to a homogeneous population of pregnant women in one city and the frequency of the children’s food consumption was reported by the mothers. Further studies are needed to examine the associations between household income and the dietary intake of children.
School lunches have been part of an important policy to improve child nutrition in Japan since World War II( Reference Tanaka and Miyoshi 28 , Reference Ishida 29 ) and school lunch coverage was 98·4 % at primary schools in fiscal year 2014( 30 ). Previous studies in the USA suggested that participation in the national school lunch programme could modify the relationship between food insecurity and weight in school-aged children( Reference Huang and Bamidge 31 – Reference Nguen, Ford and Yaroch 34 ). Thus, we hypothesized that the school lunch programme would reduce potential disparities in diet due to household income among schoolchildren in Japan. Moreover, we thought that there may be some differences in the association between household income and diet in children on the days with school lunches compared with the days without. If the effects of household income on diet were less on days with school lunches than on days without, this would suggest that school lunches have a beneficial role. This hypothesis has not been examined in previous studies.
Therefore, the objectives of the present study were to examine the association between household income and food and nutrient intake among Japanese primary-school children and to determine potential differences between days with and without school lunches.
Methods
Study participants
The participants of the present study were all 5th grade children (10–11 years old) attending nineteen primary schools in four areas (prefectures) of East Japan. There were 1498 children registered and, after excluding children who had been absent for long periods, 1447 pairs of children and guardians were included.
Among the four areas, areas A (three schools in two cities) and C (eight schools in one city and one village) were rural residential areas, and areas B (six schools in one city) and D (two schools in one city) were urban residential areas. The schools were selected from public schools after we obtained permission from the respective Boards of Education of the cities and head teachers of the primary schools. Children in 5th grade were selected because they were able to complete a questionnaire, keep dietary records and take photographs of the foods they ate. In addition, 5th grade children are less busy at school than 6th grade children. Among the 1447 pairs, 1231 (85·1 %) child–guardian pairs agreed to participate in the study. We excluded 307 pairs who did not provide information about income, eighty-two pairs who did not complete four days of dietary records, five pairs who did not provide information about body height and weight, and one pair who did not answer the frequency of exercise questions; thus 836 pairs (57·8 %) were included in the present study.
Measures
Data were collected from September to December 2013.
Household socio-economic status and demographic characteristics
The parents/guardians of the children completed questionnaires on household SES and demographic characteristics. The questionnaires were distributed to the children at school and taken home for the parents to complete. Children returned the questionnaires to schools where they were collected for the study. The SES of participants was assessed as household income and level of education of the fathers and mothers of the children.
Participants were asked about their total annual household income, including salary, social benefits, family allowance and rental income of all household members during fiscal year 2012. Participants selected one of the following choices: (i) ‘less than 1 million yen’; (ii) ‘1 to less than 2 million yen’; (iii) ‘2 to less than 3 million yen’; (iv) ‘3 to less than 4 million yen’; (v) ‘4 to less than 6 million yen’; (vi) ‘6 to less than 8 million yen’; (vii) ‘8 million yen and over’; or (viii) ‘do not wish to answer’. The choice to limit responses to these options rather than use open questions was made to increase the response rate. We used narrower ranges for incomes of less than 4 million yen to reveal the income level that affected diet among low-income households. We included the choice ‘do not wish to answer’ for ethical reasons to protect the participants’ human rights. For each response, we calculated the equivalent household income by dividing the income by the square root of the number of household members( 35 ). To calculate equivalent household incomes, the midpoints of each income category( Reference Kondo, Kawachi and Hirai 36 ) – i.e. 0·5, 1·5, 2·5, 3·5, 5, 7 and 9 million yen – were used to represent incomes, respectively. We divided the income variables into three categories as follows: low (0·2236 to 2·2361 million yen; n 319), middle (2·3333 to 2·8868 million yen; n 194) and high (3·1305 to 6·3640 million yen; n 323) income groups. The numbers of participants in each of the three income groups were not equal because many participants had the same income. The median incomes for these groups were 1·75 million yen, 2·50 million yen and 3·50 million yen, respectively.
To assess educational attainment of parents, participants were asked about their highest level of education: (i) ‘middle school’; (ii) ‘high school’; (iii) ‘vocational/junior college’; (iv) ‘university/graduate school’; and (v) ‘do not wish to answer’. Demographic variables included the number of household members who shared the household income.
Confounding factors of children
Factors related to the diets of the children were determined by questionnaire and recorded information. The children completed questionnaires in the classrooms at the schools together. The variables included sex, frequency of exercise and foods allergy/food restriction. Children were asked how frequently they exercised (including outdoor play and excluding exercise class at school) per week. The children selected one of the following choices: (i) ‘5 days/week or more’; (ii) ‘3–4 days/week’; (iii) ‘1–2 days/week’; or (vi) ‘0 days/week’. Children were asked about the presence or absence of food allergies/food restrictions as part of clinical treatment, and the details. Children’s height and body weight were measured at school according to the manual for school check-up( 37 ) in September 2013. BMI was calculated as body weight divided by the square of height (kg/m2). Degree of obesity was calculated according to school health statistics in Japan( 37 ): degree of obesity (%)=[body weight (kg) – standard body weight for height (kg)]/[standard body weight for height (kg)]×100.
Dietary records of children
Children completed four consecutive days of 24 h dietary records with photographs. Two days with school lunches (weekdays) and two days without (weekend days) were recorded. Children recorded the names and amounts of all foods and beverages consumed at breakfast, lunch, dinner and other occasions. A camera, place mat with a graded scale (with which we could confirm the size of the tableware and amount of food), dietary record sheets for four days and instructions on keeping records were distributed to each child. The trained researcher explained how to record information and how to take photographs of meals to all participating children in each classroom using the same manual for dietary records. Thereafter, the children were trained to record information and take photographs during school lunch on the training day. Children were asked to put all dishes/foods on the place mat before eating, take photographs, and record the time of the meal, dish name, names of the foods, amounts of the foods, amounts of additional helpings and amounts of leftovers. Children were asked to take photographs after eating to confirm the amount of food left over. They were also asked to weigh the foods when possible and write the commercial name of the food if the food was processed. Parents/guardians were asked to support and confirm the records kept by the children. The children took school lunches on two school days. We printed the school lunch menu for each day on the record sheets.
Analysis of the dietary records
We used Excel Eiyo-Kun version 7.0, 2nd edition software (Kenpakusya, Tokyo, Japan), which is based on published standard tables of food composition in Japan 2010, to calculate energy values and nutrients. Fortified foods were considered in the present study and the nutrients recorded were thiamin, riboflavin, vitamin C, Ca and Fe. We added the amounts of nutrients when children consumed fortified foods and we were able to find the nutrient values. We compiled a manual for the analysis of dietary records, and registered dietitians identified the food name and amount of food from the records and photographs. The registered dietitians were recruited only for analysing the dietary records and were not informed about the socio-economic and demographic data of the participants. When the amounts of foods and nutrients were calculated, cooked values of foods and weights were used for cereals, rice, noodles, breads and beans (taking into account loss of nutrients during cooking), while raw values of foods and weights were used for vegetables, meat, fish and others (did not take into account loss of nutrients during cooking).
For school lunches, we obtained data on the average amounts of foods and nutrients per person on each recorded day. Photographs of standard amounts of school lunch items were taken by the school dietitian/researcher each day at each school. These photographs were compared with the photographs taken by each child when s/he had school lunches to identify the amounts of the foods eaten by each child.
Participants who completed dietary records for four days were included in the present study. School lunch is served for about 190 d/year in public primary schools in Japan( 30 ), which is roughly half of all days in a given calendar year. Thus, the levels of intake of nutrients/foods from the averages of the four days were used to assess average intake levels for one year.
Statistical analysis
Numerical variables were expressed as means and sd or 95 % CI, and as percentages of participants for categorical variables. The χ 2 test was used to compare the characteristics of the participants as categorical data among household incomes, and one-way ANOVA was used for body weight, height and degree of obesity of the children. Differences in means for nutrient and food intakes among household incomes were analysed using linear mixed-effects models, taking into account a random effect (prefecture and school), fixed effects (sex, exercise, allergy of child, father’s and mother’s educational attainments) and a covariate (BMI of child). Amounts of differences in nutrient and food intakes between low-income and middle/high-income groups adjusted for confounding factors explained above were also calculated. Tests for trend association were performed by modelling the median value of each category of household income as a continuous variable. Two-tailed P value of <0·05 was considered statistically significant. All analyses were performed using the statistical software package IBM SPSS Statistics version 21.0.
Results
Study population
The characteristics of the participants and income levels are shown in Table 1. The study population consisted of 409 boys and 427 girls. Significant differences in demographic and socio-economic variables existed among the three income groups, such as prefecture, educational attainment of fathers and mothers, household annual income and number of family members (all P<0·001), but there were no differences in sex, frequency of exercise (excluding physical education class), presence or absence of food allergy, height, body weight, BMI and degree of obesity among the children.
* P value for χ 2 test among income levels.
† P value for one-way ANOVA test among income levels.
‡ Degree of obesity (%)=[body weight (kg) – standard body weight for height (kg)]/[standard body weight for height (kg)]×100.
Associations between household income and diet: average over 4 d
Results of the analysis of associations between household income and diet given as averages over 4 d are shown in Table 2. The results showed that the low-income group consumed less protein, animal protein, energy from protein, α-tocopherol, niacin, vitamin B12, pantothenic acid, Mg, P, Fe, Zn and Cu at the nutrient level, and sugar, green vegetables and fish/shellfish at the food group level, than those of the higher income groups after adjusting for confounding factors (all P<0·05 for P value or P for trend). Energy obtained from carbohydrates was higher in the low-income group than in the higher income groups (P<0·05 for P value and P for trend). There was no difference in total energy intake levels among the income groups.
* Mean and sd of nutrients and food intake are unadjusted values.
† Mean and 95 % CI of difference from low income are values adjusted for prefecture (A, B, C, D), school (nineteen schools), sex (male and female), BMI (continuous), exercise (>5 d/week, 3–4 d/week, 1–2 d/week and 0 d/week) and food allergy (yes, no) of the child, and parental education (less than high school, high school, college, university/graduate school, don’t know).
‡ Linear mixed-effects model was used to test difference in nutrient and food intake of children among household income levels. Adjusted for the variables above.
§ A linear trend test was used with the median value in each income group as a continuous variable in regression analysis. Adjusted for the variables above.
Associations between household income and diet on the days with/without school lunch
The results of analysis of the associations between household income and diet among averages over 2 d with school lunches are shown in Table 3. The results show that there were no significant differences in the amounts of nutrients or foods consumed among the income groups.
* Mean and sd of nutrient and food intake are unadjusted values.
† Mean and 95 % CI of difference from low income are values adjusted for prefecture (A, B, C, D), school (nineteen schools), sex (male and female), BMI (continuous), exercise (>5 d/week, 3–4 d/week, 1–2 d/week and 0 d/week) and food allergy (Yes, No) of the child, and parental education (less than high school, high school, college, university/graduate school, don’t know).
‡ Linear mixed-effects model was used to test difference in nutrient and food intake of children among household income levels. Adjusted for the variables above.
§ A linear trend test was used with the median value in each income group as a continuous variable in regression analysis. Adjusted for the variables above.
The results from analysis of the associations between household income and diet from averages over 2 d without school lunches are shown in Table 4. The results showed strong associations between household income and nutrient and food consumption. The low-income group consumed less protein, animal protein, cholesterol, energy from protein, α-tocopherol, niacin, vitamin B12, pantothenic acid, K, Mg, P, Fe, Zn and Cu at the nutrient level, and sugar, green vegetables, fish/shellfish and meats/poultry at the food group level, than those in the higher income groups (all P<0·05 for P value or P for trend). The energy obtained from carbohydrates was higher in the low-income group than in the higher-income groups (P<0·05 for P value and P for trend).
* Mean and sd of nutrients and food intake are unadjusted values.
† Mean and 95 % CI of difference from low income are values adjusted for prefecture (A, B, C, D), school (nineteen schools), sex (male and female), BMI (continuous), exercise (>5 d/week, 3–4 d/week, 1–2 d/week and 0 d/week) and food allergy (yes, no) of the child, and parental education (less than high school, high school, college, university/graduate school, don’t know).
‡ Linear mixed-effects model was used to test difference in nutrient and food intake of children among household income levels. Adjusted for the variables above.
§ A linear trend test was used with the median value in each income group as a continuous variable in regression analysis. Adjusted for the variables above.
Discussion
The findings of the present study using 4 d dietary records confirm the association between household income and amounts of foods and nutrients consumed among Japanese schoolchildren. The results suggested that the lower levels of protein and micronutrients, and higher percentage of energy from carbohydrate consumed by children from low-income households compared with higher-income groups were caused by a lower consumption of protein foods and green vegetables. This means that fewer main dishes (fish, etc.) and side dishes (vegetables, etc.) were consumed, while staple foods (rice, noodles, bread and cereals) constitute the meals of children from low-income households. An additional contribution of the study is that our data showed modification of the association between household income and the foods and nutrients consumed by children on days both with and without school lunches. No significant disparities in the foods and nutrients consumed among household income groups occurred in the days with school lunch compared with those without. This finding provides evidence to show that the school lunch programme may reduce diet disparities due to different household economic status among children.
Nutritional status
Previous studies showed an inverse relationship between household income and the prevalence of obesity in Western developed countries( Reference El-Sayed, Scarborough and Galea 10 – Reference Senese, Almeida and Fath 12 ). However, there were no significant associations between income, height, body weight and degree of obesity in the present study. Kwon et al. reported a positive relationship between household income and the prevalence of obesity in Korean boys (7–12 years old)( Reference Kwon, Oh and Park 20 ). Furthermore, Kachi et al. reported that lower household incomes and expenditures were associated with a risk of being overweight in Japanese adolescents (12–18 years old), but no significant associations were found for younger children (6–11 years old)( Reference Kachi, Otsuka and Kawada 38 ). The association between household income and the prevalence of obesity among primary-school children in Asian developed countries is not clear.
Associations between household income and diet of children
Strong associations were found between household income and nutrient intake of the children in the present study. The low-income group obtained less energy from protein and higher energy from carbohydrates than the high-income group for Japanese children. The same results were observed among Korean boys; however, as household income increased, the intake of energy from carbohydrates increased and that from protein decreased in American boys( Reference Kwon, Oh and Park 20 ). Our study also showed a strong positive association between income and micronutrient intake. The same positive association is observed between SES and micronutrient intake among children in Western countries( Reference Mark, Lambert and O’Loughlin 19 , Reference Navakovic, Cavelaars and Geelen 21 ).
Our study suggested that the lower levels of nutrients, protein and micronutrients consumed by children from low-income households compared with higher-income groups were caused by differences in the types of foods consumed, and a lower consumption of fish/shellfish and green vegetables. Similar associations between household income and diet have been observed among Japanese adults. Low-income groups among Japanese adults consume less fish and vegetables and obtain more energy from carbohydrates than higher-income groups( Reference Nishi, Horikawa and Murayama 26 ), and have a higher prevalence of obesity and diabetes( Reference Horikawa, Murayama and Ota 39 ). The association between household income and diet might affect not only adult family members, but also the children of families in Japan.
The possible pathway of this positive association might be linked to food cost. Previous studies among adults showed that positive associations exist between the cost of food and energy intake from fruits and vegetables, meat, fish and eggs, and there are inverse associations between cost of food and energy intake from fats and sugar( Reference Drewnowski and Specter 40 ). Nutrient-dense foods, such as fruits, vegetables and fish, are more expensive than energy-dense foods( Reference Darmon and Drewnowski 41 ). A more in-depth study on purchasing behaviour revealed concerns about the cost of food in low-income families and this influenced their desire to purchase healthy foods( Reference Turrell and Kavanagh 42 ). Japanese adults in low-SES subgroups consume cheaper, lower-quality foods and nutrients( Reference Okubo, Murakami and Sasaki 43 ). These studies suggest that food cost would be an intermediator between household income and dietary intake of children. Improvement in economic factors, such as income and/or food costs, might play important roles in decreasing disparities in the diet at home among Japanese children.
We compared the nutrient intake levels with Japanese dietary guidelines to assess whether such dietary disparity is acceptable in our society. The Estimated Average Requirement (EAR) and RDA of protein for children aged 10–11 years in Japan are 40 and 50 g/d, respectively( 44 ). The percentage of children who consumed less protein than the EAR in low- and high-income groups was 0·6 and 0 % based on averages over 4 d, and 3·8 and 2·2 % based on averages over 2 d without school lunch, respectively (data not shown). Most of the children of both income groups consumed more protein than the RDA based on averages over 4 d; however, some of the children in the low-income group consumed less protein than the RDA on days without school lunch. The intake levels of Fe have more serious implications. The EAR and RDA values for Fe in this age group are 7·0 and 10·0 mg/d (14·0 mg/d for girls who are menstruating), respectively( 44 ). The percentage of children who consumed less Fe than the EAR in the low- and high-income groups was 65·5 and 52·6 % from averages over 4 d and 73·0 and 60·1 % from averages over 2 d without school lunches, respectively (data not shown). More than 65 % of the children in the low-income group consumed Fe levels below the EAR and the difference between the low- and higher-income groups was more than 10 %. These disparities could indicate a public health problem in Japan.
Difference in associations of income with dietary intake between days with and without school lunch
Our results showed that the association between household income and dietary intake was reduced during days with school lunch. This might be the first evidence to show a difference between days with school lunch and days without in association with household income and dietary intake in primary-school children. One reason might be that food and nutrient composition in school lunches can make up for the differences in food and nutrient intake at home. Standard amounts of nutrients and food groups must be provided for each child per meal at school lunch in Japan and school lunches provide more than one-third of the dietary requirements per day. Moreover, food eaten on non-school days and non-school meals occurred on the weekend. Weekend meals might be more sensitive to household economic status. Low-income groups might consume less nutritious foods at the weekend compared with weekdays, and the difference in diet between weekend days and weekdays might be larger than that in high-income groups. Some reports from the USA suggest that the school lunch programme for children in low-income households is an effective way to improve dietary quality and prevent obesity( Reference Huang and Bamidge 31 – Reference Nguen, Ford and Yaroch 34 ). However, a systematic review of the effects of school-based health interventions, including school lunches, on SES disparities in health concluded that the effects are unclear( Reference Moore, Littlecott and Turley 45 ). School lunch programmes would be effective in terms of food security for children in low-SES households, but their effects on dietary behaviour and improving nutrition at home have not been clarified. Further research into whether school lunches reduce dietary disparities might be useful for deciding on additional school-based health intervention policies for children.
Study limitations
Our study has several limitations. First, sampling bias should be considered. The study participants were not selected randomly from the Japanese population. We selected several public primary schools that we contacted through the local governments and invited all 5th grade students to participate. If all these students participated in the study, low-income participants would be included at almost the same rate as in the Japanese population. We also need to improve the quality of the dietary records. We need to explain how to write a dietary record and how to take a photograph for children. If the participants are scattered, it is difficult to explain and exercise the same quality. Second, the response rate to our study was high (85·1 %), but many participants opted not to answer questions about income (307 out of 1231 participants) and the final response rate was only 57·8 %. To compare the percentage distribution of each household with primary-school children based on annual incomes of representative data in fiscal year 2012( 23 ) and the present study: less than 1 million yen, 1·4 and 0·7 %; 1 to less than 2 million yen, 6·1 and 4·3 %; 2 to less than 3 million yen, 6·8 and 7·3 %; 3 to less than 4 million yen, 8·9 and 12·3 %; 4 to less than 6 million yen, 25·7 and 32·9 %; 6 to less than 8 million yen, 21·3 and 25·8 %; 8 or more million yen, 29·8 and 16·7 %, respectively. The percentage distribution of households by household income in the present study was similar to national representative data. However, the difference between the national sample and the present study was small for low income levels, and the percentage of households that have 8 or more million yen was lower in the present study. The disparity of the diet in our study would be underestimated compared with that on a national level.
Policy implications
The present study’s results have important implications for nutritional policy. They revealed independent effects of income on diet after adjusting for related variables including parental education. This suggested that economic support might effectively increase the amount of protein and other micronutrients consumed in low-income households.
Our study revealed that income-related dietary disparities in children were not significant on days with school lunches, but significant on days without. This suggests that the school lunch programme might play an important role in reducing dietary disparity in children of low SES. However, school lunches are available for only 190 d/year at Japanese public primary schools; therefore, children do not have access to school lunches for the remaining 185 d. A meal programme for children at school or in the community might be recommended during the summer and other vacations.
Conclusion
Our study confirmed an association between household income and amounts of foods and nutrients consumed by Japanese schoolchildren using 4 d dietary records. Lower income level was associated with a lower intake of fish/shellfish, green vegetables and sugar at the food group level, a lower intake of protein and several micronutrients, and a higher energy intake from carbohydrates at the nutrient level among the children.
Moreover, we found variations in the association between household income and the food/nutrients children consumed, which depended on whether school lunches were consumed or not. The disparities in foods and nutrients consumed among the income groups were not significant on days with school lunches, while they were significant on days without school lunches. The evidence suggests the importance of school lunches for reducing disparities in the diets of children from households with various incomes.
Acknowledgements
Acknowledgements: The authors thank the volunteer subjects and their families who participated in this study. They are also greatly indebted to the primary schools and the Boards of Education in the cities. They thank Professor Naohito Tanabe for his contribution to the statistical analysis. Financial support: This research was supported by a Health and Labour Sciences Research Grant (Comprehensive Research on Life Style-related Diseases Including Cardiovascular Diseases and Diabetes Mellitus (N.M.; grant number H24-H26 Jyunkankitou-Seisyu-Ippan-006)) by the Ministry of Health Labour and Welfare, Japan; and a Grant-in-Aid for Scientific Research (KAKENHI) (N.M.; grant number 15K00728) by the Japan Society for the Promotion of Science. The funders had no role in the design, analysis or writing of this article. Conflict of interest: None. Authorship: N.M. formulated the research question, designed the study, conducted the study, analysed the data and wrote the article. H.I., T.Y. and S.H. designed the study, conducted the study and analysed the data. A.N., Y.A., M.N., Y.Y. and S.S. conducted the study and analysed the data. A.A. formulated the research question and designed the study. Ethics of human subject participation: This study was conducted according to the guidelines laid down in the Declaration of Helsinki and all procedures involving human subjects/patients were approved by the ethics committee of the University of Niigata Prefecture, Japan (number 1309). Written informed consent was obtained from all subjects/patients.