Extended Haplotype in the Tumor Necrosis Factor Gene Cluster Is Associated with Asthma and Asthma-related Phenotypes
http://www.100md.com
美国呼吸和危急护理医学 2005年第9期
Channing Laboratory, Brigham and Women's Hospital, Harvard Medical School
Department of Anesthesia, Children's Hospital and Harvard Medical School
Department of Biostatistics, Harvard School of Public Health, Boston, Massachusetts
ABSTRACT
Rationale: Tumor necrosis factor is a proinflammatory cytokine found in increased concentrations in asthmatic airways. The TNF- (TNF) and lymphotoxin- (LTA) genes belong to the TNF gene superfamily located within the human major histocompatibility complex on chromosome 6p in a region repeatedly linked to asthma. The TNF position eC308 and LTA NcoI polymorphisms are believed to influence TNF transcription and secretion, respectively. Objectives: This study sought to determine whether polymorphisms in TNF or LTA, or in TNF-LTA haplotypes, are associated with asthma and asthma phenotypes. Methods: We genotyped the TNF eC308 and LTA NcoI polymorphisms, and two other haplotype-tagging polymorphisms in the TNF and LTA genes, in 708 children with mild to moderate asthma enrolled in the Childhood Asthma Management Program and in their parents. Using an extension of the family-based association tests in the PBAT program, each polymorphism was tested for association with asthma, age at onset of asthma, and time series data on baseline FEV1 % predicted, postbronchodilator FEV1 % predicted, body mass index, and log of PC20. Measurements and Main Results: Although no associations were found for the individual single-nucleotide polymorphisms, the haplotype analysis found the LTA NcoI_G/LTA 4371T/TNF eC308G/TNF 1078G haplotype to be associated with asthma and with all five phenotype groups. Conclusions: We conclude that it is unlikely that the TNF eC308 or LTA NcoI polymorphisms influence asthma susceptibility individually, but that this haplotype of variants may be functional or may be in linkage disequilibrium with other functional single-nucleotide polymorphisms.
Key Words: asthma; haplotypes; lymphotoxin- polymorphism; tumor necrosis factor
Tumor necrosis factor is a proinflammatory cytokine that is found in increased concentrations in asthmatic airways. The TNF- (TNF) and lymphotoxin- (LTA) genes are members of the TNF gene superfamily located within the human major histocompatibility complex on chromosome 6p in a region repeatedly linked to asthma (1, 2). As shown in Table 1, association between a position eC308 guanine (G)-to-adenine (A) polymorphism in the TNF promoter and asthma has been tested in 18 published studies that included different age groups of individuals with asthma from a range of ethnic backgrounds, and an association between the A allele and asthma was reported in seven of these studies (3eC9). In nine additional studies (10eC18) listed in Table 1, however, authors reported no association between asthma and this single-nucleotide polymorphism (SNP), and two other studies showed an association between the wild-type allele and asthma (19, 20). The TNF eC308A allele is of interest because it is associated with increased in vitro transcription of TNF and with increased TNF levels in stimulated human white blood cells (21, 22). As shown in Figure 1, the LTA gene, previously called the TNF- gene, is located close upstream to the TNF gene. An NcoI polymorphism in the LTA gene has also been associated with asthma in children (19) and adults (5), although each study reported that a different allele was associated and 13 other studies (4, 7, 9, 11, 13, 14, 16, 17, 20, 23eC26) showed no association (see Table 1).
Conflicting findings in previous studies could be explained by the different populations studied, by population stratification, or by small sample sizes in case-control study designs. The TNF eC308 and LTA NcoI polymorphisms are in strong linkage disequilibrium (5, 19). Moffatt and colleagues (16) have advocated that haplotypes created by the TNF eC308 and LTA NcoI polymorphisms and surrounding innate immunity genes in the major histocompatibility complex may explain the conflicting results. A Japanese study showed significant transmission disequilibrium with haplotypes in the TNF genes (26). Using a family-based design to control for population stratification in a large cohort of children with mild to moderate asthma, we tested the TNF eC308 and LTA NcoI polymorphisms, and extended TNF/LTA haplotypes for association with asthma. Some of the results of these studies have been previously reported in the form of an abstract (27).
METHODS
Population
The Childhood Asthma Management Program (CAMP) was a multicenter, randomized, double-masked, placebo-controlled clinical trial evaluating the long-term effects of inhaled antiinflammatory medications in children with mild to moderate asthma (28). The diagnosis of asthma was based on methacholine hyperreactivity (PC20 no greater than 12.5 mg/ml) and the meeting of one or more of the following criteria for at least 6 months in the year before recruitment: (1) asthma symptoms two times or more per week; (2) at least two uses per week of an inhaled bronchodilator; and (3) daily asthma medication. See the online supplement for a description of spirometry, including methacholine challenge (29), and other testing methods.
Of the 1,041 children participating in the original CAMP study (29), DNA samples were obtained from 968 participating children and 1,518 of their parents. Complete family trios were available for 652 nuclear families that included 708 children. Some families had multiple children with asthma (49 families had two children and two families had three children). The characteristics of these children are shown in Table 2.
Choice of SNPs and Molecular Methods
Choice of SNPs and SNP locations are described in detail in the online supplement, including representation of the SNPs reported as being significantly overtransmitted to Japanese individuals with atopic asthma (26). Briefly, haplotypes were identified with the PHASE program (30) and this output was analyzed by the BEST program (31) to identify haplotype-tagging SNPs. As shown in Table 3, the four SNPs chosen for genotyping should have been able to distinguish all of the common (more than 7%) haplotypes identified in Europeans. SNP genotyping was performed with the SEQUENOM system (see the online supplement for a detailed description).
Statistical Analysis
The PedCheck program was used to assess the genotype data for pedigree errors (32). Hardy-Weinberg equilibrium was tested in the parental data for each locus, using the 2 goodness of fit test in SAS (version 8.1; SAS Institute, Cary, NC).
The primary analysis was for association of individual SNPs and haplotypes with asthma, applying the PBAT (33) program to the entire sample for the SNP analysis and to white individuals for the haplotype analysis. In the PBAT program the recessive model was used because that model had optimal power, as has been shown in previous studies (34, 35). The secondary analysis was to test each individual SNP/haplotype for association with a time series of asthma-related quantitative phenotypes, using the PBAT program (36). The time series data were obtained during multiple visits (before, at, and after study entry) to children in the CAMP population (see the online supplement for details). We tested the age at onset of asthma and the following four-time series: body mass index (BMI), logarithm of the provocative concentration of methacholine causing a 20% fall in FEV1 from baseline (LNPC20), baseline FEV1 % predicted (Pre-FEV), and postbronchodilator FEV1 % predicted (Pos-FEV). The time series data were tested with the FBAT-PC statistic for repeated measurements (36) and the time-to-onset data were tested with FBAT-LOGRANK (37). Both tests are available in PBAT for SNP and haplotype analysis.
Pairwise linkage disequilibrium between each pair of SNP loci was evaluated by a maximum likelihood method to infer phases for dual heterozygotes, expressed as r2.
RESULTS
Genotyping Results and Haplotype Evaluation
Genotypes for families with pedigree check errors were set to zero. The number of PedCheck errors for each SNP were as follows: LTA NcoI (23), LTA 4371 (9), TNF eC308 (3), and TNF 1078 (3). Genotyping success rates were 97% or greater for each SNP. Hardy-Weinberg equilibrium was confirmed for all loci in the parents (all p values were greater than 0.1). There was significant linkage disequilibrium between the TNF eC308 and LTA NcoI SNPs (r2 = 0.424). The LTA NcoI SNP and the LTA 4371 SNP had an r2 value of 0.172, and the TNF 282 SNP and LTA 4371 SNP had an r2 value of 0.083. The r2 values for all other SNP pairs were less than 0.04.
There were five four-SNP haplotypes represented in the sample at 5% or greater frequency in any ethnic group. The common haplotypes represented in the sample were the same as those with 5% or greater frequency reported in Table 3 for the Centre d'Etude du Polymorphisme Humain European sample.
Family-based Association Analysis of the TNF SNPs and Haplotypes with Asthma
As shown in Table 4, the A allele for the TNF eC308 SNP was not overtransmitted to children with asthma from their parents. In contrast to previous reports of an increased frequency of the A allele in the asthmatic cohorts, we found a trend toward undertransmission of the A allele in the entire cohort. using additive or recessive models. There was also a trend for the G allele of the LTA NcoI polymorphism to be overtransmitted, but the p value (less than 0.05) was not corrected for multiple comparisons.
We restricted the extended haplotype analysis to white subjects. As shown in Table 5, there was an increase in transmission of the LTA NcoI_G/LTA 4371T/TNF eC308G/TNF 1078G (GTGG) haplotype, using the recessive model (GTGG = 167 vs. other haplotypes = 146; relative risk, 2/0 copies GTGG = 3.1 and 1/0 copy GTGG = 2.75, p = 0.05). The p values shown in Table 5 are corrected for multiple comparisons. Table 5 also shows the results of the haplotype analysis for the GTGG haplotype for asthma-related phenotypes. The GTGG haplotype was associated with later age of onset of asthma based on the Kaplan-Meier plot (Figure 2). The GTGG haplotype was also associated with BMI, LNPC20, Pre-FEV, and Pos-FEV, using the time-series phenotype data.
DISCUSSION
Our data indicate that the TNF eC308 and LTA NcoI polymorphisms do not influence asthma susceptibility individually, but that an extended haplotype in the TNF and LTA genes is associated with asthma susceptibility and with multiple common asthma phenotypes. In contrast to seven published case-control studies that showed an increased frequency of the TNF eC308A allele in individuals with asthma (3eC9), we instead found a strong trend toward an increased frequency of transmission from parents of the G allele to their child with asthma, consistent with other investigators in two case-control studies (19, 20). The TNF and LTA genes are located in the major histocompatibility complex on chromosome 6 between many innate immunity genes (see Figure 1). We believe that previous linkage with asthma in this region may be explained by extended haplotypes across the TNF and LTA genes.
Our study has many strengths. We used a family-based genetic study design including a large number of parenteCchild trios. Our population included more subjects with asthma than have been previously enrolled in other published studies evaluating SNPs in TNF and LTA genes (see Table 1). All children had mild to moderate asthma and extensive asthma phenotype data had been collected longitudinally (29). Genotyping success rates were high and we were able to infer haplotypes from family-based data at a high rate.
Our study is limited by the number of SNPs that were genotyped. Because our population was largely white, we chose SNPs to replicate findings in published reports and chose additional SNPs to ensure good representation of common haplotypes. We cannot rule out the possibility that linkage disequilibrium between polymorphisms in the asthmatic population may be different from those in the nonasthmatic population. This is unlikely, however, because our four-SNP haplotype results are consistent with those reported for the nonasthmatic European Centre d'Etude du Polymorphisme Humain sample. With the four-haplotypeeCtagging SNPs we chose, we were able to test the majority of common haplotypes across the TNF and LTA genes for association with asthma.
Noguchi and colleagues (26) reported an association between a TNF polymorphism and asthma susceptibility in a population of Japanese children with atopic asthma in 144 families. They reported increased transmission of the C allele of the TNF eC857C/T SNP to children with asthma in 54 informative families. Further evaluation of haplotypes showed significantly increased transmission of the TNF eC1031T/eC857C (our SNPs LTA 4371T/TNF 1078G) haplotype in 74 informative families. Although we substituted the TNF 1078 SNP for their SNP labeled eC857, these two SNPs are reported to be in 100% linkage disequilibrium in nonasthmatic white individuals. The TNF eC308 SNP is not found in the Japanese population (26). Our results do not contradict the finding of Noguchi and colleagues (26), but our haplotypes are more extended. The TNF distal promoter haplotype they reported is broken into subhaplotypes in our mostly white sample, and we had a much larger population of individuals with asthma with many more informative families.
Whether the TNF eC308 locus modulates the levels of tumor necrosis factor transcription is controversial. To resolve this controversy, Knight and colleagues (38) used chromatin immunoprecipitation and mass spectrometry (haplo-CHIP analysis) to identify differential proteineCDNA binding in vivo associated with allelic variants of the TNF and LTA genes. Using Epstein-Barr viruseCtransformed B-cell lines, they evaluated RNA polymerase II loading and tumor necrosis factor mRNA expression. Their haplo-CHIP analysis of the eC308 TNF polymorphism did not find differential expression of tumor necrosis factor by the two alleles of the SNP. They then extended their haplotypes to include SNPs across the LTA gene. They identified only three major haplotypes across these two genes. Two of the haplotypes had increased transmission of tumor necrosis factor (TNF eC1031T_TNFeC308A_LTA NcoIG and TNF eC1031C_TNF eC308G _LTA NcoIA) in comparison with the third haplotype (TNF eC1031T_TNF eC308G_LTA NcoIA). These haplotypes represent the following respective haplotypes in Table 3: haplotype 4, haplotype 2 or 6, and haplotype 1 or 5. Haplotype 3 (Table 3), which we found to be positively associated with asthma, was not represented in their analysis. Extrapolation from the study by Knight and colleagues (38) is limited because three common haplotypes found in asthmatic and nonasthmatic subjects are missing from their analysis. However, it is clear from their study that the TNF eC308 SNP does not upregulate tumor necrosis factor levels and is not likely to be the important SNP as previously hypothesized. We were not able to test whether haplotype 3 (Table 3) had a regulatory effect on tumor necrosis factor levels.
In a large population of children with mild to moderate asthma, we have shown that the TNF eC308 and LTA NcoI loci are not associated with asthma susceptibility or asthma phenotypes. Our results are consistent with the majority of other case-control and family-based studies shown in Table 1 that failed to find an association between asthma and these two loci. One haplotype across the TNF and LTA genes was found to be associated with asthma susceptibility and with asthma phenotypes in our cohort. Future studies of the association of polymorphisms in the TNF and LTA genes should evaluate extended haplotypes across the TNF and LTA genes and across other nearby genes in the major histocompatibility complex (5).
Acknowledgments
The authors thank all families for their enthusiastic participation in the CAMP Genetics Ancillary Study, supported by the National Heart, Lung, and Blood Institute (NO1-HR-16049). The authors also acknowledge the CAMP investigators and research team, supported by the NHLBI, for collection of CAMP Genetic Ancillary Study data. All work on data from the CAMP Genetics Ancillary Study was conducted at the Channing Laboratory, Brigham and Women's Hospital under appropriate CAMP policies and human subject protections.
This article has an online supplement, which is accessible from this issue's table of contents at www.atsjournals.org
REFERENCES
Daniels SE, Bhattacharrya S, James A, Leaves NI, Young A, Hill MR, Faux JA, Ryan GF, le Souef PN, Lathrop GM, et al. A genome-wide search for quantitative trait loci underlying asthma. Nature 1996;383:247eC250.
Xu J, Meyers DA, Ober C, Blumenthal MN, Mellen B, Barnes KC, King RA, Lester LA, Howard TD, Solway J, et al. Genomewide screen and identification of geneeCgene interactions for asthma-susceptibility loci in three US populations: collaborative study on the genetics of asthma. Am J Hum Genet 2001;68:1437eC1446.
Chagani T, Pare PD, Zhu S, Weir TD, Bai TR, Behbehani NA, Fitzgerald JM, Sandford AJ. Prevalence of tumor necrosis factor- and angiotensin converting enzyme polymorphisms in mild/moderate and fatal/near-fatal asthma. Am J Respir Crit Care Med 1999;160:278eC282.
Li Kam Wa TC, Mansur AH, Britton J, Williams G, Pavord I, Richards K, Campbell DA, Morton N, Holgate ST, Morrison JF. Association between eC308 tumour necrosis factor promoter polymorphism and bronchial hyperreactivity in asthma. Clin Exp Allergy 1999;29:1204eC1208.
Moffatt MF, Cookson WO. Tumour necrosis factor haplotypes and asthma. Hum Mol Genet 1997;6:551eC554.
Sandford AJ, Chan HW, Wong GW, Lai CK, Chan-Yeung M. Candidate genetic polymorphisms for asthma in Chinese schoolchildren from Hong Kong. Int J Tuberc Lung Dis 2004;8:519eC527.
Wang TN, Chen WY, Wang TH, Chen CJ, Huang LY, Ko YC. GeneeCgene synergistic effect on atopic asthma: tumour necrosis factor--308 and lymphotoxin--NcoI in Taiwan's children. Clin Exp Allergy 2004;34:184eC188.
Winchester EC, Millwood IY, Rand L, Penny MA, Kessling AM. Association of the TNF--308 (GA) polymorphism with self-reported history of childhood asthma. Hum Genet 2000;107:591eC596.
Witte JS, Palmer LJ, O'Connor RD, Hopkins PJ, Hall JM. Relation between tumour necrosis factor polymorphism TNF-308 and risk of asthma. Eur J Hum Genet 2002;10:82eC85.
Beghe B, Padoan M, Moss CT, Barton SJ, Holloway JW, Holgate ST, Howell WM, Mapp CE. Lack of association of HLA class I genes and TNF -308 polymorphism in toluene diisocyanate-induced asthma. Allergy 2004;59:61eC64.
Buckova D, Holla LI, Vasku A, Znojil V, Vacha J. Lack of association between atopic asthma and the tumor necrosis factor -308 gene polymorphism in a Czech population. J Investig Allergol Clin Immunol 2002;12:192eC197.
Di Somma C, Charron D, Deichmann K, Buono C, Ruffilli A. Atopic asthma and TNF-308 alleles: linkage disequilibrium and association analyses. Hum Immunol 2003;64:359eC365.
El Bahlawan L, Christensen M, Binaei S, Murphy C, Zhang Q, Quasney M. Lack of association between the tumor necrosis factor- regulatory region genetic polymorphisms associated with elevated tumor necrosis factor- levels and children with asthma. Chest 2003;123:374SeC375S.
Lin YC, Lu CC, Su HJ, Shen CY, Lei HY, Guo YL. The association between tumor necrosis factor, HLA-DR alleles, and IgE- mediated asthma in Taiwanese adolescents. Allergy 2002;57:831eC834.
Louis R, Leyder E, Malaise M, Bartsch P, Louis E. Lack of association between adult asthma and the tumour necrosis factor -308 polymorphism gene. Eur Respir J 2000;16:604eC608.
Moffatt MF, James A, Ryan G, Musk AW, Cookson WO. Extended tumour necrosis factor/HLA-DR haplotypes and asthma in an Australian population sample. Thorax 1999;54:757eC761.
Trabetti E, Patuzzo C, Malerba G, Galavotti R, Martinati LC, Boner AL, Pignatti PF. Association of a lymphotoxin gene polymorphism and atopy in Italian families. J Med Genet 1999;36:323eC325.
Zhu S, Chan-Yeung M, Becker AB, Dimich-Ward H, Ferguson AC, Manfreda J, Watson WT, Pare PD, Sandford AJ. Polymorphisms of the IL-4, TNF-, and Fc RI genes and the risk of allergic disorders in at-risk infants. Am J Respir Crit Care Med 2000;161:1655eC1659.
Albuquerque RV, Hayden CM, Palmer LJ, Laing IA, Rye PJ, Gibson NA, Burton PR, Goldblatt J, Lesouef PN. Association of polymorphisms within the tumour necrosis factor (TNF) genes and childhood asthma. Clin Exp Allergy 1998;28:578eC584.
Shin HD, Park BL, Kim LH, Jung JH, Wang HJ, Kim YJ, Park HS, Hong SJ, Choi BW, Kim DJ, et al. Association of tumor necrosis factor polymorphisms with asthma and serum total IgE. Hum Mol Genet 2004;13:397eC403.
Messer G, Spengler U, Jung MC, Honold G, Blomer K, Pape GR, Riethmuller G, Weiss EH. Polymorphic structure of the tumor necrosis factor (TNF) locus: an NcoI polymorphism in the first intron of the human TNF- gene correlates with a variant amino acid in position 26 and a reduced level of TNF- production. J Exp Med 1991;173:209eC219.
Wilson AG, Symons JA, McDowell TL, McDevitt HO, Duff GW. Effects of a polymorphism in the human tumor necrosis factor promoter on transcriptional activation. Proc Natl Acad Sci USA 1997;94:3195eC3199.
Cardaba B, Moffatt MF, Fernandez E, Jurado A, Rojo M, Garcia M, Ansotegui IJ, Cortegano I, Arrieta I, Etxenagusia MA, et al. Allergy to Dermatophagoides in a group of Spanish gypsies: genetic restrictions. Int Arch Allergy Immunol 2001;125:297eC306.
Immervoll T, Loesgen S, Dutsch G, Gohlke H, Herbon N, Klugbauer S, et al. Fine mapping and single nucleotide polymorphism association results of candidate genes for asthma and related phenotypes. Hum Mutat 2001;18:327eC336.
Izakovicova HL, Vasku A, Izakovic V, Znojil V. The interaction of the polymorphisms in transporter of antigen peptides (TAP) and lymphotoxin (LT-) genes and atopic diseases in the Czech population. Clin Exp Allergy 2001;31:1418eC1423.
Noguchi E, Yokouchi Y, Shibasaki M, Inudou M, Nakahara S, Nogami T, Kamioka M, Yamakawa-Kobayashi K, Ichikawa K, Matsui A, et al. Association between TNFA polymorphism and the development of asthma in the Japanese population. Am J Respir Crit Care Med 2002;166:43eC46.
Randolph AG, Lange C, Silverman EK, Lazarus R, Weiss ST. Association of tumor necrosis factor gene polymorphisms with asthma and asthma related phenotypes . Am J Respir Crit Care Med 2005;2:A790.
Childhood Asthma Management Program Research Group. Long-term effects of budesonide or nedocromil in children with asthma. N Engl J Med 2000;343:1054eC1063.
Childhood Asthma Management Program Research Group. Childhood Asthma Management Program (CAMP): design, rationale, and methods. Control Clin Trials 1999;20:91eC120.
Stephens M, Smith NJ, Donnelly P. A new statistical method for haplotype reconstruction from population data. Am J Hum Genet 2001;68:978eC989.
Sebastiani P, Lazarus R, Weiss ST, Kunkel LM, Kohane IS, Ramoni MF. Minimal haplotype tagging. Proc Natl Acad Sci USA 2003;100:9900eC9905.
O'Connell JR, Weeks DE. PedCheck: a program for identification of genotype incompatibilities in linkage analysis. Am J Hum Genet 1998;63:259eC266.
Lange C, Laird NM. On a general class of conditional tests for family-based association studies in genetics: the asymptotic distribution, the conditional power, and optimality considerations. Genet Epidemiol 2002;23:165eC180.
Lazarus R, Raby BA, Lange C, Silverman EK, Kwiatkowski DJ, Vercelli D, Klimecki WJ, Martinez FD, Weiss ST. TOLL-like receptor 10 genetic variation is associated with asthma in two independent samples. Am J Respir Crit Care Med 2004;170:594eC600.
Randolph AG, Lange C, Silverman EK, Lazarus R, Silverman ES, Raby B, Brown A, Ozonoff A, Richter B, Weiss ST. The IL12B gene is associated with asthma. Am J Hum Genet 2004;75:709eC715.
Lange C, DeMeo D, Silverman EK, Weiss ST, Laird NM. PBAT: tools for family-based association studies. Am J Hum Genet 2004;74:367eC369.
Lange C, Blacker D, Laird NM. Family-based association tests for survival and times-to-onset analysis. Stat Med 2004;23:179eC189.
Knight JC, Keating BJ, Rockett KA, Kwiatkowski DP. In vivo characterization of regulatory polymorphisms by allele-specific quantification of RNA polymerase loading. Nat Genet 2003;33:469eC475.
McArthur JA, Zhang Q, Quasney MW. Association between the A/A genotype at the lymphotoxin-+250 site and increased mortality in children with positive blood cultures. Pediatr Crit Care Med 2002;3:341eC344.(Adrienne G. Randolph, Chr)
Department of Anesthesia, Children's Hospital and Harvard Medical School
Department of Biostatistics, Harvard School of Public Health, Boston, Massachusetts
ABSTRACT
Rationale: Tumor necrosis factor is a proinflammatory cytokine found in increased concentrations in asthmatic airways. The TNF- (TNF) and lymphotoxin- (LTA) genes belong to the TNF gene superfamily located within the human major histocompatibility complex on chromosome 6p in a region repeatedly linked to asthma. The TNF position eC308 and LTA NcoI polymorphisms are believed to influence TNF transcription and secretion, respectively. Objectives: This study sought to determine whether polymorphisms in TNF or LTA, or in TNF-LTA haplotypes, are associated with asthma and asthma phenotypes. Methods: We genotyped the TNF eC308 and LTA NcoI polymorphisms, and two other haplotype-tagging polymorphisms in the TNF and LTA genes, in 708 children with mild to moderate asthma enrolled in the Childhood Asthma Management Program and in their parents. Using an extension of the family-based association tests in the PBAT program, each polymorphism was tested for association with asthma, age at onset of asthma, and time series data on baseline FEV1 % predicted, postbronchodilator FEV1 % predicted, body mass index, and log of PC20. Measurements and Main Results: Although no associations were found for the individual single-nucleotide polymorphisms, the haplotype analysis found the LTA NcoI_G/LTA 4371T/TNF eC308G/TNF 1078G haplotype to be associated with asthma and with all five phenotype groups. Conclusions: We conclude that it is unlikely that the TNF eC308 or LTA NcoI polymorphisms influence asthma susceptibility individually, but that this haplotype of variants may be functional or may be in linkage disequilibrium with other functional single-nucleotide polymorphisms.
Key Words: asthma; haplotypes; lymphotoxin- polymorphism; tumor necrosis factor
Tumor necrosis factor is a proinflammatory cytokine that is found in increased concentrations in asthmatic airways. The TNF- (TNF) and lymphotoxin- (LTA) genes are members of the TNF gene superfamily located within the human major histocompatibility complex on chromosome 6p in a region repeatedly linked to asthma (1, 2). As shown in Table 1, association between a position eC308 guanine (G)-to-adenine (A) polymorphism in the TNF promoter and asthma has been tested in 18 published studies that included different age groups of individuals with asthma from a range of ethnic backgrounds, and an association between the A allele and asthma was reported in seven of these studies (3eC9). In nine additional studies (10eC18) listed in Table 1, however, authors reported no association between asthma and this single-nucleotide polymorphism (SNP), and two other studies showed an association between the wild-type allele and asthma (19, 20). The TNF eC308A allele is of interest because it is associated with increased in vitro transcription of TNF and with increased TNF levels in stimulated human white blood cells (21, 22). As shown in Figure 1, the LTA gene, previously called the TNF- gene, is located close upstream to the TNF gene. An NcoI polymorphism in the LTA gene has also been associated with asthma in children (19) and adults (5), although each study reported that a different allele was associated and 13 other studies (4, 7, 9, 11, 13, 14, 16, 17, 20, 23eC26) showed no association (see Table 1).
Conflicting findings in previous studies could be explained by the different populations studied, by population stratification, or by small sample sizes in case-control study designs. The TNF eC308 and LTA NcoI polymorphisms are in strong linkage disequilibrium (5, 19). Moffatt and colleagues (16) have advocated that haplotypes created by the TNF eC308 and LTA NcoI polymorphisms and surrounding innate immunity genes in the major histocompatibility complex may explain the conflicting results. A Japanese study showed significant transmission disequilibrium with haplotypes in the TNF genes (26). Using a family-based design to control for population stratification in a large cohort of children with mild to moderate asthma, we tested the TNF eC308 and LTA NcoI polymorphisms, and extended TNF/LTA haplotypes for association with asthma. Some of the results of these studies have been previously reported in the form of an abstract (27).
METHODS
Population
The Childhood Asthma Management Program (CAMP) was a multicenter, randomized, double-masked, placebo-controlled clinical trial evaluating the long-term effects of inhaled antiinflammatory medications in children with mild to moderate asthma (28). The diagnosis of asthma was based on methacholine hyperreactivity (PC20 no greater than 12.5 mg/ml) and the meeting of one or more of the following criteria for at least 6 months in the year before recruitment: (1) asthma symptoms two times or more per week; (2) at least two uses per week of an inhaled bronchodilator; and (3) daily asthma medication. See the online supplement for a description of spirometry, including methacholine challenge (29), and other testing methods.
Of the 1,041 children participating in the original CAMP study (29), DNA samples were obtained from 968 participating children and 1,518 of their parents. Complete family trios were available for 652 nuclear families that included 708 children. Some families had multiple children with asthma (49 families had two children and two families had three children). The characteristics of these children are shown in Table 2.
Choice of SNPs and Molecular Methods
Choice of SNPs and SNP locations are described in detail in the online supplement, including representation of the SNPs reported as being significantly overtransmitted to Japanese individuals with atopic asthma (26). Briefly, haplotypes were identified with the PHASE program (30) and this output was analyzed by the BEST program (31) to identify haplotype-tagging SNPs. As shown in Table 3, the four SNPs chosen for genotyping should have been able to distinguish all of the common (more than 7%) haplotypes identified in Europeans. SNP genotyping was performed with the SEQUENOM system (see the online supplement for a detailed description).
Statistical Analysis
The PedCheck program was used to assess the genotype data for pedigree errors (32). Hardy-Weinberg equilibrium was tested in the parental data for each locus, using the 2 goodness of fit test in SAS (version 8.1; SAS Institute, Cary, NC).
The primary analysis was for association of individual SNPs and haplotypes with asthma, applying the PBAT (33) program to the entire sample for the SNP analysis and to white individuals for the haplotype analysis. In the PBAT program the recessive model was used because that model had optimal power, as has been shown in previous studies (34, 35). The secondary analysis was to test each individual SNP/haplotype for association with a time series of asthma-related quantitative phenotypes, using the PBAT program (36). The time series data were obtained during multiple visits (before, at, and after study entry) to children in the CAMP population (see the online supplement for details). We tested the age at onset of asthma and the following four-time series: body mass index (BMI), logarithm of the provocative concentration of methacholine causing a 20% fall in FEV1 from baseline (LNPC20), baseline FEV1 % predicted (Pre-FEV), and postbronchodilator FEV1 % predicted (Pos-FEV). The time series data were tested with the FBAT-PC statistic for repeated measurements (36) and the time-to-onset data were tested with FBAT-LOGRANK (37). Both tests are available in PBAT for SNP and haplotype analysis.
Pairwise linkage disequilibrium between each pair of SNP loci was evaluated by a maximum likelihood method to infer phases for dual heterozygotes, expressed as r2.
RESULTS
Genotyping Results and Haplotype Evaluation
Genotypes for families with pedigree check errors were set to zero. The number of PedCheck errors for each SNP were as follows: LTA NcoI (23), LTA 4371 (9), TNF eC308 (3), and TNF 1078 (3). Genotyping success rates were 97% or greater for each SNP. Hardy-Weinberg equilibrium was confirmed for all loci in the parents (all p values were greater than 0.1). There was significant linkage disequilibrium between the TNF eC308 and LTA NcoI SNPs (r2 = 0.424). The LTA NcoI SNP and the LTA 4371 SNP had an r2 value of 0.172, and the TNF 282 SNP and LTA 4371 SNP had an r2 value of 0.083. The r2 values for all other SNP pairs were less than 0.04.
There were five four-SNP haplotypes represented in the sample at 5% or greater frequency in any ethnic group. The common haplotypes represented in the sample were the same as those with 5% or greater frequency reported in Table 3 for the Centre d'Etude du Polymorphisme Humain European sample.
Family-based Association Analysis of the TNF SNPs and Haplotypes with Asthma
As shown in Table 4, the A allele for the TNF eC308 SNP was not overtransmitted to children with asthma from their parents. In contrast to previous reports of an increased frequency of the A allele in the asthmatic cohorts, we found a trend toward undertransmission of the A allele in the entire cohort. using additive or recessive models. There was also a trend for the G allele of the LTA NcoI polymorphism to be overtransmitted, but the p value (less than 0.05) was not corrected for multiple comparisons.
We restricted the extended haplotype analysis to white subjects. As shown in Table 5, there was an increase in transmission of the LTA NcoI_G/LTA 4371T/TNF eC308G/TNF 1078G (GTGG) haplotype, using the recessive model (GTGG = 167 vs. other haplotypes = 146; relative risk, 2/0 copies GTGG = 3.1 and 1/0 copy GTGG = 2.75, p = 0.05). The p values shown in Table 5 are corrected for multiple comparisons. Table 5 also shows the results of the haplotype analysis for the GTGG haplotype for asthma-related phenotypes. The GTGG haplotype was associated with later age of onset of asthma based on the Kaplan-Meier plot (Figure 2). The GTGG haplotype was also associated with BMI, LNPC20, Pre-FEV, and Pos-FEV, using the time-series phenotype data.
DISCUSSION
Our data indicate that the TNF eC308 and LTA NcoI polymorphisms do not influence asthma susceptibility individually, but that an extended haplotype in the TNF and LTA genes is associated with asthma susceptibility and with multiple common asthma phenotypes. In contrast to seven published case-control studies that showed an increased frequency of the TNF eC308A allele in individuals with asthma (3eC9), we instead found a strong trend toward an increased frequency of transmission from parents of the G allele to their child with asthma, consistent with other investigators in two case-control studies (19, 20). The TNF and LTA genes are located in the major histocompatibility complex on chromosome 6 between many innate immunity genes (see Figure 1). We believe that previous linkage with asthma in this region may be explained by extended haplotypes across the TNF and LTA genes.
Our study has many strengths. We used a family-based genetic study design including a large number of parenteCchild trios. Our population included more subjects with asthma than have been previously enrolled in other published studies evaluating SNPs in TNF and LTA genes (see Table 1). All children had mild to moderate asthma and extensive asthma phenotype data had been collected longitudinally (29). Genotyping success rates were high and we were able to infer haplotypes from family-based data at a high rate.
Our study is limited by the number of SNPs that were genotyped. Because our population was largely white, we chose SNPs to replicate findings in published reports and chose additional SNPs to ensure good representation of common haplotypes. We cannot rule out the possibility that linkage disequilibrium between polymorphisms in the asthmatic population may be different from those in the nonasthmatic population. This is unlikely, however, because our four-SNP haplotype results are consistent with those reported for the nonasthmatic European Centre d'Etude du Polymorphisme Humain sample. With the four-haplotypeeCtagging SNPs we chose, we were able to test the majority of common haplotypes across the TNF and LTA genes for association with asthma.
Noguchi and colleagues (26) reported an association between a TNF polymorphism and asthma susceptibility in a population of Japanese children with atopic asthma in 144 families. They reported increased transmission of the C allele of the TNF eC857C/T SNP to children with asthma in 54 informative families. Further evaluation of haplotypes showed significantly increased transmission of the TNF eC1031T/eC857C (our SNPs LTA 4371T/TNF 1078G) haplotype in 74 informative families. Although we substituted the TNF 1078 SNP for their SNP labeled eC857, these two SNPs are reported to be in 100% linkage disequilibrium in nonasthmatic white individuals. The TNF eC308 SNP is not found in the Japanese population (26). Our results do not contradict the finding of Noguchi and colleagues (26), but our haplotypes are more extended. The TNF distal promoter haplotype they reported is broken into subhaplotypes in our mostly white sample, and we had a much larger population of individuals with asthma with many more informative families.
Whether the TNF eC308 locus modulates the levels of tumor necrosis factor transcription is controversial. To resolve this controversy, Knight and colleagues (38) used chromatin immunoprecipitation and mass spectrometry (haplo-CHIP analysis) to identify differential proteineCDNA binding in vivo associated with allelic variants of the TNF and LTA genes. Using Epstein-Barr viruseCtransformed B-cell lines, they evaluated RNA polymerase II loading and tumor necrosis factor mRNA expression. Their haplo-CHIP analysis of the eC308 TNF polymorphism did not find differential expression of tumor necrosis factor by the two alleles of the SNP. They then extended their haplotypes to include SNPs across the LTA gene. They identified only three major haplotypes across these two genes. Two of the haplotypes had increased transmission of tumor necrosis factor (TNF eC1031T_TNFeC308A_LTA NcoIG and TNF eC1031C_TNF eC308G _LTA NcoIA) in comparison with the third haplotype (TNF eC1031T_TNF eC308G_LTA NcoIA). These haplotypes represent the following respective haplotypes in Table 3: haplotype 4, haplotype 2 or 6, and haplotype 1 or 5. Haplotype 3 (Table 3), which we found to be positively associated with asthma, was not represented in their analysis. Extrapolation from the study by Knight and colleagues (38) is limited because three common haplotypes found in asthmatic and nonasthmatic subjects are missing from their analysis. However, it is clear from their study that the TNF eC308 SNP does not upregulate tumor necrosis factor levels and is not likely to be the important SNP as previously hypothesized. We were not able to test whether haplotype 3 (Table 3) had a regulatory effect on tumor necrosis factor levels.
In a large population of children with mild to moderate asthma, we have shown that the TNF eC308 and LTA NcoI loci are not associated with asthma susceptibility or asthma phenotypes. Our results are consistent with the majority of other case-control and family-based studies shown in Table 1 that failed to find an association between asthma and these two loci. One haplotype across the TNF and LTA genes was found to be associated with asthma susceptibility and with asthma phenotypes in our cohort. Future studies of the association of polymorphisms in the TNF and LTA genes should evaluate extended haplotypes across the TNF and LTA genes and across other nearby genes in the major histocompatibility complex (5).
Acknowledgments
The authors thank all families for their enthusiastic participation in the CAMP Genetics Ancillary Study, supported by the National Heart, Lung, and Blood Institute (NO1-HR-16049). The authors also acknowledge the CAMP investigators and research team, supported by the NHLBI, for collection of CAMP Genetic Ancillary Study data. All work on data from the CAMP Genetics Ancillary Study was conducted at the Channing Laboratory, Brigham and Women's Hospital under appropriate CAMP policies and human subject protections.
This article has an online supplement, which is accessible from this issue's table of contents at www.atsjournals.org
REFERENCES
Daniels SE, Bhattacharrya S, James A, Leaves NI, Young A, Hill MR, Faux JA, Ryan GF, le Souef PN, Lathrop GM, et al. A genome-wide search for quantitative trait loci underlying asthma. Nature 1996;383:247eC250.
Xu J, Meyers DA, Ober C, Blumenthal MN, Mellen B, Barnes KC, King RA, Lester LA, Howard TD, Solway J, et al. Genomewide screen and identification of geneeCgene interactions for asthma-susceptibility loci in three US populations: collaborative study on the genetics of asthma. Am J Hum Genet 2001;68:1437eC1446.
Chagani T, Pare PD, Zhu S, Weir TD, Bai TR, Behbehani NA, Fitzgerald JM, Sandford AJ. Prevalence of tumor necrosis factor- and angiotensin converting enzyme polymorphisms in mild/moderate and fatal/near-fatal asthma. Am J Respir Crit Care Med 1999;160:278eC282.
Li Kam Wa TC, Mansur AH, Britton J, Williams G, Pavord I, Richards K, Campbell DA, Morton N, Holgate ST, Morrison JF. Association between eC308 tumour necrosis factor promoter polymorphism and bronchial hyperreactivity in asthma. Clin Exp Allergy 1999;29:1204eC1208.
Moffatt MF, Cookson WO. Tumour necrosis factor haplotypes and asthma. Hum Mol Genet 1997;6:551eC554.
Sandford AJ, Chan HW, Wong GW, Lai CK, Chan-Yeung M. Candidate genetic polymorphisms for asthma in Chinese schoolchildren from Hong Kong. Int J Tuberc Lung Dis 2004;8:519eC527.
Wang TN, Chen WY, Wang TH, Chen CJ, Huang LY, Ko YC. GeneeCgene synergistic effect on atopic asthma: tumour necrosis factor--308 and lymphotoxin--NcoI in Taiwan's children. Clin Exp Allergy 2004;34:184eC188.
Winchester EC, Millwood IY, Rand L, Penny MA, Kessling AM. Association of the TNF--308 (GA) polymorphism with self-reported history of childhood asthma. Hum Genet 2000;107:591eC596.
Witte JS, Palmer LJ, O'Connor RD, Hopkins PJ, Hall JM. Relation between tumour necrosis factor polymorphism TNF-308 and risk of asthma. Eur J Hum Genet 2002;10:82eC85.
Beghe B, Padoan M, Moss CT, Barton SJ, Holloway JW, Holgate ST, Howell WM, Mapp CE. Lack of association of HLA class I genes and TNF -308 polymorphism in toluene diisocyanate-induced asthma. Allergy 2004;59:61eC64.
Buckova D, Holla LI, Vasku A, Znojil V, Vacha J. Lack of association between atopic asthma and the tumor necrosis factor -308 gene polymorphism in a Czech population. J Investig Allergol Clin Immunol 2002;12:192eC197.
Di Somma C, Charron D, Deichmann K, Buono C, Ruffilli A. Atopic asthma and TNF-308 alleles: linkage disequilibrium and association analyses. Hum Immunol 2003;64:359eC365.
El Bahlawan L, Christensen M, Binaei S, Murphy C, Zhang Q, Quasney M. Lack of association between the tumor necrosis factor- regulatory region genetic polymorphisms associated with elevated tumor necrosis factor- levels and children with asthma. Chest 2003;123:374SeC375S.
Lin YC, Lu CC, Su HJ, Shen CY, Lei HY, Guo YL. The association between tumor necrosis factor, HLA-DR alleles, and IgE- mediated asthma in Taiwanese adolescents. Allergy 2002;57:831eC834.
Louis R, Leyder E, Malaise M, Bartsch P, Louis E. Lack of association between adult asthma and the tumour necrosis factor -308 polymorphism gene. Eur Respir J 2000;16:604eC608.
Moffatt MF, James A, Ryan G, Musk AW, Cookson WO. Extended tumour necrosis factor/HLA-DR haplotypes and asthma in an Australian population sample. Thorax 1999;54:757eC761.
Trabetti E, Patuzzo C, Malerba G, Galavotti R, Martinati LC, Boner AL, Pignatti PF. Association of a lymphotoxin gene polymorphism and atopy in Italian families. J Med Genet 1999;36:323eC325.
Zhu S, Chan-Yeung M, Becker AB, Dimich-Ward H, Ferguson AC, Manfreda J, Watson WT, Pare PD, Sandford AJ. Polymorphisms of the IL-4, TNF-, and Fc RI genes and the risk of allergic disorders in at-risk infants. Am J Respir Crit Care Med 2000;161:1655eC1659.
Albuquerque RV, Hayden CM, Palmer LJ, Laing IA, Rye PJ, Gibson NA, Burton PR, Goldblatt J, Lesouef PN. Association of polymorphisms within the tumour necrosis factor (TNF) genes and childhood asthma. Clin Exp Allergy 1998;28:578eC584.
Shin HD, Park BL, Kim LH, Jung JH, Wang HJ, Kim YJ, Park HS, Hong SJ, Choi BW, Kim DJ, et al. Association of tumor necrosis factor polymorphisms with asthma and serum total IgE. Hum Mol Genet 2004;13:397eC403.
Messer G, Spengler U, Jung MC, Honold G, Blomer K, Pape GR, Riethmuller G, Weiss EH. Polymorphic structure of the tumor necrosis factor (TNF) locus: an NcoI polymorphism in the first intron of the human TNF- gene correlates with a variant amino acid in position 26 and a reduced level of TNF- production. J Exp Med 1991;173:209eC219.
Wilson AG, Symons JA, McDowell TL, McDevitt HO, Duff GW. Effects of a polymorphism in the human tumor necrosis factor promoter on transcriptional activation. Proc Natl Acad Sci USA 1997;94:3195eC3199.
Cardaba B, Moffatt MF, Fernandez E, Jurado A, Rojo M, Garcia M, Ansotegui IJ, Cortegano I, Arrieta I, Etxenagusia MA, et al. Allergy to Dermatophagoides in a group of Spanish gypsies: genetic restrictions. Int Arch Allergy Immunol 2001;125:297eC306.
Immervoll T, Loesgen S, Dutsch G, Gohlke H, Herbon N, Klugbauer S, et al. Fine mapping and single nucleotide polymorphism association results of candidate genes for asthma and related phenotypes. Hum Mutat 2001;18:327eC336.
Izakovicova HL, Vasku A, Izakovic V, Znojil V. The interaction of the polymorphisms in transporter of antigen peptides (TAP) and lymphotoxin (LT-) genes and atopic diseases in the Czech population. Clin Exp Allergy 2001;31:1418eC1423.
Noguchi E, Yokouchi Y, Shibasaki M, Inudou M, Nakahara S, Nogami T, Kamioka M, Yamakawa-Kobayashi K, Ichikawa K, Matsui A, et al. Association between TNFA polymorphism and the development of asthma in the Japanese population. Am J Respir Crit Care Med 2002;166:43eC46.
Randolph AG, Lange C, Silverman EK, Lazarus R, Weiss ST. Association of tumor necrosis factor gene polymorphisms with asthma and asthma related phenotypes . Am J Respir Crit Care Med 2005;2:A790.
Childhood Asthma Management Program Research Group. Long-term effects of budesonide or nedocromil in children with asthma. N Engl J Med 2000;343:1054eC1063.
Childhood Asthma Management Program Research Group. Childhood Asthma Management Program (CAMP): design, rationale, and methods. Control Clin Trials 1999;20:91eC120.
Stephens M, Smith NJ, Donnelly P. A new statistical method for haplotype reconstruction from population data. Am J Hum Genet 2001;68:978eC989.
Sebastiani P, Lazarus R, Weiss ST, Kunkel LM, Kohane IS, Ramoni MF. Minimal haplotype tagging. Proc Natl Acad Sci USA 2003;100:9900eC9905.
O'Connell JR, Weeks DE. PedCheck: a program for identification of genotype incompatibilities in linkage analysis. Am J Hum Genet 1998;63:259eC266.
Lange C, Laird NM. On a general class of conditional tests for family-based association studies in genetics: the asymptotic distribution, the conditional power, and optimality considerations. Genet Epidemiol 2002;23:165eC180.
Lazarus R, Raby BA, Lange C, Silverman EK, Kwiatkowski DJ, Vercelli D, Klimecki WJ, Martinez FD, Weiss ST. TOLL-like receptor 10 genetic variation is associated with asthma in two independent samples. Am J Respir Crit Care Med 2004;170:594eC600.
Randolph AG, Lange C, Silverman EK, Lazarus R, Silverman ES, Raby B, Brown A, Ozonoff A, Richter B, Weiss ST. The IL12B gene is associated with asthma. Am J Hum Genet 2004;75:709eC715.
Lange C, DeMeo D, Silverman EK, Weiss ST, Laird NM. PBAT: tools for family-based association studies. Am J Hum Genet 2004;74:367eC369.
Lange C, Blacker D, Laird NM. Family-based association tests for survival and times-to-onset analysis. Stat Med 2004;23:179eC189.
Knight JC, Keating BJ, Rockett KA, Kwiatkowski DP. In vivo characterization of regulatory polymorphisms by allele-specific quantification of RNA polymerase loading. Nat Genet 2003;33:469eC475.
McArthur JA, Zhang Q, Quasney MW. Association between the A/A genotype at the lymphotoxin-+250 site and increased mortality in children with positive blood cultures. Pediatr Crit Care Med 2002;3:341eC344.(Adrienne G. Randolph, Chr)