Copyright ©ERS Journals Ltd 2007 Genetic and environmental effects on exhaled nitric oxide and airway responsiveness in a population-based sample of twins1 Dept of Respiratory Medicine, Rikshospitalet University Hospital, and 2 Division of Epidemiology, Norwegian Institute of Public Health, Oslo, Norway. CORRESPONDENCE: M. B. Lund, Dept of Respiratory Medicine, Rikshospitalet University Hospital, N-0027 Oslo, Norway. Fax: 47 23073917. E-mail: may.brit.lund{at}rikshospitalet.no Keywords: Airway responsiveness, asthma phenotypes, epigenetics, exhaled nitric oxide, methacholine, twins
Received: April 14, 2005
Elevated levels of exhaled nitric oxide (eNO) and airway hyperresponsiveness are intermediate phenotypes of asthma. Using population-based data collected from a sample of twins, the present authors estimated the relative contribution of genes, family environment and nonshared environmental influences to variations in eNO and airway responsiveness (AR). In addition, the genetic and environmental sources of covariation between these two asthma-related phenotypes were investigated. The study population comprised a random sample of 377 adult twins identified through the Norwegian Twin Registry. The main outcome variables were eNO and AR to methacholine. Genetic effects accounted for 60% of the variation in eNO. Family environment accounted for 30% of the variation in AR, while nonshared environmental influences explained the remaining variation for both measures. For both eNO and AR, there were significant regression effects for atopy and smoking. The small, but significant association between eNO and AR was primarily explained by genetic factors. Sub-analyses restricted to atopic and nonsmoking twins strengthened the observation. In conclusion, variations in exhaled nitric oxide and airway responsiveness appear to be explained by different genetic and environmental variance structures. Variation in exhaled nitric oxide is explained by genetic and nonshared environmental effects, whereas an environmental model best explains the variation in airway responsiveness. Common genetic effects explain the small but significant association between exhaled nitric oxide and airway responsiveness. Asthma is a chronic inflammatory disorder presenting a wide spectrum of clinical manifestations. Twin and family studies report genetic and environmental influences on asthma aetiology 1, 2 and molecular approaches provide insight into the genes associated with asthma and the related phenotypes 3, 4. Regions linked to these phenotypes have been detected both in candidate gene studies and genome screens 5, 6. One approach designed to help elucidate issues of aetiology is to study the genetic and environmental sources of variation and covariation in asthma-related phenotypes. Elevated levels of exhaled nitric oxide (eNO) and airway hyperresponsiveness (AHR) are intermediate phenotypes (biomarkers) associated with asthma but neither is specific to the disease. eNO is a marker of airway inflammation recognised as a critical component of the asthmatic phenotype 7. Levels of eNO are increased in asthmatics 8 and the measurement of eNO is advocated as a tool to discriminate asthmatics from nonasthmatics 9. Levels of eNO correlate with markers of asthmatic inflammation, such as sputum and serum eosinophils 10, 11. An association between eNO and eosinophilic inflammation has been reported in bronchial biopsies of patients with atopic asthma 12. The relationship between eNO and airway responsiveness (AR) is controversial and conflicting results have been published 13, 14. Measuring AR is advocated in population studies of asthma due to its close relationship with the underlying mechanisms of the disease 15. However, there is no clear evidence of a close relationship between AHR and asthmatic inflammation 16. Autonomic dysfunction or airway wall remodelling may help explain inter-individual variability in AR 17, 18. A few genetically informative epidemiological studies have assessed AR and atopy 1921, but the heritability of eNO has not yet been explored. A key question is whether measures of eNO and AR are related and, if so, whether this relationship is explained by common genetic and/or environmental pathways. The purpose of the present study was to estimate the relative contribution of genes and environment to variation in eNO and AR using population-based data from a sample of young twin adults. A further aim was to investigate the genetic and environmental sources of covariation between these two asthma-related phenotypes.
Study population The clinical data were collected as a sub-study to the study of twins carried out at the Norwegian Institute of Public Health (Oslo, Norway) 22. All twins born in Norway during the period 19671979 were identified through the Norwegian Birth Registry. The twins who participated in the current study were sampled from respondents to a questionnaire study conducted in 1998. The questionnaire was sent to 12,701 like- and unlike-sexed twins aged 1831 yrs. Responses were obtained from 3,334 (69%) complete pairs. Zygosity was determined on the basis of seven questions previously validated to correctly categorise >97% of the twins in another Norwegian study 23. The sample and zygosity classification procedure are described elsewhere 23, 24. A random sample of twins living within a 2-h transportation radius from Oslo was invited to participate in the present study. Exclusion criteria were pregnancy and medical conditions that might interfere with carrying out the tests. In total, out of 1,446 twins (723 pairs) selected, 643 twins were willing to participate. The clinical study was completed on 377 twins comprising 171 pairs and 35 twins whose co-twin was not tested. The participants in the clinical study did not differ from the nonparticipants with respect to sex, age or zygosity. The regional medical ethics committee approved the study and written informed consent was obtained from all subjects.
Measures
eNO
Baseline lung function
AR
Atopy
Analyses The relationship between eNO and AR was then analysed using a bivariate ACE model that decomposed the covariance between these measures into correlated genetic (rG), shared environmental (rC) and nonshared environmental (rE) pathways. A full model estimating all three latent factor correlations was analysed first and subsequent models tested for the significance of each correlation by dropping it from the model.
The raw data option in Mx calculates minus twice the log-likelihood (-2LL) of the data but does not provide an overall measure of fit. Model comparisons were conducted using the likelihood ratio Chi-squared test (
AIC = which combines parsimony and fit. Prior to model specification, the pattern of MZ versus DZ intra-class correlations were inspected to determine which factors (A, C, D, E) were important to parameterise. Greater MZ values than DZ values indicate genetic effects and MZ correlations that exceed twice the DZ value suggest nonadditive (D) effects. Common environmental effects are indicated if the DZ correlation is greater than twice the MZ correlation 31.
Clinical characteristics and measures Table 1
Univariate results for eNO and AR Tests of the modelled covariates revealed no effect for age. For eNO there were significant regression effects for smoking (change in Chi-squared tests between nested models with one degree of freedom ( ![]() 2(1)) = 18.71, p<0.001) and atopy (![]() 2(1) = 9.99, p<0.01). AR regressions were significant for sex (![]() 2(1) = 14.36, p<0.001), smoking (![]() 2(1) = 16.09, p<0.001) and atopy (![]() 2(1) = 9.68, p<0.01). The intra-class correlations for the MZ and DZ twins were 0.57 and 0.28 for the log score of eNO and 0.35 and 0.32 for the log score of DRS, respectively. These values suggest that different patterns of genetic and environmental influences explain the variations in eNO and AR. Specifically, genetic and nonshared environmental effects are indicated for variation in eNO, whereas shared and nonshared environmental influences seem more important for variation in AR. Given the lack of evidence for genetic dominance, ACE models and sub-models were tested for both measures. Results including the genetic and environmental variance estimates, 95% confidence intervals (CI) and model fit statistics are presented in tables 2![]() 2(1) = 0.00, p = 0.99) or C (![]() 2(1) = 1.24, p = 0.27) did not significantly affect the model fit, but that one of these sources of shared effects (A or C) must be retained. The E-only model showed a significant reduction in fit (change in Chi-squared tests between nested models with two degrees of freedom (![]() 2(2)) = 15.35, p<0.01). The completely environmental model (CE) is the most consistent with the pattern of intra-class correlations and also yielded the lowest AIC value. Accordingly, a CE model best explains variation in AR, with common effects accounting for 30% of the variation and nonshared environmental factors explaining the remaining variation. In contrast, the results for eNO revealed highly significant genetic effects. Dropping A from the model resulted in a significant deterioration in fit (![]() 2(1) = 10.13, p 0.01). However, there was no effect of shared environment (![]() 2(1) = 0) as evidenced by no change in -2LL when dropping C from the model. The AE model showed the best fit according to AIC with genetic and nonshared environmental influences for 60 and 40% of the variation, respectively.
Bivariate analyses of the relationship between eNO and AR The association between eNO and AR is weak but significant (r = 0.14, p = 0.006). The cross-twin cross-trait correlations (eNO in twin A with AR in twin B) were 0.23 for MZ and 0.04 for DZ twins. This pattern suggests that genetic effects contribute importantly to the association between eNO and AR. Bivariate analyses revealed that a model with no covariation was unfit for the data. Therefore, tests of the importance of each source of covariation (rG, rC and rE; M1) were conducted and results are listed in table 4 ![]() 2 = 5.818, df = 1 and p = 0.016; whereas dropping rC hardly altered the -2LL (![]() 2 = 0.544, df = 1, p = 0.461). Collectively, results showing a significant change in fit only when genetic sources of covariation are excluded suggest that association between eNO and AR is primarily explained by common genetic effects.
This population-based study of young adult twins revealed that the variation in eNO and AR appeared to be explained by different genetic and environmental variance structures. Whereas variation in eNO was explained by genetic and environmental factors, variation in AR was mainly attributed to environmental effects. The small but significant association between eNO and AR was primarily explained by common genetic factors. This finding may seem inconsistent given that an environmental model best described the AR data. However, the CI were wide in the univariate results and power to detect potential genetic effects is increased in the bivariate analysis. Genetic studies using the twin-design are based upon the assumption that twins are representative of the general population for the outcomes being studied. Comparisons with prevalence rates from non-twin population-based Norwegian surveys suggest that the twin data in the current study are representative of the general population with respect to atopy, AHR 32 and smoking habits 33. The prevalence of atopy was 35%. Similar prevalence rates have been reported in other Scandinavian studies of young adults and adolescents 32, 34. In the current study, the 21% prevalence of AHR is in line with data previously reported from a population survey of Norwegian adolescents 32. The observation that AHR is more prevalent in females than in males has also been reported in other Scandinavian studies 35, 36. No evidence was found of zygosity differences in the prevalence of any of the clinical measures. Given the similarity in the reported prevalence of asthma between the sample who participated in the clinical testing and the population-based sample of twins responding to the main questionnaire, it is unlikely that selection bias presents a major problem in the present study. In recent years, eNO has become recognised as an important component of the asthmatic phenotype 7. Endogenous NO is derived from L-arginine by the enzyme NO synthase (NOS), of which at least three distinct isoforms have been identified. Two of these are constitutively expressed (NOS1 and NOS3), while a third (NOS2) is inducible and shows increased expression in inflammatory disorders. Genetic polymorphisms of all three NOS isoforms have been identified and associations have been reported between NOS1 gene polymorphisms and asthma 37. In asthmatic subjects, an association between the length of the AAT repeat polymorphism in intron 20 of the NOS1 gene and levels of eNO has been reported 37. Data on the heritability of eNO has not previously been published. In the present population-based sample of twins, genetic effects were found that accounted for 57% of the variation in eNO. Numerous studies have documented a strong association between atopy and eNO 7, 10, 32. Since current knowledge indicates that atopy is genetically influenced 36, a common heritable origin for the two phenotypes may be suggested. Since the present study was conducted outside the pollen allergy season, the only individuals likely to be tested during exposure are those with indoor allergies such as house dust mites or pets. Ideally, the modelling analyses should be stratified by twin-pair discordance and concordance for atopy. Such analyses would provide insight into whether genes and environments differentially mediate the relationship between eNO and AR contingent upon atopy. Although the current sample is not statistically powered for such stratified analyses, the association between eNO and AR is significantly larger (0.27) among those with atopy than among those without (-0.07). Furthermore, inspection of the intra-class and cross-twin cross-trait correlations (eNO in twin 1 with AR in twin 2) in groups stratified by zygosity and atopy status of the pair (data not shown) suggests that genetic effects are highly important in explaining the relationship between eNO and AR among atopic subjects. Due to sample size limitations, these results should be interpreted with caution but warrant further investigation into the differences in the nature of genetic influences affecting eNO and AR when associated with atopy. Compared with data on the genetics of atopy, data on the genetics of AR are scarce and incomplete. However, most current literature comprising data from twin studies, segregation analyses, linkage studies and genome-wide screens suggest genetic influences on AR 6, 20, 21, 38. In 1984, Hopp et al. 38 conducted a study of 107 twin pairs aged 631 yrs who were tested for AR to methacholine. Heritability, calculated as twice the difference between the MZ and DZ correlations, was 0.66. An Australian study 20 using a population-based sample of nuclear families reported a much lower heritability (0.30, SE = 12.3) for DRS to methacholine challenge than in the study by Hopp et al. 38. Finally, another Australian study 21 of twins aged 818 yrs included pairs who were concordant and discordant for current wheeze status. Although heritability was not estimated, odds ratios were consistent with genetic influences on AR as measured by responsiveness to hypertonic saline. In contrast, evidence of significant heritability for AR was not found. Heritability is a sample statistic and is expected to vary based on characteristics of the sample that are associated with the expression of genetic and environmental effects. For example, the present authors found significant main effects for smoking and atopy, but this does not reveal potential differences in the genetic and environmental variance components as a function of smoking or atopy status. Furthermore, differences in the sampling schemes, age group, prevalence and analytic methods between the Australian and Norwegian study could contribute to discrepancies regarding heritability estimates of AR. The Australian results based upon patient populations suggest that genetic effects may become expressed once asthma has developed, but may not be contributing to variation in AR when the disease is not manifest. If this is the case, then the heritability estimates may be sensitive to differences in disease prevalence in the investigated samples. The prevalence of atopy (58%) and asthma (26%) in the Australian study 21 was much higher than that observed in the general Norwegian population. Furthermore, whereas the present authors have studied young adults, children comprised a large percentage of the Australian twin samples 21, 38. AR is known to be significantly associated with age, with AHR being more prevalent in children than in adults 39. Finally, AR is also influenced by smoking 40. The present twin sample comprised 39% current and ex-smokers, while the proportion of smokers in the Australian studies (not specifically stated) would be negligible or nil. It seems conceivable that mechanisms responsible for AR in adult smokers would differ from those in allergic children. Although the mechanisms are unknown, cigarette smoking reduces eNO and increases AR, and it is reasonable to believe that cigarette smoking could influence the genetic environmental variance structure for both eNO and AR. Due to sample size restrictions, the present authors were unable to stratify the quantitative genetic analyses according to pair-wise smoking status. However, inspection of the twin intra-class correlations among the pairs, in which both twins were currently nonsmokers, revealed that similar conclusions to those reported herein would be drawn about the variance structure for NO. This is, that genetic effects may be slightly greater (rMZ = 0.70, rDZ = 0.18) but are important for variation in AR (rMZ = 0.50 and rDZ = 0.03). Furthermore, the cross-twin cross-trait correlations (0.35 among the MZ and -0.07 among the DZ nonsmoking pairs) strengthens the conclusion that a genetic correlation explains the relationship between AR and eNO. The exposure effect of cigarette smoking on NO and AR should be investigated more thoroughly using a larger sample. The present bivariate results suggest that common sets of genes account for the small but significant association between eNO and AR. Due to the fact that other studies analysing genetic covariance structures for asthma-related phenotypes have not included measures of eNO, the present results cannot be directly compared with other findings but are consistent with those based on a multivariate analysis indicating genetic effects common to asthma, atopy and AR 20. Collectively, these studies are beginning to elucidate questions of aetiology and help guide linkage studies. In conclusion, in a population-based sample of young adult twins, variation in exhaled nitric oxide was explained by genetic and nonshared environmental effects, whereas variation in airway responsiveness was best accounted for by environmental influences. Common genetic effects explained the small but significant association between exhaled nitric oxide and airway responsiveness.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||