A Replication Study of the Association between Rheumatoid Arthritis and Deletion of the Late Cornified Envelope Genes LCE3B and LCE3C

Abstract

Objective

Two recent studies, in a Spanish and a Chinese population, point to an association between rheumatoid arthritis (RA) risk and the deletion of the Late Cornified Envelope (LCE) 3B and 3C genes (LCE3C_LCE3B-del), a known risk factor for psoriasis. We aimed to replicate these studies in a large Dutch cohort.

Methods

1039 RA cases and 759 controls were genotyped for LCE3C_LCE3B-del. Association analysis was performed for the complete cohort and after stratification for the serologic markers anti-cyclic citrullinated peptide and rheumatoid factor. A meta-analysis was performed combining our data with the Spanish and Chinese datasets, resulting in an analysis including 2466 RA cases and 2438 controls.

Results

In the Dutch cohort we did not observe a significant association of LCE3C_LCE3B-del (p = 0.093) with RA risk. A stratified analysis for the serologic positive and negative group did not show an association between the genetic variant and disease risk, either. The meta-analysis, however, confirmed a significant association (p<0.0001, OR = 1.31, 95% confidence interval 1.16–1.47).

Conclusion

Our meta-analysis confirms the association of the LCE3 deletion with RA, suggesting that LCE3C_LCE3B-del is a common risk factor for (auto)immune diseases.

About the Authors:

Judith G. M. Bergboer

Affiliation: Department of Dermatology, Radboud University Nijmegen Medical Centre, Nijmegen Centre for Molecular Life Sciences, Nijmegen, The Netherlands

Maša Umićević-Mirkov

Affiliation: Department of Human Genetics, Radboud University Nijmegen Medical Centre, Nijmegen, The Netherlands

Jaap Fransen

Affiliation: Department of Rheumatology, Radboud University Nijmegen Medical Centre, Nijmegen, The Netherlands

Martin den Heijer

Affiliation: Department of Endocrinology and Department of Epidemiology and Biostatistics, Radboud University Nijmegen Medical Centre, Nijmegen, The Netherlands

Barbara Franke

Affiliation: Department of Human Genetics, Radboud University Nijmegen Medical Centre, Nijmegen, The Netherlands

Piet L. C. M. van Riel

Affiliation: Department of Rheumatology, Radboud University Nijmegen Medical Centre, Nijmegen, The Netherlands

Joost Schalkwijk

Affiliation: Department of Dermatology, Radboud University Nijmegen Medical Centre, Nijmegen Centre for Molecular Life Sciences, Nijmegen, The Netherlands

Marieke J. H. Coenen

* E-mail: [email protected]

Affiliation: Department of Human Genetics, Radboud University Nijmegen Medical Centre, Nijmegen, The Netherlands

on behalf of the Nijmegen Biomedical Study principal investigators

¶A list of principal investigators of the Nijmegen Biomedical Study is provided in the Acknowledgments.

Introduction

Rheumatoid arthritis (RA) is a chronic inflammatory joint disease affecting approximately 1% of the population. Although the causes of the disease are largely unknown, both genetic and environmental factors seem to contribute. Thus far, genome-wide association studies have revealed more than 30 risk loci explaining 16% of the disease variance.

The copy number of DNA segments in the human genome varies in normal individuals and most likely this diversity is responsible for a significant proportion of phenotypic variation. Not much is known about the role of copy number variation (CNV) in the etiology of RA, but relevance of these common genetic variants is expected. A small number of candidate gene studies have indicated that CNVs in Chemokine ligand 3-like 1 (CCL3L1), Fc gamma receptor 3B (FCGR3B) and pre-B lymphocyte 1 (VPREB1) are associated with RA disease susceptibility [1], [2], [3]. Using a genome-wide approach the Welcome Trust Case Control Consortium (WTCCC) reported that CNVs in the HLA region are significantly associated with the development of RA [4]. More recently using the WTCCC dataset, 11 rare copy number variable regions associated with RA were identified [5]. These regions included genes that could be directly linked to the development of RA. Another potential candidate CNV is the deletion of two genes of the Late Cornified Envelope gene cluster, LCE3B and LCE3C (LCE3C_LCE3B-del), which was recently found associated with psoriasis in several populations [6], as was rs4112788, a single nucleotide polymorphism (SNP) in high linkage disequilibrium with the deletion [7]. The observation that the risk of RA is increased in the offspring of individuals affected with other autoimmune disorders, including psoriasis, prompted Docampo and colleagues [8] to investigate the association between the LCE3C_LCE3B-del and RA in two independent Spanish cohorts. In their study an association was found between RA and homozygosity for LCE3C_LCE3B-del. Recently it has been demonstrated that the deletion is also associated with RA susceptibility in a non-Caucasian (Chinese-Han) population [9]. However, they were unable to confirm the Spanish results that the association between LCE3C_LCE3B-del and RA risk was restricted to anti-CCP-positive patients. We sought to replicate these studies in a similarly powered Dutch cohort. A total of 1039 RA cases and 759 controls were genotyped for LCE3C_LCE3B-del, and we investigated whether this genetic variant is associated with RA. Additionally, we combined the Spanish, Chinese and Dutch data in a meta-analysis.

Results

An overview of the demographics of the patients and controls used in this study are presented in Table 1. Table 2 provides the results of the association analysis for the LCE3C_LCE3B-del and RA risk in the Dutch cohort. There were 1039 RA patients and 779 controls successfully genotyped for LCE3C_LCE3B-del. The LCE3C_LCE3B-del and the linked SNP rs4112788 were in near-perfect linkage disequilibrium (LD) with each other (r² = 0.928 and D′ = 0.984), confirming that LCE3C_LCE3B-del was genotyped accurately.

[Figure omitted. See PDF.]

Table 1. Demographics of the rheumatoid arthritis cohort.

https://doi.org/10.1371/journal.pone.0032045.t001

[Figure omitted. See PDF.]

Table 2. No association of the LCE3C_LCE3B deletion with RA in the Dutch cohort.

https://doi.org/10.1371/journal.pone.0032045.t002

In the Dutch cohort, we did not find evidence for an association between LCE3C_LCE3B-del and RA, at the genotypic (p = 0.093, OR 1.18, 95% confidence interval (95%CI): 0.97–1.43 (Table 2)) nor the allelic level (p = 0.147, OR 1.11 95%CI: 0.97–1.27 (data not shown)). It has been hypothesized that the LCE3C_LCE3B-del might result in a defective skin barrier repair leading to an abnormal or increased exposure to environmental antigens. This could lead to an earlier age of disease onset in carriers of LCE3C_LCE3B-del compared to non-carriers. An analysis for age of onset, however, did not show association with the LCE3C_LCE3B-del status (p = 0.25, t-test).

Docampo et al. [8] demonstrated that LCE3C_LCE3B-del showed the strongest association with RA in anti-CCP and rheumatoid factor (RF) positive patients. In the study of Lu et al. [9] no associations were found in the different RA subsets, these analyses were corrected for age and gender. We performed association analyses of the anti-CCP and RF patient subsets separately, for those patients with information available. We did not find an association between LCE3C_LCE3B-del and RA in any of the subsets analyzed (Table 2).

We performed an overall meta-analysis, using a recessive logistic regression model unadjusted for age and gender, on the combined data from this study and the studies of Docampo et al. [8] and Lu et al. [9] (Table 3). Despite the lack of significance of the association in the Dutch cohort, in the meta-analysis a significant association between LCE3C_LCE3B-del and RA was found (p<0.0001, OR 1.31 95%CI: 1.16–1.47). Next to the overall analysis, we performed also a meta-analysis on the anti-CCP and RF subgroups combining the Spanish, Chinese and our data (Table 3), unfortunately not for all patients these data were available. In this subgroup analysis all previously found associations [8] between LCE3C_LCE3B-del and RA were significant, with the association apparently strongest in RF positive subset of patients (p = 0.0007, OR 1.27, 95%CI: 1.11–1.45) (Table 3).

[Figure omitted. See PDF.]

Table 3. Meta-analyses of association of the LCE3C_LCE3B deletion with RA in the Spanish, Chinese and Dutch cohort.

https://doi.org/10.1371/journal.pone.0032045.t003

Discussion

Although a previously reported association of the LCE3C_LCE3B-del with RA risk did not reach significance in our Dutch cohort, a meta-analysis combining our data and the previously reported Spanish and Chinese data [8], [9] confirmed the association of LCE3C_LCE3B-del with RA.

Our study had a power of over 95% to find an association, although we based this calculation on the previously found OR of 1.45. When we recalculate the power using the OR of the meta-analysis (OR = 1.31), the power of our study drops to 73%. A study population of 1256 samples will be needed to reach 80% power. The ORs observed in the Dutch population are in the same direction as seen in the previous studies [8], [9] and our meta-analysis showed a statistically significant association of the deletion with RA. This suggests that a larger population is necessary to prove whether the deletion is associated with RA. The contradictory findings from this study and the previous RA studies [8], [9] could also be caused by differences in the study population, e.g. disease phenotype or ethnicity. In addition, there is a discrepancy between the allele frequency of the deletion in the controls from Spain (55%) and China (54%) compared to The Netherlands (61%). This heterogeneity of LCE3C_LCE3B-del allele frequency in control groups from different ethnic backgrounds was already known from previous studies [6], [10]. Indeed, a recent meta-analysis confirming an association between LCE3C_LCE3B-del and psoriatic arthritis showed that the OR of the Spanish population was much higher than observed in an Italian population (1.66 versus 1.23) [10].

In our meta-analysis we used a recessive model, because this model fitted the LCE data in the psoriasis study best, based on the Akaike Information Criteria [11] which is in line with the two previous studies on LCE3 in patients with RA [8], [9]. Although we could not identify an association between the deletion and RA in the Dutch cohort, addition of the data to the other studies resulted in lower p-values than without this dataset, showing a significant contribution of the Dutch population to the overall results. The addition of our Dutch samples leads to a lower OR and a smaller CI than previously reported. The other two studies may have significantly overestimated the magnitude of the odds ratio, the so-called winner's curse [12]. Due to publication bias it might be possible that some studies investigating the same subject have not been published, this information will be important to confirm the association between RA and the deletion as shown in this study.

In our subgroup analysis, we observed the strongest association in the RF positive patients. Larger datasets will be necessary to clarify whether the observed association is patient subset specific. In contrast to the results from our meta-analysis, genome-wide CNV analyses using the WTCCC RA dataset did not identify association with the LCE3 region [4], [5]. One study only assessed rare CNVs (population frequency <5%) which excludes the LCE3C_LCE3B-del [5]. The other study included a probe covering the region that passed QC in their analysis (probe CNVR358.1, chr1:150,822,234–150,856,715 (UCSC genome browser hg18) [4]). It might be possible that the region is associated with RA in the WTCCC dataset at lower significance level than their genome-wide significance threshold.

Until now, LCE3C_LCE3B-del was found to be associated with psoriasis, psoriatic arthritis and RA in several populations [6], [7], [8], [9], [10], [11], [13], [14], whereas no association was found for atopic dermatitis [15]. Also in one study an association with systemic lupus erythematosus was found [9], implying that LCE3C_LCE3B-del may be a common risk factor for (auto)immune diseases. The function of the LCE genes has only been studied in skin in general [16], [17] and in relation to psoriasis [18]. From these studies it is known that the LCE proteins are likely to be incorporated in the cornified cell envelope, which is an important structure in the barrier function of skin. When comparing normal and psoriasis skin, the genes of the LCE1, 2, 5 and 6 groups are mainly expressed in normal skin, whereas the LCE3 genes are predominantly expressed in psoriasis skin. Moreover, upon barrier disruption of normal skin the LCE1, 2, 5 and 6 gene-groups are downregulated, while the expression of the LCE3 genes is upregulated. Altogether, these data imply a role in barrier repair for the LCE3 proteins and a role in barrier maintenance for the other LCEs. As hypothesized by Docampo et al. [8], the absence of LCE3B and LCE3C could compromise barrier function of the epithelia and possibly facilitate the entrance of environmental antigens or pathogens. Since the LCE genes are mainly expressed in skin and oral epithelia, these would be relevant tissues that could facilitate the entrance of antigens or pathogens, like the Epstein-Barr virus and the cytomegalovirus, thereby triggering RA. We reasoned that if this would be the case patients with LCE3C_LCE3B-del might be exposed more readily to common antigens and therefore would have an earlier age of onset due to easy access of pathogens/triggers through the skin. However, we did not observe such an association between the genotype and age of RA onset.

Our meta-analysis showed that the LCE3B and LCE3C deletion is associated with RA, though the contribution of our large Dutch sample is small. Therefore it will be necessary to test even larger patient cohorts to shed more light on the possible association of the deletion in specific RA patient subsets. An interesting next step would be to perform functional studies to unravel the mechanisms underlying this association.

Materials and Methods

Ethics statement

The ethical committee of the Radboud University Nijmegen Medical Centre (“Commissie Mensgebonden Onderzoek (CMO) Regio Arnhem Nijmegen”) approved the study (CMO number 2004/014); all patients and controls gave written informed consent.

DNA samples

All patients met the American College of Rheumatology 1987 revised criteria for RA. The patients were part of the Dutch Rheumatoid Arthritis Monitoring registry (www.dreamregistry.nl) (n = 1039) [19], [20]. The controls (n = 779) were participating in the Nijmegen Biomedical Study [21]. For cohort characteristics see Table 1.

Genotyping

DNA was extracted from blood using salt extraction. For the LCE3C_LCE3B-del genotyping a polymerase chain reaction was performed as earlier described [11]. 112 samples were genotyped in duplicate; error rate was 2.9% and assay failure rate was 0.3%. To validate the LCE3C-LCE3B-del genotype, we also genotyped the linked SNP rs4112788, which was performed using a Taqman® SNP genotyping assay (assay ID C__31910050_10) according to manufacturer's recommendations (Applied Biosystems, Nieuwerkerk aan den IJssel, The Netherlands). Five percent of the samples were analyzed in duplicate; all genotypes were concordant. The assay failure rate was 2.6%.

Statistical analysis

Descriptive statistics for quantitative values are given as mean ± standard deviation (SD). Both variants showed no deviations from Hardy-Weinberg equilibrium for either patients or controls. Co-dominant and recessive logistic regression models were used to assess the genetic effect of LCE3C_LCE3B-del on RA risk, using homozygosity for the allele without deletion as a reference category. These models were similar to those used in the previous studies [8], [9]. The statistical analyses were performed using SPSS software 16.0 (SPSS Inc., Chicago, IL, U.S.A.). A meta-analysis was performed in Review Manager 5 (Review Manager (RevMan) Version 5.0. Copenhagen: The Nordic Cochrane Centre, The Cochrane Collaboration, 2008) using the aggregated data from this studies, Docampo et al. [8] and Lu et al. [9]. Heterogeneity of odds ratios (ORs) among cohorts was calculated using the Breslow-Day method, and pooled ORs were calculated under a fixed effects model (Mantel-Haenszel meta-analysis). A power calculation was performed using the Genetic Power Calculator [22] with input values derived from Docampo et al. [8] (significance level 0.05, OR = 1.45, allele frequency 0.62, recessive model).

Acknowledgments

We would like to thank all the volunteers who participated in this study. Principal investigators of the Nijmegen Biomedical Study are L.A.L.M. Kiemeney, M. den Heijer, A.L.M. Verbeek, D.W. Swinkels en B. Franke.

Author Contributions

Conceived and designed the experiments: JB MU JF JS MC. Performed the experiments: JB MU. Analyzed the data: JB MU JF MC. Contributed reagents/materials/analysis tools: MH BF PR JS MC. Wrote the paper: JB MU MC JS. Revised the paper: JF MH BF PR. Final approval: JB MU JF MH BF PR JS MC.

Citation: Bergboer JGM, Umićević-Mirkov M, Fransen J, den Heijer M, Franke B, van Riel PLCM, et al. (2012) A Replication Study of the Association between Rheumatoid Arthritis and Deletion of the Late Cornified Envelope Genes LCE3B and LCE3C. PLoS ONE 7(2): e32045. https://doi.org/10.1371/journal.pone.0032045

References

1. McKinney C, Merriman ME, Chapman PT, Gow PJ, Harrison AA, et al. (2008) Evidence for an influence of chemokine ligand 3-like 1 (CCL3L1) gene copy number on susceptibility to rheumatoid arthritis. Ann Rheum Dis 67: 409–413.C. McKinneyME MerrimanPT ChapmanPJ GowAA Harrison2008Evidence for an influence of chemokine ligand 3-like 1 (CCL3L1) gene copy number on susceptibility to rheumatoid arthritis.Ann Rheum Dis67409413

2. McKinney C, Fanciulli M, Merriman ME, Phipps-Green A, Alizadeh BZ, et al. (2010) Association of variation in Fcgamma receptor 3B gene copy number with rheumatoid arthritis in Caucasian samples. Ann Rheum Dis 69: 1711–1716.C. McKinneyM. FanciulliME MerrimanA. Phipps-GreenBZ Alizadeh2010Association of variation in Fcgamma receptor 3B gene copy number with rheumatoid arthritis in Caucasian samples.Ann Rheum Dis6917111716

3. Yim SH, Chung YJ, Jin EH, Shim SC, Kim JY, et al. (2011) The potential role of VPREB1 gene copy number variation in susceptibility to rheumatoid arthritis. Mol Immunol 48: 1338–1343.SH YimYJ ChungEH JinSC ShimJY Kim2011The potential role of VPREB1 gene copy number variation in susceptibility to rheumatoid arthritis.Mol Immunol4813381343

4. Craddock N, Hurles ME, Cardin N, Pearson RD, Plagnol V, et al. (2010) Genome-wide association study of CNVs in 16,000 cases of eight common diseases and 3,000 shared controls. Nature 464: 713–720.N. CraddockME HurlesN. CardinRD PearsonV. Plagnol2010Genome-wide association study of CNVs in 16,000 cases of eight common diseases and 3,000 shared controls.Nature464713720

5. Uddin M, Sturge M, Rahman P, Woods MO (2011) Autosome-wide Copy Number Variation Association Analysis for Rheumatoid Arthritis Using the WTCCC High-density SNP Genotype Data. J Rheumatol 38: 797–801.M. UddinM. SturgeP. RahmanMO Woods2011Autosome-wide Copy Number Variation Association Analysis for Rheumatoid Arthritis Using the WTCCC High-density SNP Genotype Data.J Rheumatol38797801

6. Riveira-Munoz E, He SM, Escaramis G, Stuart PE, Huffmeier U, et al. (2011) Meta-Analysis Confirms the LCE3C_LCE3B Deletion as a Risk Factor for Psoriasis in Several Ethnic Groups and Finds Interaction with HLA-Cw6. J Invest Dermatol 131: 1105–1109.E. Riveira-MunozSM HeG. EscaramisPE StuartU. Huffmeier2011Meta-Analysis Confirms the LCE3C_LCE3B Deletion as a Risk Factor for Psoriasis in Several Ethnic Groups and Finds Interaction with HLA-Cw6.J Invest Dermatol13111051109

7. Zhang XJ, Huang W, Yang S, Sun LD, Zhang FY, et al. (2009) Psoriasis genome-wide association study identifies susceptibility variants within LCE gene cluster at 1q21. Nat Genet 41: 205–210.XJ ZhangW. HuangS. YangLD SunFY Zhang2009Psoriasis genome-wide association study identifies susceptibility variants within LCE gene cluster at 1q21.Nat Genet41205210

8. Docampo E, Rabionet R, Riveira-Munoz E, Escaramis G, Julia A, et al. (2010) Deletion of the late cornified envelope genes, LCE3C and LCE3B, is associated with rheumatoid arthritis. Arthritis Rheum 62: 1246–1251.E. DocampoR. RabionetE. Riveira-MunozG. EscaramisA. Julia2010Deletion of the late cornified envelope genes, LCE3C and LCE3B, is associated with rheumatoid arthritis.Arthritis Rheum6212461251

9. Lu X, Guo J, Zhou X, Li R, Liu X, et al. (2011) Deletion of LCE3C_LCE3B is associated with rheumatoid arthritis and systemic lupus erythematosus in the Chinese Han population. Ann Rheum Dis 70: 1648–1651.X. LuJ. GuoX. ZhouR. LiX. Liu2011Deletion of LCE3C_LCE3B is associated with rheumatoid arthritis and systemic lupus erythematosus in the Chinese Han population.Ann Rheum Dis7016481651

10. Docampo E, Giardina E, Riveira-Munoz E, de Cid R, Escaramis G, et al. (2011) Deletion of the LCE3C and LCE3Bgenes is a susceptibility factor for psoriatic arthritis: A study in Spanish and Italian populations and meta-analysis. Arthritis Rheum 63: 1860–1865.E. DocampoE. GiardinaE. Riveira-MunozR. de CidG. Escaramis2011Deletion of the LCE3C and LCE3Bgenes is a susceptibility factor for psoriatic arthritis: A study in Spanish and Italian populations and meta-analysis.Arthritis Rheum6318601865

11. de Cid R, Riveira-Munoz E, Zeeuwen PL, Robarge J, Liao W, et al. (2009) Deletion of the late cornified envelope LCE3B and LCE3C genes as a susceptibility factor for psoriasis. Nat Genet 41: 211–215.R. de CidE. Riveira-MunozPL ZeeuwenJ. RobargeW. Liao2009Deletion of the late cornified envelope LCE3B and LCE3C genes as a susceptibility factor for psoriasis.Nat Genet41211215

12. Young NS, Ioannidis JP, Al-Ubaydli O (2008) Why current publication practices may distort science. PLoS Med 5: e201.NS YoungJP IoannidisO. Al-Ubaydli2008Why current publication practices may distort science.PLoS Med5e201

13. Huffmeier U, Bergboer JG, Becker T, Armour JA, Traupe H, et al. (2010) Replication of LCE3C-LCE3B CNV as a risk factor for psoriasis and analysis of interaction with other genetic risk factors. J Invest Dermatol 130: 979–984.U. HuffmeierJG BergboerT. BeckerJA ArmourH. Traupe2010Replication of LCE3C-LCE3B CNV as a risk factor for psoriasis and analysis of interaction with other genetic risk factors.J Invest Dermatol130979984

14. Bowes J, Flynn E, Ho P, Aly B, Morgan AW, et al. (2010) Variants in linkage disequilibrium with the late cornified envelope gene cluster deletion are associated with susceptibility to psoriatic arthritis. Ann Rheum Dis 69: 2199–2203.J. BowesE. FlynnP. HoB. AlyAW Morgan2010Variants in linkage disequilibrium with the late cornified envelope gene cluster deletion are associated with susceptibility to psoriatic arthritis.Ann Rheum Dis6921992203

15. Bergboer JG, Zeeuwen PL, Irvine AD, Weidinger S, Giardina E, et al. (2010) Deletion of Late Cornified Envelope 3B and 3C genes is not associated with atopic dermatitis. J Invest Dermatol 130: 2057–2061.JG BergboerPL ZeeuwenAD IrvineS. WeidingerE. Giardina2010Deletion of Late Cornified Envelope 3B and 3C genes is not associated with atopic dermatitis.J Invest Dermatol13020572061

16. Jackson B, Tilli CM, Hardman MJ, Avilion AA, MacLeod MC, et al. (2005) Late cornified envelope family in differentiating epithelia–response to calcium and ultraviolet irradiation. J Invest Dermatol 124: 1062–1070.B. JacksonCM TilliMJ HardmanAA AvilionMC MacLeod2005Late cornified envelope family in differentiating epithelia–response to calcium and ultraviolet irradiation.J Invest Dermatol12410621070

17. Marshall D, Hardman MJ, Nield KM, Byrne C (2001) Differentially expressed late constituents of the epidermal cornified envelope. Proc Natl Acad Sci U S A 98: 13031–13036.D. MarshallMJ HardmanKM NieldC. Byrne2001Differentially expressed late constituents of the epidermal cornified envelope.Proc Natl Acad Sci U S A981303113036

18. Bergboer JG, Tjabringa GS, Kamsteeg M, van Vlijmen-Willems IM, Rodijk-Olthuis D, et al. (2011) Psoriasis risk genes of the late cornified envelope-3 group are distinctly expressed compared with genes of other LCE groups. Am J Pathol 178: 1470–1477.JG BergboerGS TjabringaM. KamsteegIM van Vlijmen-WillemsD. Rodijk-Olthuis2011Psoriasis risk genes of the late cornified envelope-3 group are distinctly expressed compared with genes of other LCE groups.Am J Pathol17814701477

19. Toonen EJ, Coenen MJ, Kievit W, Fransen J, Eijsbouts AM, et al. (2008) The tumour necrosis factor receptor superfamily member 1b 676T>G polymorphism in relation to response to infliximab and adalimumab treatment and disease severity in rheumatoid arthritis. Ann Rheum Dis 67: 1174–1177.EJ ToonenMJ CoenenW. KievitJ. FransenAM Eijsbouts2008The tumour necrosis factor receptor superfamily member 1b 676T>G polymorphism in relation to response to infliximab and adalimumab treatment and disease severity in rheumatoid arthritis.Ann Rheum Dis6711741177

20. Wagener FA, Toonen EJ, Wigman L, Fransen J, Creemers MC, et al. (2008) HMOX1 promoter polymorphism modulates the relationship between disease activity and joint damage in rheumatoid arthritis. Arthritis Rheum 58: 3388–3393.FA WagenerEJ ToonenL. WigmanJ. FransenMC Creemers2008HMOX1 promoter polymorphism modulates the relationship between disease activity and joint damage in rheumatoid arthritis.Arthritis Rheum5833883393

21. Hoogendoorn EH, Hermus AR, de Vegt F, Ross HA, Verbeek AL, et al. (2006) Thyroid function and prevalence of anti-thyroperoxidase antibodies in a population with borderline sufficient iodine intake: influences of age and sex. Clin Chem 52: 104–111.EH HoogendoornAR HermusF. de VegtHA RossAL Verbeek2006Thyroid function and prevalence of anti-thyroperoxidase antibodies in a population with borderline sufficient iodine intake: influences of age and sex.Clin Chem52104111

22. Purcell S, Cherny SS, Sham PC (2003) Genetic Power Calculator: design of linkage and association genetic mapping studies of complex traits. Bioinformatics 19: 149–150.S. PurcellSS ChernyPC Sham2003Genetic Power Calculator: design of linkage and association genetic mapping studies of complex traits.Bioinformatics19149150