Abstract
African-Americans have smaller lung function compared with European-Americans. The aim of this study was to disentangle the contribution of genetics from other variables on lung function.
A cohort was followed from birth to 30 years of age in Brazil. Several variables were collected: genomic analysis based on DNA; forced expiratory volume in 1 s (FEV1) and forced vital capacity (FVC) obtained by spirometry; height measured by anthropometrists; and thorax circumference evaluated by photonic scanner. Crude and adjusted linear regression models were calculated according to African ancestry.
The sample comprised 2869 participants out of 3701 members of the cohort. Males with higher African ancestry by DNA analysis had a smaller FEV1 (−0.13 L, 95% CI −0.23– −0.03 L) and FVC (−0.21 L, 95% CI −0.32– −0.09 L) compared with those with less African ancestry, having accounted for height, sitting to standing height ratio and other confounders. Similar effects were seen in females.
After adjustment, ancestry remained significantly associated with lung function, but the large effect of adjustment for confounding among males (but not females) does not allow us to exclude the possibility that residual confounding may still account for these findings.
Abstract
Smaller lung function is related to African ancestry but adjustments for confounders reduced the effect sizes http://ow.ly/Iq0Ve
Introduction
African-Americans have shown consistently smaller lung function compared with European-Americans [1–9]. In addition, forced expiratory volume in 1 s (FEV1) and forced vital capacity (FVC) in African-Americans differed from those in Caucasians by the same percentage, signifying that, for the same age and height, lung dimensions differed proportionately [6]; it should be mentioned that the information about race and ethnicity in most of the literature has been self-reported. However, according to Kumar et al. [3], in the absence of genetic typing, predicted values in self-reported African-Americans may be biased by up to 200 mL. It should be emphasised that most of the literature does not define race at all [9].
Commonly reported ethnic and racial differences in pulmonary function have been attributed to differences in body build, such as chest size relative to height (approached by the Cormic index: sitting to standing height ratio). Lung function in general is proportional to body size as indicated by standing height. However, ethnic groups with smaller lung function also have historically lower socioeconomic status (determining bodily development in early life and leading to trends in body size and pulmonary function that continue into adulthood) and possibly other environmental factors [10–19], which are important determinants of lung growth and development.
A recent systematic review about race and/or ethnicity and lung function showed that researchers have failed to explore this association in a proper way, and an appeal was made to the scientific community to clarify this relationship [9]. Of the 226 eligible papers published between 1922 and 2008, in 83.6% of the articles it was reported that “other racial and ethnic groups” had a smaller lung volume compared with “white” groups, and 94% failed to examine socioeconomic status; in those papers that reported smaller lung function in “other racial and ethnic groups”, 21.8 and 29.4% of explanations cited inherent factors and anthropometric differences, respectively, whereas 23.1% cited environmental and social factors [9].
The debate on this association has been ongoing for decades and in recent years it has been emphasised due to different interpretations of the available data [5, 7]. The present study addresses whether lung function varies according to the proportion of African genomic ancestry (based on DNA) and how this interacts with anthropometric measurements, socioeconomic status, smoking, comorbidities and birth weight.
Methods
In 1982, the maternity hospitals in Pelotas, a southern Brazilian city, were visited daily and all births were identified. Those live-borns whose family lived in the urban area of the city were examined and their mothers interviewed (n=5914). These subjects have been followed up several times [20–22].
The study protocol was approved by the Federal University of Pelotas Ethics Committee. In 2012–2013 (at 30 years of age), we invited all the cohort members to participate in a follow-up examination at the study clinic. Participants gave informed consent. Spirometry was undertaken in individuals who did not present any exclusion criteria (recent abdominal or thoracic surgery, pregnancy or possible pregnancy, cardiovascular diseases, retinal detachment or eye surgery, tuberculosis under treatment, or comprehension difficulties) using a portable, battery-operated, ultrasound transit time-based Easy-One spirometer (ndd Medical Technologies, Zurich, Switzerland). Spirometry was carried out following American Thoracic Society and European Respiratory Society recommended procedures [23], before and 15 min after inhalation of 200 µg salbutamol. The spirometric measures used in the current analysis were pre-bronchodilator FEV1 and FVC.
Skin colour was self-reported by the subjects and categorised into white, brown (or mixed), black or as “other” (52 subjects reported yellow skin colour and 45 reported as indigenous).
Standing height was assessed to the nearest 0.1 cm using a full-length wall-mounted stadiometer (SECA 240; Birmingham, UK) [24]; all anthropometric measurements were carried out twice and the mean value was used. Thorax circumference was measured by photonic scanner (TC2 3D Body Scanner; Textile Clothing and Technology Corporation, Cary, NC, USA).
The genomic ancestry analysis was based on DNA samples of 3736 cohort members who had been evaluated at 22–23 years of age. The Illumina Omni 2.5M array (Illumina, Inc., San Diego, CA, USA) was used for genotyping, and the Admixture analyses [25] were based on 370 539 single nucleotide polymorphisms (SNPs) shared by samples from the HapMap Project (http://hapmap.ncbi.nlm.nih.gov), the Human Genome Diversity Project (HGDP) (www.hagsc.org/hgdp) and the Pelotas cohort. As external panels, we used the following HapMap samples: 266 Africans, 262 Europeans (Americans from European ancestry and Italian), 77 admixed Mexican-Americans and 83 African-Americans; and 93 native Americans from the HGDP. For each subject the proportion of European, African-American and native American ancestry was estimated. By using parental populations from different public sources, we believe we included the relevant contributors to Latin-American populations: Yoruba from Western Africa (one of the main ethnic groups that contributed to the African diaspora from this region), Luyha from Kenya (East Africa), two European populations and four native American populations.
The following variables were considered as possible confounders in the multivariate analyses: standing height and sitting/standing height ratio, chest circumference, family income and maternal schooling at birth (as number of completed years), maternal smoking during pregnancy (number of cigarettes smoked per day), birth weight in grams (using paediatric scales with a precision of 10 g), household asset index in childhood (obtained through factor analysis and based on the ownership of household goods), pneumonia-related hospitalisation in the first 2 years and wheezing (defined according to the International Study of Asthma and Allergies (ISAAC) questionnaire [26]). Socioeconomic status was evaluated again during the 2012–2013 follow-up visit on the basis of number of years of achieved schooling and the possession of 17 different household assets [27]. A household asset index based on the possession of household assets was estimated using factor analysis. Prevalence of smoking and wheezing (in the past 12 months) were also collected during the 2012–2013 follow-up visit.
Statistical analyses
Descriptive analysis was performed according to the African ancestry percentiles in our cohort (0–40, 40–80, and 80–100%), which corresponded to the following categories for proportion of African ancestry: 0–5.2, 5.3–25.8, and 25.9–87.9. Analyses were based on a conceptual framework of the determinants of the association between African ancestry and pulmonary function. Crude and adjusted linear regression models for FEV1 and FVC according to self-reported skin colour and African ancestry were carried out.
The regression models were adjusted by adding new variables to African ancestry (model 1): standing height (model 2); sitting/standing height ratio and chest circumference (model 3); current asset index, current achieved schooling, current smoking, wheezing (in 2012), family income and maternal schooling at birth, birth weight, gestational age, asset index in 1984, maternal smoking during pregnancy, and history of wheezing and pneumonia from birth to 1984 (model 4).
Analyses were stratified by sex and statistical significance was a priori defined as p<0.05. All analyses were performed using R 3.0.2 (The R Project for Statistical Computing; www.r-project.org).
Results
In the 2012–2013 follow-up (participants aged 30 years old), 3701 subjects were interviewed, which added to the 366 subjects who were known to have died represents a 68.8% follow-up rate. The present analysis took into account only those individuals who had information on ancestry (from follow-up at 22–23 years of age) and on lung function (obtained at the 2012–2013 visit) (n=2869).
Table 1 shows the percentage of subjects with spirometric and ancestry data and the total number of individuals evaluated at 2012–2013 follow-up visit. More than 75% of individuals evaluated in the last follow-up visit were included in the present analysis. For those with the highest family income at birth and those born to mothers who had ≥12 years of schooling the percentage participating in the last follow-up visit was lower compared with those individuals with less favourable socioeconomic variables.
Proportion of eligible individuals among the entire 2012 follow-up# according to study variables
Table 2 shows that those who had more than one quarter African ancestry had lower birth weight, as well as income and maternal schooling at birth. By the same token, socioeconomic indicators at 30 years were unfavourable for those with higher proportion of African ancestry. Current anthropometric variables were inversely associated with African ancestry, except mean chest circumference (table 2).
Sample descriptions stratified according to proportion of African ancestry
Table 3 shows the crude and adjusted association between African ancestry with FEV1 and FVC stratified by sex. We opted for using the denominator of 927 for males and 943 for females due to missing information for some covariables. For males, those with more than 25.9% African ancestry had a smaller FEV1 (crude analysis β= −0.31 L, 95% CI −0.43 – −0.20 L; and full adjusted β (model 4)= −0.13 L; 95% CI −0.23 – −0.03 L) and FVC (crude analysis β= −0.46 L; 95% CI −0.59 – −0.32 L; and fully adjusted β (model 4)= −0.21 L; 95% CI −0.32 – −0.09 L) compared with those with <5% African ancestry. For females, the results were similar, but the strength of the association was lower compared with males.
Linear regression coefficients (in litres) for FEV1 and FVC according to African ancestry
The linear regression coefficients (height and sitting/standing height ratio in cm) and the adjusted-R2 values according to African ancestry are shown in table 4. African ancestry explains less of the variability for height than for the sitting/standing height ratio, especially in males; the R2 for the crude model was 2% for height and around 6% for the sitting/standing height ratio.
Linear regression coefficients and adjusted-R2 values for height and sitting/standing height ratio according to African ancestry
Discussion
Our results show smaller lung function volumes among those subjects with increased African ancestry after adjustment for anthropometrics and several other variables. We have shown that the sitting to standing height ratio is related to African ancestry having taken into account socioeconomic confounders, whereas overall height is not affected by African ancestry per se (table 4). This confirms the work of others showing that height, sitting/standing height ratio and thoracic circumference are linked to the association between African ancestry and lung function [2, 6, 7, 13, 28, 29]. We did not include sitting height on its own in our analysis due to its colinearity with other variables.
Having taken into account height and the sitting to standing height ratio African ancestry was still associated with significantly lower FEV1 (−150 mL in men and −200 mL in women) and FVC (−240 and −190 mL, respectively). After then accounting for all the additional confounders in table 2 these estimates were only changed by 10–30 mL. We accept that residual confounding may still be an issue and that 18% of the variance in height in our subjects was related to our measured confounders. However, the effect of African ancestry, height and sitting to standing height ratio on determining FEV1 and FVC are much larger than that of the confounders we have measured.
Socioeconomic status may affect health in a variety of ways [30], including adverse environmental exposures, and it is unlikely that specific indicators, or even a combination of indicators, can take into account all of its impact. In addition, socioeconomic status has impacts throughout life, even before birth, and current status may differ from that present at key stages of growth and development. In fact, separating genetic African ancestry from socioeconomic status may be unfeasible for some time as both have been deeply linked for centuries and it will require a considerable population of individuals with high African ancestry living their whole life in adequate socioeconomic status. A key strength of our study was the availability of socioeconomic status at more than one time point.
Although several researchers have carried out genomic ancestry analysis it should be mentioned that this method has some limitations; in the present study inferred ancestry was performed through a larger set of SNPs from different datasets and not only from Pelotas. We believe we have adjusted for a variety of indicators of ancestral history that contribute to current Latin-American populations, and therefore have partially reduced the problem of studying an unusually varied population.
The additional analysis carried out to identify whether height was more determined by genetics or by the environment shows that the inclusion of African ancestry resulted in minor increases in the R2 and adjusted-R2 values for both sexes. There was no convincing evidence supporting that the genetic differences among ethnic groups in our sample played an important role in explaining the differences in height.
The complex relationship between lung function and race and/or ethnicity [5, 7] has been addressed in several studies; a key issue in this discussion is to disentangle the contribution of “genetics” and “other variables” which is very difficult or impossible to achieve at present time. To the best of our knowledge this is the first study in a middle-income country assessing the independent role of genomic ancestry on lung function and several anthropometric measures, indicators of socioeconomic status at birth and later in life, morbidities, smoking, and birth weight in a longitudinal cohort.
Brazil is a huge country and it is common to use three skin colour categories to assess race: white, brown or “mixed”, and black [31]. In the city of Pelotas, site of the present study, 80% of the population classify themselves as white [32]; however, the population is highly admixed [33, 34]. It has been shown that those subjects who self-reported their skin colour as black had lower lung function parameters compared with those who self-reported white skin colour (table S1). Previous data from the genome-wide study in Brazil showed that self-reported skin colour has a good correlation with genomic ancestry (data not shown).
Categorising subjects into racial/ethnical groups is done on the assumption that race/ethnicity is a proxy for genetic relatedness, but this can misrepresent genetic variation and leads to confounding [8]. Yaeger et al. [35], comparing genetic ancestry and self-described race, found that self-reported race generally agreed well with the inferred genetic population cluster, but did not reveal the extent of admixture; other researchers have found a good correlation between self-defined race with inferred genetic clusters [36, 37]. However, the opposite has also been found: a poor correlation between self-reported race and genetic typing [38], with practical consequences such as lung function underestimation using predictive equations based on self-reported race alone [3].
Study strengths and limitations
Although the analysed sample was smaller than that seen at the 2012–2013 follow-up, information was available for more than 75% of the original eligible subjects. We should acknowledge the limitation of the chest circumference measurement by the photonic scanner, since this equipment measures the thorax circumference for men and women in the same way. Among women we should consider that differences in breast size and position could introduce some error in this measurement. We hope the stratification by sex in the multivariate analysis has minimised this problem.
Some strengths of this study should be mentioned: 1) the high quality of the spirometric tests; 2) the possibility of taking into account so many variables and measures due to the longitudinal design of this study and that these were measured prospectively; and 3) the very large number of markers used in genomic ancestry analysis.
The socioeconomic indicators we adjusted for in the present analysis are highly correlated, but the main objective of including them into the models was to take into account possible confounders of the association between African ancestry and lung function.
Conclusion
In summary, those individuals with high percentage of “African ancestry” presented lower lung function as compared to those with less “African ancestry”, but this was largely explained by differences in anthropometric variables (mainly height).
Footnotes
Editorial comment in Eur Respir J 2015; 45: 1529–1533 [DOI: 10.1183/09031936.00053115].
This article has supplementary material available from erj.ersjournals.com
Support statement: This article is based on data from the study “Pelotas Birth Cohort, 1982” conducted by Postgraduate Program in Epidemiology at Universidade Federal de Pelotas. The 1982 birth cohort study is currently supported by the Wellcome Trust through the programme entitled Major Awards for Latin America on Health Consequences of Population Change (Grant number: 086974/Z/08/Z). The International Development Research Center, World Health Organization, Overseas Development Administration, European Union, National Support Program for Centers of Excellence (PRONEX), the Brazilian National Research Council (CNPq), and the Brazilian Ministry of Health supported previous phases of the study. We are grateful to the Brazilian government for providing scholarships (a master's scholarship to F.P. Hartwig and a post-doctoral scholarship to G.D. Ferreira). Funding information for this article has been deposited with FundRef
Conflict of interest: Disclosures can be found alongside the online version of this article at erj.ersjournals.com
- Received June 19, 2014.
- Accepted January 2, 2015.
- Copyright ©ERS 2015
ERJ Open articles are open access and distributed under the terms of the Creative Commons Attribution Non-Commercial Licence 4.0.