UGT2B17 copy number gain in a large ankylosing spondylitis multiplex family
© Uddin et al.; licensee BioMed Central Ltd. 2013
Received: 17 September 2012
Accepted: 3 August 2013
Published: 8 August 2013
The primary objective of this study is to identify novel copy number variations (CNVs) associated with familial ankylosing spondylitis (AS). A customized genome-wide microarray was designed to detect CNVs and applied to a multiplex AS family with six (6) affected family members. CNVs were detected using the built-in DNA analytics aberration detection method-2 (ADM-2) algorithm. Gene enrichment analysis was performed to observe the segregation. Subsequent validation was performed using real time quantitative fluorescence polymerase reaction (QF-PCR). The frequency of copy number variation for the UGT2B17 gene was then performed on two well-defined AS cohorts. Fisher exact test was performed to quantify the association.
Our family-based analysis revealed ten gene-enriched CNVs that segregate with all six family members affected with AS. Based on the proposed function and the polymorphic nature of the UGT2B17 gene, the UGT2B17 gene CNV was selected for validation using real time QF-PCR with full concordance. The frequency of two copies of the UGT2B17 gene CNV was 0.41 in the Newfoundland AS cases and 0.35 in the Newfoundland controls (OR = 1.26(0.99-1.59); p < 0.05)), whereas the frequency of two (2) copies of the UGT2B17 gene CNV was 0.40 in the Alberta AS cases and 0.39 in the Alberta controls (OR = 1.05(95% CI: 0.83-1.33); p < 0.71)).
A genome-wide microarray interrogation of a large multiplex AS family revealed segregation of the UGT2B17 gene CNV among all affected family members. The association of the UGT2B17 CNV with AS is particularly interesting given the recent association of this CNV with osteoporosis and the proposed function as it encodes a key enzyme that inhibits androgens. However, two copies of the UGT2B17 gene CNV were only marginally significant in a uniplex AS cohort from Newfoundland but not in a uniplex AS cohort from Alberta.
KeywordsCopy number variants (CNV) UGT2B17 Ankylosing spondylitis Complex autoimmune disease
Ankylosing spondylitis (AS) is the prototypic inflammatory spondyloarthritis with a peak onset between 20 to 30 years [1, 2]. Extra-articular features of AS include inflammation of the eyes, skin, bowels, and more rarely the lungs and heart. There is a male preponderance as men are affected 2 to 3 times more frequently than females [3, 4].
Genetic factors are of major importance in susceptibility to AS and, in fact, genetic epidemiological studies suggest that AS represents one of the most heritable complex autoimmune diseases with a heritability greater than 90% and a sibling recurrence ratio of at least 52 [2–4]. HLA-B*27, which was first recognized to be associated with AS in 1973, remains the strongest genetic association signal with AS, and it is estimated that HLA-B*27 accounts for 23% of the genetic heritability [2–4]. Despite this strong association, only a small fraction (1-5%) of HLA-B*27 positive individuals develop AS [2–4]. Recent genome-wide association studies (GWAS) in several European populations and in the Han Chinese population have identified up to 15 high priority genes including IL23R, RUNX3, KIF21B, 2p15, IL1R2, PTGER4, ERAP1, IL12B, CARD9, TNFR1/ LTBR, TBKBP1 . These associations were primarily reported based on single nucleotide polymorphisms (SNP) analysis. Recently, an association study revealed the complex interactions between the non-HLA gene, ERAP1, with the HLA-B*27 gene [5, 6]. However, the genetic risk described by HLA or non-HLA genes suggests that other genomic variants may contribute to the risk factor for AS.
The premise that CNVs contribute to disease pathogenesis is supported by their capacity to disrupt gene expression and to interrupt functional pathways . In recent years, many CNVs have been associated with complex diseases and the majority of these associations were in autoimmune-mediated diseases, including rheumatoid arthritis, psoriasis, Crohn’s disease, and systemic lupus erythematous [8–11]. Here, we report the first validated CNV associated specifically with AS using a custom genome-wide microarray in a well-defined multiplex AS family.
Disease prevalence (case control studies)
In the case–control cohorts, real time QF-PCR analysis targeting the identical CNVR (from family analysis) revealed the presence of three (3) alleles (0 copy - homozygous deletion, 1 copy - heterozygous deletion, and 2 copies) within this population for the UGT2B17 gene.
UGT2B17 gene CNV association in ankylosing spondylitis
Copy number variations are increasingly being recognized in complex autoimmune diseases as they are capable of altering gene dosage and consequently affect gene function . The contribution of CNVs in AS pathogenesis remains to be systematically evaluated. The primary goal of this study was to identify functionally relevant CNVs segregating within a multiplex AS family and then subsequently determine the allele frequencies of CNVs in multiple case–control cohorts. The genome-wide CNV analysis performed in this study revealed that increased copy number of the UGT2B17 gene is a potential risk factor for AS.
Given that UGT2B17 is a functional gene encoding an enzyme that metabolizes steroid hormones including testosterone and selected xenobiotics [13–16], it represents an excellent candidate gene for susceptibility to AS. A gene-dosage effect is present as this CNV is associated with urine testosterone level, male insulin sensitivity, fat mass, and prostate-cancer risk (as summarized by Xue) . Importantly, copies of the UGT2B17 gene have recently been associated with osteoporosis (including hip fracture) . The dose of the UGT2B17 gene was significantly associated with bone mineral density, cortical thickness, and buckling ratio. These results support the role of UGT2B17 CNVs in the pathogenesis of osteoporosis. Although AS is characterized by new bone formation, osteoporosis is a well recognized complication of AS and the risk of clinical vertebral fractures is increased in AS patients [17, 18]. This complexity with respect to genotype-phenotype correlation requires further investigation.
From an evolutionary perspective, the UGT2B17 gene in humans is highly stratified which renders it more likely to be associated with disease. The deletion allele (i.e., 0 copies) frequency varies significantly among populations. For example, the frequency of the deletion allele is approximately 0.15 in Europeans, ~0.22 in Africans, and ~0.80-0.90 in Asians [7, 13]. In contrast, the frequency of two (2) copies of the UGT2B17 gene is rare in Asian populations (~0.01-0.03), ~0.30-0.40 (reports with varying frequencies) in Europeans, and ~0.60 in African populations.
In this study, two copies of the UGT2B17 gene clearly segregates within a large well-defined AS multiplex family. Given the nature and proposed function of the UGT2B17 gene [19, 20], the UGT2B17 CNV may contribute to susceptibility for AS within this family. Conflicting results were obtained regarding the excess prevalence of this CNV in large uniplex AS cohorts. A trend was observed in the cohort from Newfoundland. However, this was not replicated in the Alberta cohort. In this study, the deletion (i.e., 0 copies) frequency of the UGT2B17 gene in the Newfoundland population is consistent with the Alberta population. There was a noticeable frequency difference between the Alberta (0.51) and Newfoundland (0.46) populations for a single copy of the UGT2B17 gene. This stratification may represent a contributing factor for the lack of association.
The analysis of CNVs within multiplex families provides a unique opportunity to detect variants that may be family-specific. There are multiple reports in other complex diseases where family specific CNVs manifest with the disease . Although it is a difficult task to obtain more multi-generational pedigrees to follow up the segregation of this CNV, the findings observed within the investigated family suggests that segregation analysis into multiple multi-generation pedigrees is promising. Thus, it is important to assess other large multiplex families to determine if this CNV may be a cause of familial AS.
This is study applied a custom designed genome-wide microarray to investigate a large multiplex AS family revealed segregation of the UGT2B17 gene CNV among all affected family members. The association of the UGT2B17 CNV with AS is particularly interesting given the recent association of this CNV with osteoporosis and the proposed function as it encodes a key enzyme that inhibits androgens. This result provides a basis for further exploration on the role of CNV in AS pathogenesis. Further replication is required for the marginal significance of UGT2B17 association identified with AS cohort from Newfoundland.
Multiplex AS family
A multiplex AS family from Newfoundland was identified for this study. All members of the multiplex family were “native Newfoundlanders” of North European ancestry. Ten (10) members from this large family participated in the study. Each individual was assessed clinically including imaging studies. Six (6) of the ten (10) family members had AS, as defined by the 1984 modified New York criteria, while the remaining family members were not affected. Informed consent for participation in the study was obtained from participants or, where participants are children, a parent or guardian. After consent was obtained, all participating family members had their DNA collected. All samples involving human subjects reported in this manuscript have been approved by the Memorial University of Newfoundland Human Investigation Committee.
DNA was obtained from two (2) well-defined AS cohorts from the Newfoundland and Alberta populations. All cases and controls were of North European ancestry. We assessed 298 AS cases satisfying the modified New York criteria and 299 ethnically-matched controls from a homogenous population of Newfoundland. The second cohort consisted of 289 AS cases satisfying the modified New York criteria and 285 ethnically-matched controls from Alberta.
We designed a customized genome-wide microarray based on genomic hotspot breakpoints previously identified . The targeted regions are prone to produce CNVs due to their structural architecture (i.e. presence of segmental duplication, dose sensitive genes etc.). The custom design tiling microarray consisted of 2 × 1 million probes covering the genome with a mean spacing of 280 bp. Prior to CNV analysis, QC measures were applied and the derivative of log ratio spread (DLRS) <0.25 was considered as threshold. CNVs were detected using the built-in Aberration Detection Method-2 (ADM-2) algorithm DNA Analytics v.4.0.85 (Agilent Technologies) using the following criteria: 1) at least 5 probes for a CNV call on GC-corrected intensity; 2) nested filter was set to 2; and 3) log intensity >0.25 for duplication and <−0.25 for deletion. A custom script was applied to detect gene-enriched CNVs (i.e., overlaps with or constitutes a gene) that segregate (at least three cases) within affected AS family members and absent in the unaffected family members.
To validate the CNV encompassing the UGT2B17 (UDP glucuronosyltransferase 2 family, polypeptide B17) gene in the family members and in the case–control cohorts, copy number estimation of the UGT2B17 gene was performed using the Taqman copy number assay Hs03185327_cn (Life Technologies) using the manufacturer’s recommended protocol. The assay was performed in quadruplicate using 10ng genomic DNA for each sample in a 96-well plate. The 10 μl reaction mix consisted of 2 μl of 2× Taqman Genotyping Master Mix (Life Technologies), 0.5 μl of 20× copy number assay (described above), 0.5 μl of TaqMan RNAse P Copy Number Reference Assay (Life Technologies, part 4403326), 2 μl of water, and 2 μl of 5 ng/μl genomic DNA. Cycling conditions for the reaction were 95°C for 10 min, followed by 40 cycles of 95°C for 15 sec and 60°C for 1 min. Samples were analyzed using the ViiA™ 7 Real-Time PCR System (Life Technologies) and analyzed using CopyCaller Software (Life Technologies, PN 4412907). Three reference (calibrator) DNA HapMap samples (NA10851, NA15510 and NA07048 (Coriell Institute)) plus one non-template control were included with the test samples. Fisher exact test was performed to determine association of the UGT2B17 gene copy number with AS in our case–control data.
We would like to thank Genome Quebec (GQ) for helping us on microarray array analysis and The Centre for Applied Genomics (TCAG) for QF-PCR validation. This work was funded by the Canadian Arthritis Society.
- Brown MA, Laval SH, Brophy S, Calin A: A Recurrence risk modelling of the genetic susceptibility to ankylosing spondylitis. Ann. Rheum. Dis. 2000, 59: 883-886. 10.1136/ard.59.11.883.PubMed CentralView ArticlePubMedGoogle Scholar
- Brown MA, Kennedy LG, MacGregor AJ, Darke C, Duncan E, Shatford JL, Taylor A, Calin A, Wordsworth P: Susceptibility to ankylosing spondylitis in twins: the role of genes, HLA, and the environment. Arthritis Rheum. 1997, 40: 1823-1828. 10.1002/art.1780401015.View ArticlePubMedGoogle Scholar
- Calin A, Marder A, Becks E, Burns T: Genetic differences between B27 positive patients with ankylosing spondylitis and B27 positive healthy controls. Arthritis Rheum. 1983, 26 (12): 1460-1464. 10.1002/art.1780261207.View ArticlePubMedGoogle Scholar
- Linden SVD, Valkenburg H, Cats A: The risk of developing ankylosing spondylitis in HLA-B27 positive individuals: a family and population study. Br. J. Rheumatol. 1983, 22: 18-19. 10.1093/rheumatology/22.1.18.View ArticlePubMedGoogle Scholar
- International Genetics of Ankylosing Spondylitis Consortium (IGAS): Identification of multiple risk variants for ankylosing spondylitis through high-density genotyping of immune-related loci. Nat Genet. 2013, 45: 730-738. 10.1038/ng.2667.View ArticleGoogle Scholar
- Reveille JD, Sims AM, Danoy P, Evans DM, Leo P, Pointon JJ, Jin R, Zhou X, Bradbury LA, Appleton LH, Davis JC, Diekman L, Doan T, Dowling A, Duan R, Duncan EL, Farrar C, Hadler J, Harvey D, Karaderi T, Mogg R, Pomeroy E, Pryce K, Taylor J, Savage L, Deloukas P, Kumanduri V, Peltonen L, Ring SM, Whittaker P, Glazov E, Thomas GP, Maksymowych WP, Inman RD, Ward MM, Stone MA, Weisman MH, Wordsworth BP, Brown MA, Australo-Anglo-American Spondyloarthritis Consortium (TASC): Genome-wide association study of ankylosing spondylitis identifies non-MHC susceptibility loci. Nat. Genet. 2010, 42: 123-127. 10.1038/ng.513.PubMed CentralView ArticlePubMedGoogle Scholar
- Conrad DF, Pinto D, Redon R, Feuk L, Gokcumen O, Zhang Y, Aerts J, Andrews TD, Barnes C, Campbell P, Fitzgerald T, Hu M, Ihm CH, Kristiansson K, Macarthur DG, Macdonald JR, Onyiah I, Pang AW, Robson S, Stirrups K, Valsesia A, Walter K, Wei J, Tyler-Smith C, Carter NP, Lee C, Scherer SW, Hurles ME, Wellcome Trust Case Control Consortium: Origins and functional impact of copy number variation in the human genome. Nature. 2010, 464 (7289): 704-712. 10.1038/nature08516.PubMed CentralView ArticlePubMedGoogle Scholar
- McKinney C, Fanciulli M, Merriman ME, Phipps-Green A, Alizadeh BZ, Koeleman BP, Dalbeth N, Gow PJ, Harrison AA, Highton J, Jones PB, Stamp LK, Steer S, Barrera P, Coenen MJ, Franke B, van Riel PL, Vyse TJ, Aitman TJ, Radstake TR, Merriman TR: Association of variation in Fcgamma receptor 3B gene copy number with rheumatoid arthritis in Caucasian samples. Ann Rheum Dis. 2010, 69 (9): 1711-1716. 10.1136/ard.2009.123588.PubMed CentralView ArticlePubMedGoogle Scholar
- Hollox EJ, Huffmeier U, Zeeuwen PL, Palla R, Lascorz J, Rodijk-Olthuis D, van de Kerkhof PC, Traupe H, de Jongh G, den Heijer M, Reis A, Armour JA, Schalkwijk J: Psoriasis is associated with increased β-defensin genomic copy number. Nat. Genet. 2008, 40 (1): 23-25. 10.1038/ng.2007.48.PubMed CentralView ArticlePubMedGoogle Scholar
- Fellermann K, Stange DE, Schaeffeler E, Schmalzl H, Wehkamp J, Bevins CL, Reinisch W, Teml A, Schwab M, Lichter P, Radlwimmer B, Stange EF: A Chromosome 8 Gene-Cluster Polymorphism with Low Human Beta-Defensin 2 Gene Copy Number Predisposes to Crohn Disease of the Colon. Am. J. Hum. Genet. 2006, 79 (3): 439-448. 10.1086/505915.PubMed CentralView ArticlePubMedGoogle Scholar
- Molokhia M, Fanciulli M, Petretto E, Patrick AL, McKeigue P, Roberts AL, Vyse TJ, Aitman TJ: FCGR3B copy number variation is associated with systemic lupus erythematosus risk in Afro-Caribbeans. Rheumatology. 2011, 50 (7): 1206-1210. 10.1093/rheumatology/keq456.PubMed CentralView ArticlePubMedGoogle Scholar
- Henrichsen CN, Chaignat E, Reymond A: Copy number variants, diseases and gene expression. Hum. Mol. Genet. 2009, 18: R1-R8. 10.1093/hmg/ddp011.View ArticlePubMedGoogle Scholar
- Xue Y, Sun D, Daly A, Yang F, Zhou X, Zhao M, Huang N, Zerjal T, Lee C, Carter NP, Hurles ME, Tyler-Smith C: Adaptive Evolution of UGT2B17 Copy-Number Variation. Am. J. Hum. Genet. 2008, 83: 337-346. 10.1016/j.ajhg.2008.08.004.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang TL, Chen XD, Guo Y, Lei SF, Wang JT, Zhou Q, Pan F, Chen Y, Zhang ZX, Dong SS, Xu XH, Yan H, Liu X, Qiu C, Zhu XZ, Chen T, Li M, Zhang H, Zhang L, Drees BM, Hamilton JJ, Papasian CJ, Recker RR, Song XP, Cheng J, Deng HW: Genome-wide Copy-Number-Variation Study Identified a Susceptibility Gene, UGT2B17, for Osteoporosis. Am. J. Hum. Genet. 2008, 83 (6): 663-674. 10.1016/j.ajhg.2008.10.006.PubMed CentralView ArticlePubMedGoogle Scholar
- Karypidis AH, Olsson M, Andersson SO, Rane A, Ekström L: Deletion polymorphism of the UGT2B17 gene is associated with increased risk for prostate cancer and correlated to gene expression in the prostate. Pharmacogenomics J. 2008, 8 (2): 147-151. 10.1038/sj.tpj.6500449.View ArticlePubMedGoogle Scholar
- Olsson M, Lindström S, Häggkvist B, Adami HO, Bälter K, Stattin P, Ask B, Rane A, Ekström L, Grönberg H: The UGT2B17 Gene Deletion is Not Associated with Prostate Cancer Risk. Prostate. 2008, 68 (5): 571-575. 10.1002/pros.20700.View ArticlePubMedGoogle Scholar
- Cooper C, Carbone L, Michet CJ, Atkinson EJ, O'Fallon WM, Melton LJ: Fracture risk in patients with ankylosing spondylitis: a population based study. J Rheumatol. 1994, 21 (10): 1877-1882.PubMedGoogle Scholar
- Vosse D, Landewé R, van der Heijde D, van der Linden S, van Staa TP, Geusens P: Ankylosing spondylitis and the risk of fracture: results from a large primary care-based nested case–control study. Ann Rheum Dis. 2009, 68 (12): 1839-1842. 10.1136/ard.2008.100503.View ArticlePubMedGoogle Scholar
- Schulze JJ, Lundmark J, Garle M, Skilving I, Ekström L, Rane A: Doping test results dependent on genotype of uridine diphospho-glucuronosyl transferase 2B17, the major enzyme for testosterone glucuronidation. J. Clin. Endocrinol. Metab. 2008, 93: 2500-2506. 10.1210/jc.2008-0218.View ArticlePubMedGoogle Scholar
- Clarke BL, Khosla S: Androgens and Bone. Steroids. 2009, 74: 296-305. 10.1016/j.steroids.2008.10.003.PubMed CentralView ArticlePubMedGoogle Scholar
- Cook EH, Scherer SW: Copy-number variations associated with neuropsychiatric conditions. Nature. 2008, 455 (7215): 919-923. 10.1038/nature07458.View ArticlePubMedGoogle Scholar
- Uddin M, Sturge M, Peddle L, O'Rielly DD, Rahman P: Genome-Wide Signatures of ‘Rearrangement Hotspots’ within Segmental Duplications in Humans. PLoS ONE. 2011, 6: e28853-10.1371.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.