Skip to main content
Advertisement
Browse Subject Areas
?

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

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Challenges in Collating Spirometry Reference Data for South-Asian Children: An Observational Study

  • Sooky Lum ,

    s.lum@ucl.ac.uk

    Affiliation Respiratory, Critical Care & Anaesthesia section (Portex Unit), UCL, Institute of Child Health, London, United Kingdom

  • Vassiliki Bountziouka,

    Affiliation Respiratory, Critical Care & Anaesthesia section (Portex Unit), UCL, Institute of Child Health, London, United Kingdom

  • Philip Quanjer,

    Affiliation Department of Pulmonary Diseases and Department of Paediatrics-Pulmonary Diseases, Erasmus Medical Centre, Erasmus University, Rotterdam, Netherlands

  • Samatha Sonnappa,

    Affiliations Respiratory, Critical Care & Anaesthesia section (Portex Unit), UCL, Institute of Child Health, London, United Kingdom, Department of Paediatric Pulmonology, Rainbow Children’s Hospital, Bangalore, India

  • Angela Wade,

    Affiliation Clinical Epidemiology, Nutrition and Biostatistics section, UCL, Institute of Child Health, London, United Kingdom

  • Caroline Beardsmore,

    Affiliation Institute for Lung Health, NIHR Leicester Respiratory Biomedical Research Unit, and Department of Infection, Immunity & Inflammation, University of Leicester, Leicester, United Kingdom

  • Sunil K. Chhabra,

    Affiliation Department of Pulmonary Medicine, Vallabhbhai Patel Chest Institute, University of Delhi, Delhi, India

  • Rajesh K. Chudasama,

    Affiliation Community Medicine Department, PDU Medical College, Rajkot, Gujarat, India

  • Derek G. Cook,

    Affiliation Population Health Research Institute, St George’s, University of London, London, United Kingdom

  • Seeromanie Harding,

    Affiliation Diabetes & Nutritional Sciences Division, Kings College London, London, United Kingdom

  • Claudia E. Kuehni,

    Affiliation Institute of Social and Preventative Medicine, University of Bern, Switzerland

  • K. V. V. Prasad,

    Affiliation Department of Physiology, Vemana Yoga Research Institute, Hyderabad, India

  • Peter H. Whincup,

    Affiliation Population Health Research Institute, St George’s, University of London, London, United Kingdom

  • Simon Lee,

    Affiliation Respiratory, Critical Care & Anaesthesia section (Portex Unit), UCL, Institute of Child Health, London, United Kingdom

  • Janet Stocks

    Affiliation Respiratory, Critical Care & Anaesthesia section (Portex Unit), UCL, Institute of Child Health, London, United Kingdom

Abstract

Availability of sophisticated statistical modelling for developing robust reference equations has improved interpretation of lung function results. In 2012, the Global Lung function Initiative(GLI) published the first global all-age, multi-ethnic reference equations for spirometry but these lacked equations for those originating from the Indian subcontinent (South-Asians). The aims of this study were to assess the extent to which existing GLI-ethnic adjustments might fit South-Asian paediatric spirometry data, assess any similarities and discrepancies between South-Asian datasets and explore the feasibility of deriving a suitable South-Asian GLI-adjustment.

Methods

Spirometry datasets from South-Asian children were collated from four centres in India and five within the UK. Records with transcription errors, missing values for height or spirometry, and implausible values were excluded(n = 110).

Results

Following exclusions, cross-sectional data were available from 8,124 children (56.3% male; 5–17 years). When compared with GLI-predicted values from White Europeans, forced expired volume in 1s (FEV1) and forced vital capacity (FVC) in South-Asian children were on average 15% lower, ranging from 4–19% between centres. By contrast, proportional reductions in FEV1 and FVC within all but two datasets meant that the FEV1/FVC ratio remained independent of ethnicity. The ‘GLI-Other’ equation fitted data from North India reasonably well while ‘GLI-Black’ equations provided a better approximation for South-Asian data than the ‘GLI-White’ equation. However, marked discrepancies in the mean lung function z-scores between centres especially when examined according to socio-economic conditions precluded derivation of a single South-Asian GLI-adjustment.

Conclusion

Until improved and more robust prediction equations can be derived, we recommend the use of ‘GLI-Black’ equations for interpreting most South-Asian data, although ‘GLI-Other’ may be more appropriate for North Indian data. Prospective data collection using standardised protocols to explore potential sources of variation due to socio-economic circumstances, secular changes in growth/predictors of lung function and ethnicities within the South-Asian classification are urgently required.

Introduction

Lung function tests are an integral part of clinical management of respiratory disease but reliable interpretation of results relies on availability of suitable reference data to help distinguish the effects of disease from those of growth and development. Appropriate reference equations are therefore crucial both in clinical management and for interpretation of clinical trials in which lung function is a primary outcome[1,2]. In addition to the major determinants of height, age and sex, lung function is also influenced by ethnicity.[3,4] Some studies have suggested that these differences may be primarily attributed to social deprivation[5,6]. However, while there is evidence that severe deprivation may impact negatively on both growth and lung function,[7,8] recent studies in developed countries have shown that the contribution of socio-economic factors is minimal[913] and that ethnic differences in lung function persist even when such factors are taken into account.

Since publication of the most recent ERS/ATS guidelines for spirometry in 2005 [14] there have only been two publications on spirometry reference ranges for Indian adults.[15,16] By contrast, of the various publications reporting spirometry reference equations for children from the Indian sub-continent (hereafter referred to as South-Asian), seven have been published in the past 15 years. [8,1722] These have, however, been derived using simple regression techniques based on data collected in different parts of the Indian subcontinent, using different equipment in children of different age ranges and socio-economic backgrounds and may therefore not be generalisable.

Development of sophisticated statistical modelling techniques for deriving more robust reference equations has provided opportunities to improve interpretation of lung function results across the age span, with the GLI publishing the first global all-age, multi-ethnic reference equations for spirometry in 2012[3]. These equations were available for 5 distinct ethnic groups, i.e. Caucasian (hereafter referred to as “White”), African-American (hereafter referred to as “Black”), North-East Asian (e.g. North China, Korea), South-East Asian (e.g. South China, Thailand, Taiwan) and Other (consisting of groups other than the 4 main groups and those of mixed ethnic origin). Although some spirometry data South-Asians were provided to the GLI team, results from these studies were sufficiently disparate with respect both to mean results and their distribution to preclude the necessary combination of datasets to derive reliable reference equations [3]. The lack of a ‘normal range’ for South-Asian children currently limits the application and interpretation of their lung function tests.

Two recent studies of 5–12 year old South-Asian school children in London, UK and in Bangalore, India using identical equipment, protocol and quality control(QC)[7,23] showed that after adjusting for height, age and sex[3], average forced expired volume in 1 second (FEV1) and forced vital capacity (FVC) were very similar in Indian urban children residing in Bangalore city to those living in the UK; both being approximately 11% (~0.9z-score) lower than that predicted for White European children[7], with no ethnic differences in the FEV1/FVC ratio. The similarity of results from these two studies, suggested that it should be possible to derive an additional GLI-adjustment (coefficient) suitable for use in South-Asian children. The primary aim of this study was to assess the extent to which existing GLI-ethnic adjustments might fit South-Asian paediatric spirometry data. The secondary aim was to assess any similarities and discrepancies between South-Asian datasets and the feasibility of deriving a suitable GLI-adjustment if needed for interpreting such data. Some results from this study have been reported in abstract form[24].

Materials and Methods

UK centres who had recently (last 15 years) published or collected spirometry data from healthy South-Asian subjects aged 5–18 years or Indian centres who had either submitted such data to the GLI team or recently published their findings were invited to collaborate. See supporting information in S1 File, Section 1.1 for recruitment and exclusion criteria and S1 Table for further details. All centres indicating willingness to participate were requested to provide information on population characteristics (e.g. age, sex, anthropometry, ethnicity, socio-economic circumstances (SEC)) and equipment used. Data were anonymised before submission to this collaboration and came from research studies for which full local ethics approvals were obtained, i.e. from the M.S. Ramaiah Medical College and Teaching Hospitals Ethics Board, Bangalore; Vallabhbhai Patel Chest Institutional Ethics Committee, Delhi; Ethical Committee of the Government Medical College, Surat, Gujarat; Ethical Committee of the National Institute of Nutrition, Hyderabad; Multicentre Research Ethics Committee, UK (for CHASE and DASH study); Leicestershire Research Ethics Committee (for the two Leicester studies) and Research Ethics Committee: London-Hampstead: REC 10/H0720/53 (SLIC study). Parental written consent and verbal assent from each child were obtained prior to assessments.

Data

Among the 11 centres (six UK and five Indian centres) contacted, ten centres responded positively although only nine were able to submit the requested data within the available time frame. Datasets (n = 8413 initially) were available from four centres in India (Bangalore, Delhi, Hyderabad and Gujarat) [7,17,18,21,25] and five in the UK (three in London and two in Leicester; Table 1)[10,11,23,2628]. Records with transcription errors that could not be resolved, with missing values for height, FEV1 or FVC, or where lung function values were deemed implausible (i.e. FEV1/FVC ratio >1.0; FEV1 or FVC ≤0.3L) were discarded (n = 110).

thumbnail
Table 1. Summary of studies included in the collation of South-Asian data.

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

Data management and statistical analysis

Results are presented as mean(SD) for continuous variables or as frequencies(%) for categorical variables. Anthropometry was expressed as sex-specific z-scores based on Indian growth charts derived from well-nourished children[29]. Only three studies were able to provide spirometric flow-volume curves for retrospective quality check on data provided (Table 1). To ascertain overall distribution of data from each study, FEV1, FVC, and FEV1/FVC were initially expressed as z-scores using the GLI-2012 equations for White subjects as a baseline[3].

Within each centre, lung function results were then inspected to ascertain i) the relationship between the lung function z-scores against age and height to ensure that there were no trends in residuals (see S1 File, section 1.2 for details)[3] ii) the spread of data between-subjects iii) whether offsets in FEV1 and FVC relative to the ‘GLI-White’ equation were proportional, as had been observed in over 60 datasets previously submitted to the GLI-2012[30]. After assessing the South-Asian spirometry data in relation to White subjects, the exercise was repeated to assess the fit of existing GLI-ethnic adjustments [3]. The appropriateness of any given reference equation to specific datasets was also ascertained by checking the percentage of healthy subjects within each centre with results falling at or below the 5th centile (i.e. 5% lower limit of normal (LLN) ≤ -1.645 z-scores).

Similarities and discrepancies between datasets, including the potential impact of SEC in Indian-based studies, were examined. Prior to deriving any new South-Asian GLI-adjustment factors (S1 File, Section 1.2.1), datasets with non-proportional reductions in FEV1 and FVC, were excluded and the remaining data were analysed both before and after removing additional datasets that were visibly discrepant from others with respect to either mean offset or between-subject variability.

Analyses were performed using IBM SPSS Statistics version 22 (IBM Corp. Armonk, NY). One-way ANOVA was performed to compare the lung function indices (z-scores based on GLI-White) between centres. Post-hoc pairwise comparisons were subsequently performed using the Tukey’s honestly significant difference test which takes into account sample size when calculating the 95% Confidence Intervals (CI).

Results

Data for analysis were available from 8,124 children (56.3% male; age range 5–17 years), the majority of whom were of Indian origin (Tables 1 and 2).

thumbnail
Table 2. Group characteristics and spirometry results (based on GLI-White equations) according to centre.

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

Anthropometry

Based on the Indian growth charts[29], after adjusting for sex and age, South-Asian children residing in the UK were significantly taller and heavier (by 0.56z-scores and 0.84z-scores respectively) than those living in India (S2 Table). While anthropometry was similar between UK centres, mean values from Indian centres ranged from -0.63 to +0.49z-scores for height and -0.97 to +0.29z-scores for weight (Table 2), children from Northern India (Delhi-centre B) being tallest and heaviest.

Spirometry

Lung function results are summarised in Table 2. When compared with GLI-predicted values from White Europeans, FEV1 in South-Asian children were on average 15% lower, ranging from 4–19% (-0.37 to -1.57 z-scores) according to centre, with an excessive proportion of healthy children (i.e. >>5%) falling below the lower limit of normal (LLN) as defined by the 5th centile. Similar results were found for FVC (average 16% lower; range 6–23%), with Centre B (Delhi) and Centre G (Leicester) having significantly larger FEV1 and FVC than those from all other centres (p<0.01) except Centre I, in whom zFVC was similar to that in Centre G (p = 0.93). However, proportional reductions in FEV1 and FVC within all but two datasets (D and G) meant that FEV1/FVC generally remained independent of ethnicity[3]. The distributions of z-scores according to centre are shown in the Fig 1, illustrating the wide scatter of data from some centres. When data were limited to only those from India (i.e. excluding Pakistani, Bangladeshi, Sri Lankan and Other/mixed South-Asian), results were similar to those observed from all data.

thumbnail
Fig 1. Distribution of lung function data (based on GLI-White) vs. age according to centre.

Symbols: Blue denotes data from boys and Red denotes data from girls. Dashed line = mean value and dotted lines 95% limits of agreement (Mean ± 2 SD) for each dataset. Centres: A = Bangalore; B = Delhi; C = Gujarat; D = Hyderabad; E = CHASE (London); F = DASH (London); G = Leicester City; H = LRC (Leicester); I = SLIC (London). Note the different scales used on the y-axis reflecting the greater spread of FVC than FEV1 data. Note: Age was only recorded to the nearest year for data from Centres B, C and D.

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

Associations between extreme poverty and lung function

Among the UK studies, lung function differences were not explained by socio-economic factors[911]. By contrast, the association between SEC and lung function among children living in India is clearly illustrated in data from Bangalore(A) where both anthropometry and lung volumes were significantly lower in those living in rural areas or exposed to poorer SEC (Table 3)[7]. Anthropometry and lung function were also higher in those from high SEC compared to those from medium or low SEC in the Hyderabad(D) dataset, although the distinction between the latter groups was less clear. Details of SEC were not available from the remaining two India-based studies (B and C).

thumbnail
Table 3. Association between extreme poverty and lung function in children residing in India.

https://doi.org/10.1371/journal.pone.0154336.t003

Extent to which existing GLI-ethnic adjustments might fit South-Asian data

The extent to which the existing GLI-South-East Asian equations (derived from subjects from South China, Thailand and Taiwan, who are geographically closest to those originating from the Indian subcontinent) fitted the South-Asian data is summarised in Table 4. While relative offsets for zFEV1 and zFVC were smaller than when compared with GLI-White equations, the GLI-South-East Asian equations provided a poorer fit for zFEV1/FVC, and an excessive proportion of healthy children had results that fell below the 5th centile for all outcomes. While the entire South-Asian dataset fitted the predicted values for Black-African origin subjects (15% lower than White subjects) reasonably well (Table 5), marked differences in relative mean offsets and distribution of results (SD of z-scores >>1.0) between the various South-Asian datasets were still observed.

thumbnail
Table 4. Spirometry data from South-Asian subjects according to GLI-South-East Asian reference.

https://doi.org/10.1371/journal.pone.0154336.t004

thumbnail
Table 5. Spirometry data from South-Asian subjects according to GLI-Black reference.

https://doi.org/10.1371/journal.pone.0154336.t005

Although not appropriate for most of the centres, application of the ‘GLI-Other’ reference appeared to fit data from Centre B (northern South-Asian) relatively well, with mean (SD) for both FEV1 and FVC approximating 0 (1), albeit with a slight negative offset for FEV1/FVC (Table 6). Under such circumstances, use of an adjusted LLN to reflect the actual 5th centile for FEV1/FVC in that dataset (Table 7; see also S1 File, section 2.1.1) could avoid over-diagnosing abnormalities. Similarly, although the ‘Other’ equation fitted FVC data from centre G, non-proportional differences in FEV1 and FVC meant that there would be considerable under-estimation of airway disease unless an adjusted LLN to reflect the 5th centile was applied to data from this centre (Table 7).

thumbnail
Table 6. Spirometry data from South-Asian subjects according to GLI-Other reference.

https://doi.org/10.1371/journal.pone.0154336.t006

thumbnail
Table 7. Spirometry data from Delhi (Centre B) and Leicester (Centre G) according to GLI-Other reference.

https://doi.org/10.1371/journal.pone.0154336.t007

This exercise was not repeated for the GLI-North-East Asian prediction equations since on average these are only 3% lower than GLI-White, and would therefore under-estimate lung function for all South-Asian data.

Derivation of new GLI-adjustments for South-Asian children

The marked discrepancies in the mean lung function z-scores between centres (Table 2) especially when examined according to SEC (Table 3), meant that not only did none of the existing GLI-equations fit data from all of the individual centres, but precluded derivation of a single South-Asian GLI-adjustment. As data from Centres D and G had non-proportional differences in FEV1 and FVC (Table 2; Tables 46), these data were excluded from subsequent attempts to derive new South-Asian GLI-adjustments. Data with similar mean offsets were collated before modelling the data as described in the S1 File. The selection of datasets for inclusion within different models were based on the i) significantly higher anthropometric (Table 2) and spirometric indices in children from north India compared to other centres (Model 1) ii) potential impact of SEC on lung function from Indian-based studies (Model 2) iii) all remaining datasets with proportional reductions in FEV1 and FVC (Model 3).

After extensive analyses, preliminary GLI-adjustments were derived for the three models (S3 Table). Details of these preliminary models and data fit are presented in the S1 File, section 2.1 to facilitate future investigation and validation by any researchers who are able to collect high quality data from healthy South-Asian children.

Discussion

In this study we collated spirometry data from 8,124 subjects from four Indian-based and five UK-based studies of South-Asian children. After adjustment for age, sex and height using the GLI equations for White Europeans, FEV1 was, on average, 15% lower in South-Asian children, with similar results for FVC emphasising the need for an ethnic adjustment when interpreting spirometry results from such children. By contrast, proportional reductions in FEV1 and FVC in the majority of datasets meant that the FEV1/FVC ratio, which is commonly used in the diagnosis of airway obstruction, remained independent of ethnic group. Despite the smaller lung volumes in relation to height observed across all datasets, there was marked variability between centres in the magnitude of both the mean offset from predicted values in White subjects (which ranged from 4–19%) and the between-subject variability (coefficient of variation). Although the existing ‘GLI-Other’ equation fitted data from North India reasonably well while the ‘Black’ ethnic adjustment provided a better approximation for most of the other South-Asian data than the ‘White equation’, the marked disparities between centres persisted, irrespective of which equation was applied.

These findings confirm previous observations by Aggarwal and colleagues who compared spirometry equations from different subsets of the Indian adult population and found that the three sets of equations (from North, South and West India) did not yield interchangeable results, suggesting that there may be considerable heterogeneity in predicted lung function in various Indian subgroups[31]. While several paediatric reference ranges for specific regions in India have been published, the applicability of these equations to other parts of India has yet to be tested.

Strengths and limitations

Results from this study represent the largest collation of spirometry data from South-Asian children. Data were subjected to sophisticated statistical modelling which has been shown to be highly effective for developing all-age, multi-ethnic reference ranges[3]. However, the current data were collected for various purposes under a wide range of conditions, using different instrumentation and quality control with marked variability in the extent to which potential determinants of lung function had been documented. Thus, despite growing awareness of the potential effects of severe poverty [8,3234], only two of the Indian-based studies reported such details. Similarly, lack of routine documentation regarding birth weight and gestational age in many regions within India, meant that children born prematurely or of low birthweight, both of which can have long term effects on lung function, could not be excluded from such studies. Since this was an opportunistic study based on retrospective data collection, we had no access to information on the extent of migration between rural and urban areas within India, nor details regarding ethnic sub-groups within the South-Asian population from several studies. The population included in this study was predominantly of Indian-origin and the relatively small sample size within other groups such as Pakistani and Bangladeshi precluded any definitive conclusions regarding the extent to which lung function may vary between such groups.

How representative were the subjects recruited to this study?

The issue of which subjects to include in a reference population remains debateable[23,35]. We used an opportunistic approach by contacting authors who had recently submitted data from healthy South-Asian children to the GLI-initiative or who had published such data in the past 15 years, but did not attempt to collate data from South-Asian adults[5], this being beyond the remit of the current study.

All data collated were community based, with most studies specifically designed to recruit from a random sample of schools which catered to a wide range of socio-economic circumstances to ensure a representative sample. Although the Hyderabad data were collected before publication of modern standards for paediatric spirometry[14] which could have introduced some bias, as discussed below, this dataset was excluded when modelling the new GLI-adjustments, so will not have influenced our overall findings.

To assess how representative the data were, we also compared anthropometric data to the current Indian growth charts, derived from 18,666 children attending 10 affluent schools from five major geographical regions of India[29]. After adjusting for age and sex, anthropometric data from all but three centres (A, B and D) were comparable, with mean(SD) values ~0 (1), suggesting that the socio-demographics of most subjects included in this study were representative of well-nourished South-Asian children throughout India and the UK. As previously reported from a study of Punjabi children[36], South-Asian children from Northern India (Centre B) were taller and heavier than those from other centres. Whereas both anthropometry and lung function were lower in children from semi-urban/rural areas when compared with their urban counterparts in centre A, height and weight z-scores remained significantly lower in subjects from Centre D (Hyderabad) despite excluding all data from subjects from middle/low SEC, potential reasons for which could not be ascertained retrospectively.

Variability between and within centres.

One of the most striking findings in this study was the variation in mean lung function between different centres. While FEV1 and FVC were on average 15% lower than that in White children of similar age and height, this difference ranged from 4–19% across the datasets. While mean differences of up to 0.5 z-scores (~5%) in lung function can occur by chance in small samples[37], this could not account for differences between some of the larger centres included in this study. Similarly, there were marked differences in within-centre, between-subject variability. While not necessarily a reason for exclusion per se, inclusion of datasets with markedly increased coefficient of variation when deriving reference equations may decrease their sensitivity to detect lung disease[38]. Review of the literature reveals a similar degree of heterogeneity within the South-Asian population[5,7,39], reasons for which need urgent elucidation if reliable reference ranges are to be developed.

Methodological issues.

An important factor contributing to differences in both mean values between centres and scatter of results within-centres may relate to equipment used. One possible contributory explanation for differences between Centre G and other centres could be systematic overestimation by the spirometer despite daily calibration, since the relative differences in FEV1 and FVC between South-Asian and White children in this study were similar to other reports. While all centres strived to achieve international ATS/ERS standards, the quality of data will inevitably be influenced by the training, expertise and experience of the technical staff. Unfortunately, the lack of stored flow-volume curves with which to assess quality of spirometry data from six of nine centres(Table 1) meant that there were no objective criteria to separate valid from invalid measurements unless data were physiologically impossible. In contrast to some studies wherein strict QC was precluded by the type of spirometer used[5], Centres A (Bangalore) and I (SLIC study) both used a modern ultrasonic device which stores all data. This enabled immediate inspection of flow-volume data at time of collection as well as subsequent over-read, both essential features when assessing children[7,40]. Following quality control, this equipment also allows automatic upload of results into a research database, precluding transcription errors that may occur during manual entry.

In the current study, mean zFEV1/FVC from all but two centres approximated zero (Table 2). When studying healthy children, a relatively low FVC combined with a normal or high FEV1/FVC ratio, especially if combined with a wide scatter of results, may reflect poor subject cooperation, rather than lung restriction[41,42]. Additional discrepancies will have been introduced in several studies by failure to record age to at least 1 decimal year (Fig 1), the importance of which for accurate prediction of lung function was emphasised recently[3].

Impact of socio-economic circumstances.

The marked secular changes in both lung function and growth that have occurred within developing nations such as India over the past 50 years[43] may contribute to increased variability of results in South-Asians, especially if exposed to extreme social deprivation[7], but are less likely to explain discrepancies in studies undertaken in the UK[9,26]. Despite minimal impact of SEC on lung function in children living in the UK[9,11,44], we recently reported a potential threshold effect of poverty, with lung function being significantly lower in children residing in rural India than in their urban counterparts[7], which is in keeping with previous evidence[8]. In an attempt to avoid either over-diagnosis of lung disease in otherwise healthy children with evidence of growth restriction or poorer living conditions, or under-diagnosis in well-nourished children, separate GLI-adjustment factors were derived for these groups.

Regional differences? The extent to which inclusion of various South-Asian subgroups contributed to both heterogeneity of lung function within studies and differences between centres is difficult to ascertain due to sample size and inconsistency in the way in which such details were recorded. Despite inevitable mixing, people from north India have different ancestry from those in the south[45,46] which impacts on both anthropometry and lung function, and which may account for the higher lung function for data collected in Delhi[47]. Consequently, a separate preliminary equation was derived for this dataset.

While it is customary to use stature as a proxy for lung volume, this approach does not adjust for frame size including leg length in relation to torso (i.e. sit/standing height ratio). We have not reported detailed anthropometry regarding sit/stand ratio or chest dimensions in this report, as so few studies have included such measures, and those that did found that they made a relatively small contribution to observed differences between White and South-Asian children[9,26]. Nevertheless, while absolute lung volumes adjusted for stature may differ between ethnic groups, there is strong evidence that among healthy individuals, such differences in FEV1 and FVC are proportional, as observed in over 60 datasets submitted to the GLI-2012[30], resulting in an FEV1/FVC ratio that is relatively independent of ethnicity and body size[3,23]. Datasets with non-proportional reductions in FEV1 and FVC based on the ‘White equation’ were excluded prior to deriving South-Asian GLI-adjustment factors since not only may an elevated FEV1/FVC in health reflect difficulties in achieving a full exhalation, but such findings are contrary to the majority of data previously published or submitted to the GLI.

Clinical and practical implications

While this study has provided clear evidence of ethnic differences in spirometric lung function among South- Asian children, it must be emphasised that the proportionally smaller values of FEV1 and FVC observed do not reflect intrinsic differences in the functionality of the respiratory system, as reflected by the consistency of the FEV1/FVC ratio. While derivation of an ethnic GLI-adjustment would potentially be the most efficient and cost-effective way of accommodating those of South-Asian origin as has been done for those of Black-African origin, results from this study suggest that no single adjustment will currently fit all data collected from this group, since there appear to be complex factors contributing to the heterogeneity observed within the South-Asian population. While it may be possible to address some of these factors in future, for example by greater standardisation of equipment, methodology and quality control, whereas the contribution of others, such as extreme poverty, may diminish over time with improved living standards, there may be some intrinsic differences related to genetic ancestry, which need further investigation. Meanwhile, the clinical dilemma of which equations to choose when interpreting spirometry results from South-Asian children remains. Personal communication with the authors from the Indian based studies included in this initiative, revealed that within most clinical centres in India, technicians frequently rely on the unspecified and disparate reference values presented as default by the manufacturers, a practice not unknown even in the UK[1]. Under these circumstances, adoption of a common standard such as the GLI all-age equations which are already available in most commercial lung function equipment (http://www.ers-education.org/guidelines/global-lung-function-initiative.aspx) would be far preferable, even if this did entail a small offset. Results from the current study suggest that until an improved solution can be found, use of the GLI-Black equations for the majority of South-Asian children, or the GLI-other equation for those originating from North India, could provide an interim solution. For any clinical centre with access to a sizeable population (ideally at least 300[37]) of healthy children, use of such equations could be validated and if necessary the lower limit of normal adjusted to the local population, as discussed above (Table 7; see also S1 File, section 2.1.1, S6 Table).

Given that such an approach will need updating in future, we have included three alternative equations in the S1 File, in the hopes that this might stimulate further research using high quality equipment and standardised protocols to investigate the potential impact that SEC, geographical location or ethnic sub-groups may impact anthropometry and lung function within the South-Asian. When attempting to derive any new adjustment for South Asians, there needs to be a balance between limited generalisability if equations are based on a few homogenous datasets collected under ideal circumstances (which may not be achievable during routine clinical assessments) versus the lack of discrimination that will occur if such equations are derived from extremely heterogeneous data.

Conclusions and Future Directions

Variations in both the mean level of lung function and the degree of within-subject variability meant that none of the existing GLI-ethnic adjustments provided an ideal fit for all the datasets examined in this study and also precluded derivation of a single new South-Asian paediatric GLI adjustment. Given the current wide variation in equations used to interpret lung function from South-Asian children, an alternative approach would be to express results from South-Asian children using either the ‘GLI-Black’ or ‘GLI -Other’ equations, depending on region of origin, both of which are widely available in current spirometers.

Further prospective data collection among South-Asians residing within the Indian sub-continent as well as those who have migrated, using standardised protocols over a wider age range, is required to explore the potential sources of variation observed, with regular updating if there is evidence of significant secular changes. When undertaking such studies, age should be recorded in decimal years and standardised anthropometry include assessments of sitting height. All lung function assessments need to be undertaken according to ATS/ERS guidelines[14], adapted for children where necessary[40] using equipment that allows prospective QC at time of data collection, storage of all data for subsequent independent over-read and automated export of results to avoid transcription errors. Furthermore, documentation of relevant past and current medical history, environmental exposures and SEC should be as standardised as much as possible to facilitate subsequent collation or sub-division of data, and clarification of potential sources of variability (S8 Table).

In conclusion, despite some progress being made, the challenges we faced in attempting to derive spirometry reference equations for South-Asian children are perhaps best summarised by the following quote from Earl C. Kelley, Professor of Secondary Education at Wayne University (1951)[48]:

"We have not succeeded in answering all our problems—indeed we sometimes feel we have not completely answered any of them. The answers we have found have only served to raise a whole set of new questions. In some ways we feel that we are as confused as ever, but we think we are confused on a higher level and about more important things. So this report does not purport to give final answers, or to claim that we now "know how to do it". We see more need for revision than ever. But we are doing better than we did. And this is a progress report, rendered with humility because of the unsolved problems we see now which we could not see before."

Supporting Information

S1 Fig. Amendment to Excel Sheet calculator for calculation of lung function z-scores based on preliminary GLI-adjustments (for Model 3b).

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

(PDF)

S2 Fig. Data fit of the GLI-adjustment for Centre B using the smoothing function, plotted against height.

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

(PDF)

S3 Fig. FVC z-scores based on Model 2 (GLI-adjustments for A2-3 & C) according to centre.

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

(PDF)

S4 Fig. Distribution of lung function z-scores calculated using GLI-adjustment based on Model 3a.

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

(PDF)

S5 Fig. FVC z-scores calculated using Model 3a (GLI-adjustments for A1, E, F, H & I) according to centre.

https://doi.org/10.1371/journal.pone.0154336.s005

(PDF)

S6 Fig. FVC z-scores based on Model 3b (GLI-adjustments for A1, H & I) according to centre.

https://doi.org/10.1371/journal.pone.0154336.s006

(PDF)

S1 Table. Recruitment and exclusion criteria according to respective studies.

https://doi.org/10.1371/journal.pone.0154336.s008

(PDF)

S2 Table. Comparison of anthropometry between children residing in the UK and in India.

https://doi.org/10.1371/journal.pone.0154336.s009

(PDF)

S3 Table. Preliminary GLI-adjustments according to the various models.

https://doi.org/10.1371/journal.pone.0154336.s010

(PDF)

S4 Table. Lung function results based on Model 1 GLI-coefficients derived from Centre B (Delhi).

https://doi.org/10.1371/journal.pone.0154336.s011

(PDF)

S5 Table. Lung function results based on Model 2 GLI-coefficients derived from Centres A2-3 & C.

https://doi.org/10.1371/journal.pone.0154336.s012

(PDF)

S6 Table. Lung function results based on GLI-coefficients derived from Centres A1(urban), E, F, H & I (Model 3a).

https://doi.org/10.1371/journal.pone.0154336.s013

(PDF)

S7 Table. Lung function results based on Model 3b GLI-coefficients derived from Centres A1, H and I

https://doi.org/10.1371/journal.pone.0154336.s014

(PDF)

S8 Table. Data required for prospective data collection.

https://doi.org/10.1371/journal.pone.0154336.s015

(PDF)

Acknowledgments

The SLIC study team acknowledges the support of the National Institute for Health Research (NIHR), through the Comprehensive Clinical Research Network and the NIHR Biomedical Research Centre at Great Ormond Street Hospital for Children NHS Foundation Trust and University College London.

Author Contributions

Conceived and designed the experiments: S. Lum JS. Performed the experiments: S. Lum SS CB SKC RKC DC SH CEK KVVP PW. Analyzed the data: VB PQ AW. Contributed reagents/materials/analysis tools: S. Lum SS CB SKC RKC DC SH CEK KVVP PW. Wrote the paper: S. Lum PQ JS AW VB SS CB SKC RKC DC SH CEK KVVP PW S. Lee. Collated the data: S. Lum S. Lee. Interpreted the data: JS PQ S. Lum.

References

  1. 1. Stanojevic S, Wade A, Stocks J (2010) Reference values for lung function: past, present and future. Eur Respir J 36: 12–19. pmid:20595163
  2. 2. Subbarao P, Lebecque P, Corey M, Coates AL (2004) Comparison of spirometric reference values. Pediatr Pulmonol 37: 515–522. pmid:15114552
  3. 3. Quanjer PH, Stanojevic S, Cole TJ, Baur X, Hall GL, Culver BH, et al. (2012) Multi-ethnic reference values for spirometry for the 3-95-yr age range: the global lung function 2012 equations. Eur Respir J 40: 1324–1343. pmid:22743675
  4. 4. Pellegrino R, Viegi G, Brusasco V, Crapo RO, Burgos F, Casaburi R, et al. (2005) Interpretative strategies for lung function tests. Eur Respir J 26: 948–968. pmid:16264058
  5. 5. Duong M, Islam S, Rangarajan S, Teo K, O'Byrne PM, Schunemann HJ, et al. (2013) Global differences in lung function by region (PURE): an international, community-based prospective study. Lancet Respir Med 1: 599–609. pmid:24461663
  6. 6. Burney PG, Hooper RL (2012) The use of ethnically specific norms for ventilatory function in African-American and white populations. Int J Epidemiol 41: 782–790. pmid:22434864
  7. 7. Sonnappa S, Lum S, Kirkby J, Bonner R, Wade A, Subramanya V, et al. (2015) Disparities in pulmonary function in healthy children across the Indian urban-rural continuum. Am J Respir Crit Care Med 191: 79–86. pmid:25412016
  8. 8. Raju PS, Prasad KVV, Ramana YV, Balakrishna N, Murthy KJ (2005) Influence of socioeconomic status on lung function and prediction equations in Indian children. Pediatr Pulmonol 39: 528–536. pmid:15789442
  9. 9. Lum S, Bountziouka V, Sonnappa S, Wade A, Cole TJ, Harding S, et al. (2015) Lung function in children in relation to ethnicity, physique and socioeconomic factors. Eur Respir J 46: 1662–1671. pmid:26493801
  10. 10. Strippoli MP, Kuehni CE, Dogaru CM, Spycher BD, McNally T, Silverman M, et al. (2013) Etiology of ethnic differences in childhood spirometry. Pediatrics 131: e1842–1849. pmid:23713103
  11. 11. Whitrow MJ, Harding S (2008) Ethnic differences in adolescent lung function: anthropometric, socioeconomic, and psychosocial factors. Am J Respir Crit Care Med 177: 1262–1267. pmid:18323540
  12. 12. Menezes AM, Wehrmeister FC, Hartwig FP, Perez-Padilla R, Gigante DP, Barros FC, et al. (2015) African ancestry, lung function and the effect of genetics. Eur Respir J 45: 1582–1589. pmid:25700383
  13. 13. Harik-Khan RI, Fleg JL, Muller DC, Wise RA (2001) The effect of anthropometric and socioeconomic factors on the racial difference in lung function. Am J Respir Crit Care Med 164: 1647–1654. pmid:11719304
  14. 14. Miller MR, Hankinson J, Brusasco V, Burgos F, Casaburi R, Coates A, et al. (2005) Standardisation of spirometry. Eur Respir J 26: 319–338. pmid:16055882
  15. 15. Chhabra SK, Kumar R, Gupta U, Rahman M, Dash DJ (2014) Prediction equations for spirometry in adults from northern India. Indian J Chest Dis Allied Sci 56: 221–229. pmid:25962195
  16. 16. Memon MA, Sandila MP, Ahmed ST (2007) Spirometric reference values in healthy, non-smoking, urban Pakistani population. J Pak Med Assoc 57: 193–195. pmid:17489528
  17. 17. Chhabra SK, Vijayan VK, Rahman M, Mittal V, Singh PD (2012) Regression equations for spirometry in children aged 6 to 17 years in Delhi region. Indian J Chest Dis Allied Sci 54: 59–63. pmid:22779126
  18. 18. Doctor TH, Trivedi SS, Chudasama RK (2010) Pulmonary function test in healthy school children of 8 to 14 years age in south Gujarat region, India. Lung India 27: 145–148. pmid:20931033
  19. 19. Budhiraja S, Singh D, Pooni PA, Dhooria GS (2010) Pulmonary functions in normal school children in the age group of 6–15 years in north India. Iran J Pediatr 20: 82–90. pmid:23056687
  20. 20. Prasad R, Verma SK, Agrawal GG, Mathur N (2006) Prediction model for peak expiratory flow in North Indian population. Indian J Chest Dis Allied Sci 48: 103–106. pmid:16696523
  21. 21. Raju PS, Prasad KVV, Ramana YV, Murthy KJ (2004) Pulmonary function tests in Indian girls—prediction equations. Indian J Pediatr 71: 893–897. pmid:15531830
  22. 22. Vijayan VK, Reetha AM, Kuppurao KV, Venkatesan P, Thilakavathy S (2000) Pulmonary function in normal south Indian children aged 7 to 19 years. Indian J Chest Dis Allied Sci 42: 147–156. pmid:11089318
  23. 23. Lum S, Bountziouka V, Sonnappa S, Cole TJ, Bonner R, Stocks J (2015) How "healthy" should children be when selecting reference samples for spirometry? Eur Respir J 45: 1576–1581. pmid:25700391
  24. 24. Kirkby J, Lum S, Stocks J, Bonner R, Sonnappa S (2014) Adaptation of the GLI-2012 spirometry reference equations for use in Indian children. European Respiratory Journal 44: 191.
  25. 25. Raju PS, Prasad KVV, Ramana YV, Ahmed SK, Murthy KJ (2003) Study on lung function tests and prediction equations in Indian male children. Indian Pediatr 40: 705–711. pmid:12951372
  26. 26. Whittaker AL, Sutton AJ, Beardsmore CS (2005) Are ethnic differences in lung function explained by chest size? Arch Dis Child Fetal Neonatal Ed 90: F423–428. pmid:15871993
  27. 27. Whincup PH, Nightingale CM, Owen CG, Rudnicka AR, Gibb I, McKay CM, et al. (2010) Early emergence of ethnic differences in type 2 diabetes precursors in the UK: the Child Heart and Health Study in England (CHASE Study). PLoS Med 7: e1000263. pmid:20421924
  28. 28. Barone-Adesi F, Dent JE, Dajnak D, Beevers S, Anderson HR, Kelly FJ, et al. (2015) Long-Term Exposure to Primary Traffic Pollutants and Lung Function in Children: Cross-Sectional Study and Meta-Analysis. PLoS One 10: e0142565. pmid:26619227
  29. 29. Khadilkar VV, Khadilkar AV, Cole TJ, Sayyad MG (2009) Cross-sectional growth curves for height, weight and body mass index for affluent Indian children, 2007. Indian Pediatr 46: 477–489. pmid:19556658
  30. 30. Quanjer PH, Stanojevic S, Stocks J, Hall GL, Prasad KV, Cole TJ, et al. (2010) Changes in the FEV(1)/FVC ratio during childhood and adolescence: an intercontinental study. Eur Respir J 36: 1391–1399. pmid:20351026
  31. 31. Aggarwal AN, Gupta D, Jindal SK (2007) Comparison of Indian reference equations for spirometry interpretation. Respirology 12: 763–768. pmid:17875069
  32. 32. Lawlor DA, Ebrahim S, Davey Smith G (2004) Association between self-reported childhood socioeconomic position and adult lung function: findings from the British Women's Heart and Health Study. Thorax 59: 199–203. pmid:14985552
  33. 33. Jackson B, Kubzansky LD, Cohen S, Weiss S, Wright RJ (2004) A matter of life and breath: childhood socioeconomic status is related to young adult pulmonary function in the CARDIA study. Int J Epidemiol 33: 271–278. pmid:15082626
  34. 34. Stein CE, Kumaran K, Fall CH, Shaheen SO, Osmond C, Barker DJ (1997) Relation of fetal growth to adult lung function in south India. Thorax 52: 895–899. pmid:9404378
  35. 35. Stanojevic S, Wade A, Lum S, Stocks J (2007) Reference equations for pulmonary function tests in preschool children: a review. Pediatr Pulmonol 42: 962–972. pmid:17726704
  36. 36. Singh SP, Sidhu LS, Malhotra P (1987) Growth performance of Punjabi children aged 6–12 years. Ann Hum Biol 14: 169–179. pmid:3606028
  37. 37. Quanjer PH, Stocks J, Cole TJ, Hall GL, Stanojevic S (2011) Influence of secular trends and sample size on reference equations for lung function tests. Eur Respir J 37: 658–664. pmid:20817707
  38. 38. Stanojevic S, Wade A, Stocks J, Hankinson J, Coates AL, Pan H, et al. (2008) Reference ranges for spirometry across all ages: a new approach. Am J Respir Crit Care Med 177: 253–260. pmid:18006882
  39. 39. Khongsdier R (2001) Body mass index of adult males in 12 populations of northeast India. Ann Hum Biol 28: 374–383. pmid:11459235
  40. 40. Kirkby J, Welsh L, Lum S, Fawke J, Rowell V, Thomas S, et al. (2008) The EPICure study: comparison of pediatric spirometry in community and laboratory settings. Pediatr Pulmonol 43: 1233–1241. pmid:19009621
  41. 41. Swanney MP, Beckert LE, Frampton CM, Wallace LA, Jensen RL, Crapo RO (2004) Validity of the American Thoracic Society and other spirometric algorithms using FVC and forced expiratory volume at 6 s for predicting a reduced total lung capacity. Chest 126: 1861–1866. pmid:15596685
  42. 42. Aaron SD, Dales RE, Cardinal P (1999) How accurate is spirometry at predicting restrictive pulmonary impairment? Chest 115: 869–873. pmid:10084506
  43. 43. Quanjer PH, Kubota M, Kobayashi H, Omori H, Tatsumi K, Kanazawa M, et al. (2015) Secular changes in relative leg length confound height-based spirometric reference values. Chest 147: 792–797. pmid:25254426
  44. 44. Kuehni CE, Strippoli MP, Spycher BD, Silverman M, Beardsmore CS (2015) Lung function in the children of immigrant and UK-born south-Asian mothers. Eur Respir J 45: 1163–1166. pmid:25573409
  45. 45. Ali M, Liu X, Pillai EN, Chen P, Khor CC, Ong RT, et al. (2014) Characterizing the genetic differences between two distinct migrant groups from Indo-European and Dravidian speaking populations in India. BMC Genet 15: 86. pmid:25053360
  46. 46. Majumder PP (1998) People of India: Biological diversity and affinities. Evolutionary Anthropology 6: 100–110.
  47. 47. Chhabra SK (2009) Regional variations in vital capacity in adult males in India: comparison of regression equations from four regions and impact on interpretation of spirometric data. Indian J Chest Dis Allied Sci 51: 7–13. pmid:19317357
  48. 48. Kelley EC (1951) Workshop way of learning: Harper & Row.