Clonal selection and double-hit events involving tumor suppressor genes underlie relapse in myeloma.

To elucidate the mechanisms underlying relapse from chemotherapy in multiple myeloma, we performed a longitudinal study of 33 patients entered into Total Therapy protocols investigating them using gene expression profiling, high-resolution copy number arrays, and whole-exome sequencing. The study illustrates the mechanistic importance of acquired mutations in known myeloma driver genes and the critical nature of biallelic inactivation events affecting tumor suppressor genes, especially TP53, the end result being resistance to apoptosis and increased proliferation rates, which drive relapse by Darwinian-type clonal evolution. The number of copy number aberration changes and biallelic inactivation of tumor suppressor genes was increased in GEP70 high risk, consistent with genomic instability being a key feature of high risk. In conclusion, the study highlights the impact of acquired genetic events, which enhance the evolutionary fitness level of myeloma-propagating cells to survive multiagent chemotherapy and to result in relapse.

Host genetic susceptibility to Clostridium difficile infections in patients undergoing autologous stem cell transplantation: a genome-wide association study.

BACKGROUND: Clostridium difficile infection (CDI) is the most common hospital-acquired infection. Unfortunately, genes that identify CDI-susceptible patients have not been well described. We performed a genome-wide association study (GWAS) to determine genetic variants associated with the development of CDI. METHODS: A cohort study of Caucasian patients undergoing autologous stem cell transplantation for multiple myeloma was performed. Patients were genotyped using Illumina® Whole Genome Genotyping Infinium chemistry. We then compared CDI-positive to CDI-negative patients using logistic regression for baseline clinical factors and false discovery rate (FDR) for genetic factors [single nucleotide polymorphisms (SNPs)]. SNPs associated with CDI at FDR of p < 0.01 were then incorporated into a logistic regression model combining clinical and genetic factors. RESULTS: Of the 646 patients analyzed (59.7% male), 57 patients were tested CDI positive (cases) and were compared to 589 patients who were tested negative (controls). Hemoglobin, albumin, and hematocrit were lower for cases (p < 0.05). Eight SNPs on five genes (FLJ16171, GORASP2, RLBP1L1, ASPH, ATP7B) were associated with CDI at FDR p < 0.01. In the combined clinical and genetic model, low albumin and three genes RLBP1L1, ASPH, and ATP7B were associated with CDI. CONCLUSION: Low serum albumin and genes RLBP1L1 and ASPH located on chromosome 8 and ATP7B on chromosome 13 were associated with CDI. Of particular interest is ATP7B given its copper modulatory role and the sporicidal properties of copper against Clostridium difficile.

Frequent engagement of the classical and alternative NF-kappaB pathways by diverse genetic abnormalities in multiple myeloma.

Mechanisms of constitutive NF-kappaB signaling in multiple myeloma are unknown. An inhibitor of IkappaB kinase beta (IKKbeta) targeting the classical NF-kappaB pathway was lethal to many myeloma cell lines. Several cell lines had elevated expression of NIK due to genomic alterations or protein stabilization, while others had inactivating mutations of TRAF3; both kinds of abnormality triggered the classical and alternative NF-kappaB pathways. A majority of primary myeloma patient samples and cell lines had elevated NF-kappaB target gene expression, often associated with genetic or epigenetic alteration of NIK, TRAF3, CYLD, BIRC2/BIRC3, CD40, NFKB1, or NFKB2. These data demonstrate that addiction to the NF-kappaB pathway is frequent in myeloma and suggest that IKKbeta inhibitors hold promise for the treatment of this disease.

GWAS of 972 autologous stem cell recipients with multiple myeloma identifies 11 genetic variants associated with chemotherapy-induced oral mucositis.

PURPOSE: High-dose chemotherapy and autologous stem cell transplant (ASCT) to treat multiple myeloma (MM) and other cancers carries the risk of oral mucositis (OM) with sequelae including impaired nutritional and fluid intake, pain, and infectious complications. As a result of these problems, cancer treatment may have to be interrupted or delayed. In this study, we looked beyond OM's known risk factors of renal function and melphalan dose with a genome-wide association study (GWAS) to evaluate whether genetic variants in conjunction with clinical risk factors influence predisposition for OM. METHODS: Genotyping was performed using Illumina HumanOmni1-Quad v1.0 BeadChip and further assessed for data quality. We tested 892,589 germline single-nucleotide polymorphisms (SNPs) for association with OM among 972 Caucasian patients treated with high-dose melphalan and ASCT in Total Therapy clinical trials (TT2, TT3, TT4) for newly diagnosed MM. Statistical analyses included t tests, stepwise regression modeling, and logistic regression modeling to find baseline clinical factors and genotypes associated with OM. RESULTS: We found that 353 (36.3 %) patients had grades 2-4 OM. Type of treatment protocol, baseline estimated glomerular filtration rate, and melphalan dose along with baseline serum albumin and female gender predicted 43.6 % of grades 2-4 OM cases. Eleven SNPs located in or near matrix metalloproteinase 13, JPH3, DHRS7C, CEP192, CPEB1/LINC00692, FBN2, ALDH1A1, and DMRTA1/FLJ35282 were associated with grades 2-4 OM. The addition of these SNPs increased sensitivity in detecting grades 2-4 OM cases to 52 %. CONCLUSIONS: These SNPs may be important for their roles in inflammatory pathways, epithelial healing, and chemotherapy detoxification.

An intermediate-risk multiple myeloma subgroup is defined by sIL-6r: levels synergistically increase with incidence of SNP rs2228145 and 1q21 amplification.

IL-6 signaling can be enhanced through transsignaling by the soluble IL-6 receptor (sIL-6r), allowing for the pleiotropic cytokine to affect cells it would not ordinarily have an effect on. Serum levels of sIL-6r can be used as an independent prognostic indicator and further stratify the GEP 70-gene low-risk group to identify an intermediate-risk group in multiple myeloma (MM). By analyzing more than 600 MM patients with ELISA, genotyping, and gene expression profiling tools, we show how the combination of 2 independent molecular genetic events is related to synergistic increases in sIL-6r levels. We also show that the rs2228145 minor allele is related to increased expression levels of an IL-6r splice variant that purportedly codes exclusively for a sIL-6r isoform. Together, the SNP rs2228145 minor allele C and amplification of chromosome 1q21 are significantly correlated to an increase in sIL-6r levels, which are associated with lower overall survival in 70-gene low-risk disease, and aid in identification of the intermediate-risk MM group.

Prediction of cytogenetic abnormalities with gene expression profiles.

Cytogenetic abnormalities are important clinical parameters in various types of cancer, including multiple myeloma. We developed a model to predict cytogenetic abnormalities in patients with multiple myeloma using gene expression profiling and validated it by different cytogenetic techniques. The model has an accuracy rate up to 0.89. These results provide proof of concept for the hypothesis that gene expression profiling is a superior genomic method for clinical molecular diagnosis and/or prognosis.

High-resolution genomic profiles define distinct clinico-pathogenetic subgroups of multiple myeloma patients.

To identify genetic events underlying the genesis and progression of multiple myeloma (MM), we conducted a high-resolution analysis of recurrent copy number alterations (CNAs) and expression profiles in a collection of MM cell lines and outcome-annotated clinical specimens. Attesting to the molecular heterogeneity of MM, unsupervised classification using nonnegative matrix factorization (NMF) designed for array comparative genomic hybridization (aCGH) analysis uncovered distinct genomic subtypes. Additionally, we defined 87 discrete minimal common regions (MCRs) within recurrent and highly focal CNAs. Further integration with expression data generated a refined list of MM gene candidates residing within these MCRs, thereby providing a genomic framework for dissection of disease pathogenesis, improved clinical management, and initiation of targeted drug discovery for specific MM patients.

Pharmacogenomics of bortezomib test-dosing identifies hyperexpression of proteasome genes, especially PSMD4, as novel high-risk feature in myeloma treated with Total Therapy 3.

Gene expression profiling (GEP) of purified plasma cells 48 hours after thalidomide and dexamethasone test doses showed these agents' mechanisms of action and provided prognostic information for untreated myeloma patients on Total Therapy 2 (TT2). Bortezomib was added in Total Therapy 3 (TT3), and 48 hours after bortezomib GEP analysis identified 80 highly survival-discriminatory genes in a training set of 142 TT3A patients that were validated in 128 patients receiving TT3B. The 80-gene GEP model (GEP80) also distinguished outcomes when applied at baseline in both TT3 and TT2 protocols. In context of our validated 70-gene model (GEP70), the GEP80 model identified 9% of patients with a grave prognosis among those with GEP70-defined low-risk disease and 41% of patients with favorable prognosis among those with GEP70-defined high-risk disease. PMSD4 was 1 of 3 genes common to both models. Residing on chromosome 1q21, PSMD4 expression is highly sensitive to copy number. Both higher PSMD4 expression levels and higher 1q21 copy numbers affected clinical outcome adversely. GEP80 baseline-defined high risk, high lactate dehydrogenase, and low albumin were the only independent adverse variables surviving multivariate survival model. We are investigating whether second-generation proteasome inhibitors (eg, carfilzomib) can overcome resistance associated with high PSMD4 levels.

High-risk myeloma is associated with global elevation of miRNAs and overexpression of EIF2C2/AGO2.

MicroRNAs (miRNAs) are noncoding RNAs that regulate global gene expression. miRNAs often act synergistically to repress target genes, and their dysregulation can contribute to the initiation and progression of a variety of cancers. The clinical relationship between global expression of miRNA and mRNA in cancer has not been studied in detail. We used whole-genome microarray analyses of CD138-enriched plasma cells from 52 newly diagnosed cases of multiple myeloma to correlate miRNA expression profiles with a validated mRNA-based risk stratification score, proliferation index, and predefined gene sets. In stark contrast to mRNAs, we discovered that all tested miRNAs were significantly up-regulated in high-risk disease as defined by a validated 70-gene risk score (P < 0.01) and proliferation index (P < 0.05). Increased expression of EIF2C2/AGO2, a master regulator of the maturation and function of miRNAs and a component of the 70-gene mRNA risk model, is driven by DNA copy number gains in MM. Silencing of AGO2 dramatically decreased viability in MM cell lines. Genome-wide elevated expression of miRNAs in high-risk MM may be secondary to deregulation of AGO2 and the enzyme complexes that regulate miRNA maturation and function.

Microdissection Methods Utilizing Single-Cell Subtype Analysis and the Impact on Precision Medicine.

Improving the utilization of tumor tissue from diagnostic biopsies is an unmet medical need. This is especially relevant today in the rapidly evolving precision oncology field where tumor genotyping is often essential for the indication of many advanced and targeted therapies. National Comprehensive Cancer Network (NCCN) guidelines now mandate molecular testing for clinically actionable targets in certain malignancies. Utilizing advanced stage lung cancer as an example, an improved genotyping approach for solid tumors is possible. The strategy involves optimization of the microdissection process and analysis of a large number of identical target cells from formalin-fixed paraffin-embedded (FFPE) specimens sharing similar characteristics, in other words, single-cell subtype analysis. The shared characteristics can include immunostaining status, cell phenotype, and/or spatial location within a histological section. Synergy between microdissection and droplet digital PCR (ddPCR) enhances the molecular analysis. We demonstrate here a methodology that illustrates genotyping of a solid tumor from a small tissue biopsy sample in a time- and cost-efficient manner, using immunostain targeting as an example.

Expansion of Human Papillomavirus-Specific T Cells in Periphery and Cervix in a Therapeutic Vaccine Recipient Whose Cervical High-Grade Squamous Intraepithelial Lesion Regressed.

Advances in high-throughput sequencing have revolutionized the manner with which we can study T cell responses. We describe a woman who received a human papillomavirus (HPV) therapeutic vaccine called PepCan, and experienced complete resolution of her cervical high-grade squamous intraepithelial lesion. By performing bulk T cell receptor (TCR) ß deep sequencing of peripheral blood mononuclear cells before and after 4 vaccinations, 70 putatively vaccine-specific clonotypes were identified for being significantly increased using a beta-binomial model. In order to verify the vaccine-specificity of these clonotypes, T cells with specificity to a region, HPV 16 E6 91-115, previously identified to be vaccine-induced using an interferon-γ enzyme-linked immunospot assay, were sorted and analyzed using single-cell RNA-seq and TCR sequencing. HPV specificity in 60 of the 70 clonotypes identified to be vaccine-specific was demonstrated. TCR ß bulk sequencing of the cervical liquid-based cytology samples and cervical formalin-fixed paraffin-embedded samples before and after 4 vaccinations demonstrated the presence of these HPV-specific T cells in the cervix. Combining traditional and cutting-edge immunomonitoring techniques enabled us to demonstrate expansion of HPV-antigen specific T cells not only in the periphery but also in the cervix. Such an approach should be useful as a novel approach to assess vaccine-specific responses in various anatomical areas.

Exogenous phosphatidic acid reduces acetaminophen-induced liver injury in mice by activating hepatic interleukin-6 signaling through inter-organ crosstalk.

We previously demonstrated that endogenous phosphatidic acid (PA) promotes liver regeneration after acetaminophen (APAP) hepatotoxicity. Here, we hypothesized that exogenous PA is also beneficial. To test that, we treated mice with a toxic APAP dose at 0 h, followed by PA or vehicle (Veh) post-treatment. We then collected blood and liver at 6, 24, and 52 h. Post-treatment with PA 2 h after APAP protected against liver injury at 6 h, and the combination of PA and N-acetyl-l-cysteine (NAC) reduced injury more than NAC alone. Interestingly, PA did not affect canonical mechanisms of APAP toxicity. Instead, transcriptomics revealed that PA activated interleukin-6 (IL-6) signaling in the liver. Consistent with that, serum IL-6 and hepatic signal transducer and activator of transcription 3 (Stat3) phosphorylation increased in PA-treated mice. Furthermore, PA failed to protect against APAP in IL-6-deficient animals. Interestingly, IL-6 expression increased 18-fold in adipose tissue after PA, indicating that adipose is a source of PA-induced circulating IL-6. Surprisingly, however, exogenous PA did not alter regeneration, despite the importance of endogenous PA in liver repair, possibly due to its short half-life. These data demonstrate that exogenous PA is also beneficial in APAP toxicity and reinforce the protective effects of IL-6 in this model.

PHF19 inhibition as a therapeutic target in multiple myeloma.

Epigenetic deregulation is increasingly recognized as a contributing pathological factor in multiple myeloma (MM). In particular tri-methylation of H3 lysine 27 (H3K27me3), which is catalyzed by PHD finger protein 19 (PHF19), a subunit of the Polycomb Repressive Complex 2 (PRC2), has recently shown to be a crucial mediator of MM tumorigenicity. Overexpression of PHF19 in MM has been associated with worse clinical outcome. Yet, while there is mounting evidence that PHF19 overexpression plays a crucial role in MM carcinogenesis downstream mechanisms remain to be elucidated. In the current study we use a functional knock down (KD) of PHF19 to investigate the biological role of PHF19 and show that PHF19KD leads to decreased tumor growth in vitro and in vivo. Expression of major cancer players such as bcl2, myc and EGR1 were decreased upon PHF19KD further underscoring the role of PHF19 in MM biology. Additionally, our results highlighted the prognostic impact of PHF19 overexpression, which was significantly associated with worse survival. Overall, our study underscores the premise that targeting the PHF19-PRC2 complex would open up avenues for novel MM therapies.

Myeloma-derived Dickkopf-1 disrupts Wnt-regulated osteoprotegerin and RANKL production by osteoblasts: a potential mechanism underlying osteolytic bone lesions in multiple myeloma.

Multiple myeloma (MM) is characterized by osteolytic bone lesions (OBL) that arise as a consequence of osteoblast inactivation and osteoclast activation adjacent to tumor foci within bone. Wnt signaling in osteoblasts regulates osteoclastogenesis through the differential activation and inactivation of Receptor Activator of Nuclear factor Kappa B Ligand (RANKL) and osteoprotegerin (OPG), positive and negative regulators of osteoclast differentiation, respectively. We demonstrate here that MM cell-derived DKK1, a soluble inhibitor of canonical Wnt signaling, disrupted Wnt3a-regulated OPG and RANKL expression in osteoblasts. Confirmed in multiple independent assays, we show that pretreatment with rDKK1 completely abolished Wnt3a-induced OPG mRNA and protein production by mouse and human osteoblasts. In addition, we show that Wnt3a-induced OPG expression was diminished in osteoblasts cocultured with a DKK1-expressing MM cell line or primary MM cells. Finally, we show that bone marrow sera from 21 MM patients significantly suppressed Wnt3a-induced OPG expression and enhanced RANKL expression in osteoblasts in a DKK1-dependent manner. These results suggest that DKK1 may play a key role in the development of MM-associated OBL by directly interrupting Wnt-regulated differentiation of osteoblasts and indirectly increasing osteoclastogenesis via a DKK1-mediated increase in RANKL-to-OPG ratios.

A meta-analysis of genome-wide association studies of multiple myeloma among men and women of African ancestry.

Persons of African ancestry (AA) have a twofold higher risk for multiple myeloma (MM) compared with persons of European ancestry (EA). Genome-wide association studies (GWASs) support a genetic contribution to MM etiology in individuals of EA. Little is known about genetic risk factors for MM in individuals of AA. We performed a meta-analysis of 2 GWASs of MM in 1813 cases and 8871 controls and conducted an admixture mapping scan to identify risk alleles. We fine-mapped the 23 known susceptibility loci to find markers that could better capture MM risk in individuals of AA and constructed a polygenic risk score (PRS) to assess the aggregated effect of known MM risk alleles. In GWAS meta-analysis, we identified 2 suggestive novel loci located at 9p24.3 and 9p13.1 at P < 1 × 10-6; however, no genome-wide significant association was noted. In admixture mapping, we observed a genome-wide significant inverse association between local AA at 2p24.1-23.1 and MM risk in AA individuals. Of the 23 known EA risk variants, 20 showed directional consistency, and 9 replicated at P < .05 in AA individuals. In 8 regions, we identified markers that better capture MM risk in persons with AA. AA individuals with a PRS in the top 10% had a 1.82-fold (95% confidence interval, 1.56-2.11) increased MM risk compared with those with average risk (25%-75%). The strongest functional association was between the risk allele for variant rs56219066 at 5q15 and lower ELL2 expression (P = 5.1 × 10-12). Our study shows that common genetic variation contributes to MM risk in individuals with AA.

Ellipticine derivative NSC 338258 represents a potential new antineoplastic agent for the treatment of multiple myeloma.

High-risk multiple myeloma can be correlated with amplification and overexpression of the cell cycle regulator CKS1B. Herein, we used the COMPARE algorithm to correlate high expression of CKS1B mRNA in the NCI-60 cell line panel with the concentration causing 50% growth inhibition (GI(50)) of >40,000 synthetic compounds. This led to the identification of NSC 338258 (EPED3), a highly stable, hydrophilic derivative of the plant alkaloid ellipticine. In vitro, this synthetic anticancer compound exhibits dramatic cytotoxic activity against myeloma cells grown in suspension or in coculture with stromal cells. EPED3-induced cell cycle arrest and an apoptotic progression that appear to be a consequence of the instantaneous effect of the drug on cytoplasmic organelles, particularly mitochondria. Disruption of mitochondria and cytoplasmic distribution of cytochrome c initiated the intracellular proteolytic cascade through the intrinsic apoptotic pathway. EPED3 is able to induce apoptosis in myeloma cells with de novo or acquired resistance to commonly administered antimyeloma agents. Collectively, our data suggest that EPED3 targets mitochondrial function to rapidly deplete chemical energy and initiate apoptosis in myeloma cells at nanomolar concentrations while leaving stromal cells unharmed.

Author Correction: Identification of multiple risk loci and regulatory mechanisms influencing susceptibility to multiple myeloma.

The original version of this Article contained an error in the spelling of a member of the PRACTICAL Consortium, Manuela Gago-Dominguez, which was incorrectly given as Manuela Gago Dominguez. This has now been corrected in both the PDF and HTML versions of the Article. Furthermore, in the original HTML version of this Article, the order of authors within the author list was incorrect. The PRACTICAL consortium was incorrectly listed after Richard S. Houlston and should have been listed after Nora Pashayan. This error has been corrected in the HTML version of the Article; the PDF version was correct at the time of publication.

Dickkopf homolog 1 mediates endothelin-1-stimulated new bone formation.

Tumor-produced endothelin-1 (ET-1) stimulates osteoblasts to form new bone and is an important mediator of osteoblastic bone metastasis. The anabolic actions of ET-1 in osteoblasts were investigated by gene microarray analyses of murine neonatal calvarial organ cultures. Targets of ET-1 action were validated by real-time RT-PCR in murine primary osteoblast cultures. IL-6, IL-11, the CCN (CYR61, CTGF, NOV) family members cysteine-rich protein 61 and connective tissue growth factor, inhibin beta-A, serum/glucocorticoid regulated kinase, receptor activator of nuclear factor kappaB ligand, snail homolog 1, tissue inhibitor of metalloproteinase 3, and TG-interacting factor transcripts were increased by ET-1. ET-1 decreased the transcript for the Wnt signaling pathway inhibitor, dickkopf homolog 1 (Dkk1). Calvarial organ cultures treated with ET-1 had lower concentrations of DKK1 protein in conditioned media than control cultures. High DKK1 concentrations in bone marrow suppress bone formation in multiple myeloma. We hypothesized that the converse occurs in osteoblastic bone metastasis, where ET-1 stimulates osteoblast activity by reducing autocrine production of DKK1. Recombinant DKK1 blocked ET-1-mediated osteoblast proliferation and new bone formation in calvarial organ cultures, whereas a DKK1-neutralizing antibody increased osteoblast numbers and new bone formation. ET-1 directed nuclear translocation of beta-catenin in osteoblasts, indicating activation of the Wnt signaling pathway. The data suggest that ET-1 increases osteoblast proliferation and new bone formation by activating the Wnt signaling pathway through suppression of the Wnt pathway inhibitor DKK1.