Sicotte, Hugues

AbstractBackgroundGlioblastoma (GBM) exhibits profound intratumoral genetic heterogeneity. Each tumor comprises multiple genetically distinct clonal populations with different therapeutic sensitivities. This has implications for targeted therapy and genetically informed paradigms. Contrast-enhanced (CE)-MRI and conventional sampling techniques have failed to resolve this heterogeneity, particularly for nonenhancing tumor populations. This study explores the feasibility of using multiparametric MRI and texture analysis to characterize regional genetic heterogeneity throughout MRI-enhancing and nonenhancing tumor segments.MethodsWe collected multiple image-guided biopsies from primary GBM patients throughout regions of enhancement (ENH) and nonenhancing parenchyma (so called brain-around-tumor, [BAT]). For each biopsy, we analyzed DNA copy number variants for core GBM driver genes reported by The Cancer Genome Atlas. We co-registered biopsy locations with MRI and texture maps to correlate regional genetic status with spatially matched imaging measurements. We also built multivariate predictive decision-tree models for each GBM driver gene and validated accuracies using leave-one-out-cross-validation (LOOCV).ResultsWe collected 48 biopsies (13 tumors) and identified significant imaging correlations (univariate analysis) for 6 driver genes: EGFR, PDGFRA, PTEN, CDKN2A, RB1, and TP53. Predictive model accuracies (on LOOCV) varied by driver gene of interest. Highest accuracies were observed for PDGFRA (77.1%), EGFR (75%), CDKN2A (87.5%), and RB1 (87.5%), while lowest accuracy was observed in TP53 (37.5%). Models for 4 driver genes (EGFR, RB1, CDKN2A, and PTEN) showed higher accuracy in BAT samples (n = 16) compared with those from ENH segments (n = 32).ConclusionMRI and texture analysis can help characterize regional genetic heterogeneity, which offers potential diagnostic value under the paradigm of individualized oncology.

Inflammation plays a pivotal role in peripheral artery disease (PAD). Cellular adhesion proteins mediate the interaction of leukocytes with endothelial cells during inflammation. To determine the association of cellular adhesion molecules with ankle–brachial index (ABI) and ABI category (≤1.0 vs >1.0) in a diverse population, 15 adhesion proteins were measured in the Multi-Ethnic Study of Atherosclerosis (MESA). To assess multivariable associations of each protein with ABI and ABI category, linear and logistic regression was used, respectively. Among 2364 participants, 23 presented with poorly compressible arteries (ABI > 1.4) and were excluded and 261 had ABI ≤ 1.0. Adjusting for traditional risk factors, elevated levels of soluble P-selectin, hepatocyte growth factor, and secretory leukocyte protease inhibitor were associated with lower ABI ( P = .0004, .001, and .002, respectively). Per each standard deviation of protein, we found 26%, 20%, and 19% greater odds of lower ABI category ( P = .001, .01, and .02, respectively). Further investigation into the adhesion pathway may shed new light on biological mechanisms implicated in PAD.

ObjectiveAt the cellular level, how excess adiposity promotes atherogenesis is not fully understood. One pathway involves secretion of adipokines that stimulate endothelial dysfunction through increased expression of adhesion molecules. However, the relationship of adiposity to adhesion molecules that promote atherosclerosis is largely unknown.MethodsLinear regression models were used to assess the sex‐specific associations of soluble cellular adhesion molecules (sP‐ and sL‐selectin, sICAM‐1, sVCAM‐1, and sHGF) and adiposity in 5,974 adults examined as part of the Multi‐Ethnic Study of Atherosclerosis (MESA). Adiposity measures included body mass index (BMI), waist‐to‐hip‐ratio (WHR), and computed tomography measures of subcutaneous adipose tissue (SAT) and visceral adipose tissue (VAT).ResultsThe mean age was 64 years and 52% were female. In multivariable models adjusting for traditional cardiovascular risk factors, sHGF was positively associated with BMI, WHR, and VAT in both males and females, and sP‐selectin with WHR and VAT in males. sVCAM‐1 was inversely associated with VAT in females only.ConclusionsOur results showed the relation of adiposity to soluble cellular adhesion proteins was similar across adiposity measures and for both sexes. However, the relationship between adiposity and sVCAM‐1 and P‐selectin may be modified by sex and the measure used to assess adiposity.

SummaryHepatocyte growth factor (HGF) is a mesenchyme‐derived pleiotropic factor that regulates cell growth, motility, mitogenesis, and morphogenesis in a variety of cells, and increased serum levels of HGF have been linked to a number of clinical and subclinical cardiovascular disease phenotypes. However, little is currently known regarding which genetic factors influence HGF levels, despite evidence of substantial genetic contributions to HGF variation. Based upon ethnicity‐stratified single‐variant association analysis and trans‐ethnic meta‐analysis of 6201 participants of the Multi‐Ethnic Study of Atherosclerosis (MESA), we discovered five statistically significant common and low‐frequency variants: HGF missense polymorphism rs5745687 (p.E299K) as well as four variants (rs16844364, rs4690098, rs114303452, rs3748034) within or in proximity to HGFAC. We also identified two significant ethnicity‐specific gene‐level associations (A1BG in African Americans; FASN in Chinese Americans) based upon low‐frequency/rare variants, while meta‐analysis of gene‐level results identified a significant association for HGFAC. However, identified single‐variant associations explained modest proportions of the total trait variation and were not significantly associated with coronary artery calcium or coronary heart disease. Our findings indicate that genetic factors influencing circulating HGF levels may be complex and ethnically diverse.

Abstract Introduction: Triple negative breast cancers (TNBC) are characterized as lacking estrogen receptor (ER), progesterone receptor (PR) and human epidermal growth factor receptor 2 (HER2) expression and TNBC patients have higher rates of recurrence and death compared with other breast cancer subtypes. For TNBC patients who fail standard chemotherapy, there are a lack of novel drug therapies, given the absence of well-defined molecular targets. Recently, a microarray meta-analysis identified 7 triple negative subtypes, including the validation of the luminal androgen receptor (LAR) positive subtype [Lehmann, 2011]. However, microarray technology is dependent on probe-target specificity and the 7 subtypes have yet to be validated using RNA sequencing data, and the presences of recurrent genomic alterations in the 7 subtypes are unknown. METHODS: We obtained 1106 breast cancer RNA-Seq bam files from The Cancer Genome Atlas (TCGA) and aligned with Tophat v1.3. The PAM50 intrinsic gene signature was used to extract a cohort of 128 TNBC samples. Consensus clustering of genes, greater than 75th percentile of variance, was performed using Kmeans clustering in Spearman's correlation space. A nearest centroid prediction model was developed from genes differentially expressed among the clusters. Eighty independent TNBC RNA sequencing samples were obtained (British Colombia; BC) [Shah, 2012] which were calibrated to our TNBC conditional quantile normalized cohort and sub-typed by our model. RESULTS: Using RNA-Seq gene expression count data, we identified 5 clusters, all of which were stable, including the LAR cluster. Signaling pathway impact analysis (SPIA) implicated cytokine-cytokine receptor interaction, leukocyte transendothelial migration, and regulation of actin cytoskeleton pathways commonly altered in the non-LAR TNBC subtypes. In contrast, cell cycle, ECM-receptor interaction, endocrine regulated calcium reabsorption, and insulin signaling pathways were altered in the LAR versus non-LAR subtypes. Neuroactive ligand-receptor interactions were observed to be altered commonly between all sub-types. We then applied our model to the Shah, et al cohort. In this cohort, the LAR subtype was consistent with Shah's classification of ‘other’ TNBC and contained no basal samples by PAM50 intrinsic modeling. Analysis of sub-type specific mutation data from the BC cohort demonstrates an increased mutational load in ECM-related proteins, particularly the myosins, along with increased TP53 clonality in the non-LAR subtypes. CONCLUSIONS: Using TCGA RNASeq data, we have confirmed the presence of 5 major TNBC subtypes, including the LAR; which was negligible in basal composition by PAM50 intrinsic modeling. SPIA pathway analysis indicates a core set of pathways demonstrating altered expression across the TNBC sub-types and the identification of molecular targets within each subtype is ongoing. Citation Format: Kevin J. Thompson, Xiaojia Tang, Zhifu Sun, Jason P. Sinnwell, Hugues Sicotte, Douglas W. Mahoney, Steven Hart, Peter T. Vedell, Poulami Barman, Jeanette E. Eckel Passow, Eric D. Wieben, James N. Ingle, Judy C. Boughey, Liewei Wang, Richard Weinshilboum, Krishna R. Kalari, Matthew P. Goetz. Molecular classification of triple negative breast cancer via RNA-sequencing data. [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 5592. doi:10.1158/1538-7445.AM2014-5592

Abstract INTRODUCTION Next generation sequencing (NGS) of patients has significantly changed our ways to study cancer genomics as it provides precise estimates of gene expression, fusion transcripts, expressed single nucleotide variants (eSNVs), splice variants and copy number variants. The Breast Cancer Genome Guided Therapy (BEAUTY) is an ongoing clinical study in which RNA sequencing (RNAseq) and whole exome sequencing (WES) are performed prior to, during and after neoadjuvant chemotherapy. Here we report the use of these sequencing technologies to investigate gene expression levels and mutational profiles in a triple negative breast cancer (TNBC) patient enrolled in BEAUTY whose disease did not respond to neoadjuvant paclitaxel and anthracycline/cycphosphamide. METHODS Computational approaches were used to integrate WES and RNAseq data obtained before therapy (V1T), after 12 weekly paclitaxel treatments (V2T) and after anthracycline-based regimen at surgical resection (V3T) to study a single patient with persistent TNBC disease after neoadjuvant chemotherapy. RESULTS Using RNA-Seq data, we identified an inter-chromosomal fusion transcript between chromosome 20 and 22 (GNAS-TTC38) that was highly expressed at V1T, V2T and V3T. We also identified intra-chromosomal fusion transcripts that were expressed at two time points, such as fusion transcript (KANSL1-ARL17A) on chromosome 17 for V1T and V2T and fusion transcript (RBM12B-LINC00535) on chromosome 8 for V2T and V3T time points. Several gene expression changes were also observed. Gene expression analysis of V1T, V2T and V3T tumors was performed. Differential temporal gene expression profiles of 9884 genes that were significant and varying at different time points were obtained for pathway analysis. Pathway analysis of 9884 genes identified up regulation and down regulation of several transcription factors with a fold change of 2x or more. When compared to blood, DNA tumor and RNA-Seq data, we identified 81 common somatic eSNVs that were expressed in both V1T and V2T time points and we are in the process of investigating V3T data. We found alterations of key transporter domains (CD225, Coatamer_beta_C, DUF2435, Dynamitin, EI24, GLTP, LMF1, Porin_3, V-ATPase_C) in our V1T and V2T SNV data. Similar to gene expression analysis, we are in the process of obtaining the list of mutations at various time points to identify driver and passenger mutation candidate genes for this specific TNBC patient. CONCLUSIONS Our initial time-series analysis of eSNV, fusion transcripts and gene expression data demonstrate that intensive analysis for individual patients is feasible. Further investigation of drug transporters and transcription regulators may help develop personalized treatment strategies for patients with disease resistant to current regimens. Citation Format: Krishna R. Kalari, Xiaojia Tang, Kevin J. Thompson, Douglas W. Mahoney, Poulami Barman, Jason P. Sinnwell, Hugues Sicotte, Peter Vedell, Steven N. Hart, Travis J. Dockter, Katie N. Jones, Amy L. Conners, Ann M. Moyer, Daniel W. Visscher, Jia Yu, Bowen Gao, Sarah A. McLaughlin, John A. Copland, Alvaro Moreno-Aspitia, Donald W. Northfelt, Richard J. Gray, Vera J. Suman, Jeanette E. Eckel Passow, Jean-Pierre A. Kocher, Eric D. Wieben, Gianrico Farrugia, Cloann G. Schultz, James N. Ingle, Richard Weinshilboum, Matthew P. Goetz, Liewei Wang, Judy C. Boughey. Analysis of sequencing data to identify potential drug targets for an individual newly diagnosed with basal breast cancer who failed to respond to current standard neoadjuvant chemotherapy. [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 4185. doi:10.1158/1538-7445.AM2014-4185

Abstract Background: Specific germline alterations within the TERT, EGFR, CCDC26, CDKN2A/B, PHLDB1, TP53 and RTEL1 regions are associated with development of glioma. While some of these germline variants are associated with all gliomas (e.g. rs2736100 in TERT), others are associated with a specific morphologic and/or molecular subtype. For example, the RTEL1 region variants are associated with primary glioblastoma and rs55705857 in CCDC26 is associated with 1p/19q co-deleted oligodendrogliomas and with IDH mutant astrocytic gliomas. We hypothesized that germline variants in these regions will be associated with other molecular subtypes. The Cancer Genome Atlas (TCGA) and other groups have described molecular subtypes of glioma based on acquired somatic mutation patterns. Herein, we grouped glioma patients into mutation-based molecular subgroups based on three molecular alterations and determined if germline polymorphisms are associated with these mutation-based molecular subtypes. Methods: Five molecular subtypes were defined based on combinations of the following molecular alterations: TERT promoter mutation (TERTmut vs. TERTwt), IDH1/2 mutation (IDH mut vs. IDH wt) and 1p/19q co-deletion (1p19qcodel vs. 1p19q noncodel). Germline single nucleotide polymorphisms (SNPs) were evaluated by custom Illumina array analysis of blood-derived DNA and 1K genome imputation. Results: One hundred ninety-nine gliomas (24 oligodendrogliomas, 59 mixed oligoastrocytomas, 41 grade 2-3 astrocytomas and 75 glioblastomas) had both SNP data and mutation data available in order to assign them into the five molecular subgroups. The prevalence of each of the molecular subtypes was as follows: TERTmut/IDHmut/1p19qcodel (20%), TERTmut/IDHmut/1p19qnoncodel (5%), TERTmut/IDHwt/1p19qnoncodel (32%), TERTwt/IDHmut/1p19qnondel (35%) and TERTwt/IDHwt/1p19qnoncodel (8%). The TERT SNP (rs2736100) was protective for TERTmut/IDHmut/1p19qnoncodel and TERTwt/IDHmut/1p19qnoncodel gliomas. The CCDC26 SNP (rs55705857) was strongly associated with the risk of developing TERTmut/IDHmut/1p19qcodel and TERTwt/IDHmut/1p19qnoncodel gliomas. The TP53 SNP (rs7837822) was strongly associated with risk of developing TERTmut/IDHmut/1p19qnoncodel and TERTwt/IDHmut/1p19qnoncodel gliomas. Interestingly, while one RTEL1 region SNP (rs6062297) was strongly associated with glioma risk, RTEL1 SNPs rs6010620 and rs2297440 were protective for the development of TERTmut/IDHwt/1p19qnoncodel gliomas. We are currently in the process of validating these results using the TCGA glioblastoma and low grade glioma data. Conclusions: Glioma patients can be meaningfully classified according to their acquired mutation subtype. Importantly, our results suggest that there are significant associations between germline polymorphisms and mutation-based glioma molecular subtypes. Citation Format: Jeanette E. Eckel-Passow, Thomas M. Kollmeyer, Gobinda Sarkar, Anisha Chada, Paul A. Decker, Matthew L. Kosel, Alissa A. Caron, Hugues Sicotte, Kannabiran Nandakumar, Naresh Prodduturi, Brian P. O'Neill, Daniel H. Lachance, Robert B. Jenkins. The association of glioma germline risk SNPs with mutation-based molecular subgroups. [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 4714. doi:10.1158/1538-7445.AM2014-4714

Abstract Recent studies have shown that associations between common genetic variation and gene expression in breast tumors provide insight into the functional mechanisms underlying breast cancer risk loci. However, the relationship between genetic variation, gene expression, and risk of triple negative (TN) breast cancer remains largely unexplored. We performed a genome-wide expression quantitative trait locus (eQTL) analysis of TN breast cancer using 668 formalin-fixed paraffin embedded TN tumors from the Triple Negative Breast Cancer Consortium (TNBCC) to identify novel genes relevant to TN breast cancer risk and to further explore the biology underlying known TN risk loci. Cis-eQTLs were defined as correlations between single nucleotide polymorphisms (SNP) and expression of genes within 1Mb, and trans-eQTLs were defined as the remainder of all SNP-gene correlations. In a genome-wide analysis, we identified 68,012 cis-eQTLs, representing correlations between 42,193 SNPs and 2,092 genes that were significant at a 10% false discovery rate (FDR). We also identified 70 trans-eQTLs across two unique genes at the same FDR threshold. Five cis-eQTLs involved SNPs previously associated with TN breast cancer risk (p<5.0x10-5) in 3,700 TN cases and 4,700 controls from the TNBCC. Three of these, located on chromosomes 5p (cis-eQTL p=4.8x10-34; risk p=1.1x10-7), 7p (cis-eQTL p=6.1x10-10; risk p=1.3x10-7), and 17q (cis-eQTL p=1.8x10-15; risk p=2.9x10-7), have previously been reported in public databases of normal tissue. The remaining two TN breast cancer cis-eQTLs, located on chromosomes 8q (cis-eQTL p=6.3x10-6, risk p=8.3x10-6) and 14q (cis-eQTL p=3.1x10-6, risk p=1.2x10-5), have not been reported as eQTLs in other tissues. In addition, 12 of 25 known TN breast cancer risk loci (PEX14, MDM4, 2q31.1, ESR1, 11q13.1, 11q24.3, 12p13.1, PTHLH, NTN4, 12q24, 19p13.1, MLK1), contained at least one significant cis-eQTL. In summary, we have identified five cis-eQTLs that may identify novel risk loci for TN breast cancer. Citation Format: Kristen S. Purrington, Drakoulis Yannoukakos, Jane Carpenter, Heli Nevanlinna, Angela Cox, Gianluca Severi, Christine Ambrosone, Amanda Ewart Toland, Andrew K. Godwin, Hiltrud Brauch, Peter A. Fasching, Penelope Miron, Jenny Chang-Claude, Nicholas G. Martin, Grant W. Montgomery, Vessela Kristensen, Hoda Anton-Culver, Paul Goodfellow, Janet E. Olson, Hugues Sicotte, Naresh Prodduturi, Daniel W. Visscher, Jeanette E. Eckel-Passow, S. Keith Anderson, Seth Slettedahl, Curtis Olswold, Xianshu Wang, V. Shane Pankratz, Susan Slager, Wei Zheng, Arto Mannermaa, Ute Hamann, Diana M. Eccles, Celine M. Vachon, Fergus J. Couch. Expression quantitative trait locus analysis of triple negative breast cancer. [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 3266. doi:10.1158/1538-7445.AM2014-3266

Abstract INTRODUCTION Patient derived xenografts (PDX) may better reflect individual patient (pt) tumor biology; however, the feasibility of collecting PDX from percutaneous tumor biopsies (PTB) in the neoadjuvant setting is unknown. Furthermore, drug response phenotypes observed in PDX have not been prospectively compared to the corresponding pt clinical outcomes. METHODS The Breast Cancer Genome Guided Therapy Study (BEAUTY) is a prospective Mayo study of pts with high-risk breast cancer treated with neoadjuvant weekly paclitaxel (T) +/- trastuzumab followed by anthracycline based chemotherapy. PTB (at baseline) and residual surgical tissue (after all chemotherapy) are obtained for next generation sequencing (NGS) and PDX. Tumor biopsies (1-2 cores from 14 gauge needle) were implanted with Matrigel <1 hour of collection in the flanks of NOD-SCID or NSG mice. Low dose estradiol was supplemented in the drinking water. Primary outgrowth rate was defined as PDX tumor volume >50 mm3. Take rate was defined as development of at least 1 stably transplantable xenograft line/pt. To determine whether clinical T response assessed by MRI corresponded with in vivo T response, pretreatment PDX from 5 pts were injected into NOD-SCID mice (20 mice per pt PDX) and when tumors reached 100-200mm3, mice were randomized to no treatment vs T (20 mg/kg, ip. every 3-4 days). Two of these 5 patients had a MRI response defined as >30% decrease in longest lesion. RESULTS Pretreatment PTB from 81 unique pts were implanted in 251 mice (2-4 mice/pt). PDX outgrowth rates were 33.3% (27/81 pts) and 22 stable PDX were established (overall take rate 27.2%). Take rates were as follows: triple negative breast cancer (46%; 13/28); HER2 (27%; 6/22), Luminal B (13%; 3/22), and luminal A (0%; 0/9). Residual surgical tumor (after all treatment) from 17 pts was injected into 85 mice (average 5 mice/pt) and the initial outgrowth rate was 23% (4/17) with 3 stably transplantable lines established. PDX, derived from pretreatment PTB of 5 pts (2 responders and 3 non-responders), were assessed for in vivo T response. The size of the T treated group was significantly smaller than the no treatment group for the PDX derived from the 2 clinical responders, with complete disappearance of tumor by 18 days. In contrast, the PDX derived from the 3 clinical non-responders had no evidence for T response. CONCLUSIONS We have demonstrated the feasibility of using PTB to establish PDX in a prospective neoadjuvant clinical study and have demonstrated similar T drug response phenotypes in in the PDX as seen in the corresponding pt. These data suggest that PDX generated prospectively may be useful for biomarker validation and the development and individualization of new drug therapy. Funded by the Mayo Clinic Center for Individualized Medicine and the Mayo Clinic Cancer Center Citation Format: Jia Yu, Ping Yin, Bowen Gao, Jason P. Sinnwell, Ann M. Moyer, Daniel W. Visscher, Amy L. Conners, Travis J. Dockter, Krishna R. Kalari, Xiaojia Tang, Kevin J. Thompson, Hugues Sicotte, Douglas W. Mahoney, Steven N. Hart, Peter T. Vedell, Poulami Barman, Katie N. Jones, Sarah A. McLaughlin, John A. Copland, Alvaro Moreno Aspitia, Donald W. Northfelt, Richard J. Gray, Vera J. Suman, Jeanette E. Eckel Passow, Eric D. Wieben, James N. Ingle, Zhenkun Lou, Gianrico Farrugia, Richard Weinshilboum, Matthew P. Goetz, Judy C. Boughey, Liewei Wang. Feasibility of using percutaneous tumor biopsies from a prospective neoadjuvant breast cancer study to develop patient derived xenografts and assess in vivo chemotherapy sensitivity. [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 1195. doi:10.1158/1538-7445.AM2014-1195