The genetic basis of spondyloarthritis: SPARTAN/IGAS 2009

A joint meeting was held in July 2009 in Houston, Texas, of members of the Spondyloarthritis Research and Therapy Network (SPARTAN), founded in 2003 to promote research, education, and treatment of ankylosing spondylitis (AS) and related forms of spondyloarthritis (SpA), and members of International Genetics of AS (IGAS), founded in 2003 to encourage and coordinate studies internationally in the genetics of AS. The general topic was the genetic basis of SpA, with presentations on the future of human genetic studies; microbes, SpA, and innate immunity; susceptibility of AS to the major histocompatibility complex (MHC) and non-MHC; and individual discussions of the genetics of psoriasis and psoriatic arthritis, uveitis, inflammatory bowel disease, and enteropathic arthritis. Summaries of those discussions are presented.

Genetic determinants of bone density and fracture risk: State of the art and future directions

Context: Osteoporosis is a common, highly heritable condition that causes substantial morbidity and mortality, the etiopathogenesis of which is poorly understood. Genetic studies are making increasingly rapid progress in identifying the genes involved. Evidence Acquisition and Synthesis: In this review, we will summarize the current understanding of the genetics of osteoporosis based on publications from PubMed from the year 1987 onward. Conclusions: Most genes involved in osteoporosis identified to date encode components of known pathways involved in bone synthesis or resorption, but as the field progresses, new pathways are being identified. Only a small proportion of the total genetic variation involved in osteoporosis has been identified, and new approaches will be required to identify most of the remaining genes.

Investigating the genetic association between ERAP1 and ankylosing spondylitis

A strong association between ERAP1 and ankylosing spondylitis (AS) was recently identified by the Wellcome Trust Case Control Consortium and the Australo-Anglo-American Spondylitis Consortium (WTCCC-TASC) study. ERAP1 is highly polymorphic with strong linkage disequilibrium evident across the gene. We therefore conducted a series of experiments to try to identify the primary genetic association(s) with ERAP1. We replicated the original associations in an independent set of 730 patients and 1021 controls, resequenced ERAP1 to define the full extent of coding polymorphisms and tested all variants in additional association studies. The genetic association with ERAP1 was independently confirmed; the strongest association was with rs30187 in the replication set (P = 3.4 × 103). When the data were combined with the original WTCCC-TASC study the strongest association was with rs27044 (P = 1.1 × 10-9). We identified 33 sequence polymorphisms in ERAP1, including three novel and eight known non-synonymous polymorphisms. We report several new associations between AS and polymorphisms distributed across ERAP1 from the extended case-control study, the most significant of which was with rs27434 (P = 4.7 × 10-7). Regression analysis failed to identify a primary association clearly; we therefore used data from HapMap to impute genotypes for an additional 205 non-coding SNPs located within and adjacent to ERAP1. A number of highly significant associations (P < 5 × 10-9) were identified in regulatory sequences which are good candidates for causing susceptibility to AS, possibly by regulating ERAP1 expression. © 2009 The Author(s).

Genetic analyses in a sample of individuals with high or low BMD shows association with multiple Wnt pathway genes

INTRODUCTION Although the high heritability of BMD variation has long been established, few genes have been conclusively shown to affect the variation of BMD in the general population. Extreme truncate selection has been proposed as a more powerful alternative to unselected cohort designs in quantitative trait association studies. We sought to test these theoretical predictions in studies of the bone densitometry measures BMD, BMC, and femoral neck area, by investigating their association with members of the Wnt pathway, some of which have previously been shown to be associated with BMD in much larger cohorts, in a moderate-sized extreme truncate selected cohort (absolute value BMD Z-scores = 1.5-4.0; n = 344). MATERIALS AND METHODS Ninety-six tag-single nucleotide polymorphism (SNPs) lying in 13 Wnt signaling pathway genes were selected to tag common genetic variation (minor allele frequency [MAF] > 5% with an r(2) > 0.8) within 5 kb of all exons of 13 Wnt signaling pathway genes. The genes studied included LRP1, LRP5, LRP6, Wnt3a, Wnt7b, Wnt10b, SFRP1, SFRP2, DKK1, DKK2, FZD7, WISP3, and SOST. Three hundred forty-four cases with either high or low BMD were genotyped by Illumina Goldengate microarray SNP genotyping methods. Association was tested either by Cochrane-Armitage test for dichotomous variables or by linear regression for quantitative traits. RESULTS Strong association was shown with LRP5, polymorphisms of which have previously been shown to influence total hip BMD (minimum p = 0.0006). In addition, polymorphisms of the Wnt antagonist, SFRP1, were significantly associated with BMD and BMC (minimum p = 0.00042). Previously reported associations of LRP1, LRP6, and SOST with BMD were confirmed. Two other Wnt pathway genes, Wnt3a and DKK2, also showed nominal association with BMD. CONCLUSIONS This study shows that polymorphisms of multiple members of the Wnt pathway are associated with BMD variation. Furthermore, this study shows in a practical trial that study designs involving extreme truncate selection and moderate sample sizes can robustly identify genes of relevant effect sizes involved in BMD variation in the general population. This has implications for the design of future genome-wide studies of quantitative bone phenotypes relevant to osteoporosis.

Fine mapping of the MHC Class III region demonstrates association of AIF1 and rheumatoid arthritis

Objective. The heritability of RA has been estimated to be ∼55%, of which the MHC contributes about one-third. HLA-DRB1 alleles are strongly associated with RA, but it is likely that significant non-DRB1 MHC genetic susceptibility factors are involved. Previously, we identified two three-marker haplotypes in a 106-kb region in the MHC class III region immediately centromeric to TNF, which are strongly associated with RA on HLA-DRB1*0404 haplotypes. In the present study, we aimed to refine these associations further using a combination of genotyping and gene expression studies. Methods. Thirty-nine nucleotide polymorphisms (SNPs) were genotyped in 95 DRB1*0404 carrying unrelated RA cases, 125 DRB1*0404 - carrying healthy controls and 87 parent-case trio RA families in which the affected child carried HLA-DRB1*04. Quantitative RT-PCR was used to assess the expression of the positional candidate MHC class III genes APOM, BAT2, BAT3, BAT4, BAT5, AIF1, C6orf47, CSNK2β and LY6G5C, and the housekeeper genes, hypoxanthine-guanine phosphoribosyltransferase (HPRT) and β2-microglobulin (B2M) in 31 RA cases and 21 ethnically, age- and sex-matched healthy controls. Synovial membrane specimens from RA, PsA and OA cases were stained by an indirect immunoperoxidase technique using a mouse-anti-human AIF1 monoclonal antibody. Results. Association was observed between RA and single markers or two marker haplotypes involving AIF1, BAT3 and CSNK. AIF1 was also significantly overexpressed in RA mononuclear cells (1.5- to 1.9-fold difference, P = 0.02 vs HPRT, P = 0.002 vs B2M). AIF1 protein was clearly expressed by synovial macrophages in all the inflammatory synovial samples in contrast to the non-inflammatory OA samples. Conclusions. The results of the genotyping and expression studies presented here suggest a role for AIF1 in both the aetiology and pathogenesis of RA.

Breakthroughs in genetic studies of ankylosing spondylitis

Ankylosing spondylitis (AS), the prototypic seronegative arthropathy, is known to be highly heritable, with >90% of the risk of developing the disease determined genetically. As with most common heritable diseases, progress in identifying the genes involved using family-based or candidate gene approaches has been slow. The recent development of the genome-wide association study approach has revolutionized genetic studies of such diseases. Early studies in ankylosing spondylitis have produced two major breakthroughs in the identification of genes contributing roughly one third of the population attributable risk of the disease, and pointing directly to a potential therapy. These exciting findings highlight the potential of future more comprehensive genetic studies of determinants of disease risk and clinical manifestations, and are the biggest advance in our understanding of the causation of the disease since the discovery of the association with HLA-B27.

Investigation of the role of ENPP1 and TNAP genes in chondrocalcinosis

Objectives. Extracellular inorganic pyrophosphate (ePPi) inhibits certain forms of pathological mineralization while promoting others. Three molecules involved in ePPi regulation are important candidates for the development of calcium pyrophosphate dihydrate chondrocalcinosis (CPPD CC). These include ANKH, ectonucleotide pyrophosphatase (ENPP1) and TNAP. We have previously showed that genetic variation in ANKH is a cause of autosomal dominant familial CC and also some sporadic cases of CPPD CC. We now investigate the possible role of ENPP1 and TNAP in CPPD CC. Methods. Exons, untranslated regions (UTR) and exon-intron boundaries of ENPP1 and TNAP were sequenced using ABI Big Dye chemistry on automated sequencers. Sixteen variants were identified (3 in ENPP1 and 13 in TNAP) and were subsequently genotyped in 128 sporadic Caucasian CPPD CC patients and 600 healthy controls using a combination of polymerase chain reaction/restriction fragment-length polymorphism analysis or using Taqman. Allele and genotype frequencies were compared between cases and controls using the χ 2 test. Linkage disequilibrium, haplotype and the single nucleotide polymorphism-specific analyses were also performed. This study had 80% power to detect an odds ratio of 2.2 or more at these loci. Results. No difference was observed in the allele or genotype frequencies between patients and controls at either ENPP1 or TNAP. Conclusions. Polymorphisms of ENPP1 and TNAP are not major determinants of susceptibility to CC in the population studied. Further studies of the aetiology of sporadic CPPD CC are required to determine its causes.

Estrogen receptor α regulates area-adjusted bone mineral content in late pubertal girls

Context: Whether the action of estrogen in skeletal development depends on estrogen receptor α as encoded by the ESR1 gene is unknown. Objectives: The aim of this study was to establish whether the gain in area-adjusted bone mineral content (ABMC) in girls occurs in late puberty and to examine whether the magnitude of this gain is related to ESR1 polymorphisms. Design: We conducted a cross-sectional analysis. Setting: The study involved the Avon Longitudinal Study of Parents and Children (ALSPAC), a population-based prospective study. Participants: Participants included 3097 11-yr-olds with DNA samples, dual x-ray absorptiometry measurements, and pubertal stage information. Outcomes: Outcome measures included separate prespecified analyses in boys and girls of the relationship between ABMC derived from total body dual x-ray absorptiometry scans and Tanner stage and of the interaction between ABMC, Tanner stage, and ESR1 polymorphisms. Results: Total body less head and spinal ABMC were higher in girls in Tanner stages 4 and 5, compared with those in Tanner stages 1, 2, and 3. In contrast, height increased throughout puberty. No differences were observed in ABMC according to Tanner stage in boys. For rs2234693 (PvuII) and rs9340799 (XbaI) polymorphisms, differences in spinal ABMC in late puberty were 2-fold greater in girls who were homozygous for the C and G alleles, respectively (P = 0.001). For rs7757956, the difference in total body less head ABMC in late puberty was 50% less in individuals homozygous or heterozygous for the A allele (P = 0.006). Conclusions: Gains in ABMC in late pubertal girls are strongly associated with ESR1 polymorphisms, suggesting that estrogen contributes to this process via an estrogen receptor α-dependent pathway.

Effect of an estrogen receptor-α intron 4 polymorphism on fat mass in 11-year-old children

Context: in the ESR1 gene encoding estrogen receptor (ER)-α may be associated with fat mass in adults. Objectives: The objective of the study was to establish whether ESR1 polymorphisms influence fat mass in childhood. Design: This was a cross-sectional analysis after genotyping of rs9340799, rs2234693, and rs7757956 ESR1 polymorphisms. Setting: The Avon Longitudinal Study of Parents and Children (ALSPAC) was a population-based prospective study. Participants: Participants included 3097 11-yr-old children with results for ESR1 genotyping, puberty measures, and dual-energy x-ray absorptiometry results. Outcomes: Relationships between ESR1 polymorphisms and indices of body composition were measured. Results: The rs7757956 polymorphism was associated with fat mass (P = 0.002). Total body fat mass (adjusted for height) was reduced by 6% in children with TA/AA genotypes, and risk of being overweight (≥85th centile of fat mass) was decreased by 20%. This genetic effect appeared to interact with puberty in girls (P = 0.05 for interaction): in those with the TT genotype, total body fat mass (adjusted for height) was 18% higher in Tanner stages 3-5 vs. stages 1-2; the equivalent difference was 7% in those with TA/AA genotypes. Furthermore, the risk of being overweight was 36% lower in girls with TA/AA genotypes in Tanner stages 3-5, but no reduction was seen in those in stages 1-2. Neither rs9340799 nor rs2234693 polymorphisms were associated with body composition measures. Conclusions: Fat mass in 11-yr-old children was related to the rs7757956 ESR1 polymorphism. This association was strongest in girls in more advanced puberty, in whom the risk of being overweight was reduced by 36% in those with the TA/AA genotype.

Genetic, functional, and histopathological evaluation of two C-terminal BRCA1 missense variants

Background: The vast majority of BRCA1 missense sequence variants remain uncharacterised for their possible effect on protein expression and function, and therefore are unclassified in terms of their pathogenicity. BRCA1 plays diverse cellular roles and it is unlikely that any single functional assay will accurately reflect the total cellular implications of missense mutations in this gene. Objective: To elucidate the effect of two BRCA1 variants, 5236G>C (G1706A) and 5242C>A (A1708E) on BRCA1 function, and to survey the relative usefulness of several assays to direct the characterisation of other unclassified variants in BRCA genes. Methods and Results: Data from a range of bioinformatic, genetic, and histopathological analyses, and in vitro functional assays indicated that the 1708E variant was associated with the disruption of different cellular functions of BRCA1. In transient transfection experiments in T47D and 293T cells, the 1708E product was mislocalised to the cytoplasm and induced centrosome amplification in 293T cells. The 1708E variant also failed to transactivate transcription of reporter constructs in mammalian transcriptional transactivation assays. In contrast, the 1706A variant displayed a phenotype comparable to wildtype BRCA1 in these assays. Consistent with functional data, tumours from 1708E carriers showed typical BRCA1 pathology, while tumour material from 1706A carriers displayed few histopathological features associated with BRCA1 related tumours. Conclusions: A comprehensive range of genetic, bioinformatic, and functional analyses have been combined for the characterisation of BRCA1 unclassified sequence variants. Consistent with the functional analyses, the combined odds of causality calculated for the 1706A variant after multifactorial likelihood analysis (1:142) indicates a definitive classification of this variant as "benign". In contrast, functional assays of the 1708E variant indicate that it is pathogenic, possibly through subcellular mislocalisation. However, the combined odds of 262:1 in favour of causality of this variant does not meet the minimal ratio of 1000:1 for classification as pathogenic, and A1708E remains formally designated as unclassified. Our findings highlight the importance of comprehensive genetic information, together with detailed functional analysis for the definitive categorisation of unclassified sequence variants. This combination of analyses may have direct application to the characterisation of other unclassified variants in BRCA1 and BRCA2.

The v-MFG test: Investigating maternal, offspring and maternal-fetal genetic incompatibility effects on disease and viability

The MFG test is a family-based association test that detects genetic effects contributing to disease in offspring, including offspring allelic effects, maternal allelic effects and MFG incompatibility effects. Like many other family-based association tests, it assumes that the offspring survival and the offspring-parent genotypes are conditionally independent provided the offspring is affected. However, when the putative disease-increasing locus can affect another competing phenotype, for example, offspring viability, the conditional independence assumption fails and these tests could lead to incorrect conclusions regarding the role of the gene in disease. We propose the v-MFG test to adjust for the genetic effects on one phenotype, e.g., viability, when testing the effects of that locus on another phenotype, e.g., disease. Using genotype data from nuclear families containing parents and at least one affected offspring, the v-MFG test models the distribution of family genotypes conditional on offspring phenotypes. It simultaneously estimates genetic effects on two phenotypes, viability and disease. Simulations show that the v-MFG test produces accurate genetic effect estimates on disease as well as on viability under several different scenarios. It generates accurate type-I error rates and provides adequate power with moderate sample sizes to detect genetic effects on disease risk when viability is reduced. We demonstrate the v-MFG test with HLA-DRB1 data from study participants with rheumatoid arthritis (RA) and their parents, we show that the v-MFG test successfully detects an MFG incompatibility effect on RA while simultaneously adjusting for a possible viability loss.

Genetic disorders of the LRP5-Wnt signalling pathway affecting the skeleton

Osteoporosis is a common, increasingly prevalent and potentially debilitating condition of men and women. Genetic factors are major determinants of bone mass and the risk of fracture, but few genes have been definitively demonstrated to be involved. The identification of these factors will provide novel insights into the processes of bone formation and loss and thus the pathogenesis of osteoporosis, enabling the rational development of novel therapies. In this article, we present the extensive genetic and functional data indicating that the LRP5 gene and the Wnt signalling pathway are key players in bone formation and the risk of osteoporosis, and that LRP5 signalling is essential for normal morphology, developmental processes and bone health.

Genetic and genomic studies of PADI4 in rheumatoid arthritis

Objectives. Strong genetic association of rheumatoid arthritis (RA) with PADI4 (peptidyl arginine deiminase) has previously been described in Japanese, although this was not confirmed in a subsequent study in the UK. We therefore undertook a further study of genetic association between PADI4 and RA in UK Caucasians and also studied expression of PADI4 in the peripheral blood of patients with RA. Methods. Seven single-nucleotide polymorphisms (SNP) were genotyped using polymerase chain reaction (PCR)-restriction fragment length polymorphism in 111 RA cases and controls. A marker significantly associated with RA (PADI4_100, rs#2240339) in this first data set (P = 0.03) was then tested for association in a larger group of 439 RA patients and 428 controls. PADI4 transcription was also assessed by real-time quantitative PCR using RNA extracted from peripheral blood mononuclear cells from 13 RA patients and 11 healthy controls. Results. A single SNP was weakly associated with RA (P = 0.03) in the initial case-control study, a single SNP (PADI4_100) and a two marker haplotype of that SNP and the neighbouring SNP (PADI4_04) were significantly associated with RA (P = 0.02 and P = 0.03 respectively). PADI4_100 was not associated with RA in a second sample set. PADI4 expression was four times greater in cases than controls (P = 0.004), but expression levels did not correlate with the levels of markers of inflammation. Conclusion. PADI4 is significantly overexpressed in the blood of RA patients but genetic variation within PADI4 is not a major risk factor for RA in Caucasians.

Site and gender specificity of inheritance of bone mineral density

Differences in genetic control of BMD by skeletal sites and genders were examined by complex segregation analysis in 816 members of 147 families with probands with extreme low BMD. Spine BMD correlated more strongly in male-male comparisons and hip BMD in female-female comparisons, consistent with gender- and site-specificity of BMD heritability. Introduction: Evidence from studies in animals and humans suggests that the genetic control of bone mineral density (BMD) may differ at different skeletal sites and between genders. This question has important implications for the design and interpretation of genetic studies of osteoporosis. Methods: We examined the genetic profile of 147 families with 816 individuals recruited through probands with extreme low BMD (T-score < −2.5, Z-score < −2.0). Complex segregation analysis was performed using the Pedigree Analysis Package. BMD was measured by DXA at both lumbar spine (L1-L4) and femoral neck. Results: Complex segregation analysis excluded purely monogenic and environmental models of segregation of lumbar spine and femoral neck BMD in these families. Pure polygenic models were excluded at the lumbar spine when menopausal status was considered as a covariate, but not at the femoral neck. Mendelian models with a residual polygenic component were not excluded. These models were consistent with the presence of a rare Mendelian genotype of prevalence 3–19 %, causing high BMD at the hip and spine in these families, with additional polygenic effects. Total heritability range at the lumbar spine was 61–67 % and at the femoral neck was 44–67 %. Significant differences in correlation of femoral neck and lumbar spine BMD were observed between male and female relative pairs, with male-male comparisons exhibiting stronger lumbar spine BMD correlation than femoral neck, and female-female comparisons having greater femoral neck BMD correlation than lumbar spine. These findings remained true for parent-offspring correlations when menopausal status was taken into account. The recurrence risk ratio for siblings of probands of a Z-score < −2.0 was 5.4 at the lumbar spine and 5.9 at the femoral neck. Conclusions: These findings support gender- and site-specificity of the inheritance of BMD. These results should be considered in the design and interpretation of genetic studies of osteoporosis.

Genetic studies of disorders of calcium crystal deposition

In summary, although many factors are likely to be involved in regulating calcification and ossification processes, studies of the causation of articular chondrocalcinosis and disorders of spinal ossification, such as DISH and OPLL, implicate control over inorganic pyrophosphate levels as being one of the most important factors in their aetiopathogenesis. The findings of these studies may prove relevant to other rheumatic diseases in which ectopic ossification occurs, such as AS.