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.

Modification of BRCA1-associated breast cancer risk by HMMR overexpression.

Breast cancer risk for carriers of BRCA1 pathological variants is modified by genetic factors. Genetic variation in HMMR may contribute to this effect. However, the impact of risk modifiers on cancer biology remains undetermined and the biological basis of increased risk is poorly understood. Here, we depict an interplay of molecular, cellular, and tissue microenvironment alterations that increase BRCA1-associated breast cancer risk. Analysis of genome-wide association results suggests that diverse biological processes, including links to BRCA1-HMMR profiles, influence risk. HMMR overexpression in mouse mammary epithelium increases Brca1-mutant tumorigenesis by modulating the cancer cell phenotype and tumor microenvironment. Elevated HMMR activates AURKA and reduces ARPC2 localization in the mitotic cell cortex, which is correlated with micronucleation and activation of cGAS-STING and non-canonical NF-κB signaling. The initial tumorigenic events are genomic instability, epithelial-to-mesenchymal transition, and tissue infiltration of tumor-associated macrophages. The findings reveal a biological foundation for increased risk of BRCA1-associated breast cancer.

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 (IC50) 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.

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.

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.

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.

About 1% of the breast and ovarian Spanish families testing negative for BRCA1 and BRCA2 are carriers of RAD51D pathogenic variants.

RAD51D mutations have been recently identified in breast (BC) and ovarian cancer (OC) families. Although an etiological role in OC appears to be present, the association of RAD51D mutations and BC risk is more unclear. We aimed to determine the prevalence of germline RAD51D mutations in Spanish BC/OC families negative for BRCA1/BRCA2 mutations. We analyzed 842 index patients: 491 from BC/OC families, 171 BC families, 51 OC families and 129 patients without family history but with early-onset BC or OC or metachronous BC and OC. Mutation detection was performed with high-resolution melting, denaturing high-performance liquid chromatography or Sanger sequencing. Three mutations were found in four families with BC and OC cases (0.82%). Two were novel: c.1A>T (p.Met1?) and c.667+2_667+23del, leading to the exon 7 skipping and one previously described: c.674C>T (p.Arg232*). All were present in BC/OC families with only one OC. The c.667+2_667+23del cosegregated in the family with one early-onset BC and two bilateral BC cases. We also identified the c.629C>T (p.Ala210Val) variant, which was predicted in silico to be potentially pathogenic. About 1% of the BC and OC Spanish families negative for BRCA1/BRCA2 are carriers of RAD51D mutations. The presence of several BC mutation carriers, albeit in the context of familial OC, suggests an increased risk for BC, which should be taken into account in the follow-up and early detection measures. RAD51D testing should be considered in clinical setting for families with BC and OC, irrespective of the number of OC cases in the family.

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.

Low prevalence of SLX4 loss-of-function mutations in non-BRCA1/2 breast and/or ovarian cancer families.

Fanconi anemia is a genetically heterogeneous autosomal recessive disorder characterized by development abnormalities, bone marrow failure, and childhood cancers. Compelling evidence indicates a common genetic basis for FA and breast/ovarian cancer susceptibility. Recently, biallelic germ-line mutations in SLX4 have been demonstrated to cause a previously unknown FA subtype (FA-P). We address the role of SLX4/FANCP in breast/ovarian cancer susceptibility by conducting a comprehensive mutation scanning in 486 index cases from non-BRCA1/BRCA2 multiple-case breast and/or ovarian cancer families (non-BRCA1/2 families) from Spain. We detected one unequivocal loss-of-function mutation (p.Glu1517X). In addition, one missense change (p.Arg372Trp) predicted to be pathogenic by in silico analysis co-segregates with disease in one family. Overall, the study indicates that SLX4 mutation screening will have a very low impact (if any) in the genetic counseling of non-BRCA1/2 families.