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.

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.

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.

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

Functional interaction between sequestosome-1/p62 and autophagy-linked FYVE-containing protein WDFY3 in human osteoclasts.

Paget's disease of bone (PDB) is a late-onset disorder characterised by focal areas of increased bone resorption, with osteoclasts that are increased in size, multinuclearity, number and activity. PDB-causing missense and nonsense variants in the gene encoding Sequestosome-1/p62 (SQSTM1) have been identified, all of which cluster in and around the ubiquitin-associated (UBA) domain of the protein. SQSTM1 is ubiquitously expressed and there is, as yet, no clear reason why these mutations only appear to cause an osteoclast-related phenotype. Using co-immunoprecipitation and tandem mass spectrometry, we identified a novel interaction in human osteoclast-like cells between SQSTM1 and Autophagy-Linked FYVE domain-containing protein (ALFY/WDFY3). Endogenous ALFY and SQSTM1 both localised within the nuclei of osteoclasts and their mononuclear precursors. When osteoclasts were starved to induce autophagy, SQSTM1 and ALFY relocated to the cytoplasm where they formed large aggregates, with cytoplasmic relocalisation appearing more rapid in mature osteoclasts than in precursors in the same culture. Overexpression of wild-type SQSTM1 in HEK293 cells also resulted in the formation of cytoplasmic aggregates containing SQSTM1 and endogenous ALFY, as did overexpression of a PDB-causing missense mutant form of SQSTM1, indicating that this mutation does not impair the formation of SQSTM1- and ALFY-containing aggregates. Expression of ALFY in bone cells has not previously been reported, and the process of autophagy has not been studied with respect to osteoclast activity. We have identified a functional interaction between SQSTM1 and ALFY in osteoclasts under conditions of cell stress. The difference in response to starvation between mature osteoclasts and their precursors may begin to explain the cell-specific functional effects of SQSTM1 mutations in PDB.

Investigating the viability of genetic screening/testing for RA susceptibility using combinations of five confirmed risk loci.

OBJECTIVE: Five loci-the shared epitope (SE) of HLA--DRB1, the PTPN22 gene, a locus on 6q23, the STAT4 gene and a locus mapping to the TRAF1/C5 genetic region--have now been unequivocally confirmed as conferring susceptibility to RA. The largest single effect is conferred by SE. We hypothesized that combinations of susceptibility alleles may increase risk over and above that of any individual locus alone. METHODS: We analysed data from 4238 RA cases and 1811 controls, for which genotypes were available at all five loci. RESULTS: Statistical analysis identified eight high-risk combinations conferring an odds ratio >6 compared with carriage of no susceptibility variants and, interestingly, 10% population controls carried a combination conferring high risk. All high-risk combinations included SE, and all but one contained PTPN22. Statistical modelling showed that a model containing only these two loci could achieve comparable sensitivity and specificity to a model including all five. Furthermore, replacing SE (which requires full subtyping at the HLA-DRB1 gene) with DRB1*1/4/10 carriage resulted in little further loss of information (correlation coefficient between models = 0.93). CONCLUSIONS: This represents the first exploration of the viability of population screening for RA and identifies several high-risk genetic combinations. However, given the population incidence of RA, genetic screening based on these loci alone is neither sufficiently sensitive nor specific at the current time.

Identification of a major locus for Paget's disease on chromosome 10p13 in families of British descent.

UNLABELLED: Mutations of SQSTM1 are an important cause of PDB, but other genes remain to be discovered. A major susceptibility locus for PDB was identified on chromosome 10p13 by a genome-wide linkage scan in families of British descent, which accounted for the vast majority of cases not caused by SQSTM1 mutations. INTRODUCTION: Paget's disease of bone (PDB) has a strong genetic component, and several susceptibility loci have been identified by genome-wide linkage scans. We previously identified three susceptibility loci for PDB using this approach on chromosomes 5q35, 2q36, and 10p13 in 62 families of mainly British descent, but subsequently, mutations in the SQSTM1 gene were found to be the cause of PDB in 23 families from this cohort. Here we reanalyzed the results of our genome-wide search in families from this cohort who did not have SQSTM1 mutations. MATERIALS AND METHODS: The study population consisted of 210 individuals from 39 families of predominantly British descent with autosomal dominant inheritance of PDB in whom SQSTM1 mutations had been excluded by mutation screening. The average family size was 5.44 +/- 3.98 (SD) individuals (range, 2-24 individuals). Genotyping was performed using standard techniques with 382 microsatellite markers spaced at an average distance of 9.06 cM throughout the autosomes. Multipoint linkage analysis was performed using the GENEHUNTER program under models of homogeneity and heterogeneity. RESULTS: Multipoint parametric linkage analysis under a model of homogeneity and nonparametric linkage analysis under a model of heterogeneity both showed strong evidence of linkage to a single locus on chromosome 10p13 (LOD score, +4.08) close to the marker D10S1653 at 41.43cM. No evidence of linkage was detected at the chromosome 2q36 locus previously identified in this population, and linkage to other candidate loci previously implicated in the pathogenesis of PDB was excluded. CONCLUSIONS: We conclude that there is an important susceptibility gene for PDB on chromosome 10p13 in families of British descent and find no evidence to support the existence of a susceptibility locus on chromosome 2q36 or other previously identified candidate loci for PDB in this population. The gene that lies within the 10p13 locus seems to account for the development of PDB in the vast majority of families of British descent who do not carry SQSTM1 mutations.