Estimating overdiagnosis in giant cell arteritis diagnostic pathways using genetic data: genetic association study.

OBJECTIVES: Giant cell arteritis (GCA) may be confirmed by temporal artery biopsy (TAB) but false negatives can occur. GCA may be overdiagnosed in TAB-negative cases, or if neither TAB nor imaging is done. We used Human Leucocyte Antigen (HLA) genetic association of TAB-positive GCA as an "unbiased umpire" test to estimate historic overdiagnosis of GCA. METHODS: Patients diagnosed with GCA between 1990-2014 were genotyped. During this era, vascular imaging alone was rarely used to diagnose GCA. HLA region variants were jointly imputed from genome-wide genotypic data of cases and controls. Per-allele frequencies across all HLA variants with p< 1.0x1 0 -5 were compared with population control data to estimate overdiagnosis rates in cases without a positive TAB. RESULTS: Genetic data from 663 GCA patients were compared with data from 2619 population controls. TAB-negative GCA (n = 147) and GCA without TAB result (n = 160) had variant frequencies intermediate between TAB-positive GCA (n = 356) and population controls. For example, the allele frequency of HLA-DRB1*04 was 32% for TAB-positive GCA, 29% for GCA without TAB result, 27% for TAB-negative GCA and 20% in population controls. Making several strong assumptions, we estimated that around two-thirds of TAB-negative cases and one-third of cases without TAB result may have been overdiagnosed. From these data, TAB sensitivity is estimated as 88%. CONCLUSIONS: Conservatively assuming 95% specificity, TAB has a negative likelihood ratio of around 0.12. Our method for utilising standard genotyping data as an "unbiased umpire" might be used as a way of comparing the accuracy of different diagnostic pathways.

Common co-morbidities in polymyalgia rheumatica and giant cell arteritis: cross-sectional study in UK Biobank.

Objective: The aim was to determine prevalent co-morbidities in cases with PMR or GCA compared with matched controls. Methods: This was a nested, cross-sectional case-control study within the UK Biobank, which recruited participants aged 40-69 years. Case status was defined as self-reported prior diagnosis of PMR or GCA. Ten controls per case were matched for age, sex, ethnicity and assessment centre. Associations with selected self-reported co-morbidities were studied using conditional logistic regression. Results: Of PMR (n = 1036) or GCA (n = 102) cases, 72% were female, 98% White, and 58% reported current use of glucocorticoids. Mean age was 63 years. At the time of the assessment visit, compared with controls, PMR/GCA cases were more likely to report poor general health and at least several days of low mood in the past 2 weeks. PMR was associated with hypothyroidism [odds ratio (OR) = 1.34; 95% CI = 1.07, 1.67] and ever-use of HRT (OR = 1.26; 95% CI = 1.07, 1.47). Regarding common co-morbidities, PMR and GCA were both associated with hypertension (PMR: OR = 1.21; 95% CI = 1.06, 1.39; GCA: OR = 1.86; 95% CI = 1.23, 2.81) and cataract (PMR: OR = 1.51; 95% CI = 1.19, 1.93; GCA: OR = 3.84; 95% CI = 2.23, 6.60). Additionally, GCA was associated with depression (OR = 3.05; 95% CI = 1.59, 5.85). Neither condition was associated with diabetes. Conclusion: Participants with a history of PMR/GCA, including those not currently taking glucocorticoids, rated their health as poorer than matched controls. Some previously described disease associations (hypothyroidism and early menopause) were replicated. Hypertension and cataract, both of which can be exacerbated by long-term glucocorticoid therapy, were over-represented in both diseases, particularly GCA.

Comprehensive genetic and functional analyses of Fc gamma receptors influence on response to rituximab therapy for autoimmunity.

BACKGROUND: Rituximab is widely used to treat autoimmunity but clinical response varies. Efficacy is determined by the efficiency of B-cell depletion, which may depend on various Fc gamma receptor (FcγR)-dependent mechanisms. Study of FcγR is challenging due to the complexity of the FCGR genetic locus. We sought to assess the effect of FCGR variants on clinical response, B-cell depletion and NK-cell-mediated killing in rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE). METHODS: A longitudinal cohort study was conducted in 835 patients [RA = 573; SLE = 262]. Clinical outcome measures were two-component disease activity score in 28-joints (2C-DAS28CRP) for RA and British Isles Lupus Assessment Group (BILAG)-2004 major clinical response (MCR) for SLE at 6 months. B-cells were evaluated by highly-sensitive flow cytometry. Single nucleotide polymorphism and copy number variation for genes encoding five FcγRs were measured using multiplex ligation-dependent probe amplification. Ex vivo studies assessed NK-cell antibody-dependent cellular cytotoxicity (ADCC) and FcγR expression. FINDINGS: In RA, carriage of FCGR3A-158V and increased FCGR3A-158V copies were associated with greater 2C-DAS28CRP response (adjusted for baseline 2C-DAS28CRP). In SLE, MCR was associated with increased FCGR3A-158V, OR 1.64 (95% CI 1.12-2.41) and FCGR2C-ORF OR 1.93 (95% CI 1.09-3.40) copies. 236/413 (57%) patients with B-cell data achieved complete depletion. Homozygosity for FCGR3A-158V and increased FCGR3A-158V copies were associated with complete depletion in combined analyses. FCGR3A genotype was associated with rituximab-induced ADCC, and increased NK-cell FcγRIIIa expression was associated with improved clinical response and depletion in vivo. Furthermore, disease status and concomitant therapies impacted both NK-cell FcγRIIIa expression and ADCC. INTERPRETATION: FcγRIIIa is the major low affinity FcγR associated with rituximab response. Increased copies of the FCGR3A-158V allele (higher affinity for IgG1), influences clinical and biological responses to rituximab in autoimmunity. Enhancing FcγR-effector functions could improve the next generation of CD20-depleting therapies and genotyping may stratify patients for optimal treatment protocols. FUNDING: Medical Research Council, National Institute for Health and Care Research, Versus Arthritis.

Artificial Intelligence-Assisted Amphiregulin and Epiregulin IHC Predicts Panitumumab Benefit in RAS Wild-Type Metastatic Colorectal Cancer.

PURPOSE: High tumor mRNA levels of the EGFR ligands amphiregulin (AREG) and epiregulin (EREG) are associated with anti-EGFR agent response in metastatic colorectal cancer (mCRC). However, ligand RNA assays have not been adopted into routine practice due to issues with analytic precision and practicality. We investigated whether AREG/EREG IHC could predict benefit from the anti-EGFR agent panitumumab. EXPERIMENTAL DESIGN: Artificial intelligence algorithms were developed to assess AREG/EREG IHC in 274 patients from the PICCOLO trial of irinotecan with or without panitumumab (Ir vs. IrPan) in RAS wild-type mCRC. The primary endpoint was progression-free survival (PFS). Secondary endpoints were RECIST response rate (RR) and overall survival (OS). Models were repeated adjusting separately for BRAF mutation status and primary tumor location (PTL). RESULTS: High ligand expression was associated with significant PFS benefit from IrPan compared with Ir [8.0 vs. 3.2 months; HR, 0.54; 95% confidence interval (CI), 0.37-0.79; P = 0.001]; whereas low ligand expression was not (3.4 vs. 4.4 months; HR, 1.05; 95% CI, 0.74-1.49; P = 0.78). The ligand-treatment interaction was significant (P interaction = 0.02) and remained significant after adjustment for BRAF-mutation status and PTL. Likewise, RECIST RR was significantly improved in patients with high ligand expression (IrPan vs. Ir: 48% vs. 6%; P < 0.0001) but not those with low ligand expression (25% vs. 14%; P = 0.10; P interaction = 0.01). The effect on OS was similar but not statistically significant. CONCLUSIONS: AREG/EREG IHC identified patients who benefitted from the addition of panitumumab to irinotecan chemotherapy. IHC is a practicable assay that may be of use in routine practice.

Microbiome Analysis of More Than 2,000 NHS Bowel Cancer Screening Programme Samples Shows the Potential to Improve Screening Accuracy.

PURPOSE: There is potential for fecal microbiome profiling to improve colorectal cancer screening. This has been demonstrated by research studies, but it has not been quantified at scale using samples collected and processed routinely by a national screening program. EXPERIMENTAL DESIGN: Between 2016 and 2019, the largest of the NHS Bowel Cancer Screening Programme hubs prospectively collected processed guaiac fecal occult blood test (gFOBT) samples with subsequent colonoscopy outcomes: blood-negative [n = 491 (22%)]; colorectal cancer [n = 430 (19%)]; adenoma [n = 665 (30%)]; colonoscopy-normal [n = 300 (13%)]; nonneoplastic [n = 366 (16%)]. Samples were transported and stored at room temperature. DNA underwent 16S rRNA gene V4 amplicon sequencing. Taxonomic profiling was performed to provide features for classification via random forests (RF). RESULTS: Samples provided 16S amplicon-based microbial profiles, which confirmed previously described colorectal cancer-microbiome associations. Microbiome-based RF models showed potential as a first-tier screen, distinguishing colorectal cancer or neoplasm (colorectal cancer or adenoma) from blood-negative with AUC 0.86 (0.82-0.89) and AUC 0.78 (0.74-0.82), respectively. Microbiome-based models also showed potential as a second-tier screen, distinguishing from among gFOBT blood-positive samples, colorectal cancer or neoplasm from colonoscopy-normal with AUC 0.79 (0.74-0.83) and AUC 0.73 (0.68-0.77), respectively. Models remained robust when restricted to 15 taxa, and performed similarly during external validation with metagenomic datasets. CONCLUSIONS: Microbiome features can be assessed using gFOBT samples collected and processed routinely by a national colorectal cancer screening program to improve accuracy as a first- or second-tier screen. The models required as few as 15 taxa, raising the potential of an inexpensive qPCR test. This could reduce the number of colonoscopies in countries that use fecal occult blood test screening.

Germline variants are associated with increased primary melanoma tumor thickness at diagnosis.

Germline genetic variants have been identified, which predispose individuals and families to develop melanoma. Tumor thickness is the strongest predictor of outcome for clinically localized primary melanoma patients. We sought to determine whether there is a heritable genetic contribution to variation in tumor thickness. If confirmed, this will justify the search for specific genetic variants influencing tumor thickness. To address this, we estimated the proportion of variation in tumor thickness attributable to genome-wide genetic variation (variant-based heritability) using unrelated patients with measured primary cutaneous melanoma thickness. As a secondary analysis, we conducted a genome-wide association study (GWAS) of tumor thickness. The analyses utilized 10 604 individuals with primary cutaneous melanoma drawn from nine GWAS datasets from eight cohorts recruited from the general population, primary care and melanoma treatment centers. Following quality control and filtering to unrelated individuals with study phenotypes, 8125 patients were used in the primary analysis to test whether tumor thickness is heritable. An expanded set of 8505 individuals (47.6% female) were analyzed for the secondary GWAS meta-analysis. Analyses were adjusted for participant age, sex, cohort and ancestry. We found that 26.6% (SE 11.9%, P = 0.0128) of variation in tumor thickness is attributable to genome-wide genetic variation. While requiring replication, a chromosome 11 locus was associated (P < 5 × 10-8) with tumor thickness. Our work indicates that sufficiently large datasets will enable the discovery of genetic variants associated with greater tumor thickness, and this will lead to the identification of host biological processes influencing melanoma growth and invasion.

Association between a 46-SNP Polygenic Risk Score and melanoma risk in Dutch patients with familial melanoma.

BACKGROUND: Familial clustering of melanoma suggests a shared genetic predisposition among family members, but only 10%-40% of familial cases carry a pathogenic variant in a known high-risk melanoma susceptibility gene. We investigated whether a melanoma-specific Polygenic Risk Score (PRS) is associated with melanoma risk in patients with genetically unexplained familial melanoma. METHODS: Dutch familial melanoma cases (n=418) were genotyped for 46 SNPs previously identified as independently associated with melanoma risk. The 46-SNP PRS was calculated and standardised to 3423 healthy controls (sPRS) and the association between PRS and melanoma risk was modelled using logistic regression. Within the case series, possible differences were further explored by investigating the PRS in relation to (1) the number of primary melanomas in a patient and (2) the extent of familial clustering of melanoma. RESULTS: The PRS was significantly associated with melanoma risk, with a per-SD OR of 2.12 (95% CI 1.90 to 2.35, p<0.001), corresponding to a 5.70-fold increased risk (95% CI 3.93 to 8.28) when comparing the top 90th to the middle 40-60th PRS percentiles. The mean PRS was significantly higher in cases with multiple primary melanomas than in cases with a single melanoma (sPRS 1.17 vs 0.71, p=0.001). Conversely, cases from high-density melanoma families had a lower (but non-significant) mean PRS than cases from low-density families (sPRS 0.60 vs 0.94, p=0.204). CONCLUSION: Our work underlines the significance of a PRS in determining melanoma susceptibility and encourages further exploration of the diagnostic value of a PRS in genetically unexplained melanoma families.

Investigation of genetically regulated gene expression and response to treatment in rheumatoid arthritis highlights an association between IL18RAP expression and treatment response.

OBJECTIVES: In this study, we sought to investigate whether there was any association between genetically regulated gene expression (as predicted using various reference panels) and anti-tumour necrosis factor (anti-TNF) treatment response (change in erythrocyte sedimentation rate (ESR)) using 3158 European ancestry patients with rheumatoid arthritis. METHODS: The genetically regulated portion of gene expression was estimated in the full cohort of 3158 subjects (as well as within a subcohort consisting of 1575 UK patients) using the PrediXcan software package with three different reference panels. Estimated expression was tested for association with anti-TNF treatment response. As a replication/validation experiment, we also investigated the correlation between change in ESR with measured gene expression at the Interleukin 18 Receptor Accessory Protein (IL18RAP) gene in whole blood and synovial tissue, using an independent replication data set of patients receiving conventional synthetic disease modifying anti-rheumatic drugs, with directly measured (via RNA sequencing) gene expression. RESULTS: We found that predicted expression of IL18RAP showed a consistent signal of association with treatment response across the reference panels. In our independent replication data set, IL18RAP expression in whole blood showed correlation with the change in ESR between baseline and follow-up (r=-0.35, p=0.0091). Change in ESR was also correlated with the expression of IL18RAP in synovial tissue (r=-0.28, p=0.02). CONCLUSION: Our results suggest that IL18RAP expression is worthy of further investigation as a potential predictor of treatment response in rheumatoid arthritis that is not specific to a particular drug type.

Massively parallel reporter assays of melanoma risk variants identify MX2 as a gene promoting melanoma.

Genome-wide association studies (GWAS) have identified ~20 melanoma susceptibility loci, most of which are not functionally characterized. Here we report an approach integrating massively-parallel reporter assays (MPRA) with cell-type-specific epigenome and expression quantitative trait loci (eQTL) to identify susceptibility genes/variants from multiple GWAS loci. From 832 high-LD variants, we identify 39 candidate functional variants from 14 loci displaying allelic transcriptional activity, a subset of which corroborates four colocalizing melanocyte cis-eQTL genes. Among these, we further characterize the locus encompassing the HIV-1 restriction gene, MX2 (Chr21q22.3), and validate a functional intronic variant, rs398206. rs398206 mediates the binding of the transcription factor, YY1, to increase MX2 levels, consistent with the cis-eQTL of MX2 in primary human melanocytes. Melanocyte-specific expression of human MX2 in a zebrafish model demonstrates accelerated melanoma formation in a BRAFV600E background. Our integrative approach streamlines GWAS follow-up studies and highlights a pleiotropic function of MX2 in melanoma susceptibility.

Genetic associations with radiological damage in rheumatoid arthritis: Meta-analysis of seven genome-wide association studies of 2,775 cases.

BACKGROUND: Previous studies of radiological damage in rheumatoid arthritis (RA) have used candidate-gene approaches, or evaluated single genome-wide association studies (GWAS). We undertook the first meta-analysis of GWAS of RA radiological damage to: (1) identify novel genetic loci for this trait; and (2) test previously validated variants. METHODS: Seven GWAS (2,775 RA cases, of a range of ancestries) were combined in a meta-analysis. Radiological damage was assessed using modified Larsen scores, Sharp van Der Heijde scores, and erosive status. Single nucleotide polymophsim (SNP) associations with radiological damage were tested at a single time-point using regression models. Primary analyses included age and disease duration as covariates. Secondary analyses also included rheumatoid factor (RF). Meta-analyses were undertaken in trans-ethnic and European-only cases. RESULTS: In the trans-ethnic primary meta-analysis, one SNP (rs112112734) in close proximity to HLA-DRB1, and strong linkage disequilibrium with the shared-epitope, attained genome-wide significance (P = 4.2x10-8). In the secondary analysis (adjusting for RF) the association was less significant (P = 1.7x10-6). In both trans-ethnic primary and secondary meta-analyses 14 regions contained SNPs with associations reaching P<5x10-6; in the European primary and secondary analyses 13 and 10 regions contained SNPs reaching P<5x10-6, respectively. Of the previously validated SNPs for radiological progression, only rs660895 (tagging HLA-DRB1*04:01) attained significance (P = 1.6x10-5) and had a consistent direction of effect across GWAS. CONCLUSIONS: Our meta-analysis confirms the known association between the HLA-DRB1 shared epitope and RA radiological damage. The lack of replication of previously validated non-HLA markers highlights a requirement for further research to deliver clinically-useful prognostic genetic markers.

Transcriptomic Analysis Reveals Prognostic Molecular Signatures of Stage I Melanoma.

PURPOSE: Previously identified transcriptomic signatures have been based on primary and metastatic melanomas with relatively few American Joint Committee on Cancer (AJCC) stage I tumors, given difficulties in sampling small tumors. The advent of adjuvant therapies has highlighted the need for better prognostic and predictive biomarkers, especially for AJCC stage I and stage II disease. EXPERIMENTAL DESIGN: A total of 687 primary melanoma transcriptomes were generated from the Leeds Melanoma Cohort (LMC). The prognostic value of existing signatures across all the AJCC stages was tested. Unsupervised clustering was performed, and the prognostic value of the resultant signature was compared with that of sentinel node biopsy (SNB) and tested as a biomarker in three published immunotherapy datasets. RESULTS: Previous Lund and The Cancer Genome Atlas signatures predicted outcome in the LMC dataset (P = 10-8 to 10-4) but showed a significant interaction with AJCC stage (P = 0.04) and did not predict outcome in stage I tumors (P = 0.3-0.7). Consensus-based classification of the LMC dataset identified six classes that predicted outcome, notably in stage I disease. LMC class was a similar indicator of prognosis when compared with SNB, and it added prognostic value to the genes reported by Gerami and colleagues. One particular LMC class consistently predicted poor outcome in patients receiving immunotherapy in two of three tested datasets. Biological characterization of this class revealed high JUN and AXL expression and evidence of epithelial-to-mesenchymal transition. CONCLUSIONS: A transcriptomic signature of primary melanoma was identified with prognostic value, including in stage I melanoma and in patients undergoing immunotherapy.

Validity of a two-component imaging-derived disease activity score for improved assessment of synovitis in early rheumatoid arthritis.

OBJECTIVES: Imaging of joint inflammation provides a standard against which to derive an updated DAS for RA. Our objectives were to develop and validate a DAS based on reweighting the DAS28 components to maximize association with US-assessed synovitis. METHODS: Early RA patients from two observational cohorts (n = 434 and n = 117) and a clinical trial (n = 59) were assessed at intervals up to 104 weeks from baseline; all US scans were within 1 week of clinical exam. There were 899, 163 and 183 visits in each cohort. Associations of combined US grey scale and power Doppler scores (GSPD) with 28 tender joint count and 28 swollen joint count (SJC28), CRP, ESR and general health visual analogue scale were examined in linear mixed model regressions. Cross-validation evaluated model predictive ability. Coefficients learned from training data defined a re-weighted DAS28 that was validated against radiographic progression in independent data (3037 observations; 717 patients). RESULTS: Of the conventional DAS28 components only SJC28 and CRP were associated with GSPD in all three development cohorts. A two-component model including SJC28 and CRP outperformed a four-component model (R2 = 0.235, 0.392, 0.380 vs 0.232, 0.380, 0.375, respectively). The re-weighted two-component DAS28CRP outperformed conventional DAS28 definitions in predicting GSPD (Δtest log-likelihood <-2.6, P < 0.01), Larsen score and presence of erosions. CONCLUSION: A score based on SJC28 and CRP alone demonstrated stronger associations with synovitis and radiographic progression than the original DAS28 and should be considered in research on pathophysiological manifestations of early RA. Implications for clinical management of RA remain to be established.

Publisher Correction: Novel pleiotropic risk loci for melanoma and nevus density implicate multiple biological pathways.

The original version of this Article contained errors in the spelling of the authors Fan Liu and M. Arfan Ikram, which were incorrectly given as Fan Lui and Arfan M. Ikram. In addition, the original version of this Article also contained errors in the author affiliations which are detailed in the associated Publisher Correction.

Novel pleiotropic risk loci for melanoma and nevus density implicate multiple biological pathways.

The total number of acquired melanocytic nevi on the skin is strongly correlated with melanoma risk. Here we report a meta-analysis of 11 nevus GWAS from Australia, Netherlands, UK, and USA comprising 52,506 individuals. We confirm known loci including MTAP, PLA2G6, and IRF4, and detect novel SNPs in KITLG and a region of 9q32. In a bivariate analysis combining the nevus results with a recent melanoma GWAS meta-analysis (12,874 cases, 23,203 controls), SNPs near GPRC5A, CYP1B1, PPARGC1B, HDAC4, FAM208B, DOCK8, and SYNE2 reached global significance, and other loci, including MIR146A and OBFC1, reached a suggestive level. Overall, we conclude that most nevus genes affect melanoma risk (KITLG an exception), while many melanoma risk loci do not alter nevus count. For example, variants in TERC and OBFC1 affect both traits, but other telomere length maintenance genes seem to affect melanoma risk only. Our findings implicate multiple pathways in nevogenesis.

Prediction of treatment response in rheumatoid arthritis patients using genome-wide SNP data.

Although a number of treatments are available for rheumatoid arthritis (RA), each of them shows a significant nonresponse rate in patients. Therefore, predicting a priori the likelihood of treatment response would be of great patient benefit. Here, we conducted a comparison of a variety of statistical methods for predicting three measures of treatment response, between baseline and 3 or 6 months, using genome-wide SNP data from RA patients available from the MAximising Therapeutic Utility in Rheumatoid Arthritis (MATURA) consortium. Two different treatments and 11 different statistical methods were evaluated. We used 10-fold cross validation to assess predictive performance, with nested 10-fold cross validation used to tune the model hyperparameters when required. Overall, we found that SNPs added very little prediction information to that obtained using clinical characteristics only, such as baseline trait value. This observation can be explained by the lack of strong genetic effects and the relatively small sample sizes available; in analysis of simulated and real data, with larger effects and/or larger sample sizes, prediction performance was much improved. Overall, methods that were consistent with the genetic architecture of the trait were able to achieve better predictive ability than methods that were not. For treatment response in RA, methods that assumed a complex underlying genetic architecture achieved slightly better prediction performance than methods that assumed a simplified genetic architecture.

Assessing the Incremental Contribution of Common Genomic Variants to Melanoma Risk Prediction in Two Population-Based Studies.

It is unclear to what degree genomic and traditional (phenotypic and environmental) risk factors overlap in their prediction of melanoma risk. We evaluated the incremental contribution of common genomic variants (in pigmentation, nevus, and other pathways) and their overlap with traditional risk factors, using data from two population-based case-control studies from Australia (n = 1,035) and the United Kingdom (n = 1,460) that used the same questionnaires. Polygenic risk scores were derived from 21 gene regions associated with melanoma and odds ratios from published meta-analyses. Logistic regression models were adjusted for age, sex, center, and ancestry. Adding the polygenic risk score to a model with traditional risk factors increased the area under the receiver operating characteristic curve (AUC) by 2.3% (P = 0.003) for Australia and by 2.8% (P = 0.002) for Leeds. Gene variants in the pigmentation pathway, particularly MC1R, were responsible for most of the incremental improvement. In a cross-tabulation of polygenic by traditional tertile risk scores, 59% (Australia) and 49% (Leeds) of participants were categorized in the same (concordant) tertile. Of participants with low traditional risk, 9% (Australia) and 21% (Leeds) had high polygenic risk. Testing of genomic variants can identify people who are susceptible to melanoma despite not having a traditional phenotypic risk profile.

Genome-wide association study of response to methotrexate in early rheumatoid arthritis patients.

Methotrexate (MTX) monotherapy is a common first treatment for rheumatoid arthritis (RA), but many patients do not respond adequately. In order to identify genetic predictors of response, we have combined data from two consortia to carry out a genome-wide study of response to MTX in 1424 early RA patients of European ancestry. Clinical endpoints were change from baseline to 6 months after starting treatment in swollen 28-joint count, tender 28-joint count, C-reactive protein and the overall 3-component disease activity score (DAS28). No single nucleotide polymorphism (SNP) reached genome-wide statistical significance for any outcome measure. The strongest evidence for association was with rs168201 in NRG3 (p = 10-7 for change in DAS28). Some support was also seen for association with ZMIZ1, previously highlighted in a study of response to MTX in juvenile idiopathic arthritis. Follow-up in two smaller cohorts of 429 and 177 RA patients did not support these findings, although these cohorts were more heterogeneous.

Publisher correction: Functional characterization of a multi-cancer risk locus on chr5p15.33 reveals regulation of TERT by ZNF148.

This corrects the article DOI: 10.1038/ncomms15034.

Interaction between polymorphisms in aspirin metabolic pathways, regular aspirin use and colorectal cancer risk: A case-control study in unselected white European populations.

Regular aspirin use is associated with reduced risk of colorectal cancer (CRC). Variation in aspirin's chemoprevention efficacy has been attributed to the presence of single nucleotide polymorphisms (SNPs). We conducted a meta-analysis using two large population-based case-control datasets, the UK-Leeds Colorectal Cancer Study Group and the NIH-Colon Cancer Family Registry, having a combined total of 3325 cases and 2262 controls. The aim was to assess 42 candidate SNPs in 15 genes whose association with colorectal cancer risk was putatively modified by aspirin use, in the literature. Log odds ratios (ORs) and standard errors were estimated for each dataset separately using logistic regression adjusting for age, sex and study site, and dataset-specific results were combined using random effects meta-analysis. Meta-analysis showed association between SNPs rs6983267, rs11694911 and rs2302615 with CRC risk reduction (All P<0.05). Association for SNP rs6983267 in the CCAT2 gene only was noteworthy after multiple test correction (P = 0.001). Site-specific analysis showed association between SNPs rs1799853 and rs2302615 with reduced colon cancer risk only (P = 0.01 and P = 0.004, respectively), however neither reached significance threshold following multiple test correction. Meta-analysis of SNPs rs2070959 and rs1105879 in UGT1A6 gene showed interaction between aspirin use and CRC risk (Pinteraction = 0.01 and 0.02, respectively); stratification by aspirin use showed an association for decreased CRC risk for aspirin users having a wild-type genotype (rs2070959 OR = 0.77, 95% CI = 0.68-0.86; rs1105879 OR = 0.77 95% CI = 0.69-0.86) compared to variant allele cariers. The direction of the interaction however is in contrast to that published in studies on colorectal adenomas. Both SNPs showed potential site-specific interaction with aspirin use and colon cancer risk only (Pinteraction = 0.006 and 0.008, respectively), with the direction of association similar to that observed for CRC. Additionally, they showed interaction between any non-steroidal anti-inflammatory drugs (including aspirin) use and CRC risk (Pinteraction = 0.01 for both). All gene x environment (GxE) interactions however were not significant after multiple test correction. Candidate gene investigation indicated no evidence of GxE interaction between genetic variants in genes involved in aspirin pathways, regular aspirin use and colorectal cancer risk.