Spatially resolved proteomic profiling identifies tumor cell CD44 as a biomarker associated with sensitivity to PD-1 axis blockade in advanced non-small-cell lung cancer

BACKGROUND

Most patients with advanced non-small-cell lung cancer (NSCLC) fail to derive significant benefit from programmed cell death protein-1 (PD-1) axis blockade, and new biomarkers of response are needed. In this study, we aimed to discover and validate spatially resolved protein markers associated with sensitivity to PD-1 axis inhibition in NSCLC.

METHODS

We initially assessed a discovery cohort of 56 patients with NSCLC treated with PD-1 axis inhibitors at Yale Cancer Center. Using the GeoMx Digital Spatial Profiling (DSP) system, 71 proteins were measured in spatial context on each spot in a tissue microarray. We used the AQUA method of quantitative immunofluorescence (QIF) to orthogonally validate candidate biomarkers. For external independent validation, we assessed whole tissue sections derived from 128 patients with NSCLC treated with single-agent PD-1 axis inhibitors at the 12 de Octubre Hospital (Madrid) using DSP. We further analyzed two immunotherapy untreated cohorts to address prognostic significance (n=252 from Yale Cancer Center; n=124 from University Clinic of Navarra) using QIF and DSP, respectively.

RESULTS

Using continuous log-scaled data, we identified CD44 expression in the tumor compartment (pan-cytokeratin (CK)+) as a novel predictor of prolonged progression-free survival (PFS) (multivariate HR=0.68, p=0.043) in the discovery set. We validated by QIF that tumor CD44 levels assessed as continuous QIF scores were associated with longer PFS (multivariate HR=0.31, p=0.022) and overall survival (multivariate HR=0.29, p=0.038). Using DSP in an independent immunotherapy treated cohort, we validated that CD44 levels in the tumor compartment, but not in the immune compartment (panCK-/CD45+), were associated with clinical benefit (OR=1.22, p=0.018) and extended PFS under PD-1 axis inhibition using the highest tertile cutpoint (multivariate HR=0.62, p=0.03). The effect of tumor cell CD44 in predicting PFS remained significant after correcting for programmed death-ligand 1 (PD-L1) Tumor Proportion Score (TPS) in both cohorts. High tumor cell CD44 was not prognostic in the absence of immunotherapy. Using DSP data, intratumoral regions with elevated tumor cell CD44 expression showed prominent (fold change>1.5, adjusted p<0.05) upregulation of PD-L1, TIM-3, ICOS, and CD40 in two independent cohorts.

CONCLUSIONS

This work highlights CD44 as a novel indicative biomarker of sensitivity to PD-1 axis blockade that might help to improve immunotherapy strategies for NSCLC.

Abstract 2028: Spatially resolved proteomic profiling identifies tumor cell CD44 as a novel indicative biomarker of sensitivity to PD-1 axis blockade in advanced non-small cell lung cancer (NSCLC)

Introduction: Most patients with advanced NSCLC fail to develop durable responses from PD-1 axis blockade, and more robust predictive biomarkers are needed. In this study, we aimed to discover and validate spatially resolved protein markers associated with sensitivity to PD-1 axis inhibition across multiple NSCLC cohorts.

Methods: We first used the GeoMx Digital Spatial Profiling (DSP) system in a discovery cohort of 56 NSCLC patients treated with PD-1 axis inhibitors at Yale (USA). Pre-treatment tumors were represented in a tissue microarray (YTMA471) and analyzed using a 71-plex primary antibody panel. Proteins were measured from three molecularly defined tissue compartments: tumor (CK+), leukocytes (CD45+/CD68-), and macrophages (CD68+). We used quantitative immunofluorescence (QIF) to orthogonally validate candidate biomarkers. For external validation of identified predictors, we assessed pre-treatment whole tissue sections from a cohort of 128 NSCLC patients treated with single-agent PD-1 axis inhibitors at the Hospital 12 de Octubre (Madrid, Spain) using DSP (39-plex protein panel, measured from CK+ and CD45+ compartments). In addition, we analyzed two immunotherapy untreated cohorts to address prognostic significance: YTMA423 (Yale, USA; n = 252) and CIMA-CUN (UNAV, Spain; n = 124), using QIF and DSP respectively.

Results: Using continuous log-scaled data, we found 5 markers independently associated with progression-free survival (PFS) in the tumor compartment (including PD-L1, p = 0.002). Among the novel candidate predictors, tumor cell CD44, a marker of pluripotency and stemness, was associated with longer PFS (multivariate HR = 0.63, p= 0.002). Using QIF, we orthogonally validated that CD44 expression in the tumor compartment was associated with longer PFS (p &lt;0.001) and overall survival (OS) (p = 0.03). In the external validation cohort, CD44 expressed in the tumor compartment, but not in the immune compartment, was predictive of clinical benefit (OR 1.22, p = 0.018), and was significantly associated with longer PFS (first tertile cutpoint: HR 0.62, p = 0.03). In contrast, high CD44 expression was not associated with survival in the two untreated cohorts. Using DSP data from two cohorts, regions of interest with elevated expression of CD44 in tumor cells consistently showed prominent (Fold Change&gt;1.5, p&lt;0.05) upregulation of TIM-3, and PD-L1 and multiple co-stimulatory molecules (including CD40, ICOS, CD27).

Conclusion: This work highlights CD44 as a novel indicative biomarker of sensitivity to PD-1 axis blockade in NSCLC. NSCLCs with high CD44 expression in epithelial cells may be associated with an immune contexture primed for higher response to immune checkpoint blockade. Further studies are needed to understand the interplay between CD44+ cancer stem cell phenotype and mechanisms of immune evasion.

Citation Format: Myrto Moutafi, Magdalena Molero, Sandra Martinez-Morilla, Javier Baena, Ioannis Vathiotis, Niki Gavrielatou, Laura Castro, Gorka Ruiz de Garibay, Vera Adradas, Daniel Orive, Karmele Valencia, Alfonso Galvo, Luis M. Montuenga, Santiago Ponce, Kurt Schalper, Luis Paz-Ares, David L. Rimm, Jon Zugazagoitia. Spatially resolved proteomic profiling identifies tumor cell CD44 as a novel indicative biomarker of sensitivity to PD-1 axis blockade in advanced non-small cell lung cancer (NSCLC) [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 2028.

Histamine signaling and metabolism identify potential biomarkers and therapies for lymphangioleiomyomatosis

Inhibition of mTOR is the standard of care for lymphangioleiomyomatosis (LAM). However, this therapy has variable tolerability and some patients show progressive decline of lung function despite treatment. LAM diagnosis and monitoring can also be challenging due to the heterogeneity of symptoms and insufficiency of non-invasive tests. Here, we propose monoamine-derived biomarkers that provide preclinical evidence for novel therapeutic approaches. The major histamine-derived metabolite methylimidazoleacetic acid (MIAA) is relatively more abundant in LAM plasma, and MIAA values are independent of VEGF-D. Higher levels of histamine are associated with poorer lung function and greater disease burden. Molecular and cellular analyses, and metabolic profiling confirmed active histamine signaling and metabolism. LAM tumorigenesis is reduced using approved drugs targeting monoamine oxidases A/B (clorgyline and rasagiline) or histamine H1 receptor (loratadine), and loratadine synergizes with rapamycin. Depletion of Maoa or Hrh1 expression, and administration of an L-histidine analog, or a low L-histidine diet, also reduce LAM tumorigenesis. These findings extend our knowledge of LAM biology and suggest possible ways of improving disease management.

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.

Repurposing 18F-FMISO as a PET tracer for translational imaging of nitroreductase-based gene directed enzyme prodrug therapy

Nitroreductases (NTR) are a family of bacterial enzymes used in gene directed enzyme prodrug therapy (GDEPT) that selectively activate prodrugs containing aromatic nitro groups to exert cytotoxic effects following gene transduction in tumours. The clinical development of NTR-based GDEPT has, in part, been hampered by the lack of translational imaging modalities to assess gene transduction and drug cytotoxicity, non-invasively. This study presents translational preclinical PET imaging to validate and report NTR activity using the clinically approved radiotracer, ¹⁸F-FMISO, as substrate for the NTR enzyme.

Methods: The efficacy with which ¹⁸F-FMISO could be used to report NfsB NTR activity in vivo was investigated using the MDA-MB-231 mammary carcinoma xenograft model. For validation, subcutaneous xenografts of cells constitutively expressing NTR were imaged using ¹⁸F-FMISO PET/CT and fluorescence imaging with CytoCy5S, a validated fluorescent NTR substrate. Further, examination of the non-invasive functionality of ¹⁸F-FMISO PET/CT in reporting NfsB NTR activity in vivo was assessed in metastatic orthotopic NfsB NTR expressing xenografts and metastasis confirmed by bioluminescence imaging. ¹⁸F-FMISO biodistribution was acquired ex vivo by an automatic gamma counter measuring radiotracer retention to confirm in vivo results. To assess the functional imaging of NTR-based GDEPT with ¹⁸F-FMISO, PET/CT was performed to assess both gene transduction and cytotoxicity effects of prodrug therapy (CB1954) in subcutaneous models.

Results:¹⁸F-FMISO retention was detected in NTR⁺ subcutaneous xenografts, displaying significantly higher PET contrast than NTR^- xenografts (p < 0.0001). Substantial ¹⁸F-FMISO retention was evident in metastases of orthotopic xenografts (p < 0.05). Accordingly, higher ¹⁸F-FMISO biodistribution was prevalent ex vivo in NTR⁺ xenografts. ¹⁸F-FMISO NfsB NTR PET/CT imaging proved useful for monitoring in vivo NTR transduction and the cytotoxic effect of prodrug therapy.

Conclusions:¹⁸F-FMISO NfsB NTR PET/CT imaging offered significant contrast between NTR⁺ and NTR^- tumours and effective resolution of metastatic progression. Furthermore, ¹⁸F-FMISO NfsB NTR PET/CT imaging proved efficient in monitoring the two steps of GDEPT, in vivo NfsB NTR transduction and response to CB1954 prodrug therapy. These results support the repurposing of ¹⁸F-FMISO as a readily implementable PET imaging probe to be employed as companion diagnostic test for NTR-based GDEPT systems.

Selecting the optimal parameters for sonoporation of pancreatic cancer in a pre-clinical model

Pancreatic ductal adenocarcinoma (PDAC) is one of the deadliest cancers in the modern world, in part due to poor delivery of chemotherapeutics. Sonoporation can be used to enhance the efficacy of standard of care therapies for PDAC. Using xenograft models of PDAC we investigate sonoporation using four ifferent ultrasound contrast agents (UCAs) and two ultrasound regimens to identify the ideal parameters to increase therapeutic efficacy. MIA-PaCa2 xenografts in over 175 immunodeficient mice were treated with gemcitabine and paclitaxel and subjected to low or high power ultrasound (60 and 200 mW/cm² respectively) in conjunction with one of four different UCAs. The UCAs investigated were Definity®, SonoVue®, Optison™ or Sonazoid™. Tumor volumes, vascularity, hemoglobin, and oxygenation were measured and compared to controls. High power treatment in conjunction with Sonazoid sonoporation led to significantly smaller tumors when started early (tumors ~50mm³; p = .0105), while no UCAs significantly increased efficacy in the low power cohort. This trend was also found in larger tumors (~250mm³) where all four UCA agents significantly increased therapeutic efficacy in the high power group (p < .01), while only Definity and SonoVue increased efficacy in the low power cohort (p < .03). Overall, the higher power ultrasound treatment modality was more consistently effective at decreasing tumor volume and increasing vascularity characteristics. In conclusion, Sonazoid was the most consistently effective UCA at decreasing tumor volume and increasing vascularity. Thus, we are pursuing a larger phase II clinical trial to validate the increased efficacy of sonoporation in conjunction with chemotherapy in PDAC patients.

Low-Intensity Sonoporation-Induced Intracellular Signalling of Pancreatic Cancer Cells, Fibroblasts and Endothelial Cells

The use of ultrasound (US) and microbubbles (MB), usually referred to as sonoporation, has great potential to increase the efficacy of chemotherapy. However, the molecular mechanisms that mediate sonoporation response are not well-known, and recent research suggests that cell stress induced by US + MBs may contribute to the treatment benefit. Furthermore, there is a growing understanding that the effects of US + MBs are beyond only the cancer cells and involves the tumour vasculature and microenvironment. We treated pancreatic cancer cells (MIA PaCa-2) and stromal cells, fibroblasts (BJ) and human umbilical vein endothelial cells (HUVECs), with US ± MB, and investigated the extent of uptake of cell impermeable dye (calcein, by flow cytometry), viability (cell count, Annexin/PI and WST-1 assays) and activation of a number of key proteins in important intracellular signalling pathways immediately and 2 h after sonoporation (phospho flow cytometry). Different cell types responded differently to US ± MBs in all these aspects. In general, sonoporation induces immediate, transient activation of MAP-kinases (p38, ERK1/2), and an increase in phosphorylation of ribosomal protein S6 together with dephosphorylation of 4E-BP1. The sonoporation stress-response resembles cellular responses to electroporation and pore-forming toxins in membrane repair and restoring cellular homeostasis, and may be exploited therapeutically. The stromal cells were more sensitive to sonoporation than tumoural cells, and further efforts in optimising sonoporation-enhanced therapy should be targeted at the microenvironment.

Intracellular Cytidine Deaminase Regulates Gemcitabine Metabolism in Pancreatic Cancer Cell Lines

Cytidine deaminase (CDA) is a determinant of in vivo gemcitabine elimination kinetics and cellular toxicity. The impact of CDA activity in pancreatic ductal adenocarcinoma (PDAC) cell lines has not been elucidated. We hypothesized that CDA regulates gemcitabine flux through its inactivation and activation pathways in PDAC cell lines. Three PDAC cell lines (BxPC-3, MIA PaCa-2, and PANC-1) were incubated with 10 or 100 µM gemcitabine for 60 minutes or 24 hours, with or without tetrahydrouridine, a CDA inhibitor. Extracellular inactive gemcitabine metabolite (dFdU) and intracellular active metabolite (dFdCTP) were quantified with liquid chromatography tandem mass spectrometry. Cellular expression of CDA was assessed with real-time PCR and Western blot. Gemcitabine conversion to dFdU was extensive in BxPC-3 and low in MIA PaCa-2 and PANC-1, in accordance with their respective CDA expression levels. CDA inhibition was associated with low or undetectable dFdU in all three cell lines. After 24 hours gemcitabine incubation, dFdCTP was highest in MIA PaCa-2 and lowest in BxPC-3. CDA inhibition resulted in a profound dFdCTP increase in BxPC-3 but not in MIA PaCa-2 or PANC-1. dFdCTP concentrations were not higher after exposure to 100 versus 10 µM gemcitabine when CDA activities were low (MIA PaCa-2 and PANC-1) or inhibited (BxPC-3). The results suggest a regulatory role of CDA for gemcitabine activation in PDAC cells but within limits related to the capacity in the activation pathway in the cell lines. SIGNIFICANCE STATEMENT: The importance of cytidine deaminase (CDA) for cellular gemcitabine toxicity, linking a lower activity to higher toxicity, is well described. An underlying assumption is that CDA, by inactivating gemcitabine, limits the amount available for the intracellular activation pathway. Our study is the first to illustrate this regulatory role of CDA in pancreatic ductal adenocarcinoma cell lines by quantifying intracellular and extracellular gemcitabine metabolite concentrations.

Tumor xenograft modeling identifies an association between TCF4 loss and breast cancer chemoresistance

Understanding the mechanisms of cancer therapeutic resistance is fundamental to improving cancer care. There is clear benefit from chemotherapy in different breast cancer settings; however, knowledge of the mutations and genes that mediate resistance is incomplete. In this study, by modeling chemoresistance in patient-derived xenografts (PDXs), we show that adaptation to therapy is genetically complex and identify that loss of transcription factor 4 (TCF4; also known as ITF2) is associated with this process. A triple-negative BRCA1-mutated PDX was used to study the genetics of chemoresistance. The PDX was treated in parallel with four chemotherapies for five iterative cycles. Exome sequencing identified few genes with de novo or enriched mutations in common among the different therapies, whereas many common depleted mutations/genes were observed. Analysis of somatic mutations from The Cancer Genome Atlas (TCGA) supported the prognostic relevance of the identified genes. A mutation in TCF4 was found de novo in all treatments, and analysis of drug sensitivity profiles across cancer cell lines supported the link to chemoresistance. Loss of TCF4 conferred chemoresistance in breast cancer cell models, possibly by altering cell cycle regulation. Targeted sequencing in chemoresistant tumors identified an intronic variant of TCF4 that may represent an expression quantitative trait locus associated with relapse outcome in TCGA. Immunohistochemical studies suggest a common loss of nuclear TCF4 expression post-chemotherapy. Together, these results from tumor xenograft modeling depict a link between altered TCF4 expression and breast cancer chemoresistance.

Summary: By modeling chemoresistance in patient-derived breast cancer xenografts, this study shows that adaptation to therapy is genetically complex and that loss of transcription factor 4 (TCF4) is associated with this process.

Decapping protein EDC4 regulates DNA repair and phenocopies BRCA1

BRCA1 is a tumor suppressor that regulates DNA repair by homologous recombination. Germline mutations in BRCA1 are associated with increased risk of breast and ovarian cancer and BRCA1 deficient tumors are exquisitely sensitive to poly (ADP-ribose) polymerase (PARP) inhibitors. Therefore, uncovering additional components of this DNA repair pathway is of extreme importance for further understanding cancer development and therapeutic vulnerabilities. Here, we identify EDC4, a known component of processing-bodies and regulator of mRNA decapping, as a member of the BRCA1-BRIP1-TOPBP1 complex. EDC4 plays a key role in homologous recombination by stimulating end resection at double-strand breaks. EDC4 deficiency leads to genome instability and hypersensitivity to DNA interstrand cross-linking drugs and PARP inhibitors. Lack-of-function mutations in EDC4 were detected in BRCA1/2-mutation-negative breast cancer cases, suggesting a role in breast cancer susceptibility. Collectively, this study recognizes EDC4 with a dual role in decapping and DNA repair whose inactivation phenocopies BRCA1 deficiency.

Cancer network activity associated with therapeutic response and synergism

Background

Cancer patients often show no or only modest benefit from a given therapy. This major problem in oncology is generally attributed to the lack of specific predictive biomarkers, yet a global measure of cancer cell activity may support a comprehensive mechanistic understanding of therapy efficacy. We reasoned that network analysis of omic data could help to achieve this goal.

Methods

A measure of “cancer network activity” (CNA) was implemented based on a previously defined network feature of communicability. The network nodes and edges corresponded to human proteins and experimentally identified interactions, respectively. The edges were weighted proportionally to the expression of the genes encoding for the corresponding proteins and relative to the number of direct interactors. The gene expression data corresponded to the basal conditions of 595 human cancer cell lines. Therapeutic responses corresponded to the impairment of cell viability measured by the half maximal inhibitory concentration (IC₅₀) of 130 drugs approved or under clinical development. Gene ontology, signaling pathway, and transcription factor-binding annotations were taken from public repositories. Predicted synergies were assessed by determining the viability of four breast cancer cell lines and by applying two different analytical methods.

Results

The effects of drug classes were associated with CNAs formed by different cell lines. CNAs also differentiate target families and effector pathways. Proteins that occupy a central position in the network largely contribute to CNA. Known key cancer-associated biological processes, signaling pathways, and master regulators also contribute to CNA. Moreover, the major cancer drivers frequently mediate CNA and therapeutic differences. Cell-based assays centered on these differences and using uncorrelated drug effects reveals novel synergistic combinations for the treatment of breast cancer dependent on PI3K-mTOR signaling.

Conclusions

Cancer therapeutic responses can be predicted on the basis of a systems-level analysis of molecular interactions and gene expression. Fundamental cancer processes, pathways, and drivers contribute to this feature, which can also be exploited to predict precise synergistic drug combinations.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-016-0340-x) contains supplementary material, which is available to authorized users.

Study of breast cancer incidence in patients of lymphangioleiomyomatosis

Molecular evidence has linked the pathophysiology of lymphangioleiomyomatosis (LAM) to that of metastatic breast cancer. Following on this observation, we assessed the association between LAM and subsequent breast cancer. An epidemiological study was carried out using three LAM country cohorts, from Japan, Spain, and the United Kingdom. The number of incident breast cancer cases observed in these cohorts was compared with the number expected on the basis of the country-specific incidence rates for the period 2000–2014. Immunohistochemical studies and exome sequence analysis were performed in two and one tumors, respectively. All cohorts revealed breast cancer standardized incidence ratios (SIRs) ≥ 2.25. The combined analysis of all cases or restricted to pre-menopausal age groups revealed significantly higher incidence of breast cancer: SIR = 2.81, 95 % confidence interval (CI) = 1.32–5.57, P = 0.009; and SIR = 4.88, 95 % CI = 2.29–9.99, P = 0.0007, respectively. Immunohistochemical analyses showed positivity for known markers of lung metastatic potential. This study suggests the existence of increased breast cancer risk among LAM patients. Prospective studies may be warranted to corroborate this result, which may be particularly relevant for pre-menopausal women with LAM.

Lymphangioleiomyomatosis Biomarkers Linked to Lung Metastatic Potential and Cell Stemness

Lymphangioleiomyomatosis (LAM) is a rare lung-metastasizing neoplasm caused by the proliferation of smooth muscle-like cells that commonly carry loss-of-function mutations in either the tuberous sclerosis complex 1 or 2 (TSC1 or TSC2) genes. While allosteric inhibition of the mechanistic target of rapamycin (mTOR) has shown substantial clinical benefit, complementary therapies are required to improve response and/or to treat specific patients. However, there is a lack of LAM biomarkers that could potentially be used to monitor the disease and to develop other targeted therapies. We hypothesized that the mediators of cancer metastasis to lung, particularly in breast cancer, also play a relevant role in LAM. Analyses across independent breast cancer datasets revealed associations between low TSC1/2 expression, altered mTOR complex 1 (mTORC1) pathway signaling, and metastasis to lung. Subsequently, immunohistochemical analyses of 23 LAM lesions revealed positivity in all cases for the lung metastasis mediators fascin 1 (FSCN1) and inhibitor of DNA binding 1 (ID1). Moreover, assessment of breast cancer stem or luminal progenitor cell biomarkers showed positivity in most LAM tissue for the aldehyde dehydrogenase 1 (ALDH1), integrin-ß3 (ITGB3/CD61), and/or the sex-determining region Y-box 9 (SOX9) proteins. The immunohistochemical analyses also provided evidence of heterogeneity between and within LAM cases. The analysis of Tsc2-deficient cells revealed relative over-expression of FSCN1 and ID1; however, Tsc2-deficient cells did not show higher sensitivity to ID1-based cancer inhibitors. Collectively, the results of this study reveal novel LAM biomarkers linked to breast cancer metastasis to lung and to cell stemness, which in turn might guide the assessment of additional or complementary therapeutic opportunities for LAM.

BRCA1 Alternative splicing landscape in breast tissue samples

Background

BRCA1 is a key protein in cell network, involved in DNA repair pathways and cell cycle. Recently, the ENIGMA consortium has reported a high number of alternative splicing (AS) events at this locus in blood-derived samples. However, BRCA1 splicing pattern in breast tissue samples is unknown. Here, we provide an accurate description of BRCA1 splicing events distribution in breast tissue samples.

Methods

BRCA1 splicing events were scanned in 70 breast tumor samples, 4 breast samples from healthy individuals and in 72 blood-derived samples by capillary electrophoresis (capillary EP). Molecular subtype was identified in all tumor samples. Splicing events were considered predominant if their relative expression level was at least the 10% of the full-length reference signal.

Results

54 BRCA1 AS events were identified, 27 of them were annotated as predominant in at least one sample. Δ5q, Δ13, Δ9, Δ5 and ▼1aA were significantly more frequently annotated as predominant in breast tumor samples than in blood-derived samples. Predominant splicing events were, on average, more frequent in tumor samples than in normal breast tissue samples (P = 0.010). Similarly, likely inactivating splicing events (PTC-NMDs, Non-Coding, Δ5 and Δ18) were more frequently annotated as predominant in tumor than in normal breast samples (P = 0.020), whereas there were no significant differences for other splicing events (No-Fs) frequency distribution between tumor and normal breast samples (P = 0.689).

Conclusions

Our results complement recent findings by the ENIGMA consortium, demonstrating that BRCA1 AS, despite its tremendous complexity, is similar in breast and blood samples, with no evidences for tissue specific AS events. Further on, we conclude that somatic inactivation of BRCA1 through spliciogenic mutations is, at best, a rare mechanism in breast carcinogenesis, albeit our data detects an excess of likely inactivating AS events in breast tumor samples.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-015-1145-9) contains supplementary material, which is available to authorized users.

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.

Genetic polymorphisms in estrogen-related genes and the risk of breast cancer among Han Chinese women

Exposure to high levels of estrogen is considered an important risk factor for susceptibility to breast cancer. Common polymorphisms in genes that affect estrogen levels may be associated with breast cancer risk, but no comprehensive study has been performed among Han Chinese women. In the present study, 32 single-nucleotide polymorphisms (SNPs) in estrogen-related genes were genotyped using the MassARRAY IPLEX platform in 1076 Han Chinese women. Genotypic and allelic frequencies were compared between case and control groups. Unconditional logistic regression was used to assess the effects of SNPs on breast cancer risk. Associations were also evaluated for breast cancer subtypes stratified by estrogen receptor (ER) and progesterone receptor (PR) status. Case-control analysis showed a significant relation between heterozygous genotypes of rs700519 and rs2069522 and breast cancer risk (OR = 0.723, 95% CI = 0.541-0.965, p = 0.028 and OR = 1.500, 95% CI = 1.078-2.087, p = 0.016, respectively). Subgroup comparisons revealed that rs2446405 and rs17268974 were related to ER status, and rs130021 was associated with PR status. Our findings suggest that rs700519 and rs2069522 are associated with susceptibility to breast cancer among the Han Chinese population and have a cumulative effect with three other identified SNPs. Further genetic and functional studies are needed to identify additional SNPs, and to elucidate the underlying molecular mechanisms.

An A91V SNP in the perforin gene is frequently found in NK/T-cell lymphomas

NK/T-cell lymphoma (NKTCL) is the most frequent EBV-related NK/T-cell disease. Its clinical manifestations overlap with those of familial haemophagocytic lymphohistiocytosis (FHLH). Since PERFORIN (PRF1) mutations are present in FHLH, we analysed its role in a series of 12 nasal and 12 extranasal-NKTCLs. 12.5% of the tumours and 25% of the nasal-origin cases had the well-known g.272C>T(p.Ala91Val) pathogenic SNP, which confers a poor prognosis. Two of these cases had a double-CD4/CD8-positive immunophenotype, although no correlation was found with perforin protein expression. p53 was overexpressed in 20% of the tumoral samples, 80% of which were of extranasal origin, while none showed PRF1 SNVs. These results suggest that nasal and extranasal NKTCLs have different biological backgrounds, although this requires validation.

Comparison of mRNA splicing assay protocols across multiple laboratories: recommendations for best practice in standardized clinical testing

BACKGROUND

Accurate evaluation of unclassified sequence variants in cancer predisposition genes is essential for clinical management and depends on a multifactorial analysis of clinical, genetic, pathologic, and bioinformatic variables and assays of transcript length and abundance. The integrity of assay data in turn relies on appropriate assay design, interpretation, and reporting.

METHODS

We conducted a multicenter investigation to compare mRNA splicing assay protocols used by members of the ENIGMA (Evidence-Based Network for the Interpretation of Germline Mutant Alleles) consortium. We compared similarities and differences in results derived from analysis of a panel of breast cancer 1, early onset (BRCA1) and breast cancer 2, early onset (BRCA2) gene variants known to alter splicing (BRCA1: c.135-1G>T, c.591C>T, c.594-2A>C, c.671-2A>G, and c.5467+5G>C and BRCA2: c.426-12_8delGTTTT, c.7988A>T, c.8632+1G>A, and c.9501+3A>T). Differences in protocols were then assessed to determine which elements were critical in reliable assay design.

RESULTS

PCR primer design strategies, PCR conditions, and product detection methods, combined with a prior knowledge of expected alternative transcripts, were the key factors for accurate splicing assay results. For example, because of the position of primers and PCR extension times, several isoforms associated with BRCA1, c.594-2A>C and c.671-2A>G, were not detected by many sites. Variation was most evident for the detection of low-abundance transcripts (e.g., BRCA2 c.8632+1G>A Δ19,20 and BRCA1 c.135-1G>T Δ5q and Δ3). Detection of low-abundance transcripts was sometimes addressed by using more analytically sensitive detection methods (e.g., BRCA2 c.426-12_8delGTTTT ins18bp).

CONCLUSIONS

We provide recommendations for best practice and raise key issues to consider when designing mRNA assays for evaluation of unclassified sequence variants.

Capillary electrophoresis analysis of conventional splicing assays: IARC analytical and clinical classification of 31 BRCA2 genetic variants

Rare sequence variants in "high-risk" disease genes, often referred as unclassified variants (UVs), pose a serious challenge to genetic testing. However, UVs resulting in splicing alterations can be readily assessed by in vitro assays. Unfortunately, analytical and clinical interpretation of these assays is often challenging. Here, we explore this issue by conducting splicing assays in 31 BRCA2 genetic variants. All variants were assessed by RT-PCR followed by capillary electrophoresis and direct sequencing. If assays did not produce clear-cut outputs (Class-2 or Class-5 according to analytical International Agency for Research on Cancer guidelines), we performed qPCR and/or minigene assays. The latter were performed with a new splicing vector (pSAD) developed by authors of the present manuscript (patent #P201231427 CSIC). We have identified three clinically relevant Class-5 variants (c.682-2A>G, c.7617+1G>A, and c.8954-5A>G), and 27 analytical Class-2 variants (not inducing splicing alterations). In addition, we demonstrate that rs9534262 (c.7806-14T>C) is a BRCA2 splicing quantitative trait locus.

Characterization of four novel BRCA2 large genomic rearrangements in Spanish breast/ovarian cancer families: review of the literature, and reevaluation of the genetic mechanisms involved in their origin

Large genomic rearrangements (LGRs) at the BRCA2 locus explain a non-negligible proportion of hereditary breast and ovarian cancer (HBOC) syndromes. The multiplex ligation and probe amplification (MLPA) assay has permitted in recent years to identify several families carrying LGRs at this locus, but very few such alterations have been fully characterized at the molecular level. Yet, molecular characterization is essential to identify recurrent alterations, to analyze the genetic mechanisms underlying such alterations, or to investigate potential genotype/phenotype relationships. We have used MLPA to identify BRCA2 LGRs in 7 out of 813 Spanish HBOC families previously tested negative for BRCA1 and BRCA2 small genomic alterations (substitutions and indels) and BRCA1 LGRs. We used a combination of long-range PCR, restriction mapping, and cDNA analysis to characterize the alterations at the molecular level. We found that Del Exon1-Exon2, Del Exon12-Exon16 and Del Exon22-Exon24 explain one family each, while Del Exon2 appears to be a Spanish founder mutation explaining four independent families. Finally, we have combined our data with a comprehensive review of the literature to reevaluate the genetic mechanisms underlying LGRs at the BRCA2 locus. Our study substantially increases the spectrum of BRCA2 LGRs fully characterized at the molecular level. Further on, we provide data to suggest that non-allelic homologous recombination has been overestimated as a mechanism underlying these alterations, while the opposite might be true for microhomology-mediated events.