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.

Multiple alignment and sorting of peptides derived from phage-displayed random peptide libraries with polyclonal sera allows discrimination of relevant phagotopes

Biopanning of phage-displayed random peptide libraries is a powerful technique for identifying peptides that mimic epitopes (mimotopes) for monoclonal antibodies (mAbs). However, peptides derived using polyclonal antisera may represent epitopes for a diverse range of antibodies. Hence following screening of phage libraries with polyclonal antisera, including autoimmune disease sera, a procedure is required to distinguish relevant from irrelevant phagotopes. We therefore applied the multiple sequence alignment algorithm PILEUP together with a matrix for scoring amino acid substitutions based on physicochemical properties to generate guide trees depicting relatedness of selected peptides. A random heptapeptide library was biopanned nine times using no selecting antibodies, immunoglobulin G (IgG) from sera of subjects with autoimmune diseases (primary biliary cirrhosis (PBC) and type 1 diabetes) and three murine ascites fluids that contained mAbs to overlapping epitope(s) on the Ross River Virus envelope protein 2. Peptides randomly sampled from the library were distributed throughout the guide tree of the total set of peptides whilst many of the peptides derived in the absence of selecting antibody aligned to a single cluster. Moreover peptides selected by different sources of IgG aligned to separate clusters, each with a different amino acid motif. These alignments were validated by testing all of the 53 phagotopes derived using IgG from PBC sera for reactivity by capture ELISA with antibodies affinity purified on the E2 subunit of the pyruvate dehydrogenase complex (PDC-E2), the major autoantigen in PBC: only those phagotopes that aligned to PBC-associated clusters were reactive. Hence the multiple sequence alignment procedure discriminates relevant from irrelevant phagotopes and thus a major difficulty with biopanning phage-displayed random peptide libraries with polyclonal antibodies is surmounted.

Genetics of ankylosing spondylitis and rheumatoid arthritis: Where are we at currently, and how do they compare?

Both ankylosing spondylitis (AS) and rheumatoid arthritis (RA) are common, highly heritable conditions, the pathogenesis of which are incompletely understood. Gene-mapping studies in both conditions have over the last couple of years made major breakthroughs in identifying the mechanisms by which these diseases occur. Considering RA, there is an over-representation of genes involved in TNF signalling and the NFκB pathway that have been shown to influence the disease risk. There is also considerable sharing of susceptibility genes between RA and other autoimmune diseases such as systemic lupus erythematosus, type 1 diabetes, autoimmune thyroid disease and celiac disease, with thus far little overlap with AS. In AS, genes involved in response to IL12/IL23, and in endoplasmic reticulum peptide presentation, have been identified, but a full genomewide association study has not yet been reported.

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.

Fine mapping and replication of genetic risk Loci in primary sclerosing cholangitis

Background and aims. Primary sclerosing cholangitis (PSC) is a chronic cholestatic liver disease characterized by progressive inflammation and fibrosis of the bile ducts eventually leading to biliary cirrhosis. Recent genetic studies in PSC have identified associations at 2q13, 2q35, 3p21, 4q27, 13q31 and suggestive association at 10p15. The aim of this study was to further characterize and refine the genetic architecture of PSC. Methods. We analyzed previously reported associated SNPs at four of these non-HLA loci and 59 SNPs tagging the IL-2/IL-21 (4q27) and IL2RA (10p15) loci in 992 UK PSC cases and 5162 healthy UK controls. Results. The most associated SNPs identified were rs3197999 (3p21 (MST1), p = 1.9 × 10 -6, OR A vs G = 1.28, 95% CI (1.16-1.42)); rs4147359 (10p15 (IL2RA), p = 2.6 × 10 -4, OR A vs G = 1.20, 95% CI (1.09-1.33)) and rs12511287 (4q27 (IL-2/IL-21), p = 3.0 × 10 -4, OR A vs T = 1.21, 95% CI (1.09-1.35)). In addition, we performed a meta-analysis for selected SNPs using published summary statistics from recent studies. We observed genome-wide significance for rs3197999 (3p21 (MST1), P combined = 3.8 × 10 -12) and rs4147359 (10p15 (IL2RA), P combined = 1.5 × 10 -8). Conclusion. We have for the first time confirmed the association of PSC with genetic variants at 10p15 (IL2RA) locus at genome-wide significance and replicated the associations at MST1 and IL-2/IL-21 loci in a large homogeneous UK population. These results strongly implicate the role of IL-2/IL2RA pathway in PSC and provide further confirmation of MST1 association. © Informa Healthcare.

Interleukin 10 polymorphisms in ankylosing spondylitis

Genetic polymorphisms of the IL10 promoter region have been implicated in many autoimmune diseases, including seronegative spondyloarthropathies. We studied three SNPs (IL10-1087,-824, and -597) and two microsafellites(IL10R and IL10G) lying within the promoter region of IL10 for association with susceptibility to and clinical manifestations of ankylosing spondylitis (AS), a common form of spondyloarthritis. Four hundred and sixty-eight individuals from 182 Finnish families affected with AS were studied. No association between individual IL10 promoter region polymorphisms or marker haplotype was observed with susceptibility to AS, but weak association was noted between the IL10-597 and -824 SNPs and age of disease onset (P= 0.01 for each SNP). The IL10.G4 allele was associated with BASFI (corrected for disease duration) (P= 0.03). We conclude that IL10 promoter polymorphisms have no significant effect on susceptibility to AS, but may play a minor role in determining age of disease onset and disease severity. © 2003 Nature Publishing Group All rights reserved.

Polymorphism in codon 17 of the CTLA-4 gene (+49 A/G) is not associated with susceptibility to rheumatoid arthritis in British Caucasians

The role of the CTLA-4 antigen in the development of autoimmune diseases is well documented, with several autoimmune disorders showing association or linkage with the CTLA-4 locus. Its role in the aetiology of rheumatoid arthritis (RA) however, remains unclear, as the functional studies of the B7-CTLA-4 pathway in mouse models of RA and genetic studies in humans have given contrasting results. We have studied the single nucleotide polymorphism at position +49 (A/G) of the CTLA-4 gene, in a cohort of 421 RA cases and 452 healthy controls from the UK. Despite the high statistical power to detect even a weak susceptibility effect, no significant association was found. We also analysed the distribution of the allele and genotype frequencies with respect to the presence of the shared epitope (a known RA susceptibility factor) and found no statistically significant differences. We conclude that, although the importance of the B7-CTLA-4 interaction in the development of RA can not be excluded, the CTLA-4 gene is unlikely to be a predisposing factor to this disease.

DC-SIGN and alpha 2,6-sialylated IgG Fc interaction is dispensable for the anti-inflammatory activity of IVIg on human dendritic cells

Intravenous immunoglobulin (IVIg) is widely used to treat autoimmune diseases. Several mutually nonexclusive mechanisms are proposed to explain the beneficial effects of IVIg in patients (1, 2). Lately, Ravetch and colleagues (3) demonstrate that anti-inflammatory activity of IVIg is mediated mainly by antibodies that contain terminal _2,6-sialic acid linkages at the Asn297-linked glycan of Fc region.

Of tolerance in mice and men : studies in APECED and Aire

Autoimmune diseases affect 5 % of the population and come in many forms, such as diabetes, rheumatoid arthritis and MS. However, how and why autoimmune diseases arise are not yet fully resolved. In this thesis, the onset of autoimmunity was investigated using both patient samples and a mouse model of autoimmunity. Autoimmune diseases are usually complex, due to a number of different causative genes and environmental factors. However, a few monogenic autoimmune diseases have been described, which are caused by mutations in only one gene per disease. One of such disease is called APECED (autoimmune polyendocrinopathy-candidiasis-ectodermal dystrophy) and is enriched in the Finnish population. The causative gene behind APECED is named AIRE from AutoImmune REgulator. How malfunction of just one gene product can cause the multitude of disease components found in APECED is not yet resolved. This thesis sought out to find out more about the functions of AIRE, in order to reveal why APECED and other autoimmune diseases arise and what goes wrong? Usually, immune cells are taught to distinguish between self and non-self during their development. That way, immune cells can fight off bacteria and microbes while leaving the tissues and organs of the host organism itself unharmed. In APECED, the development of immune cells called αβ T cells is incomplete. The cells are not able to fully distinguish between self and non-self. This leads to autodestruction of self tissues and autoimmune disease. One of the achievements of this thesis was the finding that the development of another set of T cells called γδ T cells is not affected by AIRE in mice or in men. Instead, we found that another type of immune cell important in tolerance, called the dendritic cell is defective in APECED patients and is not able to respond to microbial stimulus in a normal fashion. Finally, we studied Aire-deficient mice and found that autoantibodies expressed in the mice were not targeted against the same molecules as those found in APECED patients. This indicates differences in the autoimmune pathology in mice and men. More work is still required before we understand the mechanisms of tolerance and autoimmunity well enough to be able to cure APECED, let alone the more complex autoimmune diseases. Yet altogether, the findings of this thesis work bring us one step closer to finding out why and how APECED and common autoimmune diseases arise.

Dietary elecrolytes and cyclosporine toxicity. Experimental studies in spontaneuosly hypertensive rats

Cyclosporine-A (CsA) is widely used after organ transplantation to prevent rejection and in the treatment of autoimmune diseases. Hypertension and nephrotoxicity are common side-effects of CsA. Studies in patients on the prevention of the side-effects of CsA are difficult to conduct because the patients often receive a combination of different drugs thus making study of the side-effects of a single drug impossible. A challenge in experimental studies has been the lack of an animal model in which the side-effects concomitantly occur. Epidemiological data show an association between sodium (Na) intake and blood pressure. There is also evidence on low dietary intake of magnesium (Mg) and potassium (K) and high blood pressure. Our study was designed to develop an experimental model to study the side-effects of CsA in spontaneously hypertensive rats (SHR). On high dietary sodium, CsA caused hypertension, left ventricular hypertrophy (LVH), narrowing of the coronary arteries, small myocardial infarctions, and proteinuria, reduced creatinine clearance and histopathological renal injury in SHR. Loss of Mg into the urine caused by CsA resulted in Mg depletion in the tissues. Renal excretion of dopamine was reduced and the renin-angiotensin-aldosterone system was activated. We investigated the effects of dietary Mg and/or K and the calcium antagonist drug, isradipine, on the prevention of CsA toxicity. Dietary supplementation of Mg alone or in combination with K prevented from the deleterious pathophysiological and histopathological changes in the kidneys and the heart. K alone had little effect. Isradipine protected better than Mg from LVH, but the combination of isradipine and Mg was the most effective. Isradipine did not, however, protect against Mg loss. In our animal model, the combination of high dietary Na and treatment with CsA accelerated the development of the cardiovascular and renal changes clinically known as the side-effects of CsA. Dietary supplementation of Mg and K and reduction of Na intake and the calcium antagonist drug isradipine prevent from the deleterious effects of CsA.

Identification of diverse full-length endogenous betaretroviruses in megabats and microbats

Background: Betaretroviruses infect a wide range of species including primates, rodents, ruminants, and marsupials. They exist in both endogenous and exogenous forms and are implicated in animal diseases such as lung cancer in sheep, and in human disease, with members of the human endogenous retrovirus-K (HERV-K) group of endogenous betaretroviruses (βERVs) associated with human cancers and autoimmune diseases. To improve our understanding of betaretroviruses in an evolutionarily distinct host species, we characterized βERVs present in the genomes and transcriptomes of mega- and microbats, which are an important reservoir of emerging viruses.Results: A diverse range of full-length βERVs were discovered in mega- and microbat genomes and transcriptomes including the first identified intact endogenous retrovirus in a bat. Our analysis revealed that the genus Betaretrovirus can be divided into eight distinct sub-groups with evidence of cross-species transmission. Betaretroviruses are revealed to be a complex retrovirus group, within which one sub-group has evolved from complex to simple genomic organization through the acquisition of an env gene from the genus Gammaretrovirus. Molecular dating suggests that bats have contended with betaretroviral infections for over 30 million years.Conclusions: Our study reveals that a diverse range of betaretroviruses have circulated in bats for most of their evolutionary history, and cluster with extant betaretroviruses of divergent mammalian lineages suggesting that their distribution may be largely unrestricted by host species barriers. The presence of βERVs with the ability to transcribe active viral elements in a major animal reservoir for viral pathogens has potential implications for public health. © 2013 Hayward et al.; licensee BioMed Central Ltd.

Genetics of T cell co-stimulatory receptors -CD28, CTLA4, ICOS and PDCD1 in immunity and transplantation

Co-stimulatory signals are essential for the activation of naïve T cells and productive immune response. Naïve T cells receive first, antigen-specific signal through T cell receptor. Co-stimulatory receptors provide the second signal which can be either activating or inhibitory. The balance between signals determines the outcome of an immune response. CD28 is crucial for T cell activation; whereas cytotoxic T lymphocyte associated antigen 4 (CTLA4) mediates critical inhibitory signal. Inducible co-stimulator (ICOS) augments cytokine expression and plays role in immunoglobulin class switching. Programmed cell death 1 (PDCD1) acts as negative regulator of T cell proliferation and cytokine responses. The co-stimulatory receptor pathways are potentially involved in self-tolerance and thus, they provide a promising therapeutic strategy for autoimmune diseases and transplantation. The genes encoding CD28, CTLA4 and ICOS are located adjacently in the chromosome region 2q33. The PDCD1 gene maps further, to the region 2q37. CTLA4 and PDCD1 are associated with the risk of a few autoimmune diseases. There is strong linkage disequilibrium (LD) on the 2q33 region; the whole gene of CD28 exists in its own LD block but CTLA4 and the 5' part of ICOS are within a same LD block. The 3' part of ICOS and PDCD1 are in their own separate LD blocks. Extended haplotypes covering the 2q33 region can be identified. This study focuses on immune related conditions like coeliac disease (CD) which is a chronic inflammatory disease with autoimmune features. Immunoglobulin A deficiency (IgAD) belongs to the group of primary antibody deficiencies characterised by reduced levels of immunoglobulins. IgAD co-occurs often with coeliac disease. Renal transplantation is needed in the end stage kidney diseases. Transplantation causes strong immune response which is tried to suppress with drugs. All these conditions are multifactorial with complex genetic background and multiple environmental factors affecting the outcome. We have screened ICOS for polymorphisms by sequencing the exon regions. We detected 11 new variants and determined their frequencies in Finnish population. We have measured linkage disequilibrium on the 2q33 region in Finnish as well as other European populations and observed conserved haplotypes. We analysed genetic association and linkage of the co-stimulatory receptor gene region aiming to study if it is a common risk locus for immune diseases. The 2q33 region was replicated to be linked to coeliac disease in Finnish population and CTLA4-ICOS haplotypes were found to be associated with CD and IgAD being the first non-HLA risk locus common for CD and immunodeficiencies. We also showed association between ICOS and the outcome of kidney transplantation. Our results suggest new evidence for CTLA4-ICOS gene region to be involved in susceptibility of coeliac disease. The earlier published contradictory association results can be explained by involvement of both CTLA4 and ICOS in disease susceptibility. The pattern of variants acting together rather than a single polymorphism may confer the disease risk. These genes may predispose also to immunodeficiencies as well as decreased graft survival and delayed graft function. Consequently, the present study indicates that like the well established HLA locus, the co-stimulatory receptor genes predispose to variety of immune disorders.

Targeting antigen to diverse APCs inactivates memory CD8+ T cells without eliciting tissue-destructive effector function

Memory T cells develop early during the preclinical stages of autoimmune diseases and have traditionally been considered resistant to tolerance induction. As such, they may represent a potent barrier to the successful immunotherapy of established autoimmune diseases. It was recently shown that memory CD8+ T cell responses are terminated when Ag is genetically targeted to steady-state dendritic cells. However, under these conditions, inactivation of memory CD8+ T cells is slow, allowing transiently expanded memory CD8+ T cells to exert tissue-destructive effector function. In this study, we compared different Ag-targeting strategies and show, using an MHC class II promoter to drive Ag expression in a diverse range of APCs, that CD8+ memory T cells can be rapidly inactivated by MHC class II+ hematopoietic APCs through a mechanism that involves a rapid and sustained downregulation of TCR, in which the effector response of CD8+ memory cells is rapidly truncated and Ag-expressing target tissue destruction is prevented. Our data provide the first demonstration that genetically targeting Ag to a broad range of MHC class II+ APC types is a highly efficient way to terminate memory CD8+ T cell responses to prevent tissue-destructive effector function and potentially established autoimmune diseases. Copyright © 2010 by The American Association of Immunologists, Inc.

An automated framework for understanding structural variations in the binding grooves of MHC class II molecules

Background: MHC/HLA class II molecules are important components of the immune system and play a critical role in processes such as phagocytosis. Understanding peptide recognition properties of the hundreds of MHC class II alleles is essential to appreciate determinants of antigenicity and ultimately to predict epitopes. While there are several methods for epitope prediction, each differing in their success rates, there are no reports so far in the literature to systematically characterize the binding sites at the structural level and infer recognition profiles from them. Results: Here we report a new approach to compare the binding sites of MHC class II molecules using their three dimensional structures. We use a specifically tuned version of our recent algorithm, PocketMatch. We show that our methodology is useful for classification of MHC class II molecules based on similarities or differences among their binding sites. A new module has been used to define binding sites in MHC molecules. Comparison of binding sites of 103 MHC molecules, both at the whole groove and individual sub-pocket levels has been carried out, and their clustering patterns analyzed. While clusters largely agree with serotypic classification, deviations from it and several new insights are obtained from our study. We also present how differences in sub-pockets of molecules associated with a pair of autoimmune diseases, narcolepsy and rheumatoid arthritis, were captured by PocketMatch(13). Conclusion: The systematic framework for understanding structuralvariations in MHC class II molecules enables large scale comparison of binding grooves and sub-pockets, which is likely to have direct implications towards predicting epitopes and understanding peptide binding preferences.

An automated framework for understanding structural variations in the binding grooves of MHC class II molecules

Background: MHC/HLA class II molecules are important components of the immune system and play a critical role in processes such as phagocytosis. Understanding peptide recognition properties of the hundreds of MHC class II alleles is essential to appreciate determinants of antigenicity and ultimately to predict epitopes. While there are several methods for epitope prediction, each differing in their success rates, there are no reports so far in the literature to systematically characterize the binding sites at the structural level and infer recognition profiles from them. Results: Here we report a new approach to compare the binding sites of MHC class II molecules using their three dimensional structures. We use a specifically tuned version of our recent algorithm, PocketMatch. We show that our methodology is useful for classification of MHC class II molecules based on similarities or differences among their binding sites. A new module has been used to define binding sites in MHC molecules. Comparison of binding sites of 103 MHC molecules, both at the whole groove and individual sub-pocket levels has been carried out, and their clustering patterns analyzed. While clusters largely agree with serotypic classification, deviations from it and several new insights are obtained from our study. We also present how differences in sub-pockets of molecules associated with a pair of autoimmune diseases, narcolepsy and rheumatoid arthritis, were captured by PocketMatch(13). Conclusion: The systematic framework for understanding structural variations in MHC class II molecules enables large scale comparison of binding grooves and sub-pockets, which is likely to have direct implications towards predicting epitopes and understanding peptide binding preferences.