Comparison of body composition parameters in the study of the association between body composition and pulmonary function

Background The excess adiposity, even in the absence of diseases, is responsible for a decline in pulmonary function, which is considered a predictor of mortality and a risk factor for diseases in several epidemiological studies. However, studies on the association between obesity and pulmonary function have found only few associations or inconclusive results. The aim of the study is to evaluate the association between body composition and spirometric parameters, comparing simple obesity measures such as body mass index (BMI) and waist circumference with more precise body composition measurements such as dual-energy X-ray absorptiometry (DXA) and air-displacement plethysmography (BOD POD). Methods This is an observational, cross-sectional study that used data from the 1978/79 Ribeirão Preto birth cohort (São Paulo, Brazil). The study included 1746 participants from the 5th follow-up of the cohort. Linear regressions were calculated to evaluate the association between BMI, waist circumference, waist–height ratio (WHtR), BOD POD- and DXA-measured fat mass percentage, and spirometric parameters FEV1, and FVC. Results For every 1-kg/m2 BMI increase, FVC decreased by 13 ml in males and by 6 ml in females and FEV1 decreased by 11 ml and 5 ml, respectively. Regarding body composition measurements, for a 1% increase in fat mass assessed by BOD POD, FVC decreased by 16 ml in males and by 8 ml in females and FEV1 decreased by 13 ml and 7 ml, respectively. Hence, negative associations between body measurements and FEV1 and FVC were observed in both genders, especially when using the fat mass measurement and were more expressive in men. Conclusion The anthropometric and body composition parameters were negatively associated with the spirometric variables FVC and FEV1. We have also observed that simple measures such as waist-height ratio were sufficient to detect the association of body composition with pulmonary function reduction. Supplementary Information The online version contains supplementary material available at 10.1186/s12890-021-01543-1.

Excess adiposity is known to interfere with the respiratory system through different mechanisms, including altered respiratory mechanics and pulmonary inflammation. Fat accumulation both in the cavities and in the thoracic and abdominal walls hinders the movement of the respiratory muscles, increases pleural pressure and reduces lung compliance, with a consequent decrease in the compliance of the respiratory system as a whole [7−11]. Furthermore, recent studies have demonstrated the influence of adipose tissue on the secretion of proinflammatory cytokines and on changes in the pulmonary immune system [12−16]. These effects are reflected in altered pulmonary function tests, including a reduction in functional residual capacity, expiratory reserve volume, forced expiratory volume in one second (FEV1), and forced vital capacity (FVC) [17,18].
However, some studies found only a weak association or conflicting results regarding the relationship between obesity and pulmonary function in men and women, mainly because of methodological differences such as the choice of parameter to define obesity [6,[18][19][20][21]. Furthermore, new body composition measurement methods have emerged in recent years, which provided new data for the study of the association between body composition and pulmonary function. These methods range from simple measures such as bioimpedance to more complex tests such as magnetic resonance imaging, computed tomography, dual-energy X-ray absorptiometry (DXA), and airdisplacement plethysmography (BOD POD). Although they have shown a stronger association between visceral or central fat and pulmonary function, such studies are still scarce. Most of them used small samples [22] or samples restricted to healthy volunteers [19,23], making generalization difficult, or they did not control for variables that could interfere with pulmonary function and body composition such as exercise and the presence of comorbidities [24,25].
We hypothesized that the relationship between body composition and pulmonary function reduction may depend on the parameter employed to characterize body composition. The present study aimed to evaluate different body composition parameters and their associations with spirometric values.

Participants
This is a cross-sectional study that used data from a birth cohort started in 1978 in the city of Ribeirão Preto, located 320 km from the capital of São Paulo (SP), Brazil. The estimated population at the time of the study was 682,302 inhabitants, with a Human Development Index of 0.800. Ribeirão Preto is a reference in retail and service activities and is one of the main biomedical research centers in the country [26].
This Brazilian birth cohort study took place between June 1978 and May 1979 in the city of Ribeirão Preto. Its initial objective was to analyze the behavior of some health indicators at birth and their associations with maternal conditions. The second phase of the study focused on the growth of the children (weight and height) during the school years and comprised the period between 1987 and 1989. In a third smaller phase conducted in 1996/97 which only involved men during recruitment for military service, the focus was to study life and health conditions at 18 years of age. The fourth phase of the study was conducted between 2002 and 2004 and was aimed at determining the association between events that had occurred from the prenatal period until the beginning of adult life [27].
In 2016, the data collection for the study "Determinants throughout the life cycle of obesity, precursors of chronic diseases, human capital and mental health" was started, which corresponds to the fifth follow-up of the cohort. The main objective of this phase was to investigate precursors of adult chronic non-communicable diseases, body composition and essential health aspects such as human capital, mental health, and other disorders. The priority was to contact subjects who had participated in the previous phase in 2002 (n = 2103), with 1117 (53.1%) subjects participating in the current data collection. There were also 658 subjects of the original cohort who had not participated in the 2002 follow-up, totaling 1775 subjects.
The present study used data collected during the fifth phase of the cohort. All participants who underwent spirometry (n = 1775) were included. Participants without anthropometric data such as weight and height and pregnant were excluded. Thus, data from 1746 subjects were used in this study (Fig. 1). Participants who did not undergo body composition measurements using BOD POD or DXA, but for whom other measures such as BMI and waist circumference were available, remained in the study because technical limitations of those tests such as the maximum weight supported by the equipment could exclude severely obese individuals (obesity class III).
The project corresponding to the fifth phase of the 1978/79 Ribeirao Preto birth cohort was approved by the Ethics Committee of the Ribeirao Preto Medical School (number 1.282.710). For this study, a new project was submitted to the Ethics Committee for ethical assessment under the terms of Resolution 466/12 of the National Health Council, requesting the waiver of a new informed consent form for the participants and authorization to use the information stored in the database through the term of commitment for data use. This study was approved under number 2.947.100. The names of the participants were not linked to the data to ensure confidentiality.

Independent variables
The following independent variables that were available in the database of the fifth follow-up of the cohort were chosen based on previous studies: waist circumference (cm), body mass index (BMI) calculated by dividing weight (kg) by height in meters squared (m 2 ), waist-toheight ratio (WHtR) (waist circumference divided by height), and fat mass percentage measured by BOD POD and DXA.
The anthropometric measures were obtained by a trained team according to standard techniques, with the subject barefoot and wearing light clothing. Weight was measured with a high-precision electronic scale to the nearest 100 g. Height was measured with a wall-mounted stadiometer to the nearest 0.1 cm. Waist circumference (WC) was measured at the midpoint between the last rib and the upper edge of the iliac crest with an inextensible measuring tape, with the subject standing and abdomen relaxed. The measurement was obtained at the end of a normal expiration.
Body fat percentage was obtained by BOD POD using the COSMED BodPod ® Gold Standard, which calculates the body volume and density by measuring the air displaced in a closed chamber. Weight (kg) was also obtained with a high-precision scale coupled to the equipment. After weighing on the scale, two sequential measurements of body volume were obtained. The subject used tight fitting shorts and a sports bra for women, as well as a swim cap to compress the hair. In addition to weight and body volume, this equipment provides estimates of total fat mass and fat-free mass using predefined equations.
For the assessment of body fat percentage by DXA, a GE Healthcare ® Lunar Prodigy DXA densitometer was used which evaluates body composition using the principle of X-ray attenuation by different body tissues, permitting the estimation of total lean mass, fat mass, and bone tissue. Each subject was subjected to a full body scan, lying completely still in the supine position on a table, with the legs together and arms along the body.

Dependent variables
The outcome variables were FVC and FEV1. Spirometry was performed with a computer-controlled spirometer (Koko Digidoser System, PDS Instrumentation, Louisville, CO, USA), which was calibrated daily. The test was performed with the subject sitting and wearing a nose clip. The technical procedures, acceptability and reproducibility criteria adopted by the American Thoracic Society for spirometry were followed [28].

Adjustment variables
A directed acyclic graph (DAG) was constructed using the DAGitty 2.3 software to verify the minimum adjustment necessary for the analysis model (Fig. 2). The DAG is a causal diagram built on already known theoretical assumptions regarding certain causal relationships [29]. Based on this diagram, the following variables were identified as potential confounders: self-reported skin color (white, black, brown, yellow, or other), educational level in years of schooling (illiterate, 1 to 4, 5 to 8, 9 to 11, and 12 or more years), occupation (first category: industry and construction, second category: less exposed, third category: no occupation), smoking (current, ex-smoker, never smoked), physical activity (classified as low, moderate and high according to the short form of the International Physical Activity Questionnaire-IPAQ [30]), and asthma (yes or no based on the following questions: Have you ever had asthma? and/or Did you experience wheezing in the last year?). Although the DAG identified sex as a confounding variable, we decided to stratify the analysis by sex instead of simply adjusting the model in order to permit comparison with other studies [18,24].

Statistical analysis
The statistical power of the sample was calculated, with a minimum number of 800 participants per group (men and women) and an R 2 of 5% resulting in a power of 99% for this sample. For comparison of the characteristics between genders, the Student t-test was used for continuous variables and the chi-squared test for categorical variables. The association between body composition parameters and pulmonary function was evaluated using multiple linear regression models adjusted for the covariates identified by DAG. Adjustment for height, which is common in studies on obesity and pulmonary function, was not identified as a confounder in the DAG; however, height-adjusted analysis was also performed and the results are reported in supplemental material (Additional file 1: Table S1, S2).
For comparison of the models, the independent variables were standardized based on their means and standard deviations to analyze the associations with Legend. Directed acyclic graph (DAG) illustrating the causal effect of obesity on pulmonary function. Minimum adjustment to estimate the total effect of obesity on pulmonary function: skin color, education level, occupation, smoking, physical activity, sex, and asthma pulmonary function. The estimates, confidence intervals and R 2 were thus used to contrast the models. The Stata v15 program was used for all analyses, adopting a level of significance of 5%. The following WHO classification of BMI categories was considered for description of the sample 1: low weight (BMI < 18.5 kg/m 2 ), eutrophic (BMI ≥ 18.5 and < 25 kg/m 2 ), overweight (BMI ≥ 25 and < 30 kg/m 2 ), and obesity (BMI > 30 kg/ m 2 ).

Results
The characteristics of the sample studied are shown in Table 1. The mean age was 38 years (38.12), with a slight predominance of female participants (52.3%). Most participants had a white skin color (79.1%) and more than 9 years of schooling (87.5%). The proportion of subjects working in industry and construction, who might be at higher risk of occupational exposure to hazardous materials, was low in men (7.1%) and women (5.1%). The prevalence of smokers did not differ between genders (p = 0.737). The presence of asthma was reported by 12.4% of men and 15.1% of women, with no significant difference (p = 0.084).
As can be seen in Table 2, mean height, weight, waist circumference and WHtR were higher in men. The body composition measurements showed a higher fat mass percentage in women. There was a higher proportion of overweight men (44.5%) than women (34.5%) (p < 0.001). The mean spirometric values (in absolute numbers, L or L/s) were higher in men, as expected (p < 0.001) ( Table 2).
Unadjusted linear regression analysis revealed a negative association of BMI, WHtR and BOD POD-and DXA-measured fat mass percentage with FEV1 in men and women. FEV1 was not associated with waist circumference in men, it was in women (Tables 3 and 4). For FVC, all parameters showed a negative association in men and women, except for waist circumference.
After adjusting the analysis for confounders, most body composition variables (BMI, WHtR and BOD POD-and DXA-measured fat mass percentage) were negatively associated with FEV1 and FVC in men and women. Only waist circumference was not associated with either FEV1 or FVC.
Considering the model adjusted for skin color, education level, occupation, smoking, physical activity, and asthma in men and women, and observing the coefficients of association, body fat measured by the BOD POD method was associated with greater negative variation than the other methods, especially in men. The observation of confidence intervals confirmed these findings that are attenuated in women.

Supplemental results
The analysis was also adjusted for height, added to all the adjusting variables described above. Results are very similar, except for the fact that waist circumference is now associated to FEV1 and FVC reduction. This means that waist circumference is a good obesity parameter if controlled by height. Even though, it is not as good as the other variables, as shown by the negative association values that are lower than those relative to fat mass measurements (Additional file 1).

Discussion
This study found a negative association of anthropometric (BMI and WHtR) and body composition parameters (DXA-and BOD POD-measured fat mass) with pulmonary function assessed by FEV1 and FVC in both genders, particularly in men. The adjusted models showed that, for every increase of 1 kg/m 2 in BMI, FVC decreased by 12 ml in men and by 6 ml in women and FEV1 decreased by 10 ml and 5 ml, respectively. Regarding the body composition parameters, for every 1% increase in BOD POD-measured fat mass, FVC decreased by 15 ml in men and by 8 ml in women and FEV1 decreased by 12 ml and 6 ml, respectively. Similar values were obtained for the DXA measurements in which FVC decreased by 11 ml in men and by 7 ml in women for every 1% increase in fat mass and FEV1 decreased by 11 ml and 6 ml, respectively (Additional file 1: Table S3). Rowe et al. [31] evaluated the association of different anthropometric measures (BMI, waist circumference, hip circumference, waist-hip ratio, WHtR, and skinfolds) with pulmonary function assessed by FVC   and FEV1. Although the authors did not report height as a confounding variable in the methods, the linear regression models were adjusted for height, except for the models using BMI and WHtR. The R 2 was used as a measure to compare the associations. The associations with the highest R 2 for both FEV1 and FVC in men and women were obtained when skinfold thickness was used, followed by waist circumference, and the worst associations were obtained using BMI and WHtR. These results suggest that height is an important factor in the study of associations and that it can be addressed in different ways. Some previous studies have shown that waist circumference is a better marker of obesity than BMI in the association with pulmonary function. This can be explained by the fact that waist circumference is a direct marker of central obesity and an indirect marker of visceral fat accumulation; however, those studies adjusted for height in their analyses [2,6,32,33]. In the present study, when height was added to the model (supplemental results), waist circumference was found to be significantly associated with FEV1 in men and women and with FVC in men. Each 1-cm increase in waist circumference decreased FEV1 and FVC by 6 ml in men and by 2 ml in women.
An increasing number of studies have demonstrated the importance of assessing obesity beyond BMI, using more accurate measures of the amount and distribution of body fat, and differentiating fat mass from lean mass. Methods such as bioimpedance, computed tomography, magnetic resonance imaging, BOD POD, and DXA have been employed in an attempt to avoid erroneous classifications of obesity based only on BMI [34]. Today, DXA is considered the gold standard for this assessment [35,36].
In all models, BOD POD-measured body fat was associated with a greater negative variation in both FEV1 and FVC in men and women when compared to the other methods. Similar values were observed for the association of pulmonary function with DXA-measured fat mass, indicating that the relationship of obesity with pulmonary function is based on the content and endocrinological function of fat mass, in addition to the respiratory mechanical injury. Recent studies using fat mass measured by computed tomography found that both visceral adipose tissue and total adipose tissue are associated with poor pulmonary function, irrespective of waist circumference [22,36,37].
Obesity is known to affect pulmonary function regardless of the presence of respiratory, cardiovascular, or metabolic diseases [4,7,8], interfering with respiratory mechanics and lung-thorax compliance. Furthermore, studies have demonstrated that adipose tissue, especially the visceral adipose tissue, is an active tissue in terms of inflammation (cytokine production) and endocrinological activity, releasing hormones that also interfere with pulmonary function [13,15,16].
This study has several strengths. It has used crosssectional data from a birth cohort, which allowed us to obtain a large and representative sample of the population studied, in addition to being subjected to standardization and quality techniques of a longitudinal study in which the responsibility to record detailed quality data is fundamental. Additionally, the study employed sophisticated and accurate methods for body composition measurement. On the other hand, the study has limitations. Its cross-sectional design does not allow to evaluate the causal relationship between body composition and pulmonary function reduction, nor how weight evolution over the time influences the maximal lung function and its decline over the years. The study used a population sample with a limited age range from the city of Ribeirão Preto, whose socioeconomic pattern differs from that of most other Brazilian cities, impairing generalization of our findings to other ages and regions. The original project also did not have this study as the central objective. It is difficult to extract data from a questionnaire that was predefined with other initial general objectives. Some anthropometric parameters such as hip circumference and skinfold thickness that could be useful for comparison with other publications were not measured.
There are other variables that can confuse the association between body composition and pulmonary function, such as diet and socioeconomic status which were not included. Another limitation is that pulmonary function was evaluated only by spirometry; static volumes measurements would certainly improve this study. It is known that obesity can change other functional parameters such as residual volume and gas exchange variables, which spirometry does not measure.
Since reduced pulmonary function and obesity have been consistently shown to be independent predictors of morbidity and mortality and previous studies have demonstrated an association between the two, identification of the method that best demonstrates this association can provide more information about the interaction between adiposity and lung physiology. The present data demonstrate the advantage of more accurate measurements of fat content (BOD POD and DXA) over anthropometric measures. The results mainly indicate the inferiority of the one-dimensional variable waist circumference, which failed to demonstrate the association that was evidenced by WHtR and BMI, probably because waist circumference is a variable that does not contain height. This understanding may be important for future studies, for clinical assessments, and for elaboration of public health policies.

Conclusion
Anthropometric measures of excess weight and increased body composition parameters were associated with lower FVC and FEV1. These associations were more expressive in men. Furthermore, in the study of the association between obesity and pulmonary function, measures that include height for their calculation, such as WHtR, are sufficient and somehow equivalent to more sophisticated body composition parameters obtained by BOD POD and DXA.