Abstract
Decreased spirometric indices are characteristic of asthma and other respiratory diseases. The aim of this study was to investigate the genetic and environmental components of variance of forced expiratory volume in one second (FEV1) and forced vital capacity (FVC) measured in adulthood in an Australian population-based sample of 468 Caucasian nuclear families. The inter-relationships of the genetic determinants of these traits with asthma and atopic rhinitis were also investigated.
Serial cross-sectional studies were conducted in the town of Busselton in Western Australia between 1966 and 1981 and follow-up of previous attendees was undertaken in 1995. Data from each subject included in this study were from a single survey in adulthood (25–60 yrs of age) when the subject was as close to age 45 yrs as possible.
Multivariate analysis suggested that FEV1 and FVC levels were associated with age, sex, height, tobacco smoke exposure, asthma and atopic rhinitis. After adjustment for relevant covariates, FEV1 levels had a narrow-sense heritability (h2N) of 38.9% (SE 9.1 %). FVC levels had an h2N of 40.6% (SE 8.9%). Extended modelling demonstrated little overlap in the genetic determinants of asthma or atopic rhinitis and either FEV1 or FVC levels.
The results of this study were consistent with the existence of important genetic determinants of adult lung function that are independent of asthma or other atopic disease, cigarette smoking, height, age or sex.
- Busselton Health Study
- familial aggregation
- forced expiratory volume
- forced vital capacity
- heritability
- lung function
This study was supported by the National Health and Medical Research Council, the Australian-American Fulbright Educational Foundation and by Healthway Western Australia.
Spirometric indices measure parameters of respiratory function in an individual which may reflect underlying pathological factors resulting in airflow obstruction. The most commonly used spirometric indices are the forced expiratory volume in one second (FEV1) and the forced vital capacity (FVC). FEV1 and FVC are strongly related to the severity of respiratory symptoms, atopy and elevated serum immunoglobulin E (IgE) levels in children and adults 1–4.
Population-based studies indicate that FEV1 and FVC are approximately normally distributed in the general population 5, suggesting that multiple factors are involved in determining these traits. Pedigree-based studies of unselected, asthmatic and chronic obstructive pulmonary disease (COPD) families give consistent evidence for familial aggregation of spirometric indices 6–11, suggesting that around 20–60% of total phenotypic variance may be accounted for by familial factors. Estimates of the narrow-sense heritability (h2N) of the most commonly investigated spirometric measure, FEV1, have ranged from 28% 12 to 47% 13. Twin studies also give consistent evidence of a genetic contribution to the variability of spirometric indices and a significantly higher concordance for these indices in monozygotic (MZ) twins 14–16, although estimates of the broad-sense heritability of spirometric indices in twin studies have ranged from just over 0% to almost 100% 8, 12, 14–16.
Few studies of general population samples with reasonable sample sizes have adequately adjusted for all known potential confounders such as sex, age, race, tobacco smoke exposure and body size. Those that have done so have reported inconsistent results 9, 17, 18.
There is substantial evidence that asthma and other atopic diseases aggregate within families and that a substantial proportion of this aggregation is due to genetic factors 19. The extent to which the close relationship of asthma with decreased spirometric indices reflects shared genetic determinants is unclear.
The aim of this study was to use variance components analysis to estimate the genetic and environmental components of lung function variance measured in adulthood in a population-based sample of nuclear families. A further aim was to investigate the inter-relationships between any genetic determinants of these traits with asthma and atopic rhinitis.
Methods
Study population
Busselton is a rural town on the coast of South-Western Australia and the population is predominantly of European origin. Six cross-sectional surveys of adults in the town of Busselton were undertaken, approximately one every three years, over the period 1966–1981 and a follow-up survey of all survivors from these cross-sectional surveys was conducted in 1994/1995. A wide range of health related data was gathered in each survey, including demographic variables, general health and lifestyle variables, health history variables, and physical, biochemical, haematological and immunological measurements. General descriptions of the cross-sectional surveys have been reported previously 20, 21.
The current study is based on an analysis of the 468 nuclear families for whom lung function and other data were available for each family member from ≥1 survey when they were aged 25–60 yrs. If a family member participated in >1 survey during this age range, then the data from the survey when their age was closest to 45 yrs was used. Data came from surveys conducted in 1966 (10.9%), 1969 (16.3%), 1972 (10.0%), 1975 (7.4%), 1978 (5.8%), 1981 (9.2%) and 1994/1995 (40.4%).
The ongoing cross-sectional studies were approved by the Human Rights Committee of the University of Western Australia, and informed personal consent was obtained from all subjects.
Questionnaire
Individual and family histories of respiratory symptoms, demographic information and smoking were completed at interview using a modified British Medical Research Council questionnaire 22. Doctor diagnosed asthma (ever) was defined as a positive response to the question: “Has your doctor ever told you that you have asthma/bronchial asthma?” Doctor diagnosed atopic rhinitis (ever) was defined as a positive response to the question: “Has your doctor ever told you that you have hay fever?”. Smoking was classified as never-smoked, exsmoker, light smoker (<15 cigarettes·day-1) or heavy smoker (≥15 cigarettes·day-1).
Spirometry
In the surveys performed 1966–1978, FEV1 and FVC were measured using a McDermott dry spirometer (Pneumoconiosis Research Unit, Penarth, Wales, UK) calibrated daily with a 3-L syringe. All values obtained were corrected to body temperature and pressure, saturated (BTPS) assuming a fixed room temperature and atmospheric pressure. In the survey performed in 1981, spirometry was measured using wedge spirometers (Vitalograph, Buckingham, UK) and in 1994/1995 using a pneumotachograph spirometer (Welch Allyn, Skaneateles Falls, USA). At all time-points, the best FEV1 and FVC was measured according to the guidelines of the American Thoracic Society 23. Results were recorded as the highest values from three maximum expiratory manoeuvres, provided that the best two recordings were within 5% of each other.
Because different methods and circumstances of lung function measurement were used in the different surveys, possible measurement biases in spirometric measures across the successive surveys were investigated. This analysis suggested a systematic study bias in the measurement of FEV1 and FVC in the 1969 and 1978 surveys (data not shown; the bias was unrelated to sex, smoking, height or other variables measured at the time of the studies). Therefore, FEV1 and FVC assessed in 1969 and 1978 were appropriately adjusted in order to correct for the apparent measurement bias. Height of subjects was measured at the time of spirometric assessment using a stadiometer.
Statistical analysis
The primary response variables modelled were FEV1 and FVC. Explanatory variables included sex, age, height, and smoking status. Asthma and atopic rhinitis status were also included as explanatory covariates in certain models. All explanatory variables except sex, asthma status, atopic rhinitis status and smoking were analysed as continuous covariates. Smoking was analysed as a categorical variable (never smoked=0; exsmoker=1; current light smoker=2; current heavy smoker=3). The marginal distributions of FEV1 and FVC levels were approximately normal. All continuous covariates were centred at or close to their mean. Bivariate analyses were performed using unpaired t-tests (two-tailed, equivariance not assumed) and analysis of variance (ANOVA) 24.
The software package FISHER (www . biomath . medsch . ucla . edu / faculty / klange / software .html) 25 was used to undertake multivariate modelling and to partition observed phenotypic variance into genetic and nongenetic components by maximum likelihood methods. Each model assumed that the distribution of the response phenotype for a family was multivariate normal, with a mean that depended upon a particular set of explanatory covariates. The mean models and specification of variance and covariance structures are given in the appendices. Similar modelling methods were used in the investigation of familial aggregation of cardiovascular risk factors in Busselton families 21. The model fitting procedure included an investigation of interaction or polynomial terms and examination of the effect of observations with high regression leverage. Definitive final models were shown to provide a valid summary of the observed data.
Phenotypic variance was partitioned into four components: 1) additive genetic effects (σ2A) (the additive effects of genes); 2) common family environment (σ2C) (environmental exposures shared by an entire family, e.g. ownership of a cat or dog or consumption of the same food); 3) common sibling environment (σ2CS) (environmental exposures unique to siblings over and above the common family environment, e.g. sharing a bed or a bedroom); and 4) the residual variance (σ2E) (which is assumed to arise from nonfamilial random factors, e.g. occupational exposure of a father to environmental pollutants).
The narrow-sense heritability (h2N) was defined as the ratio of variance due to additive genetic effects (σ2A) to the total phenotypic variance of each trait 26:
Asymptotic standard errors for h2N were obtained by reparameterizing the covariance model in terms of h2N rather than σ2A.
The statistical associations of covariates entered as fixed effects and the current response variable were assessed by removal of terms from the mean model and calculation of a likelihood ratio Chi-squared test statistic 26. The same approach was used as an approximate guide to the “significance” of a departure of the value of a variance component from its null value (zero). Statistical significance was taken as the p≤0.05.
Standard goodness-of-fit tests to check overall validity of models were performed using the FISHER program 25.
Extended modelling
In order to investigate the extent to which additive genetic effects were shared with asthma and atopic rhinitis, the original mean model for each phenotype of interest (e.g. FEV1 levels) was extended by adding (one at a time) terms for these other outcomes of interest (e.g. asthmatic status). Under such circumstances, a large reduction in the magnitude of σ2A suggested a sharing of additive genetic factors, and therefore sharing of genetic determinants.
Because this was not a standard nested model problem, an approach based upon profile likelihoods 27 was used to assess the magnitude of any change in the estimate of the genetic component of variance (σ2A). σ2Ao was defined as the estimated value of σ2A in the original model and σ2AE as its equivalent in the extended model. The extended model was then refitted and the value of σ2A constrained so it was forced to take the value σ2Ao. The ratio of the likelihoods of the free and constrained models measured the plausibility of the hypothesis that the true value of σ2A was unchanged by extending the model. In a standard nested model problem, a likelihood ratio of 6.82 approximates to p=0.05. More generally, a likelihood ratio in excess of 10 is moderate evidence that the hypothesis of shared determinants is plausible, and one in excess of 100 is strong evidence.
Results
Characteristics of families
Table 1⇓ shows the characteristics of the family members; data shown were collected at the survey when each individual was aged 25–60 yrs and was closest to age 45 yrs. The mean number of children per family was 2.0. The number of offspring ranged from 1–7 and there were a total of 938 offspring with available data. Males and females were equally represented in the study population.
Effects of age, sex, cigarette smoking and height
Average FEV1 and FVC levels were higher for males (fathers and sons) than females (mothers and daughters), and higher for offspring (sons and daughters) than parents (mothers and fathers). Conversely, average % predicted FEV1 and FVC levels were lower for males (fathers and sons) than females (mothers and daughters) (table 2⇓). The prevalence of current smoking was higher in parents than offspring and higher in males than females. Asthma was more prevalent in offspring compared with parents. Both asthma and atopic rhinitis were more prevalent in females than males (table 1⇓). Bivariate analysis by ANOVA indicated that both FEV1 % pred (F3,1870=43.40, p<0.0001) and FVC % pred (F3,1870=43.40, p<0.0001) were significantly lower in the group of heavy smokers (table 2⇓). FEV1 % pred were significantly lower in asthmatics (T213=7.10, p<0.0001) and in those with atopic rhinitis (T843=2.63, p=0.01) (table 2⇓). Similarly, FVC % pred were significantly lower in asthmatics (T216=4.82, p< 0.0001) and in those with atopic rhinitis (T859=2.06, p=0.04) (table 2⇓).
The multivariate variance components modelling confirmed these relationships, and indicated that age had a significant effect on FEV1 levels in both males and females (Chi-squared1=275.6, p<0.0001 and Chi-squared1=204.2, p<0.0001, respectively) and on FVC levels (Chi-squared1=95.2, p<0.0001 and Chi-squared1=85.3, p<0.0001 respectively). FEV1 and FVC were lower in older males and females. Height also had a significant effect in both males and females on FEV1 levels (Chi-squared1=275.5, p<0.0001 and Chi-squared1=143.9, p<0.0001, respectively) and on FVC levels (Chi-squared1=414.6, p<0.0001 and Chi-squared1= 241.6, p<0.0001, respectively). Levels were higher in taller males and females. Cigarette smoking was closely associated in both males and females with reduced levels of FEV1 (Chi-squared3=75.7, p<0.0001 and Chi-squared3=28.3, p<0.0001, respectively) and FVC (Chi-squared3=19.0, p=0.0001 and Chi-squared3=12.5, p= 0.003, respectively). The overall sex effect (see mean model in appendices) was significant for both FEV1 (Chi-squared6=225.9, p<0.0001) and FVC (Chi-squared6=317.9, p<0.0001) levels.
Association with asthma and atopic rhinitis
Doctor diagnosed asthma was associated with decreased FEV1 (Chi-squared1=77.7, p<0.0001) and FVC (Chi-squared1=23.6, p<0.0001). Doctor diagnosed atopic rhinitis was also associated with decreased FEV1 (Chi-squared1=13.3, p=0.0001) and FVC (Chi-squared1=5.1, p=0.01). These associations were independent of age, sex, smoking status, height and familial correlations.
Variance components and heritability estimates
After adjustment for all covariates, h2N of FEV1 was 38.9% (SE 9.1%), i.e. additive genetic effects (σ2A) contributed around two fifths of the total variance (fig. 1⇓ and model 1, table 3⇓). σ2A was significantly greater than zero (Chi-squared1=17.4, p<0.0001). The remaining variance was largely the result of nonfamilial environmental effects (σ2E). Familial environmental effects (σ2CS and σ2C) were not significantly different from zero.
The h2N of FVC was estimated to be 40.6% (SE 8.9%), i.e. σ2A contributed approximately two-fifths of the total variance (figure 2⇓ and model 1, table 4⇓. σ2A was significantly greater than zero (Chi-squared1=19.4, p<0.0001). Environmental effects common to families (σ2C) were not significantly different from zero and contributed only minimally, if at all, to total phenotypic variance. Environmental effects common to siblings (σ2CS) were significantly greater than zero (Chi-squared1=5.9, p=0.009). The majority of phenotypic variance was attributable to nonfamilial environmental effects (σ2E).
The extended modelling indicated that adjustment of the FEV1 for asthma status resulted in a small fall (to mean±sem 0.086±0.022), 89.6% of baseline model, see model 2 in table 3⇓) in the σ2A estimate. This was consistent with an ∼10% overlap in genetic (σ2A) determinants, and was associated with a likelihood ratio of 1.23, consistent with little or no change in σ2A. Adjustment of the FVC model for asthma status resulted in an estimated 0.8% fall (to 0.125±0.028), 99.2% of baseline model, see model 2 in table 4⇓) in the estimate of σ2A. The likelihood of the unconstrained model was only 1.12 times greater than that of the constrained model, consistent with little or no change in σ2A.
Adjustment of the FEV1 for atopic rhinitis status resulted in a small fall (to 0.094±0.023), 97.9% of baseline model, see model 3 in table 3⇓) in the σ2A estimate. This was consistent with an ∼2% overlap in genetic (σ2A) determinants, and was associated with a likelihood ratio of 1.08, consistent with little or no change in σ2A. Adjustment of the FVC model for atopic rhinitis status resulted in an estimated 0.8% fall (to 0.125±0.028), 99.2% of baseline model, see model 3 in table 4⇓) in the estimate of σ2A. The likelihood of the unconstrained model was only 1.10 times greater than that of the constrained model, consistent with little or no change in σ2A. The goodness-of-fit tests did not indicate any significant lack-of-fit problems for any of the models reported.
Discussion
The present study was designed to recruit a sample of families containing individuals assessed in adulthood and who were representative of a general Caucasian population. It has shown that adult FEV1 and FVC levels are strongly heritable traits, each with genetic determinants that are distinct from asthma and atopic rhinitis.
The prevalence rates for self-reported asthma and wheeze for parents and children in this study are similar to those previously reported in other Australian studies 28, 29. The higher prevalence of doctor diagnosed asthma in the offspring (table 1⇓) reflects the widely observed rising prevalence rates of asthma over the past two decades in developed nations 30, 31. Surveys of random samples have previously shown that the prevalence of wheeze and diagnosed asthma in the Busselton community increased substantially over the decade 1981–1990 32. The relationships of FEV1 and FVC levels to age, sex and cigarette smoking in the Busselton population were consistent with previous population−based studies of Caucasians 28–31, 33. The smoking rates for parents and offspring (table 1⇓) were similar to those derived from relevant demographic estimates for the Australian population 34.
A problem with some previous studies of spirometric indices has been inadequate adjustment for the effects of the many factors which might bias a heritability estimate. Such factors include height, sex, age, race and tobacco smoke exposure 12. Because these factors aggregate in families, it is important that they are adjusted for in genetic epidemiological analyses which attempt to define unique genetic determinants of the spirometric measures themselves. The current study adjusted for height, sex, age, and active tobacco smoke exposure. It was unnecessary to adjust for race, as all members of the Busselton population included in the current study were Caucasian. Thus, the present study provides evidence that FEV1 and FVC levels may be under some degree of significant genetic control, independent of genetic factors influencing body size, sex-specific expression and susceptibility to environmental tobacco smoke.
In the nuclear family design, estimates of the effects of common sibling environment are completely confounded with the estimated effects of genetic dominance (nonadditive effects of alleles at the same locus) 26. The heritability estimates presented in this paper are therefore conservative, as the numerator includes only additive genetic effects. However, scrutiny of tables 3 and 4⇓⇓ suggests that even the extreme assumption that all of the σ2CS estimate was in fact due to dominance genetic effects, would have made little or no difference to the substantive conclusions.
Both pedigree- and twin-based studies give consistent evidence of a genetic contribution to the variability of lung function 6–11, 13–16, 18, 35–37. Estimates of the broad-sense heritabilities of FEV1 and FVC derived from twin studies that have adjusted for body size, have ranged 0–∼75%. In agreement with previous studies reporting a substantial heritability of FEV1 and FVC levels 11, 18 the present results suggest substantial correlations in adult FEV1 and FVC levels between genetically related. Most of the observed familial correlation's in adult FEV1 and FVC were attributable to genetic rather than environmental factors. However, the majority of the total variance in adult FEV1 and FVC levels was due to nonfamilial environmental factors.
In adults and children, FEV1 and FVC levels are strongly associated with asthma 1–3. The present analysis allowed estimation of the overlap between genetic determinants of the spirometric measures studied and asthma and atopic rhinitis. In the present study, the observed close relationships of these spirometric measures to asthma were consistent with previous reports. However, little evidence was found that the association of asthma or atopic rhinitis with both FEV1 and FVC levels was due to the sharing of additive genetic determinants. Conversely, there was evidence of genetic and environmental determinants influencing FEV1 and FVC which were unshared with asthma and atopic rhinitis. The unshared genetic determinants suggest the existence of a distinct genetic pathway modulating FEV1 and FVC independently of asthma or other atopic disease (atopic rhinitis), age-, height- or sex-dependent expression or genetic susceptibility to tobacco smoke exposure. Independence of the genetic determinants of FEV1 and FVC from asthma has not been previously demonstrated.
The present study suggests the presence of important genetic determinants of two pathophysiological traits associated with asthma COPD and lower respiratory symptoms and further shows for the first time that FEV1 and FVC levels in adulthood are traits that are genetically distinct from asthma and atopic rhinitis. COPD is a disease that like asthma is currently the subject of intensive genetic investigation by many groups internationally and which makes extensive use of FEV1 as an intermediate phenotype 38, 39.
The results of this study, which suggest that immunoglobulin E mediated inflammatory processes are genetically distinct from spirometric indices, therefore also have important implications for genetic studies of chronic obstructive pulmonary disease. Programmes of gene identification in asthma will be facilitated by the recognition that forced expiratory volume in one second and forced vital capacity are not proxies for “asthma” in adults; different genetic mechanisms are likely to regulate these spirometric measures and susceptibility to asthma or atopy. The specific genes regulating these spirometric measures remain to be defined.
Components of total variance in forced expiratory volume in one second (FEV1). The following data are presented as estimated mean±sem (% of total variance): □ : residual error variance (σ2E), 0.146±0.014 (59.1%); ┘ : additive genetic variance (σ2A), 0.096±0.023 (38.9%); : common family environmental variance (σ2C), 0.005±0.011 (2.0%); note that common sibling environmental variance (σ2CS) is zero.
Components of total variance in forced vital capacity (FVC). The following data are presented as estimated mean± sem (% of total variance): □ : residual error variance (σ2E), 0.153±0.019 (49.4%); ┘ : additive genetic variance (σ2A), 0.126± 0.028 (40.6%); : common family environmental variance (σ2C), 0.004±0.014 (1.3%); └ : common sibling environmental variance (σ2CS), 0.027±0.012 (8.7%).
Characteristics of family members
Values for % predicted forced expiratory volume in one second (FEV1) and forced vital capacity (FVC) levels in asthmatics and smokers
Maximum likelihood models showing variance components estimates for forced expiratory volume in one second (FEV1) at baseline and after sequential adjustment for covariates
Maximum likelihood models showing variance components estimates for forced vital capacity (FVC) at baseline and after sequential adjustment for covariates
Appendices
1. The mean models were specified as follows:
Where βn with the nth fixed regression coefficient and δ is a binary indicator variable taking the value 1 in males and 0 in females. This model permitted independent covariate profiles for the phenotype of interest in males and females.
2. The total phenotypic variance (conditional on the mean model) was based on a conventional covariance structure 26 and was specified as:
3. The conditional covariances within a family were specified as:
(i) 0.5σ2A+σ2CS+σ2C between two siblings;
(ii) 0.5σ2A+σ2C between a parent and a child; and
(iii) σ2C between two parents.
Acknowledgments
The authors thank the people of the Busselton community for their participation in this study, the Busselton Population Medical Research Foundation and the many colleagues who assisted in the collection of this data.
- Received June 14, 2000.
- Accepted October 3, 2000.
- © ERS Journals Ltd