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.

Effects of PTPN22 C1858T polymorphism on susceptibility and clinical characteristics of British Caucasian rheumatoid arthritis patients

Objectives. To confirm the association of a functional single-nucleotide polymorphism (SNP), C1858T (rs2476601), in the PTPN22 gene of British Caucasian rheumatoid arthritis (RA) patients and to evaluate its influence on the RA phenotype. Methods. A total of 686 RA patients and 566 healthy volunteers, all of British Caucasian origin, were genotyped for C1858T polymorphism by PCR-restriction fragment length polymorphism assay. Data were analysed using SPSS software and the χ 2 test as applicable. Results. The PTPN22 1858T risk allele was more prevalent in the RA patients (13.9%) compared with the healthy controls (10.3%) (P = 0.008, odds ratio 1.4, 95% confidence interval 1.09-1.79). The association of the T allele was restricted to those with rheumatoid factor (RF)-positive disease (n = 524, 76.4%) (P = 0.004, odds ratio 1.5, 95% confidence interval 1.1-1.9). We found no association between PTPN22 and the presence of the HLA-DRB1 shared epitope or clinical characteristics. Conclusions. We confirmed the previously reported association of PTPN22 with RF-positive RA, which was independent from the HLA-DRB1 genotype.

A review of the MHC genetics of rheumatoid arthritis

Rheumatoid arthritis is a common complex genetic disease, and, despite a significant genetic element, no gene other than HLA-DRB1 has been clearly demonstrated to be involved in the disease. However, this association accounts for less than half the overall genetic susceptibility. Investigation of other candidate genes, in particular those that reside within the major histocompatibility complex, are hampered by the presence of strong linkage disequilibrium and problems with study design. © 2004 Nature Publishing Group All rights reserved.

ERAP2 functional knockout in humans does not alter surface heavy chains or HLA-B27, inflammatory cytokines or endoplasmic reticulum stress markers

Introduction Single nucleotide polymorphisms in ERAP2 are strongly associated with ankylosing spondylitis (AS). One AS-associated single nucleotide polymorphism, rs2248374, causes a truncated ERAP2 protein that is degraded by nonsense-mediated decay. Approximately 25% of the populations of European ancestry are therefore natural ERAP2 knockouts. We investigated the effect of this associated variant on HLA class I allele presentation, surface heavy chains, endoplasmic reticulum (ER) stress markers and cytokine gene transcription in AS. Methods Patients with AS and healthy controls with either AA or GG homozygous status for rs2248374 were studied. Antibodies to CD14, CD19-ECD, HLA-A-B-C, Valpha7.2, CD161, anti-HC10 and anti-HLA-B27 were used to analyse peripheral blood mononuclear cells. Expression levels of ER stress markers (GRP78 and CHOP) and proinflammatory genes (tumour necrosis factor (TNF), IL6, IL17 and IL22) were assessed by qPCR. Results There was no significant difference in HLAclass I allele presentation or major histocompatibility class I heavy chains or ER stress markers GRP78 and CHOP or proinflammatory gene expression between genotypes for rs2248374 either between cases, between cases and controls, and between controls. Discussion Large differences were not seen in HLAB27 expression or cytokine levels between subjects with and without ERAP2 in AS cases and controls. This suggests that ERAP2 is more likely to influence AS risk through other mechanisms.

Evidence for CRHR1 in multiple sclerosis using supervised machine learning and meta-analysis in 12 566 individuals

The primary genetic risk factor in multiple sclerosis (MS) is the HLA-DRB1*1501 allele; however, much of the remaining genetic contribution to MS has yet to be elucidated. Several lines of evidence support a role for neuroendocrine system involvement in autoimmunity which may, in part, be genetically determined. Here, we comprehensively investigated variation within eight candidate hypothalamic-pituitary-adrenal (HPA) axis genes and susceptibility to MS. A total of 326 SNPs were investigated in a discovery dataset of 1343 MS cases and 1379 healthy controls of European ancestry using a multi-analytical strategy. Random Forests, a supervised machine-learning algorithm, identified eight intronic SNPs within the corticotrophin-releasing hormone receptor 1 or CRHR1 locus on 17q21.31 as important predictors of MS. On the basis of univariate analyses, six CRHR1 variants were associated with decreased risk for disease following a conservative correction for multiple tests. Independent replication was observed for CRHR1 in a large meta-analysis comprising 2624 MS cases and 7220 healthy controls of European ancestry. Results from a combined meta-analysis of all 3967 MS cases and 8599 controls provide strong evidence for the involvement of CRHR1 in MS. The strongest association was observed for rs242936 (OR = 0.82, 95% CI = 0.74-0.90, P = 9.7 × 10-5). Replicated CRHR1 variants appear to exist on a single associated haplotype. Further investigation of mechanisms involved in HPA axis regulation and response to stress in MS pathogenesis is warranted. © The Author 2010. Published by Oxford University Press. All rights reserved.

A novel serogenetic approach determines the community prevalence of celiac disease and informs improved diagnostic pathways

Background: Changing perspectives on the natural history of celiac disease (CD), new serology and genetic tests, and amended histological criteria for diagnosis cast doubt on past prevalence estimates for CD. We set out to establish a more accurate prevalence estimate for CD using a novel serogenetic approach.Methods: The human leukocyte antigen (HLA)-DQ genotype was determined in 356 patients with 'biopsy-confirmed' CD, and in two age-stratified, randomly selected community cohorts of 1,390 women and 1,158 men. Sera were screened for CD-specific serology.Results: Only five 'biopsy-confirmed' patients with CD did not possess the susceptibility alleles HLA-DQ2.5, DQ8, or DQ2.2, and four of these were misdiagnoses. HLA-DQ2.5, DQ8, or DQ2.2 was present in 56% of all women and men in the community cohorts. Transglutaminase (TG)-2 IgA and composite TG2/deamidated gliadin peptide (DGP) IgA/IgG were abnormal in 4.6% and 5.6%, respectively, of the community women and 6.9% and 6.9%, respectively, of the community men, but in the screen-positive group, only 71% and 75%, respectively, of women and 65% and 63%, respectively, of men possessed HLA-DQ2.5, DQ8, or DQ2.2. Medical review was possible for 41% of seropositive women and 50% of seropositive men, and led to biopsy-confirmed CD in 10 women (0.7%) and 6 men (0.5%), but based on relative risk for HLA-DQ2.5, DQ8, or DQ2.2 in all TG2 IgA or TG2/DGP IgA/IgG screen-positive subjects, CD affected 1.3% or 1.9%, respectively, of females and 1.3% or 1.2%, respectively, of men. Serogenetic data from these community cohorts indicated that testing screen positives for HLA-DQ, or carrying out HLA-DQ and further serology, could have reduced unnecessary gastroscopies due to false-positive serology by at least 40% and by over 70%, respectively.Conclusions: Screening with TG2 IgA serology and requiring biopsy confirmation caused the community prevalence of CD to be substantially underestimated. Testing for HLA-DQ genes and confirmatory serology could reduce the numbers of unnecessary gastroscopies. © 2013 Anderson et al.; licensee BioMed Central Ltd.

Genomewide screens in ankylosing spondylitis

Efforts to identify genes other than HLA-B27 in AS have been driven by the strength of the evidence from genetic epidemiology studies indicating that HLA-B27, although a major gene in AS, is clearly not the only significant gene operating. This is the case for both genetic determinants of disease-susceptibility and phenotypic characteristics such as disease severity and associated disease features. In this chapter the genetic epidemiology of AS and the gene-mapping studies performed to date will be reviewed and the future direction of research in this field discussed.

The interleukin 1 gene cluster contains a major susceptibility locus for ankylosing spondylitis

Ankylosing spondylitis (AS) is a common and highly heritable inflammatory arthropathy. Although the gene HLA-B27 is almost essential for the inheritance of the condition, it alone is not sufficient to explain the pattern of familial recurrence of the disease. We have previously demonstrated suggestive linkage of AS to chromosome 2q13, a region containing the interleukin 1 (IL-1) family gene cluster, which includes several strong candidates for involvement in the disease. In the current study, we describe strong association and transmission of IL-1 family gene cluster single-nucleotide polymorphisms and haplotypes with AS.

Matters Arising: Zhu et al, "A novel gene variation of TNFα associated with ankylosing spondylitis: A reconfirmed study"

In a recently published study1 involving 79 ankylosing spondylitis (AS) patients and 132 unrelated healthy blood donors, Zhu et al report association of a tumour necrosis factor α (TNFα) single nucleotide polymorphism (SNP), rs1799724, with AS.1 I have many concerns with this paper and its conclusions...

Replication of association of IL1 gene complex members with ankylosing spondylitis in Taiwanese Chinese

Objective: To test the association of interleukin 1 (IL1) gene family members with ankylosing spondylitis (AS), previously reported in Europid subjects, in an ethnically remote population. Methods: 200 Taiwanese Chinese AS patients and 200 ethnically matched healthy controls were genotyped for five single nucleotide polymorphisms (SNPs) and the IL1RN.VNTR, markers previously associated with AS. Allele, genotype, and haplotype frequencies were compared between cases and controls. Results: Association of alleles and genotypes of the markers IL1F10.3, IL1RN.4, and IL1RN.VNTR was observed with AS (p<0.05). Haplotypes of pairs of these markers and of the markers IL1RN.6/1 and IL1RN.6/2 were also significantly associated with AS. The strongest associations observed were with the marker IL1RN.4, and with the two-marker haplotype IL1RN.4-IL1RN.VNTR (both p = 0.004). Strong linkage disequilibrium was observed between all marker pairs except those involving IL1B-511 (D′ 0.4 to 0.9, p<0.01). Conclusions: The IL1 gene cluster is associated with AS in Taiwanese Chinese. This finding provides strong statistical support that the previously observed association of this gene cluster with AS is a true positive finding.

Ankylosing spondylitis is associated with the anthrax toxin receptor 2 gene (ANTXR2)

Objectives ANTXR2 variants have been associated with ankylosing spondylitis (AS) in two previous genome-wide association studies (GWAS) (p∼9×10-8). However, a genome-wide significant association (p<5×10-8) was not observed. We conducted a more comprehensive analysis of ANTXR2 in an independent UK sample to confirm and refine this association. Methods A replication study was carried out with 2978 cases and 8365 controls. Then, these were combined with non-overlapping samples from the two previous GWAS in a meta-analysis. Human leukocyte antigen (HLA)-B27 stratification was also performed to test for ANTXR2-HLA-B27 interaction. Results Out of nine single nucleotide polymorphisms (SNP) in the study, five SNPs were nominally associated (p<0.05) with AS in the replication dataset. In the meta-analysis, eight SNPs showed evidence of association, the strongest being with rs12504282 (OR=0.88, p=6.7×10-9). Seven of these SNPs showed evidence for association in the HLA-B27-positive subgroup, but none was associated with HLA-B27-negative AS. However, no statistically significant interaction was detected between HLA-B27 and ANTXR2 variants. Conclusions ANTXR2 variants are clearly associated with AS. The top SNPs from two previous GWAS (rs4333130 and rs4389526) and this study (rs12504282) are in strong linkage disequilibrium (r2≥0.76). All are located near a putative regulatory region. Further studies are required to clarify the role played by these ANTXR2 variants in AS.

Axial spondyloarthritis: A new disease entity, not necessarily early ankylosing spondylitis

New classification criteria for axial spondyloarthritis have been developed with the goal of increasing sensitivity of criteria for early inflammatory spondyloarthritis. However these criteria substantially increase heterogeneity of the resulting disease group, reducing their value in both research and clinical settings. Further research to establish criteria based on better knowledge of the natural history of non-radiographic axial spondyloarthritis, its aetiopathogenesis and response to treatment is required. In the meantime the modified New York criteria for ankylosing spondylitis remain a very useful classification criteria set, defining a relatively homogenous group of cases for clinical use and research studies.

Efficacy and safety of adalimumab in patients with non-radiographic axial spondyloarthritis: Results of a randomised placebo-controlled trial (ABILITY-1)

Purpose: To evaluate the efficacy and safety of adalimumab in patients with non-radiographic axial spondyloarthritis (nr-axSpA). Methods: Patients fulfilled Assessment of Spondyloarthritis international Society (ASAS) criteria for axial spondyloarthritis, had a Bath Ankylosing Spondylitis Disease Activity Index (BASDAI) score of ≥ 4, total back pain score of ≥ 4 (10 cm visual analogue scale) and inadequate response, intolerance or contraindication to non-steroidal anti-inflammatory drugs (NSAIDs); patients fulfilling modified New York criteria for ankylosing spondylitis were excluded. Patients were randomised to adalimumab (N=91) or placebo (N=94). The primary endpoint was the percentage of patients achieving ASAS40 at week 12. Efficacy assessments included BASDAI and Ankylosing Spondylitis Disease Activity Score (ASDAS). MRI was performed at baseline and week 12 and scored using the Spondyloarthritis Research Consortium of Canada (SPARCC) index. Results: Significantly more patients in the adalimumab group achieved ASAS40 at week 12 compared with patients in the placebo group (36% vs 15%, p<0.001). Significant clinical improvements based on other ASAS responses, ASDAS and BASDAI were also detected at week 12 with adalimumab treatment, as were improvements in quality of life measures. Inflammation in the spine and sacroiliac joints on MRI significantly decreased after 12 weeks of adalimumab treatment. Shorter disease duration, younger age, elevated baseline C-reactive protein or higher SPARCC MRI sacroiliac joint scores were associated with better week 12 responses to adalimumab. The safety profile was consistent with what is known for adalimumab in ankylosing spondylitis and other diseases. Conclusions: In patients with nr-axSpA, adalimumab treatment resulted in effective control of disease activity, decreased inflammation and improved quality of life compared with placebo. Results from ABILITY-1 suggest that adalimumab has a positive benefit-risk profile in active nr-axSpA patients with inadequate response to NSAIDs.

The X-chromosome and susceptibility to ankylosing spondylitis

Objective. Ankylosing spondylitis (AS) affects 0.25-1.0% of the population, and its etiology is incompletely understood. Susceptibility to this highly familial disease (λ(s) = 58) is primarily genetically determined. There is a significant sex bias in AS, and there are differences in recurrence risk to the offspring of affected mothers and fathers, suggesting that there may be an X-linked recessive effect. We undertook an X- chromosome linkage study to determine any contribution of the X-chromosome to AS susceptibility. Methods. A linkage study of the X-chromosome using 234 affected sibling pairs was performed to investigate this hypothesis. Results. No linkage of the X-chromosome with susceptibility to AS was found. Model- free multipoint linkage analysis strongly excluded any significant genetic contribution (λ ≥1.5) to AS susceptibility encoded on the X-chromosome (logarithm of odds [LOD] <-2.0). Smaller genetic effects (A ≥1.3) were also found to be unlikely (LOD <-1.0). Conclusion. The sex bias in AS is not explained by X-chromosome-encoded genetic effects. The disease model best explaining the sex bias in occurrence and transmission of AS is a polygenic model with a higher susceptibility threshold in females.