Breast Cancer Risk Genes - Association Analysis in More than 113,000 Women.

BACKGROUND: Genetic testing for breast cancer susceptibility is widely used, but for many genes, evidence of an association with breast cancer is weak, underlying risk estimates are imprecise, and reliable subtype-specific risk estimates are lacking. METHODS: We used a panel of 34 putative susceptibility genes to perform sequencing on samples from 60,466 women with breast cancer and 53,461 controls. In separate analyses for protein-truncating variants and rare missense variants in these genes, we estimated odds ratios for breast cancer overall and tumor subtypes. We evaluated missense-variant associations according to domain and classification of pathogenicity. RESULTS: Protein-truncating variants in 5 genes (ATM, BRCA1, BRCA2, CHEK2, and PALB2) were associated with a risk of breast cancer overall with a P value of less than 0.0001. Protein-truncating variants in 4 other genes (BARD1, RAD51C, RAD51D, and TP53) were associated with a risk of breast cancer overall with a P value of less than 0.05 and a Bayesian false-discovery probability of less than 0.05. For protein-truncating variants in 19 of the remaining 25 genes, the upper limit of the 95% confidence interval of the odds ratio for breast cancer overall was less than 2.0. For protein-truncating variants in ATM and CHEK2, odds ratios were higher for estrogen receptor (ER)-positive disease than for ER-negative disease; for protein-truncating variants in BARD1, BRCA1, BRCA2, PALB2, RAD51C, and RAD51D, odds ratios were higher for ER-negative disease than for ER-positive disease. Rare missense variants (in aggregate) in ATM, CHEK2, and TP53 were associated with a risk of breast cancer overall with a P value of less than 0.001. For BRCA1, BRCA2, and TP53, missense variants (in aggregate) that would be classified as pathogenic according to standard criteria were associated with a risk of breast cancer overall, with the risk being similar to that of protein-truncating variants. CONCLUSIONS: The results of this study define the genes that are most clinically useful for inclusion on panels for the prediction of breast cancer risk, as well as provide estimates of the risks associated with protein-truncating variants, to guide genetic counseling. (Funded by European Union Horizon 2020 programs and others.).

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

Reassessment of genotype 1 hepatitis C virus subtype misclassification by LiPA 2.0: implications for direct-acting antiviral treatment.

The accuracy of LiPA 2.0 for hepatitis C virus 1 (HCV-1) subtype classification was analyzed. LiPA 2.0 genotype results from 101 HCV-1-infected patients were compared to genotype findings determined by direct core sequencing. Eleven (11%) samples were misclassified. Given the influence of the HCV-1-subtype in the anti-HCV therapy response, an alternative classification method is warranted.

Health and wellbeing status of the long-lived individuals of the Spanish LONGECYL cross-sectional study.

BACKGROUND: The increase in life expectancy and long-lived individuals is a challenge for public health and provides an opportunity to understand the determinants of longevity. However, few studies have addressed the factors associated with the health status and quality of life in a long-lived individual population. We described the perceived health, clinical status, quality of life, and dependency for activities of daily living in a representative population in Castile and Leon, Spain. METHODS: A sample of 759 long-lived individuals aged 95 years and older was studied by the Health Sentinel Network of Castile and Leon (Spain) through a health examination and a structured questionnaire covering quality of life (EQ-5D-3), lifestyle habits, diet, working life and family health. A blood sample was taken for the study of biological and genetic markers. Chi Square and logistic regression OR with 95% confidence intervals were used to analyze the determinants of the long-lived individuals' health status. The significant level for the bivariate analysis was established at 0.05. RESULTS: Perceived health was good, very good or excellent in 64.2%, while only 46.0% had a quality-of-life index above 0.5 (ranging from 0 to 1) and 44.1% maintained acceptable independence for activities of daily living. Quality-of-life index was higher in the oldest, (OR 7.98 [2,32-27.41]) above 100 years compared to those under 98, and men had better values for independence than women (OR 2.43 [1.40-4.29]). Cardiovascular diseases were the most prevalent (85.5%), but neurological and mental diseases and vision problems had the highest impact on quality of life and independence. CONCLUSION: The long-lived individuals of Castile and Leon have a relatively well-preserved health status, although the perception of health is higher than that describing their quality of life and dependence. The quality of life was higher in the oldest age group and showed differences according to sex, with a better quality of life in men. Public health policies and programs should take in account the differences by sex and age as well as the prevention and control of the main conditions related with poor quality of life or dependence. Future research must include the interaction among genetic, socioeconomic, environmental, and other clinical factors in the quality of life and disability of long-lived individuals.

A Comprehensive Analysis of 21 Actionable Pharmacogenes in the Spanish Population: From Genetic Characterisation to Clinical Impact.

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.

Exome sequencing identifies breast cancer susceptibility genes and defines the contribution of coding variants to breast cancer risk.

Linkage and candidate gene studies have identified several breast cancer susceptibility genes, but the overall contribution of coding variation to breast cancer is unclear. To evaluate the role of rare coding variants more comprehensively, we performed a meta-analysis across three large whole-exome sequencing datasets, containing 26,368 female cases and 217,673 female controls. Burden tests were performed for protein-truncating and rare missense variants in 15,616 and 18,601 genes, respectively. Associations between protein-truncating variants and breast cancer were identified for the following six genes at exome-wide significance (P < 2.5 × 10-6): the five known susceptibility genes ATM, BRCA1, BRCA2, CHEK2 and PALB2, together with MAP3K1. Associations were also observed for LZTR1, ATR and BARD1 with P < 1 × 10-4. Associations between predicted deleterious rare missense or protein-truncating variants and breast cancer were additionally identified for CDKN2A at exome-wide significance. The overall contribution of coding variants in genes beyond the previously known genes is estimated to be small.

Association of X-linked hydrocephalus and Hirschsprung disease: report of a new patient with a mutation in the L1CAM gene.

X-linked hydrocephalus (XLH) has an incidence of 1/30,000 male births and is characterized by intellectual disability, spastic paraplegia, adducted thumbs, and agenesis of corpus callosum, and/or corticospinal tract. The great proportion of cases is ascribed to loss of function mutations of L1CAM gene. Hirschsprung disease (HSCR) is characterized by the absence of ganglion cells in the submucosal and myenteric plexuses along a variable portion of the intestinal tract and has incidence of about 1/5,000. Although with several genes involved in its pathogenesis, the major HSCR gene is the RET proto-oncogene. To date only a few patients have been reported with both phenotypes and mutations in the L1CAM gene. In this report, we describe a new patient with concurrent XLH and HSCR. L1CAM mutational screening showed the presence of the G698R hemizygous mutation, which is a deleterious substitution affecting a key residue necessary for the correct folding of the protein. Moreover, the patient also carried the transcriptional enhancer RET mutation (c.73 + 9277T > C) in heterozygosis. We speculate that both the RET enhancer variant, and the L1CAM mutation may act in combination to produce the enteric phenotype, probably with the participation of other still unidentified molecular events. While it is obvious that additional studies are necessary to further delineate the association between XLH and HSCR in the presence of L1CAM mutations, the documentation of this new patient reinforces the role of this gene acting either in a direct or indirect way in the pathogenesis of Hirschsprung disease.

Genetic, environmental and life-style factors associated with longevity. Protocol and response of the LONGECYL Study.

OBJECTIVE: To describe the objectives, the methodological approach, the response rate of the Genetic, Environmental and Life-style Factors Study in Castilla y León (Spain). METHOD: The Health Sentinel Network studied a sample of long-lived individuals aged 95 or more (LLI). The study included biological samples processed with the Global Screening Array v3.0 that contains a total of 730,059 markers. Written consent was obtained before the examination. CONCLUSIONS: The LLI contacted were 944, and 760 were completed studied. The 87.4% of LLI were born in Castile and Leon and only 1% were non-native of Spain. Severe cognitive impairment was declared in 8.1% of men and 19.2% of women. Genotyping was performed in 739 LLI, the 78.3% of the contacted sample. Family doctors and nurses achieve high participation in population-based studies. DNA samples were taken from 94% of fully studied LLI, and 100% of these samples where successfully genotyped.

A Large Case-Control Study Performed in Spanish Population Suggests That RECQL5 Is the Only RECQ Helicase Involved in Breast Cancer Susceptibility.

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.

Comprehensive analysis of RET common and rare variants in a series of Spanish Hirschsprung patients confirms a synergistic effect of both kinds of events.

BACKGROUND: RET is the major gene associated to Hirschsprung disease (HSCR) with differential contributions of its rare and common, coding and noncoding mutations to the multifactorial nature of this pathology. In the present study, we have performed a comprehensive study of our HSCR series evaluating the involvement of both RET rare variants (RVs) and common variants (CVs) in the context of the disease. METHODS: RET mutational screening was performed by dHPLC and direct sequencing for the identification of RVs. In addition Taqman technology was applied for the genotyping of 3 RET CVs previously associated to HSCR, including a variant lying in an enhancer domain within RET intron 1 (rs2435357). Statistical analyses were performed using the SPSS v.17.0 to analyze the distribution of the variants. RESULTS: Our results confirm the strongest association to HSCR for the "enhancer" variant, and demonstrate a significantly higher impact of it in male versus female patients. Integration of the RET RVs and CVs analysis showed that in 91.66% of cases with both kinds of mutational events, the enhancer allele is in trans with the allele bearing the RET RV. CONCLUSIONS: A gender effect exists on both the transmission and distribution of rare coding and common HSCR causing mutations. In addition, these RET CVs and RVs seem to act in a synergistic way leading to HSCR phenotype.

POLRMT as a Novel Susceptibility Gene for Cardiotoxicity in Epirubicin Treatment of Breast Cancer Patients.

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.

Regulatory CDH4 Genetic Variants Associate With Risk to Develop Capecitabine-Induced Hand-Foot Syndrome.

Capecitabine-induced hand-foot syndrome (CiHFS) is a common dermatological adverse reaction affecting around 30% of patients with capecitabine-treated cancer, and the main cause of dose reductions and chemotherapy delays. To identify novel genetic factors associated with CiHFS in patients with cancer, we carried out an extreme-phenotype genomewide association study in 166 patients with breast and colorectal capecitabine-treated cancer with replication in a second cohort of 85 patients. We discovered and replicated a cluster of four highly correlated single-nucleotide polymorphisms associated with susceptibility to CiHFS at 20q13.33 locus (top hit = rs6129058, hazard ratio = 2.40, 95% confidence interval = 1.78-3.20; P = 1.2 × 10-8 ). Using circular chromosome conformation capture sequencing, we identified a chromatin contact between the locus containing the risk alleles and the promoter of CDH4, located 90 kilobases away. The risk haplotype was associated with decreased levels of CDH4 mRNA and the protein it encodes, R-cadherin (RCAD), which mainly localizes in the granular layer of the epidermis. In human keratinocytes, CDH4 downregulation resulted in reduced expression of involucrin, a protein of the cornified envelope, an essential structure for skin barrier function. Immunohistochemical analyses revealed that skin from patients with severe CiHFS exhibited low levels of RCAD and involucrin before capecitabine treatment. Our results uncover a novel mechanism underlying individual genetic susceptibility to CiHFS with implications for clinically relevant risk prediction.

Novel association of severe neonatal encephalopathy and Hirschsprung disease in a male with a duplication at the Xq28 region.

BACKGROUND: Hirschsprung disease (HSCR) is a neurocristopathy characterized by the absence of parasympathetic intrinsic ganglion cells in the submucosal and myenteric plexuses along a variable portion of the intestinal tract. In approximately 18% of the cases HSCR also presents with multiple congenital anomalies including recognized syndromes. METHODS: A combination of MLPA and microarray data analysis have been undertaken to refine a duplication at the Xq28 region. RESULTS: In this study we present a new clinical association of severe neonatal encephalopathy (Lubs syndrome) and HSCR, in a male patient carrying a duplication at the Xq28 region which encompasses the MECP2 and L1CAM genes. CONCLUSIONS: While the encephalopathy has been traditionally attributed to the MECP2 gene duplication in patients with Lubs syndrome, here we propose that the enteric phenotype in our patient might be due to the dosage variation of the L1CAM protein, together with additional molecular events not identified yet. This would be in agreement with the hypothesis previously forwarded that mutations in L1CAM may be involved in HSCR development in association with a predisposing genetic background.

Novel MLPA procedure using self-designed probes allows comprehensive analysis for CNVs of the genes involved in Hirschsprung disease.

BACKGROUND: Hirschsprung disease is characterized by the absence of intramural ganglion cells in the enteric plexuses, due to a fail during enteric nervous system formation. Hirschsprung has a complex genetic aetiology and mutations in several genes have been related to the disease. There is a clear predominance of missense/nonsense mutations in these genes whereas copy number variations (CNVs) have been seldom described, probably due to the limitations of conventional techniques usually employed for mutational analysis. In this study, we have looked for CNVs in some of the genes related to Hirschsprung (EDNRB, GFRA1, NRTN and PHOX2B) using the Multiple Ligation-dependent Probe Amplification (MLPA) approach. METHODS: CNVs screening was performed in 208 HSCR patients using a self-designed set of MLPA probes, covering the coding region of those genes. RESULTS: A deletion comprising the first 4 exons in GFRA1 gene was detected in 2 sporadic HSCR patients and in silico approaches have shown that the critical translation initiation signal in the mutant gene was abolished. In this study, we have been able to validate the reliability of this technique for CNVs screening in HSCR. CONCLUSIONS: The implemented MLPA based technique presented here allows CNV analysis of genes involved in HSCR that have not been not previously evaluated. Our results indicate that CNVs could be implicated in the pathogenesis of HSCR, although they seem to be an uncommon molecular cause of HSCR.

Association Between ABCB1 Genetic Variants and Persistent Chemotherapy-Induced Alopecia in Women With Breast Cancer.

Importance: Persistent chemotherapy-induced alopecia (pCIA) has been recently described in patients with breast cancer and in its most severe form occurs in up to 10% of these patients. Genetic risk factors associated with pCIA have not been adequately explored. Objective: To identify genetic variants associated with pCIA. Design, Setting, and Participants: In this genetic association study, 215 women with breast cancer treated with docetaxel-based chemotherapy with a follow-up of 1.5 to 10 years after the end of the treatment were recruited retrospectively through 3 hospital oncology units across Spain between 2005 and 2018. Severe pCIA was defined as lack of scalp hair recovery (Common Terminology Criteria for Adverse Events, version 3.0, grade 2) 18 months or more after the end of treatment. Patients with grade 2 pCIA were selected as cases, and those with no sign of residual alopecia 12 months after the end of docetaxel treatment were selected as controls. A genome-wide association study in a discovery phase was conducted, and logistic regression was used to identify variants associated with the risk to develop this adverse effect. The validity of the association was addressed through a replication phase. Exposures: Docetaxel-based chemotherapy. Main Outcomes and Measures: Genotypes of single-nucleotide variants associated with pCIA. Results: In total, 215 women with breast cancer (median age, 51.6 years; interquartile range, 44-60 years) were recruited (173 patients for the discovery phase and 42 patients for the replication phase). In the discovery phase, ABCB1 genetic variants were associated with risk to develop pCIA. In particular, single-nucleotide variation rs1202179, a regulatory variant located in an enhancer element that interacts with the ABCB1 promoter, was associated with the occurrence of pCIA. This finding was validated in the replication cohort (combined odds ratio, 4.05; 95% CI, 2.46-6.67; P = 3.946 × 10-8). This variant is associated with ABCB1 mRNA expression, and the risk allele was associated with decreased ABCB1 expression levels (P = 1.64 × 10-20). Conclusions and Relevance: This is the first study, to our knowledge, that identifies an association between a regulatory variant in the ABCB1 gene and the occurrence of pCIA in patients with breast cancer who were treated with docetaxel-based therapies. This finding suggests an important insight into the biological mechanisms underlying pCIA and opens the opportunity to explore personalized treatment of these patients.

A novel study of copy number variations in Hirschsprung disease using the multiple ligation-dependent probe amplification (MLPA) technique.

BACKGROUND: Hirschsprung disease (HSCR) is a congenital malformation of the hindgut produced by a disruption in neural crest cell migration during embryonic development. HSCR has a complex genetic etiology and mutations in several genes, mainly the RET proto-oncogene, have been related to the disease. There is a clear predominance of missense/nonsense mutations in these genes whereas copy number variations (CNVs) have been seldom described, probably due to the limitations of conventional techniques usually employed for mutational analysis. METHODS: In this study we have aimed to analyze the presence of CNVs in some HSCR genes (RET, EDN3, GDNF and ZFHX1B) using the Multiple Ligation-dependent Probe Amplification (MLPA) approach. RESULTS: Two alterations in the MLPA profiles of RET and EDN3 were detected, but a detailed inspection showed that the decrease in the corresponding dosages were due to point mutations affecting the hybridization probes regions. CONCLUSION: Our results indicate that CNVs of the gene coding regions analyzed here are not a common molecular cause of Hirschsprung disease. However, further studies are required to determine the presence of CNVs affecting non-coding regulatory regions, as well as other candidate genes.

BRAF V600E Detection in Liquid Biopsies from Pediatric Central Nervous System Tumors.

Pediatric Central Nervous System (CNS) tumors are the most fatal cancer diseases in childhood. Due to their localization and infiltrative nature, some tumor resections or biopsies are not feasible. In those cases, the use of minimally invasive methods as diagnostic, molecular marker detection, prognostic or monitoring therapies are emerging. The analysis of liquid biopsies which contain genetic information from the tumor has been much more widely explored in adults than in children. We compare the detection of BRAF V600E targetable mutation by digital-PCR from cell-free-DNA and EV-derived DNA (ctDNA) in serum, plasma and cerebrospinal fluid (CSF) isolated from a cohort of 29 CNS pediatric patients. Here we demonstrate that ctDNA isolated from serum and plasma could be successfully analyzed to obtain tumor genetic information which could be used to guide critical treatment decisions.

Corrigendum: miARma-Seq: a comprehensive tool for miRNA, mRNA and circRNA analysis.

This corrects the article DOI: 10.1038/srep25749.