A novel heterophilic graph diffusion convolutional network for identifying cancer driver genes.

Identifying cancer driver genes plays a curial role in the development of precision oncology and cancer therapeutics. Although a plethora of methods have been developed to tackle this problem, the complex cancer mechanisms and intricate interactions between genes still make the identification of cancer driver genes challenging. In this work, we propose a novel machine learning method of heterophilic graph diffusion convolutional networks (called HGDCs) to boost cancer-driver gene identification. Specifically, HGDC first introduces graph diffusion to generate an auxiliary network for capturing the structurally similar nodes in a biomolecular network. Then, HGDC designs an improved message aggregation and propagation scheme to adapt to the heterophilic setting of biomolecular networks, alleviating the problem of driver gene features being smoothed by its neighboring dissimilar genes. Finally, HGDC uses a layer-wise attention classifier to predict the probability of one gene being a cancer driver gene. In the comparison experiments with other existing state-of-the-art methods, our HGDC achieves outstanding performance in identifying cancer driver genes. The experimental results demonstrate that HGDC not only effectively identifies well-known driver genes on different networks but also novel candidate cancer genes. Moreover, HGDC can effectively prioritize cancer driver genes for individual patients. Particularly, HGDC can identify patient-specific additional driver genes, which work together with the well-known driver genes to cooperatively promote tumorigenesis.

Predicting Lipid-Rich Plaque Progression in Coronary Arteries Using Multimodal Imaging and Wall Shear Stress Signatures.

BACKGROUND: Plaque composition and wall shear stress (WSS) magnitude act as well-established players in coronary plaque progression. However, WSS magnitude per se does not completely capture the mechanical stimulus to which the endothelium is subjected, since endothelial cells experience changes in the WSS spatiotemporal configuration on the luminal surface. This study explores WSS profile and lipid content signatures of plaque progression to identify novel biomarkers of coronary atherosclerosis. METHODS: Thirty-seven patients with acute coronary syndrome underwent coronary computed tomography angiography, near-infrared spectroscopy intravascular ultrasound, and optical coherence tomography of at least 1 nonculprit vessel at baseline and 1-year follow-up. Baseline coronary artery geometries were reconstructed from intravascular ultrasound and coronary computed tomography angiography and combined with flow information to perform computational fluid dynamics simulations to assess the time-averaged WSS magnitude (TAWSS) and the variability in the contraction/expansion action exerted by WSS on the endothelium, quantifiable in terms of topological shear variation index (TSVI). Plaque progression was measured as intravascular ultrasound-derived percentage plaque atheroma volume change at 1-year follow-up. Plaque composition information was extracted from near-infrared spectroscopy and optical coherence tomography. RESULTS: Exposure to high TSVI and low TAWSS was associated with higher plaque progression (4.00±0.69% and 3.60±0.62%, respectively). Plaque composition acted synergistically with TSVI or TAWSS, resulting in the highest plaque progression (≥5.90%) at locations where lipid-rich plaque is exposed to high TSVI or low TAWSS. CONCLUSIONS: Luminal exposure to high TSVI, solely or combined with a lipid-rich plaque phenotype, is associated with enhanced plaque progression at 1-year follow-up. Where plaque progression occurred, low TAWSS was also observed. These findings suggest TSVI, in addition to low TAWSS, as a potential biomechanical predictor for plaque progression, showing promise for clinical translation to improve patient prognosis.

Plaque Ruptures Are Related to High Plaque Stress and Strain Conditions: Direct Verification by Using In Vivo OCT Rupture Data and FSI Models.

BACKGROUND: While it has been hypothesized that high plaque stress and strain may be related to plaque rupture, its direct verification using in vivo coronary plaque rupture data and full 3-dimensional fluid-structure interaction models is lacking in the current literature due to difficulty in obtaining in vivo plaque rupture imaging data from patients with acute coronary syndrome. This case-control study aims to use high-resolution optical coherence tomography-verified in vivo plaque rupture data and 3-dimensional fluid-structure interaction models to seek direct evidence for the high plaque stress/strain hypothesis. METHODS: In vivo coronary plaque optical coherence tomography data (5 ruptured plaques, 5 no-rupture plaques) were acquired from patients using a protocol approved by the local institutional review board with informed consent obtained. The ruptured caps were reconstructed to their prerupture morphology using neighboring plaque cap and vessel geometries. Optical coherence tomography-based 3-dimensional fluid-structure interaction models were constructed to obtain plaque stress, strain, and flow shear stress data for comparative analysis. The rank-sum test in the nonparametric test was used for statistical analysis. RESULTS: Our results showed that the average maximum cap stress and strain values of ruptured plaques were 142% (457.70 versus 189.22 kPa; P=0.0278) and 48% (0.2267 versus 0.1527 kPa; P=0.0476) higher than that for no-rupture plaques, respectively. The mean values of maximum flow shear stresses for ruptured and no-rupture plaques were 145.02 dyn/cm2 and 81.92 dyn/cm2 (P=0.1111), respectively. However, the flow shear stress difference was not statistically significant. CONCLUSIONS: This preliminary case-control study showed that the ruptured plaque group had higher mean maximum stress and strain values. Due to our small study size, larger scale studies are needed to further validate our findings.

Hemodynamics and Wall Mechanics of Vascular Graft Failure.

Blood vessels are subjected to complex biomechanical loads, primarily from pressure-driven blood flow. Abnormal loading associated with vascular grafts, arising from altered hemodynamics or wall mechanics, can cause acute and progressive vascular failure and end-organ dysfunction. Perturbations to mechanobiological stimuli experienced by vascular cells contribute to remodeling of the vascular wall via activation of mechanosensitive signaling pathways and subsequent changes in gene expression and associated turnover of cells and extracellular matrix. In this review, we outline experimental and computational tools used to quantify metrics of biomechanical loading in vascular grafts and highlight those that show potential in predicting graft failure for diverse disease contexts. We include metrics derived from both fluid and solid mechanics that drive feedback loops between mechanobiological processes and changes in the biomechanical state that govern the natural history of vascular grafts. As illustrative examples, we consider application-specific coronary artery bypass grafts, peripheral vascular grafts, and tissue-engineered vascular grafts for congenital heart surgery as each of these involves unique circulatory environments, loading magnitudes, and graft materials.

Patient-specific signaling signatures predict optimal therapeutic combinations for triple negative breast cancer.

Triple negative breast cancer (TNBC) is a heterogeneous group of tumors which lack estrogen receptor, progesterone receptor, and HER2 expression. Targeted therapies have limited success in treating TNBC, thus a strategy enabling effective targeted combinations is an unmet need. To tackle these challenges and discover individualized targeted combination therapies for TNBC, we integrated phosphoproteomic analysis of altered signaling networks with patient-specific signaling signature (PaSSS) analysis using an information-theoretic, thermodynamic-based approach. Using this method on a large number of TNBC patient-derived tumors (PDX), we were able to thoroughly characterize each PDX by computing a patient-specific set of unbalanced signaling processes and assigning a personalized therapy based on them. We discovered that each tumor has an average of two separate processes, and that, consistent with prior research, EGFR is a major core target in at least one of them in half of the tumors analyzed. However, anti-EGFR monotherapies were predicted to be ineffective, thus we developed personalized combination treatments based on PaSSS. These were predicted to induce anti-EGFR responses or to be used to develop an alternative therapy if EGFR was not present.In-vivo experimental validation of the predicted therapy showed that PaSSS predictions were more accurate than other therapies. Thus, we suggest that a detailed identification of molecular imbalances is necessary to tailor therapy for each TNBC. In summary, we propose a new strategy to design personalized therapy for TNBC using pY proteomics and PaSSS analysis. This method can be applied to different cancer types to improve response to the biomarker-based treatment.

A patient-based murine model recapitulates human STAT3 gain-of-function syndrome.

STAT3 gain-of-function (GOF) variants results in a heterogeneous clinical syndrome characterized by early onset immunodeficiency, multi-organ autoimmunity, and lymphoproliferation. While 191 documented cases with STAT3 GOF variants have been reported, the impact of individual variants on immune regulation and the broad clinical spectrum remains unclear. We developed a Stat3p.L387R mouse model, mirroring a variant identified in a family exhibiting common STAT3 GOF symptoms, and rare phenotypes including pulmonary hypertension and retinal vasculitis. In vitro experiments revealed increased STAT3 phosphorylation, nuclear migration, and DNA binding of the variant. Our Stat3p.L387R model displayed similar traits from previous Stat3GOF strains, such as splenomegaly and lymphadenopathy. Notably, Stat3p.L387R/+ mice exhibited heightened embryonic lethality compared to prior Stat3GOF/+ models and ocular abnormalities were observed. This research underscores the variant-specific pathology in Stat3p.L387R/+ mice, highlighting the ability to recapitulate human STAT3 GOF syndrome in patient-specific transgenic murine models. Additionally, such models could facilitate tailored treatment development.

Initial needle tracking with the first standalone combined infrared camera - CT system for brachytherapy-analysis of tracking accuracy and uncertainties.

PURPOSE: A prototype infrared camera - cone-beam computed tomography (CBCT) system for tracking in brachytherapy has recently been developed. We evaluated for the first time the corresponding tracking accuracy and uncertainties, and implemented a tracking-based prediction of needles on CBCT scans. METHODS: A marker tool rigidly attached to needles was 3D printed. The precision and accuracy of tool tracking was then evaluated for both static and dynamic scenarios. Euclidean distances between the tracked and CBCT-derived markers were assessed as well. To implement needle tracking, ground truth models of the tool attached to 200 mm and 160 mm needles were matched to the tracked positions in order to project the needles into CBCT scans. Deviations between projected and actual needle tips were measured. Finally, we put our results into perspective with simulations of the system's tracking uncertainties. RESULTS: For the stationary scenario and dynamic movements, we achieved tool-tracking precision and accuracy of 0.04 ± 0.06 mm and 0.16 ± 0.18 mm, respectively. The tracked marker positions differed by 0.52 ± 0.18 mm from the positions determined via CBCT. In addition, the predicted needle tips in air deviated from the actual tip positions by only 1.62 ± 0.68 mm (200 mm needle) and 1.49 ± 0.62 mm (160 mm needle). The simulated tracking uncertainties resulted in tip variations of 1.58 ± 0.91 mm and 1.31 ± 0.69 mm for the 200 mm and 160 mm needles, respectively. CONCLUSION: With the innovative system it was possible to achieve a high tracking and prediction accuracy of marker tool and needles. The system shows high potential for applicator tracking in brachytherapy.

Pre-operative lung ablation prediction using deep learning.

OBJECTIVE: Microwave lung ablation (MWA) is a minimally invasive and inexpensive alternative cancer treatment for patients who are not candidates for surgery/radiotherapy. However, a major challenge for MWA is its relatively high tumor recurrence rates, due to incomplete treatment as a result of inaccurate planning. We introduce a patient-specific, deep-learning model to accurately predict post-treatment ablation zones to aid planning and enable effective treatments. MATERIALS AND METHODS: Our IRB-approved retrospective study consisted of ablations with a single applicator/burn/vendor between 01/2015 and 01/2019. The input data included pre-procedure computerized tomography (CT), ablation power/time, and applicator position. The ground truth ablation zone was segmented from follow-up CT post-treatment. Novel deformable image registration optimized for ablation scans and an applicator-centric co-ordinate system for data analysis were applied. Our prediction model was based on the U-net architecture. The registrations were evaluated using target registration error (TRE) and predictions using Bland-Altman plots, Dice co-efficient, precision, and recall, compared against the applicator vendor's estimates. RESULTS: The data included 113 unique ablations from 72 patients (median age 57, interquartile range (IQR) (49-67); 41 women). We obtained a TRE ≤ 2 mm on 52 ablations. Our prediction had no bias from ground truth ablation volumes (p = 0.169) unlike the vendor's estimate (p < 0.001) and had smaller limits of agreement (p < 0.001). An 11% improvement was achieved in the Dice score. The ability to account for patient-specific in-vivo anatomical effects due to vessels, chest wall, heart, lung boundaries, and fissures was shown. CONCLUSIONS: We demonstrated a patient-specific deep-learning model to predict the ablation treatment effect prior to the procedure, with the potential for improved planning, achieving complete treatments, and reduce tumor recurrence. CLINICAL RELEVANCE STATEMENT: Our method addresses the current lack of reliable tools to estimate ablation extents, required for ensuring successful ablation treatments. The potential clinical implications include improved treatment planning, ensuring complete treatments, and reducing tumor recurrence.

Development of a patient-specific model of the human coronary system for percutaneous transluminal coronary angioplasty balloon catheter training and testing.

BACKGROUND: To treat stenosed coronary arteries, percutaneous transluminal coronary angioplasty (PTCA) balloon catheters must combine pushability, trackability, crossability, and rewrap behavior. The existing anatomic track model (ASTM F2394) for catheter testing lacks 3D morphology, vessel tortuosity, and compliance, making evaluating performance characteristics difficult. This study aimed to develop a three-dimensional patient-specific phantom (3DPSP) for device testing and safe training for interventional cardiologists. METHODS: A range of silicone materials with different shore hardnesses (00-30-45 A) and wall thicknesses (0.5 mm, 1 mm, 2 mm) were tested to determine compliance for creating coronary vessel phantoms. Compliance was assessed using optical coherence tomography (OCT) and compared to values in the literature. Stenosis was induced using multilayer casting and brushing methods, with gypsum added for calcification. The radial tensile properties of the samples were investigated, and the relationship between Young's modulus and compliance was determined. Various methods have been introduced to approximate the friction between silicone and real coronary vessel walls. Computerized tomography (CT) scans were used to obtain patient-specific anatomy from the femoral artery to the coronary arteries. Artery lumens were segmented from the CT scans to create dissolvable 3D-printed core models. RESULTS: A 15A shore hardness silicone yielded an experimental compliance of 12.3-22.4 m m 2 mmHg · 10 3 for stenosed tubes and 14.7-57.9 m m 2 mmHg · 10 3 for uniform tubes, aligning closely with the literature data (6.28-40.88 m m 2 mmHg · 10 3 ). The Young's modulus ranged from 43.2 to 75.5 kPa and 56.6-67.9 kPa for the uniform and calcified materials, respectively. The dependency of the compliance on the wall thickness, Young's modulus, and inner diameter could be shown. Introducing a lubricant reduced the silicone friction coefficient from 0.52 to 0.13. The 3DPSP was successfully fabricated, and comparative analyses were conducted among eight commercially available catheters. CONCLUSION: This study presents a novel method for crafting 3DPSPs with realistic mechanical and frictional properties. The proposed approach enables the creation of comprehensive and anatomically precise setups spanning the right femoral artery to the coronary arteries, highlighting the importance of such realistic environments for advancing medical device development and fostering safe training conditions.

An In Silico Model for Predicting the Efficacy of Edge-to-Edge Repair for Mitral Regurgitation.

In recent years, transcatheter edge-to-edge repair (TEER) has been widely adopted as an effective treatment for mitral regurgitation (MR). The aim of this study is to develop a personalized in silico model to predict the effect of edge-to-edge repair in advance to the procedure for each individual patient. For this purpose, we propose a combination of a valve deformation model for computing the mitral valve (MV) orifice area (MVOA) and a lumped parameter model for the hemodynamics, specifically mitral regurgitation volume (RVol). Although we cannot obtain detailed information on the three-dimensional flow field near the mitral valve, we can rapidly simulate the important medical parameters for the clinical decision support. In the present method, we construct the patient-specific pre-operative models by using the parameter optimization and then simulate the postoperative state by applying the additional clipping condition. The computed preclip MVOAs show good agreement with the clinical measurements, and the correlation coefficient takes 0.998. In addition, the MR grade in terms of RVol also has good correlation with the grade by ground truth MVOA. Finally, we try to investigate the applicability for the predicting the postclip state. The simulated valve shapes clearly show the well-known double orifice and the improvement of the MVOA, compared with the preclip state. Similarly, we confirmed the improved reverse flow and MR grade in terms of RVol. A total computational time is approximately 8 h by using general-purpose PC. These results obviously indicate that the present in silico model has good capability for the assessment of edge-to-edge repair.