Genome-wide association study of antidepressant treatment resistance in a population-based cohort using health service prescription data and meta-analysis with GENDEP.

Antidepressants demonstrate modest response rates in the treatment of major depressive disorder (MDD). Despite previous genome-wide association studies (GWAS) of antidepressant treatment response, the underlying genetic factors are unknown. Using prescription data in a population and family-based cohort (Generation Scotland: Scottish Family Health Study; GS:SFHS), we sought to define a measure of (a) antidepressant treatment resistance and (b) stages of antidepressant resistance by inferring antidepressant switching as non-response to treatment. GWAS were conducted separately for antidepressant treatment resistance in GS:SFHS and the Genome-based Therapeutic Drugs for Depression (GENDEP) study and then meta-analysed (meta-analysis n = 4213, cases = 358). For stages of antidepressant resistance, a GWAS on GS:SFHS only was performed (n = 3452). Additionally, we conducted gene-set enrichment, polygenic risk scoring (PRS) and genetic correlation analysis. We did not identify any significant loci, genes or gene sets associated with antidepressant treatment resistance or stages of resistance. Significant positive genetic correlations of antidepressant treatment resistance and stages of resistance with neuroticism, psychological distress, schizotypy and mood disorder traits were identified. These findings suggest that larger sample sizes are needed to identify the genetic architecture of antidepressant treatment response, and that population-based observational studies may provide a tractable approach to achieving the necessary statistical power.

Quantifying the extent to which index event biases influence large genetic association studies.

As genetic association studies increase in size to 100 000s of individuals, subtle biases may influence conclusions. One possible bias is 'index event bias' (IEB) that appears due to the stratification by, or enrichment for, disease status when testing associations between genetic variants and a disease-associated trait. We aimed to test the extent to which IEB influences some known trait associations in a range of study designs and provide a statistical framework for assessing future associations. Analyzing data from 113 203 non-diabetic UK Biobank participants, we observed three (near TCF7L2, CDKN2AB and CDKAL1) overestimated (body mass index (BMI) decreasing) and one (near MTNR1B) underestimated (BMI increasing) associations among 11 type 2 diabetes risk alleles (at P < 0.05). IEB became even stronger when we tested a type 2 diabetes genetic risk score composed of these 11 variants (-0.010 standard deviations BMI per allele, P = 5 × 10- 4), which was confirmed in four additional independent studies. Similar results emerged when examining the effect of blood pressure increasing alleles on BMI in normotensive UK Biobank samples. Furthermore, we demonstrated that, under realistic scenarios, common disease alleles would become associated at P < 5 × 10- 8 with disease-related traits through IEB alone, if disease prevalence in the sample differs appreciably from the background population prevalence. For example, some hypertension and type 2 diabetes alleles will be associated with BMI in sample sizes of >500 000 if the prevalence of those diseases differs by >10% from the background population. In conclusion, IEB may result in false positive or negative genetic associations in very large studies stratified or strongly enriched for/against disease cases.

Meta-analysis of 49†549 individuals imputed with the 1000 Genomes Project reveals an exonic damaging variant in ANGPTL4 determining fasting TG levels.

BACKGROUND: So far, more than 170 loci have been associated with circulating lipid levels through genome-wide association studies (GWAS). These associations are largely driven by common variants, their function is often not known, and many are likely to be markers for the causal variants. In this study we aimed to identify more new rare and low-frequency functional variants associated with circulating lipid levels. METHODS: We used the 1000 Genomes Project as a reference panel for the imputations of GWAS data from ∼60 000 individuals in the discovery stage and ∼90 000 samples in the replication stage. RESULTS: Our study resulted in the identification of five new associations with circulating lipid levels at four loci. All four loci are within genes that can be linked biologically to lipid metabolism. One of the variants, rs116843064, is a damaging missense variant within the ANGPTL4 gene. CONCLUSIONS: This study illustrates that GWAS with high-scale imputation may still help us unravel the biological mechanism behind circulating lipid levels.

Investigating shared aetiology between type 2 diabetes and major depressive disorder in a population based cohort.

Type II diabetes (T2D) and major depressive disorder (MDD) are often co-morbid. The reasons for this co-morbidity are unclear. Some studies have highlighted the importance of environmental factors and a causal relationship between T2D and MDD has also been postulated. In the present study we set out to investigate the shared aetiology between T2D and MDD using Mendelian randomization in a population based sample, Generation Scotland: the Scottish Family Health Study (N = 21,516). Eleven SNPs found to be associated with T2D were tested for association with MDD and psychological distress (General Health Questionnaire scores). We also assessed causality and genetic overlap between T2D and MDD using polygenic risk scores (PRS) assembled from the largest available GWAS summary statistics to date. No single T2D risk SNP was associated with MDD in the MR analyses and we did not find consistent evidence of genetic overlap between MDD and T2D in the PRS analyses. Linkage disequilibrium score regression analyses supported these findings as no genetic correlation was observed between T2D and MDD (rG = 0.0278 (S.E. 0.11), P-value = 0.79). As suggested by previous studies, T2D and MDD covariance may be better explained by environmental factors. Future studies would benefit from analyses in larger cohorts where stratifying by sex and looking more closely at MDD cases demonstrating metabolic dysregulation is possible. © 2016 The Authors. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics Published by Wiley Periodicals, Inc.

Genetic and Environmental Risk for Chronic Pain and the Contribution of Risk Variants for Major Depressive Disorder: A Family-Based Mixed-Model Analysis.

BACKGROUND: Chronic pain is highly prevalent and a significant source of disability, yet its genetic and environmental risk factors are poorly understood. Its relationship with major depressive disorder (MDD) is of particular importance. We sought to test the contribution of genetic factors and shared and unique environment to risk of chronic pain and its correlation with MDD in Generation Scotland: Scottish Family Health Study (GS:SFHS). We then sought to replicate any significant findings in the United Kingdom Biobank study. METHODS AND FINDINGS: Using family-based mixed-model analyses, we examined the contribution of genetics and shared family environment to chronic pain by spouse, sibling, and household relationships. These analyses were conducted in GS:SFHS (n = 23,960), a family- and population-based study of individuals recruited from the Scottish population through their general practitioners. We then examined and partitioned the correlation between chronic pain and MDD and estimated the contribution of genetic factors and shared environment in GS:SFHS. Finally, we used data from two independent genome-wide association studies to test whether chronic pain has a polygenic architecture and examine whether genomic risk of psychiatric disorder predicted chronic pain and whether genomic risk of chronic pain predicted MDD. These analyses were conducted in GS:SFHS and repeated in UK Biobank, a study of 500,000 from the UK population, of whom 112,151 had genotyping and phenotypic data. Chronic pain is a moderately heritable trait (heritability = 38.4%, 95% CI 33.6% to 43.9%) that is significantly concordant in spouses (variance explained 18.7%, 95% CI 9.5% to 25.1%). Chronic pain is positively correlated with depression (ρ = 0.13, 95% CI 0.11 to 0.15, p = 2.72x10-68) and shows a tendency to cluster within families for genetic reasons (genetic correlation = 0.51, 95%CI 0.40 to 0.62, p = 8.24x10-19). Polygenic risk profiles for pain, generated using independent GWAS data, were associated with chronic pain in both GS:SFHS (maximum ß = 6.18x10-2, 95% CI 2.84 x10-2 to 9.35 x10-2, p = 4.3x10-4) and UK Biobank (maximum ß = 5.68 x 10-2, 95% CI 4.70x10-2 to 6.65x10-2, p < 3x10-4). Genomic risk of MDD is also significantly associated with chronic pain in both GS:SFHS (maximum ß = 6.62x10-2, 95% CI 2.82 x10-2 to 9.76 x10-2, p = 4.3x10-4) and UK Biobank (maximum ß = 2.56x10-2, 95% CI 1.62x10-2 to 3.63x10-2, p < 3x10-4). Limitations of the current study include the possibility that spouse effects may be due to assortative mating and the relatively small polygenic risk score effect sizes. CONCLUSIONS: Genetic factors, as well as chronic pain in a partner or spouse, contribute substantially to the risk of chronic pain for an individual. Chronic pain is genetically correlated with MDD, has a polygenic architecture, and is associated with polygenic risk of MDD.

Exome-wide analysis of rare coding variation identifies novel associations with COPD and airflow limitation in MOCS3, IFIT3 and SERPINA12.

BACKGROUND: Several regions of the genome have shown to be associated with COPD in genome-wide association studies of common variants. OBJECTIVE: To determine rare and potentially functional single nucleotide polymorphisms (SNPs) associated with the risk of COPD and severity of airflow limitation. METHODS: 3226 current or former smokers of European ancestry with lung function measures indicative of Global Initiative for Chronic Obstructive Lung Disease (GOLD) 2 COPD or worse were genotyped using an exome array. An analysis of risk of COPD was carried out using ever smoking controls (n=4784). Associations with %predicted FEV1 were tested in cases. We followed-up signals of interest (p<10(-5)) in independent samples from a subset of the UK Biobank population and also undertook a more powerful discovery study by meta-analysing the exome array data and UK Biobank data for variants represented on both arrays. RESULTS: Among the associated variants were two in regions previously unreported for COPD; a low frequency non-synonymous SNP in MOCS3 (rs7269297, pdiscovery=3.08×10(-6), preplication=0.019) and a rare SNP in IFIT3, which emerged in the meta-analysis (rs140549288, pmeta=8.56×10(-6)). In the meta-analysis of % predicted FEV1 in cases, the strongest association was shown for a splice variant in a previously unreported region, SERPINA12 (rs140198372, pmeta=5.72×10(-6)). We also confirmed previously reported associations with COPD risk at MMP12, HHIP, GPR126 and CHRNA5. No associations in novel regions reached a stringent exome-wide significance threshold (p<3.7×10(-7)). CONCLUSIONS: This study identified several associations with the risk of COPD and severity of airflow limitation, including novel regions MOCS3, IFIT3 and SERPINA12, which warrant further study.

MICL controls inflammation in rheumatoid arthritis.

BACKGROUND: Myeloid inhibitory C-type lectin-like receptor (MICL, Clec12A) is a C-type lectin receptor (CLR) expressed predominantly by myeloid cells. Previous studies have suggested that MICL is involved in controlling inflammation. OBJECTIVE: To determine the role of this CLR in inflammatory pathology using Clec12A(-/-) mice. METHODS: Clec12A(-/-) mice were generated commercially and primarily characterised using the collagen antibody-induced arthritis (CAIA) model. Mechanisms and progress of disease were characterised by clinical scoring, histology, flow cytometry, irradiation bone-marrow chimera generation, administration of blocking antibodies and in vivo imaging. Characterisation of MICL in patients with rheumatoid arthritis (RA) was determined by immunohistochemistry and single nucleotide polymorphism analysis. Anti-MICL antibodies were detected in patient serum by ELISA and dot-blot analysis. RESULTS: MICL-deficient animals did not present with pan-immune dysfunction, but exhibited markedly exacerbated inflammation during CAIA, owing to the inappropriate activation of myeloid cells. Polymorphisms of MICL were not associated with disease in patients with RA, but this CLR was the target of autoantibodies in a subset of patients with RA. In wild-type mice the administration of such antibodies recapitulated the Clec12A(-/-) phenotype. CONCLUSIONS: MICL plays an essential role in regulating inflammation during arthritis and is an autoantigen in a subset of patients with RA. These data suggest an entirely new mechanism underlying RA pathogenesis, whereby the threshold of myeloid cell activation can be modulated by autoantibodies that bind to cell membrane-expressed inhibitory receptors.

Polygenic risk for alcohol dependence associates with alcohol consumption, cognitive function and social deprivation in a population-based cohort.

Alcohol dependence is frequently co-morbid with cognitive impairment. The relationship between these traits is complex as cognitive dysfunction may arise as a consequence of heavy drinking or exist prior to the onset of dependence. In the present study, we tested the genetic overlap between cognitive abilities and alcohol dependence using polygenic risk scores (PGRS). We created two independent PGRS derived from two recent genome-wide association studies (GWAS) of alcohol dependence (SAGE GWAS: n = 2750; Yale-Penn GWAS: n = 2377) in a population-based cohort, Generation Scotland: Scottish Family Health Study (GS:SFHS) (n = 9863). Data on alcohol consumption and four tests of cognitive function [Mill Hill Vocabulary (MHV), digit symbol coding, phonemic verbal fluency (VF) and logical memory] were available. PGRS for alcohol dependence were negatively associated with two measures of cognitive function: MHV (SAGE: P = 0.009, ß = -0.027; Yale-Penn: P = 0.001, ß = -0.034) and VF (SAGE: P = 0.0008, ß = -0.036; Yale-Penn: P = 0.00005, ß = -0.044). VF remained robustly associated after adjustment for education and social deprivation; however, the association with MHV was substantially attenuated. Shared genetic variants may account for some of the phenotypic association between cognitive ability and alcohol dependence. A significant negative association between PGRS and social deprivation was found (SAGE: P = 5.2 × 10(-7) , ß = -0.054; Yale-Penn: P = 0.000012, ß = -0.047). Individuals living in socially deprived regions were found to carry more alcohol dependence risk alleles which may contribute to the increased prevalence of problem drinking in regions of deprivation. Future work to identify genes which affect both cognitive impairment and alcohol dependence will help elucidate biological processes common to both disorders.

Predicting the risk of rheumatoid arthritis and its age of onset through modelling genetic risk variants with smoking.

The improved characterisation of risk factors for rheumatoid arthritis (RA) suggests they could be combined to identify individuals at increased disease risks in whom preventive strategies may be evaluated. We aimed to develop an RA prediction model capable of generating clinically relevant predictive data and to determine if it better predicted younger onset RA (YORA). Our novel modelling approach combined odds ratios for 15 four-digit/10 two-digit HLA-DRB1 alleles, 31 single nucleotide polymorphisms (SNPs) and ever-smoking status in males to determine risk using computer simulation and confidence interval based risk categorisation. Only males were evaluated in our models incorporating smoking as ever-smoking is a significant risk factor for RA in men but not women. We developed multiple models to evaluate each risk factor's impact on prediction. Each model's ability to discriminate anti-citrullinated protein antibody (ACPA)-positive RA from controls was evaluated in two cohorts: Wellcome Trust Case Control Consortium (WTCCC: 1,516 cases; 1,647 controls); UK RA Genetics Group Consortium (UKRAGG: 2,623 cases; 1,500 controls). HLA and smoking provided strongest prediction with good discrimination evidenced by an HLA-smoking model area under the curve (AUC) value of 0.813 in both WTCCC and UKRAGG. SNPs provided minimal prediction (AUC 0.660 WTCCC/0.617 UKRAGG). Whilst high individual risks were identified, with some cases having estimated lifetime risks of 86%, only a minority overall had substantially increased odds for RA. High risks from the HLA model were associated with YORA (P<0.0001); ever-smoking associated with older onset disease. This latter finding suggests smoking's impact on RA risk manifests later in life. Our modelling demonstrates that combining risk factors provides clinically informative RA prediction; additionally HLA and smoking status can be used to predict the risk of younger and older onset RA, respectively.

A genome-wide copy number association study of osteoporotic fractures points to the 6p25.1 locus.

BACKGROUND: Osteoporosis is a systemic skeletal disease characterised by reduced bone mineral density and increased susceptibility to fracture; these traits are highly heritable. Both common and rare copy number variants (CNVs) potentially affect the function of genes and may influence disease risk. AIM: To identify CNVs associated with osteoporotic bone fracture risk. METHOD: We performed a genome-wide CNV association study in 5178 individuals from a prospective cohort in the Netherlands, including 809 osteoporotic fracture cases, and performed in silico lookups and de novo genotyping to replicate in several independent studies. RESULTS: A rare (population prevalence 0.14%, 95% CI 0.03% to 0.24%) 210 kb deletion located on chromosome 6p25.1 was associated with the risk of fracture (OR 32.58, 95% CI 3.95 to 1488.89; p = 8.69 × 10(-5)). We performed an in silico meta-analysis in four studies with CNV microarray data and the association with fracture risk was replicated (OR 3.11, 95% CI 1.01 to 8.22; p = 0.02). The prevalence of this deletion showed geographic diversity, being absent in additional samples from Australia, Canada, Poland, Iceland, Denmark, and Sweden, but present in the Netherlands (0.34%), Spain (0.33%), USA (0.23%), England (0.15%), Scotland (0.10%), and Ireland (0.06%), with insufficient evidence for association with fracture risk. CONCLUSIONS: These results suggest that deletions in the 6p25.1 locus may predispose to higher risk of fracture in a subset of populations of European origin; larger and geographically restricted studies will be needed to confirm this regional association. This is a first step towards the evaluation of the role of rare CNVs in osteoporosis.

Informed conditioning on clinical covariates increases power in case-control association studies.

Genetic case-control association studies often include data on clinical covariates, such as body mass index (BMI), smoking status, or age, that may modify the underlying genetic risk of case or control samples. For example, in type 2 diabetes, odds ratios for established variants estimated from low-BMI cases are larger than those estimated from high-BMI cases. An unanswered question is how to use this information to maximize statistical power in case-control studies that ascertain individuals on the basis of phenotype (case-control ascertainment) or phenotype and clinical covariates (case-control-covariate ascertainment). While current approaches improve power in studies with random ascertainment, they often lose power under case-control ascertainment and fail to capture available power increases under case-control-covariate ascertainment. We show that an informed conditioning approach, based on the liability threshold model with parameters informed by external epidemiological information, fully accounts for disease prevalence and non-random ascertainment of phenotype as well as covariates and provides a substantial increase in power while maintaining a properly controlled false-positive rate. Our method outperforms standard case-control association tests with or without covariates, tests of gene x covariate interaction, and previously proposed tests for dealing with covariates in ascertained data, with especially large improvements in the case of case-control-covariate ascertainment. We investigate empirical case-control studies of type 2 diabetes, prostate cancer, lung cancer, breast cancer, rheumatoid arthritis, age-related macular degeneration, and end-stage kidney disease over a total of 89,726 samples. In these datasets, informed conditioning outperforms logistic regression for 115 of the 157 known associated variants investigated (P-value = 1 × 10(-9)). The improvement varied across diseases with a 16% median increase in χ(2) test statistics and a commensurate increase in power. This suggests that applying our method to existing and future association studies of these diseases may identify novel disease loci.

Common genetic variants explain the majority of the correlation between height and intelligence: the generation Scotland study.

Greater height and higher intelligence test scores are predictors of better health outcomes. Here, we used molecular (single-nucleotide polymorphism) data to estimate the genetic correlation between height and general intelligence (g) in 6,815 unrelated subjects (median age 57, IQR 49-63) from the Generation Scotland: Scottish Family Health Study cohort. The phenotypic correlation between height and g was 0.16 (SE 0.01). The genetic correlation between height and g was 0.28 (SE 0.09) with a bivariate heritability estimate of 0.71. Understanding the molecular basis of the correlation between height and intelligence may help explain any shared role in determining health outcomes. This study identified a modest genetic correlation between height and intelligence with the majority of the phenotypic correlation being explained by shared genetic influences.

Identification of BACH2 and RAD51B as rheumatoid arthritis susceptibility loci in a meta-analysis of genome-wide data.

OBJECTIVE: A recent high-density fine-mapping (ImmunoChip) study of genetic associations in rheumatoid arthritis (RA) identified 14 risk loci with validated genome-wide significance, as well as a number of loci showing associations suggestive of significance (P = 5 × 10(-5) < 5 × 10(-8)), but these have yet to be replicated. The aim of this study was to determine whether these potentially significant loci are involved in the pathogenesis of RA, and to explore whether any of the loci are associated with a specific RA serotype. METHODS: A total of 16 single-nucleotide polymorphisms (SNPs) were selected for genotyping and association analyses in 2 independent validation cohorts, comprising 6,106 RA cases and 4,290 controls. A meta-analysis of the data from the original ImmunoChip discovery cohort and from both validation cohorts was carried out, for a combined total of 17,581 RA cases and 20,160 controls. In addition, stratified analysis of patient subsets, defined according to their anti-cyclic citrullinated peptide (anti-CCP) antibody status, was performed. RESULTS: A significant association with RA risk (P < 0.05) was replicated for 6 of the SNPs assessed in the validation cohorts. All SNPs in the validation study had odds ratios (ORs) for RA susceptibility in the same direction as those in the ImmunoChip discovery study. One SNP, rs72928038, mapping to an intron of BACH2, achieved genome-wide significance in the meta-analysis (P = 1.2 × 10(-8), OR 1.12), and a second SNP, rs911263, mapping to an intron of RAD51B, was significantly associated in the anti-CCP-positive RA subgroup (P = 4 × 10(-8), OR 0.89), confirming that both are RA susceptibility loci. CONCLUSION: This study provides robust evidence for an association of RA susceptibility with genes involved in B cell differentiation (BACH2) and DNA repair (RAD51B). The finding that the RAD51B gene exhibited different associations based on serologic subtype adds to the expanding knowledge base in defining subgroups of RA.

Molecular genetic contributions to socioeconomic status and intelligence.

Education, socioeconomic status, and intelligence are commonly used as predictors of health outcomes, social environment, and mortality. Education and socioeconomic status are typically viewed as environmental variables although both correlate with intelligence, which has a substantial genetic basis. Using data from 6815 unrelated subjects from the Generation Scotland study, we examined the genetic contributions to these variables and their genetic correlations. Subjects underwent genome-wide testing for common single nucleotide polymorphisms (SNPs). DNA-derived heritability estimates and genetic correlations were calculated using the 'Genome-wide Complex Trait Analyses' (GCTA) procedures. 21% of the variation in education, 18% of the variation in socioeconomic status, and 29% of the variation in general cognitive ability was explained by variation in common SNPs (SEs ~ 5%). The SNP-based genetic correlations of education and socioeconomic status with general intelligence were 0.95 (SE 0.13) and 0.26 (0.16), respectively. There are genetic contributions to intelligence and education with near-complete overlap between common additive SNP effects on these traits (genetic correlation ~ 1). Genetic influences on socioeconomic status are also associated with the genetic foundations of intelligence. The results are also compatible with substantial environmental contributions to socioeconomic status.

Genome-wide association study meta-analysis of chronic widespread pain: evidence for involvement of the 5p15.2 region.

BACKGROUND AND OBJECTIVES: Chronic widespread pain (CWP) is a common disorder affecting ∼10% of the general population and has an estimated heritability of 48-52%. In the first large-scale genome-wide association study (GWAS) meta-analysis, we aimed to identify common genetic variants associated with CWP. METHODS: We conducted a GWAS meta-analysis in 1308 female CWP cases and 5791 controls of European descent, and replicated the effects of the genetic variants with suggestive evidence for association in 1480 CWP cases and 7989 controls. Subsequently, we studied gene expression levels of the nearest genes in two chronic inflammatory pain mouse models, and examined 92 genetic variants previously described associated with pain. RESULTS: The minor C-allele of rs13361160 on chromosome 5p15.2, located upstream of chaperonin-containing-TCP1-complex-5 gene (CCT5) and downstream of FAM173B, was found to be associated with a 30% higher risk of CWP (minor allele frequency=43%; OR=1.30, 95% CI 1.19 to 1.42, p=1.2×10(-8)). Combined with the replication, we observed a slightly attenuated OR of 1.17 (95% CI 1.10 to 1.24, p=4.7×10(-7)) with moderate heterogeneity (I2=28.4%). However, in a sensitivity analysis that only allowed studies with joint-specific pain, the combined association was genome-wide significant (OR=1.23, 95% CI 1.14 to 1.32, p=3.4×10(-8), I2=0%). Expression levels of Cct5 and Fam173b in mice with inflammatory pain were higher in the lumbar spinal cord, not in the lumbar dorsal root ganglions, compared to mice without pain. None of the 92 genetic variants previously described were significantly associated with pain (p>7.7×10(-4)). CONCLUSIONS: We identified a common genetic variant on chromosome 5p15.2 associated with joint-specific CWP in humans. This work suggests that CCT5 and FAM173B are promising targets in the regulation of pain.

Genome sequencing with gene panel-based analysis for rare inherited conditions in a publicly funded healthcare system: implications for future testing.

NHS genetics centres in Scotland sought to investigate the Genomics England 100,000 Genomes Project diagnostic utility to evaluate genome sequencing for in rare, inherited conditions. Four regional services recruited 999 individuals from 394 families in 200 rare phenotype categories, with negative historic genetic testing. Genome sequencing was performed at Edinburgh Genomics, and phenotype and sequence data were transferred to Genomics England for variant calling, gene-based filtering and variant prioritisation. NHS Scotland genetics laboratories performed interpretation, validation and reporting. New diagnoses were made in 23% cases - 19% in genes implicated in disease at the time of variant prioritisation, and 4% from later review of additional genes. Diagnostic yield varied considerably between phenotype categories and was minimal in cases with prior exome testing. Genome sequencing with gene panel filtering and reporting achieved improved diagnostic yield over previous historic testing but not over now routine trio-exome sequence tests. Re-interpretation of genomic data with updated gene panels modestly improved diagnostic yield at minimal cost. However, to justify the additional costs of genome vs exome sequencing, efficient methods for analysis of structural variation will be required and / or cost of genome analysis and storage will need to decrease.

Relationship between area-level socio-economic deprivation and autoantibody status in patients with rheumatoid arthritis: multicentre cross-sectional study.

OBJECTIVES: The aims of this study were to assess the association between area-level socio-economic deprivation and the phenotype of rheumatoid arthritis (RA), defined by rheumatoid factor (RF) and anticitrullinated peptide antibody (AC PA) status, and to determine whether any observed association can be explained by smoking. METHODS: The authors performed logistic regression analysis of 6298 patients with RA, defined by American College of Rheumatology classification criteria modified for genetic studies. Analysis was stratified by cohort/recruitment centre. Socio-economic deprivation was measured using the Townsend Index. RESULTS: Deprivation predicted RF but not ACPA positivity, independent of smoking. The ORs for trend across tertiles, adjusted for smoking, gender, period of birth and cohort/recruitment centre, were 1.14 (95% CI 1.01 to 1.29) for RF and 1.01 (95% CI 0.87 to 1.16) for ACPA. Even after adjusting for deprivation, smoking was strongly associated with ACPA positivity (OR 1.38, 95% CI 1.22 to 1.55). There was no evidence of any effect modification by the RA risk alleles (HLA-DRB1 shared epitope and PTPN22 rs2476601) that have previously been shown to modify the effect of smoking on ACPA and RF positivity. CONCLUSIONS: Among patients with RA, deprivation predicted RF positivity but not ACPA positivity. The effect of deprivation did not appear to be explained by smoking. Deprivation may be a marker for previously unrecognised, potentially modifiable environmental influences on the immunological phenotype of RA. Furthermore, given the known associations of RF positivity with prognosis and response to treatment in RA, these findings have potential implications for resource allocation and healthcare delivery.

Identification of AF4/FMR2 family, member 3 (AFF3) as a novel rheumatoid arthritis susceptibility locus and confirmation of two further pan-autoimmune susceptibility genes.

The concept of susceptibility genes common to different autoimmune diseases is now firmly established with previous studies demonstrating overlap of loci conferring susceptibility to type 1 diabetes (T1D) with both Coeliac disease and multiple sclerosis. Rheumatoid arthritis (RA) is an archetypal autoimmune disease and we, therefore, targeted putative T1D susceptibility loci for genotyping in UK RA cases and unrelated controls. A novel RA susceptibility locus at AFF3 was identified with convincing evidence for association in a combined sample cohort of 6819 RA cases and 12 650 controls [OR 1.12 95% confidence intervals (CI) 1.07-1.17, P = 2.8 x 10(-7)]. Association of two previously described loci (CTLA-4 and 4q27) with RA was also replicated (OR 0.87, 95% CI 0.82-0.94, P = 1.1 x 10(-4) and OR 0.86, 95% CI 0.79-0.94, P = 5.4 x 10(-4), respectively). These findings take the number of established RA susceptibility loci to 13, only one of which has not been associated with other autoimmune diseases.

Combined effects of three independent SNPs greatly increase the risk estimate for RA at 6q23.

The most consistent finding derived from the WTCCC GWAS for rheumatoid arthritis (RA) was association to a SNP at 6q23. We performed a fine-mapping of the region in order to search the 6q23 region for additional disease variants. 3962 RA patients and 3531 healthy controls were included in the study. We found 18 SNPs associated with RA. The SNP showing the strongest association was rs6920220 [P = 2.6 x 10(-6), OR (95% CI) 1.22 (1.13-1.33)]. The next most strongly associated SNP was rs13207033 [P = 0.0001, OR (95% CI) 0.86 (0.8-0.93)] which was perfectly correlated with rs10499194, a SNP previously associated with RA in a US/European series. Additionally, we found a number of new potential RA markers, including rs5029937, located in the intron 2 of TNFAIP3. Of the 18 associated SNPs, three polymorphisms, rs6920220, rs13207033 and rs5029937, remained significant after conditional logistic regression analysis. The combination of the carriage of both risk alleles of rs6920220 and rs5029937 together with the absence of the protective allele of rs13207033 was strongly associated with RA when compared with carriage of none [OR of 1.86 (95% CI) (1.51-2.29)]. This equates to an effect size of 1.50 (95% CI 1.21-1.85) compared with controls and is higher than that obtained for any SNP individually. This is the first study to show that the confirmed loci from the GWA studies, that confer only a modest effect size, could harbour a significantly greater effect once the effect of additional risk variants are accounted for.

Common genetic determinants of vitamin D insufficiency: a genome-wide association study.

BACKGROUND: Vitamin D is crucial for maintenance of musculoskeletal health, and might also have a role in extraskeletal tissues. Determinants of circulating 25-hydroxyvitamin D concentrations include sun exposure and diet, but high heritability suggests that genetic factors could also play a part. We aimed to identify common genetic variants affecting vitamin D concentrations and risk of insufficiency. METHODS: We undertook a genome-wide association study of 25-hydroxyvitamin D concentrations in 33 996 individuals of European descent from 15 cohorts. Five epidemiological cohorts were designated as discovery cohorts (n=16 125), five as in-silico replication cohorts (n=9367), and five as de-novo replication cohorts (n=8504). 25-hydroxyvitamin D concentrations were measured by radioimmunoassay, chemiluminescent assay, ELISA, or mass spectrometry. Vitamin D insufficiency was defined as concentrations lower than 75 nmol/L or 50 nmol/L. We combined results of genome-wide analyses across cohorts using Z-score-weighted meta-analysis. Genotype scores were constructed for confirmed variants. FINDINGS: Variants at three loci reached genome-wide significance in discovery cohorts for association with 25-hydroxyvitamin D concentrations, and were confirmed in replication cohorts: 4p12 (overall p=1.9x10(-109) for rs2282679, in GC); 11q12 (p=2.1x10(-27) for rs12785878, near DHCR7); and 11p15 (p=3.3x10(-20) for rs10741657, near CYP2R1). Variants at an additional locus (20q13, CYP24A1) were genome-wide significant in the pooled sample (p=6.0x10(-10) for rs6013897). Participants with a genotype score (combining the three confirmed variants) in the highest quartile were at increased risk of having 25-hydroxyvitamin D concentrations lower than 75 nmol/L (OR 2.47, 95% CI 2.20-2.78, p=2.3x10(-48)) or lower than 50 nmol/L (1.92, 1.70-2.16, p=1.0x10(-26)) compared with those in the lowest quartile. INTERPRETATION: Variants near genes involved in cholesterol synthesis, hydroxylation, and vitamin D transport affect vitamin D status. Genetic variation at these loci identifies individuals who have substantially raised risk of vitamin D insufficiency. FUNDING: Full funding sources listed at end of paper (see Acknowledgments).