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.

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...

Gene expression profiling reveals a downregulation in immune-associated genes in patients with AS

Objective: To identify differentially expressed genes in peripheral blood mononuclear cells (PBMCs) from patients with ankylosing spondylitis (AS) compared with healthy individuals. Methods: RNA was extracted from PBMCs collected from 18 patients with active disease and 18 gender-matched and age-matched controls. Expression profiles of these cells were determined using microarray. Candidate genes with differential expressions were confirmed in the same samples using quantitative reverse transcription-PCR (qRT-PCR). These genes were then validated in a different sample cohort of 35 patients with AS and 18 controls by qRT-PCR. Results: Microarray analysis identified 452 genes detected with 485 probes which were differentially expressed between patients with AS and controls. Underexpression of NR4A2, tumour necrosis factor AIP3 (TNFAIP3) and CD69 was confirmed. These genes were further validated in a different sample group in which the patients with AS had a wider range of disease activity. Predictive algorithms were also developed from the expression data using receiver-operating characteristic curves, which demonstrated that the three candidate genes have ∼80% power to predict AS according to their expression levels. Conclusions: The findings show differences in global gene expression patterns between patients with AS and controls, suggesting an immunosuppressive phenotype in the patients. Furthermore, downregulated expression of three immune-related genes was confirmed. These candidate genes were also shown to be strong predictive markers for 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.

Enrichment of circulating interleukin-17-secreting interleukin-23 receptor-positive γ/δ T cells in patients with active ankylosing spondylitis

Objective Ankylosing spondylitis (AS) is a common inflammatory arthritis affecting primarily the axial skeleton. IL23R is genetically associated with AS. This study was undertaken to investigate and characterize the role of interleukin-23 (IL-23) signaling in AS pathogenesis. Methods The study population consisted of patients with active AS (n = 17), patients with psoriatic arthritis (n = 8), patients with rheumatoid arthritis, (n = 9), and healthy subjects (n = 20). IL-23 receptor (IL-23R) expression in T cells was determined in each subject group, and expression levels were compared. Results The proportion of IL-23R-expressing T cells in the periphery was 2-fold higher in AS patients than in healthy controls, specifically driven by a 3-fold increase in IL-23R-positive γ/δ T cells in AS patients. The proportions of CD4+ and CD8+ cells that were positive for IL-17 were unchanged. This increased IL-23R expression on γ/δ T cells was also associated with enhanced IL-17 secretion, with no observable IL-17 production from IL-23R-negative γ/δ T cells in AS patients. Furthermore, γ/δ T cells from AS patients were heavily skewed toward IL-17 production in response to stimulation with IL-23 and/or anti-CD3/CD28. Conclusion Recently, mouse models have shown IL-17-secreting γ/δ T cells to be pathogenic in infection and autoimmunity. Our data provide the first description of a potentially pathogenic role of these cells in a human autoimmune disease. Since IL-23 is a maturation and growth factor for IL-17-producing cells, increased IL-23R expression may regulate the function of this putative pathogenic γ/δ T cell population.

The effect of HLA-DR on susceptibility to rheumatoid arthritis is influenced by the associated lymphotoxin α-tumor necrosis factor haplotype

Objective. HLA-DRB1, a major genetic determinant of susceptibility to rheumatoid arthritis (RA), is located within 1,000 kb of the gene encoding tumor necrosis factor (TNF). Because certain HLA-DRB1*04 subtypes increase susceptibility to RA, investigation of the role of the TNF gene is complicated by linkage disequilibrium (LD) between TNF and DRB1 alleles. By adequately controlling for this LD, we aimed to investigate the presence of additional major histocompatibility complex (MHC) susceptibility genes. Methods. We identified 274 HLA-DRB1*04-positive cases of RA and 271 HLA-DRB1*04-positive population controls. Each subject was typed for 6 single-nucleotide polymorphisms within a 4.5-kb region encompassing TNF and lymphotoxin a (LTA). LTA-TNF haplotypes in these unrelated individuals were determined using a combination of family data and the PHASE software program. Results. Significant differences in LTA-TNF haplotype frequencies were observed between different subtypes of HLA-DRB1*04. The LTA-TNF haplotypes observed were very restricted, with only 4 haplotypes constituting 81% of all haplotypes present. Among individuals carrying DRB1*0401, the LTA-TNF 2 haplotype was significantly underrepresented in cases compared with controls (odds ratio 0.5 [95% confidence interval 0.3-0.8], P = 0.007), while in those with DRB1*0404, the opposite effect was observed (P = 0.007). Conclusion. These findings suggest that the MHC contains genetic elements outside the LTA-TNF region that modify the effect of HLA-DRB1 on susceptibility to RA.

Serum levels of soluble receptor for advanced glycation end products and of S100 proteins are associated with inflammatory, autoantibody, and classical risk markers of joint and vascular damage in rheumatoid arthritis

Introduction: The receptor for advanced glycation end products (RAGE) is a member of the immunoglobulin superfamily of cell surface receptor molecules. High concentrations of three of its putative proinflammatory ligands, S100A8/A9 complex (calprotectin), S100A8, and S100A12, are found in rheumatoid arthritis (RA) serum and synovial fluid. In contrast, soluble RAGE (sRAGE) may prevent proinflammatory effects by acting as a decoy. This study evaluated the serum levels of S100A9, S100A8, S100A12 and sRAGE in RA patients, to determine their relationship to inflammation and joint and vascular damage. Methods: Serum sRAGE, S100A9, S100A8 and S100A12 levels from 138 patients with established RA and 44 healthy controls were measured by ELISA and compared by unpaired t test. In RA patients, associations with disease activity and severity variables were analyzed by simple and multiple linear regressions. Results: Serum S100A9, S100A8 and S100A12 levels were correlated in RA patients. S100A9 levels were associated with body mass index (BMI), and with serum levels of S100A8 and S100A12. S100A8 levels were associated with serum levels of S100A9, presence of anti-citrullinated peptide antibodies (ACPA), and rheumatoid factor (RF). S100A12 levels were associated with presence of ACPA, history of diabetes, and serum S100A9 levels. sRAGE levels were negatively associated with serum levels of C-reactive protein (CRP) and high-density lipoprotein (HDL), history of vasculitis, and the presence of the RAGE 82Ser polymorphism. Conclusions: sRAGE and S100 proteins were associated not just with RA inflammation and autoantibody production, but also with classical vascular risk factors for end-organ damage. Consistent with its role as a RAGE decoy molecule, sRAGE had the opposite effects to S100 proteins in that S100 proteins were associated with autoantibodies and vascular risk, whereas sRAGE was associated with protection against joint and vascular damage. These data suggest that RAGE activity influences co-development of joint and vascular disease in rheumatoid arthritis patients.

Cardiovascular disease is increased prior to onset of rheumatoid arthritis but not osteoarthritis: The population-based Nord-Trøndelag health study (HUNT)

Introduction: Patients with rheumatoid arthritis (RA) have increased risk of cardiovascular (CV) events. We sought to test the hypothesis that due to increased inflammation, CV disease and risk factors are associated with increased risk of future RA development. Methods: The population-based Nord-Trøndelag health surveys (HUNT) were conducted among the entire adult population of Nord-Trøndelag, Norway. All inhabitants 20 years or older were invited, and information was collected through comprehensive questionnaires, a clinical examination, and blood samples. In a cohort design, data from HUNT2 (1995-1997, baseline) and HUNT3 (2006-2008, follow-up) were obtained to study participants with RA (n = 786) or osteoarthritis (n = 3,586) at HUNT3 alone, in comparison with individuals without RA or osteoarthritis at both times (n = 33,567). Results: Female gender, age, smoking, body mass index, and history of previous CV disease were associated with self-reported incident RA (previous CV disease: odds ratio 1.52 (95% confidence interval 1.11-2.07). The findings regarding previous CV disease were confirmed in sensitivity analyses excluding participants with psoriasis (odds ratio (OR) 1.70 (1.23-2.36)) or restricting the analysis to cases with a hospital diagnosis of RA (OR 1.90 (1.10-3.27)) or carriers of the shared epitope (OR 1.76 (1.13-2.74)). History of previous CV disease was not associated with increased risk of osteoarthritis (OR 1.04 (0.86-1.27)). Conclusion: A history of previous CV disease was associated with increased risk of incident RA but not osteoarthritis.

Unique and cross-reactive T cell epitope peptides of the major bahia grass pollen allergen, pas n 1

Background Bahia grass pollen (BaGP) is a major cause of allergic rhinitis. Subcutaneous allergen-specific immunotherapy is effective for grass pollen allergy, but is unsuitable for patients with moderate to severe asthma due to the risk of anaphylaxis. T cell-reactive but IgE nonreactive peptides provide a safer treatment option. This study aimed to identify and characterize dominant CD4+ T cell epitope peptides of the major BaGP allergen, Pas n 1. Methods Pas n 1-specific T cell lines generated from the peripheral blood of BaGP-allergic subjects were tested for proliferative and cytokine response to overlapping 20-mer Pas n 1 peptides. Cross-reactivity to homologous peptides from Lol p 1 and Cyn d 1 of Ryegrass and Bermuda grass pollen, respectively, was assessed using Pas n 1 peptide-specific T cell clones. MHC class II restriction of Pas n 1 peptide T cell recognition was determined by HLA blocking assays and peptide IgE reactivity tested by dot blotting. Results Three Pas n 1 peptides showed dominant T cell reactivity; 15 of 18 (83%) patients responded to one or more of these peptides. T cell clones specific for dominant Pas n 1 peptides showed evidence of species-specific T cell reactivity as well as cross-reactivity with other group 1 grass pollen allergens. The dominant Pas n 1 T cell epitope peptides showed HLA binding diversity and were non-IgE reactive. Conclusions The immunodominant T cell-reactive Pas n 1 peptides are candidates for safe immunotherapy for individuals, including those with asthma, who are allergic to Bahia and possibly other grass pollens.

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.

Interleukin 17A is an immune marker for chlamydial disease severity and pathogenesis in the koala (Phascolarctos cinereus)

The koala (Phascolarctos cinereus) is an iconic Australian marsupial species that is facing many threats to its survival. Chlamydia pecorum infections are a significant contributor to this ongoing decline. A major limiting factor in our ability to manage and control chlamydial disease in koalas is a limited understanding of the koala’s cell-mediated immune response to infections by this bacterial pathogen. To identify immunological markers associated with chlamydial infection and disease in koalas, we used koala-specific Quantitative Real Time PCR (qrtPCR) assays to profile the cytokine responses of Peripheral Blood Mononuclear Cells (PBMCs) collected from 41 koalas with different stages of chlamydial disease. Target cytokines included the principal Th1 (Interferon gamma; IFNγ), Th2 (Interleukin 10; IL10), and pro-inflammatory cytokines (Tumor Necrosis Factor alpha; TNFα). A novel koala-specific IL17A qrtPCR assay was also developed as part of this study to quantitate the gene expression of this Th17 cytokine in koalas. A statistically significant higher IL17A gene expression was observed in animals with current chlamydial disease compared to animals with asymptomatic chlamydial infection. A modest up-regulation of pro-inflammatory cytokines, such as TNFα and IFNγ, was also observed in these animals with signs of current chlamydial disease. IL10 gene expression was not evident in the majority of animals from both groups. Future longitudinal studies are now required to confirm the role played by cytokines in pathology and/or protection against C. pecorum infection in the koala.

A comprehensive ovine model of blood transfusion

Background The growing awareness of transfusion-associated morbidity and mortality necessitates investigations into the underlying mechanisms. Small animals have been the dominant transfusion model but have associated limitations. This study aimed to develop a comprehensive large animal (ovine) model of transfusion encompassing: blood collection, processing and storage, compatibility testing right through to post-transfusion outcomes. Materials and methods Two units of blood were collected from each of 12 adult male Merino sheep and processed into 24 ovine-packed red blood cell (PRBC) units. Baseline haematological parameters of ovine blood and PRBC cells were analysed. Biochemical changes in ovine PRBCs were characterized during the 42-day storage period. Immunological compatibility of the blood was confirmed with sera from potential recipient sheep, using a saline and albumin agglutination cross-match. Following confirmation of compatibility, each recipient sheep (n = 12) was transfused with two units of ovine PRBC. Results Procedures for collecting, processing, cross-matching and transfusing ovine blood were established. Although ovine red blood cells are smaller and higher in number, their mean cell haemoglobin concentration is similar to human red blood cells. Ovine PRBC showed improved storage properties in saline–adenine–glucose–mannitol (SAG-M) compared with previous human PRBC studies. Seventy-six compatibility tests were performed and 17·1% were incompatible. Only cross-match compatible ovine PRBC were transfused and no adverse reactions were observed. Conclusion These findings demonstrate the utility of the ovine model for future blood transfusion studies and highlight the importance of compatibility testing in animal models involving homologous transfusions.

Flow evaluation and hemolysis analysis of BVAD centrifugal blood pump by computational fluids dynamics

Computational fluid dynamics (CFD) and particle image velocimetry (PIV) are commonly used techniques to evaluate the flow characteristics in the development stage of blood pumps. CFD technique allows rapid change to pump parameters to optimize the pump performance without having to construct a costly prototype model. These techniques are used in the construction of a bi-ventricular assist device (BVAD) which combines the functions of LVAD and RVAD in a compact unit. The BVAD construction consists of two separate chambers with similar impellers, volutes, inlet and output sections. To achieve the required flow characteristics of an average flow rate of 5 l/min and different pressure heads (left – 100mmHg and right – 20mmHg), the impellers were set at different rotating speeds. From the CFD results, a six-blade impeller design was adopted for the development of the BVAD. It was also observed that the fluid can flow smoothly through the pump with minimum shear stress and area of stagnation which are related to haemolysis and thrombosis. Based on the compatible Reynolds number the flow through the model was calculated for the left and the right pumps. As it was not possible to have both the left and right chambers in the experimental model, the left and right pumps were tested separately.

Patient reported adherence to hypertension treatment: A revalidation study

Background Adherence to hypertension management in patients with hypertension is known to influence their blood pressure control. It is important to measure patients’ adherence behaviours to assist with designing appropriate interventions to improve blood pressure control. Aims The purposes of this study were to use confirmatory factor analysis to revalidate the Therapeutic Adherence Scale for Hypertensive Patients (TASHP), and to calculate the cut-off score for classifying adherence behaviours into two groups: satisfactory and low adherence behaviours. Methods Systematic random sampling was used to recruit patients with hypertension in China. Demographic characteristics, the TASHP and blood pressure were collected. The psychometric tests of the TASHP included: construct validity, criteria-related validity, internal reliability, and split-half reliability. The area under the receiver operating characteristics curve and Youden index were used to identify the cut-off score of the TASHP for blood pressure control. Results This study involved 366 patients. Confirmatory factor analysis supported the four-component structure of the TASHP proposed in the original scale development study. The TASHP has a satisfactory internal reliability (Cronbach’s α > 0.7) and a satisfactory split-half reliability (Spearman–Brown coefficients > 0.7). The patients with overall scores of the TASHP ⩾ 109 points were considered to have satisfactory adherence behaviours. Conclusion The TASHP is a validated and reliable instrument to measure the adherence to hypertension management in Chinese patients with hypertension. The cut-off score of 109 points can be considered as an effective measure to classify the level of adherence into satisfactory and low adherence behaviours.

Cardiometabolic effects of genetic upregulation of the interleukin 1 receptor antagonist: A Mendelian randomisation analysis

Background To investigate potential cardiovascular and other effects of long-term pharmacological interleukin 1 (IL-1) inhibition, we studied genetic variants that produce inhibition of IL-1, a master regulator of inflammation. Methods We created a genetic score combining the effects of alleles of two common variants (rs6743376 and rs1542176) that are located upstream of IL1RN, the gene encoding the IL-1 receptor antagonist (IL-1Ra; an endogenous inhibitor of both IL-1α and IL-1β); both alleles increase soluble IL-1Ra protein concentration. We compared effects on inflammation biomarkers of this genetic score with those of anakinra, the recombinant form of IL-1Ra, which has previously been studied in randomised trials of rheumatoid arthritis and other inflammatory disorders. In primary analyses, we investigated the score in relation to rheumatoid arthritis and four cardiometabolic diseases (type 2 diabetes, coronary heart disease, ischaemic stroke, and abdominal aortic aneurysm; 453 411 total participants). In exploratory analyses, we studied the relation of the score to many disease traits and to 24 other disorders of proposed relevance to IL-1 signalling (746 171 total participants). Findings For each IL1RN minor allele inherited, serum concentrations of IL-1Ra increased by 0·22 SD (95% CI 0·18–0·25; 12·5%; p=9·3 × 10−33), concentrations of interleukin 6 decreased by 0·02 SD (−0·04 to −0·01; −1·7%; p=3·5 × 10−3), and concentrations of C-reactive protein decreased by 0·03 SD (−0·04 to −0·02; −3·4%; p=7·7 × 10−14). We noted the effects of the genetic score on these inflammation biomarkers to be directionally concordant with those of anakinra. The allele count of the genetic score had roughly log-linear, dose-dependent associations with both IL-1Ra concentration and risk of coronary heart disease. For people who carried four IL-1Ra-raising alleles, the odds ratio for coronary heart disease was 1·15 (1·08–1·22; p=1·8 × 10−6) compared with people who carried no IL-1Ra-raising alleles; the per-allele odds ratio for coronary heart disease was 1·03 (1·02–1·04; p=3·9 × 10−10). Per-allele odds ratios were 0·97 (0·95–0·99; p=9·9 × 10−4) for rheumatoid arthritis, 0·99 (0·97–1·01; p=0·47) for type 2 diabetes, 1·00 (0·98–1·02; p=0·92) for ischaemic stroke, and 1·08 (1·04–1·12; p=1·8 × 10−5) for abdominal aortic aneurysm. In exploratory analyses, we observed per-allele increases in concentrations of proatherogenic lipids, including LDL-cholesterol, but no clear evidence of association for blood pressure, glycaemic traits, or any of the 24 other disorders studied. Modelling suggested that the observed increase in LDL-cholesterol could account for about a third of the association observed between the genetic score and increased coronary risk. Interpretation Human genetic data suggest that long-term dual IL-1α/β inhibition could increase cardiovascular risk and, conversely, reduce the risk of development of rheumatoid arthritis. The cardiovascular risk might, in part, be mediated through an increase in proatherogenic lipid concentrations. Funding UK Medical Research Council, British Heart Foundation, UK National Institute for Health Research, National Institute for Health Research Cambridge Biomedical Research Centre, European Research Council, and European Commission Framework Programme 7.