Dissection of class III major histocompatibility complex haplotypes associated with rheumatoid arthritis

Objective. We have previously identified a single-nucleotide polymorphism (SNP) haplotype involving the lymphotoxin α (LTA) and tumor necrosis factor (TNF) loci (termed haplotype LTA-TNF2) on chromosome 6 that shows differential association with rheumatoid arthritis (RA) on HLA-DRB1*0404 and *0401 haplotypes, suggesting the presence of additional non-HLA-DRB1 RA susceptibility genes on these haplotypes. To refine this association, we performed a case-control association study using both SNPs and microsatellite markers in haplotypes matched either for HLA-DRB1*0404 or for HLA-DRB1*0401. Methods. Fourteen SNPs lying between HLA-DRB1 and LTA were genotyped in 87 DRB1*04-positive families. High-density microsatellite typing was performed using 24 markers spanning 2,500 kb centered around the TNF gene in 305 DRB1*0401 or *0404 cases and 400 DRB1*0401 or *0404 controls. Single-marker, 2-marker, and 3-marker minihaplotypes were constructed and their frequencies compared between the DRB1*0401 and DRB1*0404 matched case and control haplotypes. Results. Marked preservation of major histocompatibility complex haplotypes was seen, with chromosomes carrying LTA-TNF2 and either DRB1*0401 or DRB1*0404 both carrying an identical SNP haplotype across the 1-Mb region between TNF and HLA-DRB1. Using microsatellite markers, we observed two 3-marker minihaplotypes that were significantly overrepresented in the DRB1*0404 case haplotypes (P = 0.00024 and P = 0.00097). Conclusion. The presence of a single extended SNP haplotype between LTA-TNF2 and both DRB1*0401 and DRB1*0404 is evidence against this region harboring the genetic effects in linkage disequillbrium with LTA-TNF2. Two RA-associated haplotypes on the background of DRB1*0404 were identified in a 126-kb region surrounding and centromeric to the TNF locus.

Novel risk loci for rheumatoid arthritis in Han Chinese and congruence with risk variants in Europeans

Objective To investigate differences in genetic risk factors for rheumatoid arthritis (RA) in Han Chinese as compared with Europeans. Methods A genome-wide association study was conducted in China with 952 patients and 943 controls, and 32 variants were followed up in 2,132 patients and 2,553 controls. A transpopulation meta-analysis with results from a large European RA study was also performed to compare the genetic architecture across the 2 ethnic remote populations. Results Three non-major histocompatibility complex (non-MHC) loci were identified at the genome-wide significance level, the effect sizes of which were larger in anti-citrullinated protein antibody (ACPA)-positive patients than in ACPA-negative patients. These included 2 novel variants, rs12617656, located in an intron of DPP4 (odds ratio [OR] 1.56, P = 1.6 × 10 -21), and rs12379034, located in the coding region of CDK5RAP2 (OR 1.49, P = 1.1 × 10-16), as well as a variant at the known CCR6 locus, rs1854853 (OR 0.71, P = 6.5 × 10-15). The analysis of ACPA-positive patients versus ACPA-negative patients revealed that rs12617656 at the DPP4 locus showed a strong interaction effect with ACPAs (P = 5.3 × 10-18), and such an interaction was also observed for rs7748270 at the MHC locus (P = 5.9 × 10-8). The transpopulation meta-analysis showed genome-wide overlap and enrichment in association signals across the 2 populations, as confirmed by prediction analysis. Conclusion This study has expanded the list of alleles that confer risk of RA, provided new insight into the pathogenesis of RA, and added empirical evidence to the emerging polygenic nature of complex trait variation driven by common genetic variants. Copyright © 2014 by the American College of Rheumatology.

Is there a higher genetic load of susceptibility loci in familial ankylosing spondylitis?

Objective Several genetic risk variants for ankylosing spondylitis (AS) have been identified in genome-wide association studies. Our objective was to examine whether familial AS cases have a higher genetic load of these susceptibility variants. Methods Overall, 502 AS patients were examined, consisting of 312 patients who had first-degree relatives (FDRs) with AS (familial) and 190 patients who had no FDRs with AS or spondylarthritis (sporadic). All patients and affected FDRs fulfilled the modified New York criteria for AS. The patients were recruited from 2 US cohorts (the North American Spondylitis Consortium and the Prospective Study of Outcomes in Ankylosing Spondylitis) and from the UK-Oxford cohort. The frequencies of AS susceptibility loci in IL-23R, IL1R2, ANTXR2, ERAP-1, 2 intergenic regions on chromosomes 2p15 and 21q22, and HLA-B27 status as determined by the tag single-nucleotide polymorphism (SNP) rs4349859 were compared between familial and sporadic cases of AS. Association between SNPs and multiplex status was assessed by logistic regression controlling for sibship size. Results HLA-B27 was significantly more prevalent in familial than sporadic cases of AS (odds ratio 4.44 [95% confidence interval 2.06, 9.55], P = 0.0001). Furthermore, the AS risk allele at chromosome 21q22 intergenic region showed a trend toward higher frequency in the multiplex cases (P = 0.08). The frequency of the other AS risk variants did not differ significantly between familial and sporadic cases, either individually or combined. Conclusion HLA-B27 is more prevalent in familial than sporadic cases of AS, demonstrating higher familial aggregation of AS in patients with HLA-B27 positivity. The frequency of the recently described non-major histocompatibility complex susceptibility loci is not markedly different between the sporadic and familial cases of AS.

Progress in spondylarthritis. Progress in studies of the genetics of ankylosing spondylitis

The advent of high-throughput SNP genotyping methods has advanced research into the genetics of common complex genetic diseases such as ankylosing spondylitis (AS) rapidly in recent times. The identification of associations with the genes IL23R and ERAP1 have been robustly replicated, and advances have been made in studies of the major histocompatibility complex genetics of AS, and of KIR gene variants and the disease. The findings are already being translated into increased understanding of the immunological pathways involved in AS, and raising novel potential therapies. The current studies in AS remain underpowered, and no full genomewide association study has yet been reported in AS; such studies are likely to add to the significant advances that have already been made.

Mapping genes for osteoporosis-Old dogs and new tricks

In stark contrast to its horticultural origins, modern genetics is an extremely technology-driven field. Almost all the major advances in the field over the past 20 years have followed technological developments that have permitted change in study designs. The development of PCR in the 1980s led to RFLP mapping of monogenic diseases. The development of fluorescent-tagged genotyping methods led to linkage mapping approaches for common diseases that dominated the 1990s. The development of microarray SNP genotyping has led to the genome-wide association study era of the new millennium. And now the development of next-generation sequencing technologies is about to open up a new era of gene-mapping, enabling many potential new study designs. This review aims to present the strengths and weaknesses of the current approaches, and present some new ideas about gene-mapping approaches that are likely to advance our knowledge of the genes involved in heritable bone traits such as bone mineral density (BMD) and fracture.

Genome-wide association study for radiographic vertebral fractures: A potential role for the 16q24 BMD locus

Vertebral fracture risk is a heritable complex trait. The aim of this study was to identify genetic susceptibility factors for osteoporotic vertebral fractures applying a genome-wide association study (GWAS) approach. The GWAS discovery was based on the Rotterdam Study, a population-based study of elderly Dutch individuals aged >55years; and comprising 329 cases and 2666 controls with radiographic scoring (McCloskey-Kanis) and genetic data. Replication of one top-associated SNP was pursued by de-novo genotyping of 15 independent studies across Europe, the United States, and Australia and one Asian study. Radiographic vertebral fracture assessment was performed using McCloskey-Kanis or Genant semi-quantitative definitions. SNPs were analyzed in relation to vertebral fracture using logistic regression models corrected for age and sex. Fixed effects inverse variance and Han-Eskin alternative random effects meta-analyses were applied. Genome-wide significance was set at p<5×10-8. In the discovery, a SNP (rs11645938) on chromosome 16q24 was associated with the risk for vertebral fractures at p=4.6×10-8. However, the association was not significant across 5720 cases and 21,791 controls from 14 studies. Fixed-effects meta-analysis summary estimate was 1.06 (95% CI: 0.98-1.14; p=0.17), displaying high degree of heterogeneity (I2=57%; Qhet p=0.0006). Under Han-Eskin alternative random effects model the summary effect was significant (p=0.0005). The SNP maps to a region previously found associated with lumbar spine bone mineral density (LS-BMD) in two large meta-analyses from the GEFOS consortium. A false positive association in the GWAS discovery cannot be excluded, yet, the low-powered setting of the discovery and replication settings (appropriate to identify risk effect size >1.25) may still be consistent with an effect size <1.10, more of the type expected in complex traits. Larger effort in studies with standardized phenotype definitions is needed to confirm or reject the involvement of this locus on the risk for vertebral fractures.

PTHR1 polymorphisms influence BMD variation through effects on the growing skeleton

We investigated whether polymorphisms in PTHR1 are associated with bone mineral density (BMD), to determine whether the association of this gene with BMD was due to effects on attainment of peak bone mass or effects on subsequent bone loss. The PTHR1 gene, including its 14 exons, their exon-intron boundaries, and 1,500 bp of its promoter region, was screened for polymorphisms by denaturing high-performance liquid chromatography (dHPLC) and sequencing in 36 osteoporotic cases. Eleven single-nucleotide polymorphisms (SNPs), one tetranucleotide repeat, and one tetranucleotide deletion were identified. A cohort of 634 families, including 1,236 men (39%) and 1,926 women (61%) ascertained with probands with low BMD (Z< -2.0) and the Children in Focus subset of the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort (785 unrelated individuals, mean age 118 months), were genotyped for the five most informative SNPs (minor allele frequency >5%) and the tetranucleotide repeat. In our osteoporosis families, association was noted between lumbar spine BMD and alleles of a known functional tetranucleotide repeat (U4) in the PTHR1 promoter region (P = 0.042) and between two and three marker haplotypes of PTHR1 polymorphisms with lumbar spine, femoral neck, and total hip BMD (P = 0.021-0.047). This association was restricted to the youngest tertile of the population (age 16-39 years, P = 0.013-0.048). A similar association was found for the ALSPAC cohort: two marker haplotypes of SNPs A48609T and C52813T were associated with height (P = 0.006) and total body less head BMD (P = 0.02), corrected for age and gender, confirming the family findings. These findings suggest a role for PTHR1 variation in determining peak BMD.

Evidence of genetic association between TNFRSF1A encoding the p55 tumour necrosis factor receptor, and ankylosing spondylitis in UK Caucasians

Objectives: To replicate the possible genetic association between ankylosing spondylitis (AS) and TNFRSF1A. Methods: TNFRSF1A was re-sequenced in 48 individuals with AS to identify novel polymorphisms. Nine single nucleotide polymorphisms (SNPs) in TNFRSF1A and 5 SNPs in the neighbouring gene SCNN1A were genotyped in 1604 UK Caucasian individuals with AS and 1019 matched controls. An extended study was implemented using additional genotype data on 8 of these SNPs from 1400 historical controls from the 1958 British Birth Cohort. A meta-analysis of previously published results was also undertaken. Results: One novel variant in intron 6 was identified but no new coding variants. No definite associations were seen in the initial study but in the extended study there were weak associations with rs4149576 (p=0.04) and rs4149577 (p=0.007). In the metaanalysis consistent, somewhat stronger associations were seen with rs4149577 (p=0.002) and rs4149578 (p=0.006). Conclusions: These studies confirm the weak genetic associations between AS and TNFRSF1A. In view of the previously reported associations of TNFRSF1A with AS, in Caucasians and Chinese, and the biological plausibility of this candidate gene, replication of this finding in well powered studies is clearly indicated.

Association of IL23R and ERAP1 genes with ankylosing spondylitis in a Portuguese population

Objective: Association between ankylosing spondylitis (AS) and two genes, ERAP1 and IL23R, has recently been reported in North American and British populations. The population attributable risk fraction for ERAP1 in this study was 25%, and for IL23R, 9%. Confirmation of these findings to ERAP1 in other ethnic groups has not yet been demonstrated. We sought to test the association between single nucleotide polymorphisms (SNPs) in these genes and susceptibility to AS among a Portuguese population. We also investigated the role of these genes in clinical manifestations of AS, including age of symptom onset, the Bath Ankylosing Spondylitis Disease Activity, Metrology and Functional Indices, and the modified Stoke Ankylosing Spondylitis Spinal Score. Methods: The study was conducted on 358 AS cases and 285 ethnically matched Portuguese healthy controls. AS was defined according to the modified New York Criteria. Genotyping of IL23R and ERAP1 allelic variants was carried out with TaqMan allelic discrimination assays. Association analysis was performed using the Cochrane-Armitage and linear regression tests of genotypes as implemented in PLINK for dichotomous and quantitative variables respectively. A meta-analysis for Portuguese and previously published Spanish IL23R data was performed using the StatsDirect® Statistical tools, by fixed and random effects models. Results: A total of 14 nsSNPs markers (8 for IL23R, 5 for ERAPl, 1 for LN-PEP) were analysed. Three markers (2 for IL23R and 1 for ERAP1) showed significant single-locus disease associations, confirming that the association of these genes with AS in the Portuguese population. The strongest associated SNP in IL23R was rs1004819 (OR=1.4, p=0.0049), and in ERAP1 was rs30187 (OR=1.26, p=0.035). The population attributable risk fractions in the Portuguese population for these SNPs are 11% and 9.7% respectively. No association was seen with any SNP in LN-PEP, which flanks ERAP1 and was associated with AS in the British population. No association was seen with clinical manifestations of AS. Conclusions: These results show that IL23R and ERAP1 genes are also associated with susceptibility to AS in the Portuguese population, and that they contribute a significant proportion of the population risk for this disease.

Cuff tear arthropathy: Evidence of functional variation in pyrophosphate metabolism genes

We investigated the role of two genes, ANKH and TNAP, in patients with cuff tear arthropathy. These genes encode proteins which regulate the extracellular concentration of inorganic pyrophosphate, fluctuations of which can lead to calcium crystal formation. Variants were detected by direct sequencing of DNA and their frequencies compared with healthy controls. The effect of variants on protein function was further studied by in vitro approaches. Variant genotypes were observed more frequently in the cases when compared with controls in ANKH (45% and 20%) and TNAP (32% and 9%). Variants in ANKH altered inorganic pyrophosphate (PPi) concentrations in transfected human chondrocytes. There was a higher mean serum concentration of TNAP detected in female patients compared with normal ranges. Cuff tear arthropathy is associated with variants in ANKH and TNAP that alter extracellular inorganic pyrophosphate concentrations causing calcium crystal deposition. This supports a theory that genetic variants predispose patients to primary crystal deposition which when combined with a massive rotator cuff tear leads to the development of arthritis.

Genetics and the pathogenesis of ankylosing spondylitis

Purpose of review The field of genetic research in ankylosing spondylitis (AS) is advancing rapidly. The purpose of this review is to outline recent findings, particularly, in regard to genetic studies of the major histocompatibility complex (MHC) and the non-MHC genes IL23R, ERAP1, and killer cell immunologlobulin-like receptor (KIR) complex, in AS. Recent findings: Convincing evidence has been reported for the existence of further non-B27 MHC genes involved in AS. Strong, replicated association has been reported with IL23R and ERAP1 and AS. The IL23R finding strongly implicates the TH17 lymphocyte system in AS aetiopathogenesis. Suggestive evidence of a role for KIR gene polymorphism in AS exists, but definitive findings are awaited. Summary: The findings suggest that further genome-wide studies in large case-control cohorts are likely to be very productive in this disease. The IL23R findings and subsequent immunological investigations suggest that targeted intervention in the TH17 system is likely to have major therapeutic benefit, as it does in the genetically related diseases, inflammatory bowel disease and psoriasis.

Genetic determinants of heel bone properties: Genome-wide association meta-analysis and replication in the GEFOS/GENOMOS consortium

Quantitative ultrasound of the heel captures heel bone properties that independently predict fracture risk and, with bone mineral density (BMD) assessed by X-ray (DXA), may be convenient alternatives for evaluating osteoporosis and fracture risk. We performed a meta-analysis of genome-wide association (GWA) studies to assess the genetic determinants of heel broadband ultrasound attenuation (BUA; n 5 14 260), velocity of sound (VOS; n 5 15 514) and BMD (n 5 4566) in 13 discovery cohorts. Independent replication involved seven cohorts with GWA data (in silico n 5 11 452) and new genotyping in 15 cohorts (de novo n 5 24 902). In combined random effects, meta-analysis of the discovery and replication cohorts, nine single nucleotide polymorphisms (SNPs) had genome-wide significant (P < 5 3 108) associations with heel bone properties. Alongside SNPs within or near previously identified osteoporosis susceptibility genes including ESR1 (6q25.1: rs4869739, rs3020331, rs2982552), SPTBN1 (2p16.2: rs11898505), RSPO3 (6q22.33: rs7741021), WNT16 (7q31.31: rs2908007), DKK1 (10q21.1: rs7902708) and GPATCH1 (19q13.11: rs10416265), we identified a new locus on chromosome 11q14.2 (rs597319 close to TMEM135, a gene recently linked to osteoblastogenesis and longevity) significantly associated with both BUA and VOS (P < 8.23 3 1014). In meta-analyses involving 25 cohorts with up to 14 985 fracture cases, six of 10 SNPs associated with heel bone properties at P < 5 3 106 also had the expected direction of association with any fracture (P < 0.05), including threeSNPswithP < 0.005: 6q22.33 (rs7741021), 7q31.31 (rs2908007) and 10q21.1 (rs7902708). In conclusion, thisGWAstudy reveals the effect of several genescommon to central DXA-derivedBMDand heel ultrasound/DXAmeasures and points to anewgenetic locus with potential implications for better understanding of osteoporosis pathophysiology.

Utility of temporal artery biopsy samples for genome-wide analysis of giant cell arteritis

Giant Cell Arteritis (GCA) is the most common vasculitis affecting the elderly. Archived formalin-fixed paraffin-embedded (FFPE) temporal artery biopsy (TAB) specimens potentially represent a valuable resource for large-scale genetic analysis of this disease. FFPE TAB samples were obtained from 12 patients with GCA. Extracted TAB DNA was assessed by real time PCR before restoration using the Illumina HD FFPE Restore Kit. Paired FFPE-blood samples were genotyped on the Illumina OmniExpress FFPE microarray. The FFPE samples that passed stringent quality control measures had a mean genotyping success of >97%. When compared with their matching peripheral blood DNA, the mean discordant heterozygote and homozygote single nucleotide polymorphisms calls were 0.0028 and 0.0003, respectively, which is within the accepted tolerance of reproducibility. This work demonstrates that it is possible to successfully obtain high-quality microarray-based genotypes FFPE TAB samples and that this data is similar to that obtained from peripheral blood.

Genetic associations and functional characterization of M1 aminopeptidases and immune-mediated diseases

Endosplasmic reticulum aminopeptidase 1 (ERAP1), endoplasmic reticulum aminopeptidase 2 (ERAP2) and puromycin-sensitive aminopeptidase (NPEPPS) are key zinc metallopeptidases that belong to the oxytocinase subfamily of M1 aminopeptidase family. NPEPPS catalyzes the processing of proteosome-derived peptide repertoire followed by trimming of antigenic peptides by ERAP1 and ERAP2 for presentation on major histocompatibility complex (MHC) Class I molecules. A series of genome-wide association studies have demonstrated associations of these aminopeptidases with a range of immune-mediated diseases such as ankylosing spondylitis, psoriasis, Behçet's disease, inflammatory bowel disease and type I diabetes, and significantly, genetic interaction between some aminopeptidases and HLA Class I loci with which these diseases are strongly associated. In this review, we highlight the current state of understanding of the genetic associations of this class of genes, their functional role in disease, and potential as therapeutic targets.

Schizophrenia Genetic Variants are not Associated with Intelligence

Background Schizophrenia is associated with lower pre-morbid intelligence (IQ) in addition to (pre-morbid) cognitive decline. Both schizophrenia and IQ are highly heritable traits. Therefore, we hypothesized that genetic variants associated with schizophrenia, including copy number variants (CNVs) and a polygenic schizophrenia (risk) score (PSS), may influence intelligence. Method IQ was estimated with the Wechsler Adult Intelligence Scale (WAIS). CNVs were determined from single nucleotide polymorphism (SNP) data using the QuantiSNP and PennCNV algorithms. For the PSS, odds ratios for genome-wide SNP data were calculated in a sample collected by the Psychiatric Genome-Wide Association Study (GWAS) Consortium (8690 schizophrenia patients and 11 831 controls). These were used to calculate individual PSSs in our independent sample of 350 schizophrenia patients and 322 healthy controls. Results Although significantly more genes were disrupted by deletions in schizophrenia patients compared to controls (p = 0.009), there was no effect of CNV measures on IQ. The PSS was associated with disease status (R 2 = 0.055, p = 2.1 × 10 -7) and with IQ in the entire sample (R 2 = 0.018, p = 0.0008) but the effect on IQ disappeared after correction for disease status. Conclusions Our data suggest that rare and common schizophrenia-associated variants do not explain the variation in IQ in healthy subjects or in schizophrenia patients. Thus, reductions in IQ in schizophrenia patients may be secondary to other processes related to schizophrenia risk. © Cambridge University Press 2013.