The effect of applying ethnicity-specific spirometric reference equations to Asian migrant workers in Korea
Annals of Occupational and Environmental Medicine volume 27, Article number: 14 (2015)
Asian migrant workers in Korea have various ethnicities. The purpose of this study was to evaluate the difference in spirometric interpretation made using the set of third National Health and Nutrition Examination Survey (NHANES III) reference equations and the relevant ethnicity-specific reference sets.
Spirometry was performed on 166 migrant and 498 Korean male workers between March and November 2012. We analyzed the spirometric data of healthy never-smokers. Spirometric patterns were evaluated using the NHANES III reference set and some relevant ethnicity-specific equations (Eom’s equation for Koreans, Ip’s equation for East Asians, Crapo’s equation for Central Asians, Memon’s equation for South Asians, and Gnanou’s equation for Southeast Asian people).
In all migrant groups except the Central Asian group, the forced expiratory volume in 1 second percentage (FEV1%) and forced vital capacity percentage (FVC%) calculated using each of the ethnicity-specific reference equations considered were significantly higher than those calculated using the NHANES III reference set. This study showed that in the evaluation of the spirometric result of subjects from Southeast Asia or South Asia, the percentage of cases with an abnormal FEV1 or FVC increased when the NHANES III set of equations was used as compared to when the ethnicity-specific equations were used.
We found that the spirometric results of all ethnic groups were higher and the composition of the abnormal spirometric result was lower when the ethnicity-specific reference equations were used instead of the NHANES III reference set.
The Occupational Safety and Health Act in Korea requires that a spirometric test be done as a regular health checkup for employees exposed to pulmonary hazards in the workplace, such as dust, metal fumes, oil mist, or organic solvents that can lead to some ventilatory disorders like occupational asthma, chronic obstructive pulmonary disease (COPD), and interstitial lung disease. Some parameters measured during a spirometric test include forced vital capacity (FVC), forced expiratory volume in 1 second (FEV1), FEV1/FVC, and peak expiratory flow (PEF) . A spirometric test is an effective and reliable diagnostic tool for assessing obstructive ventilatory disorders such as asthma or COPD by using FEV1/FVC and PEF, and restrictive ventilator disorders such as lung parenchymal disease by using FVC. For the interpretation of a spirometric test, a normal pulmonary function reference value is necessary for determining the normal range, and this value is calculated using spirometric reference equations. The job fitness of employees and work-relatedness of reduced pulmonary function are judged with a simultaneous consideration of the spirometric result, work-related symptoms, and workplace respiratory hazards [2,3].
The number of migrant workers in Korea has increased. In 2012, 1,445,103 non-Korean citizens were employed in Korea and 230,237 of them had an E-9 visa for non-professional work . China and Southeast Asia are the most common regions of origin of the migrant laborers in Korea and account for 79% of the non-Koreans in Korea. Further, since many immigrants, mainly from Southeast Asia, remain undocumented, the actually ratio of migrant workers from Asian countries is known to be higher than the calculated value . Considering the diversity of language families in Asia, it is supposed that Asians have a number of ethnic origins .
Respiratory impairment is influenced by ethnic differences . The variation of ethnic spirometric reference equations has been previously reported in the literature [8-10], and official recommendations have been made. The most recent guidelines, published in 2005 by the American Thoracic Society (ATS)/European Respiratory Society (ERS), recommend ethnicity-specific reference standards. The ERS/Global Lung Initiative task force announced spirometric reference equations (ERS/GLI 2012) derived from data collected from healthy nonsmokers in the age group of 3–95 years from 33 countries. The ERS/GLI 2012 equations provided multi-ethnic values and the lower limit of normal (LLN) for spirometry [11,12]. However, many countries still have no choice but to use pulmonary function reference equations that have not been developed for the local ethnic groups. Most spirometric tests use European descent as the base population for a prediction equation of normal reference values . Therefore, a spirometric interpretation, that is, a spirometric reference equation developed considering the ethnicity of subjects, is needed [14,15]. Spirometric values of Asians in general are known to be lower than those of Caucasians, but spirometric values of Koreans in particular are not lower than those of Caucasians [16-19].
In this study, we compared the spirometric results of migrant workers and Korean workers in Busan metropolitan city and compared the interpretations of the spirometric data determined by ethnicity-specific reference equations with those determined by the third National Health and Nutrition Examination Survey (NHANES III) reference equation set.
Materials and methods
The subjects in this study were 3,067 male workers who underwent spirometry as part of a regular health examination in the occupational health examination center of a hospital located in Busan between March and November 2012. All subjects were Korean and non-ethnic-Korean migrant workers. In this study, we considered 1658 subjects in the age group of 20 to 40 years because most migrant workers belong to this age group. We excluded 24 subjects (22 Korean workers and 2 migrant workers) with an abnormal chest radiograph or unconfirmed nationality, or with a history of significant lung or heart disease, or with respiratory symptom such as persistent cough, persistent phlegm production, etc. Subjects with a satisfactory result with respect to suitability and reproducibility were selected for inclusion by analyzing the spirometry data using the technical guidelines of the spirometric test and interpretation of the Korea Occupational Safety and Health Agency . 166 migrant workers and 1,434 Korean workers showed satisfactory results. Three Korean subjects were matched with each migrant subject randomly selected with respect to age and smoking history. Finally, a total of 498 Korean and 166 migrant workers were enrolled.
Ethnicity-specific spirometric equations
The migrant workers were classified into the following four groups on the basis of the subjects’ country of origin: East Asia, Southeast Asia, Central Asia, and South Asia. East Asia consisted of China, Japan, and Taiwan. Southeast Asia consisted of Brunei, Cambodia, Indonesia, Laos, Malaysia, Myanmar, Philippines, Singapore, Thailand, and Vietnam. Central Asia included Kazakhstan, Kyrgyzstan, Tajikistan, Turkmenistan, Uzbekistan, and Mongolia. South Asia included Afghanistan, Bangladesh, Bhutan, India, Maldives, Nepal, Pakistan, and Sri Lanka.
We selected an ethnicity-specific reference equation for each of the abovementioned subregions on the basis of a literature review. We searched for studies of spirometric reference equations that had been established for each of the subregions and for Korea that met the following criteria: First, spirometry had been performed on the basis of the ATS guidelines. Second, spirometric equations for the FEV1 and FVC had been devised. Third, the studies needed to have been based on subjects belonging to the age group of 20 to 40 years. Fourth, the number of subjects was higher than 500. Fifth, the studies must have been published in recent years, at least in or after 1999. We could not find the equations developed for Central Asia that met the fourth criterion, and therefore, the selected study for Central Asia had 344 subjects. The selected spirometric reference equations of each subregion were as follows: Eom’s equation for Koreans, Ip’s equation from Hong Kong for East Asians, Gnanou’s equation from Malaysia for Southeast Asians, Crapo’s equation from Mongolia for Central Asians, and Memon’s equation from Pakistan for South Asians. Although the most limited definition of Central Asia does not include Mongolia, no other spirometric reference has been developed recently in Central Asia strictly defined. The former definition of Central Asia by United Nations Educational, Scientific, and Cultural Organization (UNESCO) included Mongolia . The reference equation used for the group of Central Asians was developed in Mongolia [9,21-25] (Table 1).
The subjects’ basic data such as nationality, workplace, job position, smoking history, age, height, and weight were collected. A chest radiograph was performed. The participants were asked about their smoking history through a self-administered questionnaire before the spirometry.
For the spirometry, FVC and FEV1 were measured. Spirometry was performed in subjects in stable condition after explaining its purpose and the required position. The spirometry was performed with the examinee in a seated position after one practice test. The better result of three attempts was selected. A FlowScreen testing device (Cardinal Health, Inc., Germany) was used for the screening. Standardized pulmonary function tests were performed according to the ATS/ERS guidelines .
Each migrant worker’s country of origin was stratified into four groups: East Asia, Southeast Asia, Central Asia, and South Asia. By applying the NHANES III reference set, we compared the spirometric results of each group with those of the Korean workers. Then, the results of applying ethnicity-specific spirometric references were compared with those of applying the NHANES III reference set. To acquire values with percentages, the actual measured FEV1 and FVC were divided into the calculated FEV1 and FVC from each equation (FEV1% = measured FEV1/calculated FEV1, FVC% = measured FVC/calculated FVC). Abnormal FVC and FEV1 values were defined as values that were less than 80% of the reference value . Subjects were divided into three groups according to their smoking status: never smoker (n = 183 for Korean workers, n = 88 for all migrant workers), ex-smoker (n = 56 for Korean workers, n = 11 for all migrant workers), and current smokers at the time the spirometry tests were conducted (n = 259 for Korean workers, n = 67 for all migrant workers). As cigarette smoke affects lung function and may be considered a potentially confounding factor, current and ex-smokers were excluded from the statistical analysis.
The differences in the values of the FEV1% and FVC% between results obtained by applying NHANES III and the ethnicity-specific equations for the never smokers were analyzed using a paired t-test. The differences in the values of the FEV1% and FVC% between never-smoker Korean and migrant workers obtained by applying the NHANES III set of equations and the ethnicity-specific sets were analyzed using a t-test. The differences in the frequency of abnormal FEV1 and FVC were analyzed with McNemar’s test. The evaluation agreement for abnormal FVC and FEV1 obtained using the reference equations and the concordance coefficients were analyzed using the Kappa distribution. IBM SPSS Statistics version 21.0 (IBM, Corp., Armonk, NY, USA) was used for the analysis.
Among the 166 migrant workers grouped by region of origin, the group of Southeast Asians was the largest, with 109 subjects. The smallest group was that of East Asians with 11 subjects. Among individual countries, 55 subjects were from Indonesia; on the other hand, there was only one subject from Kazakhstan and one from Taiwan. Most of the subjects, 476 Korean workers (95.6%) and 165 migrant workers (99.4%), were working in the manufacturing sector. The mean age of the migrant workers was 29.0 ± 5.3 years, and the mean age of the Korean workers was 29.3 ± 5.3 years. The highest mean age group was that of the East Asians, with the mean age of 33.9 ± 4.8 years. The lowest mean age group was that of the Southeast Asians, with the mean age of 28.1 ± 5.3 years. The Koreans were significantly taller and heavier than the migrants. The proportion of never smokers was 36.7% in Korean workers and 53% in migrant workers (Figure 1, Table 2).
Comparing the results of applying NHANES III and ethnicity-specific reference equations to the never-smoker group
In all the migrant groups except the Central Asian group, the FEV1% and FVC% values calculated using the ethnicity-specific equations were significantly higher than those calculated using the NHANES III reference set. The value of the FEV1% of the Central Asian group calculated using the NHANES III reference set and that calculated using the corresponding ethnicity-specific set did not differ significantly. Some differences between the spirometric measurements calculated using both the equations were more than 20%; these are the differences in FVC% in the Southeast Asian group and in FEV1% in the South Asian group. The FEV1% and FVC% values calculated using the NHANES III reference set were higher for Korean workers than for all the migrant workers in aggregate (p = 0.000 for FEV1%, p = 0.000 for FVC%). In contrast, the result of the FEV1% and FVC% calculated using the ethnicity-specific equations for Korean workers were lower than those for all the migrant workers (p = 0.015 for FEV1%, p = 0.000 for FVC%) (Tables 3 and 4).
In the evaluation of the normal or abnormal FEV1 and FVC, the concordance coefficient between applying NHANES III and applying ethnicity-specific equations in the case of the never-smoker Korean workers was 0.378 and 0.433, respectively, which implies moderate agreement. In the case of the never-smoker migrant workers, the concordance coefficient between applying NHANES III and applying ethnicity-specific equations to evaluate the FEV1 and FVC was 0.095 and 0.159, respectively, which implies poor agreement. In the evaluation of the abnormal ratio of FEV1 and FVC, the abnormal ratio obtained by applying NHANES III was greater than that obtained by applying the ethnicity-specific equations for Korean and all migrants except those of East Asian and Central Asian origins (Tables 5 and 6).
In this study, we made an assumption that Asian migrant workers from the same subregion belonged to the same ethnic classification. Terms such as ethnicity or race are used for the classification of human beings. Ethnicity involves belonging to cultural, sociological aspect of region of origin, and the race classification is based on physical characteristics . However, there are genetic differences between ethnic groups or races; therefore, clinical practices involving racial classification should be carried out carefully . For an accurate spirometric interpretation, it is necessary to know the ethnicity of the subjects tested . Some countries in Asia have a population composed of several ethnic groups. For example, the Singaporean population consists of the major ethnic groups of Chinese, Malays, Indians, and Eurasians . While the country of origin cannot be used to determine race or genetic characteristics with accuracy, it can be used for representing ethnicity.
We could not find reference equations for all the countries of origin of the subjects of this study. Therefore, reference equations developed for people from the same subregions were used. We excluded reference equations developed before 1999 in the selection of equations developed for individuals from a similar spirometric guidelines or socioeconomic status of subjects. Although the ethnicity-specific reference equations applied in this study were developed in Asia, there are some differences among these equations.
Our results demonstrate that in almost all migrant groups, the FEV1% and FVC% values calculated using the ethnicity-specific equations and the NHANES III reference set were significantly different. The FEV1% and FVC% of Southeast Asians calculated using the ethnicity-specific equation were higher than those calculated using the NHANES III reference set by more than 19%. Further, by comparing the pulmonary function classification of individuals as “normal” or “abnormal” using the ethnicity-specific equations and the commonly used NHANES III set of equations, we found that the proportion of cases with an abnormal FEV1 and FVC increased when using the NHANES III set of equations as compared to when using the ethnicity-specific sets. This finding agrees with previous data from other immigrant studies performed in the United States and Israel [30,31]. Further, in a study conducted in India, when the spirometric reference of North Indians was applied to South Indians, the proportion of abnormal values increased . The proportion of abnormal spirometric values of the East Asian and Central Asian groups did not differ significantly between when the NHANES III set of reference equations was applied and when the ethnicity-specific sets were applied. This could be attributed to the limited number of subjects or to the ethnic characteristics of the subjects.
Morris equation  and Choi’s equation , as well as NHANES III set of reference equation are also widely used as the references for interpretation of spirometry results in Korean. So, although not shown earlier as a result, we compared two reference equations with the relevant ethnicity-specific equations. We found that all the migrant groups except the Central Asian group, the FEV1% and FVC % values calculated using the relevant ethnicity-specific equations were significantly higher than those calculated using the Morris equation and Choi’s equation. In the evaluation of the spirometric results of subjects from Southeast Asia or South Asia, the percentage of abnormal FEV1 increased when the Morris equation and Choi’s equation were used instead of the relevant ethnicity-specific equations. These results were very similar to compared NHANES III reference set with the relevant ethnicity-specific equations.
Ethnic origin affects pulmonary prediction values. The ERS/GLI 2012 equation provides age, height, sex, and ethnicity-specific reference equations and the LLN for spirometry . However, thus far, no definite valid explanation for differences in the normal lung function parameters attributed to ethnicity has been reported. Nevertheless, the literature suggests that inherited differences, anthropometric differences, environmental and social factors among ethnic groups play a role .
In 2012, in Korea, there were 128,274 migrant male workers from Southeast Asia, 47,552 from South Asia, 22,773 from Central Asia, and 6,199 from East Asia. With respect to the country of origin, the largest number of migrant workers (52,084) were from Vietnam, followed by 24,701 from Indonesia, and then, 19,830 from Sri Lanka, 15,116 from Uzbekistan, and 14,619 from Nepal. Similarly, according to the number of migrant workers, our subject groups were from Southeast Asia, South Asia, Central Asia, and East Asia . Pulmonary function of the migrant workers are currently being calculated using the same equation as that is used for Koreans. Therefore, we should pay attention to their normal lung function predictive values.
This study has several limitations and possible confounders. Firstly, all subjects were employees undergoing a routine annual health examinations at a health exam center located in Busan city. This may have caused a selection bias due to the “healthy worker effect”. Further, the material hazardous to the pulmonary system by occupational exposure was not considered. However, the migrants’ working conditions were almost all similar. Another possible limitation is the fact that the selection of subjects was not representative of all migrant workers in Korea. Nevertheless, we think that the results of this study should be typical of spirometry results from migrant workers throughout Korea.
This study showed that the spirometric results of migrant workers from Southeast Asia or South Asia can be calculated as too low when obtained using the same equation as that used for Koreans. Laws requiring regular health check-ups of workers would support the health management of workers exposed to hazardous materials at the workplace. However, an overdiagnosis of the pulmonary dysfunction of migrant workers can influence the assessment of their work fitness and consequently, limit the occupational options of migrant workers. Therefore, an appropriate standard for evaluating the spirometric results of migrant workers should be established.
The number of migrant workers in Korea is increasing. According to Korean law, employees exposed to pulmonary hazards in the workplace must undergo regular health check-ups including spirometry. This study showed that all the migrant groups except the Central Asian group, the FEV1% and FVC% values calculated using the ethnicity-specific equations were significantly higher than those calculated using the NHANES III reference set. And in the evaluation of the spirometric results of subjects from Southeast Asia or South Asia, the percentage of abnormal FEV1 or FVC increased when the NHANES III reference set was used instead of the relevant ethnicity-specific equations.
Ministry of Employment and Labor: Enforcement regulation of the occupational safety and health act. http://english.kosha.or.kr/english/cmsTiles.do?url=/cms/board/board/Board.jsp?communityKey=B0488&menuId=1240&searchType=ALL&searchWord=&pageNum=1&pageSize=&boardId=29&act=VIEW (2012). Accessed 30 April 2015
Korea occupational Safety & Health Agency: Technical guideline of spirometric test and interpretation. http://www.kosha.or.kr/cms/generate/FileDownload.jsp?content_id=358431&category_id=4828&version=1.0&file_name=358431_2014_2.pdf (2014). Accessed 30 April 2015. (Korean, translated by Kim NM).
The National Institute for Occupational Safety and Health (NIOSH). Spirometry. http://www.cdc.gov/niosh/topics/spirometry/spirola.html (2012). Accessed 11 October 2014.
Korean Statistical Information Service: Populaltion of registered foreigners by nationality or types of visa. http://kosis.kr/statHtml/statHtml.do?orgId=111&tblId=DT_1B040A5&conn_path=I2 (2014). Accessed 15 October 2014. (Korean, translated by Kim NM).
Korea Institute for International Economic Policy. Southeast Asian Labor Migration to Korea: Origin-Country Factors and Policy Implications. http://www.kiep.go.kr/include/filedown.jsp?fname=WEU-13-35.pdf&fpath=EngWorld&NO=58&FNO=49 (2014). Accessed 30 April 2015.
Office for Coordination of Humanitarian Affairs. Major Language Families in Asia-Pacific. http://reliefweb.int/sites/reliefweb.int/files/resources/OCHA_ROAP_Language_v6_110519.pdf (2011). Accessed 19 October 2014.
Vaz Fragoso CA, McAvay G, Gill TM, Concato J, Quanjer PH, Van Ness PH. Ethnic differences in respiratory impairment. Thorax. 2014;69:55–62.
Fulambarker A, Copur AS, Javeri A, Jere S, Cohen ME. Reference values for pulmonary function in Asian Indians living in the united states. Chest. 2004;126:1225–33.
Ip MS, Ko FW, Lau AC, Yu WC, Tang KS, Choo K, et al. Updated spirometric reference value for adult Chinese in hong Kong and implications on clinical utilization. Chest. 2006;129:384–92.
Torre-Bouscoulet L, Perez-Padilla R. Adjustment of several spirometric reference equations to a population-based sample in Mexico. Salud Publ Mex. 2006;48:466–73.
Miller MR, Hankinson J, Brusasco V, Burgos F, Casaburi R, Coates A, et al. Standardisation of spirometry. Eur Respir J. 2005;26:319–38.
Quanjer PH, Stanojevic S, Cole TJ, Baur X, Hall GL, Culver BH, et al. ERS global lung function initiative. Multi-ethnic reference values for spirometry for the 3 e 95-yr age range: the global lung function 2012 equations. Eur Respir J. 2012;40:1324–43.
Ghio JA, Crapo RO, Elliot G. Reference equations used to predict pulmonary function. Reference equations used to predict pulmonary function. Survey at institutions with respiratory disease training programs in the united states and Canada. Chest. 1990;97:400–3.
Pellegrino R, Viegi G, Brusasco V, Crapo RO, Burgos F, Casaburi R, et al. Interpertative strategies for lung function tests. Eur Respir J. 2005;26:948–68.
Park JH, Lee JS, Huh JH, Oh YM, Lee SD, Lee SW, et al. Proposal of new criteria for assessing respiratory impairment. Tuberc Respir Dis. 2011;70:199–205. Korean with English abstract.
Hankinson JL, Kawut SM, Shahar E, Smith LJ, Stukovsky KH, Barr RG. Performance of american thoracic society-recommended spirometry reference values in a multiethnic sample of adults: the multi-ethnic study of atherosclerosis (MESA) lung study. Chest. 2010;137:138–45.
Sharp DS, Enright PL, Chiu D, Burchfiel CM, Rodriguez BL, Curb JD. Reference values for pulmonary function tests of Japanese-American men aged 71 to 90 years. Am J Respir Crit Care Med. 1996;153:805–11.
Lee CH, Lee JY, Jang EJ, Park KY, Han HW. New predictive equations for spirometric reference values and comparison with Morris equation in a Korean population. Respirology. 2008;13:365–71.
Eom SY, Moon SI, Yim DH, Lee CH, Kim GB, Kim YD, et al. Goodness-of-fits of the spirometric reference values for Koreans and USA Caucasians to spirometry data from residents of a region within chungbuk province. Tuberc Respir Dis. 2012;72:302–9. Korean with English abstract.
Eom GH. Research trend about central Asia in Korea. Asia Review. 2012;2:121–52. Korean, translated by Kim NM.
Hankinson JL, Odencrantz JR, Fedan KB. Spirometric reference values from a sample of the general U.S. Population. Am J Respir Crit Care Med. 1999;159:179–87.
Eom SY, Kim H. Reference values for the pulmonary function of Korean adults using the data Korea national health and nutrition examination survey IV (2007–2009). J Korean Med Sci. 2013;28:424–30.
Gnanou J, Caszo B, Mohamad WHW, Nawawi H, Yusoff K, Ismail T. Prediction equation for lung function in health, life time never-smoking Malaysian population. Southeast Asian J Trop Med Public Health. 2011;42:965–76.
Crapo RO, Jensen RL, Oyunchimeg M, Tsh T, Schmidt CD. Differences in spirometry reference values: a statisca comparison of a Mongolian and a Caucasian study. Eur Respir J. 1999;13:606–9.
Memon MA, Sandila MP, Ahmed ST. Spirometric reference values in healthy, non-smoking, urban Pakistani population. JPMA. 2007;57:193–5.
Ministry of Employment and Labor. Act on the prevention of pneumoconiosis and protection, etc., of Pneumoconiosis workers. http://www.ilo.org/dyn/natlex/docs/ELECTRONIC/95009/111687/F-964324601/KOR95009.pdf (2010). Accessed 30 April 2015.
Sankar P. MEDLINE definitions of race and ethnicity and their application to genetic research. Nat Genet. 2003;34:119.
Jorde LB, Wooding SP. Genetic variation, classification and ‘race’. Nat Genet. 2004;36:28–33.
Singapore tourism board: About Singapore: Culture, Language and People. http://www.yoursingapore.com/content/traveller/en/browse/aboutsingapore/people-lang-culture.html. Accessed 8 December 2014.
Steinvil A, Fireman E, Wolach O, Rebhun U, Cohen M, Shapira I, et al. The effect of ethnic origin on pulmonary prediction equations in a Jewish immigrant population. Respir Med. 2008;102:919–26.
Collen J, Greenburg D, Holley A, King C, Roop S, Hnatiuk O. Racial discordance in spirometry comparing four commonly used reference equations to the national health and nutrition examination study III. Respir Med. 2010;104:705–11.
Aggarwal AN, Gupta D, Jindal SK. Comparison of Indian reference equations for spirometry Interpretation. Respirology. 2007;12:763–8.
Morris JF, Koski A, Johnson LC. Spirometric standards for healthy nonsmoking adults. Am Rev Respir Dis. 1971;103:57–67.
Choi JK, Paek D, Lee JO. Normal predictive values of spirometry in Korean population. Tuberc Respir Dis. 2005;58:230–42. Korean, translated by Song YJ.
Braun L, Wolfgang M, Dickersin K. Defining race/ethnicity and explaining difference in research studies on lung function. Eur Respir J. 2013;41:1362–70.
The authors declare that they have no competing interests.
NMK conceived this study and participated in its design. NMK and YJS performed the statistical analyses and drafted the article. SYK devised and supervised the entire process. CHS, KHK, JHK, BCS, CKL, and JTL reviewed and revised the article. All authors read and approved the final manuscript.
About this article
Cite this article
Kim, N., Kim, SY., Song, Y. et al. The effect of applying ethnicity-specific spirometric reference equations to Asian migrant workers in Korea. Ann of Occup and Environ Med 27, 14 (2015). https://doi.org/10.1186/s40557-015-0065-0