|
|
||||||||
1 Dept of Biostatistics and Research Epidemiology, Henry Ford Health System, Detroit, Michigan and 2 Division of Pulmonary, Critical Care, and Sleep Medicine, Mount Sinai Medical Center, New York, New York, USA
CORRESPONDENCE: B.A. Rybicki, Dept of Biostatistics and Research Epidemiology, Henry Ford Health System, One Ford Place, 3E, Detroit, MI 48202, USA. Fax: 1 3138746730. E-mail: brybick1@hfhs.org
Keywords: Calcitriol receptor, genetic predisposition, peptidyl-dipeptidase A, polymorphism
Received: January 15, 2004
Accepted April 9, 2004
This study was supported by a National Institutes of Health grant R01 HL54306.
| Abstract |
|---|
|
|
|---|
In the present study, using a family-based study design, which controls for population stratification, the authors attempted to replicate previously reported associations between sarcoidosis and three attractive candidate genes studied primarily in case-control samples.
In 225 nuclear families, ascertained through African Americans with a history of sarcoidosis, no evidence was found for an association between sarcoidosis susceptibility and polymorphisms in the angiotensin converting enzyme, vitamin D receptor and tumour necrosis factor-
genes. Further analyses of chronic and acute disease phenotypes failed to reveal any notable associations. Assuming an underlying inheritance model with an additive allelic effect on disease risk, the current study had
8090% statistical power to detect a 3-fold increased risk associated with the putative risk allele of the polymorphisms under study.
The present authors conclude that in African-Americans, the angiotensin converting enzyme, vitamin D receptor, and tumour necrosis factor-
genes are not significant risk factors for sarcoidosis susceptibility.
Sarcoidosis is a systemic granulomatous disease of unknown aetiology that can affect any organ system, but most frequently includes the lungs (in
90% of cases), skin and eyes. A genetic susceptibility to sarcoidosis has been supported by studies demonstrating familial clustering 1 and a higher prevalence of disease in first-degree relatives 2. Clustering of cases within specific ethnic groups also suggests a genetic susceptibility. In the USA, African Americans have a higher disease incidence and morbidity compared with Caucasians 3. The pathological hallmark of sarcoidosis is multinucleated noncaseating granulomas in affected organs. Granuloma burden is indicative of disease activity and long-term active disease can lead to fibrosis and organ dysfunction. Key genes that play a role in granuloma formation and persistence include angiotensin converting enzyme (ACE), tumour necrosis factor (TNF)-
, and the vitamin D receptor (VDR). The epithelioid cells in sarcoid-granulomas secrete >40 different cytokines and other mediators, including ACE 4. Epithelioid cells of sarcoid granulomas contain ACE, and ACE synthesis in monocyte cell culture can be modulated by T-lymphocytes obtained from sarcoid patients 5. Serum ACE levels in healthy individuals 6 and sarcoidosis patients 7 are influenced by a 287 base pair Alu sequence either present (insertion (I)) or absent (deletion (D)) in intron 16 ofthe ACE gene. Association studies of the ACE I/D polymorphism and sarcoidosis have produced mixed results. The ACE D allele, which is associated with higher serum ACE levels 6, has been found associated with increased sarcoidosis risk in African Americans 8, Japanese women 7 and in sarcoid-affected family members from a German population 9. Conversely, other studies in both Caucasians 10, 11 and Japanese 12 have failed to find an association between the ACE I/D polymorphism and sarcoidosis susceptibility or progression.
TNF-
is a cytokine that probably has a key role in granuloma formation 13. TNF-
is increased at sites of disease activity 14 and anti-TNF-
agents have shown promising therapeutic value in chronic sarcoidosis 15. At least eleven polymorphisms in the promoter region of the TNF-
gene have been reported 16 but of these a G to A single base pair (bp) transition at position 307 is considered the most likely to affect gene transcription 16. Several studies have examined this polymorphism in relation to sarcoidosis susceptibility 1719, with a single case-control study in African Americans having null findings 20.
The active form of vitamin D, 1,25-dihydroxyvitamin D3 (1,25(OH)2D3), a promoter of multinucleated giant cell formation 21, is produced at sites of sarcoid granulomas 22, and is elevated in sarcoidosis patients 23. This hormone acts through the binding of intranuclear VDRs 24 and has been shown to interact directly with interferon gamma to promote the granulomatous process 25. Several polymorphisms with potential functional significance have been identified in the VDR gene 26. The Bsm I and Apa I polymorphisms are located in noncoding regions, whereas the Taq I polymorphism is in an exon but leads to a silent codon change with the nucleotides ATT and ATC both coding for isoleucine at codon 352. The b allele of the Bsm I polymorphism, which has a reported association with sarcoidosis susceptibility 27, is in strong linkage disequilibrium with the absence of a Taq I cleavage site (T allele) in both Caucasians and African-Americans 28. Another DNA sequence variant in the VDR gene is a thymine/cytosine (T/C) polymorphism in the first of two potential start (ATG) codons detected by restriction fragment length polymorphism using the endonuclease Fok I. The functional significance of the Fok I polymorphism is uncertain. However, most data indicate that the F allele, which lacks the first ATG start codon resulting in a transcribed protein three amino acids shorter than the variant f allele, is more effective in transactivation of the 1,25(OH)2D3 signal 29. One study has reported an association of the f allele with tuberculosis 30.
Since conflicting case-control reports of sarcoidosis associations with variants in the ACE, VDR and TNF-
genes exist, the present authors attempted to replicate these associations in a family-based study of sarcoidosis in African-American nuclear families. The strength of the family-based design is that it controls for genetic background unrelated to disease. While case-control samples are generally easier to collect and are statistically more powerful than family-based samples under some circumstances 31, case-control samples can produce spurious associations because of ethnic admixture or other problems relating to inappropriate controls 32. Therefore, the results presented here should help clarify the role of the ACE (I/D), VDR (TaqI and Fok I) and TNF-
(307) genetic polymorphisms with regard to sarcoidosis risk in African Americans.
| Methods |
|---|
|
|
|---|
2 yrs with no other medical condition that could explain the clinical course.
From a total of 623 eligible probands, 359 (58%) were enrolled with one or more first-degree family member(s). Of these, 234 had two or more parents or sibs that donated a blood sample. An additional 10 African-American sarcoidosis families were recruited outside the Henry Ford Health System, resulting in 244 total families. Both parents were genotyped when available. When one or both parents were unavailable for genotyping, all available full sibs were genotyped. Using a panel of unlinked markers, 20 families (8.2%) with Mendelian segregation inconsistencies and two families (0.8%), where the proband, after further follow-up, was not found to have sarcoidosis, were excluded from the analysis. Three of the remaining 222 families had two separate nuclear families, which were analysed separately, resulting in 225 nuclear families or a total of 704 individuals for study. The three basic family configurations are shown in table 1
. While the majority of family members with a sarcoidosis history were probands through whom the family was ascertained, there were also 22 sibs and four additional family members (who were treated as index cases in the analysis) with a history of sarcoidosis based on self-report of a physician diagnosis of this condition. In 18 of the 22 sib cases, and the other four additional family member cases, the diagnosis was confirmed by biopsy. In the remaining four sib cases, a medical record review revealed a chest radiograph report of bilateral hilar adenopathy at the time of diagnosis. Most families (90%) had a single affected offspring and >80% of families had additional unaffected sibs.
|
3 yrs, despite therapy.
Genotyping
High molecular weight DNA was isolated from anticoagulated blood by detergent lysis and organic extraction. The diallelic ACE polymorphism was genotyped by enzymatic amplification of a segment within intron 16 with the 190 and 490 base pair alleles resolved by 6% polyacrylamide gelelectrophoresis. The primers were designed as previously described 8. Amplification was performed using 100 ng of genomic DNA and 10 picomol (pmol) of each primer at 94°C for 1 min, 58°C for 1 min and 72°C for 2 min, for 30 cycles. To confirm that genotypes with only a D band were true homozygotes rather than mistyped ID heterozygotes due to preferential amplification of the D allele, a second reaction was performed using the sense primer labelled with gamma [32P] ATP. Alleles were separated by acrylamide gel electrophoresis, stained with ethidium bromide to display molecular weight markers and autoradiographed to detect alleles.
The VDR Taq I and Fok I polymorphisms were typed by amplification and restriction fragment length polymorphism as previously described 30. Taq I identifies a silent T to C base change at nucleotide 352. Fok I identifies a C to T transition that creates an alternative initiation site three codons upstream from a second ATG start site. Alleles were separated by agarose gel electrophoresis, and detected with ethidium bromide. Each DNA sample was assayed at least twice per marker, to confirm genotype.
The G to A transition at position 307 of the TNF-
gene was determined by amplification of a 413 bp region incorporating two polymorphisms. Amplification was performed using the sense primer 5'-TTCCTGCATCCTGTCTGGAA-3' and antisense primer, 5'-AGGGAGCGTCTGA-TGGCTG-3' at 1 pmol per reaction, 2 mM deoxynucleotidetriphosphates, 10 ng DNA, and 0.62 unit Taq DNA polymerase (Promega, Madison, WI, USA). Amplification conditions were 94°C for 60 s, 62°C for 60 s, 72°C for 60 s, with a final extension step of 6 min at 72°C. Products were sequenced in both directions on a Beckman CEQ 2000XL, by the CEQ Dye Terminator Cycle Sequencing kit, according to the manufacturer's protocol (Beckman Coulter, Fullerton, CA, USA). All genotypes were performed in duplicate with genotypes in disagreement (
10%) run a third time to verify results.
Statistical methodology
To determine whether one or more alleles at the locus of interest was associated with the sarcoidosis phenotype, a Family-Based Association Test (FBAT) statistic, S, calculated using the FBAT software was used 34. In short, this test statistic determines whether a certain allele in a nuclear family is transmitted to affected offspring more (or less) often than would be expected under Mendelian segregation. Test statistics for both dominant and additive inheritance models were computed. The outcome for initial analyses was a dichotomous affection status variable. More in depth phenotype models were also examined in which the three phenotypic levels were treated as either a dichotomous or an ordinal variable. For the former parameterisation, a multivariate model was used where three 0, 1 dummy variables were coded corresponding to absence or presence of the phenotype in question. For the latter parameterisation, affected sibs were coded as 1, 2 or 3 with a higher number corresponding to a more severe phenotypic level. For both types of parameterisations, unaffected sibs were coded as zero.
The statistical power of the study sample for detecting an association between the putative risk allele of each polymorphism and sarcoidosis was calculated with the power calculation of Family-Based Test software 35. Power calculations were conditional on the frequency of each risk allele, an estimated lifetime risk of 1% 3, the study sample configuration depicted in table 1
and a range of genotype relative risks (
s). Using these parameters, expected genotypic-specific penetrance values, f2 (two risk alleles), f1 (one risk allele), and f0 (no risk alleles), under additive (f0=2f1f2), multiplicative (f0=f12/f2) and dominant (f1=f2) inheritance models could be calculated. The routine used to compute power estimates for the family-based analysis was based on a second-order Taylor expansion, which has been shown to have good precision for samples sizes of
100 families 35.
| Results |
|---|
|
|
|---|
|
| (001) |
|
|
|
) ranging 15 were tested at the type I error (
) levels of 5 and 1% (fig. 1
s of
3 were detected at the
=0.05 level with
80% statistical power under an additive inheritance model, and with 90% statistical power under a multiplicative inheritance model. At the
=0.01 level for theadditive and multiplicative models,
s of
4 could be detected with
80% power depending on the inheritance model. Under a dominant inheritance model, the statistical power only approached 60% for a
value of 5 at
=0.05. For the VDR Fok I*f and Taq I*t alleles,
s of
3 were detected at the
=0.05 level with
90% statistical power under all three inheritance models tested. Likewise
s of
3 could be detected with
70% statistical power under all three inheritance models tested at the
=0.01 level. Statistical power for
s of
4 could be detected with
100% power under all three inheritance models tested at both
levels (0.01 and 0.05) tested. Statistical power estimates for the TNF-307*A allele were virtually identical to the two VDR alleles under the dominant and additive inheritance models and slightly less under the multiplicative inheritance model.
|
| Discussion |
|---|
|
|
|---|
1720, 42 genes (table 5
3 under an additive inheritance model, no significant associations were found.
|
17, 42 polymorphisms, then a disproportionate mix of disease types among cases may bias an odds ratio measuring sarcoidosis susceptibility. In the present study population most cases were phenotyped, so the present authors were able to examine the relationship between the four polymorphisms under study and disease severity. No association was found between either the ACE I/D or TNF-307 polymorphisms and disease phenotype, which would suggest that neither of these polymorphisms are related to disease progression in African Americans and that the case mix in the present study population did not influence its susceptibility results. If the TNF-307 polymorphism is actually a marker for Lofgrens syndrome and erythema nodosum in sarcoidosis, as several studies report 17, 42, then these conditions being rare in African-American sarcoidosis cases 43 may explain the negative findings with respect to this polymorphism in the current study. However, the extent to which the current authors could exclude a disease progression gene was limited by a smaller sample size due to missing and incomplete phenotype data, which did not allow the present study to consider the full range of phenotypic expressions the genes under study could exhibit. Genetic and environmental background can also affect case-control genetic association results. All four polymorphisms that were studied are observed in Caucasian, Japanese and African-American populations but their effect on disease susceptibility may be altered by other genetic or environmental factors. If these factors vary among ethnic groups, then ethnic-specific associations of a given polymorphism with sarcoidosis susceptibility are possible.
In the study by Niimi et al. 27, the Bsm I B allele was overrepresented in sarcoidosis patients but Guleva and Seitzer 41 subsequently reported a failure to replicate this finding in a German Caucasian population. In the current study population, the corresponding t allele for the Taq I polymorphism showed no preferential transmission to affected offspring. The current study also examined the Fok I VDR polymorphism, which is in linkage equilibrium with the Taq I/Bsm I allele group 30, and found that this polymorphism was not associated with sarcoidosis susceptibility but the variant f allele may have a modest association with chronic sarcoidosis. The potential functional relevance of the Fok I polymorphism 29 and the association of VDR genotypes with other granulomatous disorders 44, 45 would suggest that the VDR gene warrants further investigation in sarcoidosis. However, in terms of sarcoidosis susceptibility both VDR alleles tested in the present study had an 80% power to detect a 2-fold increased risk under a dominant model, which is the risk level for the heterozygote VDR Bsm I polymorphism previously reported 27.
The extent to which population stratification biases results of case-control genetic association studies is a debated topic. While in theory, selection of a control sample from the same population that cases are drawn should avoid population stratification, in practice uniform case/control sampling is difficult to achieve. Statistical methods to adjust for population stratification using population-specific alleles have been proposed 46 but these methods have yet to be applied widely and may be prone to overcorrection of potential biases. The family-based study design used in the present study provides a sampling procedure that avoids the potential bias of case-control sampling without making any elaborate statistical corrections.
In summary, the present study found no evidence for an association between sarcoidosis susceptibility in African Americans and three key genes, angiotensin converting enzyme, vitamin D receptor and tumour necrosis factor-
, involved in granuloma formation. Overall results from past case-control studies of polymorphisms in these genes are equivocal with regard to sarcoidosis susceptibility. The results of the present study should not be subject to population stratification since potential confounding by genetic background was controlled in the study design. The current study was also sufficiently powered to detect levels of association with disease susceptibility that have been previously reported. The presented results suggest that genetic control of the granulomatous reaction is not related to sarcoidosis susceptibility in African Americans. Alternatively, the angiotensin converting enzyme, vitamin D receptor and tumour necrosis factor-
genes may not be optimal candidates for studying the sarcoidosis granulomatous reaction in relation to disease susceptibility in African Americans.
| Acknowledgements |
|---|
|
|
|---|
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. C. Iannuzzi, B. A. Rybicki, and A. S. Teirstein Sarcoidosis N. Engl. J. Med., November 22, 2007; 357(21): 2153 - 2165. [Full Text] [PDF] |
||||
![]() |
M. C. Iannuzzi Thomas A. Neff Lecture. Advances in the Genetics of Sarcoidosis Proceedings of the ATS, August 15, 2007; 4(5): 457 - 460. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. J. Plant, S. Ghani, M. J. O'Mahony, L. Morgan, C. M. O'Connor, K. Morgan, J. A. Baugh, and S. C. Donnelly Sarcoidosis and MIF gene polymorphism: a case-control study in an Irish population Eur. Respir. J., February 1, 2007; 29(2): 325 - 329. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. C. Iannuzzi and B. A. Rybicki Genetics of Sarcoidosis: Candidate Genes and Genome Scans Proceedings of the ATS, January 1, 2007; 4(1): 108 - 116. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |