Hung-Ching Chang1
Antony M. Gitau2
Siri Kothapalli3
Danny R. Welch4,5
Mihaela E. Sardiu5,6,7
Matthew D. McCoy8,9*- 1Department of Biostatistics, University of Pittsburgh, Pittsburgh, PA, United States
- 2Department of Electrical and Electronics Engineering, Kenyatta University, Nairobi, Kenya
- 3Department of Engineering and Computer Science, Baylor University, Waco, TX, United States
- 4Department of Cancer Biology, University of Kansas Medical Center, Kansas City, KS, United States
- 5The University of Kansas Cancer Center, Kansas City, KS, United States
- 6Department of Biostatistics and Data Science, University of Kansas Medical Center, Kansas City, KS, United States
- 7Kansas Institute for Precision Medicine, University of Kansas Medical Center, Kansas City, KS, United States
- 8Innovation Center for Biomedical Informatics, Department of Oncology, Georgetown University Medical Center, Washington, DC, United States
- 9Lombardi Comprehensive Cancer Center, Washington, DC, United States
Digital twins are made of a real-world component where data is measured and a virtual component where those measurements are used to parameterize computational models. There is growing interest in applying digital twins-based approaches to optimize personalized treatment plans and improve health outcomes. The integration of artificial intelligence is critical in this process, as it enables the development of sophisticated disease models that can accurately predict patient response to therapeutic interventions. There is a unique and equally important application of AI to the real-world component of a digital twin when it is applied to medical interventions. The patient can only be treated once, and therefore, we must turn to the experience and outcomes of previously treated patients for validation and optimization of the computational predictions. The physical component of a digital twins instead must utilize a compilation of available data from previously treated cancer patients whose characteristics (genetics, tumor type, lifestyle, etc.) closely parallel those of a newly diagnosed cancer patient for the purpose of predicting outcomes, stratifying treatment options, predicting responses to treatment and/or adverse events. These tasks include the development of robust data collection methods, ensuring data availability, creating precise and dependable models, and establishing ethical guidelines for the use and sharing of data. To successfully implement digital twin technology in clinical care, it is crucial to gather data that accurately reflects the variety of diseases and the diversity of the population.
1. Introduction
Digital Twins (DT) have gained widespread adoption as a crucial component of the industrial engineering approach to optimize manufacturing processes and drive advancements in modern industrial design. However, the transformative potential of DT extends beyond the industrial realm and finds numerous impactful applications in healthcare (Madhavan et al., 2021; Stahlberg et al., 2022). In particular, DT have emerged as powerful tools in precision medicine, clinical trials, and drug development (Haleem et al., 2023; Sun et al., 2023). A notable example is Unlearn. AI, a pioneering startup that secured an impressive $17 million in funding to develop DT specifically for trials. In a groundbreaking study by Venkatesh et al. (2022), Unlearn. AI predicted that Health Digital Twins (HDT) could potentially reduce clinical trial expenses by a staggering 25%. This groundbreaking finding underscores the immense value and potential of DT in revolutionizing the healthcare landscape and highlights their ability to drive substantial cost savings while enhancing patient outcomes.
In its most general form, a DT consists of a real-world object and its virtual counterpart, which work together to enhance each other through an iterative process. Data is collected under different conditions to refine the computational simulations, and the virtual representation is used to predict how to optimize the real-world system. By testing these predictions, we can identify their success and failure, and further improve the virtual model. As a result, through each DT iteration, the real-world system is continuously updated and improved.
In order to avoid altering the state of a real-world system, successful DT applications often rely on having identical copies of the real-world counterpart available for testing the predicted response. This allows for testing multiple, potentially destructive conditions, and, unfortunately, does not apply to HDT. In the absence of identical copies of an individual, we must look to other patients when training predictive models and evaluating the predicted outcomes of various therapeutic interventions. These patients are closely matched with respect to their biological characteristics, and longitudinal data detailing their intervention and response, and the outcome can serve to validate model predictions.
The incorporation of artificial intelligence (AI) plays a pivotal role in the realm of DT and can be dissected into three key components: (1) identify patient data from individuals that closely mimic the biological system of the target patient to evaluate predicted outcomes; (2) predict the patient outcome to optimize an individualized therapeutic approach and reduce the cost associated with healthcare interventions; (3) unravel the underlying biological mechanisms that drive diseases. It is worth noting that components (2) and (3) would greatly benefit from an efficient method for identifying real-world DT, as the increased information on patient disease and response would expedite the development of disease models and drive further understanding disease mechanisms.
To collect real-world datasets that closely align with an individual patient’s DT, we present a novel approach that capitalizes on the public availability of high throughput “omics” data and represents an early iteration of a workflow for identifying matched datasets to incorporate into patient DT. In the 21st century, many studies (McClellan and King, 2010; Tam et al., 2019) propose that human disease can be characterized by marked genetic heterogeneity. While some studies apply clinical feature maps to represent a real-world patient virtually (Venkatesh et al., 2022), we propose a pioneering concept that characterizes clinical and molecular features (Figure 1) from The Cancer Genome Atlas (Weinstein et al., 2013) and identifies a subset of individuals to incorporate into a given patient’s DT. The outcomes and experience of this subset of “like” individuals can then be used to evaluate the robustness of the computational predictions.
Figure 1. Pipeline of patient DT identification. The pipeline is divided into three stages: (i) Data Integration for integrating enriched omics and clinical data, (ii) Digital Twins identification based on Euclidean distance in UMAP space, and (iii) Prediction of clinical outcomes for better patient stratification.
Our methodology begins by constructing a low-dimensional embedding of clinical and molecular features, enriched by an external annotation database. This low-dimensional feature map serves as the foundation for identifying a subset of individuals most representative of the DT through Euclidean similarity, and the results are visually presented using UMAP. Consequently, these embeddings offer a better representation of each patient by projecting those with similar genetic patterns and clinical features to the same region of the plot. The selected DT have the potential to predict the survival outcomes of new patients and even forecast drug responses. Furthermore, this pipeline can be extended to incorporate multi-omics data and environmental factors, enhancing its adaptability and comprehensive analytical capabilities.
2. Challenges to data collection
When collecting biological or clinical data for DT, there are various challenges that need to be addressed to ensure accurate and valuable representations of individuals. DT data is collected at various scales, capturing different types of information related to lifestyle, clinical and routine diagnostic data, as well as research data. For instance, they can gather data on individual lifestyle factors like physical activity, sleep patterns, nutrition, and stress levels, using wearable devices, mobile apps, or self-reporting. Additionally, they can collect clinical data on medication usage, laboratory test results, and vital signs, which is typically sourced from electronic health records (EHRs), health monitoring devices, or patient self-reporting. Processing and analyzing large amounts of data in digital twin computational models requires considerable computational power. To effectively handle the computational demands of these models, it is crucial to have a high-performance computing infrastructure, consisting of powerful processors, memory, and storage systems.
In addition to the many challenges associated with the data collection, the fact that the data comes from human patients presents additional considerations. These include data accessibility, data quality, standardization, data complexity, and ethical and privacy considerations. For example, getting the required biological data can be difficult because of privacy regulations, restricted entry to medical records or research databases, and the need for agreement from participants. It is essential to access diverse and representative datasets to create thorough and widely applicable DT. Furthermore, collecting and utilizing biological data for DT raises ethical and privacy concerns, such as informed consent, data anonymization, and data sharing. Ensuring compliance with ethical guidelines and implementing robust data security measures are essential to protect individual privacy and maintain public trust.
To tackle these obstacles, researchers, clinicians, data scientists, and policymakers need to work together. They must establish consistent protocols, create innovative analytics tools, and navigate ethical concerns. By overcoming these challenges, we can achieve more precise and trustworthy DT to benefit personalized medicine and healthcare.
3. Modeling for breast cancer
When applied to cancer patients, DT offer a forward-thinking approach to healthcare by creating virtual representations of patients. These representations use mathematical and computational models to simulate personalized diagnosis and treatment options and are increasingly integrating AI in several ways. In addition to directly learning the relationships between high throughput omics and clinical datasets and outcomes, it also serves a role in mapping patient data directly to model parameters. By utilizing DT, healthcare professionals can analyze the potential effects of different treatment approaches, leading to more informed decision-making and personalized patient care (Hernandez-Boussard et al., 2021).
Applying the digital twin paradigm to human subjects creates yet another application of AI, identifying previously treated patients that can bolster the data in the real-world component of the DT. It is impossible to test and validate multiple therapeutic interventions, and data on how patients who are matched on both their individual and tumor biologies must be used to evaluate the robustness of the virtual component of the DT.
Shown in Figure 1, a workflow was developed as part of a community hackathon event with the goal of utilizing public data (specifically, the BRCA (Breast Invasive Carcinoma) dataset from The Cancer Genome Atlas (TCGA), which provides genomic, histological, and clinical data on breast cancer patients) to build a real-world cohort of patients from a query patient. To identify additional patient data as input to a patient’s DT, the data analysis process is divided into three stages: (a) Data Integration, (b) DT identification, and (c) Prediction of Clinical Outcomes. Our objective is to conduct an integrative analysis of BRCA data utilizing omics techniques, such as enriched variant scoring (Ioannidis et al., 2016), and clinical features, such as histological subtypes for outcome prediction and DT identification (Liu et al., 2018; Thennavan et al., 2021). The BRCA dataset from TCGA consists of 970 samples and 12,908 variants. We utilize annotation files from REVEL to enhance our analysis by assessing the likelihood of a variant being pathogenic. Our final feature matrix is enriched with variant scores and supported clinical features such as PAM50, ER, PR, HER2, and triple negative subtyping data. We then utilized the UMAP clustering approach to identify real world examples that closely matched data from the input DT. This involved visualizing the entire population and calculating similarity by measuring the Euclidean distance between new patients and each sample using the first two low-dimensional representations in UMAP space. We selected k (default 50) samples with the smallest dissimilarity score as the set of DT and visualized them in the scatter plot. Through our research, we have discovered that DT can capture survival heterogeneity, making them a valuable resource for predicting the clinical outcome of new patients. It is important to note that the various sets of DT may have different survival outcomes, either better or worse than the overall population. Our survival data is obtained from UCSC Xena (Goldman et al., 2020).
Our approach can benefit from incorporating various other data types and computational modelling, including data that can be utilized for drug response prediction and prioritization of therapeutic options (McCoy and Madhavan, 2018; Kumar et al., 2019; Adam et al., 2020; Chang et al., 2020; Baptista et al., 2021; McCoy et al., 2021; Thiemeyer et al., 2021; Babbitt et al., 2022; Wang et al., 2022). DT can greatly aid breast cancer research by allowing for better patient stratification, improved treatment decision-making, and faster drug development. This technique has the capability to revolutionize breast cancer management by promoting personalized medicine and optimizing clinical practices.
4. Discussion
Artificial intelligence will play a crucial role in deploying DT to cancer patients. Aside from the direct application as models that predict the outcomes of different therapeutic interventions, they will also map patient data to parameters in mechanistic models of disease progression. However, applying AI to supplement the real-world component of the DT is equally essential. The progression and outcomes of the computational models of divergent intervention strategies can be validated by directly comparing patient trajectories and experiences. Also, the subset of patients that most closely match the real-world component of the DT can be used to evaluate the robustness of model predictions by providing additional inputs. For any given model prediction, convergence of its output for the subset of real-world counterparts can boost confidence in the predicted outcomes.
The utilization of DT has the potential to greatly benefit patients and their caregivers by providing tailored information, tools, and resources. Patients can leverage their DT to gain insights into their health status, monitor their progress, and make informed decisions about their care. Patients can engage in meaningful conversations about their health and treatment options by sharing the DT data with their healthcare team. This collaborative approach prioritizes the patient’s voice, ensuring their preferences and goals are considered in the decision-making process. Moreover, the data collected from DT can help identify population-level trends, patterns, and insights, which are instrumental in supporting public health initiatives, advancing research studies, and informing policy-making.
Artificial Intelligence needs data to learn from, and development in all its applications to DT is dependent on widespread data collection. Unfortunately, collecting the data required for DT is invasive and disruptive. Educating patient populations about the importance and potential of collecting is critical to realizing the clinical application of DT.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found at: https://portal.gdc.cancer.gov/projects/TCGA-BRCA TCGA BRCA dataset.
Author contributions
H-CC: Investigation, Methodology, Software, Writing – original draft, Writing – review & editing. AG: Data curation, Investigation, Visualization, Writing – review & editing. SK: Data curation, Investigation, Visualization, Writing – review & editing. DW: Writing – review & editing. MS: Data curation, Investigation, Methodology, Project administration, Supervision, Writing – original draft, Writing – review & editing. MM: Conceptualization, Data curation, Investigation, Methodology, Project administration, Software, Supervision, Visualization, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
We gratefully acknowledge the technical support Susan G. Komen and UT Southwestern provided during the 2023 Susan G. Komen Hackathon, and to PIVOT (Patient and Investigator Voices Organizing Together) at The University of Kansas Cancer Center for their valuable insight into patient advocacy.
Conflict of interest
The authors declare that the research was conducted without any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Adam, G., Rampášek, L., Safikhani, Z., Smirnov, P., Haibe-Kains, B., and Goldenberg, A. (2020). Machine learning approaches to drug response prediction: challenges and recent progress. NPJ Precis. Oncol. 4:19. doi: 10.1038/s41698-020-0122-1
Babbitt, G. A., Lynch, M. L., McCoy, M., Fokoue, E. P., and Hudson, A. O. (2022). Function and evolution of B-Raf loop dynamics relevant to cancer recurrence under drug inhibition. J. Biomol. Struct. Dyn. 40, 468–483. doi: 10.1080/07391102.2020.1815578
Baptista, D., Ferreira, P. G., and Rocha, M. (2021). Deep learning for drug response prediction in cancer. Brief. Bioinform. 22, 360–379. doi: 10.1093/bib/bbz171
Chang, H.-C., Chu, C.-P., Lin, S.-J., and Hsiao, C. K. (2020). Network hub-node prioritization of gene regulation with intra-network association. BMC Bioinformat. 21:101. doi: 10.1186/s12859-020-3444-7
Goldman, M. J., Craft, B., Hastie, M., Repečka, K., McDade, F., Kamath, A., et al. (2020). Visualizing and interpreting cancer genomics data via the Xena platform. Nat. Biotechnol. 38, 675–678. doi: 10.1038/s41587-020-0546-8
Haleem, A., Javaid, M., Singh, R. P., and Suman, R. (2023). Exploring the revolution in healthcare systems through the applications of digital twin technology. Biomed. Technol. 4, 28–38. doi: 10.1016/j.bmt.2023.02.001
Hernandez-Boussard, T., Macklin, P., Greenspan, E. J., Gryshuk, A. L., Stahlberg, E., Syeda-Mahmood, T., et al. (2021). Digital twins for predictive oncology will be a paradigm shift for precision cancer care. Nat. Med. 27, 2065–2066. doi: 10.1038/s41591-021-01558-5
Ioannidis, N. M., Rothstein, J. H., Pejaver, V., Middha, S., McDonnell, S. K., Baheti, S., et al. (2016). REVEL: an ensemble method for predicting the pathogenicity of rare missense variants. Am. J. Hum. Genet. 99, 877–885. doi: 10.1016/j.ajhg.2016.08.016
Kumar, D., Lisok, A., Dahmane, E., McCoy, M., Shelake, S., Chatterjee, S., et al. (2019). Peptide-based PET quantifies target engagement of PD-L1 therapeutics. J. Clin. Invest. 129, 616–630. doi: 10.1172/JCI122216
Liu, J., Lichtenberg, T., Hoadley, K. A., Poisson, L. M., Lazar, A. J., Cherniack, A. D., et al. (2018). An integrated TCGA pan-cancer clinical data resource to drive high-quality survival outcome analytics. Cells 173, 400–16.e11. doi: 10.1016/j.cell.2018.02.052
Madhavan, S., Beckman, R. A., McCoy, M. D., Pishvaian, M. J., Brody, J. R., and Macklin, P. (2021). Envisioning the future of precision oncology trials. Nat. Cancer 2, 9–11. doi: 10.1038/s43018-020-00163-8
McClellan, J., and King, M. C. (2010). Genetic heterogeneity in human disease. Cells 141, 210–217. doi: 10.1016/j.cell.2010.03.032
McCoy, M. D., Hamre, J., Klimov, D. K., and Jafri, M. S. (2021). Predicting genetic variation severity using machine learning to interpret molecular simulations. Biophys. J. 120, 189–204. doi: 10.1016/j.bpj.2020.12.002
McCoy, M. D., and Madhavan, S. (2018). A computational approach for prioritizing selection of therapies targeting drug resistant variation in anaplastic lymphoma kinase. AMIA Jt. Summits Transl. Sci. Proc. 2017:160:167.
Stahlberg, E. A., Abdel-Rahman, M., Aguilar, B., Asadpoure, A., Beckman, R. A., Borkon, L. L., et al. (2022). Exploring approaches for predictive cancer patient digital twins: opportunities for collaboration and innovation. Front. Digit. Health 4:1007784. doi: 10.3389/fdgth.2022.1007784
Sun, T., He, X., and Li, Z. (2023). Digital twin in healthcare: recent updates and challenges. Digit. Health 9:20552076221149651. doi: 10.1177/20552076221149651
Tam, V., Patel, N., Turcotte, M., Bossé, Y., Paré, G., and Meyre, D. (2019). Benefits and limitations of genome-wide association studies. Nat. Rev. Genet. 20, 467–484. doi: 10.1038/s41576-019-0127-1
Thennavan, A., Beca, F., Xia, Y., Recio, S. G., Allison, K., Collins, L. C., et al. (2021). Molecular analysis of TCGA breast cancer histologic types. Cell Genom. 1:100067. doi: 10.1016/j.xgen.2021.100067
Thiemeyer, H., Taher, L., Schille, J. T., Packeiser, E.-M., Harder, L. K., Hewicker-Trautwein, M., et al. (2021). An RNA-Seq-based framework for characterizing canine prostate cancer and prioritizing clinically relevant biomarker candidate genes. Int. J. Mol. Sci. 22:11481. doi: 10.3390/ijms222111481
Venkatesh, K. P., Raza, M. M., and Kvedar, J. C. (2022). Health digital twins as tools for precision medicine: considerations for computation, implementation, and regulation. NPJ Digit. Med. 5:150. doi: 10.1038/s41746-022-00694-7
Wang, H.-Y., Hsieh, T.-T., Chung, C.-R., Chang, H.-C., Horng, J.-T., Lu, J.-J., et al. (2022). Efficiently predicting vancomycin resistance of Enterococcus faecium from MALDI-TOF MS spectra using a deep learning-based approach. Front. Microbiol. 13:821233. doi: 10.3389/fmicb.2022.821233
Keywords: digital twins, cancer, precision medicine, omics, data integration
Citation: Chang H-C, Gitau AM, Kothapalli S, Welch DR, Sardiu ME and McCoy MD (2023) Understanding the need for digital twins’ data in patient advocacy and forecasting oncology. Front. Artif. Intell. 6:1260361. doi: 10.3389/frai.2023.1260361
Edited by:
Qi Wang, University of South Carolina, United StatesReviewed by:
Gregory R. Hart, Institute for Disease Modeling (IDM), United StatesCopyright © 2023 Chang, Gitau, Kothapalli, Welch, Sardiu and McCoy. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Matthew D. McCoy, bWF0dGhldy5tY2NveUBnZW9yZ2V0b3duLmVkdQ==