Geographic variation of mutagenic exposures in kidney cancer genomes.

International differences in the incidence of many cancer types indicate the existence of carcinogen exposures that have not yet been identified by conventional epidemiology make a substantial contribution to cancer burden1. In clear cell renal cell carcinoma, obesity, hypertension and tobacco smoking are risk factors, but they do not explain the geographical variation in its incidence2. Underlying causes can be inferred by sequencing the genomes of cancers from populations with different incidence rates and detecting differences in patterns of somatic mutations. Here we sequenced 962 clear cell renal cell carcinomas from 11 countries with varying incidence. The somatic mutation profiles differed between countries. In Romania, Serbia and Thailand, mutational signatures characteristic of aristolochic acid compounds were present in most cases, but these were rare elsewhere. In Japan, a mutational signature of unknown cause was found in more than 70% of cases but in less than 2% elsewhere. A further mutational signature of unknown cause was ubiquitous but exhibited higher mutation loads in countries with higher incidence rates of kidney cancer. Known signatures of tobacco smoking correlated with tobacco consumption, but no signature was associated with obesity or hypertension, suggesting that non-mutagenic mechanisms of action underlie these risk factors. The results of this study indicate the existence of multiple, geographically variable, mutagenic exposures that potentially affect tens of millions of people and illustrate the opportunities for new insights into cancer causation through large-scale global cancer genomics.

Genome-wide association meta-analysis identifies pleiotropic risk loci for aerodigestive squamous cell cancers.

Squamous cell carcinomas (SqCC) of the aerodigestive tract have similar etiological risk factors. Although genetic risk variants for individual cancers have been identified, an agnostic, genome-wide search for shared genetic susceptibility has not been performed. To identify novel and pleotropic SqCC risk variants, we performed a meta-analysis of GWAS data on lung SqCC (LuSqCC), oro/pharyngeal SqCC (OSqCC), laryngeal SqCC (LaSqCC) and esophageal SqCC (ESqCC) cancers, totaling 13,887 cases and 61,961 controls of European ancestry. We identified one novel genome-wide significant (Pmeta<5x10-8) aerodigestive SqCC susceptibility loci in the 2q33.1 region (rs56321285, TMEM273). Additionally, three previously unknown loci reached suggestive significance (Pmeta<5x10-7): 1q32.1 (rs12133735, near MDM4), 5q31.2 (rs13181561, TMEM173) and 19p13.11 (rs61494113, ABHD8). Multiple previously identified loci for aerodigestive SqCC also showed evidence of pleiotropy in at least another SqCC site, these include: 4q23 (ADH1B), 6p21.33 (STK19), 6p21.32 (HLA-DQB1), 9p21.33 (CDKN2B-AS1) and 13q13.1(BRCA2). Gene-based association and gene set enrichment identified a set of 48 SqCC-related genes rel to DNA damage and epigenetic regulation pathways. Our study highlights the importance of cross-cancer analyses to identify pleiotropic risk loci of histology-related cancers arising at distinct anatomical sites.

Evaluation of a clinical pharmacogenetics model to predict methotrexate response in patients with rheumatoid arthritis.

Variability of response to treatment hinders successful management of rheumatoid arthritis (RA). Consequently, a clinical pharmacogenetics model for predicting response to methotrexate (CP-MTX) has been previously proposed that includes four clinical variables (disease activity, sex, the presence of rheumatoid factor and smoking status) and four SNPs (rs2236225, rs17602729, rs1127354, and rs2372536) in genes of the folate pathway. It showed good performance, but failed to attract attention, likely, in relation with lack of clear clinical benefit. Here, we have revised the value of the CP-MTX model directly addressing its clinical benefit by focusing on the expected benefit-cost of the predictions. In addition, our study included a much larger number of RA patients (n = 720) in MTX monotherapy than previous studies. Benefit of CP-MTX prediction was defined as the patients that would have received combination therapy as first treatment because they were correctly predicted as non-responders to MTX monotherapy. In contrast, cost of CP-MTX prediction was defined as the responder patients that were wrongly predicted as non-responders. Application of CP-MTX predictions to our patients showed a good benefit-cost relationship, with half of the 66.7% non-responders to MTX monotherapy rightly directed to alternative treatments (a benefit of 33.3%) at the cost of 8.5% wrongly predicted non-responders. These benefits-costs were consistent with reanalysis of the previously published studies. Therefore, predictions of CP-MTX showed a good benefit-cost relationship for informing MTX prescription.

Multi-ancestry genome-wide association study of kidney cancer identifies 63 susceptibility regions.

Here, in a multi-ancestry genome-wide association study meta-analysis of kidney cancer (29,020 cases and 835,670 controls), we identified 63 susceptibility regions (50 novel) containing 108 independent risk loci. In analyses stratified by subtype, 52 regions (78 loci) were associated with clear cell renal cell carcinoma (RCC) and 6 regions (7 loci) with papillary RCC. Notably, we report a variant common in African ancestry individuals ( rs7629500 ) in the 3' untranslated region of VHL, nearly tripling clear cell RCC risk (odds ratio 2.72, 95% confidence interval 2.23-3.30). In cis-expression quantitative trait locus analyses, 48 variants from 34 regions point toward 83 candidate genes. Enrichment of hypoxia-inducible factor-binding sites underscores the importance of hypoxia-related mechanisms in kidney cancer. Our results advance understanding of the genetic architecture of kidney cancer, provide clues for functional investigation and enable generation of a validated polygenic risk score with an estimated area under the curve of 0.65 (0.74 including risk factors) among European ancestry individuals.

The Mutographs biorepository: A unique genomic resource to study cancer around the world.

Large-scale biorepositories and databases are essential to generate equitable, effective, and sustainable advances in cancer prevention, early detection, cancer therapy, cancer care, and surveillance. The Mutographs project has created a large genomic dataset and biorepository of over 7,800 cancer cases from 30 countries across five continents with extensive demographic, lifestyle, environmental, and clinical information. Whole-genome sequencing is being finalized for over 4,000 cases, with the primary goal of understanding the causes of cancer at eight anatomic sites. Genomic, exposure, and clinical data will be publicly available through the International Cancer Genome Consortium Accelerating Research in Genomic Oncology platform. The Mutographs sample and metadata biorepository constitutes a legacy resource for new projects and collaborations aiming to increase our current research efforts in cancer genomic epidemiology globally.

Penile Cancers Attributed to Human Papillomavirus Are Associated with Improved Survival for Node-positive Patients. Findings from a Norwegian Cohort Study Spanning 50 Years.

BACKGROUND: Human papillomavirus (HPV) infection is a risk factor for the development of penile squamous cell carcinoma (PSCC). It remains inconclusive whether HPV-related PSCC has a different prognosis from non-HPV-related PSCC. OBJECTIVE: To investigate the relationship between HPV status and survival as well as temporal changes in the proportion of HPV-related PSCC. DESIGN, SETTING, AND PARTICIPANTS: A retrospective cohort of 277 patients treated in Norway between 1973 and 2022 was investigated for HPV DNA in tumor tissue. Clinicopathological variables and disease course were registered. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Kaplan-Meier curves and Cox regression were used to investigate the determinants of cancer-specific survival (CSS). The chi-square test for trend in proportions enabled investigation of temporal changes in the HPV-related proportion of PSCC patients treated in Western Norway (n = 211). RESULTS AND LIMITATIONS: HPV DNA was detected in tumor tissue from 131 (47%) patients. Stratified by HPV status, 5-yr CSS did not differ between groups (p = 0.37). When investigating only node-positive patients, however, presence of HPV DNA was an independent predictor of better survival in multivariable Cox regression after adjustment for age, nodal stage, and adjuvant therapy (hazard ratio 0.54, 95% confidence interval: [0.30-0.99], p = 0.04). In cases from Western Norway, an increasing proportion of HPV-related cases over time was found (p = 0.01). The main limitation is the retrospective study design. CONCLUSIONS: HPV DNA in tumor tissue was associated with significantly better CSS for node-positive patients. The proportion of HPV DNA-positive PSCC has increased significantly in Western Norway over the past 50 yr. PATIENT SUMMARY: We investigated the impact of human papillomavirus (HPV) on the survival of penile cancer patients treated over a 50-yr period in Norway. We found that for patients with lymph node metastasis, survival was better for HPV-related cases. We also found that the proportion of cases due to HPV has increased in Western Norway.

Nasopharyngeal carcinoma patients from Norway show elevated Epstein-Barr virus IgA and IgG antibodies prior to diagnosis.

BACKGROUND: IgA antibodies against few Epstein-Barr virus (EBV) proteins are established serological markers for nasopharyngeal carcinoma (NPC). We recently validated a novel, comprehensive EBV marker panel and showed that IgA, but also IgG antibodies against multiple EBV proteins are highly sensitive and specific for EBV-positive NPC at diagnosis. However, data about these novel biomarkers as prospective markers for NPC are sparse. METHODS: This study included 30 incident NPC cases and 60 matched controls from the Norwegian Janus Serum Bank. For 21 NPCs, molecular EBV and human papillomavirus (HPV) status were assessed by EBER-ISH and HPV DNA/RNA testing by PCR, respectively. IgA and IgG serum antibodies against 17 EBV antigens were analyzed in prediagnostic sera of cases (median lead time 14 years) and controls using multiplex serology. Sensitivities were calculated using receiver operating characteristic analysis pre-specified to yield 90% specificity in the control group. From 10 cases, serial samples were available. RESULTS: Quantitative EBV antibody levels were significantly elevated among all cases (p < 0.05) for three IgA and six IgG antibodies. The highest sensitivities for defining 12 EBER-ISH-positive NPCs were observed for BGLF2 IgA (67%) and BGLF2 IgG (83%). Increased IgA and IgG antibody levels between the first and last draw before diagnosis were observed for EBER-ISH positive, but not for EBER-ISH negative NPCs. Among 21 molecularly analyzed NPCs, 4 EBER-ISH negative NPCs showed concomitant positivity to HPV type-specific DNA and RNA; 3 NPCs were HPV16 and 1 NPC was HPV18 positive. CONCLUSION: Both, EBV IgA and IgG antibody levels are significantly elevated many years before diagnosis of EBV-positive NPCs in Norway, an NPC low-incidence region. This study provides insights into one of the largest available prospective sample collections of NPCs in a non-endemic country.

Absolute Risk of Oropharyngeal Cancer After an HPV16-E6 Serology Test and Potential Implications for Screening: Results From the Human Papillomavirus Cancer Cohort Consortium.

PURPOSE: Seropositivity for the HPV16-E6 oncoprotein is a promising marker for early detection of oropharyngeal cancer (OPC), but the absolute risk of OPC after a positive or negative test is unknown. METHODS: We constructed an OPC risk prediction model that integrates (1) relative odds of OPC for HPV16-E6 serostatus and cigarette smoking from the human papillomavirus (HPV) Cancer Cohort Consortium (HPVC3), (2) US population risk factor data from the National Health Interview Survey, and (3) US sex-specific population rates of OPC and mortality. RESULTS: The nine HPVC3 cohorts included 365 participants with OPC with up to 10 years between blood draw and diagnosis and 5,794 controls. The estimated 10-year OPC risk for HPV16-E6 seropositive males at age 50 years was 17.4% (95% CI, 12.4 to 28.6) and at age 60 years was 27.1% (95% CI, 19.2 to 45.4). Corresponding 5-year risk estimates were 7.3% and 14.4%, respectively. For HPV16-E6 seropositive females, 10-year risk estimates were 3.6% (95% CI, 2.5 to 5.9) at age 50 years and 5.5% (95% CI, 3.8 to 9.2) at age 60 years and 5-year risk estimates were 1.5% and 2.7%, respectively. Over 30 years, after a seropositive result at age 50 years, an estimated 49.9% of males and 13.3% of females would develop OPC. By contrast, 10-year risks among HPV16-E6 seronegative people were very low, ranging from 0.01% to 0.25% depending on age, sex, and smoking status. CONCLUSION: We estimate that a substantial proportion of HPV16-E6 seropositive individuals will develop OPC, with 10-year risks of 17%-27% for males and 4%-6% for females age 50-60 years in the United States. This high level of risk may warrant periodic, minimally invasive surveillance after a positive HPV16-E6 serology test, particularly for males in high-incidence regions. However, an appropriate clinical protocol for surveillance remains to be established.

Germline determinants of humoral immune response to HPV-16 protect against oropharyngeal cancer.

Although several oropharyngeal cancer (OPC) susceptibility loci have been identified, most previous studies lacked detailed information on human papillomavirus (HPV) status. We conduct a genome-wide analysis by HPV16 serology status in 4,002 oral cancer cases (OPC and oral cavity cancer (OCC)) and 5,256 controls. We detect four susceptibility loci pointing to a distinct genetic predisposition by HPV status. Our most notable finding in the HLA region, that is now confirmed to be specific of HPV(+)OPC risk, reveal two independent loci with strong protective effects, one refining the previously reported HLA class II haplotype association. Antibody levels against HPV16 viral proteins strongly implicate the protective HLA variants as major determinants of humoral response against L1 capsid protein or E6 oncoprotein suggesting a natural immune response against HPV(+)OPC promoted by HLA variants. This indicates that therapeutic vaccines that target E6 and attenuate viral response after established HPV infections might protect against HPV(+)OPC.

Evaluation of 12 GWAS-drawn SNPs as biomarkers of rheumatoid arthritis response to TNF inhibitors. A potential SNP association with response to etanercept.

Research in rheumatoid arthritis (RA) is increasingly focused on the discovery of biomarkers that could enable personalized treatments. The genetic biomarkers associated with the response to TNF inhibitors (TNFi) are among the most studied. They include 12 SNPs exhibiting promising results in the three largest genome-wide association studies (GWAS). However, they still require further validation. With this aim, we assessed their association with response to TNFi in a replication study, and a meta-analysis summarizing all non-redundant data. The replication involved 755 patients with RA that were treated for the first time with a biologic drug, which was either infliximab (n = 397), etanercept (n = 155) or adalimumab (n = 203). Their DNA samples were successfully genotyped with a single-base extension multiplex method. Lamentably, none of the 12 SNPs was associated with response to the TNFi in the replication study (p > 0.05). However, a drug-stratified exploratory analysis revealed a significant association of the NUBPL rs2378945 SNP with a poor response to etanercept (B = -0.50, 95% CI = -0.82, -0.17, p = 0.003). In addition, the meta-analysis reinforced the previous association of three SNPs: rs2378945, rs12142623, and rs4651370. In contrast, five of the remaining SNPs were less associated than before, and the other four SNPs were no longer associated with the response to treatment. In summary, our results highlight the complexity of the pharmacogenetics of TNFi in RA showing that it could involve a drug-specific component and clarifying the status of the 12 GWAS-drawn SNPs.

Replication study of polymorphisms associated with response to methotrexate in patients with rheumatoid arthritis.

About 70 genetic studies have already addressed the need of biomarkers to predict the response of patients with rheumatoid arthritis (RA) to methotrexate (MTX) treatment. However, no genetic biomarker has yet been sufficiently validated. Here, we aimed to replicate a selection of 25 SNPs in the largest collection of patients up to date, which consisted of 915 patients treated with MTX. The change in disease activity (measured as ΔDAS28) from baseline was considered the primary outcome. In addition, response according to widely used criteria (EULAR) was taken as secondary outcome. We considered consistency between outcomes, P values accounting for the number of SNPs, and independence from potential confounders for interpretation of the results. Only the rs1801394 SNP in MTRR fulfilled the high association standards. Its minor allele was associated with less improvement than the major allele according to ΔDAS28 (p = 0.0016), and EULAR response (p = 0.004), with independence of sex, age, baseline DAS28, smoking, seropositivity, concomitant corticosteroid use or previous treatments. In addition, previous evidence suggests the association of this SNP with response to MTX in another autoimmune disease, juvenile idiopathic arthritis, and with high intracellular folate levels, which could contribute to poor response.

Fine mapping of MHC region in lung cancer highlights independent susceptibility loci by ethnicity.

Lung cancer has several genetic associations identified within the major histocompatibility complex (MHC); although the basis for these associations remains elusive. Here, we analyze MHC genetic variation among 26,044 lung cancer patients and 20,836 controls densely genotyped across the MHC, using the Illumina Illumina OncoArray or Illumina 660W SNP microarray. We impute sequence variation in classical HLA genes, fine-map MHC associations for lung cancer risk with major histologies and compare results between ethnicities. Independent and novel associations within HLA genes are identified in Europeans including amino acids in the HLA-B*0801 peptide binding groove and an independent HLA-DQB1*06 loci group. In Asians, associations are driven by two independent HLA allele sets that both increase risk in HLA-DQB1*0401 and HLA-DRB1*0701; the latter better represented by the amino acid Ala-104. These results implicate several HLA-tumor peptide interactions as the major MHC factor modulating lung cancer susceptibility.

Lack of replication of interactions between polymorphisms in rheumatoid arthritis susceptibility: case-control study.

INTRODUCTION: Approximately 100 loci have been definitively associated with rheumatoid arthritis (RA) susceptibility. However, they explain only a fraction of RA heritability. Interactions between polymorphisms could explain part of the remaining heritability. Multiple interactions have been reported, but only the shared epitope (SE) × protein tyrosine phosphatase nonreceptor type 22 (PTPN22) interaction has been replicated convincingly. Two recent studies deserve attention because of their quality, including their replication in a second sample collection. In one of them, researchers identified interactions between PTPN22 and seven single-nucleotide polymorphisms (SNPs). The other showed interactions between the SE and the null genotype of glutathione S-transferase Mu 1 (GSTM1) in the anti-cyclic citrullinated peptide-positive (anti-CCP+) patients. In the present study, we aimed to replicate association with RA susceptibility of interactions described in these two high-quality studies. METHODS: A total of 1,744 patients with RA and 1,650 healthy controls of Spanish ancestry were studied. Polymorphisms were genotyped by single-base extension. SE genotypes of 736 patients were available from previous studies. Interaction analysis was done using multiple methods, including those originally reported and the most powerful methods described. RESULTS: Genotypes of one of the SNPs (rs4695888) failed quality control tests. The call rate for the other eight polymorphisms was 99.9%. The frequencies of the polymorphisms were similar in RA patients and controls, except for PTPN22 SNP. None of the interactions between PTPN22 SNPs and the six SNPs that met quality control tests was replicated as a significant interaction term--the originally reported finding--or with any of the other methods. Nor was the interaction between GSTM1 and the SE replicated as a departure from additivity in anti-CCP+ patients or with any of the other methods. CONCLUSIONS: None of the interactions tested were replicated in spite of sufficient power and assessment with different assays. These negative results indicate that whether interactions are significant contributors to RA susceptibility remains unknown and that strict standards need to be applied to claim that an interaction exists.

Lack of validation of genetic variants associated with anti-tumor necrosis factor therapy response in rheumatoid arthritis: a genome-wide association study replication and meta-analysis.

INTRODUCTION: In this study, our aim was to elucidate the role of four polymorphisms identified in a prior large genome-wide association study (GWAS) in which the investigators analyzed the responses of patients with rheumatoid arthritis (RA) to treatment with tumor necrosis factor inhibitors (TNFi). The authors of that study reported that the four genetic variants were significantly associated. However, none of the associations reached GWAS significance, and two subsequent studies failed to replicate these associations. METHODS: The four polymorphisms (rs12081765, rs1532269, rs17301249 and rs7305646) were genotyped in a total of 634 TNFi-treated RA patients of Spanish Caucasian origin. Four outcomes were evaluated: changes in the Disease Activity Score in 28 joints (DAS28) after 6 and 12 months of treatment and classification according to the European League Against Rheumatism (EULAR) response criteria at the same time points. Association with DAS28 changes was assessed by linear regression using an additive genetic model. Contingency tables of genotype and allele frequencies between EULAR responder and nonresponder patients were compared. In addition, we combined our data with those of previously reported studies in a meta-analysis including 2,998 RA patients. RESULTS: None of the four genetic variants showed an association with response to TNFi in any of the four outcomes analyzed in our Spanish patients. In addition, only rs1532269 yielded a suggestive association (P = 0.0033) with the response to TNFi when available data from previous studies were combined in the meta-analysis. CONCLUSION: Our data suggest that the rs12081765, rs1532269, rs17301249 and rs7305646 genetic variants do not have a role as genetic predictors of TNFi treatment outcomes.

Ultrasonographic assessment of enthesitis in HLA-B27 positive patients with rheumatoid arthritis, a matched case-only study.

INTRODUCTION: HLA-B27 has a modifier effect on the phenotype of multiple diseases, both associated and non-associated with it. Among these effects, an increased frequency of clinical enthesitis in patients with Rheumatoid Arthritis (RA) has been reported but never explored again. We aimed to replicate this study with a sensitive and quantitative assessment of enthesitis by using standardized ultrasonography (US). METHODS: The Madrid Sonography Enthesitis Index (MASEI) was applied to the US assessment of 41 HLA-B27 positive and 41 matched HLA-B27 negative patients with longstanding RA. Clinical characteristics including explorations aimed to evaluate spondyloarthrtitis and laboratory tests were also done. RESULTS: A significant degree of abnormalities in the entheses of the patients with RA were found, but the MASEI values, and each of its components including the Doppler signal, were similar in HLA-B27 positive and negative patients. An increase of the MASEI scores with age was identified. Differences in two clinical features were found: a lower prevalence of rheumatoid factor and a more common story of low back pain in the HLA-B27 positive patients than in the negative. The latter was accompanied by radiographic sacroiliitis in two HLA-B27 positive patients. No other differences were detected. CONCLUSION: We have found that HLA-B27 positive patients with RA do not have more enthesitis as assessed with US than the patients lacking this HLA allele. However, HLA-B27 could be shaping the RA phenotype towards RF seronegativity and axial involvement.

Variants within STAT genes reveal association with anticitrullinated protein antibody-negative rheumatoid arthritis in 2 European populations.

OBJECTIVE: STAT3 and 4 are, among other factors, critical for the interleukin 12 (IL-12)-mediated Th1 response, for transfer of IL-23 signals, and for survival and expansion of Th17 cells. We investigated the association of STAT3 and STAT4 polymorphisms with serologically distinct subgroups of rheumatoid arthritis (RA). METHODS: A total of 41 single-nucleotide polymorphisms (SNP) within STAT3 and STAT1-STAT4 loci were investigated in a Swedish cohort of 2043 RA cases and 1115 controls. Nine of the associated SNP were tested in a Spanish cohort of 1223 RA cases and 1090 controls. RESULTS: Fourteen SNP in the STAT3 and STAT1-STAT4 loci were associated with anticitrullinated protein antibody (ACPA)-negative RA in the Swedish cohort. Three of the SNP in STAT4 and 2 SNP in STAT3 remained associated with ACPA-negative RA after considering the Spanish results. In addition, rs7574865 and rs10181656, in STAT4, were associated with ACPA-positive RA in the Swedish study. One of these SNP, rs7574865, showed a similar pattern of the association in serologically distinct subgroups of RA in a metaanalysis of all 7 published studies. CONCLUSION: Our findings suggest that variants in STAT genes may contribute differentially to susceptibility to RA in seropositive and in seronegative patients.