Serie, Daniel

Abstract Introduction: Invasive lobular carcinoma (ILC) comprises approximately ∼10-20% of breast cancers (BC) and appears to have a distinct biology from invasive ductal carcinoma (IDC). ILC is known to have a higher incidence of multifocal (MF) disease when compared to IDC (Biglia N, EJSO, 2013) and in general, MF BC has been associated with an increased rate of regional lymph node metastases (Cabioglu N, JACS, 2009). Treatment decisions are most often defined by the largest of the MF tumors without targeting potentially relevant biologic heterogeneity of the disease. We conducted a genomic analysis of MF ILC to help address this knowledge void. Methods: We characterized gene expression in multiple foci from 11 patients (pts) with MF ILC (all ER+ and HER2-, mean number of foci = 2.62) and 13 benign breast tissues. RNA was extracted from 3×1.5mm cores from each foci and gene expression was measured in 730 known cancer genes using the NanoString PanCancer pathways panel. Differential gene expression analyses were performed between 1) MF tumors and benign tissue (BT) accounting for pairing of pts and 2) different foci of MF tumors from the same ILC pt by ANOVA. The association of clinical covariates (age, number of tumors, grade and stage) in the top genes that were differentially expressed between multiple foci was assessed by linear mixed models accounting for pt and tumor. Results: A total of 20 genes (including MET and EGFR, top 5 genes: SPRY2, SOX17, FIGF, PDGFA, IL6) were differentially expressed (p<0.05) between malignant foci from the same pt, after correction by false discovery rate (FDR). EGFR expression was also negatively associated with the number of positive lymph nodes (N1a vs N0 absolute fold change = 0.51, p = 0.043). We observed an excess of genes that were down-regulated in MF ILC compared to BT. 277 genes were differentially expressed at p<0.05 after correction by FDR. The top 5 differentially expressed genes (by fold change) between ILC and BT were COL11A1, COL1A1, WIF1, FIGF and PROM1. CDH1, a known marker of ILC, was down-regulated in ILC compared to benign (absolute fold change = 2.54, p = 6.02×10−6), ranked only 39/277 genes by fold change. MET and EGFR ranked 29 and 42/277 genes (p = 5.99×10−9 and 3.32×10−8 respectively). Conclusion: We observed significant heterogeneity for 20/730 genes between individual foci of MF ILC within the same pt suggesting that MF synchronous ILC lesions may have different biologic behavior. This in turn may have important prognostic implications with potential therapeutic utility. Our analyses also shows significant down-regulation of gene expression between ILC tumors and BT for many genes other than CDH1 suggesting that additional biomarkers can be used to further characterize ILC Study funded in part by Mayo Clinic Cancer Focus Research Team, the 26.2 with Donna Foundation and the Serene M. and Frances C. Durling endowment. Citation Format: Pooja Advani, Nadine Norton, Daniel J. Serie, Xochiquetzal Geiger, Julia E. Crook, Aubrey E. Thompson, Edith A. Perez, Alvaro Moreno Aspitia. Characterization of individual foci of multifocal invasive lobular breast cancer using gene expression profiling. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4761. doi:10.1158/1538-7445.AM2015-4761

Abstract Genome-wide association studies (GWAS) of multiple myeloma (MM) in Northern Europeans have identified seven novel risk loci (2p23.3, 3p22.1, 3q26.2, 6p21.33, 7p15.3, 17p11.2, 22q13.1). We performed a multiethnic meta-analysis of these regions in 1,274 MM patients and 1,486 controls of European ancestry (EA) and 1,049 MM patients and 7,080 controls of African ancestry (AA), leveraging the differential linkage-disequilibrium of these populations in order to better localize the putative functional variants. We observed directionally consistent effects for all seven index SNPs in both populations, with four significantly associated (p<0.05) with risk in EAs (3p22.1, 7p15.3, 17p11.2, 22q13.1), and two significantly associated with risk in AAs (7p15.3 and 22q13.1). In a fixed effects meta-analysis of six regions (excluding the HLA region on chromosome 6), variation in five of the regions (2p33.3, 3p22.1, 7p15.3, 17p11.2, 22q13.1) had statistically significant associations with risk (Table 1). In one region, the index variant had the strongest association [rs4487645 at 7p15.3, (OR = 1.30, p = 8.7×10−8)]. Five of the six most significantly associated variants identified in the multiethnic analyses overlapped with biologically relevant features indicating regulatory activity based on CD20+ (B lymphocyte) cells, showing evidence of potential function; those included a missense variant in (17p11.2, rs34562254, Pro251Leu) in TNFRSF13B, which encodes a lymphocyte-specific protein in the tumor necrosis factor receptor family that interacts with the NF-kb pathway. Our study shows that these regions are important in MM risk across ethnicities and further supports the use of multiple ethnic groups in genetic studies to enhance identification of risk variants. Table 1.Most significant associations for each region in the multiethnic meta-analysis.Individuals of European AncestryIndividuals of African AncestryMulitethnic Metar2 with IndexcCHRSNPRAaFreqbORP-valueFreqbORP-valueORP-valueP-het2rs732075G0.591.222.0×10−30.621.122.0×10−21.162.6×10−40.280.09/0.283rs73069394A0.191.243.0×10−30.621.181.5×10−21.201.3×10−50.550.77/0.963rs12637184G0.761.136.0×10−20.921.192.7×10−11.151.0×10−20.640.94/1.007rs4487645C0.671.237.0×10−40.891.485.5×10−51.308.7×10−80.07-d17rs34562254A0.121.452.4×10−50.131.212.2×10−31.312.5×10−60.120.33/0.9022rs139400T0.491.224.0×10−40.531.172.1×10−31.191.2×10−60.630.63/0.96aRisk allelebFrequency of the risk allele in European and African ancestry studiescr2 metrics based on 1000 Genomes Project AFR/EUR populationsdIndex SNP Citation Format: Kristin A. Rand, Chi Song, Eric Dean, Daniel Serie, Karen Curtin, Dennis Hazelett, Amie E. Hwang, Xin Sheng, Alex Stram, David J. Van Den Berg, Carol Ann Huff, Leon Bernal-Mizrachi, Michael H. Tomasson, Sikander Ailawadhi, Anneclaire De Roos, Seema Singhal, Karen Pawlish, Edward Peters, Catherine Bock, David V. Conti, Graham Colditz, Todd Zimmerman, Scott Huntsman, John Graff, African Ancestry Prostate Cancer GWAS Consortium,African Ancestry Breast Cancer GWAS Consortium, Stephen J. Chanock, Michael Lieber, Jayesh Mehta, Eric A. Klein, Nalini Janakiraman, Richard K. Severson, Angela R. Brooks-Wilson, Vincent Rajkumar, Elizabeth E. Brown, Laurence Kolonel, Susan Slager, Brian E. Henderson, Graham G. Giles, John J. Spinelli, Brian Chiu, Kenneth C. Anderson, Jeffrey Zonder, Robert Z. Orlowski, Sagar Lonial, Nicola Camp, Celine Vachon, Elad Ziv, Dan O. Stram, Christopher A. Haiman, Wendy Cozen. Multiple myeloma susceptibility loci examined in African and European ancestry populations. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4629. doi:10.1158/1538-7445.AM2015-4629

Abstract Background: Genome wide association studies (GWAS) have recently identified 8 loci identified with multiple myeloma (MM) susceptibility. We hypothesized that germ line genetic variants may also affect MM progression. We performed the first GWAS of MM survival, by conducting a meta-analysis of two existing studies at University of California San Francisco (UCSF) and the Mayo Clinic, and replicating findings in the International Multiple Myeloma rESEarch (IMMEnSE) consortium. Methods: The UCSF study included 353 MM patients genotyped on an Illumina Omni5 and Illumina 660 array. The Mayo clinic study included 239 MM patients genotyped on an Affymetrix 6.0 array. We imputed missing genotypes using the 1000 Genomes dataset and performed a GWAS for survival using proportional hazards models adjusting for age and gender and genetic ancestry using principal components analysis. We conducted a meta-analysis of the results from the two GWAS. We replicated the top SNPs in the IMMEnSE cohort which included 772 patients with survival data from 7 European countries and one North American site. We used the dataset from Grundberg et al (Nat Genetics 2012) to test for an association between SNPs associated with survival and gene expression in immortalized lymphocytes. We used gene expression data from Zhan et al (Blood 2006) to analyze the association between expression and survival. Results: In meta-analysis, we found a genome wide significant association between SNPs on a region at 16p13 and survival (rs72773978; p=2.9x10e-10). Patients with the minor allele were at increased risk for mortality (HR 2.73; 95% CI: 1.99 - 3.73). The top associated SNPs were all in strong linkage disequilibrium and were centered over the gene FOPNL. We replicated the association in the IMMEnSE cohort and found a significant association between the top SNP from the discovery dataset (rs72773978 p=0.037). We used previously published data to determine the association between gene expression and the top associated SNPs from the GWAS. The top SNP (rs72773978) was not in previously published eQTL databases, so we used rs7201759 which was also associated with survival in our data (p=3.8x10e-10) as a proxy. The minor allele for this SNP was associated with increased levels of expression of FOPNL (p=9x10e-5), but not with expression of any other genes within a 1 megabase region. Finally, we also observed an association between increased expression of FOPNL and shorter survival (HR: 2.62 per standard deviation.; 95% CI: 1.08 - 6.32). Conclusion: Germ line variants at 16p13 near the FOPNL gene are associated with survival among MM patients. FOPNL is thought to be involved in centrosomal function, and centrosomal activity has previously been associated with survival in MM patients. Our results support this previous observation. Furthermore, they suggest that germ line variants could also be useful in assessing prognosis for patients with MM. Citation Format: Elad Ziv, Eric Dean, Donglei Hu, Alessandro Martino, Daniel Serie, Daniele Campa, Blake Aftab, Paige Bracci, Gabriele Buda, Jennifer Caswell, Charles Dumontet, Marek Dudziński, Laura Fejerman, Alexandra Greenberg, Scott Hunstman, Artur Jurczyszyn, Krzysztof Jamroziak, Shaji Kumar, Herlander Marques, Thomas Martin, Joaquin Martinez-Lopez, Vincent Rajkumar, Juan Sainz, Annette Juul Vangsted, Marzena Watek, Jeffrey Wolf, Susan Slager, Federico Canzian, Celine Vachon. Genome wide association study identifies variants at 16p13 associated with survival in multiple myeloma patients. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 5078. doi:10.1158/1538-7445.AM2014-5078

BACKGROUND The majority of patients diagnosed with clear cell renal cell carcinoma (ccRCC) have low‐risk disease with a < 10% chance of ccRCC‐specific death. DNA sequencing revealed that mutations in BAP1 (BRCA1 associated protein‐1) occur in 5% to 15% of ccRCC cases and are associated with poor outcomes. The vast majority of BAP1 mutations abolish protein expression. In this study, we used a highly sensitive and specific immunohistochemistry (IHC) assay to test whether BAP1 expression is an independent marker of ccRCC‐specific survival, particularly in patients with low‐risk disease. METHODS BAP1 expression was assessed, using IHC, in 1479 patients who underwent nephrectomy to treat clinically localized ccRCC. A centralized pathologist dichotomized patients as either BAP1‐positive or BAP1‐negative. The authors employed Kaplan‐Meier and Cox regression models to associate BAP1 expression with cancer‐specific survival. RESULTS A total of 10.5% of tumors were BAP1‐negative, 84.8% of tumors were BAP1‐positive, and 4.6% of tumors had ambiguous staining for BAP1. Patients with BAP1‐negative tumors have an increased risk of ccRCC‐related death (hazard ratio [HR] = 3.06; 95% confidence interval [CI] = 2.28‐4.10; P = 6.77 × 10 −14 ). BAP1 expression remained an independent marker of prognosis after adjusting for the UCLA integrated staging system (UISS) (HR = 1.67; 95% CI = 1.24‐2.25; P < .001). Finally, BAP1 was an independent prognostic marker in low‐risk patients with a Mayo Clinic stage, size, grade, and necrosis (SSIGN) score of ≤ 3 (HR = 3.24; 95% CI = 1.26‐8.33; P = .015). CONCLUSIONS This study used a large patient cohort to demonstrate that BAP1 expression is an independent marker of prognosis in patients with low‐risk (SSIGN≤ 3) ccRCC. Cancer 2014;1059–1067 . © 2014 American Cancer Society .

Context: The BRAF V600E mutation (BRAF-MUT) confers an aggressive phenotype in papillary thyroid carcinoma, but unidentified additional genomic abnormalities may be required for full phenotypic expression. Objective: RNA sequencing (RNA-Seq) was performed to identify genes differentially expressed between BRAF-MUT and BRAF wild-type (BRAF-WT) tumors and to correlate changes to patient clinical status. Design: BRAF-MUT and BRAF-WT tumors were identified in patients with T1N0 and T2–3N1 tumors evaluated in a referral medical center. Gene expression levels were determined (RNA-Seq) and fusion transcripts were detected. Multiplexed capture/detection and digital counting of mRNA transcripts (nCounter, NanoString Technologies) validated RNA-Seq data for immune system-related genes. Patients: BRAF-MUT patients included nine women, three men; nine were TNM stage I and three were stage III. Three (25%) had tumor infiltrating lymphocytes. BRAF-WT included five women, three men; all were stage I, and five (62.5%) had tumor infiltrating lymphocytes. Results: RNA-Seq identified 560 of 13 085 genes differentially expressed between BRAF-MUT and BRAF-WT tumors. Approximately 10% of these genes were related to MetaCore immune function pathways; 51 were underexpressed in BRAF-MUT tumors, whereas 4 (HLAG, CXCL14, TIMP1, IL1RAP) were overexpressed. The four most differentially overexpressed immune genes in BRAF-WT tumors (IL1B; CCL19; CCL21; CXCR4) correlated with lymphocyte infiltration. nCounter confirmed the RNA-Seq expression level data. Eleven different high-confidence fusion transcripts were detected (four interchromosomal; seven intrachromosomal) in 13 of 20 tumors. All in-frame fusions were validated by RT-PCR. Conclusion: BRAF-MUT papillary thyroid cancers have reduced expression of immune/inflammatory response genes compared with BRAF-WT tumors and correlate with lymphocyte infiltration. In contrast, HLA-G and CXCL14 are overexpressed in BRAF-MUT tumors. Sixty-five percent of tumors had between one and three fusion transcripts. Functional studies will be required to determine the potential role of these newly identified genomic abnormalities in contributing to the aggressiveness of BRAF-MUT and BRAF-WT tumors.

AbstractBackgroundGenetic variants at the CLU, CR1, and PICALM loci associate with risk for late‐onset Alzheimer's disease (LOAD) in genomewide association studies. In this study, our aim was to determine whether the LOAD risk variants at these three loci influence memory endophenotypes in black and white subjects.MethodsWe pursued an association study between single nucleotide polymorphism genotypes at the CLU, CR1, and PICALM loci and memory endophenotypes. We assessed black subjects (AA series: 44 with LOAD and 224 control subjects) recruited at Mayo Clinic Florida and whites recruited at Mayo Clinic Minnesota (RS series: 372 with LOAD and 1690 control subjects) and Florida (JS series: 60 with LOAD and 529 control subjects). Single nucleotide polymorphisms at the LOAD risk loci CLU (rs11136000), CR1 (rs6656401, rs3818361), and PICALM (rs3851179) were genotyped and tested for association with Logical Memory immediate recall, Logical Memory delayed recall, Logical Memory percent retention, Visual Reproduction immediate recall, Visual Reproduction delayed recall, and Visual Reproduction percent retention scores from the Wechsler Memory Scale–Revised using multivariable linear regression analysis, adjusting for age at exam, sex, education, and apolipoprotein E ε4 dosage.ResultsWe identified nominally significant or suggestive associations between the LOAD‐risky CR1 variants and worse Logical Memory immediate recall scores in blacks (P = .068–.046, β = −2.7 to −1.2). The LOAD‐protective CLU variant is associated with better logical memory endophenotypes in white subjects (P = .099–.027, β = 0.31–0.93). The CR1 associations persisted when the control subjects from the AA series were assessed separately. The CLU associations appeared to be driven by one of the white series (RS) and were also observed when the control subset from RS was analyzed.ConclusionThese results suggest for the first time that LOAD risk variants at CR1 may influence memory endophenotypes in blacks. In addition, the CLU LOAD‐protective variant may confer enhanced memory in whites. Although these results would not remain significant after stringent corrections for multiple testing, they need to be considered in the context of the LOAD associations with which they have biological consistency. They also provide estimates for effect sizes on memory endophenotypes that could guide future studies. The detection of memory effects for these variants in clinically normal subjects, implies that these LOAD risk loci might modify memory prior to clinical diagnosis of AD.