Research | Open | Published:
Body composition and physical activity in Italian university students
Journal of Translational Medicinevolume 12, Article number: 120 (2014)
Increased sedentary lifestyle and prevalence of overweight/obesity are common in western countries. The purposes of this study were (i) to assess the main health-related anthropometric characteristics in a sample of students in relation to sex, amount of physical activity and sport discipline, and (ii) to investigate the accuracy of the Body Mass Index (BMI) and Waist-to-Stature Ratio (WSR) as indicators of body fat percentage (%F) in young adults.
734 university students, both sexes, participated in the present research. A self-administered questionnaire acquired socio-demographic information (sex, age) and sport participation (hours/week, sport discipline). Anthropometric measurements and grip strength values were acquired according to standardized procedures. Body composition was assessed by means of the skinfold method.
Most students had normal BMI, WSR and %F. There were significant statistical differences in all anthropometric traits between the two sexes. One-way ANOVAs within sex showed statistically significant differences in biceps skinfold, waist circumference (WC), WSR, body density (BD), %F and fat mass (FM) among different levels of physical activity in males; and in weight, BMI, arm girths and fat free mass (FFM) in females. One-way ANOVAs within sex showed statistically significant differences in arm girths, grip strength and FFM among different sport disciplines in males, and in height, weight, BMI, WC, relaxed arm girth, grip strength, FM and FFM in females. Despite the significant and positive correlation of BMI and WSR with %F both indices had poor sensitivity.
Physical activity plays an important role in body composition parameters: the most active males had the least amount of FM and the most active females had the greatest amount of FFM. BMI and WSR are not accurate indices of adiposity in young adults.
Body composition assessment is used to monitor performance and training in the athletic community, and to verify the health status of the population in general. The Body Mass Index (BMI) is often used to evaluate the weight status, even if it does not discriminate between different components of the overall body mass by definition (BMI = weight/height2). Therefore, the adoption of BMI as a predictor of adiposity and of consequent health risk should be used with caution, especially with physically active individuals, who usually have a higher body density and fat free mass (FFM) than the general population[2–4]. Body fat percentage (%F) instead is directly correlated with increased health risk, especially for metabolic and cardiovascular diseases[5–9]. Waist-to-Stature Ratio (WSR) and Waist Circumference (WC) are supposed to have greater discriminatory power compared to BMI[10, 11] and are more sensitive than BMI as an early predictor of health-related risks. In particular, WSR is probably the most sensitive anthropometric index for the screening of the metabolic syndrome in Mediterranean populations, compared to both BMI and WC.
Low levels of physical activity may place individuals at increased risk of obesity and cardiovascular diseases. On the other side, physical activity has been suggested as a means to reduce and control body fatness. More in general, regular physical activity has proved to effectively reduce diverse health risk factors, especially those related to cardiovascular diseases and the metabolic syndrome[15, 16]. In particular, the American College of Sports Medicine recommends that adults engage in at least 150 min∙wk-1 of moderate intensity cardiovascular exercise and at least 75 min∙wk-1 of vigorous intensity training, in order to maintain a sufficient level of cardio-respiratory fitness. Resistance training is also suggested 2–3 d∙wk-1. We can therefore assume that these recommendations amount for a total of more than 4 h∙wk-1 of moderate-to-intense physical activity.
The purpose of this study was to assess the main health-related anthropometric characteristics of a group of university students, in order to evaluate their relationship with quantity and type of physical activity according to sex. In particular, FFM, %F, WC, WSR, BMI and grip strength were taken into consideration. Furthermore, the accuracy of BMI and of WSR as predictors of %F was evaluated.
Participants and study design
This was a cross-sectional study carried out on a total of 734 university students, 354 females aged 21.5 ± 2.9 yrs (mean ± standard deviation) and 380 males aged 22.1 ± 3.6 yrs, of the School of Sport Science (Faculty of Medicine, University of Ferrara, Italy) who volunteered for the study. The sample was composed of North Italian students (mainly coming from the regions of Emilia Romagna and Veneto). Body image perception was previously assessed on the same sample. The research protocol was approved by the Ethic Committee for Biomedical Research of the University of Ferrara, and all participants provided written informed consent.
A questionnaire on training and physical activity patterns was administered to participants. The mean weekly amount of physical activity was 6.7 ± 4.2 hrs for males and 4.2 ± 3.8 hrs for females; 28 males (7.4% of the total male sub-sample) and 83 females (23.4% of the total female sub-sample) did not practice any sport activity.
All measures were taken in the Anthropometry Laboratory at the University of Ferrara, during the tutorials for the students of the course of Anthropometry and Ergonomics in the second year of the School of Sport Science.
Standing (H, cm) and sitting heights (SH, cm) were measured to the nearest 0.1 cm using a wall-mounted stadiometer (Magnimeter, Raven Equipment Limited, UK). Weight (W, kg) was measured to the nearest 0.1 kg using a calibrated electronic scale. BMI was calculated as W/H2 (kg/m2). Skinfold thicknesses at biceps (B Sk) and triceps (T Sk) were measured to the nearest 0.1 cm using a Lange caliper (Beta Technology Inc.). All girths (Waist Circumference WC, Contracted Arm Girth CAG, Relaxed Arm Girth RAG) were measured to the nearest 0.1 cm using a non-metallic and non-stretchable tape. WSR was calculated as WC/H.
All measurements were taken on the left side of the body, according to standardized procedures. During the anthropometric measurements, all participants were barefoot and clothed appropriately.
Left and right hand grip strength was measured to the nearest 0.5 kg by means of a Takei dynamometer (T.K.K. 5001 grip-A Takei scientific instruments Co., LTD, Japan). The highest value of two trials was recorded, after an adequate period of rest between attempts, for each hand.
Assessment of body composition
Body density (BD) was calculated using Durnin & Womersley equations with two skinfolds (biceps and triceps), according to sex and age of the student. %F was calculated from BD using Siri equation. Fat Mass (FM, kg) was calculated as (%F*W)/100 and FFM (kg) as W-FM.
Indices and classifications
According to the World Health Organization, underweight was defined as BMI < 18.5 kg/m2, normal weight as 18.5 kg/m2 ≤ BMI < 25 kg/m2, overweight as 25 kg/m2 ≤ BMI < 30 kg/m2, and obesity as a BMI ≥ 30 kg/m2. Because of the small number of students with a BMI ≥ 30 kg/m2 (only one female and 15 males), they were included in the overweight group for further elaboration. Even if there is widespread consensus on cut-points for weight status, this is not the case for what concern fatness. According to Gallagher et al., %F ≥ 20% (males) and %F ≥ 33% (females) are the cut-points adopted to define overfatness, corresponding to overweight classification using BMI in a population of young adults.
According to the National Institute for Health and Clinical Excellence guidelines, WC ≥102 cm for men and ≥88 cm for women are prerequisite risk factors for the diagnosis of the metabolic syndrome, as WSR ≥ 0.5 for both males and females.
All variables were checked for normality and logarithmically (10-based) transformed where necessary (skinfold at biceps and triceps). Results were expressed as mean ± standard deviation. Comparisons between sexes were carried out using a two-sample Student’s t-test for continuous data and a chi-square (χ2) test for categorical data.
Subsequently, both females and males were divided into 3 tertiles, according to their level of weekly physical activity: low (≤3 h∙wk-1 for females, ≤5 h∙wk-1 for males), medium (3 < h∙wk-1 < 6 for females, 5 < h∙wk-1 < 8 for males) and high (≥6 h∙wk-1 for females, ≥8 h∙wk-1 for males). One-way ANOVAs were used to assess the differences in anthropometric variables and grip strength among the 3 groups and post hoc comparisons were performed using Tukey test.
Correlation analysis between total weekly hours of physical activity and anthropometric variables was carried out.
In order to assess the anthropometric differences among subjects practicing different activities, one-way ANOVA was performed on sports with at least 10 participants: soccer, body building, basketball, swimming and volleyball in males; gymnastics, other gym activities (O.G.A.), ballet, volleyball, swimming and jogging in females. When a significant F value was obtained, post-hoc comparisons were performed by means of Tukey test.
To determine the accuracy of BMI as a measure of overfatness - and therefore of poor health status - participants were classified into one of four categories: 1) overweight and overfat (True Positive, TP), 2) overweight and normal fat (False Positive, FP), 3) normal weight and overfat (False Negative, FN), and 4) normal weight and normal fat (True Negative, TN). The sensitivity, specificity, and predictive values of BMI were calculated for each group. Sensitivity was calculated as the proportion of overfat individuals who were identified as overweight by BMI (i.e. TP/(TP + FN)). Specificity was calculated as the proportion of normal fat individuals who were identified as normal weight by BMI (i.e. TN/(TN + FP)). Positive Predictive Value (PPV) was calculated as the probability that a participant identified as overweight by BMI was truly overfat: PPV = TP/(TP + FP). Negative Predictive Value (NPV) was calculated as the probability that a participant who was identified as normal weight by BMI was normal fat: NPV = TN/(TN + FN). Test accuracy increases as the total number of FP and FN decreases.
To test the accuracy of WSR as a measure of overfatness, the same procedure was adopted, substituting the overweight category with excessive WSR.
The statistical significance was set at p < 0.05. All analyses were performed using “Statistica” for Windows, Version 11.0 (StatSoft Italia srl, Padua, Italy).
There were significant differences among all anthropometric traits between sexes (Table 1). Males were on average heavier, taller, leaner and stronger than females and had wider girths. Females had thicker skinfolds than males, as expected[20, 25], therefore they had lower BD and higher %F. 72% of males and 89% of females were normal fat, while 27.3% of males and 10% of females were overfat.
Only 4 females (1.2% of the sub-sample) had WC ≥ 88 cm and 7 males (2.0% of the sub-sample) had WC ≥ 102 cm; 5% of females and 13% of males had WSR ≥ 0.5.
BMI mean values were in the normal range according to WHO weight status categories. χ2 test proved there was a significant difference (p < 0.001) between sexes in weight status distribution. No male student was underweight, compared to 5.6% of females who fell into this category. Most males (71.7%) and females (80.9%) were normal weight. Males were more overweight (24.2%) and obese (4.2%), than females (13.2% overweight and only 0.3% obese).
ANOVAs within male sub-sample with different levels of physical activity (Table 2) show significant statistical differences in biceps skinfold, WC, WSR, BD, %F and FM, supporting the positive effects of physical activity on health-related anthropometric traits. Tukey post-hoc test shows significant differences only between the low and high level groups.
ANOVAs within female sub-sample with different levels of physical activity (Table 3) show significant statistical differences in weight, BMI, contracted and relaxed arm girths and FFM, supporting the positive effects of physical activity, particularly on FFM. Tukey post-hoc test shows significant differences between the high level group and the other two.
Statistical correlations between hours of physical activity and BMI, triceps and biceps skinfolds, WC, WSR, BD, %F and FM were significant (p < 0.05) in males, and biceps skinfold, left and right had grip strength, BD, %F and FFM in females (Table 4).
ANOVAs between sport disciplines with more than 10 participants in males (Table 5) show significant statistical differences in relaxed and contracted arm girths, left and right hand grip strength and FFM. Tukey post-hoc test shows significant differences between body building and other sports, especially soccer, for all the traits above.
Body builders had the highest BMI – similar to that of volleyball players -, arm girths, right and left hand grip strength, BD and FFM, and the lowest H, skinfold thicknesses, %F and FM. Soccer players had the lowest W, arm girths and hand grip. Volleyball players had the highest W, WC and FM, and the lowest SH. Basketball players had the highest H and SH, and the lowest BMI. Swimmers had the thickest skinfolds and the highest %F, the lowest WC, BD, and FFM.
ANOVAs between sport discipline with more than 10 participants in females (Table 6) show significant statistical differences in height, weight, BMI, WC, RAG, left and right hand grip strength, FM and FFM. Tukey post-hoc test shows significant differences between volleyball players, gymnasts and dancers for the traits above. Volleyball players had the highest H, SH, W, BMI, triceps skinfold, girths, hand grip strength, %F and FFM. Gymnasts were the shortest and lightest and had the greatest BD, the lowest SH, skinfold thickness, WC, %F, FM and FFM. Dancers had the smallest arm girths (RAG values being similar to those of gymnasts) and grip strength.
A significant positive correlation between BMI and %F was found in males (r = 0.476, p < 0.001), but it did not reach significance in basketball players (p = 0.300) and body builders (p = 0.906). In fact, one third of the subjects who were classified as overweight according to BMI, but who were actually normal fat, practiced body building. Twelve percent of total participants fell within the FP quadrant and 10% in the FN one (Figure 1[a]). Sensitivity was 0.62 and specificity was 0.83, while PPV was 0.58 and NPV was 0.85. A significant positive correlation between BMI and %F was found in females (r = 0.622, p < 0.001) but it did not reach significance in gymnasts (p = 0.752). Seven percent were classified as FP and 4% as FN (Figure 1[b]). Sensitivity was 0.59 and specificity was 0.92, while PPV was 0.45 and NPV was 0.95. Therefore, sensitivity was poor for both sexes, reflecting the fact that the individuals who were at the same time classified as overfat (according to their %F) and overweight (according to their BMI) were only a small proportion of those who were actually overfat. Also PPV of BMI was poor, because really fat individuals were about a half of those who were classified as overweight.
A significant positive correlation between WSR and %F was found in males (r = 0.439, p < 0.001). Four percent of total participants fell in the FP quadrant and 17% in the FN one (Figure 1[c]). Sensitivity was 0.36 and specificity was 0.95, while PPV was 0.73 and NPV was 0.80. A significant positive correlation between WSR and %F was found in females (r = 0.527, p < 0.001). Three percent of total participants fell within the FP quadrant and 8% in the FN one (Figure 1[d]). Sensitivity was 0.24 and specificity was 0.97, while PPV was 0.47 and NPV was 0.92. Therefore, sensitivity was poor for both sexes, and PPV was poor especially in females.
In the present study, we found a different trend in the two sexes in relation to training volume: female students performing more hours of weekly physical activity had a significantly higher amount of FFM compared to the less active individuals, while male students showed a lower %F and FM. A study by Westerterp et al. found a negative correlation between energy expenditure and %F in males, but not in females. Also, a negative correlation between physical activity and FM was found in males, but not in females. Even if FM can be reduced by means of increased physical activity, females seem to compensate for excess energy expenditure with added energy intake. Since women tend to preserve their energy balance more than men, FM loss can be not significant. Increased caloric intake could also justify added FFM, as in the present research.
The different behaviours may be consistent with both sex-related differences and sport preferences. The examined females are more often than males engaged in individual sports and in disciplines with a relevant aesthetic component (gym activities, ballet, gymnastics). Males are more often than females engaged in team sports (soccer, basketball, volleyball) or in strength-related activities, like body building, both involving intense repeated efforts, which have been positively correlated to fat loss[29, 30]. So the different adaptation could be sport-related. Body building determines an evident increase on muscle hypertrophy, which is significant in comparison to soccer. This fact may contribute to the limited accuracy of BMI as an index of body fatness and general health status, since body builders have a high BMI, close to overweight, even if they have the lowest %F in our sample. Also, females are less physically active, on average, therefore it can be hypothesized that physical adaptation in response to moderate physical activity can be correlated to increased muscle mass, and, vice-versa, that physical adaptation in response to high volume of weekly physical activity can be correlated to reduced %F.
The variance in weekly hours of physical activity within the sample determined significant differences in body composition, and showed the limits of BMI and WSR as indices of adiposity. Intersecting BMI values with %F, we have obtained important indications on its limited applicability in a sample of young adults with different levels of weekly training hours.
The analysis of specificity and sensitivity showed that neither BMI nor WSR can be considered accurate indices of the health status of the population of young adults because they are not consistent measures of body fatness. In fact, both BMI and WSR had good specificity versus %F, but low sensitivity, suggesting that a significant percentage of overfat individuals were classified as normal according to BMI or WSR.
A possible limiting factor of the present study is that physical activity assessment (weekly training hours and type of sport) was based only on a self-reported questionnaire. Also, the training volume does not account for training intensity and quality (mainly aerobic, anaerobic etc.). A lower volume of weekly training hours involving a strenuous practice may have more significant outcomes than a higher volume with a less intense effort, in particular for what concern body composition. Moreover, it must be highlighted that the skinfold-thickness technique is an indirect method for assessing body composition, based on population-specific predictive equations. Although relatively inexpensive, non-invasive, and widely used in sportsmen, its accuracy cannot be granted especially in individuals with adipose tissue that is not well separated from the underlying muscle.
This study examined a large sample of Italian university students from the same geographical area by means of rigorous anthropometric procedures.
In conclusion, a different behaviour was highlighted in the two sexes in relation to weekly amount of physical activity: males mainly showed a decrease in %F, whereas females showed an increase in FFM, which could be explained by a stronger tendency to maintain energetic balance or by different sport preferences.
BMI and WSR have been suggested as indirect measures of %F, because of the ease with which they can be collected. The present study confirms their low accuracy. In fact, in females, misclassification (FP + FN) was 11% for both BMI and WSR. In males, misclassification was 22% for BMI and 21% for WSR. Therefore, regardless of the fact that WSR has been proposed as a better index of adiposity than BMI, both indices show similar low accuracy and they cannot be considered reliable predictors of body fatness, especially in young males. Greater accuracy can be found in females, possibly because of lower overall FFM compared to males. In fact, high FFM contributes to increased BMI, without any real detrimental effect (e.g. in body builders).
The present study confirms that an active lifestyle, including regular weekly physical activity, is significantly correlated to body composition parameters.
De Lorenzo A, Bianchi A, Maroni P, Iannarelli A, Di Daniele N, Iacopino L, Di Renzo L: Adiposity rather than BMI determines metabolic risk. Int J Cardiol. 2013, 166 (1): 111-117. 10.1016/j.ijcard.2011.10.006.
Zaccagni L, Onisto N, Gualdi-Russo E: Biological characteristics and ageing in former elite volleyball players. J Sci Med Sport. 2009, 12 (6): 667-672. 10.1016/j.jsams.2008.06.006.
Barbieri D, Zaccagni L, Cogo A, Gualdi Russo E: Body composition and somatotype of experienced mountain climbers. High Alt Med Biol. 2012, 13 (1): 46-50. 10.1089/ham.2011.1062.
Klungland Torstveit M, Sundgot-Borgen J: Are under- and overweight female elite athletes thin and fat? A controlled study. Med Sci Sports Exerc. 2012, 44 (5): 949-957. 10.1249/MSS.0b013e31823fe4ef.
Tanaka S, Togashi K, Rankinen T, Pérusse L, Leon AS, Rao DC, Skinner JS, Wilmore JH, Bouchard C: Is adiposity at normal body weight relevant for cardiovascular disease risk?. Int J Obes Relat Metab Disord. 2002, 26 (2): 176-183. 10.1038/sj.ijo.0801880.
Cho YG, Song HJ, Kim JM, Park KH, Paek YJ, Cho JJ, Caterson I, Kang JG: The estimation of cardiovascular risk factors by body mass index and body fat percentage in Korean male adults. Metabolism. 2009, 58 (6): 765-771. 10.1016/j.metabol.2009.01.004.
Onisto N, Teofoli F, Zaccagni L, Gualdi Russo E: Anthropometric traits and aging: a cross-sectional survey in diabetic elderly women. Arch Gerontol Geriatr. 2009, 48 (2): 197-200. 10.1016/j.archger.2008.01.007.
Gokulakrishnan K, Deepa M, Monickaraj F, Mohan V: Relationship of body fat with insulin resistance and cardiometabolic risk factors among normal glucose-tolerant subjects. J Postgrad Med. 2011, 57 (3): 184-188. 10.4103/0022-3859.85200.
Chuang HH, Li WC, Sheu BF, Liao SC, Chen JY, Chang KC, Tsai YW: Correlation between body composition and risk factors for cardiovascular disease and metabolic syndrome. Biofactors. 2012, 38: 284-291. 10.1002/biof.1027.
Gualdi-Russo E, Zironi A, Dallari GV, Toselli S: Migration and Health in Italy: A Multiethnic Adult Sample. J Travel Med. 2009, 16 (2): 88-95. 10.1111/j.1708-8305.2008.00280.x.
Ashwell M, Gunn P, Gibson S: Waist-to-height ratio is a better screening tool than waist circumference and BMI for adult cardiometabolic risk factors: systematic review and meta-analysis. Obes Rev. 2012, 13 (3): 275-286. 10.1111/j.1467-789X.2011.00952.x.
Ashwell M, Hsieh SD: Six reasons why the waist-to-height ratio is a rapid and effective global indicator for health risks of obesity and how its use could simplify the international public health message on obesity. Int J Food Sci Nutr. 2005, 56 (5): 303-307. 10.1080/09637480500195066.
Mombelli G, Zanaboni AM, Gaito S, Sirtori CR: Waist-to-height ratio is a highly sensitive index for the metabolic syndrome in a Mediterranean population. Metab Syndr Relat Disord. 2009, 7 (5): 477-484. 10.1089/met.2008.0101.
Sacheck JM, Kuder JF, Economos CD: Physical fitness, adiposity, and metabolic risk factors in young college students. Med Sci Sports Exerc. 2010, 42 (6): 1039-1044.
Reimers CD, Knapp G, Reimers AK: Does physical activity increase life expectancy? A review of the literature. J Aging Res. 2012, 2012: 243958-doi:10.1155/2012/243958. Epub 2012 Jul 1
Wagner A, Dallongeville J, Haas B, Ruidavets JB, Amouyel P, Ferrières J, Simon C, Arveiler D: Sedentary behaviour physical activity and dietary patterns are independently associated with the metabolic syndrome. Diabetes Metab. 2012, 38 (5): 428-435. 10.1016/j.diabet.2012.04.005.
Garber CE, Blissmer B, Deschenes MR, Franklin BA, Lamonte MJ, Lee IM, Nieman DC, Swain DP, American College of Sports Medicine: Quantity and quality of exercise for developing and maintaining cardiorespiratory, musculoskeletal and neuromotor fitness in apparently healthy adults: guidance for prescribing exercise. Med Sci Sports Exerc. 2011, 43 (7): 1334-1359. 10.1249/MSS.0b013e318213fefb.
Zaccagni L, Masotti S, Donati R, Mazzoni G, Gualdi-Russo E: Body image and weight perceptions in relation to actual measurements by means of a new index and level of physical activity in Italian university students. J Transl Med. 2014, 12: 42-10.1186/1479-5876-12-42.
Weiner JS, Lourie JA: Practical Human Biology. 1981, London: Academic Press
Durnin JV, Womersley J: Body fat assessed from total body density and its estimation from skinfold thickness: measurements on 481 men and women aged from 16 to 72 years. Br J Nutr. 1974, 32 (1): 77-97. 10.1079/BJN19740060.
Siri WE: The gross composition of the body. Adv Biol Med Phys. 1956, 4: 239-280.
James PT, Leach R, Kalamara E, Shayeghi M: The worldwide obesity epidemic. Obes Res. 2001, 9: S228-S233. 10.1038/oby.2001.123.
Gallagher D, Heymsfield SB, Heo M, Jebb SA, Murgatroyd PR, Sakamoto Y: Healthy percentage body fat ranges: an approach for developing guidelines based on body mass index. Am J Clin Nutr. 2000, 72 (3): 694-701.
McNeil B, Keeler E, Adelstein J: Primer on certain elements of clinical decision making. N Engl J Med. 1975, 293: 211-215. 10.1056/NEJM197507312930501.
Gualdi Russo E, Gruppioni G, Gueresi P, Belcastro MG, Marchesini V: Skinfolds and body composition of sports participants. J Sports Med Phys Fitness. 1992, 32: 303-313.
Westerterp KR, Goran MI: Relationship between physical activity related energy expenditure and body composition: a gender difference. Int J Obes Relat Metab Disord. 1997, 21: 184-188.
Westerterp KR, Meijer GA, Kester AD, Wouters L, Ten Hoor F: Fat-free mass as a function of fat mass and habitual activity level. Int J Sports Med. 1992, 13 (2): 163-166. 10.1055/s-2007-1021249.
Westerterp KR, Meijer GA, Janssen EM, Saris WH, Ten Hoor F: Long-term effect of physical activity on energy balance and body composition. Br J Nutr. 1992, 68: 21-30. 10.1079/BJN19920063.
Tremblay A, Simoneau JA, Bouchard C: Impact of exercise intensity on body fatness and skeletal muscle metabolism. Metabolism. 1994, 43 (7): 814-818. 10.1016/0026-0495(94)90259-3.
Sijie T, Hainai Y, Fengying Y, Jianxiong W: High intensity interval exercise training in overweight young women. J Sports Med Phys Fitness. 2012, 52 (3): 255-262.
Selkow NM, Pietrosimone BG, Saliba SA: Subcutaneous thigh fat assessment: a comparison of skinfold calipers and ultrasound imaging. J Athl Train. 2011, 46: 50-54. 10.4085/1062-6050-46.1.50.
We would like to thank all the students who volunteered for the study.
The authors declare that they have no competing interests.
EG-R and LZ have planned the study, ZL and BD participated to the data collection and did all the statistical analysis. All authors were involved with data interpretation, critical revisions of the paper. All authors read and approved the final manuscript.