Patient-Derived Multiple Myeloma 3D Models for Personalized Medicine—Are We There Yet?

Despite the wide variety of existing therapies, multiple myeloma (MM) remains a disease with dismal prognosis. Choosing the right treatment for each patient remains one of the major challenges. A new approach being explored is the use of ex vivo models for personalized medicine. Two-dimensional culture or animal models often fail to predict clinical outcomes. Three-dimensional ex vivo models using patients’ bone marrow (BM) cells may better reproduce the complexity and heterogeneity of the BM microenvironment. Here, we review the strengths and limitations of currently existing patient-derived ex vivo three-dimensional MM models. We analyze their biochemical and biophysical properties, molecular and cellular characteristics, as well as their potential for drug testing and identification of disease biomarkers. Furthermore, we discuss the remaining challenges and give some insight on how to achieve a more biomimetic and accurate MM BM model. Overall, there is still a need for standardized culture methods and refined readout techniques. Including both myeloma and other cells of the BM microenvironment in a simple and reproducible three-dimensional scaffold is the key to faithfully mapping and examining the relationship between these players in MM. This will allow a patient-personalized profile, providing a powerful tool for clinical and research applications.


Introduction
Multiple myeloma (MM) is one of the most common hematological malignancies worldwide [1]. It is characterized by bone marrow (BM) infiltration by monoclonal plasma cells and overproduction of monoclonal immunoglobulins, leading to end-organ damage and significant patient morbidity [2,3]. Despite the introduction of several novel drugs, MM remains incurable, with more than 90% of relapsed patients developing drug resistance [4].
Myeloma cells grow within the BM, a hematopoietic organ in which various cell types establish a complex interplay. This involves adhesion molecules and soluble factors that play a critical role in disease progression, making the microenvironment more permissive for immune evasion, proliferation, survival and migration of myeloma cells [5][6][7][8]. Presently, only 5% of the new antitumor molecules end up obtaining clinical approval [9,10], reinforcing the importance of unveiling molecular mechanisms associated with the disease and identifying therapies with less toxicity and meaningful efficacy.
In order to develop the fittest and safest treatment for each individual, a deeper understanding of the underlying disparities between MM patients [11] and of the complex

The Importance of MM Bone Marrow Microenvironment
MM development is conditioned and supported by a wide range of cell types, namely BM stromal cells (BMSC), adipocytes, endothelial cells, osteoclasts, neurons and immune cells. Through soluble factors and/or direct contact, myeloma cells can modulate the phenotype of surrounding BM cells to evade immune surveillance, building a protective environment that confers growth advantages and resistance to therapeutic agents [5,19,20]. Hence, many approved therapeutic strategies can act directly on myeloma cells, or indirectly, by targeting other cell types of the microenvironment, including immune cells, osteoclasts or endothelial cells ( Figure 1).
Immune evasion is a hallmark of MM [21,22]. This includes a decrease in B cells with altered differentiation and antibody response [23], and an upregulation of inhibitory molecules either by myeloma cells, antigen-presenting cells (e.g., programmed deathligand 1 (PD-L1)), or lymphocytes (e.g., PD-1) [24,25]. Furthermore, the presence of tumorassociated macrophages with higher levels of CD163 or CC motif chemokine receptor 2 (CCR2) has also been associated with worse prognosis and drug resistance [26,27]. The impairment in immune homeostasis led to the development of novel therapies that significantly improved MM patients' survival in the last decade by inducing both tumoricidal and immune-boosting effects [28]. These agents include immunomodulators (IMiD) and, more recently, cereblon E3 ligase modulators (CelMoDs), proteasome inhibitors (PI), monoclonal antibodies (mAb), antibody drug conjugates and chimeric antigen receptor (CAR) T cells.
Although not yet approved in MM, other agents disrupting the crosstalk between myeloma cells and the BM microenvironment are also being tested, such as immune checkpoint inhibitors and CAR-NK cells [29][30][31].
The inhibition of osteoblast function and the activation of bone-resorbing osteoclasts are also found in the MM BM milieu [32,33]. This imbalance can be mediated directly by myeloma cells via the receptor activator of NF-B ligand (RANKL), hepatocyte growth factor (HGF) or interleukin (IL)-3, or indirectly via BMSC through the upregulation of RANKL and secretion of IL-6, IL-1β or tumor necrosis factor (TNF)-α [34]. While PIs such as bortezomib have shown to both inhibit osteoclast differentiation and stimulate osteoblast formation [35][36][37], the effect of IMiD is mostly known to affect osteoclastogenesis [38][39][40][41]. Nonetheless, existing reports suggest that combined therapies (e.g., lenalidomide, bortezomib and dexamethasone) may help to induce bone formation [42][43][44]. Moreover, both zoledronic acid, which is a bisphosphonate, and the mAb denosumab are approved for the prevention of osteolytic lesions [45].
The secretion of pro-angiogenic factors has also been associated with MM progression and poor prognosis. Indeed, several research groups have reported an increase in the secretion of the vascular endothelial growth factor (VEGF), hypoxia-inducible transcription factor (HIF)-1α, fibroblast growth factor (FGF), osteopontin (OPN) or angiopoietin (Ang)-1 by myeloma cells. These signals can happen between myeloma cells, BMSC and endothelial cells. Although pre-clinical tests using anti-angiogenic agents against these specific factors have shown promising results, their translation is still not a reality due to relapse and a limited increase in progression-free survival [46][47][48]. Nonetheless, other already approved MM therapies are also reported to target angiogenesis, such as IMiD, PIs or bisphosphonates [49][50][51][52][53].
New axes of development involve the CXC chemokine receptor (CXCR) 4 [54]. For instance, the therapeutic inhibition of CXCR4 (AMD3100) disrupts the interaction between myeloma cells and the BMSC, re-sensitizing myeloma cells to the PI bortezomib, leading to tumor reduction in vivo [55].
Adipocytes may represent another critical regulator of myelomagenesis by protecting myeloma cells from chemotherapy-induced apoptosis [56,57]. The production of fatty acid-binding protein (FABP) 4 by adipocytes upregulates the energy and metabolism of myeloma cells, making the BM microenvironment more pro-tumoral [58]. Preclinical studies strongly suggest that targeting adipocytes or adipose-derived soluble molecules is a worthwhile strategy to tackle MM and/or overcome chemotherapy resistance.
The development of new therapeutic strategies is challenged by serious side effects and/or unreproducible responses due to immune signature, microbiota, specific mutations and patient inter-variability [59][60][61][62]. The discrepancy between the results obtained in pre-clinical trials and clinical outcomes can be attributed to the limitations of the classic tissue culture models, which fail to successfully recreate the whole MM BM niche [8,63,64]. Indeed, MM should be considered a complex ensemble of cells, and over-simplified models are unlikely to recapitulate the complexity of the disease. in the next section, we will address the main limitations associated with current 2D cellular and animal models.

Historical Evolution, Limitations of 2D Cell Culture and Animal Models
Cell culture represents an indispensable tool for studying cancer cell biology and drug response and has come a long way since the first HeLa cell line was established in 1951 [65]. Scientists can now culture many cell types, including immortalized cancer cell lines, immune cells and even primary human cells [66,67].
Conventional 2D cell culture offers many benefits, including low cost, high-throughput screening capability, standardization and reproducibility [68,69]. However, this technique grows cells on plastic surfaces, resulting in remodeling of the cytoskeleton, altered gene expression and protein synthesis [70]. Studies in planar cell culture have shown how cells become progressively flatter, divide abnormally and lose their differentiated phenotype [71,72]. Two-dimensional MM models mostly use cell line monocultures, which have lost their BM dependence. These do not accurately represent tissue architecture or the interplay between cells, and lack the personal heterogeneity aspects of the disease [4,13]. Transcriptomic studies demonstrated that myeloma cell lines only reflect a limited segment of the whole in vivo profile of the tumor, making them unreliable preclinical models [73]. These also fail to reproduce nutrients and oxygen gradients, limiting the ability of precisely predicting drug sensitivity [74]. The aberrant metabolic profile of MM cells is also a key feature that should be considered. Metabolic pathways provide the energy required for cell proliferation and tumor growth, and are involved in drug resistance events [75,76]. Moreover, the extracellular environment has a major impact on cell structure, mechanic transduction and signaling [77]. In 2D culture, there is a lack of natural extracellular matrix (ECM) proteins, chemokines, growth factors and sites for cellular adhesion, which are crucial for cells to interact with adjacent cells and surroundings, preserving the specificity and homeostasis of the original tissue and its regular functioning [78,79].
Animal models, mostly murine models, provide advantages over 2D models, as the tumor is surrounded by a microenvironment, reproducing tumoral complexity. MM mouse models have shown to be relevant for the recreation of several 3D features of the BM. The 5TMM mice model has been a milestone for the comprehension of MM pathogenesis [80], and the implantation of human fetal bone chips in SCID mice has allowed an even better resemblance of BM physiology and the microenvironment importance in MM [81]. However, mouse models fail to accurately recreate human disease conditions. For example, immunodeficient mice impair the study of immune and cancer cells' interactions in the microenvironment. Moreover, murine and human tumor microenvironment and immune systems vary significantly. Discrepancies in both innate and adaptive immunity include the balance of leukocyte and immunoglobulin subsets, B-cell and T-cell signaling pathway components, cytokines and cytokine receptors, Th1/Th2 differentiation, the antigen-presenting function of endothelial cells, and chemokine and chemokine receptor expression [82]. Humanized mice could solve this problem, allowing reconstitution of the human immune landscape and cancer-stroma interactions [12]. Nonetheless, studies have highlighted troubles with graft-versus-host disease [83]. Interestingly, Meyer et al. developed patient-derived tumor xenograft models of human acute lymphoblastic leukemia and acute myeloid leukemia, and reported that the murine environment selected specific subclones, resulting in a number of different sub-models [84]. Moreover, this is an expensive time-and resource-consuming approach. Taking into account the timeframe for drug testing, results may be obtained after the patient has suffered mutations and/or developed metastasis, compromising the effectiveness of the selected treatment regimen [85]. The use of poorly validated animal models and the non-reflection of the genetic heterogeneity of human cancers may also constitute some of the reasons that lead to reduced scientific validity and reproducibility of the in vivo studies in biomedical research [69,86]. Finally, with the increasing emphasis on animal welfare, murine models are also known to involve ethical issues. Accordingly, 3D cell culture models may help to decrease the use of laboratory animals in drug testing, going in line with the 3R principles: Replacing animals with alternative methods, Reducing the use of animals in research, and Refining discomfort, whenever possible [87].
All these aspects reinforce the need for a personalized, biologically relevant, complete and accurate preclinical MM model. The idea of 3D culture was first conceived in 1912 by Alexis Carrel, when this scientist cultured an explant from a chick embryo and maintained it over a period of 3 months [88]. Later, aiming to create a cell culture environment that resembled the human body more closely, Hamburger and Salmon developed a model using a soft agar solution [89]. Over the last two decades, 3D culture methods have been widely used to study biology in multiple cancer types (e.g., lung, liver, brain), giving light to new discoveries in the areas of metastasis, hypoxia, angiogenesis and drug screening [90][91][92][93]. Three-dimensional cell culture models enable cells to maintain their natural morphology and tumor architecture, with a proliferating zone and a quiescent region, with limited oxygen, nutrient and growth factor distribution, which might influence drug response [69,70]. With respect to MM, it has been proven that, for example, cytokine production (e.g., IL-6, IL-11 or HGF) is higher in 3D versus 2D cultures [94]. De la Puente et al. showed that their 3D model promoted myeloma cell proliferation better than 2D systems [95]. Spelat and colleagues' 3D gel was able to maintain pluripotent stem cells in a G0 quiescent state without inducing proliferation or differentiation commitment, but maintaining the potential to do it [96].
Thus, 3D models are proposed as viable and biomimetic alternatives. They might allow for a deeper understanding of molecular disease events and of multidrug resistance mechanisms, improving the predictive value of drug efficacy and safety. The following section depicts 3D MM BM models that are currently available or under development.

The Reconstructed Endosteum-Bone Marrow Model
In 2008, a key progress in 3D MM modelling was achieved by Kirshner et al. with the reconstruction of the human MM BM microenvironment in the "rEnd-rBM" model (Table 1) [97]. Twenty-four-well plates were pre-treated with fibronectin and collagen type I, creating the reconstructed endosteum-marrow junction (rEnd) compartment, and this layer was covered with patient BM mononuclear cells (BMMC) suspended in a gel mixture of Matrigel and fibronectin, creating the recombinant BM (rBM) compartment [97]. Cells were cultured in growth medium supplemented with the patient's own plasma and kept viable for up to 21 days. Curiously, plasma from healthy donors did not support myeloma cell growth, which reinforces the importance of patient specificity. When compared with the patient BM, the rBM environment showed similarities with natural BM niche architecture. Moreover, this breakthrough enabled the proliferation of myeloma cells, including putative stem cells. This is a very relevant feature, as this cell subset seems to be one of the main subsets responsible for drug resistance and relapse [98,99]. Interestingly, an increase in the total number of clonotypic cells was observed, and cells retained the same chromosomal abnormalities present in vivo, as observed after 15 days of culture. Obtaining the entire stromal compartment is a marked improvement of this model over standard culture methods. It allows the testing of several drugs with different targets and the understanding of how these drugs act, not only on myeloma cells, but also on non-tumoral cells. Furthermore, the interaction with stroma may be enough by itself to influence molecular mechanisms related to drug sensitivity and resistance. Results showed the inability of the system to support populations of CD34+ hematopoietic stem cells, CD20+ B cells, and CD138+ plasma cells without the stromal compartment of the BM, highlighting its importance [100].
More recently, based on the "rEnd-rBM" model, Kirshner established zPredicta Inc., commercializing Reconstructed Bone (r-Bone™) technology. This customizable organspecific 3D platform includes both hematopoietic and stromal cells, as well as extracellular components, to sustain primary human BM cell survival, proliferation and variety [101]. BMMC from 21 newly diagnosed or relapsed patients with MM were set up in r-Bone and cultures were treated according to the clinical single-agent or combination regimen received by the patients. Subsequently, myeloma cell death was evaluated by flow cytometry. The system successfully predicted the patient's clinical outcome in 90% of the cases [101]. This technology has already been applied in various medical contexts, growing primary BM cells from healthy donors and from patients with MM and amyloidosis.
Huang et al. adapted the rBM concept with the goal of studying the activation of the signal transducer and activator of the transcription (STAT) 3 pathway in MM in 3D versus conventional 2D cultures (Table 1) [102]. The STAT3 pathway was activated when cells were cultured in 3D but remained inactive in conventional 2D cultures, showing how some MM mechanisms are dependent on the 3D structure. Moreover, inhibition of the STAT3 pathway, using the pharmacological selective inhibitor Stattic, significantly decreased the viability of myeloma cells [103]. Furthermore, susceptibility to the PI bortezomib was increased [102], providing insight on targeting STAT3 in MM treatment. Overall, the obtained results suggest that the 3D environment is crucial for recapitulation of ECM proteins and cytokines' interactions with myeloma cells. This is valuable because MM is usually characterized by deregulated levels of both pro-and anti-inflammatory cytokines [104].
Adapting Huang's work, Caillot et al. tested whether the inhibition or overproduction of reactive oxygen species (ROS) could sensitize myeloma cells to bortezomib [105], as these cells are known to produce high levels of ROS [106]. Upon modulation of the redox balance, apoptosis and autophagy was increased in myeloma cells [105], suggesting that this may be useful in relapsed/refractory patients, partially reversing tumor microenvironmentmediated drug resistance. Nonetheless, the inclusion of other components of the BM microenvironment would be of interest. Moreover, further characterization of the primary myeloma cells embedded is needed, as most assays were conducted using MM cell lines.

3D Matrigel Matrix-Embedded Models
Using Matrigel and fibronectin to form a 3D matrix, Perez et al. studied myeloidderived suppressor cells [107], whose frequency is commonly increased in myeloma niches [108]. BM samples from MM patients were cultured, and granulocytic cell subsets were identified by flow cytometry before and after exposure to the mAb daratumumab [107]. Cultures were maintained for 10 days, and no differences were seen in vitro in the percentage of granulocytic subsets between timepoints of drug exposure.
Cucè and colleagues used a Matrigel-spheroid model to evaluate apoptosis, cell cycle, and changes in cytokine production and release, both in MM cell lines and patient-derived primary myeloma cells (Table 1) [109]. Cultures were exposed to increasing concentrations of trabectedin both in 2D and 3D monoculture or co-culture with monocytes from healthy donors to understand the impact on the nucleotide excision repair pathway [109], previously described to be dysregulated in MM [110]. Trabectedin was able to alter the pro-inflammatory cytokine network, reducing monocyte chemoattractant protein (MCP)-1, VEGF and IL-10 [109]. The authors verified that the presence of monocytes was of utmost importance to recreate the protective effect of the microenvironment on myeloma cells, since 3D Matrigelspheroid co-culture promoted myeloma cell viability and proliferation in the presence of trabectedin. Spheroids are of particular interest to cancer researchers because they can include heterogeneous cell populations with both proliferating and quiescent cells, owing to limited oxygen and nutrient transport [9,70]. Being scaffold-free, these models take advantage of the natural ability of the cell types to self-aggregate and secrete their own ECM over time [69]. Nevertheless, spheroids have important limitations because they grow as independent cellular aggregates and have reduced interactions with the extracellular setting [111]. Still, if combined with more recent techniques, such as microfluidics or 3D bioprinting, scientists could be able to develop more physiologically relevant 3D models of the disease. These approaches are further described in Sections 4.3 and 4.7 Very recently, Khan and others put their efforts into developing a 3D platform with the best homology possible to the native hematopoietic tissues, thus allowing the engraftment and growth of primary cells from patients with blood cancers, including MM (Table 1) [112]. For that, differentiated human induced pluripotent stem cells (iPSC) were embedded in a hydrogel with an optimized composition-Matrigel and type I and type IV collagen-in order to support differentiation of the 3D BM perivascular niche [112]. Organoids were obtained and contained hematopoietic stem and progenitor cells, as well as stromal and myeloid cellular subtypes, with homology to human BM cellular subtypes, as confirmed by comparison of the molecular profiles identified by single-cell RNA sequencing. Moreover, cell-cell communication and a vascular network were observed. Interestingly, myeloma cells isolated from BM aspirates of MM patients engrafted into organoids presented a viability of more than 90% at day 10 of culture, whereas in 2D culture, they started to die after only 2 days. These results confirm that patient-derived MM cells require a supportive niche to proliferate and stay viable in ex vivo conditions. This system addresses a major problem of most current models: the vascularization. However, it does not include key populations, such as lymphoid or osteoid cells. Mesodermal, endothelial and hematopoietic lineages were obtained by differentiation of a human episomal cell line and not primary differentiated cells, which can lead to a suboptimal mimicking of the BM physiology.
Using Matrigel as the scaffold also, Braham et al. cultured primary myeloma cells with mesenchymal stem cells (MSC) and endothelial progenitor cells, to facilitate their survival and proliferation, over a 28-day period (Table 1) [113]. Their aim was to investigate the use of a novel class of αβ-engineered T cells to express a defined γδ T-cell receptor (TCR), called TEG, in targeting and eliminating primary myeloma cells. These TEGs migrated inside the Matrigel matrix, finding and killing their targets after 48 h [113]. The 3D model was far more effective when compared to a similar 2D approach. The same 3D model was later employed to test liposomal delivery of drugs (e.g., doxorubicin, bortezomib) [114]. Additionally, it was adapted by other groups to test the efficacy of γδ TCR anti-CD3 bispecific molecules in redirecting T lymphocytes against myeloma cells, while leaving healthy tissues intact [115]. The authors exploited this system to provide pre-clinical in vitro testing of different therapies on primary BM samples from seven relapsed/refractory MM patients [116]. Interestingly, the responses of each donor cell to the given therapies were analyzed according to two different readout parameters (percentage of dead myeloma cells and number of live myeloma cells), resulting in different outcomes. This reflects the importance of selecting the best predictor of clinical response. Nevertheless, the model showed a very poor agreement between in vitro and clinical treatment responses to IMiD. This is most likely due to the lack of events' recapitulation within the BM microenvironment, such as angiogenic and inflammatory processes, possibly because of incomplete representation of BM cellular subpopulations. Braham et al. also used a bioprinting method with bioactive and biocompatible calcium phosphate cement disks to create separate but interacting endosteal and perivascular subniches of the MM BM, representing an interesting option for studying myeloma-bone interactions. The 3D BM model with combined subniches significantly increased the proliferation of CD138+ myeloma cells [117]. However, even though this is a more complete model, it still does not fully represent the whole BM microenvironment as it lacks, for instance, immune system interactions.

Bioprinted BM MM Models
Cell bioprinting is an emerging approach for 3D cancer cell patterning that simplifies the control of spatial and temporal distribution of cells in a biocompatible material [118]. In 2019, Rodriguez et al. developed a bioprinted MM organoid model using a hydrogel scaffold that combined fibronectin, denatured collagen (Gelin-S), a photo-initiator and a crosslinker, with primary BM aspirates from MM patients (Table 1) [119]. Threedimensional bioprinters were programmed to extrude 200,000 cells per well in a 96-well plate. Organoids were maintained with a viability above 70% long enough to assess chemosensitivity after one week [119]. This was the first MM organoid model produced using high-throughput bioprinting.
More recently, Wu et al. co-cultured myeloma cells with HS-5 stromal cells in a coaxial extrusion bioprinted construct that consisted of a stiff mineral outer alginate layer and a soft hydrogel core (Table 1) [120]. The authors showed that patient-derived myeloma cells were maintained with good cell viability for up to 7 days, compared to 5 days in the 2D environment [120]. Interestingly, 2D-cultured cells showed higher toxicity to bortezomib than 3D cultures, pointing to the potential improvement in physiological relevance of 3D platforms. Although promising, this model encompasses a low cellular diversity, hindering real mimicking of BM cellular interactions. Common limitations of this method include slow printing speed and the challenge of developing nontoxic bioinks. Although this is still an underexplored approach in the field of MM, its microscale resolution, high precision in forming 3D constructs, and ability to use multiple materials point to the potential of the technique, especially if using the right biomimetic compatible matrix and/or combined with techniques such as microfluidics, addressed further ahead [121,122].

The 3D Tissue-Engineered BM Model
Contributing to another key advancement, de la Puente et al. developed a 3D tissueengineered BM (3DTEBM ® ) model that consisted of a mixture of myeloma, endothelial and stromal cells from the BM of MM patients (Table 1) [95]. Fibrinogen naturally found in the human plasma and calcium chloride were added for promoting clotting and the crosslinking reaction. Tranexamic acid was added to provide stability to the structure for up to 7 days [95]. Proliferation rates of myeloma cells increased 250% when co-cultured with stromal and endothelial cells, compared to monoculture. Proliferation rates of myeloma cells in the 3DTEBM ® were also compared with conventional 2D multi-culture and commercially available 3D models, but these failed to sustain myeloma cell growth, showing the importance of recreating a 3D biomimetic environment. Notably, 3DTEBM ® allowed the recreation of drug and oxygen gradients, as suggested by the upregulation of HIF-1α and proto-oncogene serine/threonine-protein kinase (PIM) hypoxia markers, downregulation of CD138 and overexpression of CXCR4 in deeper areas of the scaffold. In addition, the drug uptake measured in cells demonstrated an inverse correlation with depth. Thus, this model seems to be able to recreate the conditions observed in vivo, with two distinct regions, where cells in proximity to vasculature are more proliferative and sensitive to therapy than those in the endosteal niche. This is an important achievement, as it can influence the prediction of resistance to delivered drugs [123][124][125].
Due to its impressive results in resembling the BM niche, the team further explored this model. For instance, they tested the potential of targeting CD47 for immunotherapy [126]. Results showed that when CD47-expressing MM cells were co-cultured with macrophages and treated with an anti-CD47 mAb, a significant killing of MM cells was seen in the 3DTEBM ® , but not in classic 2D cultures. Extensive motility of macrophages during the phagocytosis process was observed thanks to the hydrogel-like structure of the 3DTEBM ® , which faithfully simulates the in vivo MM BM niche, in opposition to the adherent nature of 2D cultures. Most recently, the 3DTEBM ® model was used to assess different therapies on 19 primary patient samples and found a 89% predictable rate when comparing in vitro results to clinical outcomes [127]. These remarkable drug efficacy findings obtained using the 3DTEBM ® technology suggest that this is a highly promising approach for studying primary BM malignancies and those of hematologic origin, such as MM and acute myeloid leukemia.

Hydrogel-Based MM BM Models
Jakubikova et al. used PuraMatrix™ (Becton Dickinson, Franklin Lakes, NJ, USA), a synthetic self-assembling peptide hydrogel (Table 1) [128]. BMMC were obtained from BM aspirates and the adherent fraction containing MSC was expanded for several days before seeding to promote a higher degree of similarity with the in vivo tissue [128]. Cultures were maintained for up to 21 days. A higher expression of both cytokines and ECM molecules known to promote in vivo stem cell expansion and MM survival was observed in the 3D MSC model. These included IL-6, IL-8, MCP-1, VEGF, and collagen I, collagen IV, fibronectin and laminin, respectively. Protein expression analysis comparing 3D and 2D models revealed activation of osteogenesis and osteoclastogenic differentiation, as shown by the upregulation of matrix metallopeptidase 13, osteopontin, and matrix Gla protein and by the downregulation of calmodulin 1 and CXCR4. They also tested a wide range of novel and conventional anti-MM drugs. Particularly, using BM samples from MM patients, sensitivity results between 3D and 2D cultures showed that the drug response in the 3D co-culture model paralleled clinical resistance. Although promising, a higher number of patients are needed to validate these results. Authors report that prospective trials are being conducted to assess the model value in testing personalized immune therapies in MM [128].
Regarding a different approach, Waldschmidt et al. established a conical agarose 3D platform for the propagation of primary BMMC ex vivo (Table 1) [129]. First, myeloma cells growth in the system were compared to a 2D monolayer culture. In 3D, myeloma cells showed a slower initial proliferation, but maintained a more stable growth after 12 days of culture [129]. When comparing 3D monoculture and co-culture conditions, primary myeloma cells did not expand in 3D monoculture, requiring co-culture support by a human stromal cell line (e.g., HS-5, MSP-1). Similar to other groups, detailed cytokine characterization within the developed 3D model revealed a high expression of pro-angiogenic and pro-inflammatory cytokines, such as IL-6, IL-1β, IL-8 and TNF-α [95,128]. This secretion seemed to be stimulated by the interaction of myeloma cells with a stromal co-culture partner, providing the optimal conditions for tumoral cell proliferation and allowing the maintenance of cells in culture for 21 days. Obtained results are in line with studies where newly diagnosed MM patients' cytokine profiling revealed high levels of IL-6, IL-8 and TNF-α when compared to healthy controls [130]. Therefore, such cytokines seem to be appealing targets for treating MM, and tests can be conducted using the aforementioned platforms. Drug sensitivity tests with bortezomib and auranofin also showed that less toxicity was induced under 3D versus 2D condition and in co-versus monoculture [129,131].
Together, these results point to the important role of stromal co-culture for the in vitro modeling of drug resistance events and to correctly determine the efficacy of therapeutic compounds. From a technical point of view, this conical agarose 3D model shows some advantages, as it allows for clear-cut observation and monitoring by confocal microscopy, flow cytometry or Western blot [132].
Hyaluronic acid (hyaluronan or HA) hydrogels also represent a reliable alternative. HA is an increasingly popular biologically derived matrix, due to its abundant presence in the natural ECM. HA is easier to chemically manipulate compared to other matrix molecules and it allows for the activation of a variety of pathways, including those affecting cell adhesion and motility, inflammation and drug sensitivity [133]. Although being tested in several types of cancer, namely those of the brain, breast and liver [134][135][136][137], few data are available in MM. In 2014, Narayanan et al. developed a 3D HA-based model and showed how matrix composition and stiffness can impact results obtained (Table 1) [138]. In this study, HA-based hydrogels were developed for encapsulating BMSC and myeloma cells. The percentage of survival of myeloma cells grown on mediumstiffness (≥90%) hydrogels was higher than those grown on low-or high-stiffness hydrogels (70-80%) [138]. Proliferation was also higher in medium-stiffness HA hydrogels when compared to 2D or Matrigel cultures. As is well-accepted, during cancer, cells are subjected to a variety of mechanical changes in their microenvironment, resultant in part from ECM remodeling, with events such as increased crosslinking activity and deficient matrix degradation [139]. These events lead to increased matrix stiffness, which in turn activates transcription factors that promote tumor cell proliferation and growth [140,141]. A study conducted with monoclonal gammopathy of undetermined significance and MM patients showed that altered ECM protein expression (e.g., osteopontin, periostin, entactin and fibulin) is associated with decreased overall survival [142]. Therefore, using 3D models for studying proteins that are involved in matrix crosslinking and stiffness management may help identify new targets, if found to be dysregulated in MM.  (Table 1) [143]. This rotary cell culture system allowed the creation of a dynamic model in which cells are in a "free fall" form, favoring mass transfer of nutrients and low shear stress conditions [143,144]. On a first approach, this was designed for culture tissue explants from MM patients, with well-preserved tissue architecture and cell viability. Through histology, it was possible to observe that myeloma cells were kept viable within their microenvironment, as well as the vessels that appeared disrupted in static culture. Assessment of β2-microglobulin levels in supernatants from bortezomib-treated samples and in patients' sera after bortezomib-based therapy showed an overall agreement in the drug response ex vivo and in vivo. However, there were issues with tissue harvesting and reproducibility: (1) the obtained biological samples are very sensitive to manipulation, requiring careful handling to preserve integrity; (2) the number and size of the tissue fragments need to be thoroughly determined to counterbalance the hydrodynamic forces generated in the culture chamber and maintain the free fall condition; and (3) for different conditions tested in parallel, each culture vessel needs to have the exact amount of samples and comparable weight and volume [144].

Dynamic MM BM Models
In a follow-up study, the same group focused on isolated primary myeloma cells (Table 1) [145]. These were co-cultured with allogeneic BMSC and human umbilical endothelial vein cells (HUVEC), supporting the survival of myeloma cells for up to 7 days, in the bioreactor. Of note, genomic analysis was performed in a MM patient and compared to the bioreactor culture, showing that it paralleled the expansion of the clone that dominated in vivo [145]. Interestingly, signaling pathways associated with tumor survival, proliferation and drug resistance were upregulated in the 3D bioreactor model when com-pared with conventional 2D models. This indicates that this 3D model can be useful for studying drug resistance and its association with specific genomic alterations. Overall, this construct allowed for the reproduction of several MM characteristics and tumor-stroma interactions. The use of allogeneic BMSC in co-culture with patient-derived myeloma cells can be considered a strength of this study when compared to the models developed using stromal cell lines (e.g., HS-5 [105,120,131]). However, one could wonder if using both autologous BMSC and HUVEC would allow for a better recapitulation of the patient's BM niche specificity.

3D Models of Myeloma-Bone Interactions
Despite carrying advantages, the rotary system model requires sophisticated devices [146]. A more affordable approach includes microfluidic devices. Microfluidics seems to be a promising technology for 3D cell culture as it brings the possibility of creating biomimetic structures in which biological processes can be recapitulated [147], including circulation and cellular compartmentalization, in a controlled physiological environment [148,149]. in the last number of years, this micromanipulation technique has been employed in efforts to develop a system that allows for the easy placing of cells and matrices into the culture chambers of the device and supports replication and characterization of tissues and tumor microenvironments at a low cost [150].
Zhang et al. focused on the continuous nourishment of the cells, attempting to replicate the body's circulation in an 8-chamber microfluidic device (Table 1) [151]. Primary BMMC were introduced in each chamber after generation of an ossified tissue scaffold from a human osteoblastic cell line (line hFOB 1.19), resembling the endosteal surface, to mimic interactions between the two cellular populations [151]. Real-time monitoring confirmed that myeloma cells expanded after 21 days in culture and were drawn towards the osteoblast layer. The group also compared cell viability and proliferation of cells in their model to a 3D static culture, showing that perfusion in the endosteal niche is a key factor of the microenvironment and improves maintenance of long-term culture of myeloma cells [152,153]. However, the authors showed concerns regarding the microfluidic device material polydimethylsiloxane, as this has been reported to leach hydrophobic components from cell culture media, including drugs, antibodies and growth factors, potentially affecting the system's reproducibility, which is essential for validation steps and clinical testing [153]. Moreover, this model is based on a simple co-culture of myeloma and bone-derived cells. Because a monolayer of osteoblasts is cultured on the surface of the microfluidic device, and then myeloma cells are pumped into the chambers, this methodology raises the question of whether a 3D scaffold should be considered here. The addition of multiple patient-derived BM components would be of interest to increase the quality and complexity of the model, paving the way for relevant studies of cell adhesion-mediated drug resistance, for example.
In this context, Reagan et al. attempted to recreate the BM microenvironment by stimulating MSC to undergo osteogenic differentiation inside porous silk scaffolds prior to the inclusion of myeloma cells (Table 1) [154]. Cell viability after bortezomib treatment was determined, showing drug resistance in their 3D culture, when compared to hydrogel and 2D cultures [154]. These findings suggest that the inclusion of bone cells and recreation of a bone-like component are important factors for cell adhesion-mediated drug resistance. Additionally, the authors were able to use this model to identify a new microRNA signature associated with dysfunctional osteogenesis of MSC in contact with MM cells [154], showing that 3D models can be used to identify novel biomarkers and potential new therapeutic targets [155]. These 3D silk scaffolds have also been employed to develop the first 3D tissueengineered BM adipose tissue (BMAT) model. Briefly, MSC were seeded in silk scaffolds and cultured with adipogenic media, driving differentiation of adipocytes. With this approach, Fairfield et al. were able to shed light on complex interactions between BMAT and tumor cells, even if no patient-derived cells were used [156]. Compared with 2D cultures, pathways of DNA replication, metabolic and proliferation events were upregulated in BM adipocytes in this model. Of note, a study combining in vitro and in vivo methods showed that BM of MM patients contained, in fact, an increased number of both preadipocytes and mature adipocytes, pointing to its potentiality as a target and to the value of using a 3D model that includes cells of the adipose tissue in culture [157].

Computational and Mathematical MM BM Models
Over the years, mathematical models have been developed to simulate tumor growth and to study tumor dynamics in response to therapy and drug resistance [158,159]. Although these models are very powerful tools for the analysis of complex interactions within the tumor microenvironment, they are often derived from experimental conditions that are very controlled, hindering their translation to the clinical practice.
In 2014, Khin et al. combined in vitro microfluidics and in silico approaches to assess drug response in MM. Myeloma and stromal cells from 10 MM patients were cultured in 3D culture slides with collagen type I, and the effect of bortezomib and melphalan was tested [160]. An algorithm of digital image analysis detected the number of live myeloma cells by the motion of the cell membrane for different drug concentrations at different timepoints. The viability measurements obtained in vitro were fitted to the mathematical model of chemosensitivity, representing one or two subpopulations, with specific size, doubling time and level of sensitivity to the drug. Data derived from these experiments were then used to parameterize mathematical models for clinical outcome simulation [160].
This model was improved through the development of an ex vivo platform capable of making a three-month prediction regarding clinical response [63]. The Ex vivo Mathematical Myeloma Advisor (EMMA) combined a digital image analysis algorithm, mathematical models and pharmacokinetic data to predict the effect of 31 chemotherapeutic agents and IMiD. Primary myeloma cells were seeded with human BMSC in medium supplemented with patient-derived plasma. A digital microscope snapped live images, which were gathered every 30 min over 4 days. Nonetheless, the described system was unable to quantify cellularity of adherent cancer cells [161]. For the characterization of tumor heterogeneity with different degrees of chemosensitivity, the group considered and tested four hypotheses: (1) the existence of only one clonal population in which all tumoral cells presented the same degree of sensitivity to a particular drug; (2) the existence of only one population, but with a chemosensitivity following a normal distribution; (3) the existence of two clonal subpopulations; or (4) the existence of two normal distributions [63]. This may be advantageous because patients presenting heterogeneous tumors can exhibit clones that are more resistant to therapy. Even though more cells of the BM microenvironment are not present in the model (e.g., immune cells), the results obtained with this approach were very promising.
Nonetheless, the use of EMMA for drug combinations revealed some flaws: predictions of response to therapy were generated by simulating each of the drugs independently and combining all responses, assuming additivity [63]. In 2020, the same group extended EMMA into a synergy-augmented model (SAM), which captured interactions between drugs, transforming the fixed combined effects of the in vitro drug concentrations into clinically relevant, time-varying ones [64]. The modified damages caused by a drug when in the presence of another drug was considered and deemed to accumulate over time. Cell death happened when either drug exceeded the corresponding tumor-specific threshold. To choose the model that best represents data of tumor heterogeneity, the Akaike information criterion is applied, balancing the goodness of fit with the complexity of the model [162]. Then, tumor drug-specific parameters are combined with pharmacokinetic data, simulating patient-specific clinical responses to each drug. To do so, tumor growth must be estimated. "Median lethal dose" (LD50) and "area under the curve" are metrics that cannot handle pharmacokinetic/pharmacodynamic complex relations [63,64]. For that reason, in this context, the best response is predicted based on clinical response parameters (from EMMA and SAM models), combined with pharmacokinetic data. Still, we cannot fail to note that this model was developed without the use of a validation cohort [63]. Although the methodology used does not require a different group, it would be a great asset to validate the results in a completely independent cohort.
EMMA was built considering a cohort of 52 MM patients. This model correctly classified 50 out of 52 patients (96%) according to response/no response to treatment, and demonstrated that 60% of the patients received at least one drug with no clinical efficacy [63]. For the improved SAM model, patient-specific parameters were estimated by fitting ex vivo drug sensitivity data of 203 patients with pharmacokinetic data [64]. The twodrug combination effect was evaluated for 46 different combinations and their synergistic or antagonistic relation was assessed, and only 4 combinations were considered to be clinically synergistic [64]. Remarkably, this result shows that the synergism observed in the bench does not always translate into the bedside, referring to the benefit of combinations over single-agent regimens mainly due to independent drug action.
Although there is no simple formula for building the perfect model, there are some characteristics that we consider essential in obtaining a system that is both representative and easily applied in clinical practice. In general terms, the mathematical/computational approach needs to extrapolate in vitro predictions into a clinical outcome. To this end, the effect that different drugs have at the cellular level, not only singularly but also when combined, needs to be accounted for. Drug concentration, drug-induced damage and tumor growth equations are required. When in the presence of a combined therapy regimen, defining a dose-response function becomes more complex, since the synergetic effect must be considered, requiring, in addition to the effects of single agents, the inclusion of functions that define the effect due to the interaction of the drugs. Importantly, most patients are treated with triple therapies [163]. Therefore, testing the combined effect of more than two agents should be the next step to obtain a more accurate model. Nonetheless, this will significantly increase the complexity of the mathematical model. Additionally, parameters such as drug-induced repair rate and the drug-induced cell damage effect must also be included, as well as microenvironment-specific parameters regarding the immune compartment, matrix stiffness, and assessment of changes/effects in cells other than myeloma cells. Finally, the inclusion of -omics parameters (e.g., proteomics, genomics, transcriptomics), also possible to assess in vitro [164], must be evaluated, as they allow for even greater customization of the model to each patient's characteristics. Table 1. Available 3D patient-derived MM models. n of BM samples, number of samples obtained from patient biopsy and used for cell culture in the 3D model; n of patients, number of patients who entered a study of correlation of the obtained laboratory results in 3D culture versus clinical outcome. Abbreviations: BM, bone marrow; BMMC, bone marrow mononuclear cells; BMSC, bone marrow stromal cells; CaCl 2 , calcium chloride; CFSE, carboxy fluorescein succinimidyl ester; CFU, colony-forming unit; ECM, extracellular matrix; FISH, fluorescence in situ hybridization; GelMA, gelatin methacryloyl; HARV, high-aspect-ratio vessel; HUVEC, human umbilical vein endothelial cells; IHC, immunohistochemistry; iPSC, induced pluripotent stem cells; me-HA, methacrylated hyaluronic acid; MM, multiple myeloma; MSC, mesenchymal stem cells; nHA, nanohydroxyapatite; PBMC, peripheral blood mononuclear cells; PDMS, polydimethylsiloxane; PI, proteasome inhibitors; rBM, reconstructed bone marrow; RCCS, rotary cell culture system; rEnd, reconstructed endosteum; scRNA seq, single-cell RNA sequencing; WB, Western blot. * Approved patent.

An Eye into the Future-Challenges and Solutions
While 3D cell culture allows a more accurate representation of the natural environment of cells, these models also present new difficulties. One of the most important challenges to overcome concerns the ability to obtain and grow primary cells, maintaining their viability, integrity and phenotype. Primary cells usually take time to adapt and easily become senescent [165]. Therefore, time is a key factor, as tests must be conducted before the growth potential of primary cells is exhausted. Preferably, the culture should be carried out using freshly harvested biopsy samples. However, they are not always available, even for research groups that have a close link to medical institutes, which makes it difficult to develop, optimize and validate new assay platforms.
The obtained sample usually needs processing, either for isolation of specific types of cells or for removal of others (e.g., erythrocytes), introducing additional variations in tissue architecture, microenvironment and sample composition. Hence, we believe that it is of the utmost importance to take benefit of the whole sample, or as much as possible, to avoid intense manipulation and ensure the presence of all the different cell types, including stem, immune and bone cells. This will enable the study of their molecular mechanisms and drug susceptibility.
Co-cultures increase the complexity of the model and, therefore, require optimization of both cell ratios and cell media composition (including medium supplements specific to the needs of the cells in culture, such as growth factors or patient-derived plasma) to allow the proper growth of all cell types in terms of biological and disease relevance. Renewal of the medium is also an important aspect to consider: On one hand, the procedure must be done carefully to avoid disruption of the scaffold. On the other hand, the amount and frequency of medium renewal need to be carefully considered and adjusted according to each case, since the process might remove factors secreted by the cells that are most likely essential for their maintenance and can interfere with physiological behavior. Interestingly, Kirshner's model can maintain cell viability for 30 days without changing the medium [100]. The authors advise to only change half of the total medium volume at a time if medium renewal is required [100]. Due to the complexity of the immune system, which includes several subpopulations and activation states, co-culturing immune cells in a 3D MM model is a major challenge. Given the increasing interest in cancer immunotherapy, efforts are being made to create models that allow tumor-immune cell interaction [166,167]. It is of utmost importance to make sure that immune cell effector responses are not affected in 3D culture. Immune cell infiltration (e.g., lymphocytes, macrophages and myeloid cells), activation and cytokine secretion must be evaluated, making sure that cells are able to display their antitumor potential, which is mediated by immunotherapeutic agents [168][169][170]. Interestingly, the work of van Diest et al., using γδ TCR anti-CD3 bispecific molecules in a 3D BM niche model with stromal, endothelial and myeloma cells in Matrigel, revealed the successful infiltration of T cells and selective elimination of myeloma cells in culture [115]. The presence of functional tumor-stroma interactions in a 3D MM model may provide valuable information on therapeutic efficacy and adverse effects of immunotherapies [126,145].
The ideal model should provide the ability to identify key factors regulating tumor development, such as cell-cell interactions, cell-matrix interactions and cell receptors [9]. For that, it must comprise the following aspects: good recapitulation of biochemical cues, adequate tissue stiffness, permeability for penetration of drugs, nutrients and growth factors, recreation of oxygen gradients, sites for cell adhesion, compartmentalization and vascularization, and must allow proliferation, differentiation and migration. Previous experiments have shown how culturing the same cells in different biomaterials leads to changes in the obtained results [138,154]. Hence, selecting the scaffold to be used is of key importance in the experimental process, and should not be disregarded. Although several new options have emerged in the past 30 years, scaffolds of natural origin, namely Matrigel [171], collagen [172] and HA [173], are still the gold standard in cancer research. These are naturally recognized and remodeled by the cells, and possess cytokines and growth factors, which allow for the improvement of viability, while cells carry out their functions [174,175]. However, most of the ECM gels are obtained from animals or cultured cells, leading to batch-to-batch variation that reduces reliability and reproducibility [176]. To avoid such inconsistencies, some alternatives are recommended, such as growth-factorreduced Matrigel instead of the standard Matrigel, for instance. In addition, these gels often require fast handling and control of temperature or crosslinking activity in order to avoid premature polymerization of the matrix, a challenge in itself [177,178]. Recently, xeno-free hydrogel options are emerging that do not require crosslinking agents and can be manipulated at room temperature [179,180].
Improvements in imaging techniques are also fundamental in making the most out of the third dimension. Confocal microscopy enables visualization of cells in situ but, for now, it is fairly difficult to visualize 3D samples on microscopes due to their thickness and poor light penetration [181]. More transparent matrices are preferable for more straightforward monitoring using non-destructive methods. For example, Waldschmidt et al. state that their agarose hydrogel allows for easy examinations and does not absorb any of the fluorescent dyes incubated to stain cells for fluorescent microscopy, providing a clear background in imaging [129,132]. On the other hand, most 3D models can be fixed and embedded in paraffin, allowing slicing into several sections to be examined. However, even though immunohistochemistry provides a wide range of markers for visualizing cell behavior, it preserves limited spatial characteristics, reducing the 3D model down to 2D [182]. Fixation of the samples is sometimes demanding because some matrices tend to depolymerize with fixative agents and the addition of reagents that minimize this action can increase the imaging background [183]. The right balance between these agents must be found, to avoid compromising the value of the model. Attention should also be paid to the processing of the sample for histological cuts, as certain scaffolds are too sensitive and may tear apart. Notably, the 3D culture of Huang et al. was formed in a well with HistoGel TM that could be easily fixed, allowing for histologic processing and immunocytochemical studies [102,184].
Analyzing 3D samples by flow cytometry can be challenging, as samples must be in cell suspension, requiring dissociation by single or combined enzymes (e.g., trypsin, hyaluronidase, collagenase, dispase, cell recovery solution), which can affect cell integrity and viability. In addition, spatial distribution is lost. Spelat et al. introduced a nice advantage regarding this matter, as their worm-like micelle gel can be converted into freeflowing spheres simply by incubating the culture in a 4 • C solution, therefore avoiding enzymatic degradation, as required by many commercially available protein-based matrix gels [96]. Thanks to technological advances, new flow cytometers are being developed, aiming to allow the analysis and sorting of intact structures up to 1500 µm, instead of cell suspension [185,186].
As for drug screening, the major challenge regards the high tumor heterogeneity, not only at the inter-patient heterogeneity level, as analyses reveal high genomic heterogeneity across patients [187,188]; but also at the temporal and spatial intra-patient heterogeneity levels [74,189]. Anti-tumoral drugs might work out differently in cells at different stages, impacting the guidelines for treatment, and the whole setting of a patient's disease may not be represented simply by the portion of biologic material collected in the biopsy [190]. Time is also an obstacle because the model needs to be grown and validated before screening. A defined repertoire of drugs to be screened is also required for setting a realistic and convenient timeframe. Interestingly, Braham et al. showed that drug resistance in their model was only seen when administrating treatments 5 to 7 days after assembling 3D cultures. This time was necessary for the supporting cells in Matrigel to spread and form networks, pointing to the need to allow the establishment of myeloma-BM interactions before testing therapeutics [114]. The lack of drug-metabolizing enzymes in an ex vivo context can also impact the evaluation of anti-tumor drug effects in the 3D model. Considering the importance of metabolizing compounds over different tissues of the human body, including the liver and kidneys [191,192], studies using multi-organ systems have emerged. These allow for the establishment of organoids of different origin, linked and communicating on a chip, for integrated responses and more predictive assessments [193,194].
For now, there is still a lack of standardized procedures for creating reliable platforms to be used in personalized medicine. Most of the developed models are complex and recreating them is a major task for other groups, hampering the transition to the drug discovery industry. The ideal 3D MM model should be easy to perform, reproducible, timesaving, cost-effective, and should allow for high-throughput screening while still mimicking the tumor cells' behavior and their crosstalk with the BM microenvironment. As a result, a greater knowledge of the mechanisms and pathways involved in MM is expected and, consequently, developing drugs capable of engaging with them should become easier [9]. From our point of view, there are some paramount concerns to take into account when trying to assess cells' behavior and the therapeutic response: (1) Do we have all the necessary cells in culture, both in type and number? (2) Are all the BM compartments replicated? (3) What are the fundamental tumor niche components affecting disease progression? (4) What are the possible drug targets that must be present in culture? (5) What is the endpoint of interest and how do we consider it?
Including emerging technologies in pre-clinical studies such as microfluidics and computational modeling will most likely push 3D models forward in clinical translation ( Figure 2 For now, there is still a lack of standardized procedures for creating reliable platforms to be used in personalized medicine. Most of the developed models are complex and recreating them is a major task for other groups, hampering the transition to the drug discovery industry. The ideal 3D MM model should be easy to perform, reproducible, timesaving, cost-effective, and should allow for high-throughput screening while still mimicking the tumor cells' behavior and their crosstalk with the BM microenvironment. As a result, a greater knowledge of the mechanisms and pathways involved in MM is expected and, consequently, developing drugs capable of engaging with them should become easier [9]. From our point of view, there are some paramount concerns to take into account when trying to assess cells' behavior and the therapeutic response: (1) Do we have all the necessary cells in culture, both in type and number? (2) Are all the BM compartments replicated? (3) What are the fundamental tumor niche components affecting disease progression? (4) What are the possible drug targets that must be present in culture? (5) What is the endpoint of interest and how do we consider it?
Including emerging technologies in pre-clinical studies such as microfluidics and computational modeling will most likely push 3D models forward in clinical translation ( Figure 2).

Figure 2.
Key points for obtaining the ideal 3D MM model. Even though efforts to create an authentic and reliable MM model are increasing, there are still issues that need to be addressed and improved. Here, we sum up the main factors and variables that are essential and must be considered before, during and after developing a 3D MM model. They include sample processing, culture establishment, biological assays, validation and prediction of clinical outcomes. Created with BioRender.com. Abbreviations: BM, bone marrow; ECM, extracellular matrix; ELISA, enzyme-linked immunosorbent assay; IHC, immunohistochemistry; RT-PCR, reverse transcriptase polymerase chain reaction; SEM, scanning electron microscope; TME, tumor microenvironment; WB, Western blot.

Conclusions
The scientific community is becoming more aware of the importance of having a reliable and representative model of the in vivo scenario of MM. The success of the 3D model's establishment and therapy response prediction depends heavily on the ability of the model to recapitulate tumors' biology and physiology as well as their microenvironment. While there are still many challenges to overcome, including finding the right balance between the necessity to achieve the perfect model and the ease of its usage, it seems clear that the ultimate model should be based on patient samples, as every tumor has its own individuality. This would provide a powerful tool for patient-tailored profiling in pathophysiology studies and drug response, potentiating its translation from bench to bedside.