Novel genes and sex differences in COVID-19 severity.

Here, we describe the results of a genome-wide study conducted in 11 939 coronavirus disease 2019 (COVID-19) positive cases with an extensive clinical information that were recruited from 34 hospitals across Spain (SCOURGE consortium). In sex-disaggregated genome-wide association studies for COVID-19 hospitalization, genome-wide significance (P < 5 × 10-8) was crossed for variants in 3p21.31 and 21q22.11 loci only among males (P = 1.3 × 10-22 and P = 8.1 × 10-12, respectively), and for variants in 9q21.32 near TLE1 only among females (P = 4.4 × 10-8). In a second phase, results were combined with an independent Spanish cohort (1598 COVID-19 cases and 1068 population controls), revealing in the overall analysis two novel risk loci in 9p13.3 and 19q13.12, with fine-mapping prioritized variants functionally associated with AQP3 (P = 2.7 × 10-8) and ARHGAP33 (P = 1.3 × 10-8), respectively. The meta-analysis of both phases with four European studies stratified by sex from the Host Genetics Initiative (HGI) confirmed the association of the 3p21.31 and 21q22.11 loci predominantly in males and replicated a recently reported variant in 11p13 (ELF5, P = 4.1 × 10-8). Six of the COVID-19 HGI discovered loci were replicated and an HGI-based genetic risk score predicted the severity strata in SCOURGE. We also found more SNP-heritability and larger heritability differences by age (<60 or ≥60 years) among males than among females. Parallel genome-wide screening of inbreeding depression in SCOURGE also showed an effect of homozygosity in COVID-19 hospitalization and severity and this effect was stronger among older males. In summary, new candidate genes for COVID-19 severity and evidence supporting genetic disparities among sexes are provided.

A crowdsourcing database for the copy-number variation of the Spanish population.

BACKGROUND: Despite being a very common type of genetic variation, the distribution of copy-number variations (CNVs) in the population is still poorly understood. The knowledge of the genetic variability, especially at the level of the local population, is a critical factor for distinguishing pathogenic from non-pathogenic variation in the discovery of new disease variants. RESULTS: Here, we present the SPAnish Copy Number Alterations Collaborative Server (SPACNACS), which currently contains copy number variation profiles obtained from more than 400 genomes and exomes of unrelated Spanish individuals. By means of a collaborative crowdsourcing effort whole genome and whole exome sequencing data, produced by local genomic projects and for other purposes, is continuously collected. Once checked both, the Spanish ancestry and the lack of kinship with other individuals in the SPACNACS, the CNVs are inferred for these sequences and they are used to populate the database. A web interface allows querying the database with different filters that include ICD10 upper categories. This allows discarding samples from the disease under study and obtaining pseudo-control CNV profiles from the local population. We also show here additional studies on the local impact of CNVs in some phenotypes and on pharmacogenomic variants. SPACNACS can be accessed at: http://csvs.clinbioinfosspa.es/spacnacs/ . CONCLUSION: SPACNACS facilitates disease gene discovery by providing detailed information of the local variability of the population and exemplifies how to reuse genomic data produced for other purposes to build a local reference database.

CSVS, a crowdsourcing database of the Spanish population genetic variability.

The knowledge of the genetic variability of the local population is of utmost importance in personalized medicine and has been revealed as a critical factor for the discovery of new disease variants. Here, we present the Collaborative Spanish Variability Server (CSVS), which currently contains more than 2000 genomes and exomes of unrelated Spanish individuals. This database has been generated in a collaborative crowdsourcing effort collecting sequencing data produced by local genomic projects and for other purposes. Sequences have been grouped by ICD10 upper categories. A web interface allows querying the database removing one or more ICD10 categories. In this way, aggregated counts of allele frequencies of the pseudo-control Spanish population can be obtained for diseases belonging to the category removed. Interestingly, in addition to pseudo-control studies, some population studies can be made, as, for example, prevalence of pharmacogenomic variants, etc. In addition, this genomic data has been used to define the first Spanish Genome Reference Panel (SGRP1.0) for imputation. This is the first local repository of variability entirely produced by a crowdsourcing effort and constitutes an example for future initiatives to characterize local variability worldwide. CSVS is also part of the GA4GH Beacon network. CSVS can be accessed at: http://csvs.babelomics.org/.

Recurrent Germline DLST Mutations in Individuals with Multiple Pheochromocytomas and Paragangliomas.

Pheochromocytomas and paragangliomas (PPGLs) provide some of the clearest genetic evidence for the critical role of metabolism in the tumorigenesis process. Approximately 40% of PPGLs are caused by driver germline mutations in 16 known susceptibility genes, and approximately half of these genes encode members of the tricarboxylic acid (TCA) cycle. Taking as a starting point the involvement of the TCA cycle in PPGL development, we aimed to identify unreported mutations that occurred in genes involved in this key metabolic pathway and that could explain the phenotypes of additional individuals who lack mutations in known susceptibility genes. To accomplish this, we applied a targeted sequencing of 37 TCA-cycle-related genes to DNA from 104 PPGL-affected individuals with no mutations in the major known predisposing genes. We also performed omics-based analyses, TCA-related metabolite determination, and 13C5-glutamate labeling assays. We identified five germline variants affecting DLST in eight unrelated individuals (∼7%); all except one were diagnosed with multiple PPGLs. A recurrent variant, c.1121G>A (p.Gly374Glu), found in four of the eight individuals triggered accumulation of 2-hydroxyglutarate, both in tumors and in a heterologous cell-based assay designed to functionally evaluate DLST variants. p.Gly374Glu-DLST tumors exhibited loss of heterozygosity, and their methylation and expression profiles are similar to those of EPAS1-mutated PPGLs; this similarity suggests a link between DLST disruption and pseudohypoxia. Moreover, we found positive DLST immunostaining exclusively in tumors carrying TCA-cycle or EPAS1 mutations. In summary, this study reveals DLST as a PPGL-susceptibility gene and further strengthens the relevance of the TCA cycle in PPGL development.

Human Genomic Diversity Where the Mediterranean Joins the Atlantic.

Throughout the past few years, a lively debate emerged about the timing and magnitude of the human migrations between the Iberian Peninsula and the Maghreb. Several pieces of evidence, including archaeological, anthropological, historical, and genetic data, have pointed to a complex and intermingled evolutionary history in the western Mediterranean area. To study to what extent connections across the Strait of Gibraltar and surrounding areas have shaped the present-day genomic diversity of its populations, we have performed a screening of 2.5 million single-nucleotide polymorphisms in 142 samples from southern Spain, southern Portugal, and Morocco. We built comprehensive data sets of the studied area and we implemented multistep bioinformatic approaches to assess population structure, demographic histories, and admixture dynamics. Both local and global ancestry inference showed an internal substructure in the Iberian Peninsula, mainly linked to a differential African ancestry. Western Iberia, from southern Portugal to Galicia, constituted an independent cluster within Iberia characterized by an enriched African genomic input. Migration time modeling showed recent historic dates for the admixture events occurring both in Iberia and in the North of Africa. However, an integrative vision of both paleogenomic and modern DNA data allowed us to detect chronological transitions and population turnovers that could be the result of transcontinental migrations dating back from Neolithic times. The present contribution aimed to fill the gaps in the modern human genomic record of a key geographic area, where the Mediterranean and the Atlantic come together.

Hsa-miR-139-5p is a prognostic thyroid cancer marker involved in HNRNPF-mediated alternative splicing.

It is critical to identify biomarkers and functional networks associated with aggressive thyroid cancer to anticipate disease progression and facilitate personalized patient management. We performed miRNome sequencing of 46 thyroid tumors enriched with advanced disease patients with a median follow-up of 96 months. MiRNome profiles correlated with tumor-specific histopathological and molecular features, such as stromal cell infiltration and tumor driver mutation. Differential expression analysis revealed a consistent hsa-miR-139-5p downexpression in primary carcinomas from patients with recurrent/metastatic disease compared to disease-free patients, sustained in paired local metastases and validated in publicly available thyroid cancer series. Exogenous expression of hsa-miR-139-5p significantly reduced migration and proliferation of anaplastic thyroid cancer cells. Proteomic analysis indicated RICTOR, SMAD2/3 and HNRNPF as putative hsa-miR-139-5p targets in our cell system. Abundance of HNRNPF mRNA, encoding an alternative splicing factor involved in cryptic exon inclusion/exclusion, inversely correlated with hsa-miR-139-5p expression in human tumors. RNA sequencing analysis revealed 174 splicing events differentially regulated upon HNRNPF repression in our cell system, affecting genes involved in RTK/RAS/MAPK and PI3K/AKT/MTOR signaling cascades among others. These results point at the hsa-miR-139-5p/HNRNPF axis as a novel regulatory mechanism associated with the modulation of major thyroid cancer signaling pathways and tumor virulence.

RECQL5: Another DNA helicase potentially involved in hereditary breast cancer susceptibility.

There is still around 50% of the familial breast cancer (BC) cases with an undefined genetic cause, here we have used next-generation sequencing (NGS) technology to identify new BC susceptibility genes. This approach has led to the identification of RECQL5, a member of RECQL-helicases family, as a new BC susceptibility candidate, which deserves further study. We have used a combination of whole exome sequencing in a family negative for mutations in BRCA1/2 throughout (BRCAX), in which we found a probably deleterious variant in RECQL5, and targeted NGS of the complete coding regions and exon-intron boundaries of the candidate gene in 699 BC Spanish BRCAX families and 665 controls. Functional characterization and in silico inference of pathogenicity were performed to evaluate the deleterious effect of detected variants. We found at least seven deleterious or likely deleterious variants among the cases and only one in controls. These results prompt us to propose RECQL5 as a gene that would be worth to analyze in larger studies to explore its possible implication in BC susceptibility.

Evidence that the 5p12 Variant rs10941679 Confers Susceptibility to Estrogen-Receptor-Positive Breast Cancer through FGF10 and MRPS30 Regulation.

Genome-wide association studies (GWASs) have revealed increased breast cancer risk associated with multiple genetic variants at 5p12. Here, we report the fine mapping of this locus using data from 104,660 subjects from 50 case-control studies in the Breast Cancer Association Consortium (BCAC). With data for 3,365 genotyped and imputed SNPs across a 1 Mb region (positions 44,394,495-45,364,167; NCBI build 37), we found evidence for at least three independent signals: the strongest signal, consisting of a single SNP rs10941679, was associated with risk of estrogen-receptor-positive (ER+) breast cancer (per-g allele OR ER+ = 1.15; 95% CI 1.13-1.18; p = 8.35 × 10-30). After adjustment for rs10941679, we detected signal 2, consisting of 38 SNPs more strongly associated with ER-negative (ER-) breast cancer (lead SNP rs6864776: per-a allele OR ER- = 1.10; 95% CI 1.05-1.14; p conditional = 1.44 × 10-12), and a single signal 3 SNP (rs200229088: per-t allele OR ER+ = 1.12; 95% CI 1.09-1.15; p conditional = 1.12 × 10-05). Expression quantitative trait locus analysis in normal breast tissues and breast tumors showed that the g (risk) allele of rs10941679 was associated with increased expression of FGF10 and MRPS30. Functional assays demonstrated that SNP rs10941679 maps to an enhancer element that physically interacts with the FGF10 and MRPS30 promoter regions in breast cancer cell lines. FGF10 is an oncogene that binds to FGFR2 and is overexpressed in ∼10% of human breast cancers, whereas MRPS30 plays a key role in apoptosis. These data suggest that the strongest signal of association at 5p12 is mediated through coordinated activation of FGF10 and MRPS30, two candidate genes for breast cancer pathogenesis.

Whole exome sequencing identifies PLEC, EXO5 and DNAH7 as novel susceptibility genes in testicular cancer.

Testicular germ cell tumors (TGCTs) are a clinically and pathologically heterogeneous disease, and little is known of its genetic basis. Only low susceptibility risk loci have been identified for both sporadic and familial cases. Therefore, we tried to identify new susceptibility genes responsible for familial testicular cancer that may contribute to increasing our knowledge about the genetic basis of the disease. Nineteen Spanish families with at least two affected individuals with TGCT were selected. WES was performed on those individuals using an Illumina Hiseq2000 sequencing platform. Data were analyzed under a monogenic and polygenic model of inheritance, and candidate variants were evaluated in a case-control association study performed on 391 Spanish sporadic cases and 1,170 healthy Spanish controls. Results were replicated in a second series consisting of 101 TGCTs from the Cancer Genome Atlas (TGCA) and 27,000 controls from the Exome Aggregation Consortium (ExAC) database. Logistic regression was carried out to analyze the association strength (risk) of candidate variants obtained among cases and controls in different populations. Despite the sample size, we detected a significant earlier age of onset in familial TGCT (28y) than sporadic cases (33y), using a Mann-Whitney U test. We identified significant variants in the comparative study of TGCT cases (391) versus controls (almost 1,170), and three of them [PLEC (OR = 6.28, p = 6.42 × 10-23 ) (p.Arg2016Trp), EXO5 (OR = 3.37, p = 4.82 × 10-09 ) (p.Arg344AlafsTer10) and DNAH7 (OR = 1.64, p = 0.048)] were replicated as potential candidates that may contribute to explaining the genetic basis of TGCT.

Pharmacogenetic variants and response to neoadjuvant single-agent doxorubicin or docetaxel: a study in locally advanced breast cancer patients participating in the NCT00123929 phase 2 randomized trial.

OBJECTIVES: Taxanes and anthracyclines are widely used in the treatment of breast cancer, although the benefit is limited to a proportion of patients and predictive biomarkers for clinical outcome remain elusive. PATIENTS AND METHODS: We carried out a pharmacogenetic study in 181 patients with locally advanced breast cancer enrolled in a phase 2 randomized clinical trial (NCT00123929), where patients were randomly assigned to receive neoadjuvant single-agent docetaxel 100 mg/m(2) (n=84) or doxorubicin 75 mg/m(2) (n=97). We studied the association of 226 single nucleotide polymorphisms (SNPs) in 15 key drug biotransformation genes with neoadjuvant pathological tumor response residual cancer burden index to docetaxel and to doxorubicin. RESULTS: We identified a significant association for rs162561, an intronic SNP located in the cytochrome P450 family 1 subfamily B member 1 (CYP1B1) gene, with tumor response in patients treated with single-agent docetaxel (dominant model: ß=1.02, 95% confidence interval=0.49-1.55; P=1.77×10(-4)), and for rs717620, an SNP located in the promoter of the ATP-binding cassette subfamily C member 2 (ABCC2) gene, in patients treated with neoadjuvant doxorubicin (recessive model: ß=1.67; 95% confidence interval=0.26-3.11; P=0.02). CONCLUSION: We identified two polymorphisms in CYP1B1 and ABCC2 associated with tumor pathological response following docetaxel or doxorubicin neoadjuvant monotherapy, respectively. Although further validation is required, these variants could be potential predictive genetic markers for treatment outcome in breast cancer patients.

Exome array analysis identifies ETFB as a novel susceptibility gene for anthracycline-induced cardiotoxicity in cancer patients.

PURPOSE: Anthracyclines are widely used chemotherapeutic drugs that can cause progressive and irreversible cardiac damage and fatal heart failure. Several genetic variants associated with anthracycline-induced cardiotoxicity (AIC) have been identified, but they explain only a small proportion of the interindividual differences in AIC susceptibility. METHODS: In this study, we evaluated the association of low-frequency variants with risk of chronic AIC using the Illumina HumanExome BeadChip array in a discovery cohort of 61 anthracycline-treated breast cancer patients with replication in a second independent cohort of 83 anthracycline-treated pediatric cancer patients, using gene-based tests (SKAT-O). RESULTS: The most significant associated gene in the discovery cohort was ETFB (electron transfer flavoprotein beta subunit) involved in mitochondrial ß-oxidation and ATP production (P = 4.16 × 10-4) and this association was replicated in an independent set of anthracycline-treated cancer patients (P = 2.81 × 10-3). Within ETFB, we found that the missense variant rs79338777 (p.Pro52Leu; c.155C > T) made the greatest contribution to the observed gene association and it was associated with increased risk of chronic AIC in the two cohorts separately and when combined (OR 9.00, P = 1.95 × 10-4, 95% CI 2.83-28.6). CONCLUSIONS: We identified and replicated a novel gene, ETFB, strongly associated with chronic AIC independently of age at tumor onset and related to anthracycline-mediated mitochondrial dysfunction. Although experimental verification and further studies in larger patient cohorts are required to confirm our finding, we demonstrated that exome array data analysis represents a valuable strategy to identify novel genes contributing to the susceptibility to chronic AIC.

Gain-of-function mutations in DNMT3A in patients with paraganglioma.

PURPOSE: The high percentage of patients carrying germline mutations makes pheochromocytomas/paragangliomas the most heritable of all tumors. However, there are still cases unexplained by mutations in the known genes. We aimed to identify the genetic cause of disease in patients strongly suspected of having hereditary tumors. METHODS: Whole-exome sequencing was applied to the germlines of a parent-proband trio. Genome-wide methylome analysis, RNA-seq, CRISPR/Cas9 gene editing, and targeted sequencing were also performed. RESULTS: We identified a novel de novo germline mutation in DNMT3A, affecting a highly conserved residue located close to the aromatic cage that binds to trimethylated histone H3. DNMT3A-mutated tumors exhibited significant hypermethylation of homeobox-containing genes, suggesting an activating role of the mutation. CRISPR/Cas9-mediated knock-in in HeLa cells led to global changes in methylation, providing evidence of the DNMT3A-altered function. Targeted sequencing revealed subclonal somatic mutations in six additional paragangliomas. Finally, a second germline DNMT3A mutation, also causing global tumor DNA hypermethylation, was found in a patient with a family history of pheochromocytoma. CONCLUSION: Our findings suggest that DNMT3A may be a susceptibility gene for paragangliomas and, if confirmed in future studies, would represent the first example of gain-of-function mutations affecting a DNA methyltransferase gene involved in cancer predisposition.

Identification and characterization of novel associations in the CASP8/ALS2CR12 region on chromosome 2 with breast cancer risk.

Previous studies have suggested that polymorphisms in CASP8 on chromosome 2 are associated with breast cancer risk. To clarify the role of CASP8 in breast cancer susceptibility, we carried out dense genotyping of this region in the Breast Cancer Association Consortium (BCAC). Single-nucleotide polymorphisms (SNPs) spanning a 1 Mb region around CASP8 were genotyped in 46 450 breast cancer cases and 42 600 controls of European origin from 41 studies participating in the BCAC as part of a custom genotyping array experiment (iCOGS). Missing genotypes and SNPs were imputed and, after quality exclusions, 501 typed and 1232 imputed SNPs were included in logistic regression models adjusting for study and ancestry principal components. The SNPs retained in the final model were investigated further in data from nine genome-wide association studies (GWAS) comprising in total 10 052 case and 12 575 control subjects. The most significant association signal observed in European subjects was for the imputed intronic SNP rs1830298 in ALS2CR12 (telomeric to CASP8), with per allele odds ratio and 95% confidence interval [OR (95% confidence interval, CI)] for the minor allele of 1.05 (1.03-1.07), P = 1 × 10(-5). Three additional independent signals from intronic SNPs were identified, in CASP8 (rs36043647), ALS2CR11 (rs59278883) and CFLAR (rs7558475). The association with rs1830298 was replicated in the imputed results from the combined GWAS (P = 3 × 10(-6)), yielding a combined OR (95% CI) of 1.06 (1.04-1.08), P = 1 × 10(-9). Analyses of gene expression associations in peripheral blood and normal breast tissue indicate that CASP8 might be the target gene, suggesting a mechanism involving apoptosis.

Exome array analysis identifies GPR35 as a novel susceptibility gene for anthracycline-induced cardiotoxicity in childhood cancer.

OBJECTIVES: Pediatric cancer survivors are a steadily growing population; however, chronic anthracycline-induced cardiotoxicity (AIC) is a serious long-term complication leading to considerable morbidity. We aimed to identify new genes and low-frequency variants influencing the susceptibility to AIC for pediatric cancer patients. PATIENTS AND METHODS: We studied the association of variants on the Illumina HumanExome BeadChip array in 83 anthracycline-treated pediatric cancer patients. In addition to single-variant association tests, we carried out a gene-based analysis to investigate the combined effects of common and low-frequency variants to chronic AIC. RESULTS: Although no single-variant showed an association with chronic AIC that was statistically significant after correction for multiple testing, we identified a novel significant association for G protein-coupled receptor 35 (GPR35) by gene-based testing, a gene with potential roles in cardiac physiology and pathology (P=7.0×10), which remained statistically significant after correction for multiple testing (PFDR=0.03). The greatest contribution to this observed association was made by rs12468485, a missense variant (p.Thr253Met, c.758C>T, minor allele frequency=0.04), with the T allele associated with an increased risk of chronic AIC and more severe symptomatic cardiac manifestations at low anthracycline doses. CONCLUSION: Using exome array data, we identified GPR35 as a novel susceptibility gene associated with chronic AIC in pediatric cancer patients.

DNA glycosylases involved in base excision repair may be associated with cancer risk in BRCA1 and BRCA2 mutation carriers.

Single Nucleotide Polymorphisms (SNPs) in genes involved in the DNA Base Excision Repair (BER) pathway could be associated with cancer risk in carriers of mutations in the high-penetrance susceptibility genes BRCA1 and BRCA2, given the relation of synthetic lethality that exists between one of the components of the BER pathway, PARP1 (poly ADP ribose polymerase), and both BRCA1 and BRCA2. In the present study, we have performed a comprehensive analysis of 18 genes involved in BER using a tagging SNP approach in a large series of BRCA1 and BRCA2 mutation carriers. 144 SNPs were analyzed in a two stage study involving 23,463 carriers from the CIMBA consortium (the Consortium of Investigators of Modifiers of BRCA1 and BRCA2). Eleven SNPs showed evidence of association with breast and/or ovarian cancer at p<0.05 in the combined analysis. Four of the five genes for which strongest evidence of association was observed were DNA glycosylases. The strongest evidence was for rs1466785 in the NEIL2 (endonuclease VIII-like 2) gene (HR: 1.09, 95% CI (1.03-1.16), p = 2.7 × 10(-3)) for association with breast cancer risk in BRCA2 mutation carriers, and rs2304277 in the OGG1 (8-guanine DNA glycosylase) gene, with ovarian cancer risk in BRCA1 mutation carriers (HR: 1.12 95%CI: 1.03-1.21, p = 4.8 × 10(-3)). DNA glycosylases involved in the first steps of the BER pathway may be associated with cancer risk in BRCA1/2 mutation carriers and should be more comprehensively studied.

Common non-synonymous SNPs associated with breast cancer susceptibility: findings from the Breast Cancer Association Consortium.

Candidate variant association studies have been largely unsuccessful in identifying common breast cancer susceptibility variants, although most studies have been underpowered to detect associations of a realistic magnitude. We assessed 41 common non-synonymous single-nucleotide polymorphisms (nsSNPs) for which evidence of association with breast cancer risk had been previously reported. Case-control data were combined from 38 studies of white European women (46 450 cases and 42 600 controls) and analyzed using unconditional logistic regression. Strong evidence of association was observed for three nsSNPs: ATXN7-K264R at 3p21 [rs1053338, per allele OR = 1.07, 95% confidence interval (CI) = 1.04-1.10, P = 2.9 × 10(-6)], AKAP9-M463I at 7q21 (rs6964587, OR = 1.05, 95% CI = 1.03-1.07, P = 1.7 × 10(-6)) and NEK10-L513S at 3p24 (rs10510592, OR = 1.10, 95% CI = 1.07-1.12, P = 5.1 × 10(-17)). The first two associations reached genome-wide statistical significance in a combined analysis of available data, including independent data from nine genome-wide association studies (GWASs): for ATXN7-K264R, OR = 1.07 (95% CI = 1.05-1.10, P = 1.0 × 10(-8)); for AKAP9-M463I, OR = 1.05 (95% CI = 1.04-1.07, P = 2.0 × 10(-10)). Further analysis of other common variants in these two regions suggested that intronic SNPs nearby are more strongly associated with disease risk. We have thus identified a novel susceptibility locus at 3p21, and confirmed previous suggestive evidence that rs6964587 at 7q21 is associated with risk. The third locus, rs10510592, is located in an established breast cancer susceptibility region; the association was substantially attenuated after adjustment for the known GWAS hit. Thus, each of the associated nsSNPs is likely to be a marker for another, non-coding, variant causally related to breast cancer risk. Further fine-mapping and functional studies are required to identify the underlying risk-modifying variants and the genes through which they act.

A large-scale assessment of two-way SNP interactions in breast cancer susceptibility using 46,450 cases and 42,461 controls from the breast cancer association consortium.

Part of the substantial unexplained familial aggregation of breast cancer may be due to interactions between common variants, but few studies have had adequate statistical power to detect interactions of realistic magnitude. We aimed to assess all two-way interactions in breast cancer susceptibility between 70,917 single nucleotide polymorphisms (SNPs) selected primarily based on prior evidence of a marginal effect. Thirty-eight international studies contributed data for 46,450 breast cancer cases and 42,461 controls of European origin as part of a multi-consortium project (COGS). First, SNPs were preselected based on evidence (P < 0.01) of a per-allele main effect, and all two-way combinations of those were evaluated by a per-allele (1 d.f.) test for interaction using logistic regression. Second, all 2.5 billion possible two-SNP combinations were evaluated using Boolean operation-based screening and testing, and SNP pairs with the strongest evidence of interaction (P < 10(-4)) were selected for more careful assessment by logistic regression. Under the first approach, 3277 SNPs were preselected, but an evaluation of all possible two-SNP combinations (1 d.f.) identified no interactions at P < 10(-8). Results from the second analytic approach were consistent with those from the first (P > 10(-10)). In summary, we observed little evidence of two-way SNP interactions in breast cancer susceptibility, despite the large number of SNPs with potential marginal effects considered and the very large sample size. This finding may have important implications for risk prediction, simplifying the modelling required. Further comprehensive, large-scale genome-wide interaction studies may identify novel interacting loci if the inherent logistic and computational challenges can be overcome.

A candidate gene study of capecitabine-related toxicity in colorectal cancer identifies new toxicity variants at DPYD and a putative role for ENOSF1 rather than TYMS.

OBJECTIVE: Capecitabine is an oral 5-fluorouracil (5-FU) pro-drug commonly used to treat colorectal carcinoma and other tumours. About 35% of patients experience dose-limiting toxicity. The few proven genetic biomarkers of 5-FU toxicity are rare variants and polymorphisms, respectively, at candidate loci dihydropyrimidine dehydrogenase (DPYD) and thymidylate synthase (TYMS). DESIGN: We investigated 1456 polymorphisms and rare coding variants near 25 candidate 5-FU pathway genes in 968 UK patients from the QUASAR2 clinical trial. RESULTS: We identified the first common DPYD polymorphisms to be consistently associated with capecitabine toxicity, rs12132152 (toxicity allele frequency (TAF)=0.031, OR=3.83, p=4.31×10(-6)) and rs12022243 (TAF=0.196, OR=1.69, p=2.55×10(-5)). rs12132152 was particularly strongly associated with hand-foot syndrome (OR=6.1, p=3.6×10(-8)). The rs12132152 and rs12022243 associations were independent of each other and of previously reported DPYD toxicity variants. Next-generation sequencing additionally identified rare DPYD variant p.Ala551Thr in one patient with severe toxicity. Using functional predictions and published data, we assigned p.Ala551Thr as causal for toxicity. We found that polymorphism rs2612091, which lies within an intron of ENOSF1, was also associated with capecitabine toxicity (TAF=0.532, OR=1.59, p=5.28×10(-6)). ENSOF1 is adjacent to TYMS and there is a poorly characterised regulatory interaction between the two genes/proteins. Unexpectedly, rs2612091 fully explained the previously reported associations between capecitabine toxicity and the supposedly functional TYMS variants, 5'VNTR 2R/3R and 3'UTR 6 bp ins-del. rs2612091 genotypes were, moreover, consistently associated with ENOSF1 mRNA levels, but not with TYMS expression. CONCLUSIONS: DPYD harbours rare and common capecitabine toxicity variants. The toxicity polymorphism in the TYMS region may actually act through ENOSF1.