Identity-by-Descent Mapping to Detect Rare Variants Conferring Susceptibility to Multiple Sclerosis

Genome-wide association studies (GWAS) have identified around 60 common variants associated with multiple sclerosis (MS), but these loci only explain a fraction of the heritability of MS. Some missing heritability may be caused by rare variants that have been suggested to play an important role in the aetiology of complex diseases such as MS. However current genetic and statistical methods for detecting rare variants are expensive and time consuming. 'Population-based linkage analysis' (PBLA) or so called identity-by-descent (IBD) mapping is a novel way to detect rare variants in extant GWAS datasets. We employed BEAGLE fastIBD to search for rare MS variants utilising IBD mapping in a large GWAS dataset of 3,543 cases and 5,898 controls. We identified a genome-wide significant linkage signal on chromosome 19 (LOD = 4.65; p = 1.9×10-6). Network analysis of cases and controls sharing haplotypes on chromosome 19 further strengthened the association as there are more large networks of cases sharing haplotypes than controls. This linkage region includes a cluster of zinc finger genes of unknown function. Analysis of genome wide transcriptome data suggests that genes in this zinc finger cluster may be involved in very early developmental regulation of the CNS. Our study also indicates that BEAGLE fastIBD allowed identification of rare variants in large unrelated population with moderate computational intensity. Even with the development of whole-genome sequencing, IBD mapping still may be a promising way to narrow down the region of interest for sequencing priority. © 2013 Lin et al.

Polymorphisms in the receptor tyrosine kinase MERTK gene are associated with Multiple Sclerosis susceptibility

Multiple sclerosis (MS) is a debilitating, chronic demyelinating disease of the central nervous system affecting over 2 million people worldwide. The TAM family of receptor tyrosine kinases (TYRO3, AXL and MERTK) have been implicated as important players during demyelination in both animal models of MS and in the human disease. We therefore conducted an association study to identify single nucleotide polymorphisms (SNPs) within genes encoding the TAM receptors and their ligands associated with MS. Analysis of genotype data from a genome-wide association study which consisted of 1618 MS cases and 3413 healthy controls conducted by the Australia and New Zealand Multiple Sclerosis Genetics Consortium (ANZgene) revealed several SNPs within the MERTK gene (Chromosome 2q14.1, Accession Number NG_011607.1) that showed suggestive association with MS. We therefore interrogated 28 SNPs in MERTK in an independent replication cohort of 1140 MS cases and 1140 healthy controls. We found 12 SNPs that replicated, with 7 SNPs showing p-values of less than 10-5 when the discovery and replication cohorts were combined. All 12 replicated SNPs were in strong linkage disequilibrium with each other. In combination, these data suggest the MERTK gene is a novel risk gene for MS susceptibility. © 2011 Ma et al.

A "Candidate-Interactome" Aggregate Analysis of Genome-Wide Association Data in Multiple Sclerosis

Though difficult, the study of gene-environment interactions in multifactorial diseases is crucial for interpreting the relevance of non-heritable factors and prevents from overlooking genetic associations with small but measurable effects. We propose a "candidate interactome" (i.e. a group of genes whose products are known to physically interact with environmental factors that may be relevant for disease pathogenesis) analysis of genome-wide association data in multiple sclerosis. We looked for statistical enrichment of associations among interactomes that, at the current state of knowledge, may be representative of gene-environment interactions of potential, uncertain or unlikely relevance for multiple sclerosis pathogenesis: Epstein-Barr virus, human immunodeficiency virus, hepatitis B virus, hepatitis C virus, cytomegalovirus, HHV8-Kaposi sarcoma, H1N1-influenza, JC virus, human innate immunity interactome for type I interferon, autoimmune regulator, vitamin D receptor, aryl hydrocarbon receptor and a panel of proteins targeted by 70 innate immune-modulating viral open reading frames from 30 viral species. Interactomes were either obtained from the literature or were manually curated. The P values of all single nucleotide polymorphism mapping to a given interactome were obtained from the last genome-wide association study of the International Multiple Sclerosis Genetics Consortium & the Wellcome Trust Case Control Consortium, 2. The interaction between genotype and Epstein Barr virus emerges as relevant for multiple sclerosis etiology. However, in line with recent data on the coexistence of common and unique strategies used by viruses to perturb the human molecular system, also other viruses have a similar potential, though probably less relevant in epidemiological terms. © 2013 Mechelli et al.

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.

Five years of GWAS discovery

The past five years have seen many scientific and biological discoveries made through the experimental design of genome-wide association studies (GWASs). These studies were aimed at detecting variants at genomic loci that are associated with complex traits in the population and, in particular, at detecting associations between common single-nucleotide polymorphisms (SNPs) and common diseases such as heart disease, diabetes, auto-immune diseases, and psychiatric disorders. We start by giving a number of quotes from scientists and journalists about perceived problems with GWASs. We will then briefly give the history of GWASs and focus on the discoveries made through this experimental design, what those discoveries tell us and do not tell us about the genetics and biology of complex traits, and what immediate utility has come out of these studies. Rather than giving an exhaustive review of all reported findings for all diseases and other complex traits, we focus on the results for auto-immune diseases and metabolic diseases. We return to the perceived failure or disappointment about GWASs in the concluding section. © 2012 The American Society of Human Genetics.

Analysis of killer immunoglobulin-like receptor genes in ankylosing spondylitis

Objectives: To assess the possible association of killer immunoglobulin-like receptor (KIR) genes, specifically KIR3DL1, KIR3DS1 and KIR3DL2, with ankylosing spondylitis (AS). Methods: 14 KIR genes were genotyped in 200 UK patients with AS and 405 healthy controls using multiplex polymerase chain reaction. Sequence-specific oligonucleotide probes were used to subtype 368 cases with AS and 366 controls for 12 KIR3DL2 alleles. Differences in KIR genotypes and KIR3DL2 allele frequencies were assessed using the χp2p test. Results: KIR3DL1 and KIR3DS1 gene frequencies were very similar in cases with AS and controls (odds ratio = 1.5, 95% confidence interval 0.8 to 3.0, and odds ratio = 1.02, 95% confidence interval 0.2 to 5.3, respectively). KIR3DL2 allele frequencies were not significantly different between cases with AS and controls. Conclusions: Neither the KIR gene content of particular KIR haplotypes nor KIR3DL2 polymorphisms contribute to AS.

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.

The chromosome 16q region associated with ankylosing spondylitis includes the candidate gene tumour necrosis factor receptor type 1-associated death domain (TRADD)

Objective: To replicate and refine the reported association of ankylosing spondylitis (AS) with two nonsynonymous single nucleotide polymorphisms (nsSNPs) on chromosome 16q22.1. Methods: Firstly, 730 independent UK patients with AS were genotyped for rs9939768 and rs6979 and allele frequencies were compared with 2879 previously typed historic disease controls. Secondly, the two data sets were combined in meta-analyses. Finally, 5 tagging SNPs, located between rs9939768 and rs6979, were analysed in 1604 cases and 1020 controls. Results: The association of rs6979 with AS was replicated, p=0.03, OR=1.14 (95% CI 1.01 to 1.28), and a trend for association with rs9939768 detected, p=0.06, OR=1.25 (95% CI 0.99 to 1.57). Meta-analyses revealed association of both SNPs with AS, p=0.0008, OR=1.31 (95% CI 1.12 to 1.54) and p=0.0009, OR=1.15 (95% CI 1.06 to 1.23) for rs9939768 and rs6979, respectively. New associations with rs9033 and rs868213 (p=0.00002, OR=1.23 (95% CI 1.12 to 1.36) and p=0.00002 OR=1.45 (95% CI 1.22 to 1.72), respectively, were identified. Conclusions: The region on chromosome 16 that has been replicated in the present work is interesting as the highly plausible candidate gene, tumour necrosis factor receptor type 1 (TNFR1)-associated death domain (TRADD), is located between rs9033 and rs868213. It will require additional work to identify the primary genetic association(s) with AS.

Role of NOD2 variants in spondylarthritis

Objective. To investigate the role of the gene NOD2 in susceptibility to, and clinical manifestations of, ankylosing spondylitis (AS). Methods. A case-control study of NOD2 polymorphisms known to be associated with Crohn's disease (CD) (Pro268 Ser, Arg702 Trp, GlY908 Arg, and Len1007fsinsC) was performed in 229 cases of primary AS with no diagnosed inflammatory bowel disease (IBD), 197 cases of AS associated with IBD (referred to as colitic spondylarthritis; comprising 78 with CD and 119 with ulcerative colitis [UC]), and 229 ethnically matched, healthy controls. Associations between NOD2 polymorphisms and several clinical features of AS, including disease severity assessed by questionnaire and age at spondylarthritis onset, were also investigated. Exclusion linkage mapping of chromosome 16 was performed in a separate group of 185 multicase families with AS. Results. An association was identified between Gly908 Arg and UC spondylarthritis (P = 0.016, odds ratio [OR] 4.6, 95% confidence interval [95% CI] 1.316), and a nonsignificant trend with a similar magnitude was observed in association with CD spondylarthritis (P = 0.08, OR 3.9, 95% CI 0.8-18). The Pro268Ser variant was inversely associated with UC spondylarthritis (P = 0.003, OR 0.55, 95% CI 0.37-0.82), but not with CD spondylarthritis. No association was demonstrated between NOD2 variants and primary AS, or between other variants of NOD2 and either UC or CD spondylarthritis. Carriage of the Pro268 Ser polymorphism was associated with greater disease activity as measured by the Bath Ankylosing Spondylitis Disease Activity Index (P = 0.002). Although patients with CD had a younger age at spondylarthritis onset than did those with UC (22.4 years versus 26.4 years; P = 0.01), no association was noted between the NOD2 variants linked with CD and age at spondylarthritis onset. In primary AS, the presence of a gene with a magnitude of association >2.0 was excluded (exclusion logarithm of odds score less than -2.0), and no association was observed with the microsatellite D16S3136. Conclusion. NOD2 variants do not significantly affect the risk of developing primary AS, but may influence susceptibility to, and clinical manifestations of, colitic spondylarthritis.

High-throughput sequencing of IL23R reveals a low-frequency, nonsynonymous single-nucleotide polymorphism that is associated with ankylosing spondylitis in a Han Chinese population

Objective Ankylosing spondylitis (AS) is a highly heritable common inflammatory arthritis that targets the spine and sacroiliac joints of the pelvis, causing pain and stiffness and leading eventually to joint fusion. Although previous studies have shown a strong association of IL23R with AS in white Europeans, similar studies in East Asian populations have shown no association with common variants of IL23R, suggesting either that IL23R variants have no role or that rare genetic variants contribute. The present study was undertaken to screen IL23R to identify rare variants associated with AS in Han Chinese. Methods A 170-kb region containing IL23R and its flanking regions was sequenced in 50 patients with AS and 50 ethnically matched healthy control subjects from a Han Chinese population. In addition, the 30-kb region of peak association in white Europeans was sequenced in 650 patients with AS and 1,300 healthy controls. Validation genotyping was undertaken in 846 patients with AS and 1,308 healthy controls. Results We identified 1,047 variants, of which 729 were not found in the dbSNP genomic build 130. Several potentially functional rare variants in IL23R were identified, including one nonsynonomous single-nucleotide polymorphism (nsSNP), Gly149Arg (position 67421184 GA on chromosome 1). Validation genotyping showed that the Gly149Arg variant was associated with AS (odds ratio 0.61, P = 0.0054). Conclusion This is the first study to implicate rare IL23R variants in the pathogenesis of AS. The results identified a low-frequency nsSNP with predicted loss-of-function effects that was protectively associated with AS in Han Chinese, suggesting that decreased function of the interleukin-23 (IL-23) receptor protects against AS. These findings further support the notion that IL-23 signaling has an important role in the pathogenesis of AS.

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.

Promise and pitfalls of the Immunochip

Genomewide association studies (GWAS) have proven a powerful hypothesis-free method to identify common disease-associated variants. Even quite large GWAS, however, have only at best identified moderate proportions of the genetic variants contributing to disease heritability. To provide cost-effective genotyping of common and rare variants to map the remaining heritability and to fine-map established loci, the Immunochip Consortium has developed a 200,000 SNP chip that has been produced in very large numbers for a fraction of the cost of GWAS chips. This chip provides a powerful tool for immunogenetics gene mapping.

The window of opportunity: A relevant concept for axial spondyloarthritis

The window of opportunity is a concept critical to rheumatoid arthritis treatment. Early treatment changes the outcome of rheumatoid arthritis treatment, in that response rates are higher with earlier disease-modifying anti-rheumatic drug treatment and damage is substantially reduced. Axial spondyloarthritis is an inflammatory axial disease encompassing both nonradiographic axial spondyloarthritis and established ankylosing spondylitis. In axial spondyloarthritis, studies of magnetic resonance imaging as well as tumor necrosis factor inhibitor treatment and withdrawal studies all suggest that early effective suppression of inflammation has the potential to reduce radiographic damage. This potential would suggest that the concept of a window of opportunity is relevant not only to rheumatoid arthritis but also to axial spondyloarthritis. The challenge now remains to identify high-risk patients early and to commence treatment without delay. Developments in risk stratification include new classification criteria, identification of clinical risk factors, biomarkers, genetic associations, potential antibody associations and an ankylosing spondylitis-specific microbiome signature. Further research needs to focus on the evidence for early intervention and the early identification of high-risk individuals.

Non-major-histocompatibility-complex genetics of ankylosing spondylitis

There is strong evidence from twin and family studies indicating that a substantial proportion of the heritability of susceptibility to ankylosing spondylitis (AS) and its clinical manifestations is encoded by non-major-histocompatibility-complex genes. Efforts to identify these genes have included genomewide linkage studies and candidate gene association studies. One region, the interleukin (IL)-1 gene complex on chromosome 2, has been repeatedly associated with AS in both Caucasians and Asians. It is likely that more than one gene in this complex is involved in AS, with the strongest evidence to date implicating IL-1A. Identifying the genes underlying other linkage regions has been difficult due to the lack of obvious candidates and the low power of most studies to date to identify genes of the small to moderate magnitude that are likely to be involved. The field is moving towards genomewide association analysis, involving much larger datasets of unrelated cases and controls. Early successes using this approach in other diseases indicates that it is likely to identify genes in common diseases like AS, but there remains the risk that the common-variant, common-disease hypothesis will not hold true in AS. Nonetheless, it is appropriate for the field to be cautiously optimistic that the next few years will bring great advances in our understanding of the genetics of this condition.