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.

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.

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.

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.

Whole blood transcriptional profiling in ankylosing spondylitis identifies novel candidate genes that might contribute to the inflammatory and tissue-destructive disease aspects

Introduction: A number of genetic-association studies have identified genes contributing to ankylosing spondylitis (AS) susceptibility but such approaches provide little information as to the gene activity changes occurring during the disease process. Transcriptional profiling generates a 'snapshot' of the sampled cells' activity and thus can provide insights into the molecular processes driving the disease process. We undertook a whole-genome microarray approach to identify candidate genes associated with AS and validated these gene-expression changes in a larger sample cohort. Methods: A total of 18 active AS patients, classified according to the New York criteria, and 18 gender- and age-matched controls were profiled using Illumina HT-12 whole-genome expression BeadChips which carry cDNAs for 48,000 genes and transcripts. Class comparison analysis identified a number of differentially expressed candidate genes. These candidate genes were then validated in a larger cohort using qPCR-based TaqMan low density arrays (TLDAs). Results: A total of 239 probes corresponding to 221 genes were identified as being significantly different between patients and controls with a P-value <0.0005 (80% confidence level of false discovery rate). Forty-seven genes were then selected for validation studies, using the TLDAs. Thirteen of these genes were validated in the second patient cohort with 12 downregulated 1.3- to 2-fold and only 1 upregulated (1.6-fold). Among a number of identified genes with well-documented inflammatory roles we also validated genes that might be of great interest to the understanding of AS progression such as SPOCK2 (osteonectin) and EP300, which modulate cartilage and bone metabolism. Conclusions: We have validated a gene expression signature for AS from whole blood and identified strong candidate genes that may play roles in both the inflammatory and joint destruction aspects of the disease.

The genetics of osteoporosis: Future diagnostic possibilities

The risk of developing osteoporosis is determined by the interaction of several mostly unknown genes and environmental factors. Genetic studies in osteoporosis have largely focussed on association studies of a small number of candidate genes, with few linkage studies performed, and large areas of the genome remaining unexplored. Identifying the genes involved in osteoporosis would be a major advance in our understanding of the causation of the disease, and lead to advances in diagnosis, risk prediction, and potentially preventive and therapeutic measures.

Genetic susceptibility to ankylosing spondylitis

Ankylosing spondylitis is a highly heritable, common rheumatic condition, primarily affecting the axial skeleton. The association with HLA-B27 has been demonstrated worldwide, and evidence for a role of HLA-B27 in disease comes from linkage and association studies in humans, and transgenic animal models. However, twin studies indicate that HLA-B27 contributes only 16% of the total genetic risk for disease. Furthermore, there is compelling evidence that non-B27 genes, both within and outwith the major histocompatability complex, are involved in disease aetiology. In this post-genomic era we have the tools to help elicit the genetic basis of disease. This review describes methods for genetic investigation of ankylosing spondylitis, and summarises the status of current research in this exciting area.

Genetics of ankylosing spondylitis

Purpose of Review Over the past 3 years, several new genes and gene deserts have been identified that are associated with ankylosing spondylitis (AS). The purpose of this review is to discuss the major findings of these studies, and the answers they provide and questions they raise about the pathogenesis of this common condition. Recent Findings: Five genes/genetic regions have now definitively been associated with AS [the major histocompatibility complex (MHC), IL23R, ERAP1, 2p15 and 21q22]. Strong evidence to support association with the disease has been demonstrated for the genes IL1R2, ANTXR2, TNFSF15, TNFR1 and a region on chromosome 16q including the gene TRADD. There is an overrepresentation of genes involved in Th17 lymphocyte differentiation/activation among genes associated with AS and the related diseases inflammatory bowel disease and psoriasis, pointing strongly to this pathway as playing a major causative role in the disease. Increasing information about differential association of HLA-B27 subtypes with disease suggests a hierarchy of strength of association of those alleles with AS, providing a useful test as to the validity of different potential mechanisms of association of HLA-B27 with AS. The mechanism underlying the association of the gene deserts, 2p15 and 21q22, suggests the involvement of noncoding RNA in AS etiopathogenesis. Summary: The increasing list of genes identified as being definitely involved in AS provides a useful platform for hypothesis-driven research in the field, providing a potential alternative route to determining the underlying mechanisms involved in the disease to research focusing on HLA-B27 alone.

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.

Single nucleotide polymorphisms in clinics: Fantasy or reality for cancer?

Single nucleotide polymorphisms (SNPs) have been classically used for dissecting various human complex disorders using candidate gene studies. During the last decade, large scale SNP analysis i.e. genome-wide association studies (GWAS) have provided an agnostic approach to identify possible genetic loci associated with heterogeneous disease such as cancer susceptibility, prognosis of survival or drug response. Further, the advent of new technologies, including microarray based genotyping as well as high throughput next generation sequencing has opened new avenues for SNPs to be used in clinical practice. It is speculated that the utility of SNPs to understand the mechanisms, biology of variable drug response and ultimately treatment individualization based on the individual’s genome composition will be indispensable in the near future. In the current review, we discuss the advantages and disadvantages of the clinical utility of genetic variants in disease risk-prediction, prognosis, clinical outcome and pharmacogenomics. The lessons and challenges for the utility of SNP based biomarkers are also discussed, including the need for additional functional validation studies.

Genetic insights into common pathways and complex relationships among iImmune-mediated diseases

Shared aetiopathogenic factors among immune-mediated diseases have long been suggested by their co-familiality and co-occurrence, and molecular support has been provided by analysis of human leukocyte antigen (HLA) haplotypes and genome-wide association studies. The interrelationships can now be better appreciated following the genotyping of large immune disease sample sets on a shared SNP array: the 'Immunochip'. Here, we systematically analyse loci shared among major immune-mediated diseases. This reveals that several diseases share multiple susceptibility loci, but there are many nuances. The most associated variant at a given locus frequently differs and, even when shared, the same allele often has opposite associations. Interestingly, risk alleles conferring the largest effect sizes are usually disease-specific. These factors help to explain why early evidence of extensive 'sharing' is not always reflected in epidemiological overlap. © 2013 Macmillan Publishers Limited. All rights reserved.

Common variants in the HLA-DRB1-HLA-DQA1 HLA class II region are associated with susceptibility to visceral leishmaniasis

To identify susceptibility loci for visceral leishmaniasis, we undertook genome-wide association studies in two populations: 989 cases and 1,089 controls from India and 357 cases in 308 Brazilian families (1,970 individuals). The HLA-DRB1-HLA-DQA1 locus was the only region to show strong evidence of association in both populations. Replication at this region was undertaken in a second Indian population comprising 941 cases and 990 controls, and combined analysis across the three cohorts for rs9271858 at this locus showed P combined = 2.76 × 10 -17 and odds ratio (OR) = 1.41, 95% confidence interval (CI) = 1.30-1.52. A conditional analysis provided evidence for multiple associations within the HLA-DRB1-HLA-DQA1 region, and a model in which risk differed between three groups of haplotypes better explained the signal and was significant in the Indian discovery and replication cohorts. In conclusion, the HLA-DRB1-HLA-DQA1 HLA class II region contributes to visceral leishmaniasis susceptibility in India and Brazil, suggesting shared genetic risk factors for visceral leishmaniasis that cross the epidemiological divides of geography and parasite species. © 2013 Nature America, Inc. All rights reserved.

Identification of 15 new psoriasis susceptibility loci highlights the role of innate immunity

To gain further insight into the genetic architecture of psoriasis, we conducted a meta-analysis of 3 genome-wide association studies (GWAS) and 2 independent data sets genotyped on the Immunochip, including 10,588 cases and 22,806 controls. We identified 15 new susceptibility loci, increasing to 36 the number associated with psoriasis in European individuals. We also identified, using conditional analyses, five independent signals within previously known loci. The newly identified loci shared with other autoimmune diseases include candidate genes with roles in regulating T-cell function (such as RUNX3, TAGAP and STAT3). Notably, they included candidate genes whose products are involved in innate host defense, including interferon-mediated antiviral responses (DDX58), macrophage activation (ZC3H12C) and nuclear factor (NF)-κB signaling (CARD14 and CARM1). These results portend a better understanding of shared and distinctive genetic determinants of immune-mediated inflammatory disorders and emphasize the importance of the skin in innate and acquired host defense. © 2012 Nature America, Inc. All rights reserved.

Using genome-Wwide complex trait analysis to quantify 'missing heritability' in Parkinson's disease

Genome-wide association studies (GWASs) have been successful at identifying single-nucleotide polymorphisms (SNPs) highly associated with common traits; however, a great deal of the heritable variation associated with common traits remains unaccounted for within the genome. Genome-wide complex trait analysis (GCTA) is a statistical method that applies a linear mixed model to estimate phenotypic variance of complex traits explained by genome-wide SNPs, including those not associated with the trait in a GWAS. We applied GCTA to 8 cohorts containing 7096 case and 19 455 control individuals of European ancestry in order to examine the missing heritability present in Parkinson's disease (PD). We meta-analyzed our initial results to produce robust heritability estimates for PD types across cohorts. Our results identify 27% (95% CI 17-38, P = 8.08E - 08) phenotypic variance associated with all types of PD, 15% (95% CI -0.2 to 33, P = 0.09) phenotypic variance associated with early-onset PD and 31% (95% CI 17-44, P = 1.34E - 05) phenotypic variance associated with late-onset PD. This is a substantial increase from the genetic variance identified by top GWAS hits alone (between 3 and 5%) and indicates there are substantially more risk loci to be identified. Our results suggest that although GWASs are a useful tool in identifying the most common variants associated with complex disease, a great deal of common variants of small effect remain to be discovered. © Published by Oxford University Press 2012.

Genetic risk and a primary role for cell-mediated immune mechanisms in multiple sclerosis

Multiple sclerosis is a common disease of the central nervous system in which the interplay between inflammatory and neurodegenerative processes typically results in intermittent neurological disturbance followed by progressive accumulation of disability. Epidemiological studies have shown that genetic factors are primarily responsible for the substantially increased frequency of the disease seen in the relatives of affected individuals, and systematic attempts to identify linkage in multiplex families have confirmed that variation within the major histocompatibility complex (MHC) exerts the greatest individual effect on risk. Modestly powered genome-wide association studies (GWAS) have enabled more than 20 additional risk loci to be identified and have shown that multiple variants exerting modest individual effects have a key role in disease susceptibility. Most of the genetic architecture underlying susceptibility to the disease remains to be defined and is anticipated to require the analysis of sample sizes that are beyond the numbers currently available to individual research groups. In a collaborative GWAS involving 9,772 cases of European descent collected by 23 research groups working in 15 different countries, we have replicated almost all of the previously suggested associations and identified at least a further 29 novel susceptibility loci. Within the MHC we have refined the identity of the HLA-DRB1 risk alleles and confirmed that variation in the HLA-A gene underlies the independent protective effect attributable to the class I region. Immunologically relevant genes are significantly overrepresented among those mapping close to the identified loci and particularly implicate T-helper-cell differentiation in the pathogenesis of multiple sclerosis. © 2011 Macmillan Publishers Limited. All rights reserved.