Ability of a polygenic risk score to refine colorectal cancer risk in Lynch syndrome.

BACKGROUND: Polygenic risk scores (PRSs) have been used to stratify colorectal cancer (CRC) risk in the general population, whereas its role in Lynch syndrome (LS), the most common type of hereditary CRC, is still conflicting. We aimed to assess the ability of PRS to refine CRC risk prediction in European-descendant individuals with LS. METHODS: 1465 individuals with LS (557 MLH1, 517 MSH2/EPCAM, 299 MSH6 and 92 PMS2) and 5656 CRC-free population-based controls from two independent cohorts were included. A 91-SNP PRS was applied. A Cox proportional hazard regression model with 'family' as a random effect and a logistic regression analysis, followed by a meta-analysis combining both cohorts were conducted. RESULTS: Overall, we did not observe a statistically significant association between PRS and CRC risk in the entire cohort. Nevertheless, PRS was significantly associated with a slightly increased risk of CRC or advanced adenoma (AA), in those with CRC diagnosed <50 years and in individuals with multiple CRCs or AAs diagnosed <60 years. CONCLUSION: The PRS may slightly influence CRC risk in individuals with LS in particular in more extreme phenotypes such as early-onset disease. However, the study design and recruitment strategy strongly influence the results of PRS studies. A separate analysis by genes and its combination with other genetic and non-genetic risk factors will help refine its role as a risk modifier in LS.

Altered regulation of BRCA1 exon 11 splicing is associated with breast cancer risk in carriers of BRCA1 pathogenic variants.

Germline pathogenic variants in BRCA1 confer a high risk of developing breast and ovarian cancer. The BRCA1 exon 11 (formally exon 10) is one of the largest exons and codes for the nuclear localization signals of the corresponding gene product. This exon can be partially or entirely skipped during pre-mRNA splicing, leading to three major in-frame isoforms that are detectable in most cell types and tissue, and in normal and cancer settings. However, it is unclear whether the splicing imbalance of this exon is associated with cancer risk. Here we identify a common genetic variant in intron 10, rs5820483 (NC_000017.11:g.43095106_43095108dup), which is associated with exon 11 isoform expression and alternative splicing, and with the risk of breast cancer, but not ovarian cancer, in BRCA1 pathogenic variant carriers. The identification of this genetic effect was confirmed by analogous observations in mouse cells and tissue in which a loxP sequence was inserted in the syntenic intronic region. The prediction that the rs5820483 minor allele variant would create a binding site for the splicing silencer hnRNP A1 was confirmed by pull-down assays. Our data suggest that perturbation of BRCA1 exon 11 splicing modifies the breast cancer risk conferred by pathogenic variants of this gene.

High-sensitivity microsatellite instability assessment for the detection of mismatch repair defects in normal tissue of biallelic germline mismatch repair mutation carriers.

INTRODUCTION: Lynch syndrome (LS) and constitutional mismatch repair deficiency (CMMRD) are hereditary cancer syndromes associated with mismatch repair (MMR) deficiency. Tumours show microsatellite instability (MSI), also reported at low levels in non-neoplastic tissues. Our aim was to evaluate the performance of high-sensitivity MSI (hs-MSI) assessment for the identification of LS and CMMRD in non-neoplastic tissues. MATERIALS AND METHODS: Blood DNA samples from 131 individuals were grouped into three cohorts: baseline (22 controls), training (11 CMMRD, 48 LS and 15 controls) and validation (18 CMMRD and 18 controls). Custom next generation sequencing panel and bioinformatics pipeline were used to detect insertions and deletions in microsatellite markers. An hs-MSI score was calculated representing the percentage of unstable markers. RESULTS: The hs-MSI score was significantly higher in CMMRD blood samples when compared with controls in the training cohort (p<0.001). This finding was confirmed in the validation set, reaching 100% specificity and sensitivity. Higher hs-MSI scores were detected in biallelic MSH2 carriers (n=5) compared with MSH6 carriers (n=15). The hs-MSI analysis did not detect a difference between LS and control blood samples (p=0.564). CONCLUSIONS: The hs-MSI approach is a valuable tool for CMMRD diagnosis, especially in suspected patients harbouring MMR variants of unknown significance or non-detected biallelic germline mutations.

Integrating gene expression and epidemiological data for the discovery of genetic interactions associated with cancer risk.

Dozens of common genetic variants associated with cancer risk have been identified through genome-wide association studies (GWASs). However, these variants only explain a modest fraction of the heritability of disease. The missing heritability has been attributed to several factors, among them the existence of genetic interactions (G × G). Systematic screens for G × G in model organisms have revealed their fundamental influence in complex phenotypes. In this scenario, G × G overlap significantly with other types of gene and/or protein relationships. Here, by integrating predicted G × G from GWAS data and complex- and context-defined gene coexpression profiles, we provide evidence for G × G associated with cancer risk. G × G predicted from a breast cancer GWAS dataset identified significant overlaps [relative enrichments (REs) of 8-36%, empirical P values < 0.05 to 10(-4)] with complex (non-linear) gene coexpression in breast tumors. The use of gene or protein data not specific for breast cancer did not reveal overlaps. According to the predicted G × G, experimental assays demonstrated functional interplay between lipoma-preferred partner and transforming growth factor-ß signaling in the MCF10A non-tumorigenic mammary epithelial cell model. Next, integration of pancreatic tumor gene expression profiles with pancreatic cancer G × G predicted from a GWAS corroborated the observations made for breast cancer risk (REs of 25-59%). The method presented here can potentially support the identification of genetic interactions associated with cancer risk, providing novel mechanistic hypotheses for carcinogenesis.

VAV3 mediates resistance to breast cancer endocrine therapy.

INTRODUCTION: Endocrine therapies targeting cell proliferation and survival mediated by estrogen receptor α (ERα) are among the most effective systemic treatments for ERα-positive breast cancer. However, most tumors initially responsive to these therapies acquire resistance through mechanisms that involve ERα transcriptional regulatory plasticity. Herein we identify VAV3 as a critical component in this process. METHODS: A cell-based chemical compound screen was carried out to identify therapeutic strategies against resistance to endocrine therapy. Binding to ERα was evaluated by molecular docking analyses, an agonist fluoligand assay and short hairpin (sh)RNA-mediated protein depletion. Microarray analyses were performed to identify altered gene expression. Western blot analysis of signaling and proliferation markers, and shRNA-mediated protein depletion in viability and clonogenic assays, were performed to delineate the role of VAV3. Genetic variation in VAV3 was assessed for association with the response to tamoxifen. Immunohistochemical analyses of VAV3 were carried out to determine its association with therapeutic response and different tumor markers. An analysis of gene expression association with drug sensitivity was carried out to identify a potential therapeutic approach based on differential VAV3 expression. RESULTS: The compound YC-1 was found to comparatively reduce the viability of cell models of acquired resistance. This effect was probably not due to activation of its canonical target (soluble guanylyl cyclase), but instead was likely a result of binding to ERα. VAV3 was selectively reduced upon exposure to YC-1 or ERα depletion, and, accordingly, VAV3 depletion comparatively reduced the viability of cell models of acquired resistance. In the clinical scenario, germline variation in VAV3 was associated with the response to tamoxifen in Japanese breast cancer patients (rs10494071 combined P value = 8.4 × 10-4). The allele association combined with gene expression analyses indicated that low VAV3 expression predicts better clinical outcome. Conversely, high nuclear VAV3 expression in tumor cells was associated with poorer endocrine therapy response. Based on VAV3 expression levels and the response to erlotinib in cancer cell lines, targeting EGFR signaling may be a promising therapeutic strategy. CONCLUSIONS: This study proposes VAV3 as a biomarker and a rationale for its use as a signaling target to prevent and/or overcome resistance to endocrine therapy in breast cancer.

Interplay between BRCA1 and RHAMM regulates epithelial apicobasal polarization and may influence risk of breast cancer.

Differentiated mammary epithelium shows apicobasal polarity, and loss of tissue organization is an early hallmark of breast carcinogenesis. In BRCA1 mutation carriers, accumulation of stem and progenitor cells in normal breast tissue and increased risk of developing tumors of basal-like type suggest that BRCA1 regulates stem/progenitor cell proliferation and differentiation. However, the function of BRCA1 in this process and its link to carcinogenesis remain unknown. Here we depict a molecular mechanism involving BRCA1 and RHAMM that regulates apicobasal polarity and, when perturbed, may increase risk of breast cancer. Starting from complementary genetic analyses across families and populations, we identified common genetic variation at the low-penetrance susceptibility HMMR locus (encoding for RHAMM) that modifies breast cancer risk among BRCA1, but probably not BRCA2, mutation carriers: n = 7,584, weighted hazard ratio ((w)HR) = 1.09 (95% CI 1.02-1.16), p(trend) = 0.017; and n = 3,965, (w)HR = 1.04 (95% CI 0.94-1.16), p(trend) = 0.43; respectively. Subsequently, studies of MCF10A apicobasal polarization revealed a central role for BRCA1 and RHAMM, together with AURKA and TPX2, in essential reorganization of microtubules. Mechanistically, reorganization is facilitated by BRCA1 and impaired by AURKA, which is regulated by negative feedback involving RHAMM and TPX2. Taken together, our data provide fundamental insight into apicobasal polarization through BRCA1 function, which may explain the expanded cell subsets and characteristic tumor type accompanying BRCA1 mutation, while also linking this process to sporadic breast cancer through perturbation of HMMR/RHAMM.

Serum 25-hydroxyvitamin D3 levels and vitamin D receptor variants in melanoma patients from the Mediterranean area of Barcelona.

BACKGROUND: Serum 25-hydroxyvitamin D3 (Vitamin D) insufficiency and single-nucleotide polymorphisms (SNPs) on its receptor, Vitamin D receptor (VDR), have been reported to be involved in melanoma susceptibility in populations mostly from northern countries. OBJECTIVE: To investigate 25-hydroxyvitamin D3 levels and VDR SNPs in melanoma patients from sunny area of Barcelona, two studies were carried out. The first study evaluated the levels of Vitamin D at time of melanoma diagnosis and the second one analyzed the association between VDR genetic variants and risk of having a high nevus number, the strongest phenotypic risk factor for melanoma. METHODS: The levels of 25-hydroxyvitamin D3 in 81 melanoma patients at diagnosis were measured. In a second group of melanoma patients, including 150 with low and 113 with high nevus number, 11 VDR SNPs were analyzed for their association with nevus number. RESULTS: In the first study, 68% of patients had insufficient levels of 25-hydroxyvitamin D3 (<25 ng/ml). Autumn-winter months and fair phototype were associated with 25-hydroxyvitamin D3 insufficiency; after multivariate analysis, season of sampling remained the only independent predictor of 25-hydroxyvitamin D3 levels. In the second study, VDR variant rs2189480 (P = 0.006) was associated with risk of high nevus number whereas rs2239179 (P = 0.044) and rs7975128 (P = 0.0005) were protective against high nevus number. After Bonferroni adjustment only rs7975128 remained significant. In stratified analysis, SNP rs7975128 was found protective against multiple melanomas (P = 0.021). CONCLUSION: This study showed that even in Barcelona, a sunny Mediterranean area, 25-hydroxyvitamin D3 levels were sub-optimal in the majority of melanoma patients at diagnosis. The involvement of VDR in nevi and, in turn, in melanoma susceptibility has also been suggested. Larger studies are needed to confirm our findings.

Oxidative stress-induced FAK activation contributes to uterine serous carcinoma aggressiveness.

Uterine serous carcinoma (USC) is an aggressive form of endometrial cancer (EC), characterized by its high propensity for metastases. In fact, while endometrioid endometrial carcinoma (EEC), which accounts for 85% of EC, presents a good prognosis, USC is the most frequently fatal. Herein, we used for the first time a peptide-based tyrosine-kinase-activity profiling approach to quantify the changes in tyrosine kinase activation between USC and EEC. Among the tyrosine kinases highly activated in USC, we identified focal adhesion kinase (FAK). We conducted mechanistic studies using cellular models. In a USC cell line, targeting FAK either by inhibitors PF-573228 and defactinib (VS-6063) or by gene silencing limits 3D cell growth and reduces cell migration. Moreover, results from our studies suggest that oxidative stress is increased in USC tumors compared to EEC ones. Reactive oxygen species (ROS) induce tyrosine phosphorylation of FAK and a concomitant tyrosine phosphorylation of paxillin, a mediator of FAK signal transduction. Mechanistically, by tracking hundreds of individual cells per condition, we show that ROS increased cell distance and migration velocity, highlighting the role of ROS-FAK-PAX signaling in cell migration. Both defactinib and ROS scavenger N-acetylcysteine (NAC) revert this effect, pointing toward ROS as potential culprits for the increase in USC cell motility. A proof of concept of the role of FAK in controlling cell growth was obtained in in vivo experiments using cancer-tissue-originated spheroids (CTOS) and a patient-derived orthotopic xenograft model (orthoxenograft/PDOX). Defactinib reduces cell proliferation and protein oxidation, supporting a pro-tumoral antioxidant role of FAK, whereas antioxidant NAC reverts FAK inhibitor effects. Overall, our data points to ROS-mediated FAK activation in USC as being responsible for the poor prognosis of this tumor type and emphasize the potential of FAK inhibition for USC treatment.

MTHFR C677T and A1298C polymorphism's effect on risk of colorectal cancer in Lynch syndrome.

Lynch syndrome (LS) is characterised by an increased risk of developing colorectal cancer (CRC) and other extracolonic epithelial cancers. It is caused by pathogenic germline variants in DNA mismatch repair (MMR) genes or the EPCAM gene, leading to a less functional DNA MMR system. Individuals diagnosed with LS (LS individuals) have a 10-80% lifetime risk of developing cancer. However, there is considerable variability in the age of cancer onset, which cannot be attributed to the specific MMR gene or variant alone. It is speculated that multiple genetic and environmental factors contribute to this variability, including two single nucleotide polymorphisms (SNPs) in the methylenetetrahydrofolate reductase (MTHFR) gene: C677T (rs1801133) and A1298C (rs1801131). By decreasing MTHFR activity, these SNPs theoretically reduce the silencing of DNA repair genes and increase the availability of nucleotides for DNA synthesis and repair, thereby protecting against early-onset cancer in LS. We investigated the effect of these SNPs on LS disease expression in 2,723 LS individuals from Australia, Poland, Germany, Norway and Spain. The association between age at cancer onset and SNP genotype (risk of cancer) was estimated using Cox regression adjusted for gender, country and affected MMR gene. For A1298C (rs1801131), both the AC and CC genotypes were significantly associated with a reduced risk of developing CRC compared to the AA genotype, but no association was seen for C677T (rs1801133). However, an aggregated effect of protective alleles was seen when combining the alleles from the two SNPs, especially for LS individuals carrying 1 and 2 alleles. For individuals with germline pathogenic variants in MLH1, the CC genotype of A1298C was estimated to reduce the risk of CRC significantly by 39% (HR = 0.61, 95% CI 0.42, 0.89, p = 0.011), while for individuals with pathogenic germline MSH2 variants, the AC genotype (compared to AA) was estimated to reduce the risk of CRC by 26% (HR = 0.66, 95% CI 0.53, 0.83, p = 0.01). In comparison, no association was observed for C677T (rs1801133). In conclusion, our study suggests that combining the MMR gene information with the MTHFR genotype, including the aggregated effect of protective alleles, could be useful in developing an algorithm that estimates the risk of CRC in LS individuals.

Metabolomic Analysis Points to Bioactive Lipid Species and Acireductone Dioxygenase 1 (ADI1) as Potential Therapeutic Targets in Poor Prognosis Endometrial Cancer.

Metabolomic profiling analysis has the potential to highlight new molecules and cellular pathways that may serve as potential therapeutic targets for disease treatment. In this study, we used an LC-MS/MS platform to define, for the first time, the specific metabolomic signature of uterine serous carcinoma (SC), a relatively rare and aggressive variant of endometrial cancer (EC) responsible for 40% of all endometrial cancer-related deaths. A metabolomic analysis of 31 ECs (20 endometrial endometrioid carcinomas (EECs) and 11 SCs) was performed. Following multivariate statistical analysis, we identified 232 statistically different metabolites among the SC and EEC patient samples. Notably, most of the metabolites identified (89.2%) were lipid species and showed lower levels in SCs when compared to EECs. In addition to lipids, we also documented metabolites belonging to amino acids and purine nucleotides (such as 2-Oxo-4-methylthiobutanoic acid, synthesised by acireductone dioxygenase 1 (ADI1) enzyme), which showed higher levels in SCs. To further investigate the role of ADI1 in SC, we analysed the expression protein levels of ADI1 in 96 ECs (67 EECs and 29 SCs), proving that the levels of ADI1 were higher in SCs compared to EECs. We also found that ADI1 mRNA levels were higher in p53 abnormal ECs compared to p53 wild type tumours. Furthermore, elevated ADI1 mRNA levels showed a statistically significant negative correlation with overall survival and progression-free survival among EEC patients. Finally, we tested the ability of ADI1 to induce migration and invasion capabilities in EC cell lines. Altogether, these results suggest that ADI1 could be a potential therapeutic target in poor-prognosis SCs and other Ecs with abnormal p53 expression.

A compound directed against S6K1 hampers fat mass expansion and mitigates diet-induced hepatosteatosis.

The ribosomal protein S6 kinase 1 (S6K1) is a relevant effector downstream of the mammalian target of rapamycin complex 1 (mTORC1), best known for its role in the control of lipid homeostasis. Consistent with this, mice lacking the S6k1 gene have a defect in their ability to induce the commitment of fat precursor cells to the adipogenic lineage, which contributes to a significant reduction of fat mass. Here, we assess the therapeutic blockage of S6K1 in diet-induced obese mice challenged with LY2584702 tosylate, a specific oral S6K1 inhibitor initially developed for the treatment of solid tumors. We show that diminished S6K1 activity hampers fat mass expansion and ameliorates dyslipidemia and hepatic steatosis, while modifying transcriptome-wide gene expression programs relevant for adipose and liver function. Accordingly, decreased mTORC1 signaling in fat (but increased in the liver) segregated with defective epithelial-mesenchymal transition and the impaired expression of Cd36 (coding for a fatty acid translocase) and Lgals1 (Galectin 1) in both tissues. All these factors combined align with reduced adipocyte size and improved lipidomic signatures in the liver, while hepatic steatosis and hypertriglyceridemia were improved in treatments lasting either 3 months or 6 weeks.

ARID1A-deficient cells require HDAC6 for progression of endometrial carcinoma.

AT-rich interactive domain-containing protein 1A (ARID1A) loss-of-function mutation accompanied by a loss of ARID1A protein expression is frequently observed in endometrial carcinomas. However, the molecular mechanisms linking these genetic changes to the altered pathways regulating tumour initiation, maintenance and/or progression remain poorly understood. Thus, the main aim of this study was to analyse the role of ARID1A loss of function in endometrial tumorigenesis. Here, using different endometrial in vitro and in vivo models, such as tumoral cell lines, 3D primary cultures and metastatic or genetically modified mouse models, we show that altered expression of ARID1A is not enough to initiate endometrial tumorigenesis. However, in an established endometrial cancer context, ARID1A loss of function accelerates tumoral progression and metastasis through the disruption of the G2/M cell cycle checkpoint and ATM/ATR-mediated DNA damage checkpoints, increases epithelial cell proliferation rates and induces epithelial mesenchymal transition through the activation of histone deacetylase 6 (HDAC6). Next, we demonstrated that the inhibition of HDAC6 function, using the HDAC6-specific inhibitor ACY1215 or by transfection with HDAC6 short hairpin RNA (shRNA), can reverse the migratory and invasive phenotype of ARID1A-knockdown cells. Further, we also show that inhibition of HDAC6 activity causes an apoptotic vulnerability to etoposide treatments in ARID1A-deficient cells. In summary, the findings exposed in this work indicate that the inhibition of HDAC6 activity is a potential therapeutic strategy for patients suffering from ARID1A-mutant endometrial cancer diagnosed in advanced stages.

Generation and Integrated Analysis of Advanced Patient-Derived Orthoxenograft Models (PDOX) for the Rational Assessment of Targeted Therapies in Endometrial Cancer.

Clinical management of endometrial cancer (EC) is handicapped by the limited availability of second line treatments and bona fide molecular biomarkers to predict recurrence. These limitations have hampered the treatment of these patients, whose survival rates have not improved over the last four decades. The advent of coordinated studies such as The Cancer Genome Atlas Uterine Corpus Endometrial Carcinoma (TCGA_UCEC) has partially solved this issue, but the lack of proper experimental systems still represents a bottleneck that precludes translational studies from successful clinical testing in EC patients. Within this context, the first study reporting the generation of a collection of endometrioid-EC-patient-derived orthoxenograft (PDOX) mouse models is presented that is believed to overcome these experimental constraints and pave the way toward state-of-the-art precision medicine in EC. The collection of primary tumors and derived PDOXs is characterized through an integrative approach based on transcriptomics, mutational profiles, and morphological analysis; and it is demonstrated that EC tumors engrafted in the mouse uterus retain the main molecular and morphological features from analogous tumor donors. Finally, the molecular properties of these tumors are harnessed to assess the therapeutic potential of trastuzumab, a human epidermal growth factor receptor 2 (HER2) inhibitor with growing interest in EC, using patient-derived organotypic multicellular tumor spheroids and in vivo experiments.

A microRNA Cluster Controls Fat Cell Differentiation and Adipose Tissue Expansion By Regulating SNCG.

The H19X-encoded miR-424(322)/503 cluster regulates multiple cellular functions. Here, it is reported for the first time that it is also a critical linchpin of fat mass expansion. Deletion of this miRNA cluster in mice results in obesity, while increasing the pool of early adipocyte progenitors and hypertrophied adipocytes. Complementary loss and gain of function experiments and RNA sequencing demonstrate that miR-424(322)/503 regulates a conserved genetic program involved in the differentiation and commitment of white adipocytes. Mechanistically, it is demonstrated that miR-424(322)/503 targets γ-Synuclein (SNCG), a factor that mediates this program rearrangement by controlling metabolic functions in fat cells, allowing adipocyte differentiation and adipose tissue enlargement. Accordingly, diminished miR-424(322) in mice and obese humans co-segregate with increased SNCG in fat and peripheral blood as mutually exclusive features of obesity, being normalized upon weight loss. The data unveil a previously unknown regulatory mechanism of fat mass expansion tightly controlled by the miR-424(322)/503 through SNCG.

Exploring the link between MORF4L1 and risk of breast cancer.

INTRODUCTION: Proteins encoded by Fanconi anemia (FA) and/or breast cancer (BrCa) susceptibility genes cooperate in a common DNA damage repair signaling pathway. To gain deeper insight into this pathway and its influence on cancer risk, we searched for novel components through protein physical interaction screens. METHODS: Protein physical interactions were screened using the yeast two-hybrid system. Co-affinity purifications and endogenous co-immunoprecipitation assays were performed to corroborate interactions. Biochemical and functional assays in human, mouse and Caenorhabditis elegans models were carried out to characterize pathway components. Thirteen FANCD2-monoubiquitinylation-positive FA cell lines excluded for genetic defects in the downstream pathway components and 300 familial BrCa patients negative for BRCA1/2 mutations were analyzed for genetic mutations. Common genetic variants were genotyped in 9,573 BRCA1/2 mutation carriers for associations with BrCa risk. RESULTS: A previously identified co-purifying protein with PALB2 was identified, MRG15 (MORF4L1 gene). Results in human, mouse and C. elegans models delineate molecular and functional relationships with BRCA2, PALB2, RAD51 and RPA1 that suggest a role for MRG15 in the repair of DNA double-strand breaks. Mrg15-deficient murine embryonic fibroblasts showed moderate sensitivity to γ-irradiation relative to controls and reduced formation of Rad51 nuclear foci. Examination of mutants of MRG15 and BRCA2 C. elegans orthologs revealed phenocopy by accumulation of RPA-1 (human RPA1) nuclear foci and aberrant chromosomal compactions in meiotic cells. However, no alterations or mutations were identified for MRG15/MORF4L1 in unclassified FA patients and BrCa familial cases. Finally, no significant associations between common MORF4L1 variants and BrCa risk for BRCA1 or BRCA2 mutation carriers were identified: rs7164529, Ptrend = 0.45 and 0.05, P2df = 0.51 and 0.14, respectively; and rs10519219, Ptrend = 0.92 and 0.72, P2df = 0.76 and 0.07, respectively. CONCLUSIONS: While the present study expands on the role of MRG15 in the control of genomic stability, weak associations cannot be ruled out for potential low-penetrance variants at MORF4L1 and BrCa risk among BRCA2 mutation carriers.

Evidence for a link between TNFRSF11A and risk of breast cancer.

Intracellular signaling mediated by the receptor activator of nuclear factor-κB [Rank, encoded by the tumor necrosis factor receptor superfamily, member 11a (Tnfrsf11a) gene] is fundamental for mammary gland development in mice, regulating the expansion of stem and progenitor cell compartments. Conversely, Rank overexpression in mice promotes abnormal proliferation and impairs differentiation, leading to an increased incidence of tumorigenesis. Here, we show that a common genetic variant near the 5'-end of TNFRSF11A, rs7226991, is associated with breast cancer risk in the general population and among carriers of mutations in the breast cancer 2, early onset (BRCA2) gene. Akin to the results of the Cancer and Genetics Markers of Susceptibility initiative, combined analysis of rs7226991 in two Spanish case-control studies (1,365 controls and 1,323 cases in total) revealed a significant association with risk: odds ratio (OR) = 0.88, 95% confidence interval (CI) 0.78-0.98, P (trend) = 0.025. Subsequent examination of BRCA1 (n = 1,017) and BRCA2 (n = 885) mutation carriers revealed a consistent association in the latter group: weighted hazard ratio ((w)HR) = 0.70; 95% CI 0.55-0.88; and P (trend) = 0.003; compared to BRCA1 mutation carriers, (w)HR = 0.91; 95% CI 0.76-1.10; and P (trend) = 0.33. The results of this study need to be replicated in other populations and with larger numbers of BRCA1/2 mutation carriers.

Gene set-based analysis of polymorphisms: finding pathways or biological processes associated to traits in genome-wide association studies.

Genome-wide association studies have become a popular strategy to find associations of genes to traits of interest. Despite the high-resolution available today to carry out genotyping studies, the success of its application in real studies has been limited by the testing strategy used. As an alternative to brute force solutions involving the use of very large cohorts, we propose the use of the Gene Set Analysis (GSA), a different analysis strategy based on testing the association of modules of functionally related genes. We show here how the Gene Set-based Analysis of Polymorphisms (GeSBAP), which is a simple implementation of the GSA strategy for the analysis of genome-wide association studies, provides a significant increase in the power testing for this type of studies. GeSBAP is freely available at http://bioinfo.cipf.es/gesbap/.

Non-Lynch Familial and Early-Onset Colorectal Cancer Explained by Accumulation of Low-Risk Genetic Variants.

A large proportion of familial and/or early-onset cancer patients do not carry pathogenic variants in known cancer predisposing genes. We aimed to assess the contribution of previously validated low-risk colorectal cancer (CRC) alleles to familial/early-onset CRC (fCRC) and to serrated polyposis. We estimated the association of CRC with a 92-variant-based weighted polygenic risk score (wPRS) using 417 fCRC patients, 80 serrated polyposis patients, 1077 hospital-based incident CRC patients, and 1642 controls. The mean wPRS was significantly higher in fCRC than in controls or sporadic CRC patients. fCRC patients in the highest (20th) wPRS quantile were at four-fold greater CRC risk than those in the middle quantile (10th). Compared to low-wPRS fCRC, a higher number of high-wPRS fCRC patients had developed multiple primary CRCs, had CRC family history, and were diagnosed at age ≥50. No association with wPRS was observed for serrated polyposis. In conclusion, a relevant proportion of mismatch repair (MMR)-proficient fCRC cases might be explained by the accumulation of low-risk CRC alleles. Validation in independent cohorts and development of predictive models that include polygenic risk score (PRS) data and other CRC predisposing factors will determine the implementation of PRS into genetic testing and counselling in familial and early-onset CRC.

Evolutionary Changes after Translational Challenges Imposed by Horizontal Gene Transfer.

Genes acquired by horizontal gene transfer (HGT) may provide the recipient organism with potentially new functions, but proper expression level and integration of the transferred genes in the novel environment are not granted. Notably, transferred genes can differ from the receiving genome in codon usage preferences, leading to impaired translation and reduced functionality. Here, we characterize the genomic and proteomic changes undergone during experimental evolution of Escherichia coli after HGT of three synonymous versions, presenting very different codon usage preference, of an antibiotic resistance gene. The experimental evolution was conducted with and without the corresponding antibiotic and the mutational patterns and proteomic profiles after 1,000 generations largely depend on the experimental growth conditions (e.g., mutations in antibiotic off-target genes), and on the synonymous gene version transferred (e.g., mutations in genes responsive to translational stress). The transfer of an exogenous gene extensively modifies the whole proteome, and these proteomic changes are different for the different version of the transferred gene. Additionally, we identified conspicuous changes in global regulators and in intermediate metabolism, confirmed the evolutionary ratchet generated by mutations in DNA repair genes and highlighted the plasticity of bacterial genomes accumulating large and occasionally transient duplications. Our results support a central role of HGT in fuelling evolution as a powerful mechanism promoting rapid, often dramatic genotypic and phenotypic changes. The profound reshaping of the pre-existing geno/phenotype allows the recipient bacteria to explore new ways of functioning, far beyond the mere acquisition of a novel function.

Assessing associations between the AURKA-HMMR-TPX2-TUBG1 functional module and breast cancer risk in BRCA1/2 mutation carriers.

While interplay between BRCA1 and AURKA-RHAMM-TPX2-TUBG1 regulates mammary epithelial polarization, common genetic variation in HMMR (gene product RHAMM) may be associated with risk of breast cancer in BRCA1 mutation carriers. Following on these observations, we further assessed the link between the AURKA-HMMR-TPX2-TUBG1 functional module and risk of breast cancer in BRCA1 or BRCA2 mutation carriers. Forty-one single nucleotide polymorphisms (SNPs) were genotyped in 15,252 BRCA1 and 8,211 BRCA2 mutation carriers and subsequently analyzed using a retrospective likelihood approach. The association of HMMR rs299290 with breast cancer risk in BRCA1 mutation carriers was confirmed: per-allele hazard ratio (HR) = 1.10, 95% confidence interval (CI) 1.04-1.15, p = 1.9 x 10(-4) (false discovery rate (FDR)-adjusted p = 0.043). Variation in CSTF1, located next to AURKA, was also found to be associated with breast cancer risk in BRCA2 mutation carriers: rs2426618 per-allele HR = 1.10, 95% CI 1.03-1.16, p = 0.005 (FDR-adjusted p = 0.045). Assessment of pairwise interactions provided suggestions (FDR-adjusted pinteraction values > 0.05) for deviations from the multiplicative model for rs299290 and CSTF1 rs6064391, and rs299290 and TUBG1 rs11649877 in both BRCA1 and BRCA2 mutation carriers. Following these suggestions, the expression of HMMR and AURKA or TUBG1 in sporadic breast tumors was found to potentially interact, influencing patients' survival. Together, the results of this study support the hypothesis of a causative link between altered function of AURKA-HMMR-TPX2-TUBG1 and breast carcinogenesis in BRCA1/2 mutation carriers.