Abstract Glioblastoma (GBM), also known as grade IV astrocytoma, is the most common and devastating primary cancer of the central nervous system in adults. Despite the most advanced therapeutic strategies combining surgery, chemotherapy, and radiation, survival expectancy of GBM patients averages 14 months, with a recurrence rate of almost 100%. Remarkably, tumor location has great implications in the prognosis of GBM. Approximately 50-60% of GBM tumors infiltrates or contacts the subventricular zone (SVZ). These patients have significantly worse outcomes, in terms of median overall survival, time to progression and reoccurrence. The reason for the worse prognosis of SVZ-infiltrating tumors is unknown. The SVZ is the largest neurogenic niche in the adult brain, here the differentiation and migration of SVZ cells is regulated by direct contact with the cerebrospinal fluid (CSF). We suggest that CSF contribute to enhance GBM malignant characteristics in SVZ-infiltrating tumors. This study evaluated the impact of CSF exposure on GBM cells proliferation, migration in vitro and in vivo, as well as changes in gene expression profile. We demonstrated that CSF affects the malignancy of GBM cells by increasing their ability to proliferate and migrate. GBM cells stimulated with cancer CSF, non-cancer CSF, or control conditions at different time points (24-72 hrs) were evaluated for viability (by Alamar Blue assay) and proliferation rate (Ki67+ immunostaining), as well as for migratory capabilities (by means of transwell assay and time-lapse microscopy). We recapitulated these CSF effects in vivo using an in vivo model in which both female and male nude mice were implanted with GBM cells in suspension with human cancer or non-cancer CSF in an hydrogel as vehicle (n=8 mice per group). Tumor growth was followed via IVIS imaging and then studied by measuring tumor area on H&E slides and immunohistochemistry staining for Ki67, 21 days after tumor implantation. Animals receiving cancer CSF presented with greater tumor size and Ki67 proliferation index. To further dissect the mechanisms behind these observations, we performed a transcriptome analysis of CSF-treated GBM cells. At the molecular level we observed that cancer CSF induced changes in pathways regulating apoptosis, survival, integrins signaling, angiogenesis, glucose metabolism and response to inflammatory stimuli, all of which might be responsible for the increase in GBM malignancy. Taken together, this study provides direct evidence that CSF is an important player in determining tumor growth and invasion through the activation of complex gene expression patterns characteristic of a malignant phenotype. Understanding the interaction of brain tumors and CSF may lead to new therapeutic strategies that will target the CSF components in GBM patients. Citation Format: Anna Carrano, Jordan Phillipps, Montserrat Lara-Velazquez, Natanael Zarco, Paola Suarez-Meade, Kaisorn Chaichana, Alfredo Quinones-Hinojosa, Yan Asmann, Hugo Guerrero-Cázares. The role of human cerebrospinal fluid in glioblastoma malignancy [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 1574.
Impact of Variant-level Batch Effects on Identification of Genetic Risk Factors in Large Sequencing Studies
No container title available | August 13, 2020 | Volume: No volume available | Issue: No issue available | Pages: No page available
Abstract Background: Genetic studies have shifted to sequencing-based rare variants discovery after decades of success in identifying common disease variants by Genome-Wide Association Studies using Single Nucleotide Polymorphism chips. Sequencing-based studies require large sample sizes for statistical power but often inadvertently introduce batch effects because samples are typically collected, processed, and sequenced at multiple centers. Conventionally, batch effects are first detected and visualized using Principal Components Analysis and then controlled by including batch covariates in the disease association models. For sequencing-based genetic studies, because all variants included in the association analyses have passed quality control measures, this conventional approach treats every variant as equal and ignores the substantial differences still remaining in variant qualities and characteristics such as genotype quality scores, alternative allele fractions (fraction of reads supporting alternative allele at a variant position) and sequencing depths. Results: In the Alzheimer’s Disease Sequencing Project (ADSP) exome dataset of 9,904 cases and controls, we discovered hidden variant-level differences between sample batches of three sequencing centers and two exome capture kits. Although sequencing centers were included as a covariate in our association models, we observed differences at the variant level in genotype quality and alternative allele fraction between samples processed by different exome capture kits that significantly impacted both the confidence of variant detection and the identification of disease-associated variants. Furthermore, we found that the association signals of a subset of top disease risk variants came exclusively from samples processed by one exome capture kit that was more effective at capturing the alternative alleles compared to the other kit. Conclusions: Our findings highlight the importance of additional variant-level quality control for large sequencing-based genetic studies. More importantly, we demonstrate that automatically filtering out variants with batch differences may lead to false negatives if the batch discordance came largely from quality differences and if the variants from one batch had better quality scores.
Immune Responses and Risk of Triple-negative Breast Cancer: Implications for Higher Rates among African American Women
Cancer Prevention Research | August 5, 2020 | Volume: 13 | Issue: 11 | Pages: 901-910
AbstractThe etiology of triple-negative breast cancers (TNBC) is poorly understood. As many TNBCs develop prior to the initiation of breast cancer screening or at younger ages when the sensitivity of mammography is comparatively low, understanding the etiology of TNBCs is critical for discovering novel prevention approaches for these tumors. Furthermore, the higher incidence rate of estrogen receptor–negative breast cancers, and specifically, of TNBCs, among young African American women (AAW) versus white women is a source of racial disparities in breast cancer mortality. Whereas immune responses to TNBCs have received considerable attention in relation to prognosis and treatment, the concept that dysregulated immune responses may predispose to the development of TNBCs has received limited attention. We present evidence that dysregulated immune responses are critical in the pathogenesis of TNBCs, based on the molecular biology of the cancers and the mechanisms proposed to mediate TNBC risk factors. Furthermore, proposed risk factors for TNBC, especially childbearing without breastfeeding, high parity, and obesity, are more prevalent among AAW than white women. Limited data suggest genetic differences in immune responses by race, which favor a stronger Thr type 2 (Th2) immune response among AAW than white women. Th2 responses contribute to wound-healing processes, which are implicated in the pathogenesis of TNBCs. Accordingly, we review data on the link between immune responses and TNBC risk and consider whether the prevalence of risk factors that result in dysregulated immunity is higher among AAW than white women.
Loss of TMEM106B leads to myelination deficits: implications for frontotemporal dementia treatment strategies
Abstract Genetic variants that define two distinct haplotypes at the TMEM106B locus have been implicated in multiple neurodegenerative diseases and in healthy brain ageing. In frontotemporal dementia (FTD), the high expressing TMEM106B risk haplotype was shown to increase susceptibility for FTD with TDP-43 inclusions (FTD-TDP) and to modify disease penetrance in progranulin mutation carriers (FTD-GRN). To elucidate the biological function of TMEM106B and determine whether lowering TMEM106B may be a viable therapeutic strategy, we performed brain transcriptomic analyses in 8-month-old animals from our recently developed Tmem106b−/− mouse model. We included 10 Tmem106b+/+ (wild-type), 10 Tmem106b+/− and 10 Tmem106−/− mice. The most differentially expressed genes (153 downregulated and 60 upregulated) were identified between Tmem106b−/− and wild-type animals, with an enrichment for genes implicated in myelination-related cellular processes including axon ensheathment and oligodendrocyte differentiation. Co-expression analysis also revealed that the most downregulated group of correlated genes was enriched for myelination-related processes. We further detected a significant loss of OLIG2-positive cells in the corpus callosum of Tmem106b−/− mice, which was present already in young animals (21 days) and persisted until old age (23 months), without worsening. Quantitative polymerase chain reaction revealed a reduction of differentiated but not undifferentiated oligodendrocytes cellular markers. While no obvious changes in myelin were observed at the ultrastructure levels in unchallenged animals, treatment with cuprizone revealed that Tmem106b−/− mice are more susceptible to cuprizone-induced demyelination and have a reduced capacity to remyelinate, a finding which we were able to replicate in a newly generated Tmem106b CRISPR/cas9 knock-out mouse model. Finally, using a TMEM106B HeLa knock-out cell line and primary cultured oligodendrocytes, we determined that loss of TMEM106B leads to abnormalities in the distribution of lysosomes and PLP1. Together these findings reveal an important function for TMEM106B in myelination with possible consequences for therapeutic strategies aimed at lowering TMEM106B levels.
Novel Mutations in Types 2 & 3 von Willebrand Disease and Correlation with von Willebrand Factor Multimer Patterns.
Abstract Background & Objectives: Diagnosis and classification of VWD are currently based on integration of personal and family bleeding histories and results of protein-based diagnostic tests (VWF:Ag, VWF:RCo, FVIII:C, VWF:CB and VWF multimers) which have performance limitations. Genetic testing is emerging as a complementary diagnostic tool. We have identified mutations and correlated VWF multimer patterns in our patients (pt) with types 2 & 3 VWD, and report novel candidate mutations. Patients & Methods: Unrelated pt with type 3 (n=6) and subtypes 2 (n=22) VWD, from Mayo Comprehensive Hemophilia Center, consented to this IRB approved study. PCR amplification, from genomic DNA, of VWF gene (splice junctions, coding and promoter regions), avoiding pseudogene amplification, followed by ABI® sequencing and analysis (Mutation Surveyor: SoftGenetic®) were performed. Mutations were compared to VWD (ISTH), NCBI NR nucleotide and DV SNP databases. Comparison with available VWF sequences from other vertebrate species was performed. Selected regions not known to contain mutations remain to be sequenced. VWF multimer analysis was performed using a novel in-gel immunostaining and infrared imaging system. Results: Type 2A VWD (n=10), a novel mutation in exon 28 (E28) (V1524G: 2 pt) with characteristic VWF multimer abnormalities (decreased HMWM, increased satellite banding, suggesting enhanced VWF proteolysis). Type 2B VWD (n=6), 2 pt were compound heterozygous for 2 previously reported mutations, the first in E28 (R1306W)/E48 (T2647M) and the second in E28 (V1316M)/E20 (R854Q). Type 2M VWD (n=4): an 84 yr female with personal (spontaneous and post-surgical) and family bleeding histories, VWF:RCo 72%, VWF:Ag 118% (RCo:Ag ratio 0.61), normal FVIII and platelets, and aberrant VWF multimer banding pattern without substantive reduction of HMWM, had a novel mutation in E52 (S2775C; conserved in mouse, rat, chicken and dog); a 14 yr female with a history of bleeding, VWF:RCo 34%, VWF:Ag 68% (RCo:Ag ratio 0.5) and normal FVIII and platelets, was compound heterozygous for a known mutation (R1399C subtype not classified in VWD database) and two novel mutations: E30 (P1725S) (conserved in mouse, rat, dog, cow and chimp) & E49 (T2666M) (conserved in mouse, dog, rat not in chicken). VWF multimers demonstrated abnormal banding pattern without substantive reduction of HMWM; 2 unrelated pt and families with ultra-large and smeary multimers, previously classified as Vicenza variant, had R1374C mutations and were reclassified as Type 2M VWD. Type 2N VWD (n=1) previously reported mutation E20 (R854Q) was confirmed. Type 3 VWD, 2 pt had novel mutations, the first in E22 (E950X) and the second was a compound heterozygote for E43 (R2478Q)/E40 (E2322V) mutations. In general, patients with similar mutations had similar multimer patterns. Conclusions: We report novel candidate mutations in types 2 & 3 VWD and demonstrate novel VWF multimer patterns with our in-gel staining system. Multimer patterns correlated well with the underlying genotype. The novel mutations are likely causative, given their occurrence at highly conserved residues, but this needs to be confirmed with expression studies. Integration of genetic testing into the diagnostic workup of VWD will potentially lead to more accurate diagnosis and subtyping of VWD, and may further refine protein based testing and provide additional biological insight into VWD.
Skeletal muscle gene transcript changes in type 1 diabetic patients following insulin deprivation
The FASEB Journal | April 10, 2020 | Volume: 21 | Issue: 5 | Pages: No page available
Insulin resistance is associated with skeletal muscle mitochondrial (mito) dysfunction in many conditions including type 2 diabetes (T2D). Poor glycemic control in T2D is associated with muscle mito gene transcript alterations. It is unknown whether muscle mito changes are related to reduced insulin action or vice versa. We obtained vastus lateralis muscle biopsies from 7 type 1 diabetic patients (T1D) following 8–12 hrs of insulin deprivation (I−) and replacement (I+) and measured skeletal muscle mito ATP production (MATP) and gene transcripts. I− and I+ changed glucose and glucagon levels (glucose=291±15mg/dl in I− vs 97±3 in I+ and glucagon=79.8±15.3 pg/mL in I− vs 47.3±5.9 in I+, p<0.04). Whole body O 2 consumption increased during I− by 15% (p<0.05) but maximal MATP decreased by 21.3–27.9% (p<0.05). Gene transcript profiling in skeletal muscle analyzed using Affymetrix GeneChip arrays showed alterations in 2355 transcripts. Pathway analysis showed down‐regulation of oxidative phosphorylation and arachidonic acid metabolism, and up‐regulation of VEGF and actin cytoskeleton signaling in I−. We conclude that the reduced insulin action that occurred during I− in T1D resulted in substantial changes in gene transcript levels of muscle mito genes and reduced maximal MATP, supporting that insulin is key regulator of muscle mito function by exerting its effect at various levels. Supported by NIHRO1DK41973.
HLA class-I and class-II restricted neoantigen loads predict overall survival in breast cancer
OncoImmunology | April 1, 2020 | Volume: 9 | Issue: 1 | Pages: No page available
7021 Background: Case-control studies (CCS) have identified germline variants (GVs) & genes associated with AML risk and/or outcome. However, whether these genes also harbor somatic mutations (SMs) in AML and whether these SMs are associated with clinical phenotypes and outcome is unknown. We systemically interrogated PubMed (terms "AML genetic polymorphisms"), yielding 612 reports (1981-2018), from which we identified unique variants/single nucleotide polymorphisms (SNPs) in 100 genes associated with AML risk and/or outcome in CCS. To determine the prevalence and relevance of SMs in these genes in an unselected AML population, we performed an analysis using WES of cases in the highly-annotated MCAEC (Finn, Ca Epidemiol 2015). Methods: We performed WES using remnant cytogenetic pellets from AML diagnosis in 98 MCAEC patients with available samples, sequenced at a depth of ~100 million paired-end 100bp reads using Agilent SureSelectXT Human All Exon V5 + UTRs target enrichment kit. Sequencing reads were aligned to human reference genome, and SMs including non-synonymous (NS) and truncating single nucleotide variants (SNVs) and small INDELs were identified and filtered using Exome Sequencing Project, 1000 genome, HapMap, & Mayo Clinic (MC) internal Biobank GV database. Copy number aberrations (CNAs) were identified & filtered using public CN polymorphism databases. For this analysis we evaluated the prevalence of somatic SNVs, INDELs, & CNAs in the 100 nominated genetic risk genes in collaboration with MC Cancer Center to determine their association with clinical phenotype and outcome. Results: Of the 100 unique genes, there were 40 with SM in our dataset, with 126 total mutations. Most were mutated in a very small number of patients, allowing insufficient statistical power to perform association analyses. However, an NS mutation in CYP2B6 (c.785A > G, p.Lys262Arg) had a minor allele frequency of 12.0%; this variant has been reported in association with breast cancer risk, but not with AML. Conclusions: We studied SM’s in AML genetic risk genes using WES of tumor samples. While representing a smaller sample size, this suggests their contribution to AML development may be limited in a broader clinical population. These results & prevalence of the novel CYP2B6 rs2279343 variant in AML must be confirmed in larger cohorts.
The relationship between quantitative HER2 gene expression by the 21-gene RT-PCR assay and adjuvant trastuzumab (H) benefit in NCCTG (Alliance) N9831.
Journal of Clinical Oncology | March 11, 2020 | Volume: 31 | Issue: 15_suppl | Pages: 520-520
520 Background: There is considerable interest in developing HER2 testing criteria for adjuvant H. We used the 21-gene assay to examine the relationship of HER2 mRNA to benefit from H. Methods: N9831 compared adjuvant chemotherapy AC-T to concurrent chemotherapy-trastuzumab AC-TH in stage I-III HER2+ breast cancer. Recurrence Score (RS) and HER2 mRNA expression were determined by Oncotype DX (neg<10.7, equiv 10.7 to <11.5, and pos ≥11.5 log2 expression units). Cox regression was used to assess the association of HER2 expression with H benefit for distant recurrence. Results: Median follow-up: 7.4 yrs. Of 1,940 total pts, 901 had consent and sufficient tissue. HER2 by RT-PCR was neg in 130 (14%), equiv in 85 (9%), and pos in 686 (76%) pts. Concordance between HER2 assessments was 95% for RT-PCR vs central IHC (>10% + cells = +), 91% for RT-PCR vs central FISH (≥2.0 = pos) and 94% for central IHC vs central FISH. In the primary analysis, the association of HER2 expression with H benefit was marginally non-significant (P=0.057). In hormone receptor pos pts (local IHC) the association was significant (P=0.002). The association was nonlinear with the greatest estimated benefit at lower and higher HER2 mRNA expression levels. The observed treatment benefit in low HER2 pts was not due to imbalance between arms in RS and individual gene expression values. Conclusions: Concordance among HER2 assessments by central IHC, FISH, and RT-PCR was high. Association of HER2 mRNA expression with H benefit was marginally non-significant. A consistent benefit of trastuzumab irrespective of mHER2 levels was observed in the pts with either IHC+ or FISH+ tumors. Benefit was observed in pts with high HER2 by RT-PCR but also observed for the small groups of pts with negative results by quantitative RT-PCR or FISH (Table). Plausible mechanisms for this observation will be discussed. [Table: see text]
Monosomic Loss of MIR15A/MIR16-1 Is a Driver of Multiple Myeloma Proliferation and Disease Progression
Blood Cancer Discovery | March 3, 2020 | Volume: 1 | Issue: 1 | Pages: 68-81
Abstract The most common genetic abnormality in multiple myeloma is the deletion of chromosome 13, seen in almost half of newly diagnosed patients. Unlike chronic lymphocytic leukemia, where a recurrent minimally deleted region including MIR15A/MIR16-1 has been mapped, the deletions in multiple myeloma predominantly involve the entire chromosome and no specific driver gene has been identified. Additional candidate loci include RB1 and DIS3, but while biallelic deletion of RB1 is associated with disease progression, DIS3 is a common essential gene and complete inactivation is not observed. The Vk*MYC transgenic mouse model of multiple myeloma spontaneously acquires del(14), syntenic to human chromosome 13, and Rb1 complete inactivation, but not Dis3 mutations. Taking advantage of this model, we explored the role in multiple myeloma initiation and progression of two candidate loci on chromosome 13: RB1 and MIR15A/MIR16-1. Monoallelic deletion of Mir15a/Mir16-1, but not Rb1, was sufficient to accelerate the development of monoclonal gammopathy in wild-type mice and the progression of multiple myeloma in Vk*MYC mice, resulting in increased expression of Mir15a/Mir16-1 target genes and plasma cell proliferation, which was similarly observed in patients with multiple myeloma. Significance: In the absence of a defined, minimally deleted region the significance of del(13) in multiple myeloma has remained controversial. Here we show that haploinsufficiency of Mir15a/Mir16-1, but not Rb1, upregulates the cell cycle–regulatory network, inducing monoclonal gammopathy in mice and promoting multiple myeloma progression in both mice and men. See related commentary by Walker, p. 16. This article is highlighted in the In This Issue feature, p. 5