RNA-seq is currently considered the most powerful, robust and adaptable technique for measuring gene expression and transcription activation at genome-wide level. As the analysis of RNA-seq data is complex, it has prompted a large amount of research on algorithms and methods. This has resulted in a substantial increase in the number of options available at each step of the analysis. Consequently, there is no clear consensus about the most appropriate algorithms and pipelines that should be used to analyse RNA-seq data. In the present study, 192 pipelines using alternative methods were applied to 18 samples from two human cell lines and the performance of the results was evaluated. Raw gene expression signal was quantified by non-parametric statistics to measure precision and accuracy. Differential gene expression performance was estimated by testing 17 differential expression methods. The procedures were validated by qRT-PCR in the same samples. This study weighs up the advantages and disadvantages of the tested algorithms and pipelines providing a comprehensive guide to the different methods and procedures applied to the analysis of RNA-seq data, both for the quantification of the raw expression signal and for the differential gene expression.

Circulating myeloma tumor cells (CTCs) as defined by the presence of peripheral blood (PB) clonal plasma cells (PCs) are a powerful prognostic marker in multiple myeloma (MM). However, the biological features of CTCs and their pathophysiological role in MM remains unexplored. Here, we investigate the phenotypic, cytogenetic, and functional characteristics as well as the circadian distribution of CTCs vs paired bone marrow (BM) clonal PCs from MM patients. Our results show that CTCs typically represent a unique subpopulation of all BM clonal PCs, characterized by downregulation (P < .05) of integrins (CD11a/CD11c/CD29/CD49d/CD49e), adhesion (CD33/CD56/CD117/CD138), and activation molecules (CD28/CD38/CD81). Fluorescence in situ hybridization analysis of fluorescence-activated cell sorter–sorted CTCs also unraveled different cytogenetic profiles vs paired BM clonal PCs. Moreover, CTCs were mostly quiescent and associated with higher clonogenic potential when cocultured with BM stromal cells. Most interestingly, CTCs showed a circadian distribution which fluctuates in a similar pattern to that of CD34⁺ cells, and opposite to stromal cell–derived factor 1 plasma levels and corresponding surface expression of CXC chemokine receptor 4 on clonal PCs, suggesting that in MM, CTCs may egress to PB to colonize/metastasize other sites in the BM during the patients’ resting period.

Genetic events mediating transformation from premalignant monoclonal gammopathies (MG) to multiple myeloma (MM) are unknown. To obtain a comprehensive genomic profile of MG from the early to late stages, we performed high-resolution analysis of purified plasma cells from 20 MGUS, 20 smoldering MM (SMM) and 34 MM by high-density 6.0 SNP array. A progressive increase in the incidence of copy number abnormalities (CNA) from MGUS to SMM and to MM (median 5, 7.5 and 12 per case, respectively) was observed (P=0.006). Gains on 1q, 3p, 6p, 9p, 11q, 19p, 19q and 21q along with 1p, 16q and 22q deletions were significantly less frequent in MGUS than in MM. Although 11q and 21q gains together with 16q and 22q deletions were apparently exclusive of MM status, we observed that these abnormalities were also present in minor subclones in MGUS. Overall, a total of 65 copy number-neutral LOH (CNN-LOH) were detected. Their frequency was higher in active MM than in the asymptomatic entities (P=0.047). A strong association between genetic lesions and fragile sites was also detected. In summary, our study shows an increasing genomic complexity from MGUS to MM and identifies new chromosomal regions involved in CNA and CNN-LOH.

Although information about the molecular pathogenesis of Waldenström's Macroglobulinemia (WM) has significantly advanced, the precise cell of origin and the mechanisms behind WM transformation from IgM MGUS remain undetermined. Here, we undertook an integrative phenotypic, molecular and genomic approach to study clonal B-cells from newly-diagnosed patients with IgM MGUS (n=22), smoldering (n=16), and symptomatic WM (n=11). Through principal-component-analysis of multidimensional flow cytometry data, we demonstrated highly overlapping phenotypic profiles for clonal B-cells from IgM MGUS, smoldering and symptomatic WM patients. Similarly, virtually no genes were significantly deregulated between FACS-sorted clonal B-cells from the three disease groups. Interestingly, the transcriptome of the Waldenström's B-cell clone was highly different than that of normal CD25^-CD22⁺ B-cells, whereas significantly less genes were differentially expressed and specific WM pathways normalized once the transcriptome of the Waldenström's B-cell clone was compared with its normal phenotypic (CD25⁺CD22^+low) B-cell counterpart. The frequency of specific copy number abnormalities [+4, del(6q23.3-6q25.3), +12, and +18q11-18q23] progressively increased from IgM MGUS and smoldering WM vs. symptomatic WM (18% vs. 20% and 73%, respectively; P=.008), suggesting a multistep transformation of clonal B-cells that albeit benign (i.e.: IgM MGUS and smoldering WM), already harbor the phenotypic and molecular signatures of the malignant Waldenström's clone.

Graft-versus-host disease (GvHD) remains the major obstacle to successful allogeneic hematopoietic stem cell transplantation, despite of the immunosuppressive regimens administered to control T cell alloreactivity. PI3K/AKT/mTOR pathway is crucial in T cell activation and function and, therefore, represents an attractive therapeutic target to prevent GvHD development. Recently, numerous PI3K inhibitors have been developed for cancer therapy. However, few studies have explored their immunosuppressive effect. The effects of a selective PI3K inhibitor (BKM120) and a dual PI3K/mTOR inhibitor (BEZ235) on human T cell proliferation, expression of activation-related molecules, and phosphorylation of PI3K/AKT/mTOR pathway proteins were analyzed. Besides, the ability of BEZ235 to prevent GvHD development in mice was evaluated. Simultaneous inhibition of PI3K and mTOR was efficient at lower concentrations than PI3K specific targeting. Importantly, BEZ235 prevented naïve T cell activation and induced tolerance of alloreactive T cells, while maintaining an adequate response against cytomegalovirus, more efficiently than BKM120. Finally, BEZ235 treatment significantly improved the survival and decreased the GvHD development in mice. These results support the use of PI3K inhibitors to control T cell responses and show the potential utility of the dual PI3K/mTOR inhibitor BEZ235 in GvHD prophylaxis.

Persistence of chemoresistant minimal residual disease (MRD) plasma cells (PCs) is associated with inferior survival in multiple myeloma (MM). Thus, characterization of the minor MRD subclone may represent a unique model to understand chemoresistance, but to our knowledge, the phenotypic and genetic features of the MRD subclone have never been investigated. Here, we compared the antigenic profile of MRD vs diagnostic clonal PCs in 40 elderly MM patients enrolled in the GEM2010MAS65 study and showed that the MRD subclone is enriched in cells overexpressing integrins (CD11a/CD11c/CD29/CD49d/CD49e), chemokine receptors (CXCR4), and adhesion molecules (CD44/CD54). Genetic profiling of MRD vs diagnostic PCs was performed in 12 patients; 3 of them showed identical copy number alterations (CNAs), in another 3 cases, MRD clonal PCs displayed all genetic alterations detected at diagnosis plus additional CNAs that emerged at the MRD stage, whereas in the remaining 6 patients, there were CNAs present at diagnosis that were undetectable in MRD clonal PCs, but also a selected number of genetic alterations that became apparent only at the MRD stage. The MRD subclone showed significant downregulation of genes related to protein processing in endoplasmic reticulum, as well as novel deregulated genes such as ALCAM that is prognostically relevant in MM and may identify chemoresistant PCs in vitro. Altogether, our results suggest that therapy-induced clonal selection could be already present at the MRD stage, where chemoresistant PCs show a singular phenotypic signature that may result from the persistence of clones with different genetic and gene expression profiles. This trial was registered at www.clinicaltrials.gov as #NCT01237249.