Search

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 (13)C(5)-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.

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.

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.

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.

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

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.

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.

Mammographic density (MD) is an intermediate phenotype for breast cancer. Previous studies have identified genetic variants associated with MD; however, much of the genetic contribution to MD is unexplained. We conducted a two-stage genome-wide association analysis among the participants in the "Determinants of Density in Mammographies in Spain" study, together with a replication analysis in women from the Australian MD Twins and Sisters Study. Our discovery set covered a total of 3,351 Caucasian women aged 45 to 68 years, recruited from Spanish breast cancer screening centres. MD was blindly assessed by a single reader using Boyd's scale. A two-stage approach was employed, including a feature selection phase exploring 575,374 SNPs in 239 pairs of women with extreme phenotypes and a verification stage for the 183 selected SNPs in the remaining sample (2,873 women). Replication was conducted in 1,786 women aged 40 to 70 years old recruited via the Australian Twin Registry, where MD were measured using Cumulus-3.0, assessing 14 SNPs with a p value

Exceptional longevity (EL) is a rare phenotype that can cluster in families, and co-segregation of genetic variation in these families may point to candidate genes that could contribute to extended lifespan. In this study, for the first time, we have sequenced a total of seven exomes from exceptionally long-lived siblings (probands ≥ 103 years and at least one sibling ≥ 97 years) that come from three separate families. We have focused on rare functional variants (RFVs) which have ≤ 1% minor allele frequency according to databases and that are likely to alter gene product function. Based on this, we have identified one candidate longevity gene carrying RFVs in all three families, APOB. Interestingly, APOB is a component of lipoprotein particles together with APOE, and variants in the genes encoding these two proteins have been previously associated with human longevity. Analysis of nonfamilial EL cases showed a trend, without reaching statistical significance, toward enrichment of APOB RFVs. We have also identified candidate longevity genes shared between two families (5-13) or within individual families (66-156 genes). Some of these genes have been previously linked to longevity in model organisms, such as PPARGC1A, NRG1, RAD52, RAD51, NCOR1, and ADCY5 genes. This work provides an initial catalog of genes that could contribute to exceptional familial longevity.

1540 Background: Individual susceptibility to carcinogens may depend on the genetic background. We characterized the constitutional exome of individuals presenting extreme phenotypes of high sensitivity and resistance to develop tobacco induced LUAD, correlating the results to KRAS status. Methods: From an identification cohort (n=3,631) we selected 100 caucasian heavy smokers that either developed LUAD at early age (cancer cohort, n=50) or did not develop LUAD or other tumors at advanced age (cancer free cohort, n=50). We sequenced their germline DNA with the Agilent Human Exome Capture v5 (21,522 genes, 357,999 exons). Using logistic regression we selected the most significant variants between both cohorts and correlated them with KRAS mutation status of LUAD patients. Results: mean ages for the cancer and cancer free cohorts were 50 (range 34-55) and 78 years (72-90). Mean tobacco consumptions were 44 (range 6-72) and 55 pack-years (20-124). Median coverage was 96% at >10X; median depth was 97X. Table shows the most significant variants. rs7240666 ( ALPK2) achieved top significance (p=8.14x10-5, OR 0.18). rs78898229 ( ANKRD36C) and rs74866537 ( PTPN4) were predominantly represented in patients with KRAS+ tumors, OR: 16 (3.3-78) and 11.9 (3.3-43); and rs12426243 (CCDC41) in KRAS- tumors (OR: 13 (3-53). Conclusions: Our study characterizes for the first time the genotype of individuals presenting extreme phenotypes of high and low risk to develop tobacco-induced LUAD according to KRAS status. Our results warrants further study to assess their value to screen these clinically relevant phenotypes; and to identify mechanisms of high susceptibility and resistance to carcinogens. [Table: see text]

Mammographic density (MD) is regarded as an intermediate phenotype in breast cancer development. This association study investigated the influence of 14 breast cancer susceptibility loci identified through previous genome-wide association studies on MD among the participants in the "Determinants of Density in Mammographies in Spain" (DDM-Spain) study.Our study covered a total of 3348 Caucasian women aged 45-68years, recruited from seven Spanish breast cancer screening centres having DNA available. Mammographic density was blindly assessed by a single reader using a semiquantitative scale. Ordinal logistic models, adjusted for age, body mass index and menopausal status, were used to estimate the association between each genotype and MD.Evidence of association with MD was found for variant rs3803662 (TOX3) (Odds Ratio (OR)=1.13, 95% Confidence Interval (CI)=1.03-1.25), and marginal evidence of association for susceptibility loci rs3817198 (LSP1) (OR=1.09, 95% CI=1.00-1.20) and rs2981582 (FGFR2) (OR=0.92, 95% CI=0.84-1.01). Two other loci were associated with MD solely among pre-menopausal women, namely, rs4973768 (SLC4A7) (OR=0.83, 95% CI=0.70-1.00) and rs4415084 (MEPS30) (OR=1.22, 95% CI=1.00-1.49).Our findings lend some support to the hypothesis which links these susceptibility loci to MD.

Genome-wide association studies have found SNPs at 17q22 to be associated with breast cancer risk. To identify potential causal variants related to breast cancer risk, we performed a high resolution fine-mapping analysis that involved genotyping 517 SNPs using a custom Illumina iSelect array (iCOGS) followed by imputation of genotypes for 3,134 SNPs in more than 89,000 participants of European ancestry from the Breast Cancer Association Consortium (BCAC). We identified 28 highly correlated common variants, in a 53 Kb region spanning two introns of the STXBP4 gene, that are strong candidates for driving breast cancer risk (lead SNP rs2787486 (OR = 0.92; CI 0.90–0.94; P = 8.96 × 10−15)) and are correlated with two previously reported risk-associated variants at this locus, SNPs rs6504950 (OR = 0.94, P = 2.04 × 10−09, r2 = 0.73 with lead SNP) and rs1156287 (OR = 0.93, P = 3.41 × 10−11, r2 = 0.83 with lead SNP). Analyses indicate only one causal SNP in the region and several enhancer elements targeting STXBP4 are located within the 53 kb association signal. Expression studies in breast tumor tissues found SNP rs2787486 to be associated with increased STXBP4 expression, suggesting this may be a target gene of this locus., published version, peerReviewed

Skeletal muscle regeneration by muscle satellite cells is a physiological mechanism activated upon muscle damage and regulated by Notch signaling. In a family with autosomal recessive limbgirdle muscular dystrophy, we identified a missense mutation in POGLUT1 (protein O-glucosyltransferase 1), an enzyme involved in Notch posttranslational modification and function. In vitro and in vivo experiments demonstrated that the mutation reduces Oglucosyltransferase activity on Notch and impairs muscle development. Muscles from patients revealed decreased Notch signaling, dramatic reduction in satellite cell pool and a muscle-specific adystroglycan hypoglycosylation not present in patients’ fibroblasts. Primary myoblasts from patients showed slow proliferation, facilitated differentiation, and a decreased pool of quiescent PAX7+ cells. A robust rescue of the myogenesis was demonstrated by increasing Notch signaling. None of these alterations were found in muscles from secondary dystroglycanopathy patients. These data suggest that a key pathomechanism for this novel form of muscular dystrophy is Notch-dependent loss of satellite cells. Junta de Andalucía PI-0017-2014

High-throughput single nucleotide polymorphism (SNP)-array technologies allow to investigate copy number variants (CNVs) in genome-wide scans and specific calling algorithms have been developed to determine CNV location and copy number. We report the results of a reliability analysis comparing data from 96 pairs of samples processed with CNVpartition, PennCNV, and QuantiSNP for Infinium Illumina Human 1Million probe chip data. We also performed a validity assessment with multiplex ligation-dependent probe amplification (MLPA) as a reference standard. The number of CNVs per individual varied according to the calling algorithm. Higher numbers of CNVs were detected in saliva than in blood DNA samples regardless of the algorithm used. All algorithms presented low agreement with mean Kappa Index (KI) 0.62). Our results indicate that the current calling algorithms should be improved for high performance CNVanalysis in genome-wide scans. Further refinement is required to assess CNVs as risk factors in complex diseases.

Please cite this paper as: Genetic analysis of three important genes in pigmentation and melanoma susceptibility: CDKN2A, MC1R and HERC2/OCA2. Experimental Dermatology 2010; 19: 836–844. Abstract: The CDKN2A gene is regarded as the major familial malignant melanoma (MM) susceptibility gene. Human pigmentation is one of the main modulators of individual risk of developing MM. Therefore, the genes involved in the determination of skin colour and tanning response are potentially implicated in MM predisposition and may be useful predictors of MM risk in the general population. The human melanocortin-1 receptor gene (MC1R) plays a crucial role in pigmentation and also appears to be important in MM. The OCA2 gene has emerged as a new and significant determinant of human iris colour variation. We present a case–control study in Spanish population including 390 consecutive patients with melanoma and 254 control subjects. Sequence analysis of the entire coding region and genotyping of 5 tag-SNPs in the genomic region of MC1R was performed. We identified 27 variants, two reaching statistical significance [R160W (OR: 4.18, 95% CI: 1.24–14.04, P = 0.02) and D294H (OR: 3.10, 95% CI: 1.37–7.01, P = 0.01)] and we detected two novel non-synonymous changes: V92L and T308M. Odds ratio for carrying two functional variants was 4.25 (95% CI: 2.30–7.84, P = 3.63 × 10−6). Haplotypes of the entire MC1R region have been established, and we observed an enrichment of a rare European haplotype similar to African values carrying variants V92M and I155T. In addition, three potentially functional SNPs were selected in p16/CDKN2A and in the promoter region of OCA2/HERC2. Our data for CDKN2A gene did not reach statistically significant results for any of the two studied alleles. We found that the variant allele A > G of OCA2/HERC2 (rs12913832) was associated with pigmentation features: eye, hair and skin colour; P-values = 1.8 × 10−29, 9.2 × 10−16, 1.1 × 10−3, respectively, validating previous results.

The search for novel oncogenes is important because they could be the target of future specific anticancer therapies. In the present paper we report the identification of novel amplified genes in lung cancer by means of global gene expression analysis. To screen for amplicons, we aligned the gene expression data according to the position of transcripts in the human genome and searched for clusters of over-expressed genes. We found several clusters with gene over-expression, suggesting an underlying genomic amplification. FISH and microarray analysis for DNA copy number in two clusters, at chromosomes 11q12 and 13q34, confirmed the presence of amplifications spanning about 0.4 and 1 Mb for 11q12 and 13q34, respectively. Amplification at these regions each occurred at a frequency of 3%. Moreover, quantitative RT-PCR of each individual transcript within the amplicons allowed us to verify the increased in gene expression of several genes. The p120ctn and DP1 proteins, encoded by two candidate oncogenes, CTNND1 and TFDP1, at 11q12 and 13q amplicons, respectively, showed very strong immunostaining in lung tumours with gene amplification. We then focused on the 13q34 amplicon and in the TFDP1 candidate oncogene. To further determine the oncogenic properties of DP1, we searched for lung cancer cell lines carrying TFDP1 amplification. Depletion of TFDP1 expression by small interference RNA in a lung cancer cell line (HCC33) with TFDP1 amplification and protein over-expression reduced cell viability by 50%. In conclusion, we report the identification of two novel amplicons, at 13q34 and 11q12, each occurring at a frequency of 3% of non-small cell lung cancers. TFDP1, which encodes the E2F-associated transcription factor DP1 is a candidate oncogene at 13q34. The data discussed in this publication have been deposited in NCBIs Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo/) and are accessible through GEO Series Accession No. GSE21168.

Despite all the research efforts made during the last few decades, most of the cases of families with breast cancer remain unexplained. Mutations in BRCA1 and BRCA2, and in other breast-cancer-susceptibility genes, account for about 25% of familial breast cancer. Linkage studies have failed to identify other breast-cancer-susceptibility genes. The selection criteria of the families, differences in the population background, or clinical and genetic heterogeneity, among other factors, might determine the power to detect the linkage signal. We have performed a SNP-based linkage scan with a total of 6000 SNP markers across the genome in 41 breast-cancer Spanish families, with an average of four breast-cancer cases per family not associated with BRCA1 or BRCA2 germline mutations. In addition, we have included three BRCA-positive families to test the power in linkage detection from a low-complexity family in which a high-penetrance mutation segregates. We have identified three regions of interest, located on 3q25, 6q24, and 21q22. The two former regions showed a suggestive linkage signal (HLOD scores 3.01 and 2.26, respectively), and the latter region showed a significant linkage signal (HLOD score 3.55). Moreover, we found that a subset of 13 families with bilateral breast cancer presented a HLOD of 3.13 on the 3q25 region. Our results suggest that several variables must be taken into account before performing a linkage study in familial breast cancer because of the high heterogeneity within non-BRCA1/2 families. Phenotypic and geographic homogeneity could be the most important factors.

GWAS have identified a breast cancer susceptibility locus on 2q35. Here we report the fine mapping of this locus using data from 101,943 subjects from 50 case-control studies. We genotype 276 SNPs using the 'iCOGS' genotyping array and impute genotypes for a further 1,284 using 1000 Genomes Project data. All but two, strongly correlated SNPs (rs4442975 G/T and rs6721996 G/A) are excluded as candidate causal variants at odds against >100:1. The best functional candidate, rs4442975, is associated with oestrogen receptor positive (ER+) disease with an odds ratio (OR) in Europeans of 0.85 (95% confidence interval=0.84-0.87; P=1.7 × 10(-43)) per t-allele. This SNP flanks a transcriptional enhancer that physically interacts with the promoter of IGFBP5 (encoding insulin-like growth factor-binding protein 5) and displays allele-specific gene expression, FOXA1 binding and chromatin looping. Evidence suggests that the g-allele confers increased breast cancer susceptibility through relative downregulation of IGFBP5, a gene with known roles in breast cell biology. ispartof: Nature communications vol:5 pages:4999- ispartof: location:England status: published

OBJECTIVE A primary challenge in identifying replicable pharmacogenomic markers from clinical genomewide association study (GWAS) trials in oncology is the difficulty in performing a second large clinical trial with the same drugs and dosage regimen. We sought to overcome this challenge by incorporating GWAS results from cell-based studies using the same chemotherapy as a clinical cohort. METHODS In this study, we test whether the overlap between genetic variants identified in a preclinical study and a clinical study on capecitabine is more than expected by chance. A GWAS of capecitabine-induced cytotoxicity was performed in 164 lymphoblastoid cell lines derived from the CEU HapMap population and compared with a GWAS of hand-foot syndrome (HFS), the most frequent capecitabine-induced adverse drug reaction, in Spanish breast and colorectal cancer patients (n=160) treated with capecitabine. RESULTS We observed an overlap of 16 single nucleotide polymorphisms associated with capecitabine-induced cytotoxicity (P

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.

Background Peripheral neuropathy is the dose limiting toxicity of paclitaxel, a chemotherapeutic drug widely used to treat solid tumours. This toxicity exhibits great inter-individual variability of unknown origin. The present study aimed to identify genetic variants associated with paclitaxel induced neuropathy via a whole genome approach. Methods A genome-wide association study (GWAS) was performed in 144 white European patients uniformly treated with paclitaxel/carboplatin and for whom detailed data on neuropathy was available. Per allele single nucleotide polymorphism (SNP) associations were assessed by Cox regression, modelling the cumulative dose of paclitaxel up to the development of grade 2 sensory neuropathy. Results The strongest evidence of association was observed for the ephrin type A receptor 4 ( EPHA4 ) locus (rs17348202, p=1.0×10–6), and EPHA6 and EPHA5 were among the top 25 and 50 hits (rs301927, p=3.4×10–5 and rs1159057, p=6.8×10–5), respectively. A meta-analysis of EPHA5 -rs7349683, the top marker for paclitaxel induced neuropathy in a previous GWAS (r2=0.79 with rs1159057), gave a hazard ratio (HR) estimate of 1.68 (p=1.4×10−9). Meta-analysis of the second hit of this GWAS, XKR4 -rs4737264, gave a HR of 1.71 (p=3.1×10−8). Imputed SNPs at LIMK2 locus were also strongly associated with this toxicity (HR=2.78, p=2.0×10−7). Conclusions This study provides independent support of EPHA5 -rs7349683 and XKR4 -rs4737264 as the first markers of risk of paclitaxel induced neuropathy. In addition, it suggests that other EPHA genes also involved in axonal guidance and repair following neural injury, as well as LIMK2 locus, may play an important role in the development of this toxicity. The identified SNPs could form the basis for individualised paclitaxel chemotherapy.

Journal article TERT-locus SNPs and leukocyte telomere measures are reportedly associated with risks of multiple cancers. Using the Illumina custom genotyping array iCOGs, we analyzed ~480 SNPs at the TERT locus in breast (n = 103,991), ovarian (n = 39,774) and BRCA1 mutation carrier (n = 11,705) cancer cases and controls. Leukocyte telomere measurements were also available for 53,724 participants. Most associations cluster into three independent peaks. The minor allele at the peak 1 SNP rs2736108 associates with longer telomeres (P = 5.8 × 10!-7), lower risks for estrogen receptor (ER)-negative (P = 1.0 × 10!-8) and BRCA1 mutation carrier (P = 1.1 × 10!-5) breast cancers and altered promoter assay signal. The minor allele at the peak 2 SNP rs7705526 associates with longer telomeres (P = 2.3 × 10!-14), higher risk of low-malignant-potential ovarian cancer (P = 1.3 × 10!-15) and greater promoter activity. The minor alleles at the peak 3 SNPs rs10069690 and rs2242652 increase ER-negative (P = 1.2 × 10!-12) and BRCA1 mutation carrier (P = 1.6 × 10!-14) breast and invasive ovarian (P = 1.3 × 10!-11) cancer risks but not via altered telomere length. The cancer risk alleles of rs2242652 and rs10069690, respectively, increase silencing and generate a truncated TERT splice variant. European Commission Seventh Framework Programme (agreement 223175-HEALTH-F2-2009-223175); Cancer Research UK (C1287/A10118 and C1287/A12014) peer-reviewed

Common genetic variants contribute to the observed variation in breast cancer risk for BRCA2 mutation carriers; those known to date have all been found through population-based genome-wide association studies (GWAS). To comprehensively identify breast cancer risk modifying loci for BRCA2 mutation carriers, we conducted a deep replication of an ongoing GWAS discovery study. Using the ranked P-values of the breast cancer associations with the imputed genotype of 1.4 M SNPs, 19,029 SNPs were selected and designed for inclusion on a custom Illumina array that included a total of 211,155 SNPs as part of a multi-consortial project. DNA samples from 3,881 breast cancer affected and 4,330 unaffected BRCA2 mutation carriers from 47 studies belonging to the Consortium of Investigators of Modifiers of BRCA1/2 were genotyped and available for analysis. We replicated previously reported breast cancer susceptibility alleles in these BRCA2 mutation carriers and for several regions (including FGFR2, MAP3K1, CDKN2A/B, and PTHLH) identified SNPs that have stronger evidence of association than those previously published. We also identified a novel susceptibility allele at 6p24 that was inversely associated with risk in BRCA2 mutation carriers (rs9348512; per allele HR = 0.85, 95% CI 0.80–0.90, P = 3.9×10−8). This SNP was not associated with breast cancer risk either in the general population or in BRCA1 mutation carriers. The locus lies within a region containing TFAP2A, which encodes a transcriptional activation protein that interacts with several tumor suppressor genes. This report identifies the first breast cancer risk locus specific to a BRCA2 mutation background. This comprehensive update of novel and previously reported breast cancer susceptibility loci contributes to the establishment of a panel of SNPs that modify breast cancer risk in BRCA2 mutation carriers. This panel may have clinical utility for women with BRCA2 mutations weighing options for medical prevention of breast cancer., Author Summary Women who carry BRCA2 mutations have an increased risk of breast cancer that varies widely. To identify common genetic variants that modify the breast cancer risk associated with BRCA2 mutations, we have built upon our previous work in which we examined genetic variants across the genome in relation to breast cancer risk among BRCA2 mutation carriers. Using a custom genotyping platform with 211,155 genetic variants known as single nucleotide polymorphisms (SNPs), we genotyped 3,881 women who had breast cancer and 4,330 women without breast cancer, which represents the largest possible, international collection of BRCA2 mutation carriers. We identified that a SNP located at 6p24 in the genome was associated with lower risk of breast cancer. Importantly, this SNP was not associated with breast cancer in BRCA1 mutation carriers or in a general population of women, indicating that the breast cancer association with this SNP might be specific to BRCA2 mutation carriers. Combining this BRCA2-specific SNP with 13 other breast cancer risk SNPs also known to modify risk in BRCA2 mutation carriers, we were able to derive a risk prediction model that could be useful in helping women with BRCA2 mutations weigh their risk-reduction strategy options.

Lorente-Galdós, Belén et al., Association studies are the choice approach in the discovery of the genomic basis of complex traits. To carry out such analysis, researchers frequently need to (1) select optimally informative sets of Single Nucleotide Polymorphisms (SNPs) in candidate regions and (2) annotate the results of associations found by means of genome-wide SNP arrays. These are complex tasks, since many criteria have to be considered, including the SNPs' functional properties, technological information and haplotype frequencies in given populations. SYSNPs implements algorithms that allow for efficient and simultaneous consideration of all the relevant criteria to obtain sets of SNPs that properly cover arbitrarily large lists of genes or genomic regions. Complementarily, SYSNPs allows for comprehensive functional annotation of SNPs linked to any given marker SNP. SYSNPs dramatically reduces the effort needed for SNP selection from days of searching various databases to a few minutes using a simple browser. © 2012 Inderscience Enterprises Ltd.

Background: Standard treatment for osteosarcoma patients consists of a combination of cisplatin, adriamycin, and methotrexate before surgical resection of the primary tumour, followed by postoperative chemotherapy including vincristine and cyclophosphamide. Unfortunately, many patients still relapse or suffer adverse events. We examined whether common germline polymorphisms in chemotherapeutic transporter and metabolic pathway genes of the drugs used in standard osteosarcoma treatment may predict treatment response. Methodology/Principal Findings: In this study we screened 102 osteosarcoma patients for 346 Single Nucleotide Polymorphisms (SNPs) and 2 Copy Number Variants (CNVs) in 24 genes involved in the metabolism or transport of cisplatin, adriamycin, methotrexate, vincristine, and cyclophosphamide. We studied the association of the genotypes with tumour response and overall survival. We found that four SNPs in two ATP-binding cassette genes were significantly associated with overall survival: rs4148416 in ABCC3 (per-allele HR = 8.14, 95%CI = 2.73-20.2, p-value = 5.1×10 -5), and three SNPs in ABCB1, rs4148737 (per-allele HR = 3.66, 95%CI = 1.85-6.11, p-value = 6.9×10 -5), rs1128503 and rs10276036 (r 2 = 1, per-allele HR = 0.24, 95%CI = 0.11-0.47 p-value = 7.9×10 -5). Associations with these SNPs remained statistically significant after correction for multiple testing (all corrected p-values [permutation test] ≤0.03). Conclusions: Our findings suggest that these polymorphisms may affect osteosarcoma treatment efficacy. If these associations are independently validated, these variants could be used as genetic predictors of clinical outcome in the treatment of osteosarcoma, helping in the design of individualized therapy, This work was supported by the AECC (Asociación Española contra el Cáncer), FIS (Fondo de Investigación Sanitaria-Instituto de Salud Carlos III) and the ‘‘Inocente Inocente’’ Foundation

10 p. As the incidence of Malignant Melanoma (MM) reflects an interaction between skin colour and UV exposure, variations in genes implicated in pigmentation and tanning response to UV may be associated with susceptibility to MM. In this study, 363 SNPs in 65 gene regions belonging to the pigmentation pathway have been successfully genotyped using a SNP array. Five hundred and ninety MM cases and 507 controls were analyzed in a discovery phase I. Ten candidate SNPs based on a p-value threshold of 0.01 were identified. Two of them, rs35414 (SLC45A2) and rs2069398 (SILV/CKD2), were statistically significant after conservative Bonferroni correction. The best six SNPs were further tested in an independent Spanish series (624 MM cases and 789 controls). A novel SNP located on the SLC45A2 gene (rs35414) was found to be significantly associated with melanoma in both phase I and phase II (P < 0.0001). None of the other five SNPs were replicated in this second phase of the study. However, three SNPs in TYR, SILV/CDK2 and ADAMTS20 genes (rs17793678, rs2069398 and rs1510521 respectively) had an overall p-value < 0.05 when considering the whole DNA collection (1214 MM cases and 1296 controls). Both the SLC45A2 and the SILV/CDK2 variants behave as protective alleles, while the TYR and ADAMTS20 variants seem to function as risk alleles. Cumulative effects were detected when these four variants were considered together. Furthermore, individuals carrying two or more mutations in MC1R, a well-known low penetrance melanoma-predisposing gene, had a decreased MM risk if concurrently bearing the SLC45A2 protective variant. To our knowledge, this is the largest study on Spanish sporadic MM cases to date. This study was supported by grants from the Ministerio de Educacion y Ciencia (MEC) [SAF2007-65542-C02-01, SAF2009_07072E; BES-2008-009234 to M.I.V. and CP08-0069 “Miquel Servet SNS Contract” to G.R.]; Ministerio de Ciencia e Innovacion [CGL2008-04066/BOS]; Network RD06-0020-1060; the Fundacion Mutua Madrileña and the Fundacion Hospital Provincial de Castellon, Spain. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

International audience; Molecular studies in a patient with Beckwith-Wiedemann syndrome phenotype who developed two different tumors revealed an unexpected observation of almost complete loss of heterozygosity of all chromosomes. We demonstrate, by means of numerous molecular methods that the absence of maternal contribution in somatic cells is due to high-degree (85%) genome-wide paternal uniparental disomy. Our observations indicate the genome-wide UPD results from diploidization, and have important implications for genetic counseling and tumor surveillance for the growing number of UPD-associated imprinting disorders.

Fanconi anemia (FA) family of proteins participates in the DNA repair pathway by homologous recombination, and it is currently formed by 13 genes. Some of these proteins also confer susceptibility to hereditary breast and ovarian cancer (HBOC), since FANCD1 is the BRCA2 breast cancer susceptibility gene, and FANCN/PALB2 and FANCJ/BRIP1 explain 2% of non-BRCA1/2 HBOC families. Thus, there is an important connection between FA and BRCA pathways. In a previous case–control association study analysing FANCA, FANCD2 and FANCL, we reported an association between FANCD2 and sporadic breast cancer (BC) risk (OR = 1.35). In order to know whether variants in other FA genes could also be involved in this association, we have extended our study with the rest of FA genes and some others implicated in the BRCA pathway. We have also analyzed the correlation with survival, nodal metastasis and hormonal receptors (ER− and PR−). A total of 61 SNPs in ten FA genes (FANC-B, -C, -D1, -E, -F, -G, -I, -J, -M, -N) and five FA related genes (ATM, ATR, BRCA1, H2AX and USP1) were studied in a total of 547 consecutive and nonrelated sporadic BC cases and 552 unaffected controls from the Spanish population. Association analyses reported marginal statistically significant results with the minor allele of intronic SNPs in three genes: BRCA1, BRCA2/FANCD1, and ATM. Survival association with SNPs on FANCC and BRCA2/FANCD1 genes were also reported. Sub-group analyses revealed associations between SNPs on FANCI and ATM and nodal metastasis status and between FANCJ/BRIP1 and FANCN/PALB2 and PR− status.

Human pigmentation appears to be one of the main modulators of individual risk of developing malignant melanoma (MM). A large number of genes are known to be involved in rare pigmentary disorders and explain most of the variation in pigmentation phenotypes seen in human populations. This Spanish case-control study included 205 patients with melanoma and 245 control subjects. Thirty-one single nucleotide polymorphisms (SNPs) in genes that had been mainly associated with congenital pigmentation syndromes (ADTB3A, ATRN, CHS1, EDNRB, HPS, KIT, MGRN1, MITF, MLANA, MYO5A, MYO7A, OA1, OCA2, PAX3 and SOX10) were selected. We found that the variant allele of OCA2 R419Q (rs1800407) was associated with increased risk of MM (OR 1.55, 95% CI 1.04-2.31, P = 0.03). This effect on melanoma risk appeared to be stronger among individuals with solar lentigines, or at least 50 nevi. We also describe, for the first time, an association with the variant S1666C (rs2276288) in the MYO7A gene (OR 1.35; 95% CI 1.04-1.76; P = 0.03). Again, this association appeared to be stronger in several phenotypic groups such as individuals with fair skin and those with childhood sunburns. We also found that several variants in the pigmentation genes considered were associated with intermediate phenotypic characteristics. Our findings highlight the potential importance of pigmentation genes in sporadic MM susceptibility.

The presence of familial history in pheochromocytoma/paraganglioma patients, including syndromic antecedents, leads in the majority of cases to a positive genetic testing for mutations in one of the major susceptibility genes described so far. Furthermore, it has been reported that in the absence of familial antecedents, about 11-24% of patients also carry a mutation in one of these related genes. In these cases, other clinical aspects like bilaterality, multiplicity, location of the tumors, or age at onset can help to recognize the underlying genes involved.The objective of the study was to discuss clinical criteria helpful in the genetic diagnosis, placing special emphasis on apparently sporadic cases.Two hundred thirty-seven nonrelated probands were analyzed for the major susceptibility genes: VHL, RET, SDHB, SDHC, and SDHD. Genetic characterization included both point mutation analysis and gross deletions in the SDH genes performed by multiplex PCR.As expected, all syndromic probands were genetically diagnosed with a mutation affecting either RET or VHL. A total of 79.1% (19 of 24) and 18.4% (31 of 168) of patients presenting with either nonsyndromic familial antecedents or apparently sporadic presentation were found to carry a mutation in one of the susceptibility genes. Finally, we found a Spanish founder effect for two mutations: SDHB c.166_170delCCTCA and SDHD c.129GA.Germline mutations are rare in apparently sporadic probands diagnosed after age 40 yr (3.9% in our series) and mainly involve SDHB. Therefore, we recommend prioritizing SDHB genetic testing in patients developing isolated tumors at any age, especially those with extraadrenal location or malignant behavior.

Background Vitamin D serum levels have been found to be related to sun exposure and diet, together with cell differentiation, growth control and consequently, cancer risk. Vitamin D receptor (VDR) genotypes may influence cancer risk; however, no epidemiological studies in sporadic breast cancer (BC) or malignant melanoma (MM) have been performed in a southern European population. In this study, the VDR gene has been evaluated in two epithelial cancers BC and MM. Methods We have conducted an analysis in 549 consecutive and non-related sporadic BC cases and 556 controls, all from the Spanish population, and 283 MM cases and 245 controls. Genotyping analyses were carried out on four putatively functional SNPs within the VDR gene. Results An association with the minor allele A of the non-synonymous SNP rs2228570 (rs10735810, FokI, Met1Thr) was observed for BC, with an estimated odds ratio (OR) of 1.26 (95% CI = 1.02–1.57; p = 0.036). The synonymous variant rs731236 (TaqI) appeared to be associated with protection from BC (OR = 0.80, 95%CI = 0.64–0.99; p = 0.047). No statistically significant associations with MM were observed for any SNP. Nevertheless, sub-group analyses revealed an association between rs2228570 (FokI) and absence of childhood sunburns (OR = 0.65, p = 0.003), between the 3'utr SNP rs739837 (BglI) and fair skin (OR = 1.31, p = 0.048), and between the promoter SNP rs4516035 and the more aggressive tumour location in head-neck and trunk (OR = 1.54, p = 0.020). Conclusion In summary, we observed associations between SNPs in the VDR gene and BC risk, and a comprehensive analysis using clinical and tumour characteristics as outcome variables has revealed potential associations with MM. These associations required confirmation in independent studies.

Familial breast cancer represents up to 5% of all breast cancer cases. Recently, our group has performed a new SNP-based linkage study in 19 non-BRCA1/2 families. We found that a single family was linked to regions in two different chromosomes (11q13 and 14q21), and observed a non-parametric LOD score of 11.5 in both regions. In the present study, we ruled out any possible translocation between the chromosomes. We also used both a panel of STRs and an indirect approach based on HapMap data to narrow down these regions from 28 to 7 Mb in chromosome 11 and from 14.5 to 8.5 Mb in chromosome 14. We performed a mutational screening on candidate genes in 11q13 (NUMA1, FGF3, CCND1, RAD9A, RNF121, FADD and hsa-mir-192), and on FOXA1 in 14q21. Although we have not found any deleterious mutations in the coding region of these genes, data from STR markers confirm 11q13 as a candidate region to contain a breast cancer susceptibility gene.

Aim: We analyzed the genome of individuals presenting extreme phenotypes of sensitivity and resistance to develop tobacco induced NSCLC to identify SNPs associated with these phenotypes. For this purpose, we used one of the most powerful GWAS platforms available. We hypothesized that SNPs may modulate individual susceptibility to carcinogens and that selection of extreme phenotypes would enrich the frequencies of alleles that contribute to the trait, thus increasing the power to identify them. Methods: From an identification cohort (n=3631) we selected caucasian heavy smokers that either developed NSCLC at an early age (cancer-cohort) or that did not present NSCLC at an advanced age (cancer-free cohort). We analyzed their genomic DNA using the array Illumina HumanOmni 2.5 Quad that includes over 2 million powerful markers selected from the 1000 Genomes Project, targeting genetic variation down to 1% minor allele frequency. Statistical significance of SNPs was calculated using logistic regression and Fisher's test. Results: 96 patients (48 per cohort) were selected. Mean age for the cancer and cancer-free cohorts was 49 years (range 38-55) and 76 years (72-84). Mean tobacco consumption was 41 pack-years (range 18-99) and 68 pack-years (40-120). GWAS identified 8 SNPs differentially expressed by logistic regression and 54 SNP by Fisher's test (p Function Gene Oncogenes MSX2 SOX11 Tumor supressors CSMD1 FOXF1 Tobacco induced NSCLC ABCB5 Regulators of transcription DROSHA HDAC9 KIAA0947 Regulators of inflammation PellinoE3 TRIM9 Related with cancer ABHD6 GRIK1 RAB40B SCN1A SLC24A2 SLC25A26 ZFYVE26 Conclusions: We identified candidate SNPs associated with the risk of developing tobacco-induced NSCLC in individuals with extreme phenotypes. Several identified SNPs were located within or near genes that constitute potentially relevant targets for modulation of cancer risk. Validation of the most significant SNPs in an independent set of individuals with similar phenotypes, selected from the EPIC-Spain project (www.epic-spain.com, n=40,000) is ongoing. Disclosure: All authors have declared no conflicts of interest.

11046 Background: SNPs may modulate individual susceptibility to carcinogens. We used GWAS to identify SNPs associated with individuals presenting extreme phenotypes of sensitivity and resistance t...

analyzed by gene array, immunofluorescence, flow cytometry and thymidine incorporation. Results Anti-TNF induced regulatory macrophages displayed increased numbers of autophagosomes as well as an increased expression of autophagy related transcripts including atg5, atg7, atg9 and atg16l2, suggesting induction of autophagy by IFX treatment. Of all donors, 7 were homozygously carrying the CD associated risk allele, 14 were heterozygous and 7 were homozygous for the WT allele. The number of CD14+ regulatory macrophages correlated significantly with the number of WT alleles present in each individual culture, with the largest number of macrophages found in cultures containing 3 or 4 WT alleles (2 WT donors or 1 WT and 1 heterozygous). Similarly, expression of CD206, which is associated with the immunosuppressive function of macrophages, positively correlated with the number of WT alleles present. To confirm the functional consequences of these findings, IFX mediated suppression of T cell proliferation was determined. Again, the level of suppression correlated significantly with the number of WT alleles present in the respective donor combinations. Conclusion Induction of regulatory macrophages by IFX is associated with induction of autophagy. In vitro, the number and function of IFX induced macrophages correlated with the number of WT alleles for the autophagy related gene ATG16L1. This suggests that an intact autophagy pathway is important for effectiveness of anti-TNF therapy.

BACKGROUND: The classical candidate-gene approach has failed to identify novel breast cancer susceptibility genes. Nowadays, massive parallel sequencing technology allows the development of studies unaffordable a few years ago. However, analysis protocols are not yet sufficiently developed to extract all information from the huge amount of data obtained. METHODOLOGY/PRINCIPAL FINDINGS: In this study, we performed high throughput sequencing in two regions located on chromosomes 3 and 6, recently identified by linkage studies by our group as candidate regions for harbouring breast cancer susceptibility genes. In order to enrich for the coding regions of all described genes located in both candidate regions, a hybrid-selection method on tiling microarrays was performed. CONCLUSIONS/SIGNIFICANCE: We developed an analysis pipeline based on SOAP aligner to identify candidate variants with a high real positive confirmation rate (0.89), with which we identified eight variants considered candidates for functional studies. The results suggest that the present strategy might be a valid second step for identifying high penetrance genes.

Background The recent development of new high-throughput technologies for SNP genotyping has opened the possibility of taking a genome-wide linkage approach to the search for new candidate genes involved in heredity diseases. The two major breast cancer susceptibility genes BRCA1 and BRCA2 are involved in 30% of hereditary breast cancer cases, but the discovery of additional breast cancer predisposition genes for the non-BRCA1/2 breast cancer families has so far been unsuccessful. Results In order to evaluate the power improvement provided by using SNP markers in a real situation, we have performed a whole genome screen of 19 non-BRCA1/2 breast cancer families using 4720 genomewide SNPs with Illumina technology (Illumina's Linkage III Panel), with an average distance of 615 Kb/SNP. We identified six regions on chromosomes 2, 3, 4, 7, 11 and 14 as candidates to contain genes involved in breast cancer susceptibility, and additional fine mapping genotyping using microsatellite markers around linkage peaks confirmed five of them, excluding the region on chromosome 3. These results were consistent in analyses that excluded SNPs in high linkage disequilibrium. The results were compared with those obtained previously using a 10 cM microsatellite scan (STR-GWS) and we found lower or not significant linkage signals with STR-GWS data compared to SNP data in all cases. Conclusion Our results show the power increase that SNPs can supply in linkage studies.