Multiomic Analysis Identifies a High-Risk Metabolic and TME Depleted Signature that Predicts Early Clinical Failure in DLBCL

PURPOSE

60-70% of newly diagnosed diffuse large B-cell lymphoma (DLBCL) patients avoid events within 24 months of diagnosis (EFS24) and the remainder have poor outcomes. Recent genetic and molecular classification of DLBCL has advanced our knowledge of disease biology, yet were not designed to predict early events and guide anticipatory selection of novel therapies. To address this unmet need, we used an integrative multiomic approach to identify a signature at diagnosis that will identify DLBCL at high risk of early clinical failure.

PATIENTS AND METHODS

Tumor biopsies from 444 newly diagnosed DLBCL were analyzed by WES and RNAseq. A combination of weighted gene correlation network analysis and differential gene expression analysis followed by integration with clinical and genomic data was used to identify a multiomic signature associated with high risk of early clinical failure.

RESULTS

Current DLBCL classifiers are unable to discriminate cases who fail EFS24. We identified a high risk RNA signature that had a hazard ratio (HR, 18.46 [95% CI 6.51-52.31] P < .001) in a univariate model, which did not attenuate after adjustment for age, IPI and COO (HR, 20.8 [95% CI, 7.14-61.09] P < .001). Further analysis revealed the signature was associated with metabolic reprogramming and a depleted immune microenvironment. Finally, WES data was integrated into the signature and we found that inclusion of ARID1A mutations resulted in identification of 45% of cases with an early clinical failure which was validated in external DLBCL cohorts.

CONCLUSION

This novel and integrative approach is the first to identify a signature at diagnosis that will identify DLBCL at high risk for early clinical failure and may have significant implications for design of therapeutic options.

AT-101 Enhances the Antitumor Activity of Lenalidomide in Patients with Multiple Myeloma

Bcl-2 and Mcl-1 proteins play a role in multiple myeloma (MM) cell survival, for which targeted inhibitors are being developed. AT-101 is an oral drug, which disrupts Bcl-2 and Mcl-1 function, impedes mitochondrial bioenergetic processes and induces apoptosis in MM cells. When combined with lenalidomide and dexamethasone (Rd), AT-101 significantly reduced tumor burden in an in vivo xenograft model of MM. These data provided rationale for a phase I/II study to establish the effective dose of AT-101 in combination with Rd (ARd regimen) in relapsed/refractory MM. A total of 10 patients were enrolled, most with high-risk cytogenetics (80%) and prior stem cell transplant (70%). Three patients were lenalidomide-refractory, 2 were bortezomib-refractory and 3 were daratumumab-refractory. The ARd combination was well tolerated with most common grade 3/4 adverse events being cytopenia's. The overall response rate was 40% and clinical benefit rate was 90%. The median progression free survival was 14.9 months (95% CI 7.1-NE). Patients responsive to ARd showed a decrease in Bcl-2:Bim or Mcl-1:Noxa protein complexes, increased CD8+ T and NK cells and depletion of T and B-regulatory cells. The ARd regimen demonstrated an acceptable safety profile and promising efficacy in patients with relapsed/refractory MM prompting further investigation in additional patients.

Loss of Tmem106b leads to cerebellum Purkinje cell death and motor deficits

TMEM106B has been recently implicated in multiple neurodegenerative diseases. Here, Rademakers et al. report a late-onset cerebellar Purkinje cell loss and progressive decline in motor function and gait deficits in a conventional Tmem106b-/- mouse model. By using high-power microscopy and bulk RNA sequencing, the authors further identify lysosomal and immune dysfunction as potential underlying mechanisms of the Purkinje cell loss.

Targeting of inflammatory pathways with R2CHOP in high-risk DLBCL

Diffuse large B-cell lymphoma (DLBCL) is the most common lymphoma, and front line therapies have not improved overall outcomes since the advent of immunochemotherapy. By pairing DNA and gene expression data with clinical response data, we identified a high-risk subset of non-GCB DLBCL patients characterized by genomic alterations and expression signatures capable of sustaining an inflammatory environment. These mutational alterations (PIM1, SPEN, and MYD88 [L265P]) and expression signatures (NF-κB, IRF4, and JAK-STAT engagement) were associated with proliferative signaling, and were found to be enriched in patients treated with RCHOP that experienced unfavorable outcomes. However, patients with these high-risk mutations had more favorable outcomes when the immunomodulatory agent lenalidomide was added to RCHOP (R2CHOP). We are the first to report the genomic validation of a high-risk phenotype with a preferential response towards R2CHOP therapy in non-GCB DLBCL patients. These conclusions could be translated to a clinical setting to identify the ~38% of non-GCB patients that could be considered high-risk, and would benefit from alternative therapies to standard RCHOP based on personalized genomic data.

Abstract 1574: The role of human cerebrospinal fluid in glioblastoma malignancy

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.

Abstract 3780: Cerebrospinal fluid-responsive factor SERPINA3 induces an increase in the malignancy of glioblastoma

Glioblastoma (GBM) is the most common and aggressive primary brain tumor in adults with a median survival of 14 months despite multimodal treatment. Patients with GBM tumors exhibit a post-resection recurrence rate of almost 100%. The tumor recurrence and treatment resistance properties of GBM reside in a stem-like cell population called Brain Tumor Initiating Cells (BTIC). Among primary GBMs, periventricular tumors exhibit worse survival than tumors distal to the lateral ventricles (LV), the cause for this worse outcome is not known. A potential explanation is the proximity of these tumors to the cerebrospinal fluid (CSF) and the neurogenic niche in the sub-ventricular zone (SVZ). Chemotactic signals present in the cerebrospinal fluid (CSF) regulate neural progenitor cells migratory response and differentiation profile. Our preliminary data show that cancer-derived CSF induces an increase in the migration and proliferation of GBM cells when compared to non-cancer CSF. We evaluated the transcriptome of human GBM cells in response to cancer and non-cancer CSF and identified alpha-1-antichymotrypsin (SERPINA3) as one of the genes that show higher over expression upon cancer CSF-stimulation. SERPINA3 is a circulating anti-protease that is produced and released in different tissues such as liver, prostate, lungs, endometrium and brain that is over-expressed in multiple cancers. We hypothesize that the overexpression of SERPINA3 in response to CSF is a contributing factor to the increase malignancy of periventricular GBM tumors. We tested this hypothesis in silico, in vitro and in vivo. Using TCGA and Rembrandt databases we observed a higher expression of SERPINA3 in brain tumor tissue compared to controls with no differences among GBM subtypes. In GBM patients, expression of SERPINA3 negatively affects survival expectancy. In vitro, Using intraoperative samples, we confirmed that SERPINA3 is overexpressed in brain cancer tissue by western-blot and qRT-PCR. Interestingly, we observed a higher nuclear localization of SERPINA3 in cancer tissue and cells. To determine if SERPINA3 has an effect on the malignancy of GBM cells, we silenced its expression using shRNA. SERPINA3 silencing induced a decrease in cell proliferation (p&lt;0.05), as well as migration (p&lt;0.05), and invasion (p&lt;0.05). Remarkably CD133 expression and stem cell markers were decreased when SERPINA3 is silenced. In contrast, the addition of recombinant SERPINA3 induced an increase in cell migration (p&lt;0.05) of GBM cells. Moreover, cancer CSF is no longer able to induce an increase in GBM cell migration in SERPINA3-KD cells (p&lt;0.05). In vivo, mice orthotopically-injected with SERPINA3-KD cells showed a prolonged overall survival (p&lt;0.05). In summary, our results support the hypothesis that SERPINA3 plays a key role in GBM malignancy in response to CSF and its inhibition results in a better outcome using GBM preclinical models. Targeting the CSF-GBM interaction is a potential therapeutic approach that remains to be explored.

Citation Format: Montserrat A. Lara Velazquez, Natanael Zarco, Anna Carrano, Jordan Phillips, Paula Schiapparelli, Emily S. Norton, Stephanie Jeanneret, Teresita Corona, Jose Segovia, Kaisorn Chaichana, Yan Asmann, Alfredo Quinones Hinojosa, Hugo Guerrero-Cazares. Cerebrospinal fluid-responsive factor SERPINA3 induces an increase in the malignancy of glioblastoma [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 3780.

Loss of TMEM106B leads to myelination deficits: implications for frontotemporal dementia treatment strategies

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.

Inherited variants at 3q13.33 and 3p24.1 are associated with risk of diffuse large B-cell lymphoma and implicate immune pathways

We previously identified five single nucleotide polymorphisms (SNPs) at four susceptibility loci for diffuse large B-cell lymphoma (DLBCL) in individuals of European ancestry through a large genome-wide association study (GWAS). To further elucidate genetic susceptibility to DLBCL, we sought to validate two loci at 3q13.33 and 3p24.1 that were suggestive in the original GWAS with additional genotyping. In the meta-analysis (5662 cases and 9237 controls) of the four original GWAS discovery scans and three replication studies, the 3q13.33 locus (rs9831894; minor allele frequency [MAF] = 0.40) was associated with DLBCL risk [odds ratio (OR) = 0.83, P = 3.62 × 10-13]. rs9831894 is in linkage disequilibrium (LD) with additional variants that are part of a super-enhancer that physically interacts with promoters of CD86 and ILDR1. In the meta-analysis (5510 cases and 12 817 controls) of the four GWAS discovery scans and four replication studies, the 3p24.1 locus (rs6773363; MAF = 0.45) was also associated with DLBCL risk (OR = 1.20, P = 2.31 × 10-12). This SNP is 29 426-bp upstream of the nearest gene EOMES and in LD with additional SNPs that are part of a highly lineage-specific and tumor-acquired super-enhancer that shows long-range interaction with AZI2 promoter. These loci provide additional evidence for the role of immune function in the etiology of DLBCL, the most common lymphoma subtype.

Managing genomic variant calling workflows with Swift/T

Bioinformatics research is frequently performed using complex workflows with multiple steps, fans, merges, and conditionals. This complexity makes management of the workflow difficult on a computer cluster, especially when running in parallel on large batches of data: hundreds or thousands of samples at a time. Scientific workflow management systems could help with that. Many are now being proposed, but is there yet the “best” workflow management system for bioinformatics? Such a system would need to satisfy numerous, sometimes conflicting requirements: from ease of use, to seamless deployment at peta- and exa-scale, and portability to the cloud. We evaluated Swift/T as a candidate for such role by implementing a primary genomic variant calling workflow in the Swift/T language, focusing on workflow management, performance and scalability issues that arise from production-grade big data genomic analyses. In the process we introduced novel features into the language, which are now part of its open repository. Additionally, we formalized a set of design criteria for quality, robust, maintainable workflows that must function at-scale in a production setting, such as a large genomic sequencing facility or a major hospital system. The use of Swift/T conveys two key advantages. (1) It operates transparently in multiple cluster scheduling environments (PBS Torque, SLURM, Cray aprun environment, etc.), thus a single workflow is trivially portable across numerous clusters. (2) The leaf functions of Swift/T permit developers to easily swap executables in and out of the workflow, which makes it easy to maintain and to request resources optimal for each stage of the pipeline. While Swift/T’s data-level parallelism eliminates the need to code parallel analysis of multiple samples, it does make debugging more difficult, as is common for implicitly parallel code. Nonetheless, the language gives users a powerful and portable way to scale up analyses in many computing architectures. The code for our implementation of a variant calling workflow using Swift/T can be found on GitHub at https://github.com/ncsa/Swift-T-Variant-Calling, with full documentation provided at http://swift-t-variant-calling.readthedocs.io/en/latest/.

Abstract 651: Analysis of the immune microenvironment to advance breast cancer risk prediction and prevention

Existing risk models underperform among women who have undergone a benign breast disease (BBD) biopsy (&gt; one million performed annually in the US) and with respect to estimating the risk of aggressive BCs. We've shown that molecular pathologic analysis of BBD biopsies can improve individual risk prediction, compared with standard risk models relying on self-reported factors, and provide insights into mechanisms mediating BC risk. BC risk factors such as obesity and ethanol use are proposed to increase production of cytokines and chemokines, resulting in chronic inflammation, and leading to production of DNA damaging free radicals and growth factors (e.g., VEGF, IGFs) that activate pro-carcinogenic pathways (e.g. NF-KB and JAK/STAT). The effects of most BC risk factors on immunity are poorly understood. Further, whereas experimental models implicate immunity throughout carcinogenesis, knowledge of immune markers and mechanisms related to BC development among women is limited. Therefore, we aim to define tissue immune cell content throughout BC development, and particularly at the earliest stages, to improve risk assessment and discover immune-based prevention strategies.

Our study combines novel resources and technologies to define the immune landscape in: 1) normal breast tissues in relation to BC risk factors, 2) BBD with respect to BC risk, and 3) invasive BC in relation to neoantigen expression and molecular subtype. We will perform NanoString GeoMxTM Digital Spatial Profiling (DSP) to quantitatively map expression of key immune proteins in healthy breast tissues donated to the Komen Tissue Bank (KTB) in relation to BC risk factors and in BBD biopsies from two cohorts to predict BC risk. We will re-analyze top prognostic markers in BBD biopsies using Vectra multiplex IF staining with machine learning algorithms to refine how immune cell content and organization impacts prognosis, and assess critical immune pathways related to BBD progression with the NanoString PanCancer IO 360 panel, which provides RNA expression of 770 immune genes (13 signatures). Finally, we will evaluate immune responses in BCs categorized for predicted neoantigen load and underlying mutation type with our novel bioinformatics pipeline that enables accurate prediction of MHC class I and class II missense, fusion and frameshift mutation neoantigens generated through faulty DNA repair, aberrant DNA and RNA splicing, insertions, and deletions.

Our project demonstrates a novel approach to defining the immune processes in early stages of BC development which may improve BC risk prediction and prevention. We present a plan to comprehensively analyze immunity throughout BC development and define changes that accompany transitions from normal to BBD and BBD to BC.

Citation Format: Doug Hinerfeld, Keith Knutson, Derek Radisky, E. Aubrey Thompson, Yan Asmann, Kim McCoy, Amy Degnim, Jodi Carter, Stacey Winham, Michele Cote, Jeroen van der Laak, Mark Sherman. Analysis of the immune microenvironment to advance breast cancer risk prediction and prevention [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 651.

Genome-wide analyses as part of the international FTLD-TDP whole-genome sequencing consortium reveals novel disease risk factors and increases support for immune dysfunction in FTLD

Frontotemporal lobar degeneration with neuronal inclusions of the TAR DNA-binding protein 43 (FTLD-TDP) represents the most common pathological subtype of FTLD. We established the international FTLD-TDP whole-genome sequencing consortium to thoroughly characterize the known genetic causes of FTLD-TDP and identify novel genetic risk factors. Through the study of 1131 unrelated Caucasian patients, we estimated that C9orf72 repeat expansions and GRN loss-of-function mutations account for 25.5% and 13.9% of FTLD-TDP patients, respectively. Mutations in TBK1 (1.5%) and other known FTLD genes (1.4%) were rare, and the disease in 57.7% of FTLD-TDP patients was unexplained by the known FTLD genes. To unravel the contribution of common genetic factors to the FTLD-TDP etiology in these patients, we conducted a two-stage association study comprising the analysis of whole-genome sequencing data from 517 FTLD-TDP patients and 838 controls, followed by targeted genotyping of the most associated genomic loci in 119 additional FTLD-TDP patients and 1653 controls. We identified three genome-wide significant FTLD-TDP risk loci: one new locus at chromosome 7q36 within the DPP6 gene led by rs118113626 (p value = 4.82e - 08, OR = 2.12), and two known loci: UNC13A, led by rs1297319 (p value = 1.27e - 08, OR = 1.50) and HLA-DQA2 led by rs17219281 (p value = 3.22e - 08, OR = 1.98). While HLA represents a locus previously implicated in clinical FTLD and related neurodegenerative disorders, the association signal in our study is independent from previously reported associations. Through inspection of our whole-genome sequence data for genes with an excess of rare loss-of-function variants in FTLD-TDP patients (n ≥ 3) as compared to controls (n = 0), we further discovered a possible role for genes functioning within the TBK1-related immune pathway (e.g., DHX58, TRIM21, IRF7) in the genetic etiology of FTLD-TDP. Together, our study based on the largest cohort of unrelated FTLD-TDP patients assembled to date provides a comprehensive view of the genetic landscape of FTLD-TDP, nominates novel FTLD-TDP risk loci, and strongly implicates the immune pathway in FTLD-TDP pathogenesis.

PCNT point mutations and familial intracranial aneurysms

OBJECTIVE

To identify novel genes involved in the etiology of intracranial aneurysms (IAs) or subarachnoid hemorrhages (SAHs) using whole-exome sequencing.

METHODS

We performed whole-exome sequencing in 13 individuals from 3 families with an autosomal dominant IA/SAH inheritance pattern to look for candidate genes for disease. In addition, we sequenced PCNT exon 38 in a further 161 idiopathic patients with IA/SAH to find additional carriers of potential pathogenic variants.

RESULTS

We identified 2 different variants in exon 38 from the PCNT gene shared between affected members from 2 different families with either IA or SAH (p.R2728C and p.V2811L). One hundred sixty-four samples with either SAH or IA were Sanger sequenced for the PCNT exon 38. Five additional missense mutations were identified. We also found a second p.V2811L carrier in a family with a history of neurovascular diseases.

CONCLUSION

The PCNT gene encodes a protein that is involved in the process of microtubule nucleation and organization in interphase and mitosis. Biallelic loss-of-function mutations in PCNT cause a form of primordial dwarfism (microcephalic osteodysplastic primordial dwarfism type II), and ≈50% of these patients will develop neurovascular abnormalities, including IAs and SAHs. In addition, a complete Pcnt knockout mouse model (Pcnt ^-/-) published previously showed general vascular abnormalities, including intracranial hemorrhage. The variants in our families lie in the highly conserved PCNT protein-protein interaction domain, making PCNT a highly plausible candidate gene in cerebrovascular disease.

Differences in genomic abnormalities among African individuals with monoclonal gammopathies using calculated ancestry

Multiple myeloma (MM) is two- to three-fold more common in African Americans (AAs) compared to European Americans (EAs). This striking disparity, one of the highest of any cancer, may be due to underlying genetic predisposition between these groups. There are multiple unique cytogenetic subtypes of MM, and it is likely that the disparity is associated with only certain subtypes. Previous efforts to understand this disparity have relied on self-reported race rather than genetic ancestry, which may result in bias. To mitigate these difficulties, we studied 881 patients with monoclonal gammopathies who had undergone uniform testing to identify primary cytogenetic abnormalities. DNA from bone marrow samples was genotyped on the Precision Medicine Research Array and biogeographical ancestry was quantitatively assessed using the Geographic Population Structure Origins tool. The probability of having one of three specific subtypes, namely t(11;14), t(14;16), or t(14;20) was significantly higher in the 120 individuals with highest African ancestry (≥80%) compared with the 235 individuals with lowest African ancestry (<0.1%) (51% vs. 33%, respectively, p value = 0.008). Using quantitatively measured African ancestry, we demonstrate a major proportion of the racial disparity in MM is driven by disparity in the occurrence of the t(11;14), t(14;16), and t(14;20) types of MM.

Comprehensive Screening for Disease Risk Variants in Early-Onset Alzheimer's Disease Genes in African Americans Identifies Novel PSEN Variants

We conducted a comprehensive screening of rare coding variants in an African American cohort to identify novel pathogenic mutations within the early-onset Alzheimer's disease (EOAD) genes (APP, PSEN1, and PSEN2) in this understudied population. Whole-exome sequencing of 238 African American subjects identified 6 rare missense variants within the EOAD genes, which were observed in AD cases but never among controls. These variants were analyzed in an independent cohort of 300 African American subjects in which PSEN2:NM_000447:exon5:c.T331C:p.Phe111Leu and PSEN1-minilin rs777923890 variants were again not observed, indicating that these novel rare variants, may contribute to AD risk in this population.

Lipocalin-2 protects the brain during inflammatory conditions

Sepsis is a prevalent health issue that can lead to central nervous system (CNS) inflammation with long-term behavioral and cognitive alterations. Using unbiased proteomic profiling of over 100 different cytokines, we found that Lipocalin-2 (LCN2) was the most substantially elevated protein in the CNS after peripheral administration of lipopolysaccharide (LPS). To determine whether the high level of LCN2 in the CNS is protective or deleterious, we challenged Lcn2^-/- mice with peripheral LPS and determined effects on behavior and neuroinflammation. At a time corresponding to peak LCN2 induction in wild-type (WT) mice injected with LPS, Lcn2^-/- mice challenged with LPS had exacerbated levels of pro-inflammatory cytokines and exhibited significantly worsened behavioral phenotypes. To determine the extent of global inflammatory changes dependent upon LCN2, we performed an RNAseq transcriptomic analysis. Compared with WT mice injected with LPS, Lcn2^-/- mice injected with LPS had unique transcriptional profiles and significantly elevated levels of multiple pro-inflammatory molecules. Several LCN2-dependent pathways were revealed with this analysis including, cytokine and chemokine signaling, nucleotide-binding oligomerization domain-like receptor signaling and Janus kinase-signal transducer and activator of transcription signaling. These findings demonstrate that LCN2 serves as a potent protective factor in the CNS in response to systemic inflammation and may be a potential candidate for limiting sepsis-related CNS sequelae.

High somatic mutation and neoantigen burden are correlated with decreased progression-free survival in multiple myeloma

Tumor-specific mutations can result in immunogenic neoantigens, both of which have been correlated with responsiveness to immune checkpoint inhibitors in highly mutagenic cancers. However, early results of single-agent checkpoint inhibitors in multiple myeloma (MM) have been underwhelming. Therefore, we sought to understand the relationship between mutation and neoantigen landscape of MM patients and responsiveness to therapies. Somatic mutation burden, neoantigen load, and response to therapy were determined using interim data from the MMRF CoMMpass study (NCT01454297) on 664 MM patients. In this population, the mean somatic and missense mutation loads were 405.84(s=608.55) and 63.90(s=95.88) mutations per patient, respectively. There was a positive linear relationship between mutation and neoantigen burdens (R²=0.862). The average predicted neoantigen load was 23.52(s=52.14) neoantigens with an average of 9.40(s=26.97) expressed neoantigens. Survival analysis revealed significantly shorter progression-free survival (PFS) in patients with greater than average somatic missense mutation load (N=163, 0.493 vs 0.726 2-year PFS, P=0.0023) and predicted expressed neoantigen load (N=214, 0.555 vs 0.729 2-year PFS, P=0.0028). This pattern is maintained when stratified by disease stage and cytogenetic abnormalities. Therefore, high mutation and neoantigen load are clinically relevant risk factors that negatively impact survival of MM patients under current standards of care.

BAP1 dependent expression of long non-coding RNA NEAT-1 contributes to sensitivity to gemcitabine in cholangiocarcinoma

BACKGROUND

Genetic alterations in chromatin modulators such as BRCA-1 associated protein-1 (BAP1) are the most frequent genetic alteration in intrahepatic cholangiocarcinomas (CCA). We evaluated the contribution of BAP1 expression on tumor cell behavior and therapeutic sensitivity to identify rationale therapeutic strategies.

METHODS

The impact of BAP1 expression on sensitivity to therapeutic agents was evaluated in CCA cells with a 7-fold difference in BAP1 expression (KMBC-low, HuCCT1-high) and genetically engineered haplo-insufficient BAP1 knockout cells. We also identified long non-coding RNA genes associated with loss of BAP1 and their role in therapeutic sensitivity.

RESULTS

Sensitivity to gemcitabine was greater in low BAP1 expressing or BAP1 knockout cells compared with the high BAP1 expressing cells or control haplo-insufficient cells respectively. Similar results were observed with TSA, olaparib, b-AP15 but not with GSK126. A differential synergistic effect was observed in combinations of gemcitabine with olaparib or GSK126 in KMBC cells and TSA or bAP15 in HuCCT1 cells, indicating BAP1 dependent target-specific synergism and sensitivity to gemcitabine. A BAP1 dependent alteration in expression of lncRNA NEAT-1 was identified by RT-PCR based lncRNA expression profiling, and an inverse relationship between this lncRNA and BAP1 was observed in analysis of the Tumor Cancer Genome Atlas cholangiocarcinoma dataset. Exogenous modulation of NEAT-1 and/or BAP1 expression altered tumor cell phenotype and modulated sensitivity to gemcitabine.

CONCLUSIONS

NEAT-1 is a downstream effector of gemcitabine sensitivity in CCA. The expression of BAP1 is a determinant of sensitivity to therapeutic drugs that can be exploited to enhance responses through combination strategies.

Human whole genome genotype and transcriptome data for Alzheimer's and other neurodegenerative diseases

Previous genome-wide association studies (GWAS), conducted by our group and others, have identified loci that harbor risk variants for neurodegenerative diseases, including Alzheimer's disease (AD). Human disease variants are enriched for polymorphisms that affect gene expression, including some that are known to associate with expression changes in the brain. Postulating that many variants confer risk to neurodegenerative disease via transcriptional regulatory mechanisms, we have analyzed gene expression levels in the brain tissue of subjects with AD and related diseases. Herein, we describe our collective datasets comprised of GWAS data from 2,099 subjects; microarray gene expression data from 773 brain samples, 186 of which also have RNAseq; and an independent cohort of 556 brain samples with RNAseq. We expect that these datasets, which are available to all qualified researchers, will enable investigators to explore and identify transcriptional mechanisms contributing to neurodegenerative diseases.

Gene expression, methylation and neuropathology correlations at progressive supranuclear palsy risk loci

To determine the effects of single nucleotide polymorphisms (SNPs) identified in a genome-wide association study of progressive supranuclear palsy (PSP), we tested their association with brain gene expression, CpG methylation and neuropathology. In 175 autopsied PSP subjects, we performed associations between seven PSP risk variants and temporal cortex levels of 20 genes in-cis, within ±100 kb. Methylation measures were collected using reduced representation bisulfite sequencing in 43 PSP brains. To determine whether SNP/expression associations are due to epigenetic modifications, CpG methylation levels of associated genes were tested against relevant variants. Quantitative neuropathology endophenotypes were tested for SNP associations in 422 PSP subjects. Brain levels of LRRC37A4 and ARL17B were associated with rs8070723; MOBP with rs1768208 and both ARL17A and ARL17B with rs242557. Expression associations for LRRC37A4 and MOBP were available in an additional 100 PSP subjects. Meta-analysis revealed highly significant associations for PSP risk alleles of rs8070723 and rs1768208 with higher LRRC37A4 and MOBP brain levels, respectively. Methylation levels of one CpG in the 3' region of ARL17B associated with rs242557 and rs8070723. Additionally, methylation levels of an intronic ARL17A CpG associated with rs242557 and that of an intronic MOBP CpG with rs1768208. MAPT and MOBP region risk alleles also associated with higher levels of neuropathology. Strongest associations were observed for rs242557/coiled bodies and tufted astrocytes; and for rs1768208/coiled bodies and tau threads. These findings suggest that PSP variants at MAPT and MOBP loci may confer PSP risk via influencing gene expression and tau neuropathology. MOBP, LRRC37A4, ARL17A and ARL17B warrant further assessment as candidate PSP risk genes. Our findings have implications for the mechanism of action of variants at some of the top PSP risk loci.

Frequent occurrence of large duplications at reciprocal genomic rearrangement breakpoints in multiple myeloma and other tumors

Using a combination of array comparative genomic hybridization, mate pair and cloned sequences, and FISH analyses, we have identified in multiple myeloma cell lines and tumors a novel and recurrent type of genomic rearrangement, i.e. interchromosomal rearrangements (translocations or insertions) and intrachromosomal inversions that contain long (1-4000 kb; median ∼100 kb) identical sequences adjacent to both reciprocal breakpoint junctions. These duplicated sequences were generated from sequences immediately adjacent to the breakpoint from at least one-but sometimes both-chromosomal donor site(s). Tandem duplications had a similar size distribution suggesting the possibility of a shared mechanism for generating duplicated sequences at breakpoints. Although about 25% of apparent secondary rearrangements contained these duplications, primary IGH translocations rarely, if ever, had large duplications at breakpoint junctions. Significantly, these duplications often contain super-enhancers and/or oncogenes (e.g. MYC) that are dysregulated by rearrangements during tumor progression. We also found that long identical sequences often were identified at both reciprocal breakpoint junctions in six of eight other tumor types. Finally, we have been unable to find reports of similar kinds of rearrangements in wild-type or mutant prokaryotes or lower eukaryotes such as yeast.

Whole exome sequencing of a patient with metastatic hidradenocarcinoma and review of the literature

Hidradenocarcinoma is a rare malignancy of the sweat glands with only a few cases reported in literature. The management of these tumors is based on the extent of disease with local disease managed with surgical resection. These can tumors carry a high potential of lymphatic and vascular spread and local and distant metastases are not uncommon. Given the rarity of the tumor and lack of genetic and clinical data about these tumors, there is no consensus on the proper management of metastatic disease. Here in we report the first case of metastatic hidradenocarcinoma with detailed molecular profiling including whole exome sequencing. We identified mutations in multiple genes including two that are potentially targetable: PTCH1 and TCF7L1. Further work is necessary to not only confirm the presence of these mutations but also to confirm the clinical significance.

Analysis of COQ2 gene in multiple system atrophy

Background

Loss of function COQ2 mutations results in primary CoQ10 deficiency. Recently, recessive mutations of the COQ2 gene have been identified in two unrelated Japanese families with multiple system atrophy (MSA). It has also been proposed that specific heterozygous variants in the COQ2 gene may confer susceptibility to sporadic MSA. To assess the frequency of COQ2 variants in patients with MSA, we sequenced the entire coding region and investigated all exonic copy number variants of the COQ2 gene in 97 pathologically-confirmed and 58 clinically-diagnosed MSA patients from the United States.

Results

We did not find any homozygous or compound heterozygous pathogenic COQ2 mutations including deletion or multiplication within our series of MSA patients. In two patients, we identified two heterozygous COQ2 variants (p.S54W and c.403 + 10G > T) of unknown significance, which were not observed in 360 control subjects. We also identified one heterozygous carrier of a known loss of function p.S146N substitution in a severe MSA-C pathologically-confirmed patient.

Conclusions

The COQ2 p.S146N substitution has been previously reported as a pathogenic mutation in primary CoQ10 deficiency (including infantile multisystem disorder) in a recessive manner. This variant is the third primary CoQ10 deficiency mutation observed in an MSA case (p.R387X and p.R197H). Therefore it is possible that in the heterozygous state it may increase susceptibility to MSA. Further studies, including reassessing family history in patients of primary CoQ10 deficiency for the possible occurrence of MSA, are now warranted to resolve the role of COQ2 variation in MSA.

Electronic supplementary material

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

Integrated genomic characterization reveals novel, therapeutically relevant drug targets in FGFR and EGFR pathways in sporadic intrahepatic cholangiocarcinoma

Advanced cholangiocarcinoma continues to harbor a difficult prognosis and therapeutic options have been limited. During the course of a clinical trial of whole genomic sequencing seeking druggable targets, we examined six patients with advanced cholangiocarcinoma. Integrated genome-wide and whole transcriptome sequence analyses were performed on tumors from six patients with advanced, sporadic intrahepatic cholangiocarcinoma (SIC) to identify potential therapeutically actionable events. Among the somatic events captured in our analysis, we uncovered two novel therapeutically relevant genomic contexts that when acted upon, resulted in preliminary evidence of anti-tumor activity. Genome-wide structural analysis of sequence data revealed recurrent translocation events involving the FGFR2 locus in three of six assessed patients. These observations and supporting evidence triggered the use of FGFR inhibitors in these patients. In one example, preliminary anti-tumor activity of pazopanib (in vitro FGFR2 IC₅₀≈350 nM) was noted in a patient with an FGFR2-TACC3 fusion. After progression on pazopanib, the same patient also had stable disease on ponatinib, a pan-FGFR inhibitor (in vitro, FGFR2 IC₅₀≈8 nM). In an independent non-FGFR2 translocation patient, exome and transcriptome analysis revealed an allele specific somatic nonsense mutation (E384X) in ERRFI1, a direct negative regulator of EGFR activation. Rapid and robust disease regression was noted in this ERRFI1 inactivated tumor when treated with erlotinib, an EGFR kinase inhibitor. FGFR2 fusions and ERRFI mutations may represent novel targets in sporadic intrahepatic cholangiocarcinoma and trials should be characterized in larger cohorts of patients with these aberrations.

Cholangiocarcinoma is a cancer that affects the bile ducts. Unfortunately, many patients diagnosed with cholangiocarcinoma have disease that cannot be treated with surgery or has spread to other parts of the body, thus severely limiting treatment options. New advances in drug treatment have enabled treatment of these cancers with “targeted therapy” that exploits an error in the normal functioning of a tumor cell, compared to other cells in the body, thus allowing only tumor cells to be killed by the drug. We sought to identify changes in the genetic material of cholangiocarcinoma patient tumors in order to identify potential errors in cellular functioning by utilizing cutting edge genetic sequencing technology. We identified three patient tumors possessing an FGFR2 gene that was aberrantly fused to another gene. Two of these patients were able to receive targeted therapy for FGFR2 with resulting tumor shrinkage. A fourth tumor contained an error in a gene that controls a very important cellular mechanism in cancer, termed epidermal growth factor pathway (EGFR). This patient received therapy targeting this mechanism and also demonstrated response to treatment. Thus, we have been able to utilize cutting edge technology with targeted drug treatment to personalize medical treatment for cancer in cholangiocarcinoma patients.

Insulin fails to enhance mTOR phosphorylation, mitochondrial protein synthesis, and ATP production in human skeletal muscle without amino acid replacement

Systemic insulin administration causes hypoaminoacidemia by inhibiting protein degradation, which may in turn inhibit muscle protein synthesis (PS). Insulin enhances muscle mitochondrial PS and ATP production when hypoaminoacidemia is prevented by exogenous amino acid (AA) replacement. We determined whether insulin would stimulate mitochondrial PS and ATP production in the absence of AA replacement. Using l-[1,2-¹³C]leucine as a tracer, we measured the fractional synthetic rate of mitochondrial as well as sarcoplasmic and mixed muscle proteins in 18 participants during sustained (7-h) insulin or saline infusion (n = 9 each). We also measured muscle ATP production, mitochondrial enzyme activities, mRNA levels of mitochondrial genes, and phosphorylation of signaling proteins regulating protein synthesis. The concentration of circulating essential AA decreased during insulin infusion. Mitochondrial, sarcoplasmic, and mixed muscle PS rates were also lower during insulin (2-7 h) than during saline infusions despite increased mRNA levels of selected mitochondrial genes. Under these conditions, insulin did not alter mitochondrial enzyme activities and ATP production. These effects were associated with enhanced phosphorylation of Akt but not of protein synthesis activators mTOR, p70(S6K), and 4EBP1. In conclusion, sustained physiological hyperinsulinemia without AA replacement did not stimulate PS of mixed muscle or protein subfractions and did not alter muscle mitochondrial ATP production in healthy humans. These results support that insulin and AA act in conjunction to stimulate muscle mitochondrial function and mitochondrial protein synthesis.

Germline copy number variation and ovarian cancer survival

Copy number variants (CNVs) have been implicated in many complex diseases. We examined whether inherited CNVs were associated with overall survival among women with invasive epithelial ovarian cancer. Germline DNA from 1,056 cases (494 deceased, average of 3.7 years follow-up) was interrogated with the Illumina 610 quad genome-wide array containing, after quality control exclusions, 581,903 single nucleotide polymorphisms (SNPs) and 17,917 CNV probes. Comprehensive analysis capitalized upon the strengths of three complementary approaches to CNV classification. First, to identify small CNVs, single markers were evaluated and, where associated with survival, consecutive markers were combined. Two chromosomal regions were associated with survival using this approach (14q31.3 rs2274736 p = 1.59 × 10⁻⁶, p = 0.001; 22q13.31 rs2285164 p = 4.01 × 10⁻⁵, p = 0.009), but were not significant after multiple testing correction. Second, to identify large CNVs, genome-wide segmentation was conducted to characterize chromosomal gains and losses, and association with survival was evaluated by segment. Four regions were associated with survival (1q21.3 loss p = 0.005, 5p14.1 loss p = 0.004, 9p23 loss p = 0.002, and 15q22.31 gain p = 0.002); however, again, after correcting for multiple testing, no regions were statistically significant, and none were in common with the single marker approach. Finally, to evaluate associations with general amounts of copy number changes across the genome, we estimated CNV burden based on genome-wide numbers of gains and losses; no associations with survival were observed (p > 0.40). Although CNVs that were not well-covered by the Illumina 610 quad array merit investigation, these data suggest no association between inherited CNVs and survival after ovarian cancer.

LIM domain only 2 protein expression, LMO2 germline genetic variation, and overall survival in diffuse large B-cell lymphoma in the pre-rituximab era

Both LMO2 (LIM domain only 2) mRNA and protein expression in diffuse large B-cell lymphoma (DLBCL) have been associated with superior survival. However, a role for germline genetic variation in LMO2 has not been previously reported. Immunohistochemistry (IHC) for LMO2 was conducted on tumor tissue from diagnostic biopsies, and 20 tag single nucleotide polymorphisms (SNPs) from LMO2 were genotyped from germline DNA. LMO2 IHC positivity was associated with superior survival (hazard ratio [HR] = 0.55; 95% confidence interval [CI] 0.31-0.97). Four LMO2 SNPs (rs10836127, rs941940, rs750781, rs1885524) were associated with survival after adjusting for LMO2 IHC and clinical factors (p < 0.05), and one of these SNPs (rs941940) was also associated with IHC positivity (p = 0.02). Compared to a model with clinical factors only (c-statistic = 0.676), adding the four SNPs (c-statistic = 0.751) or LMO2 IHC (c-statistic = 0.691) increased the predictive ability of the model, while inclusion of all three factors (c-statistic = 0.754) did not meaningfully add predictive ability above a model with clinical factors and the four SNPs. In conclusion, germline genetic variation in LMO2 was associated with DLBCL prognosis and provided slightly stronger predictive ability relative to LMO2 IHC status.

Abstract 4975: Next generation sequencing reveals a connection between KRAS mutation and the NFkB pathway in lung adenocarcinoma samples

Motivation: KRAS is commonly mutated in a variety of cancers including lung cancer. The KRAS gene is frequently mutated at codons 12 and 13 in lung adenocarcinomas in patients with a history of smoking. Tumors harboring an activating KRAS mutation are aggressive and are often resistant to available therapies. In the present study, we set out to identify pathways that are specifically altered in lung cancer patients whose tumors harbor a KRAS mutation using next generation sequencing technology.

Methods: We performed 50nt paired-end RNA-sequencing in 15 lung adenocarcinoma samples (8 and 7 samples with and without KRAS mutation). Reads were aligned to genome and exon junctions using Illumina's alignment tool ELAND_RNA. CASAVA and Genome Studio data analysis software were used to obtain read counts for genes, exons and exon junctions. We applied a variety of computational methods (Casper R package, Bowtie, TopHat) and softwares (Partek, R statistical software, JMP, Ingenuity Pathway Anlaysis, DAVID) to carry out our analyses.

Results: Differential gene expression and splicing analysis were performed between the two groups. We identified 115 genes that were differentially expressed and 112 genes that consist of splicing variants with at least 2 fold changes and p-value &lt; 0.01. We randomly selected 6/15 samples and performed real-time qRT-PCR for 6 differentially expressed genes. Significant correlations were observed ranging from 0.44 to 0.83, when qRT-PCR results were compared with RNA-sequencing expression data. Pathway analysis with 115 differentially expressed genes and 112 splicing variants revealed that the most significant pathway is composed entirely of NFκB focus genes, suggesting that there is a direct connection between oncogenic KRAS and activation of the NFκB signaling pathway.

Conclusions: NFκB has been implicated in KRas-mediated formation of early stage lung adenomas in the LSL-Kas mouse model. Our data indicate that this connection also occurs in human lung cancer and establish NFκB activation as a key manifestation of KRAS mutation in human lung cancer. Hence drugs targeting NFκB and NFκB related genes may potentially be helpful for the treatment of patients with oncogenic KRAS mutations. (Supported in part by grants from the 26.2 with Donna Foundation and the Florida Department of Health Bankhead/Coley program.)

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 4975. doi:10.1158/1538-7445.AM2011-4975

Genome-wide linkage analysis for uric acid in families enriched for hypertension

BACKGROUND

Uric acid is heritable and associated with hypertension and insulin resistance. We sought to identify genomic regions influencing serum uric acid in families in which two or more siblings had hypertension.

METHODS

Uric acid levels and microsatellite markers were assayed in the Genetic Epidemiology Network of Arteriopathy (GENOA) cohort (1075 whites and 1333 blacks) and the Hypertension Genetic Epidemiology Network (HyperGEN) cohort (1542 whites and 1627 blacks). Genome-wide linkage analyses of uric acid and bivariate linkage analyses of uric acid with an additional surrogate of insulin resistance were completed. Pathway analysis explored gene sets enriched at loci influencing uric acid.

RESULTS

In the GENOA white cohort, loci influencing uric acid were identified on chromosome 8 at 135 cM [multipoint logarithm of odds score (MLS) = 2.4], on chromosome 9 at 113 cM (MLS = 3.7) and on chromosome 16 at 93 cM (MLS = 2.3), but did not replicate in HyperGEN. At these loci, there was evidence of pleiotropy with other surrogates of insulin resistance and genes in the fructose and mannose metabolism pathway were enriched. In the HyperGEN-black cohort, there was some evidence of a locus for uric acid on chromosome 4 at 135 cM (MLS = 2.4) that had modest replication in GENOA (MLS = 1.2).

CONCLUSIONS

Several novel loci linked to uric acid were identified but none showed clear replication. Widespread diuretic use, a medication that raises uric acid levels, was an important study limitation. Bivariate linkage analyses and pathway analysis were consistent with genes regulating insulin resistance and fructose metabolism contributing to the heritability of uric acid.