Lee LX, Li SC. Hunting down the dominating subclone of cancer stem cells as a potential new therapeutic target in multiple myeloma: An artificial intelligence perspective. World J Stem Cells 2020; 12(8): 706-720
Corresponding Author of This Article
Shengwen Calvin Li, PhD, Professor, Neuro-oncology and Stem Cell Research Laboratory, CHOC Children's Research Institute, Children's Hospital of Orange County, 1201 W La Veta Ave., Orange, CA 92868, United States. firstname.lastname@example.org
Checklist of Responsibilities for the Scientific Editor of This Article
This article is an open-access article which was selected by an in-house editor and fully peer-reviewed by external reviewers. It is distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Author contributions: Lee LX and Li SC conceived and wrote the manuscript; both revised the manuscript; and all authors approved the final version submitted.
Supported bythe CHOC-UCI Joint Research Award (in part).
Conflict-of-interest statement: The authors declare no conflict of interest.
Open-Access: This article is an open-access article that was selected by an in-house editor and fully peer-reviewed by external reviewers. It is distributed in accordance with the Creative Commons Attribution NonCommercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Corresponding author: Shengwen Calvin Li, PhD, Professor, Neuro-oncology and Stem Cell Research Laboratory, CHOC Children's Research Institute, Children's Hospital of Orange County, 1201 W La Veta Ave., Orange, CA 92868, United States. email@example.com
Received: May 18, 2020 Peer-review started: May 18, 2020 First decision: June 5, 2020 Revised: July 8, 2020 Accepted: August 14, 2020 Article in press: August 14, 2020 Published online: August 26, 2020
The development of single-cell subclones, which can rapidly switch from dormant to dominant subclones, occur in the natural pathophysiology of multiple myeloma (MM) but is often "pressed" by the standard treatment of MM. These emerging subclones present a challenge, providing reservoirs for chemoresistant mutations. Technological advancement is required to track MM subclonal changes, as understanding MM's mechanism of evolution at the cellular level can prompt the development of new targeted ways of treating this disease. Current methods to study the evolution of subclones in MM rely on technologies capable of phenotypically and genotypically characterizing plasma cells, which include immunohistochemistry, flow cytometry, or cytogenetics. Still, all of these technologies may be limited by the sensitivity for picking up rare events. In contrast, more incisive methods such as RNA sequencing, comparative genomic hybridization, or whole-genome sequencing are not yet commonly used in clinical practice. Here we introduce the epidemiological diagnosis and prognosis of MM and review current methods for evaluating MM subclone evolution, such as minimal residual disease/multiparametric flow cytometry/next-generation sequencing, and their respective advantages and disadvantages. In addition, we propose our new single-cell method of evaluation to understand MM's mechanism of evolution at the molecular and cellular level and to prompt the development of new targeted ways of treating this disease, which has a broad prospect.
Core tip: Current methods for determining prognosis in multiple myeloma are limited. The prototype device called Multi-Phase Laser-cavitation Single Cell Analyzer can perform reverse transcriptase polymerase chain reaction (RT-PCR) on single cells in a one-step microfluidics chip platform. The ability of the microfluidics chip platform to enrich plasma cell content by depleting CD45+ white blood cells has been demonstrated. Further studies will need to combine single-cell selection with RT-PCR to further enhance the diagnostic capabilities of this technology. This platform has the potential to be used for clinical risk stratification in multiple myeloma as well as minimal residual disease monitoring and selection of therapies to modulate the development of resistance.
Citation: Lee LX, Li SC. Hunting down the dominating subclone of cancer stem cells as a potential new therapeutic target in multiple myeloma: An artificial intelligence perspective. World J Stem Cells 2020; 12(8): 706-720
With approximately 31000 new cases of multiple myeloma (MM) diagnosed in the United States (US) per year, the impact of this incurable disease on individual patients and society as a whole is profound. The median age at diagnosis is 70 years old. All diagnoses of MM are believed to be preceded by a state of clonal expansion of plasma cells (PCs), including monoclonal gammopathy of unknown significance (MGUS) and smoldering myeloma (SM). The duration of these precursor conditions of MM has been demonstrated to be present up to 15 years prior to the diagnosis of MM.
Diagnosis and disease prognostication
The current diagnosis of MM requires a bone marrow biopsy and aspirate, which is used to enumerate plasma cell content and to characterize PCs by immunohistochemical staining, cytogenetics, and flow cytometry. Detection of cytogenetic alterations, in particular, are paramount to provide prognostication and direct therapy and have been incorporated into the standardized staging system for MM. For example, the presence of high-risk cytogenetics, including del17p, t(4,14), and t(14;16) prognosticates for survival 1/5th that of standard-risk cytogenetics. However, the identification of such cytogenetic features may be used to guide therapy such as in patients with t(4;14), who have traditionally had significantly inferior outcome may be able to have an overall survival (OS) similar to patients with standard-risk MM when treated with bortezomib-containing regimens and autologous stem cell transplantation.
CURRENT SOLUTIONS TO OVERCOME THERAPEUTIC RESISTANCE
Initial treatment incorporating conventional drugs such as Dexamethasone (Dex) effectively induces MM cell death; however, prolonged drug exposures result in the development of chemoresistance. Thus, individual patients' survival within a risk category remains variable, and the patients relapse despite achieving a “complete response,” reflecting persistent disease that cannot be detected using the currently recommended disease evaluation techniques. It is becoming apparent that static cytogenetic categories alone are not sufficient to define subclone formation and stage. Several methods are being evaluated to enhance further our ability to individualize treatment.
First, response assessment using minimal residual disease (MRD) at varying time points in a patient's disease process can further fine-tune response-adapted treatment strategies. MRD negativity at any given time point is closely correlated with more prolonged progression-free survival (PFS). It has been incorporated into the International Myeloma Working Group recommendation for response assessment, and ongoing studies are studying adaptive treatment strategies based on achieving MRD negativity. Current methods for minimal residual disease testing include flow cytometry or next-generation sequencing. Multiparametric flow cytometry (MFC) MRD testing in MM has quite low sensitivity, detecting one cell in 104 cells and requires at least 2 × 106, preferably greater than 5 × 106 bone marrow cells to be measured, as recommended by the International Clinical Cytometry Society and the European Society for Clinical Cell Analysis. In addition to low sensitivity, MFC may be unable to differentiate between the dominant clone and various subclones. In addition, there is heterogeneity between laboratories (cross-platform flow cytometry), which depends on instrumentation used and initial gating parameters (CD38, CD138, CD45, forward, and sideward light scatter), within the same aliquot and is therefore entirely subjective. Next-generation sequencing (NGS) of immunoglobulin gene sequences is an alternative method for MRD assessment. While a more sensitive technique compared with MFC, detecting one tumor cell in 106, NGS cannot detect mutations that are present within individual cells. Therefore, based on our current technologies for MRD detection, we can only say whether a patient is positive or negative without genuinely understanding the temporal and spatial heterogeneity within a given plasma cell population.
Second, we recognize the temporal and spatial heterogeneity in MM, as clinical observations revealed that several subclones of PCs exist at diagnosis and that there is selective therapeutic pressure for the evolution of individual subclones. This phenomenon can be tracked utilizing the whole-genome sequencing of paired tumor/normal samples. In one study, from 203 MM patients revealed frequent mutations in KRAS, NRAS, BRAF, FAM46C, TP53, and DIS3. Mutations were often present in subclonal populations, and multiple mutations within the same pathway (e.g., KRAS, NRAS, and BRAF) were observed in the same patient". However, a more recent study utilizing plasma samples found mutations in the KRAS-MAPK pathway in 70% of samples in addition to multiple mutations within subclones including a notable mutation in PIK3CA signifying perhaps the relative insensitivity of a one site biopsy in addition to the development of more mutated clones with the escape of PCs from the bone marrow microenvironment. Liquid biopsy utilizing cfDNA can provide us with information on targetable mutations, but is a way to study spatial and temporal heterogeneity present. The drawbacks for above current testing mistaken population phenomena for real physiological events happening only within a single-cell (i.e., subclone) – mutations exist in different cells may not cross-talk - thereby not being able to manifest as clinically treatable phenotypes – which would not give early insight into the evolution of a given patient's MM.
Third, the identification of chemoresistant biomarkers offers a trace to the subclones, e.g., the oligonucleotide array analysis demonstrates that heat shock protein-27 (Hsp27) is upregulated in Dex-resistant, but not in Dex-sensitive MM cells. Proteomics analysis of Hsp27-immunocomplexes revealed the presence of actin in Dex-resistant, but not in Dex-sensitive cells. The activator protein-1 transcription factor family (JUNB) driving the JunB-mediated phenotype in MM cells: knockdown of JUNB restored the response to dexamethasone in dexamethasone-resistant MM cells. When JunB-ER fusion protein in dexamethasone-sensitive MM cells is activated by 4-hydroxytamoxifen, Dex-sensitive cells become to be resistant to dexamethasone- and bortezomib-induced cytotoxicity.
Thus, the ability to track mutations within a single cell subclone lends to the study of mechanisms of drug resistance, possibly leading to a better selection of targeted therapies. To that end, new technology must be developed and raised its sensitivity sufficient to evaluate the low burden of MM cells, which is currently being investigated as a way to detect pre-biochemical relapse. We propose to develop a technique that combines the detection of low-frequency events combines with the in-depth characterization of the remaining subclones.
DEVELOPMENT OF AN INNOVATIVE SINGLE-CELL MOLECULAR PROFILING PLATFORM
While single-cell proteomics is still uncertain, single-cell RNA-seq is a widespread practice in research laboratories now. Many microfluidic devices, including ours, have been developed for single-cell transcriptome analysis but clinical application of single-cell transcriptome is still not common, especially in cancer characterization and classification. US Food and Drug Administration approved CellSearch™ (Janssen, Raritan, NJ) of circulating tumor cells (CTC) for predicting PFS and OS in metastatic breast cancer for clinical use in 2005. However, CellSearch® data correlate negatively with survival in patients with metastatic breast, colorectal, or prostate cancer. CellSearch™, along with several other CTC enrichment techniques, relies on only fluorescent imaging analysis.
At least 50 competitor circulating tumor cell platforms exist (Table 1)[19,20], only one, “CellSearch™” (Janssen, Raritan, NJ) was cleared by the US Food and Drug Administration in October 2005. As we have a prototype Multi-Phase Laser-cavitation Single Cell Analyzer (MLSCA) device for single-cell transcriptome analysis (Figure 1), we develop a necessary commercialization 510k, and the CAP/CLIA component of this proposal must be side by side comparison of the applicant's technology to the CellSearch™ technology. This feature differentiates it from numerous competitor platforms used for single-cell counting. CellSearch™ can be used for predicting PFS and OS in metastatic breast cancer, however; CellSearch™ cannot generate consistent results for routine clinical use, due to the limit of the sensitivity of these devices. CellSearch™ captures CTCs from blood using magnetic particles coated with anti-EpCAM (CD326; 17-1A antigen) antibodies, which relies on an antibody that binds to the protein EpCAM (epithelial cell adhesion molecule), present on the surface of malignant epithelial cells but not of blood cells. As CellSearch®, along with several other CTC enrichment techniques, relies on the presence of epithelial cell markers, CTCs that do not express EpCAM, such as those that have undergone an epithelial to mesenchymal transition may be missed[24,25].
Table 1 Multi-Phase Laser-cavitation Single Cell Analyzer can perform both circulating tumor cells enumeration and single-cell molecular characterization.
Figure 1 A prototype LSCAT device for single-cell transcriptome analysis.
Thus, CellSearch™ cannot generate consistent results for routine clinical use, and there is an unmet need to combine image analysis with RNA-seq for cancer classification. To fill this unmet need, we developed a microfluidic prototype device to connect FACS and imaging analysis (FISH) with molecular analysis (e.g., single-cell transcriptomes) - whose prototype device is called MLSCA (Figure 1).
Our MLSCA can overcome the limitation of “CellSearch™" [multiple-step processing, i.e., separated reverse transcriptase polymerase chain reaction (RT-PCR)], and the strength is the ability of our MLSCA to perform RT-PCR on a single cell isolated on the chip (one-step processing), eliminating process errors. Our microfluidic system is equipped with both single-cell isolation and cDNA synthesis capabilities. Thus, our MLSCA enables (1) microscale fluorescence-activated cell selection for separation of rare cell subpopulations; and (2) generation of the high-quality single-cell transcriptome with nano-droplets. We have conducted and published a small Phase I clinical trial with our devices in myeloma risk stratification of MM. As a proof-of-concept for instrumentation of our prototype device, we published the clinical trial using MLSCA/MF-CD45-TACs on MM 48 patients, which shed new light for scale-up applications in clinics.
Conventional single-cell isolation techniques with microliter carry-over volumes cannot be used for sensitive nano-liter RT directly. Unlike PCR, which is a repetitive event that itself may introduce bias, RT is a single biochemical reaction for which starting mRNA concentration is critical. Further, the resultant cDNA population is thought to be unbiased. The innovations of our device include: (1) The single cell analyzer (MLSCA) combines single-cell fluorescent-activated cell selection, reverse transcript synthesis of high-quality cDNA, transcriptome analysis - all in 0.1 nano-liter droplets (50 pico-liter), thereby reducing procedure errors to improve single-cell cDNA quality and to yield reliable single-cell transcriptomes down-stream; (2) The ability to create single artificial cells, a highly controllable test system, for evaluation of platform performance, using homogeneous droplets of known low abundant RNA; (3) Demonstration of such MLSCA analysis clinical specimens; (4) Development of statistical methods for defining and quantifying molecular heterogeneity in populations of cells; and (5) Linkage to an RNAseq platform (Helicos and Ion Torrent) that does not require PCR amplification of the input cDNA is crucial for ensuring the fidelity of the measured transcriptome: read counts can be directly related to RNA abundance, with no possible distortion due to differential PCR efficiencies. Thus, the strength is the ability of the MLSCA to perform RT-PCR on a single cell on the chip. This strength differentiates it from numerous competitor platforms only used for CTC counting. We tested the MLSCA by addressing MM heterogeneity to identify underlying tumor initiation, and relapse biomarkers in MM, which was inspired by that fact that single-cell transcriptomic analysis in medulloblastomas led to mapping oncogenic networks including HIPPO-YAP/TAZ and AURORA-A/MYCN pathways and to identify pathways of drug resistance.
We applied this MLSCA for a clinical trial of risk stratification of MM. As our MF-CD45-TACs differs from CellSearch™, we expected single-cell features (FACS, biomarkers) of CellSearch™ differ from those of our MF-CD45-TACs (Figure 2). Specifically, our MF-CD45-TACs-based technology distinguished CD45- cells from MM PCs, which improved the detection of rare genetic alternation in PCs, which was a significant improvement over direct flow cytometry and FISH, and led to more precise diagnosis and prognosis of MM.
Figure 2 Processing of single-cell droplets.
A: A single-cell in an oil droplet travels to a trapping module of the 10-nl reactor (green, in inset). Other cells are forced to bypass to the next unit; B: FEOET push the droplet (with a cell) past the open valve (orange bar), which closes, locking the droplet into the 10-nl ring with RT/RT-PCR master mix (green); and C: The ring’s peristaltic pump breaks the droplet to mix the cell with master-mix. When the reaction is finished, oil (blue) pushes the product (cDNA) out of the ring in the form of a droplet (10-nL) for downstream molecular analysis (Refer to for details).
This attribute is of significance as MM is an incurable neoplasm of PCs that affects more than 20000 people annually in the United States. Risk stratification, primarily based on cytogenetic abnormalities, has emerged as essential for its management. Thalidomide, lenalidomide, and pomalidomide, first to the third generation of immunomodulatory drugs (IMiDs), respectively, are used for maintenance therapy of MM. Cytogenetic alterations are the base of risk stratification for MM and the selection of which IMiDs for MM therapy.
The rarity and sporadic distribution of PCs in bone marrow often lead to false-negative results of FACS and cytogenetic detection performed directly on a bone marrow biopsy sample. Target cell enrichment could overcome the rarity and sporadic distribution of PCs in the bone marrow. Density gradient centrifugation and magnetically labeled antibodies [DG/magnetic-activated cell sorting (MACS)] with MACS have been widely used to enrich target cells in blood samples. MACS enrichment of CD138+ cells for FISH in MM diagnosis has been reported. However, MM cells with low levels of CD138 have also been associated with poor prognosis. Therefore, a better enrichment method is needed. Here, we report a novel microfluidic approach, combining microfluidic size selection and CD45 depletion with tetrameric antibody complexes (TACs) for the enrichment of MM cells (MF-CD45-TACs) in bone marrow samples. Our study showed that this approach significantly improves the detection of rare genetic alternations in PCs. Parallel diagnosis performed for 48 patients (Figure 3) showed that the microfluidic enrichment approach represents a significant improvement over direct flow cytometry and FISH and leads to more precise diagnosis and prognosis. Implementation of this modified diagnostic assay in clinics could improve the current clinical outcomes of MM (Figure 4 and Figure 5).
Figure 3 Schematic designs of the proposed workflow.
A: Consolidated Standards of Reporting Trials diagram. A total of 63 patients were screened for eligibility. Only 48 patients were newly diagnosed with multiple myeloma before receiving any treatment. These patients were enrolled, and their bone marrow obtained at diagnosis was divided into two aliquots: One aliquot underwent traditional flow cytometry and FISH analysis, and the other aliquot was subjected to microfluidic selection for enrichment of CD45-PCs, then subjected to flow cytometry and FISH analysis. Results from both methods were compared; B: Comparison of traditional method to microfluidic method (MF-CD45-TACs). MF-CD45-TACs significantly enrich plasma cells for flow cytometry and FISH assays and improve the accuracy of these assays; C: This proposed workflow (Note that we can use both bone marrow and circulating multiple myeloma cells).
Figure 4 Improved clinical outcomes with microfluidic CD45 depletion (Patient 1).
A: Bone marrow smear at the time of initial diagnosis; B: Bone marrow smear after effective treatment (complete remission); C: Flow-cytometry without microfluidic enrichment. Plasma cells (CD38+/CD138+) is only 4.24% and no FISH cytogenetic abnormalities were found (low risk); D: After microfluidic enrichment (center inset) plasma cells (CD38+/CD138+) increased to 36.1%; E: FISH on enriched plasma cells show 17p- (Red: D13S319; Green: P53); F: FISH showed del(13q14) in enriched PC (Red: D13S319; Green: RB1). With enriched plasma cell for FISH, the patient was reclassified and treated as high-risk multiple myeloma which leads to complete remission (Refer to for details).
Figure 5 Microfluidic risk-stratification improves clinical outcomes of multiple myeloma (Patient 2).
A: Bone marrow at diagnosis. Active granulocyte hyperplasia; B: Partial remission was achieved with revised risk-stratification; C: At diagnosis, bone marrow plasma cell (PC) abnormalities included clustered and scattered distribution of primitive and immature PCs, with large cell body, fine chromatin, visible nucleolus, and abundant cytoplasm; D: After treatment for high-risk multiple myeloma, PCs were rare and had normal morphology; no typical abnormal PCs were observed; E: At diagnosis without microfluidic enrichment, PCs (CD38+/CD138+) were only 1.84 %; F: After microfluidic enrichment, PCs (CD38+/CD138+) increased to 40.79%; G: Without microfluidic enrichment, FISH showed IgH rearrangement and del(13q14), leading to classification as intermediate risk (Red: D13S319; Green: P53); H and I: After microfluidic enrichment, in addition to del(13q14), FISH showed t (4,14) fusion (yellow dots) and 17p- (Red: D13S319; Green: P53), patient was reclassified and treated as high-risk, which led to efficacious treatment (Refer to for details).
APPLICATION OF SINGLE CELL SUBCLONE TRACKING TO THE FUTURE TREATMENT OF MM
At the MGUS/Smoldering stage
While all MM is preceded by an asymptomatic MGUS/smoldering myeloma stage, only a fraction of these individuals will evolve to symptomatic MM. Currently, some high-risk features such as high bone-marrow plasma cell burden, light chain ratios, and predicates the development of symptomatic MM. Still, we do not understand the oncogenesis of MM and, therefore, cannot accurately determine who will progress and who will remain asymptomatic. Using circulating tumor cell technology, we could track the occurrence of trigger genetic events in pre-symptomatic patients without the need to perform repeated invasive procedures and potentially intervene to eradicate these emerging malignant subclones using targeted therapies prior to the development of symptomatic disease or the acquisition of additional potent genetic mutations.
MRD monitoring has become one of the most relevant prognostic factors for MM. It has been shown that persistent MRD is associated with improved progression-free survival and overall survival. Not only is single cell tracking methodology as described in this article a sensitive method to detect MRD, but the characteristics of the residual cells will also be able to be elucidated. As such, one can detect the emergence of a “dangerous” clone.
Determining the sequencing of therapies including immunotherapy
Currently, we have more than 14 unique treatments for the MM, which, when used in combination, yields dozens of combination options for patients. Clinical trials using antibody drug-conjugate and bispecific antigen-directed CD3 T-cell engager targeting, by checkpoint inhibitors and an anti-T-cell immunoglobulin and ITIM domains antibody are currently underway and has the potential to further prolong survival times. Chimeric antigen receptor T cell therapy targeting B-cell maturation antigen, immunoglobulin kappa chain, SLAM family member 7, or G-protein coupled receptor family C group 5 member D, the activated integrin beta7 is a promising treatment modality which can often give long progression-free survival in heavily treated patients[35,36]. Despite this arsenal of treatments, the elusive cure for MM has yet to be found, and the current approaches to the treatment of refractory disease produce progressively short-lived efficacy. Perhaps this is because the sequencing of these treatments is often borne out of trial and error and do not take into consideration the temporal changes and spatial relationships of MM subclones. The development of chemoresistance, leading to shorter progression-free with each subsequent treatment and overall survival, lies in understanding how therapies drive the evolution of subclones. Using sensitive methods for detection and characterization of subclones, we can understand whether therapy induces molecular alterations within myeloma cells or selects for the survival of specific clones over others. Knowing how the treatments we use drive the process of evolution can allow clinicians to choose combinations of therapies that will modulate the development of chemoresistance.
We envision that single-cell technology will innovate cancer stem cell subclonal evolution on time-space landscaping of heterogeneity and imply the lineage-tracking pathway-based prediction of therapeutic efficacies of cancer treatment. Accurately, temporal development and spatial distribution of quantitative subclonal measurement of MM will reveal therapeutic sensitizing mutations, thereby moving closer to developing a therapeutic window of cancer in the advent of new, more productive, and less toxic therapies. We hypothesize that subclonal evolution, in conjunction with current standard care, will improve outcomes in patients with heterogeneous pathologies (Figure 6). Circulating MM counts and Cav-1 molecules early during radiotherapy are independently predictive of recurrence in MM. Physicians assert that every time that there is a reference to, visual or spoken, the patient view of the landscape of an MM diagnostics that they claim to predict the outcomes legitimately, it has to be as comprehensive and individualized as possible given the data package generated from AI-Med algorithms. Recently, we applied our microfluidic devices to myeloma risk stratification. However, like many current microfluidic devices, the device only enhanced and improved traditional FACS and FISH. The gap of integrating genomic profiles into cancer characterization is still not filled. Therefore, it is a logical and necessary step for us to integrate our single-cell RNA-seq technology into cancer characterization, specifically molecular classification of MM in cancer genome landscape such as the Pan-Cancer Analysis of Whole Genomes consortium. "Timing analyses suggest that driver mutations often precede diagnosis by many years, if not decades. Together, these results determine the evolutionary trajectories of cancer and highlight opportunities for early cancer detection". All of these must rely on "A deep learning system accurately classifies primary and metastatic cancers using passenger mutation patterns" for integration of dynamic space-time changes. One such integrated platform was to scaffold the diverse datasets together, allowing them to interface not only across single-cell transcriptomics (scRNA-seq), but also across distinct cellular modalities – e.g., a bone marrow atlas to characterize lymphocyte populations – to better understand cellular identity and function beyond the taxonomic listing of clusters of cellular heterogeneity.
Figure 6 Blockade of the dominating subclonal switchboard signals in cancer stem cells as a new therapeutic strategy to suppress the dominating subclone shift to control cancer progression and post-treatment cancer recurrence.
Showed is the proposed new treatment paradigm that should target the subclonal-switchboard signals (SSS). Blocking the dominating subclonal SSS leads to subclonal quiescence, so keeping tumors alive but small and manageable (dormant/quiescent subclone). Note that SSS as mechanisms for leading to shifting dominating subclones as triggered by environmental cues (stress) for cancer progression and post-treatment. A cancer subclone may gain a mutation that, in the appropriate environment cue, leads to dominating subclonal activation due to positive selection. Showed lettering and lines/ arrows in the black color is the current concept of a treatment strategy for cancer- dominant subclonal cells (cancer stem cells) that may acquire a mutation in a suitable environment, triggering to dominating subclonal expansion and growth. When this dominating subclone is explicitly destroyed, it sends out dominating subclonal-SSS to a dormant/quiescent subclonal cell, which gets activated for dominating subclonal expansion and growth (adopted from).
Implementation of this modified diagnostic device in clinics proven to improve clinical outcomes of MM (Figure 3C). Our microfluidic-assisted stratification of single cancer cells may help understand the mechanisms underlying the temporal and spatial heterogeneities in solid tumors like brain cancer as testing is underway[44-46], thereby holding promise for using the single-cell analysis to guide treatment for targeted therapy (Table 2), as governed by artificial intelligence-based integration of genome, epigenome, and pathological measurements.
Table 2 Therapeutics targets and corresponding agents in multiple myeloma and artificial intelligence medicine.
Mechanism of action (clinical phenotype)
Resistance to chemotherapy in MM
Bcl-2/Bcl-X(L)/Bcl-w (antiapoptotic proteins)
Inhibitor ABT-737 (with bortezomib-, dexamethasone-(Dex) and thalidomide)
Preventive medicine using pulse oximetry screening
Pulse transit time (PTT) is the time it takes a pulse wave to travel between two arterial sites R-wave-gated photo-plethysmography as of measurement of PTT as a surrogate for intra-thoracic pressure changes in obstructive sleep apnea)
Respiratory sinus arrhythmia (RSA) with an algorithm for quantifying instantaneous RSA as applied to heartbeat interval and respiratory recordings to track changes in cardiorespiratory interactions elicited during meditation, otherwise not evidenced in control resting states)
Cited Literature. MM: Multiple myeloma; AP-1: Activator protein-1; ICU: Intensive care unit.
We thank Maria Minon, MD; Brent A Dethlefs; Mustafa H Kabeer, MD; William G Loudon, MD, PhD; Leonard S Sender, MD; Anthony Christopher Chang, MD, MBA, MPH; Edward Nelson, MD; Richard A van Etten, MD, PhD; Dan Cooper, MD; and Jiang F Zhong, PhD; for their support and enthusiasm.
P-Reviewer: Lin JM, Vatansever M S-Editor: Wang JL L-Editor: A P-Editor: Wu YXJ
Kazandjian D. Multiple myeloma epidemiology and survival: A unique malignancy.Semin Oncol. 2016;43:676-681.
Weiss BM, Abadie J, Verma P, Howard RS, Kuehl WM. A monoclonal gammopathy precedes multiple myeloma in most patients.Blood. 2009;113:5418-5422.
Fonseca R, Bergsagel PL, Drach J, Shaughnessy J, Gutierrez N, Stewart AK, Morgan G, Van Ness B, Chesi M, Minvielle S, Neri A, Barlogie B, Kuehl WM, Liebisch P, Davies F, Chen-Kiang S, Durie BG, Carrasco R, Sezer O, Reiman T, Pilarski L, Avet-Loiseau H; International Myeloma Working Group. International Myeloma Working Group molecular classification of multiple myeloma: spotlight review.Leukemia. 2009;23:2210-2221.
Hebraud B, Magrangeas F, Cleynen A, Lauwers-Cances V, Chretien ML, Hulin C, Leleu X, Yon E, Marit G, Karlin L, Roussel M, Stoppa AM, Belhadj K, Voillat L, Garderet L, Macro M, Caillot D, Mohty M, Facon T, Moreau P, Attal M, Munshi N, Corre J, Minvielle S, Avet-Loiseau H. Role of additional chromosomal changes in the prognostic value of t(4;14) and del(17p) in multiple myeloma: the IFM experience.Blood. 2015;125:2095-2100.
Siegel DS, Dimopoulos M, Jagannath S, Goldschmidt H, Durrant S, Kaufman JL, Leleu X, Nagler A, Offner F, Graef T, Eid JE, Houp J, Gause C, Vuocolo S, Anderson KC. VANTAGE 095: An International, Multicenter, Open-Label Study of Vorinostat (MK-0683) in Combination With Bortezomib in Patients With Relapsed and Refractory Multiple Myeloma.Clin Lymphoma Myeloma Leuk. 2016;16:329-334.e1.
Bochtler T, Merz M, Hielscher T, Granzow M, Hoffmann K, Krämer A, Raab MS, Hillengass J, Seckinger A, Kimmich C, Dittrich T, Müller-Tidow C, Hose D, Goldschmidt H, Hegenbart U, Jauch A, Schönland SO. Cytogenetic intraclonal heterogeneity of plasma cell dyscrasia in AL amyloidosis as compared with multiple myeloma.Blood Adv. 2018;2:2607-2618.
Kumar S, Paiva B, Anderson KC, Durie B, Landgren O, Moreau P, Munshi N, Lonial S, Bladé J, Mateos MV, Dimopoulos M, Kastritis E, Boccadoro M, Orlowski R, Goldschmidt H, Spencer A, Hou J, Chng WJ, Usmani SZ, Zamagni E, Shimizu K, Jagannath S, Johnsen HE, Terpos E, Reiman A, Kyle RA, Sonneveld P, Richardson PG, McCarthy P, Ludwig H, Chen W, Cavo M, Harousseau JL, Lentzsch S, Hillengass J, Palumbo A, Orfao A, Rajkumar SV, Miguel JS, Avet-Loiseau H. International Myeloma Working Group consensus criteria for response and minimal residual disease assessment in multiple myeloma.Lancet Oncol. 2016;17:e328-e346.
Soh KT, Tario JD, Wallace PK. Diagnosis of Plasma Cell Dyscrasias and Monitoring of Minimal Residual Disease by Multiparametric Flow Cytometry.Clin Lab Med. 2017;37:821-853.
Berger N, Kim-Schulze S, Parekh S. Minimal Residual Disease in Multiple Myeloma: Impact on Response Assessment, Prognosis and Tumor Heterogeneity.Adv Exp Med Biol. 2018;1100:141-159.
Arroz M, Came N, Lin P, Chen W, Yuan C, Lagoo A, Monreal M, de Tute R, Vergilio JA, Rawstron AC, Paiva B. Consensus guidelines on plasma cell myeloma minimal residual disease analysis and reporting.Cytometry B Clin Cytom. 2016;90:31-39.
Bai Y, Orfao A, Chim CS. Molecular detection of minimal residual disease in multiple myeloma.Br J Haematol. 2018;181:11-26.
Lohr JG, Stojanov P, Carter SL, Cruz-Gordillo P, Lawrence MS, Auclair D, Sougnez C, Knoechel B, Gould J, Saksena G, Cibulskis K, McKenna A, Chapman MA, Straussman R, Levy J, Perkins LM, Keats JJ, Schumacher SE, Rosenberg M; Multiple Myeloma Research Consortium, Getz G, Golub TR. Widespread genetic heterogeneity in multiple myeloma: implications for targeted therapy.Cancer Cell. 2014;25:91-101.
Hocking J, Mithraprabhu S, Kalff A, Spencer A. Liquid biopsies for liquid tumors: emerging potential of circulating free nucleic acid evaluation for the management of hematologic malignancies.Cancer Biol Med. 2016;13:215-225.
Terpos E, Suzan F, Goldschmidt H. Going the distance: Are we losing patients along the multiple myeloma treatment pathway?Crit Rev Oncol Hematol. 2018;126:19-23.
Fan F, Bashari MH, Morelli E, Tonon G, Malvestiti S, Vallet S, Jarahian M, Seckinger A, Hose D, Bakiri L, Sun C, Hu Y, Ball CR, Glimm H, Sattler M, Goldschmidt H, Wagner EF, Tassone P, Jaeger D, Podar K. The AP-1 transcription factor JunB is essential for multiple myeloma cell proliferation and drug resistance in the bone marrow microenvironment.Leukemia. 2017;31:1570-1581.
He R, Yang N, Zhang P, Liu J, Li J, Zhou F, Zhang W. Identification and expression of MMSA-8, and its clinical significance in multiple myeloma.Oncol Rep. 2017;37:3235-3243.
Zeng Y, Gao L, Luo X, Chen Y, Kabeer MH, Chen X, Stucky A, Loudon WG, Li SC, Zhang X, Zhong JF. Microfluidic enrichment of plasma cells improves treatment of multiple myeloma.Mol Oncol. 2018;12:1004-1011.
Miller MC, Doyle GV, Terstappen LW. Significance of Circulating Tumor Cells Detected by the CellSearch System in Patients with Metastatic Breast Colorectal and Prostate Cancer.J Oncol. 2010;2010:617421.
Lianidou ES, Markou A. Circulating tumor cells in breast cancer: detection systems, molecular characterization, and future challenges.Clin Chem. 2011;57:1242-1255.
Parkinson DR, Dracopoli N, Petty BG, Compton C, Cristofanilli M, Deisseroth A, Hayes DF, Kapke G, Kumar P, Lee JSh, Liu MC, McCormack R. Mikulski S, Nagahara L, Pantel K, Pearson-White S, Punnoose EA, Roadcap LT, Schade AE, Scher HI, Sigman CC, Kelloff GJ. Considerations in the development of circulating tumor cell technology for clinical use.J Transl Med. 2012;10:138.
Liu MC. By the numbers: does circulating tumor cell enumeration have a role in metastatic breast cancer?J Clin Oncol. 2014;32:3479-3482.
Münz M, Murr A, Kvesic M, Rau D, Mangold S, Pflanz S, Lumsden J, Volkland J, Fagerberg J, Riethmüller G, Rüttinger D, Kufer P, Baeuerle PA, Raum T. Side-by-side analysis of five clinically tested anti-EpCAM monoclonal antibodies.Cancer Cell Int. 2010;10:44.
Gorges TM, Tinhofer I, Drosch M, Röse L, Zollner TM, Krahn T, von Ahsen O. Circulating tumour cells escape from EpCAM-based detection due to epithelial-to-mesenchymal transition.BMC Cancer. 2012;12:178.
Königsberg R, Obermayr E, Bises G, Pfeiler G, Gneist M, Wrba F, de Santis M, Zeillinger R, Hudec M, Dittrich C. Detection of EpCAM positive and negative circulating tumor cells in metastatic breast cancer patients.Acta Oncol. 2011;50:700-710.
Chen X, Wen Q, Stucky A, Zeng Y, Gao S, Loudon WG, Ho HW, Kabeer MH, Li SC, Zhang X, Zhong JF. Relapse pathway of glioblastoma revealed by single-cell molecular analysis.Carcinogenesis. 2018;39:931-936.
Zhang L, He X, Liu X, Zhang F, Huang LF, Potter AS, Xu L, Zhou W, Zheng T, Luo Z, Berry KP, Pribnow A, Smith SM, Fuller C, Jones BV, Fouladi M, Drissi R, Yang ZJ, Gustafson WC, Remke M, Pomeroy SL, Girard EJ, Olson JM, Morrissy AS, Vladoiu MC, Zhang J, Tian W, Xin M, Taylor MD, Potter SS, Roussel MF, Weiss WA, Lu QR. Single-Cell Transcriptomics in Medulloblastoma Reveals Tumor-Initiating Progenitors and Oncogenic Cascades during Tumorigenesis and Relapse.Cancer Cell. 2019;36:302-318.e7.
Bertrand KC, Faria CC, Skowron P, Luck A, Garzia L, Wu X, Agnihotri S, Smith CA, Taylor MD, Mack SC, Rutka JT. A functional genomics approach to identify pathways of drug resistance in medulloblastoma.Acta Neuropathol Commun. 2018;6:146.
Mikhael JR, Dingli D, Roy V, Reeder CB, Buadi FK, Hayman SR, Dispenzieri A, Fonseca R, Sher T, Kyle RA, Lin Y, Russell SJ, Kumar S, Bergsagel PL, Zeldenrust SR, Leung N, Drake MT, Kapoor P, Ansell SM, Witzig TE, Lust JA, Dalton RJ, Gertz MA, Stewart AK, Rajkumar SV, Chanan-Khan A, Lacy MQ; Mayo Clinic. Management of newly diagnosed symptomatic multiple myeloma: updated Mayo Stratification of Myeloma and Risk-Adapted Therapy (mSMART) consensus guidelines 2013.Mayo Clin Proc. 2013;88:360-376.
Nathwani N, Larsen JT, Kapoor P. Consolidation and Maintenance Therapies for Newly Diagnosed Multiple Myeloma in the Era of Novel Agents.Curr Hematol Malig Rep. 2016;11:127-136.
Chen L, Li J, Xu W, Qiu H, Zhu Y, Zhang Y, Duan L, Qian S, Lu H. Molecular cytogenetic aberrations in patients with multiple myeloma studied by interphase fluorescence in situ hybridization.Exp Oncol. 2007;29:116-120.
Kawano Y, Fujiwara S, Wada N, Izaki M, Yuki H, Okuno Y, Iyama K, Yamasaki H, Sakai A, Mitsuya H, Hata H. Multiple myeloma cells expressing low levels of CD138 have an immature phenotype and reduced sensitivity to lenalidomide.Int J Oncol. 2012;41:876-884.
Landgren O, Kyle RA, Pfeiffer RM, Katzmann JA, Caporaso NE, Hayes RB, Dispenzieri A, Kumar S, Clark RJ, Baris D, Hoover R, Rajkumar SV. Monoclonal gammopathy of undetermined significance (MGUS) consistently precedes multiple myeloma: a prospective study.Blood. 2009;113:5412-5417.
Tamura H, Ishibashi M, Sunakawa M, Inokuchi K. Immunotherapy for Multiple Myeloma.Cancers (Basel). 2019;11:2009.
Hosen N. Chimeric Antigen Receptor T-Cell Therapy for Multiple Myeloma.Cancers (Basel). 2019;11:2024.
Hosen N. Chimeric antigen receptor T-cell therapy for multiple myeloma.Int J Hematol. 2020;111:530-534.
Li SC, Han YP, Dethlefs BA, Loudon WG. Therapeutic window, a critical developmental stage for stem cell therapies.Curr Stem Cell Res Ther. 2010;5:297-293.
Li SC, Lee KL, Luo J. Control dominating subclones for managing cancer progression and posttreatment recurrence by subclonal switchboard signal: implication for new therapies.Stem Cells Dev. 2012;21:503-506.
ICGC/TCGA Pan-Cancer Analysis of Whole Genomes Consortium. Pan-cancer analysis of whole genomes.Nature. 2020;578:82-93.
Gerstung M, Jolly C, Leshchiner I, Dentro SC, Gonzalez S, Rosebrock D, Mitchell TJ, Rubanova Y, Anur P, Yu K, Tarabichi M, Deshwar A, Wintersinger J, Kleinheinz K, Vázquez-García I, Haase K, Jerman L, Sengupta S, Macintyre G, Malikic S, Donmez N, Livitz DG, Cmero M, Demeulemeester J, Schumacher S, Fan Y, Yao X, Lee J, Schlesner M, Boutros PC, Bowtell DD, Zhu H, Getz G, Imielinski M, Beroukhim R, Sahinalp SC, Ji Y, Peifer M, Markowetz F, Mustonen V, Yuan K, Wang W, Morris QD; PCAWG Evolution & Heterogeneity Working Group, Spellman PT, Wedge DC, Van Loo P; PCAWG Consortium. The evolutionary history of 2,658 cancers.Nature. 2020;578:122-128.
Jiao W, Atwal G, Polak P, Karlic R, Cuppen E; PCAWG Tumor Subtypes and Clinical Translation Working Group, Danyi A, de Ridder J, van Herpen C, Lolkema MP, Steeghs N, Getz G, Morris Q, Stein LD; PCAWG Consortium. A deep learning system accurately classifies primary and metastatic cancers using passenger mutation patterns.Nat Commun. 2020;11:728.
Stuart T, Butler A, Hoffman P, Hafemeister C, Papalexi E, Mauck WM, Hao Y, Stoeckius M, Smibert P, Satija R. Comprehensive Integration of Single-Cell Data.Cell. 2019;177:1888-1902.e21.
Cao J, Packer JS, Ramani V, Cusanovich DA, Huynh C, Daza R, Qiu X, Lee C, Furlan SN, Steemers FJ, Adey A, Waterston RH, Trapnell C, Shendure J. Comprehensive single-cell transcriptional profiling of a multicellular organism.Science. 2017;357:661-667.
Li Z, Hao P, Wu Q, Li F, Zhao J, Wu K, Qu C, Chen Y, Li M, Chen X, Stucky A, Zhong J, Li L, Zhong JF. Genetic mutations associated with metastatic clear cell renal cell carcinoma.Oncotarget. 2016;7:16172-16179.
Chen X, Chakravarty T, Zhang Y, Li X, Zhong JF, Wang C. Single-cell transcriptome and epigenomic reprogramming of cardiomyocyte-derived cardiac progenitor cells.Sci Data. 2016;3:160079.
Zhang Y, Zhong JF, Qiu H, MacLellan WR, Marbán E, Wang C. Epigenomic Reprogramming of Adult Cardiomyocyte-Derived Cardiac Progenitor Cells.Sci Rep. 2015;5:17686.
Nagrath S, Sequist LV, Maheswaran S, Bell DW, Irimia D, Ulkus L, Smith MR, Kwak EL, Digumarthy S, Muzikansky A, Ryan P, Balis UJ, Tompkins RG, Haber DA, Toner M. Isolation of rare circulating tumour cells in cancer patients by microchip technology.Nature. 2007;450:1235-1239.
Stott SL, Hsu CH, Tsukrov DI, Yu M, Miyamoto DT, Waltman BA, Rothenberg SM, Shah AM, Smas ME, Korir GK, Floyd FP, Gilman AJ, Lord JB, Winokur D, Springer S, Irimia D, Nagrath S, Sequist LV, Lee RJ, Isselbacher KJ, Maheswaran S, Haber DA, Toner M. Isolation of circulating tumor cells using a microvortex-generating herringbone-chip.Proc Natl Acad Sci USA. 2010;107:18392-18397.
Zheng S, Lin H, Liu JQ, Balic M, Datar R, Cote RJ, Tai YC. Membrane microfilter device for selective capture, electrolysis and genomic analysis of human circulating tumor cells.J Chromatogr A. 2007;1162:154-161.
Simpson SJ, Vachula M, Kennedy MJ, Kaizer H, Coon JS, Ghalie R, Williams S, Van Epps D. Detection of tumor cells in the bone marrow, peripheral blood, and apheresis products of breast cancer patients using flow cytometry.Exp Hematol. 1995;23:1062-1068.
He W, Wang H, Hartmann LC, Cheng JX, Low PS. In vivo quantitation of rare circulating tumor cells by multiphoton intravital flow cytometry.Proc Natl Acad Sci USA. 2007;104:11760-11765.
Pachmann K, Clement JH, Schneider CP, Willen B, Camara O, Pachmann U, Höffken K. Standardized quantification of circulating peripheral tumor cells from lung and breast cancer.Clin Chem Lab Med. 2005;43:617-627.
Weight RM, Dale PS, Viator JA. Detection of circulating melanoma cells in human blood using photoacoustic flowmetry.Conf Proc IEEE Eng Med Biol Soc. 2009;2009:106-109.
Marrinucci D, Bethel K, Bruce RH, Curry DN, Hsieh B, Humphrey M, Krivacic RT, Kroener J, Kroener L, Ladanyi A, Lazarus NH, Nieva J, Kuhn P. Case study of the morphologic variation of circulating tumor cells.Hum Pathol. 2007;38:514-519.
Attard G, Swennenhuis JF, Olmos D, Reid AH, Vickers E, A'Hern R, Levink R, Coumans F, Moreira J, Riisnaes R, Oommen NB, Hawche G, Jameson C, Thompson E, Sipkema R, Carden CP, Parker C, Dearnaley D, Kaye SB, Cooper CS, Molina A, Cox ME, Terstappen LW, de Bono JS. Characterization of ERG, AR and PTEN gene status in circulating tumor cells from patients with castration-resistant prostate cancer.Cancer Res. 2009;69:2912-2918.
Yang L, Lang JC, Balasubramanian P, Jatana KR, Schuller D, Agrawal A, Zborowski M, Chalmers JJ. Optimization of an enrichment process for circulating tumor cells from the blood of head and neck cancer patients through depletion of normal cells.Biotechnol Bioeng. 2009;102:521-534.
Alix-Panabières C, Vendrell JP, Pellé O, Rebillard X, Riethdorf S, Müller V, Fabbro M, Pantel K. Detection and characterization of putative metastatic precursor cells in cancer patients.Clin Chem. 2007;53:537-539.
Müller V, Stahmann N, Riethdorf S, Rau T, Zabel T, Goetz A, Jänicke F, Pantel K. Circulating tumor cells in breast cancer: correlation to bone marrow micrometastases, heterogeneous response to systemic therapy and low proliferative activity.Clin Cancer Res. 2005;11:3678-3685.
Gascoyne PR, Noshari J, Anderson TJ, Becker FF. Isolation of rare cells from cell mixtures by dielectrophoresis.Electrophoresis. 2009;30:1388-1398.
Lu J, Fan T, Zhao Q, Zeng W, Zaslavsky E, Chen JJ, Frohman MA, Golightly MG, Madajewicz S, Chen WT. Isolation of circulating epithelial and tumor progenitor cells with an invasive phenotype from breast cancer patients.Int J Cancer. 2010;126:669-683.
Leary RJ, Kinde I, Diehl F, Schmidt K, Clouser C, Duncan C, Antipova A, Lee C, McKernan K, De La Vega FM, Kinzler KW, Vogelstein B, Diaz LA, Velculescu VE. Development of personalized tumor biomarkers using massively parallel sequencing.Sci Transl Med. 2010;2:20ra14.
Xi L, Nicastri DG, El-Hefnawy T, Hughes SJ, Luketich JD, Godfrey TE. Optimal markers for real-time quantitative reverse transcription PCR detection of circulating tumor cells from melanoma, breast, colon, esophageal, head and neck, and lung cancers.Clin Chem. 2007;53:1206-1215.
Wu CH, Lin SR, Yu FJ, Wu DC, Pan YS, Hsieh JS, Huang SY, Wang JY. Development of a high-throughput membrane-array method for molecular diagnosis of circulating tumor cells in patients with gastric cancers.Int J Cancer. 2006;119:373-379.
Chauhan D, Velankar M, Brahmandam M, Hideshima T, Podar K, Richardson P, Schlossman R, Ghobrial I, Raje N, Munshi N, Anderson KC. A novel Bcl-2/Bcl-X(L)/Bcl-w inhibitor ABT-737 as therapy in multiple myeloma.Oncogene. 2007;26:2374-2380.
Chauhan D, Li G, Auclair D, Hideshima T, Podar K, Mitsiades N, Mitsiades C, Chen LB, Munshi N, Saxena S, Anderson KC. 2-Methoxyestardiol and bortezomib/proteasome-inhibitor overcome dexamethasone-resistance in multiple myeloma cells by modulating Heat Shock Protein-27.Apoptosis. 2004;9:149-155.
Raje N, Hideshima T, Mukherjee S, Raab M, Vallet S, Chhetri S, Cirstea D, Pozzi S, Mitsiades C, Rooney M, Kiziltepe T, Podar K, Okawa Y, Ikeda H, Carrasco R, Richardson PG, Chauhan D, Munshi NC, Sharma S, Parikh H, Chabner B, Scadden D, Anderson KC. Preclinical activity of P276-00, a novel small-molecule cyclin-dependent kinase inhibitor in the therapy of multiple myeloma.Leukemia. 2009;23:961-970.
Chen L, Ju SG, Wang ZY, Li J, Yuan YQ, Fu JX. [Sensitivity to bortezomib of RPIM8226 cells after co-cultured with down-regulated Cav-1 expression HUVECs].Zhonghua Xue Ye Xue Za Zhi. 2013;34:946-951.
Krug B, Kugel H, Harnischmacher U, Heindel W, Schmidt R, Krings F. MR pulsatility measurements in peripheral arteries: preliminary results.Magn Reson Med. 1995;34:698-705.
Naschitz JE, Bezobchuk S, Mussafia-Priselac R, Sundick S, Dreyfuss D, Khorshidi I, Karidis A, Manor H, Nagar M, Peck ER, Peck S, Storch S, Rosner I, Gaitini L. Pulse transit time by R-wave-gated infrared photoplethysmography: review of the literature and personal experience.J Clin Monit Comput. 2004;18:333-342.
Oddie S, McGuire W. Response to the Letter "RE: Commentary on 'Pulse Oximetry Screening for Critical Congenital Heart Defects'".Neonatology. 2019;116:392.
Kodituwakku S, Lazar SW, Indic P, Chen Z, Brown EN, Barbieri R. Point process time-frequency analysis of dynamic respiratory patterns during meditation practice.Med Biol Eng Comput. 2012;50:261-275.
Kodituwakku S, Lazar SW, Indic P, Brown EN, Barbieri R. Point process time-frequency analysis of respiratory sinus arrhythmia under altered respiration dynamics.Conf Proc IEEE Eng Med Biol Soc. 2010;2010:1622-1625.
Ye C, Murano E, Stone M, Prince JL. A Bayesian approach to distinguishing interdigitated tongue muscles from limited diffusion magnetic resonance imaging.Comput Med Imaging Graph. 2015;45:63-74.
Marmar CR, Brown AD, Qian M, Laska E, Siegel C, Li M, Abu-Amara D, Tsiartas A, Richey C, Smith J, Knoth B, Vergyri D. Speech-based markers for posttraumatic stress disorder in US veterans.Depress Anxiety. 2019;36:607-616.
Chen C, Wen T, Liao W. Neurally adjusted ventilatory assist versus pressure support ventilation in patient-ventilator interaction and clinical outcomes: a meta-analysis of clinical trials.Ann Transl Med. 2019;7:382.
Swennenhuis JF, van Dalum G, Zeune LL, Terstappen LW. Improving the CellSearch® system.Expert Rev Mol Diagn. 2016;16:1291-1305.