Deep MRD profiling defines outcome and unveils different modes of treatment resistance in standard- and high-risk myeloma.

Patients with multiple myeloma (MM) carrying standard- or high-risk cytogenetic abnormalities (CAs) achieve similar complete response (CR) rates, but the later have inferior progression-free survival (PFS). This questions the legitimacy of CR as a treatment endpoint and represents a biological conundrum regarding the nature of tumor reservoirs that persist after therapy in high-risk MM. We used next-generation flow (NGF) cytometry to evaluate measurable residual disease (MRD) in MM patients with standard- vs high-risk CAs (n = 300 and 90, respectively) enrolled in the PETHEMA/GEM2012MENOS65 trial, and to identify mechanisms that determine MRD resistance in both patient subgroups (n = 40). The 36-month PFS rates were higher than 90% in patients with standard- or high-risk CAs achieving undetectable MRD. Persistent MRD resulted in a median PFS of âˆ¼3 and 2 years in patients with standard- and high-risk CAs, respectively. Further use of NGF to isolate MRD, followed by whole-exome sequencing of paired diagnostic and MRD tumor cells, revealed greater clonal selection in patients with standard-risk CAs, higher genomic instability with acquisition of new mutations in high-risk MM, and no unifying genetic event driving MRD resistance. Conversely, RNA sequencing of diagnostic and MRD tumor cells uncovered the selection of MRD clones with singular transcriptional programs and reactive oxygen species-mediated MRD resistance in high-risk MM. Our study supports undetectable MRD as a treatment endpoint for patients with MM who have high-risk CAs and proposes characterizing MRD clones to understand and overcome MRD resistance. This trial is registered at www.clinicaltrials.gov as #NCT01916252.

Tumor cells in light-chain amyloidosis and myeloma show distinct transcriptional rewiring of normal plasma cell development.

Although light-chain amyloidosis (AL) and multiple myeloma (MM) are characterized by tumor plasma cell (PC) expansion in bone marrow (BM), their clinical presentation differs. Previous attempts to identify unique pathogenic mechanisms behind such differences were unsuccessful, and no studies have investigated the differentiation stage of tumor PCs in patients with AL and MM. We sought to define a transcriptional atlas of normal PC development in secondary lymphoid organs (SLOs), peripheral blood (PB), and BM for comparison with the transcriptional programs (TPs) of tumor PCs in AL, MM, and monoclonal gammopathy of undetermined significance (MGUS). Based on bulk and single-cell RNA sequencing, we observed 13 TPs during transition of normal PCs throughout SLOs, PB, and BM. We further noted the following: CD39 outperforms CD19 to discriminate newborn from long-lived BM-PCs; tumor PCs expressed the most advantageous TPs of normal PC differentiation; AL shares greater similarity to SLO-PCs whereas MM is transcriptionally closer to PB-PCs and newborn BM-PCs; patients with AL and MM enriched in immature TPs had inferior survival; and protein N-linked glycosylation-related TPs are upregulated in AL. Collectively, we provide a novel resource to understand normal PC development and the transcriptional reorganization of AL and other monoclonal gammopathies.

Biological and clinical significance of dysplastic hematopoiesis in patients with newly diagnosed multiple myeloma.

Risk of developing myelodysplastic syndrome (MDS) is significantly increased in both multiple myeloma (MM) and monoclonal gammopathy of undetermined significance, suggesting that it is therapy independent. However, the incidence and sequelae of dysplastic hematopoiesis at diagnosis are unknown. Here, we used multidimensional flow cytometry (MFC) to prospectively screen for the presence of MDS-associated phenotypic alterations (MDS-PA) in the bone marrow of 285 patients with MM enrolled in the PETHEMA/GEM2012MENOS65 trial (#NCT01916252). We investigated the clinical significance of monocytic MDS-PA in a larger series of 1252 patients enrolled in 4 PETHEMA/GEM protocols. At diagnosis, 33 (11.6%) of 285 cases displayed MDS-PA. Bulk and single-cell-targeted sequencing of MDS recurrently mutated genes in CD34+ progenitors (and dysplastic lineages) from 67 patients revealed clonal hematopoiesis in 13 (50%) of 26 cases with MDS-PA vs 9 (22%) of 41 without MDS-PA; TET2 and NRAS were the most frequently mutated genes. Dynamics of MDS-PA at diagnosis and after autologous transplant were evaluated in 86 of 285 patients and showed that in most cases (69 of 86 [80%]), MDS-PA either persisted or remained absent in patients with or without MDS-PA at diagnosis, respectively. Noteworthy, MDS-associated mutations infrequently emerged after high-dose therapy. Based on MFC profiling, patients with MDS-PA have altered hematopoiesis and T regulatory cell distribution in the tumor microenvironment. Importantly, the presence of monocytic MDS-PA at diagnosis anticipated greater risk of hematologic toxicity and was independently associated with inferior progression-free survival (hazard ratio, 1.5; P = .02) and overall survival (hazard ratio, 1.7; P = .01). This study reveals the biological and clinical significance of dysplastic hematopoiesis in newly diagnosed MM, which can be screened with moderate sensitivity using cost-effective MFC.

Immunogenomic identification and characterization of granulocytic myeloid-derived suppressor cells in multiple myeloma.

Granulocytic myeloid-derived suppressor cells (G-MDSCs) promote tumor growth and immunosuppression in multiple myeloma (MM). However, their phenotype is not well established for accurate monitoring or clinical translation. We aimed to provide the phenotypic profile of G-MDSCs based on their prognostic significance in MM, immunosuppressive potential, and molecular program. The preestablished phenotype of G-MDSCs was evaluated in bone marrow samples from controls and MM patients using multidimensional flow cytometry; surprisingly, we found that CD11b+CD14-CD15+CD33+HLADR- cells overlapped with common eosinophils and neutrophils, which were not expanded in MM patients. Therefore, we relied on automated clustering to unbiasedly identify all granulocytic subsets in the tumor microenvironment: basophils, eosinophils, and immature, intermediate, and mature neutrophils. In a series of 267 newly diagnosed MM patients (GEM2012MENOS65 trial), only the frequency of mature neutrophils at diagnosis was significantly associated with patient outcome, and a high mature neutrophil/T-cell ratio resulted in inferior progression-free survival (P < .001). Upon fluorescence-activated cell sorting of each neutrophil subset, T-cell proliferation decreased in the presence of mature neutrophils (0.5-fold; P = .016), and the cytotoxic potential of T cells engaged by a BCMA×CD3-bispecific antibody increased notably with the depletion of mature neutrophils (fourfold; P = .0007). Most interestingly, RNA sequencing of the 3 subsets revealed that G-MDSC-related genes were specifically upregulated in mature neutrophils from MM patients vs controls because of differential chromatin accessibility. Taken together, our results establish a correlation between the clinical significance, immunosuppressive potential, and transcriptional network of well-defined neutrophil subsets, providing for the first time a set of optimal markers (CD11b/CD13/CD16) for accurate monitoring of G-MDSCs in MM.

Identification of Recurrent Mutations in the microRNA-Binding Sites of B-Cell Lymphoma-Associated Genes in Follicular Lymphoma.

Follicular lymphoma (FL) is a common indolent B-cell lymphoma that can transform into the more aggressive transformed FL (tFL). However, the molecular process driving this transformation is uncertain. In this work, we aimed to identify microRNA (miRNA)-binding sites recurrently mutated in follicular lymphoma patients, as well as in transformed FL patients. Using whole-genome sequencing data from FL tumors, we discovered 544 mutations located in bioinformatically predicted microRNA-binding sites. We then studied these specific regions using targeted sequencing in a cohort of 55 FL patients, found 16 recurrent mutations, and identified a further 69 variants. After filtering for QC, we identified 21 genes with mutated miRNA-binding sites that were also enriched for B-cell-associated genes by Gene Ontology. Over 40% of mutations identified in these genes were present exclusively in tFL patients. We validated the predicted miRNA-binding sites of five of the genes by luciferase assay and demonstrated that the identified mutations in BCL2 and EZH2 genes impaired the binding efficiency of miR-5008 and miR-144 and regulated the endogenous levels of messenger RNA (mRNA).

New Concepts in Cancer Biomarkers: Circulating miRNAs in Liquid Biopsies.

The effective and efficient management of cancer patients relies upon early diagnosis and/or the monitoring of treatment, something that is often difficult to achieve using standard tissue biopsy techniques. Biological fluids such as blood hold great possibilities as a source of non-invasive cancer biomarkers that can act as surrogate markers to biopsy-based sampling. The non-invasive nature of these "liquid biopsies" ultimately means that cancer detection may be earlier and that the ability to monitor disease progression and/or treatment response represents a paradigm shift in the treatment of cancer patients. Below, we review one of the most promising classes of circulating cancer biomarkers: microRNAs (miRNAs). In particular, we will consider their history, the controversy surrounding their origin and biology, and, most importantly, the hurdles that remain to be overcome if they are really to become part of future clinical practice.

The circulating transcriptome as a source of non-invasive cancer biomarkers: concepts and controversies of non-coding and coding RNA in body fluids.

The gold standard for cancer diagnosis remains the histological examination of affected tissue, obtained either by surgical excision, or radiologically guided biopsy. Such procedures however are expensive, not without risk to the patient, and require consistent evaluation by expert pathologists. Consequently, the search for non-invasive tools for the diagnosis and management of cancer has led to great interest in the field of circulating nucleic acids in plasma and serum. An additional benefit of blood-based testing is the ability to carry out screening and repeat sampling on patients undergoing therapy, or monitoring disease progression allowing for the development of a personalized approach to cancer patient management. Despite having been discovered over 60 years ago, the clear clinical potential of circulating nucleic acids, with the notable exception of prenatal diagnostic testing, has yet to translate into the clinic. The recent discovery of non-coding (nc) RNA (in particular micro(mi)RNAs) in the blood has provided fresh impetuous for the field. In this review, we discuss the potential of the circulating transcriptome (coding and ncRNA), as novel cancer biomarkers, the controversy surrounding their origin and biology, and most importantly the hurdles that remain to be overcome if they are really to become part of future clinical practice.

Targeted next-generation sequencing and non-coding RNA expression analysis of clear cell papillary renal cell carcinoma suggests distinct pathological mechanisms from other renal tumour subtypes.

Clear cell tubulopapillary renal cell carcinoma (CCPRCC) is a recently described rare renal malignancy that displays characteristic gross, microscopic and immunohistochemical differences from other renal tumour types. However, CCPRCC remains a very poorly understood entity. We therefore sought to elucidate some of the molecular mechanisms involved in this neoplasm by carrying out targeted next-generation sequencing (NGS) to identify associated mutations, and in addition examined the expression of non-coding (nc) RNAs. We identified multiple somatic mutations in CCPRCC cases, including a recurrent [3/14 cases (21%)] non-synonymous T992I mutation in the MET proto-oncogene, a gene associated with epithelial-to-mesenchymal transition (EMT). Using a microarray approach, we found that the expression of mature (n = 1105) and pre-miRNAs (n = 1105), as well as snoRNA and scaRNAs (n = 2214), in CCPRCC cases differed from that of clear cell renal cell carcinoma (CCRCC) or papillary renal cell carcinoma (PRCC) tumours. Surprisingly, and unlike other renal tumour subtypes, we found that all five members of the miR-200 family were over-expressed in CCPRCC cases. As these miRNAs are intimately involved with EMT, we stained CCPRCC cases for E-cadherin, vimentin and ß-catenin and found that the tumour cells of all cases were positive for all three markers, a combination rarely reported in other renal tumours that could have diagnostic implications. Taken together with the mutational analysis, these data suggest that EMT in CCPRCC tumour cells is incomplete or blocked, consistent with the indolent clinical course typical of this malignancy. In summary, as well as describing a novel pathological mechanism in renal carcinomas, this study adds to the mounting evidence that CCPRCC should be formally considered a distinct entity. Microarray data have been deposited in the GEO database [GEO accession number (GSE51554)].

Differential chemosensitivity to antifolate drugs between RAS and BRAF melanoma cells.

BACKGROUND: The importance of the genetic background of cancer cells for the individual susceptibility to cancer treatments is increasingly apparent. In melanoma, the existence of a BRAF mutation is a main predictor for successful BRAF-targeted therapy. However, despite initial successes with these therapies, patients relapse within a year and have to move on to other therapies. Moreover, patients harbouring a wild type BRAF gene (including 25% with NRAS mutations) still require alternative treatment such as chemotherapy. Multiple genetic parameters have been associated with response to chemotherapy, but despite their high frequency in melanoma nothing is known about the impact of BRAF or NRAS mutations on the response to chemotherapeutic agents. METHODS: Using cell proliferation and DNA methylation assays, FACS analysis and quantitative-RT-PCR we have characterised the response of a panel of NRAS and BRAF mutant melanoma cell lines to various chemotherapy drugs, amongst them dacarbazine (DTIC) and temozolomide (TMZ) and DNA synthesis inhibitors. RESULTS: Although both, DTIC and TMZ act as alkylating agents through the same intermediate, NRAS and BRAF mutant cells responded differentially only to DTIC. Further analysis revealed that the growth-inhibitory effects mediated by DTIC were rather due to interference with nucleotide salvaging, and that NRAS mutant melanoma cells exhibit higher activity of the nucleotide synthesis enzymes IMPDH and TK1. Importantly, the enhanced ability of RAS mutant cells to use nucleotide salvaging resulted in resistance to DHFR inhibitors. CONCLUSION: In summary, our data suggest that the genetic background in melanoma cells influences the response to inhibitors blocking de novo DNA synthesis, and that defining the RAS mutation status could be used to stratify patients for the use of antifolate drugs.

Characterization of the Common Genetic Variation in the Spanish Population of Navarre.

Large-scale genomic studies have significantly increased our knowledge of genetic variability across populations. Regional genetic profiling is essential for distinguishing common benign variants from disease-causing ones. To this end, we conducted a comprehensive characterization of exonic variants in the population of Navarre (Spain), utilizing whole genome sequencing data from 358 unrelated individuals of Spanish origin. Our analysis revealed 61,410 biallelic single nucleotide variants (SNV) within the Navarrese cohort, with 35% classified as common (MAF > 1%). By comparing allele frequency data from 1000 Genome Project (excluding the Iberian cohort of Spain, IBS), Genome Aggregation Database, and a Spanish cohort (including IBS individuals and data from Medical Genome Project), we identified 1069 SNVs common in Navarre but rare (MAF ≤ 1%) in all other populations. We further corroborated this observation with a second regional cohort of 239 unrelated exomes, which confirmed 676 of the 1069 SNVs as common in Navarre. In conclusion, this study highlights the importance of population-specific characterization of genetic variation to improve allele frequency filtering in sequencing data analysis to identify disease-causing variants.

Preclinical models for prediction of immunotherapy outcomes and immune evasion mechanisms in genetically heterogeneous multiple myeloma.

The historical lack of preclinical models reflecting the genetic heterogeneity of multiple myeloma (MM) hampers the advance of therapeutic discoveries. To circumvent this limitation, we screened mice engineered to carry eight MM lesions (NF-κB, KRAS, MYC, TP53, BCL2, cyclin D1, MMSET/NSD2 and c-MAF) combinatorially activated in B lymphocytes following T cell-driven immunization. Fifteen genetically diverse models developed bone marrow (BM) tumors fulfilling MM pathogenesis. Integrative analyses of ∼500 mice and ∼1,000 patients revealed a common MAPK-MYC genetic pathway that accelerated time to progression from precursor states across genetically heterogeneous MM. MYC-dependent time to progression conditioned immune evasion mechanisms that remodeled the BM microenvironment differently. Rapid MYC-driven progressors exhibited a high number of activated/exhausted CD8+ T cells with reduced immunosuppressive regulatory T (Treg) cells, while late MYC acquisition in slow progressors was associated with lower CD8+ T cell infiltration and more abundant Treg cells. Single-cell transcriptomics and functional assays defined a high ratio of CD8+ T cells versus Treg cells as a predictor of response to immune checkpoint blockade (ICB). In clinical series, high CD8+ T/Treg cell ratios underlie early progression in untreated smoldering MM, and correlated with early relapse in newly diagnosed patients with MM under Len/Dex therapy. In ICB-refractory MM models, increasing CD8+ T cell cytotoxicity or depleting Treg cells reversed immunotherapy resistance and yielded prolonged MM control. Our experimental models enable the correlation of MM genetic and immunological traits with preclinical therapy responses, which may inform the next-generation immunotherapy trials.

Large T cell clones expressing immune checkpoints increase during multiple myeloma evolution and predict treatment resistance.

Tumor recognition by T cells is essential for antitumor immunity. A comprehensive characterization of T cell diversity may be key to understanding the success of immunomodulatory drugs and failure of PD-1 blockade in tumors such as multiple myeloma (MM). Here, we use single-cell RNA and T cell receptor sequencing to characterize bone marrow T cells from healthy adults (n = 4) and patients with precursor (n = 8) and full-blown MM (n = 10). Large T cell clones from patients with MM expressed multiple immune checkpoints, suggesting a potentially dysfunctional phenotype. Dual targeting of PD-1 + LAG3 or PD-1 + TIGIT partially restored their function in mice with MM. We identify phenotypic hallmarks of large intratumoral T cell clones, and demonstrate that the CD27- and CD27+ T cell ratio, measured by flow cytometry, may serve as a surrogate of clonal T cell expansions and an independent prognostic factor in 543 patients with MM treated with lenalidomide-based treatment combinations.

An N-glycosylation hotspot in immunoglobulin κ light chains is associated with AL amyloidosis.

Immunoglobulin light chain (AL) amyloidosis is caused by a small, minimally proliferating B-cell/plasma-cell clone secreting a patient-unique, aggregation-prone, toxic light chain (LC). The pathogenicity of LCs is encrypted in their sequence, yet molecular determinants of amyloidogenesis are poorly understood. Higher rates of N-glycosylation among clonal κ LCs from patients with AL amyloidosis compared to other monoclonal gammopathies indicate that this post-translational modification is associated with a higher risk of developing AL amyloidosis. Here, we exploited LC sequence information from previously published amyloidogenic and control clonal LCs and from a series of 220 patients with AL amyloidosis or multiple myeloma followed at our Institutions to define sequence and spatial features of N-glycosylation, combining bioinformatics, biochemical, proteomics, structural and genetic analyses. We found peculiar sequence and spatial pattern of N-glycosylation in amyloidogenic κ LCs, with most of the N-glycosylation sites laying in the framework region 3, particularly within the E strand, and consisting mainly of the NFT sequon, setting them apart with respect to non-amyloidogenic clonal LCs. Our data further support a potential role of N-glycosylation in determining the pathogenic behavior of a subset of amyloidogenic LCs and may help refine current N-glycosylation-based prognostic assessments for patients with monoclonal gammopathies.

A Machine Learning Model Based on Tumor and Immune Biomarkers to Predict Undetectable MRD and Survival Outcomes in Multiple Myeloma.

PURPOSE: Undetectable measurable residual disease (MRD) is a surrogate of prolonged survival in multiple myeloma. Thus, treatment individualization based on the probability of a patient achieving undetectable MRD with a singular regimen could represent a new concept toward personalized treatment, with fast assessment of its success. This has never been investigated; therefore, we sought to define a machine learning model to predict undetectable MRD at the onset of multiple myeloma. EXPERIMENTAL DESIGN: This study included 487 newly diagnosed patients with multiple myeloma. The training (n = 152) and internal validation cohorts (n = 149) consisted of 301 transplant-eligible patients with active multiple myeloma enrolled in the GEM2012MENOS65 trial. Two external validation cohorts were defined by 76 high-risk transplant-eligible patients with smoldering multiple myeloma enrolled in the Grupo Español de Mieloma(GEM)-CESAR trial, and 110 transplant-ineligible elderly patients enrolled in the GEM-CLARIDEX trial. RESULTS: The most effective model to predict MRD status resulted from integrating cytogenetic [t(4;14) and/or del(17p13)], tumor burden (bone marrow plasma cell clonality and circulating tumor cells), and immune-related biomarkers. Accurate predictions of MRD outcomes were achieved in 71% of cases in the GEM2012MENOS65 trial (n = 214/301) and 72% in the external validation cohorts (n = 134/186). The model also predicted sustained MRD negativity from consolidation onto 2 years maintenance (GEM2014MAIN). High-confidence prediction of undetectable MRD at diagnosis identified a subgroup of patients with active multiple myeloma with 80% and 93% progression-free and overall survival rates at 5 years. CONCLUSIONS: It is possible to accurately predict MRD outcomes using an integrative, weighted model defined by machine learning algorithms. This is a new concept toward individualized treatment in multiple myeloma. See related commentary by Pawlyn and Davies, p. 2482.

Preneoplastic somatic mutations including MYD88L265P in lymphoplasmacytic lymphoma.

Normal cell counterparts of solid and myeloid tumors accumulate mutations years before disease onset; whether this occurs in B lymphocytes before lymphoma remains uncertain. We sequenced multiple stages of the B lineage in elderly individuals and patients with lymphoplasmacytic lymphoma, a singular disease for studying lymphomagenesis because of the high prevalence of mutated MYD88. We observed similar accumulation of random mutations in B lineages from both cohorts and unexpectedly found MYD88L265P in normal precursor and mature B lymphocytes from patients with lymphoma. We uncovered genetic and transcriptional pathways driving malignant transformation and leveraged these to model lymphoplasmacytic lymphoma in mice, based on mutated MYD88 in B cell precursors and BCL2 overexpression. Thus, MYD88L265P is a preneoplastic event, which challenges the current understanding of lymphomagenesis and may have implications for early detection of B cell lymphomas.

Characterization of herpes simplex virus 1 strains as platforms for the development of oncolytic viruses against liver cancer.

BACKGROUND: Diverse oncolytic viruses (OV) are being designed for the treatment of cancer. The characteristics of the parental virus strains may influence the properties of these agents. AIMS: To characterize two herpes simplex virus 1 strains (HSV-1 17syn(+) and HFEM) as platforms for virotherapy against liver cancer. METHODS: The luciferase reporter gene was introduced in the intergenic region 20 locus of both HSV-1 strains, giving rise to the Cgal-Luc and H6-Luc viruses. Their properties were studied in hepatocellular carcinoma (HCC) cells in vitro. Biodistribution was monitored by bioluminescence imaging (BLI) in athymic mice and immune-competent Balb/c mice. Immunogenicity was studied by MHC-tetramer staining, in vivo killing assays and determination of specific antibody production. Intratumoural transgene expression and oncolytic effect were studied in HuH-7 xenografts. RESULTS: The H6-Luc virus displayed a syncytial phenotype and showed higher cytolytic effect on some HCC cells. Upon intravenous or intrahepatic injection in mice, both viruses showed a transient transduction of the liver with rapid relocalization of bioluminescence in adrenal glands, spinal cord, uterus and ovaries. No significant differences were observed in the immunogenicity of these viruses. Local intratumoural administration caused progressive increase in transgene expression during the first 5 days and persisted for at least 2 weeks. H6-Luc achieved faster amplification of transgene expression and stronger inhibition of tumour growth than Cgal-Luc, although toxicity of these non-attenuated viruses should be reduced to obtain a therapeutic effect. CONCLUSIONS: The syncytial H6-Luc virus has a strong oncolytic potential on human HCC xenografts and could be the basis for potent OV.

The Urinary Transcriptome as a Source of Biomarkers for Prostate Cancer.

Prostate cancer (PCa) is the second most common cancer of men and is typically slow-growing and asymptomatic. The use of blood PSA as a screening method has greatly improved PCa diagnosis, but high levels of false positives has raised much interest in alternative biomarkers. We used next-generation sequencing (NGS) to elucidate the urinary transcriptome of whole urine collected from high-stage and low-stage PCa patients as well as from patients with the confounding diagnosis of benign hyperplasia (BPH). We identified and validated five differentially expressed protein-coding genes (FTH1 BRPF1, OSBP, PHC3, and UACA) in an independent validation cohort of small-volume (1 mL) centrifuged urine (n = 94) and non-centrifuged urine (n = 84) by droplet digital (dd)PCR. These biomarkers were able to discriminate between BPH and PCa patients and healthy controls using either centrifuged or non-centrifuged whole urine samples, suggesting that the urinary transcriptome is a valuable source of non-invasive biomarkers for PCa that warrants further investigation.

Circulating tumor cells for comprehensive and multiregional non-invasive genetic characterization of multiple myeloma.

Multiple myeloma (MM) patients undergo repetitive bone marrow (BM) aspirates for genetic characterization. Circulating tumor cells (CTCs) are detectable in peripheral blood (PB) of virtually all MM cases and are prognostic, but their applicability for noninvasive screening has been poorly investigated. Here, we used next-generation flow (NGF) cytometry to isolate matched CTCs and BM tumor cells from 53 patients and compared their genetic profile. In eight cases, tumor cells from extramedullary (EM) plasmacytomas were also sorted and whole-exome sequencing was performed in the three spatially distributed tumor samples. CTCs were detectable by NGF in the PB of all patients with MM. Based on the cancer cell fraction of clonal and subclonal mutations, we found that ~22% of CTCs egressed from a BM (or EM) site distant from the matched BM aspirate. Concordance between BM tumor cells and CTCs was high for chromosome arm-level copy number alterations (≥95%) though not for translocations (39%). All high-risk genetic abnormalities except one t(4;14) were detected in CTCs whenever present in BM tumor cells. Noteworthy, ≥82% mutations present in BM and EM clones were detectable in CTCs. Altogether, these results support CTCs for noninvasive risk-stratification of MM patients based on their numbers and genetic profile.

Transcriptional profiling of circulating tumor cells in multiple myeloma: a new model to understand disease dissemination.

The reason why a few myeloma cells egress from the bone marrow (BM) into peripheral blood (PB) remains unknown. Here, we investigated molecular hallmarks of circulating tumor cells (CTCs) to identify the events leading to myeloma trafficking into the bloodstream. After using next-generation flow to isolate matched CTCs and BM tumor cells from 32 patients, we found high correlation in gene expression at single-cell and bulk levels (r ≥ 0.94, P = 10-16), with only 55 genes differentially expressed between CTCs and BM tumor cells. CTCs overexpressed genes involved in inflammation, hypoxia, or epithelial-mesenchymal transition, whereas genes related with proliferation were downregulated in CTCs. The cancer stem cell marker CD44 was overexpressed in CTCs, and its knockdown significantly reduced migration of MM cells towards SDF1-α and their adhesion to fibronectin. Approximately half (29/55) of genes differentially expressed in CTCs were prognostic in patients with newly-diagnosed myeloma (n = 553; CoMMpass). In a multivariate analysis including the R-ISS, overexpression of CENPF and LGALS1 was significantly associated with inferior survival. Altogether, these results help understanding the presence of CTCs in PB and suggest that hypoxic BM niches together with a pro-inflammatory microenvironment induce an arrest in proliferation, forcing tumor cells to circulate in PB and seek other BM niches to continue growing.

The Circulating Transcriptome as a Source of Biomarkers for Melanoma.

The circulating transcriptome is a valuable source of cancer biomarkers, which, with the exception of microRNAs (miRNAs), remains relatively unexplored. To elucidate which RNAs are present in plasma from melanoma patients and which could be used to distinguish cancer patients from healthy individuals, we used next generation sequencing (NGS), and validation was carried out by qPCR and/or ddPCR. We identified 442 different microRNAs in samples, eleven of which were differentially expressed (p < 0.05). Levels of miR-134-5p and miR-320a-3p were significantly down-regulated (p < 0.001) in melanoma samples (n = 96) compared to healthy controls (n = 28). Differentially expressed protein-encoding mRNA 5'-fragments were enriched for the angiopoietin, p21-activated kinase (PAK), and EIF2 pathways. Levels of ATM1, AMFR, SOS1, and CD109 gene fragments were up-regulated (p < 0.001) in melanoma samples (n = 144) compared to healthy controls (n = 41) (AUC = 0.825). Over 40% of mapped reads were YRNAs, a class of non-coding RNAs that to date has been little explored. Expression levels of RNY3P1, RNY4P1, and RNY4P25 were significantly higher in patients with stage 0 disease than either healthy controls or more advanced stage disease (p < 0.001). In conclusion, we have identified a number of novel RNA biomarkers, which, most importantly, we validated in multi-center retrospective and prospective cohorts, suggesting potential diagnostic use of these RNA species.