Search

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

Abstract Background Genome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor (ER) status, but how these variants relate with other tumor features and intrinsic molecular subtypes is unclear. Methods Among 106,571 invasive breast cancer cases and 95,762 controls of European ancestry with data on 173 breast cancer variants identified in previous GWAS, we used novel two-stage polytomous logistic regression models to evaluate variants in relation to multiple tumor features (ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and grade) adjusting for each other, and to intrinsic-like subtypes. Results Eighty-five of 173 variants were associated with at least one tumor feature (false discovery rate

Cardiotoxicity due to anthracyclines (CDA) affects cancer patients, but we cannot predict who may suffer from this complication. CDA is a complex trait with a polygenic component that is mainly unidentified. We propose that levels of intermediate molecular phenotypes (IMPs) in the myocardium associated with histopathological damage could explain CDA susceptibility, so variants of genes encoding these IMPs could identify patients susceptible to this complication. Thus, a genetically heterogeneous cohort of mice (n = 165) generated by backcrossing were treated with doxorubicin and docetaxel. We quantified heart fibrosis using an Ariol slide scanner and intramyocardial levels of IMPs using multiplex bead arrays and QPCR. We identified quantitative trait loci linked to IMPs (ipQTLs) and cdaQTLs via linkage analysis. In three cancer patient cohorts, CDA was quantified using echocardiography or Cardiac Magnetic Resonance. CDA behaves as a complex trait in the mouse cohort. IMP levels in the myocardium were associated with CDA. ipQTLs integrated into genetic models with cdaQTLs account for more CDA phenotypic variation than that explained by cda-QTLs alone. Allelic forms of genes encoding IMPs associated with CDA in mice, including AKT1, MAPK14, MAPK8, STAT3, CAS3, and TP53, are genetic determinants of CDA in patients. Two genetic risk scores for pediatric patients (n = 71) and women with breast cancer (n = 420) were generated using machine-learning Least Absolute Shrinkage and Selection Operator (LASSO) regression. Thus, IMPs associated with heart damage identify genetic markers of CDA risk, thereby allowing more personalized patient management.

Background Comprehensive molecular studies on tumours are needed to delineate immortalization process steps and identify sensitive prognostic biomarkers in thyroid cancer. Methods and Results In this study, we extensively characterize telomere‐related alterations in a series of 106 thyroid tumours with heterogeneous clinical outcomes. Using a custom‐designed RNA‐seq panel, we identified five telomerase holoenzyme‐complex genes upregulated in clinically aggressive tumours compared to tumours from long‐term disease‐free patients, being TERT and TERC denoted as independent prognostic markers by multivariate regression model analysis. Characterization of alterations related to TERT re‐expression revealed that promoter mutations, methylation and/or copy gains exclusively co‐occurred in clinically aggressive tumours. Quantitative‐FISH (fluorescence in situ hybridization) analysis of telomere lengths showed a significant shortening in these carcinomas, which matched with a high proliferative rate measured by Ki‐67 immunohistochemistry. RNA‐seq data analysis indicated that short‐telomere tumours exhibit an increased transcriptional activity in the 5‐Mb‐subtelomeric regions, site of several telomerase‐complex genes. Gene upregulation enrichment was significant for specific chromosome‐ends such as the 5p, where TERT is located. Co‐FISH analysis of 5p‐end and TERT loci showed a more relaxed chromatin configuration in short telomere‐length tumours compared to normal telomere‐length tumours. Conclusions Overall, our findings support that telomere shortening leads to a 5p subtelomeric region reorganization, facilitating the transcription and accumulation of alterations at TERT‐locus.

The implementation of pharmacogenetics (PGx) is a main milestones of precision medicine nowadays in order to achieve safer and more effective therapies. Nevertheless, the implementation of PGx diagnostics is extremely slow and unequal worldwide, in part due to a lack of ethnic PGx information. We analysed genetic data from 3006 Spanish individuals obtained by different high-throughput (HT) techniques. Allele frequencies were determined in our population for the main 21 actionable PGx genes associated with therapeutical changes. We found that 98% of the Spanish population harbours at least one allele associated with a therapeutical change and, thus, there would be a need for a therapeutical change in a mean of 3.31 of the 64 associated drugs. We also identified 326 putative deleterious variants that were not previously related with PGx in 18 out of the 21 main PGx genes evaluated and a total of 7122 putative deleterious variants for the 1045 PGx genes described. Additionally, we performed a comparison of the main HT diagnostic techniques, revealing that after whole genome sequencing, genotyping with the PGx HT array is the most suitable solution for PGx diagnostics. Finally, all this information was integrated in the Collaborative Spanish Variant Server to be available to and updated by the scientific community.

Around 50% of the familial breast cancer (BC) cases are estimated to be caused by germline variants in known low-, moderate-, and high-risk susceptibility genes, while the other half is of unknown genetic origin. In the present study, we wanted to evaluate the role of the RECQ helicases, some of which have been studied in the past as candidates, with unclear results about their role in the disease. Using next-generation sequencing (NGS) technology, we analyzed the whole coding sequence of BLM, RECQL1, RECQL4, RECQL5, and WRN in almost 2000 index cases from BC Spanish families that had previously tested negative for the known BC susceptibility genes (BRCAX) and compared the results with the controls extracted from gnomAD. Our results suggest that BLM, RECQL1, RECQL4, and WRN do not play a major role in BC susceptibility. However, in the combined analysis, joining the present results with those previously reported in a series of 1334 BC Spanish patients and controls, we found a statistically significant association between Loss of Function (LoF) variants in RECQL5 and BC risk, with an OR of 2.56 (p = 0.009; 95% CI, 1.18–4.98). Our findings support our previous work and places the RECQL5 gene as a new moderate-risk BC gene.

Abstract Single nucleotide polymorphisms (SNPs) may modulate individual susceptibility to carcinogens. We designed a genome‐wide association study to characterize individuals presenting extreme phenotypes of high and low risk to develop tobacco‐induced non‐small cell lung cancer (NSCLC), and we validated our results. We hypothesized that this strategy would enrich the frequencies of the alleles that contribute to the observed traits. We genotyped 2.37 million SNPs in 95 extreme phenotype individuals, that is: heavy smokers that either developed NSCLC at an early age (extreme cases); or did not present NSCLC at an advanced age (extreme controls), selected from a discovery set (n = 3631). We validated significant SNPs in 133 additional subjects with extreme phenotypes selected from databases including >39,000 individuals. Two SNPs were validated: rs12660420 (pcombined = 5.66 × 10−5; ORcombined = 2.80), mapping to a noncoding transcript exon of PDE10A; and rs6835978 (pcombined = 1.02 × 10−4; ORcombined = 2.57), an intronic variant in ATP10D. We assessed the relevance of both proteins in early‐stage NSCLC. PDE10A and ATP10DmRNA expressions correlated with survival in 821 stage I–II NSCLC patients (p = 0.01 and p

Anthracyclines are among the most used chemotherapeutic agents in breast cancer (BC). However their use is hampered by anthracycline-induced cardiotoxicity (AIC). The currently known clinical and genetic risk factors do not fully explain the observed inter-individual variability and only have a limited ability to predict which patients are more likely to develop this severe toxicity. To identify novel predictive genes, we conducted a two-stage genome-wide association study in epirubicin-treated BC patients. In the discovery phase, we genotyped over 700,000 single nucleotide variants in a cohort of 227 patients. The most interesting finding was rs62134260, located 4kb upstream of POLRMT (OR = 5.76, P = 2.23 × 10−5). We replicated this association in a validation cohort of 123 patients (P = 0.021). This variant regulates the expression of POLRMT, a gene that encodes a mitochondrial DNA-directed RNA polymerase, responsible for mitochondrial gene expression. Individuals harbouring the risk allele had a decreased expression of POLRMT in heart tissue that may cause an impaired capacity to maintain a healthy mitochondrial population in cardiomyocytes under stressful conditions, as is treatment with epirubicin. This finding suggests a novel molecular mechanism involved in the development of AIC and may improve our ability to predict patients who are at risk.

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

Invasive lobular breast cancer (ILC) accounts for 10-15% of all invasive breast carcinomas. It is generally ER positive (ER+) and often associated with lobular carcinoma in situ (LCIS). Genome-wide association studies have identified more than 70 common polymorphisms that predispose to breast cancer, but these studies included predominantly ductal (IDC) carcinomas. To identify novel common polymorphisms that predispose to ILC and LCIS, we pooled data from 6,023 cases (5,622 ILC, 401 pure LCIS) and 34,271 controls from 36 studies genotyped using the iCOGS chip. Six novel SNPs most strongly associated with ILC/LCIS in the pooled analysis were genotyped in a further 516 lobular cases (482 ILC, 36 LCIS) and 1,467 controls. These analyses identified a lobular-specific SNP at 7q34 (rs11977670, OR (95%CI) for ILC = 1.13 (1.09-1.18), P = 6.0 × 10(-10); P-het for ILC vs IDC ER+ tumors = 1.8 × 10(-4)). Of the 75 known breast cancer polymorphisms that were genotyped, 56 were associated with ILC and 15 with LCIS at P

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.

The identification of the two most prevalent susceptibility genes in breast cancer, BRCA1 and BRCA2, was the beginning of a sustained effort to uncover new genes explaining the missing heritability in this disease. Today, additional high, moderate and low penetrance genes have been identified in breast cancer, such as P53, PTEN, STK11, PALB2 or ATM, globally accounting for around 35 percent of the familial cases. In the present study we used massively parallel sequencing to analyze 7 BRCA1/BRCA2 negative families, each having at least 6 affected women with breast cancer (between 6 and 10) diagnosed under the age of 60 across generations. After extensive filtering, Sanger sequencing validation and co-segregation studies, variants were prioritized through either control-population studies, including up to 750 healthy individuals, or case-control assays comprising approximately 5300 samples. As a result, a known moderate susceptibility indel variant (CHEK2 1100delC) and a catalogue of 11 rare variants presenting signs of association with breast cancer were identified. All the affected genes are involved in important cellular mechanisms like DNA repair, cell proliferation and survival or cell cycle regulation. This study highlights the need to investigate the role of rare variants in familial cancer development by means of novel high throughput analysis strategies optimized for genetically heterogeneous scenarios. Even considering the intrinsic limitations of exome resequencing studies, our findings support the hypothesis that the majority of non-BRCA1/BRCA2 breast cancer families might be explained by the action of moderate and/or low penetrance susceptibility alleles.

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

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.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⁻⁵), and three SNPs in ABCB1, rs4148737 (per-allele HR = 3.66, 95%CI = 1.85-6.11, p-value = 6.9×10⁻⁵), rs1128503 and rs10276036 (r² = 1, per-allele HR = 0.24, 95%CI = 0.11-0.47 p-value = 7.9×10⁻⁵). Associations with these SNPs remained statistically significant after correction for multiple testing (all corrected p-values [permutation test] ≤ 0.03).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.

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.

In order to identify genetic factors related to thyroid cancer susceptibility, we adopted a candidate gene approach. We studied tag- and putative functional SNPs in genes involved in thyroid cell differentiation and proliferation, and in genes found to be differentially expressed in thyroid carcinoma. A total of 768 SNPs in 97 genes were genotyped in a Spanish series of 615 cases and 525 controls, the former comprising the largest collection of patients with this pathology from a single population studied to date. SNPs in an LD block spanning the entire FOXE1 gene showed the strongest evidence of association with papillary thyroid carcinoma susceptibility. This association was validated in a second stage of the study that included an independent Italian series of 482 patients and 532 controls. The strongest association results were observed for rs1867277 (OR[per-allele] = 1.49; 95%CI = 1.30-1.70; P = 5.9x10(-9)). Functional assays of rs1867277 (NM_004473.3:c.-283G>A) within the FOXE1 5' UTR suggested that this variant affects FOXE1 transcription. DNA-binding assays demonstrated that, exclusively, the sequence containing the A allele recruited the USF1/USF2 transcription factors, while both alleles formed a complex in which DREAM/CREB/alphaCREM participated. Transfection studies showed an allele-dependent transcriptional regulation of FOXE1. We propose a FOXE1 regulation model dependent on the rs1867277 genotype, indicating that this SNP is a causal variant in thyroid cancer susceptibility. Our results constitute the first functional explanation for an association identified by a GWAS and thereby elucidate a mechanism of thyroid cancer susceptibility. They also attest to the efficacy of candidate gene approaches in the GWAS era.