Comprehensive cardiac magnetic resonance imaging in patients with idiopathic premature ventricular contractions: role of feature-tracking strain and T1/T2 mapping in detecting subclinical myocardial dysfunction
Original Article

Comprehensive cardiac magnetic resonance imaging in patients with idiopathic premature ventricular contractions: role of feature-tracking strain and T1/T2 mapping in detecting subclinical myocardial dysfunction

Jia Yang# ORCID logo, Yi Zhang# ORCID logo, Qinqin Yu ORCID logo, Xinxing Hu ORCID logo, Qian Luo ORCID logo, Lei Zhang ORCID logo

Department of Radiology, Shanghai General Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China

Contributions: (I) Conception and design: L Zhang; (II) Administrative support: L Zhang, Q Luo; (III) Provision of study materials or patients: J Yang, Y Zhang, X Hu, Q Luo; (IV) Collection and assembly of data: J Yang, Y Zhang, Q Yu, X Hu; (V) Data analysis and interpretation: J Yang, Y Zhang, L Zhang; (VI) Manuscript writing: All authors; (VII) Final approval of manuscript: All authors.

#These authors contributed equally to this work.

Correspondence to: Lei Zhang, MD, PhD; Qian Luo, MD. Department of Radiology, Shanghai General Hospital, Shanghai Jiao Tong University School of Medicine, 650 Xinsongjiang Road, Songjiang District, Shanghai 201620, China. Email: lei.zhang2@shgh.cn; sisiandlb@outlook.com.

Background: Although standard cardiac magnetic resonance (CMR) findings are typically normal in patients with idiopathic premature ventricular contractions (PVCs), it remains unclear whether a comprehensive protocol incorporating left ventricular (LV) strain and tissue characterization can detect subclinical myocardial abnormalities. This study aimed to determine if advanced CMR parameters can reveal a burden-dependent myocardial phenotype in patients with idiopathic PVCs.

Methods: In this retrospective study, 72 patients with idiopathic PVCs and 25 controls underwent comprehensive 3T CMR with feature-tracking (FT) strain analysis and native T1/T2 mapping. Patients were stratified by 24-hour PVC burden into three groups: <500 PVC/24 h (n=24), 500–10,000 PVC/24 h (n=23), and >10,000 PVC/24 h (n=25).

Results: Compared to controls, patients with >10,000 PVC/24 h (high burden) exhibited significantly impaired myocardial strain after multiple-comparison correction, including reductions in two-dimensional (2D) global radial (27.4%±7.0% vs. 33.9%±6.2%; P=0.005), 2D global circumferential (−16.7%±2.8% vs. −19.2%±2.1%; P=0.001), three-dimensional (3D) global circumferential (−17.5%±3.7% vs. −20.6%±2.0%; P=0.003), 2D apical radial (28.1%±11.0% vs. 40.2%±12.2%; P=0.002), and 2D apical circumferential strain (−16.7%±4.4% vs. −20.8%±3.8%; P=0.001). Notably, the most pronounced strain impairment was observed in the apical segments. These strain parameters were also significantly impaired in the high-burden compared to the low-burden (<500 PVC/24 h) group (all P values <0.05). Moreover, the high-burden demonstrated elevated global native T1 (1,344.5±40.0 vs. 1,289.4±16.6 ms in controls; P<0.001) and T2 (54.4±3.7 vs. 49.4±3.1 ms in controls; P<0.001) values, with the most pronounced alterations localized to the basal segments (T1 and T2; P<0.001 for comparison with controls). A significant dose-response relationship was confirmed for both strain and mapping parameters across the burden groups.

Conclusions: As a hypothesis-generating study, our work, based on comprehensive CMR assessment, identified a novel, burden-dependent myocardial phenotype in patients with idiopathic PVC. The discovery of this phenotype, characterized by distinct patterns of systolic dysfunction and interstitial alterations, should be examined in terms of its clinical and prognostic significance in subsequent longitudinal studies.

Keywords: Cardiac magnetic resonance (CMR); premature ventricular contractions (PVCs); feature-tracking (FT); T1 mapping; T2 mapping


Submitted May 22, 2025. Accepted for publication Nov 21, 2025. Published online Jan 09, 2026.

doi: 10.21037/qims-2025-1190


Introduction

Premature ventricular contractions (PVCs) are common both in individuals with and those without structural heart disease (1-4). Although frequent occurrence of PVCs is associated with an elevated risk of sudden and total cardiac death even in the absence of demonstrable underlying disease (5-9), whether an underlying structural substrate is present in these patients—and its characteristics—remains unclear (10-13).

Evaluation of patients with PVCs typically begins with transthoracic echocardiography (TTE), as recommended by major societal guidelines (14-19). However, conventional TTE has limited sensitivity for detecting subtle arrhythmogenic substrates, leaving a controversy as whether advanced imaging is needed in select cases (16). Cardiac magnetic resonance (CMR), particularly late gadolinium enhancement (LGE), has proven invaluable for identifying focal scar and stratifying arrhythmic risk (12,20). However, a significant proportion (54–85%) of patients experiencing PVCs show no LGE or obvious structural abnormalities on routine imaging (21-24).

Comprehensive CMR extends beyond LGE, incorporating techniques for the quantification of diffuse tissue pathology through native T1 and T2 mapping, which can detect interstitial fibrosis and edema (25). Although these techniques show promise in the characterization of the arrhythmic substrate in known structural heart disease (26-28), their role for idiopathic PVCs remains less defined. Similarly, CMR feature-tracking (FT) strain analysis offers a sensitive measure of systolic function but remains underinvestigated in the context of ventricular arrhythmias (VAs) (13,29-32).

Given the potential of comprehensive CMR to detect early tissue alterations, we conducted a systematic evaluation in patients with idiopathic PVC and normal CMR findings. This exploratory study aimed to determine whether advanced CMR parameters can reveal subtle myocardial abnormalities and to investigate their potential association with PVC burden. We present this article in accordance with the STROBE reporting checklist (available at https://qims.amegroups.com/article/view/10.21037/qims-2025-1190/rc).


Methods

Ethical considerations

This study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments and was approved by the Ethics Review Committee of Shanghai General Hospital (No. 2023364). The requirement for informed consent was waived due to the retrospective nature of the study.

Participant population and study design

In this single-center, retrospective study conducted at Shanghai General Hospital, eligible patients (October 2017–December 2024) were identified through a systematic query of the institutional electronic medical records. The search encompassed a broad spectrum of patients evaluated at our tertiary center, which included both referred individuals and those presenting directly for evaluation. The inclusion criteria were (I) documented PVCs on 24-hour Holter monitoring and (II) CMR examination performed within 7 days of Holter monitoring. Meanwhile, the exclusion criteria were as follows: (I) any form of structural heart disease (including but not limited to cardiomyopathies, significant valvular heart disease, congenital heart disease, or any other cardiac structural abnormality), ischemic cardiomyopathy, and mild-to-severe coronary atherosclerosis (any coronary artery stenosis ≥25%), or left ventricular ejection fraction (LVEF) <50%; (II) myocarditis or pericardial disease; (III) any systemic diseases that could affect the heart, such as hypertension, diabetes mellitus, connective tissue disease, and hyperthyroidism; (IV) electrocardiogram findings of multifocal PVC or atrioventricular blocks; (V) delayed myocardial contrast enhancement in CMR imaging; (VI) structural or functional abnormalities detected by cardiac ultrasound; (VII) previous administration of radiofrequency or alcohol ablation; (VIII) poor or incomplete images; and (IX) a family history of cardiomyopathy or premature sudden death or a history of prior cardiac arrest or persistent ventricular tachycardia. Based on these classification criteria, 72 adult patients were included in the final analysis. The patients were divided into three groups based on the number of PVCs within a 24-hour period (<500 PVCs, 500–10,000 PVCs, and >10,000 PVCs). The control group comprised 25 individuals who underwent CMR (e.g., for evaluation of atypical chest pain or palpitations) and were confirmed to have completely normal cardiac structure, function, absence of LGE, and normal T1/T2 mapping values. Moreover, none of the controls demonstrated evidence of PVCs on 24-hour Holter monitoring or met any of the study’s exclusion criteria.

CMR studies

A standardized CMR protocol was maintained throughout the study period to ensure technical consistency. All examinations were performed on a 3.0-T Ingenia scanner (Philips Healthcare, Amsterdam, the Netherlands) with a dStream torso coil and vector-electrocardiographic gating. The protocol included (I) cine imaging for ventricular function and strain analysis; (II) native T1 and T2 mapping for tissue characterization; and (III) LGE for focal fibrosis assessment. Cine images were acquired with a balanced steady-state free precession (SSFP) sequence during repeated breath-holds at end expiration, which covered the long-axis (two-, three-, and four-chamber) and short-axis (stacked from base to apex) planes. The typical sequence parameters were as follows: slice thickness =8 mm, repetition time/echo time =2.8–3.2/1.4–1.5 ms, flip angle =45°; matrix size =[160–176]×[138–192], field of view =300×300 to 350×350 mm2, and temporal resolution =15–25 ms. For LGE imaging, acquisitions commenced 10–15 minutes after intravenous administration of 0.2 mmol/kg of gadopentetic acid with a phase-sensitive inversion recovery sequence. Native T1 mapping was performed precontrast with a modified Look-Locker inversion recovery sequence in short-axis views. T2 mapping was acquired with a T2-prepared SSFP sequence with T2-preparation times of 0, 30, and 55 ms.

All image analyses were performed via CVI42 software version 5.17.1 (Circle Cardiovascular Imaging, Calgary, Canada) by an experienced observer blinded to clinical data. Left ventricular (LV) volumes, ejection fraction, and mass were quantified from short-axis cine stacks. For FT strain analysis, endocardial and epicardial contours were semiautomatically traced on cine images, with manual adjustments being applied as needed. Global and segmental peak systolic strain parameters (longitudinal, circumferential, and radial) were derived in two-dimensional (2D) and three-dimensional (3D). For tissue characterization, T1 and T2 maps were generated, and values were automatically sampled from the 16 American Heart Association segments (excluding the apex), with global and regional (basal, mid, and apical) averages calculated for analysis.

Inter- and intraobserver reproducibility

To assess interobserver reproducibility, 20 randomly selected participants (10 patients and 10 controls) were reanalyzed by a second blinded observer. To assess intraobserver reproducibility, the primary observer reanalyzed the same 20 participants after a washout period of at least 4 weeks.

Statistical analysis

Statistical analyses were performed with SPSS version 25 (IBM Corp., Armonk, NY, USA). Missing data for specific parameters were handled via exclusion of cases on an analysis-by-analysis basis; each test was performed with all available data for the relevant variables. Normally distributed continuous data are presented as the mean ± standard deviation, nonnormally distributed data as the median with interquartile range, and categorical data as frequencies (percentages). Group comparisons were performed via one-way analysis of variance or the Kruskal-Wallis test, as appropriate, with a two-tiered approach to control for multiple comparisons: first, false-discovery rate (FDR) correction was applied to P values from omnibus tests within prespecified parameter families (strain parameters and T1/T2 mapping parameters), with an FDR-corrected P value (Q value) <0.05 being considered significant. Second, for parameters surpassing this threshold, post hoc pairwise comparisons were conducted via the Tukey test (parametric) or Dunn test with Bonferroni adjustment (nonparametric). Associations between PVC burden and CMR parameters were assessed with Spearman rank correlation and partial correlation analyses (with adjustments for age, sex, and LV volumes); meanwhile, burden-dependent trends were evaluated with the Jonckheere-Terpstra test, with a two-sided P value (or Q value) <0.05 being considered statistically significant.


Results

Patient demographics

Among the 288 patients who initially met the inclusion criteria, 216 were excluded base on the predefined criteria, resulting in a final cohort of 72 patients with idiopathic PVC and 25 controls (Figure 1). The baseline characteristics are summarized in Table 1. Groups were well-balanced in terms of age, sex, and clinical comorbidities. The only exception was LV end-systolic volume (ESV), which was significantly higher in both the medium- (10,000–500 PVC/24 h) and high-burden (>10,000 PVC/24 h) groups than in the control group (P<0.05). Consistent with this, Spearman correlation analysis indicated weak positive correlations between PVC burden and both LV-ESV (r=0.258) and LV end-diastolic volume (EDV; r=0.261) (Figure S1).

Figure 1 Patient flowchart. The flowchart shows the enrollment of patients and healthy controls. CMR, cardiac magnetic resonance; PVC, premature ventricular contraction.

Table 1

Baseline characteristics of all participants

Characteristic Controls (n=25) <500 PVC/24 h (n=24) 500–10,000 PVC/24 h (n=23) >10,000 PVC/24 h (n=25) P value (ANOVA)
Clinical variable
   Age (years) 40±13 46±13 41±14 47±14 0.228
   Sex, female 13 (52.0) 11 (45.8) 12 (52.2) 15 (60.0) 0.802
   BMI (kg/m2) 22.5 (21.3, 24.6) 23.2 (21.3, 26.1) 22.8 (20.1, 25.9) 23.8 (22.5, 25.2) 0.605
   BSA (m2) 1.72±0.19 1.79±0.20 1.75±0.23 1.72±0.19 0.614
   Heart rate (bpm) 70±11 67±10 73±13 72±11 0.207
   SBP (mmHg) 123±9 124±9 120±15 120±13 0.576
   DBP (mmHg) 77±6 77±8 75±9 71±10 0.138
   CAD 2 (8.3) 1 (4.3) 2 (8.0) >0.99
   PVC burden 0.01 (0.00, 0.11) 3.13 (1.30, 6.30) 20.10 (14.80, 26.53) <0.001
   Antiarrhythmic drug 0 0 4 (16.0) 0.031
   EP examination 0 2 (13.3) 21 (84.0) <0.001
Cardiovascular MRI variable
   LVEF (%) 63±4 62±5 60±4 60±7 0.081
   LVCO (L/min) 4.8±1.3 5.1±1.1 5.5±1.4 5.1±0.8 0.431
   LV mass/BSA (g/m2) 38.5 (33.0, 43.7) 41.6 (36.7, 48.6) 41.9 (38.6, 50.8) 37.9 (32.6, 42.3) 0.056
   LVEDV/BSA (mL/m2) 66.1±9.6 68.6±12.5 73.6±12.3 72.0±11.9 0.082
   LVESV/BSA (mL/m2) 24.4±3.6 26.0±6.2 29.1±5.5 29.5±8.5 0.003
   LVSV/BSA (mL/m2) 41.6±7.3 42.6±8.0 44.4±8.4 42.7±6.6 0.641

Data are presented as mean ± standard deviation, absolute number (percentage), or median (interquartile range). , coronary artery stenosis <25%; , P<0.05 vs. the control group. “/BSA” indicates variables indexed for body surface area. ANOVA, analysis of variance; BMI, body mass index; BSA, body surface area; CAD, coronary artery disease; DBP, diastolic blood pressure; EP, electrophysiology; LV, left ventricular; LVCO, left ventricular cardiac output; LVEDV, left ventricular end-diastolic volume; LVEF, left ventricular ejection fraction; LVESV, left ventricular end-systolic volume; LVSV, left ventricular stroke volume index; MRI, magnetic resonance imaging; PVC, premature ventricular contraction; SBP, systolic blood pressure.

Notably, the four patients in the high-burden group who were on antiarrhythmic medication at the time of assessment still exhibited a high PVC burden, indicating ineffective control. Following CMR, 2 patients in the medium-burden group and 21 in the high-burden group underwent a cardiac electrophysiological examination. Analysis of available Holter records (n=54) indicated that PVCs predominantly presented with a left bundle branch block morphology (43/54, 79.6%), consistent with a right ventricular origin. In the subgroup with electrophysiological data (n=23), the right ventricular outflow tract (RVOT) was the most common site of origin (16/23, 69.6%).

Global and segmental strain analysis in patients with PVC and healthy controls

As shown in Figure 2 and Table S1, comprehensive strain analysis revealed a distinctive pattern of myocardial dysfunction in patients with a high PVC burden (>10,000 PVC/24 h). After FDR correction, significant impairments were first identified at the global level, with this group exhibiting reduced 2D global radial strain (GRS), 2D global circumferential strain (GCS), and 3D-GCS compared to both the control group (P<0.01) and the low-burden (<500 PVC/24 h) group (P<0.05). Segmental analysis further localized the most prominent alterations to the apical segments, with 2D apical radial strain (ARS) and 2D apical circumferential strain (ACS) were being reduced (P<0.01 and P<0.001, respectively), highlighting the apex as a key region of myocardial involvement.

Figure 2 Comparison of 2D left ventricular strain parameters. (A) GRS; (B) BRS; (C) MRS; (D) ARS; (E) GCS; (F) BCS; (G) MCS; (H) ACS; (I) GLS; (J) BLS; (K) MLS; (L) ALS. Data are shown for healthy controls and groups stratified by PVC burden (<500, 500–10,000, >10,000 PVC/24 h). *, P<0.05; **, P<0.01; ***, P<0.001 (Jonckheere-Terpstra test). ACS, apical circumferential strain; ALS, apical longitudinal strain; ARS, apical radial strain; BCS, basal circumferential strain; BLS, basal longitudinal strain; BRS, basal radial strain; GCS, global circumferential strain; GLS, global longitudinal strain; GRS, global radial strain; MCS, mid-circumferential strain; MLS, mid-longitudinal strain; MRS, mid-radial strain; PVC, premature ventricular contraction.

The biological relevance of these findings was strongly supported by a consistent, burden-dependent decrease across the vast majority of strain parameters, as evidenced by significant Jonckheere-Terpstra tests (P values ranging from <0.001 to <0.05). This trend was also present in several longitudinal strain parameters [e.g., 2D global longitudinal strain (GLS) and 3D-GLS], which while not meeting the significance threshold for specific intergroup comparisons after FDR correction, demonstrated a clear progression of worsening function with increasing PVC burden.

Spearman correlation and partial correlation analyses investigated the association between global and partial correlation parameters and PVC burden. The results of the simple correlation analysis indicated a moderate correlation (correlation coefficient >0.3) between 2D-GRS, 2D-ARS, 2D-GCS, 2D-ACS, 2D-GLS, 3D-GCS, 3D basal circumferential strain (BCS), 3D mid-circumferential strain (MCS), and PVC burden. This correlation remained significant after age and sex were controlled for. When the EDV and ESV related to PVC burden were excluded, the correlation coefficients were slightly reduced; however, 2D-ARS, 2D-ACS, 3D-GCS, and 3D-BCS still exhibited a moderate significant correlation with PVC burden (Figure 3 and Figure S1).

Figure 3 Simple and partial Spearman correlation between strain parameters and PVC burden. (A) Simple correlation coefficients. (B,C) Partial correlation coefficients. *, P<0.05; **, P<0.01; ***, P<0.001. 2D, two-dimensional; ACS, apical circumferential strain; ALS, apical longitudinal strain; ARS, apical radial strain; BCS, basal circumferential strain; BLS, basal longitudinal strain; BRS, basal radial strain; EDV, end-diastolic volume; ESV, end-systolic volume; GCS, global circumferential strain; GLS, global longitudinal strain; GRS, global radial strain; MCS, mid-circumferential strain; MLS, mid-longitudinal strain; MRS, mid-radial strain; PVC, premature ventricular contraction.

Global and segmental LV myocardial native T1 and T2 values

Native T1 values were available for 48 patients with PVC, while T2 values were available for 55 patients with PVC. The mean global and segmental LV myocardial T1 and T2 values are presented in Tables 2,3 and Figure 4. The bull’s-eye plots in Figure 5 illustrate the average septal and free-wall segments of the LV myocardial T1 and T2 values, while Figure 6 presents a typical example.

Table 2

Global and regional LV myocardial T1 values for the healthy controls and PVC cohorts

Native T1 (ms) PVC/24 h Controls (n=20) P value
<500 (n=18) 500–10,000 (n=15) >10,000 (n=15) ANOVA >10,000 vs. controls >10,000 vs. <500 >10,000 vs. 500–10,000 JT test
Apical 1,301.7±42.8 1,312.3±53.9 1,339.9±36.5†‡ 1,300.0±28.3 0.020 0.008 0.044 0.374 0.015
   Apical-FW 1,300.5±44.6 1,309.1±55.8 1,338.2±36.1 1,297.6±31.0 0.036 0.031 0.061 0.243 0.015
   Apical-IVS 1,313.5
(1,264.7, 1,359.0)
1,309.1
(1,273.4, 1,352.6)
1,338.7
(1,317.2, 1,381.7)
1,310.2
(1,274.2, 1,323.0)
0.129 0.044
Mid 1,279.0±46.5 1,298.1±52.7 1,319.6±47.8 1,273.5±22.2 0.022 0.013 0.088 0.648 0.007
   Mid-FW 1,281.6±54.3 1,294.2±64.9 1,346.6±66.5 1,272.0±25.5 0.035 0.026 0.112 0.705 0.011
   Mid-IVS 1,279.6
(1,249.1, 1,318.2)
1,284.5
(1,257.5, 1,338.8)
1,322.0
(1,290.6, 1,349.6)
1,272.5
(1,261.6, 1,291.5)
0.035 0.023 0.099 0.500 0.007
Basal 1,283.8±49.2 1,305.6±65.7 1,349.3±55.7†‡ 1,277.2±21.7 0.003 0.001 0.007 0.198 <0.001
   Basal-FW 1,297.9
(1,218.5, 1,332.3)
1,287.4
(1,226.6, 1,338.5)
1,327.4
(1,290.6, 1,428.7)
1,277.2
(1,251.2, 1,295.5)
0.013 0.003 0.052 0.248 0.001
   Basal-IVS 1,288.4±49.8 1,325.8±73.6 1,358.3±57.2†‡ 1,287.8±36.5 0.002 0.002 0.005 0.541 0.001
Global 1,300.4±41.9 1,314.4±44.8 1,344.5±40.0†‡ 1,289.4±16.6 0.002 <0.001 0.021 0.234 <0.001

Data are presented as the mean ± standard deviation or median (interquartile range). , P<0.05 vs. the control group; , P<0.05 vs. the <500 PVC/24 group. ANOVA, analysis of variance; FW, free wall; IVS, interventricular septal; LV, left ventricular; JT, Jonckheere-Terpstra; PVC, premature ventricular contraction.

Table 3

Global and regional LV myocardial T2 values for the healthy controls and PVC cohorts

T2 (ms) PVC/24 h Controls (n=20) P value
<500 (n=20) 500–10,000 (n=16) >10,000 (n=19) ANOVA >10,000 vs. controls >10,000 vs. <500 >10,000 vs. 500–10,000 JT test
Apical 51.1±4.3 52.8±5.5 53.7±3.1 51.3±4.2 0.210 0.298 0.232 0.917 0.026
   Apical-FW 51.0±4.2 52.4±5.3 53.2±3.1 50.8±3.8 0.229 0.274 0.358 0.940 0.045
   Apical-IVS 51.1 (48.2, 53.4) 52.9 (49.9, 55.5) 53.0 (51.2, 55.1) 51.0 (48.3, 54.0) 0.215 0.051
Mid 49.7±3.3 50.5±4.0 52.6±4.4 48.0±3.9 0.013 0.002 0.074 0.243 0.001
   Mid-FW 49.3 (47.7, 52.6) 49.7 (47.9, 52.0) 50.7 (49.6, 55.3) 47.8 (45.3, 50.8) 0.018 0.003 0.323 0.895 0.001
   Mid-IVS 49.1±3.5 49.2±3.2 51.6±4.2 48.2±3.0 0.038 0.016 0.125 0.180 0.004
Basal 51.5±3.4 52.1±4.7 54.5±4.6 48.8±4.4 0.005 <0.001 0.126 0.345 <0.001
   Basal-FW 52.2±5.3 52.6±6.1 55.0±6.7 49.0±5.5 0.044 0.012 0.446 0.633 0.004
   Basal-IVS 50.6 (46.0, 53.3) 49.9 (46.4, 53.3) 53.8 (49.0, 55.5) 48.0 (46.2, 50.3) 0.047 0.020 0.409 0.924 0.004
Global 50.8±3.2 51.6±3.1 54.4±3.7†‡ 49.4±3.1 0.001 <0.001 0.007 0.068 <0.001

Data are presented as mean ± standard deviation or median (interquartile range). , P<0.05 vs. the control group; , P<0.05 vs. the <500 PVC/24 h group. ANOVA, analysis of variance; FW, free wall; IVS, interventricular septal; LV, left ventricular; JT, Jonckheere-Terpstra; PVC, premature ventricular contraction.

Figure 4 Global and segmental native T1 and T2 values of LV. (A) Line charts of native T1 values. (B) Line charts of T2 values. (C) Bar charts of native T1 values. (D) Bar charts of T2 values. Lines represent mean values and error bars represent standard deviation. *, P<0.05; **, P<0.01; ***, P<0.001. LV, left ventricular.
Figure 5 Bull’s-eye plots displaying the mean regional LV myocardial T1 and T2 values for the control and PVC cohorts. The native T1 values (ms) are shown in red, while the native T2 values are in blue. The plots include the septal and free-wall regions of the basal (outer ring), mid-cavity (middle ring), and apical (central ring) planes. LV, left ventricular; PVC, premature ventricular contraction.
Figure 6 CMR mapping and measurement. Top row: a healthy control (a 61-year-old female patient). Bottom row: a patient with PVCs (a 66-year-old female patient; burden 25%; 34,510/24 h). (A,E) T1 mapping. (B,F) Bull’s-eye plots of T1 mapping. (C,G) T2 mapping. (D,H) Bull’s-eye plots of T2 mapping. CMR, cardiac magnetic resonance; PVC, premature ventricular contraction.

Patients with >10,000 PVC/24 h exhibited a statistically significant increase in LV global T1 values (P<0.001). Furthermore, segmental analysis of LV T1 values indicated a statistically significant increase across all segments, with the exception of the interventricular apical segments. Notably, within the interventricular basal segments, the T1 values in the high-burden group showed a significant increase compared to both the control and low-burden groups (both P values <0.01). Similarly, LV global myocardial T2 values exhibited a significant increase, with robust statistical significance (P<0.001). Segmental analysis of T2 values indicated a statistically significant increase in the basal and mid-segments, covering both septal and free wall segments (P=0.002–0.020).

The Jonckheere-Terpstra test indicated a significant upward tendency in both T1 and T2 values with increased PVC frequency (P≤0.011 for global, basal, and midventricular segments). Apical segments showed a similar but weaker tendency (P=0.015–0.045).

Data reproducibility

All measures of strain and mapping demonstrated moderate-to-excellent reliability for intra-and interobserver variability [intraclass correlation coefficient (ICC) 0.626–0.993]. Table S2 summarizes the ICC values for both intra- and interobserver reproducibility.


Discussion

This study identified a novel, burden-dependent myocardial phenotype in patients with idiopathic PVCs and normal CMR findings. The phenotype is characterized by (I) impaired LV strain, most consistently in global/apical circumferential and radial directions; and (II) elevated native T1 and T2 values, predominantly in the basal and midventricular segments. Collectively, these patterns suggest that the observed alterations may represent concomitant subclinical systolic dysfunction and interstitial changes.

CMR-FT-based strain parameters

Research on the association between CMR-FT strain analysis and PVCs remains limited and inconclusive. A previous study involving 42 patients with PVCs (>10/h) or nonsustained ventricular tachycardia (NSVT) and structurally normal hearts reported altered segmental radial and circumferential strains in the basal and mid-LV segments, alongside reduced longitudinal strain in the mid-LV segments (32). Another CMR-FT analysis reported impaired LV GRS and GCS in patients with frequent PVCs (≥500/24 h) or NSVT, although no statistically significant differences in LV GLS were observed between patients and controls (13). The heterogeneity in these findings likely stems from differences in the patient selection criteria and the multifactorial nature of myocardial strain alterations. Our study sought to clarify this picture by implementing rigorous inclusion criteria and stratifying patients based on 24-hour PVC burden, consistent with established clinical guidelines (4,18).

In this context, our findings both confirm and extend previous observations. We found impaired global LV 2D strain (2D-GRS and GCS) and 3D strain (3D-GCS) in patients with >10,000 PVC/24 h than in controls, which aligns with previous findings of GRS and GCS impairment in patients with frequent PVCs (13). A key point of consistency is the relative sparing of GLS, for which we also observed no significant differences between groups, only a tendency toward a decrease, consistent with previous work (33). This recurring pattern across studies suggests a particular susceptibility of patients with circumferential and radial deformation to PVC-induced dyssynchrony. The mechanism underlying this relationship may involve the intrinsic mechanics of ventricular contraction. In the normal heart, the majority of the myocardium exhibits relatively low longitudinal deformation (34). It has been proposed that asynchronous electrical activation, such as that from PVCs, could introduce disproportional alterations in the timing and pattern of wall thickening and shortening, potentially affecting circumferential and radial strains to a greater extent than longitudinal strains (35-37). Therefore, while strain analysis broadly demonstrates an ability to detect subtle contractile impairments in high-burden PVC patients, the choice of specific strain parameters is critical.

Our analysis identified two potential pathways through which PVC burden may influence myocardial strain. First, we observed that higher PVC burden was weakly correlated with increased LV volume (EDV and ESV), and this volumetric change appeared to contribute to a generalized reduction in myocardial strain. Second, and more notably, the correlation between PVC burden and certain strain parameters—particularly global strains and apical radial/circumferential strain (ARS and ACS)—remained significant after LV volumes were accounted for. This indicates a direct relationship between PVC burden and myocardial dysfunction that is independent of gross structural remodeling. The most pronounced impairments were consistently localized to the apical segment, a region of complex mechanics that may be uniquely vulnerable to the dyssynchronous activation patterns of PVCs, due to its shape, curvature, and fiber architecture (38-40).

Myocardial tissue characteristics of patients with PVCs

In patients with a high PVC burden, the observed reduction in strain accompanied by increased ventricular volumes suggests the presence of subclinical myocardial injury. Such impaired cardiac function has previously been linked to increased myocardial stiffness (41,42). Although LV myocardial replacement fibrosis, detectable by LGE, is a recognized marker of maladaptive remodeling (43), our study focused on the earlier, diffuse tissue changes that precede overt scar formation. Myocardial fibrosis and inflammatory edema are interrelated processes that drive LV remodeling and dysfunction (44). Although elevated T1 or T2 mapping values have been associated with VA in established cardiomyopathies (27,28), the myocardial tissue characteristics in patients with idiopathic PVCs but otherwise normal hearts remain largely unexplored.

In our study, patients with >10,000 PVC/24 h exhibited higher native T1 and T2 values, particularly in the basal segments, as compared to healthy controls. In the context of normal conventional imaging, these variations are consistent with the presence of myocardial interstitial pathology, such as edema or fibrosis. Thus, beyond functional strain analysis, CMR mapping provides preliminary evidence that myocardial involvement in patients with a high burden of PVCs might extend to the microstructural level, potentially revealing subclinical interstitial changes. Although this interpretation requires validation through future studies incorporating histological correlation or ECV measurement, our findings highlight the potential of CMR mapping to identify a novel tissue phenotype in this patient population. The myocardial interstitium, a dynamic compartment comprising stromal cells and extracellular matrix, is critical to postinjury repair. Both reparative and reactive myocardial interstitial fibrosis contribute to cardiac dysfunction and poor outcomes in individuals with nonischemic disease (45). Despite further investigation and histopathological validation being needed to determine the exact nature of the underlying tissue change, our findings nevertheless demonstrate the presence of a quantifiable interstitial abnormality.

Limitations

This study involved several limitations that should be addressed. First, we employed a retrospective, single-center design and vendor-specific methodology, which may limit the generalizability of the findings and require external validation. Second, although the stringent patient selection enhanced the internal validity of the study, the resulting cohort and modest subgroup sizes may not fully represent the broader PVC population. Third, the absence of postablation CMR data precluded assessment of causality, and the potential influence of PVC origin could not be examined due to limited mapping data. Finally, the modest sample size led us to employ partial correlation rather than more complex multivariable models for adjusting confounders. Although this approach supports independence, the stability of more comprehensive models requires validation in larger cohorts.


Conclusions

This exploratory study identified a novel, burden-dependent myocardial phenotype among patients with idiopathic PVCs, characterized by CMR-derived strain impairment and elevated native T1/T2 values. The robust dose-response relationship is consistent with the subclinical myocardial involvement directly related to PVC burden. Given the study’s limitations, including its sample size, the findings are inherently hypothesis-generating. Future longitudinal studies with histological correlation are essential to definitively establish the clinical and prognostic significance of this phenotype.


Acknowledgments

None.


Footnote

Reporting Checklist: The authors have completed the STROBE reporting checklist. Available at https://qims.amegroups.com/article/view/10.21037/qims-2025-1190/rc

Data Sharing Statement: Available at https://qims.amegroups.com/article/view/10.21037/qims-2025-1190/dss

Funding: None.

Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://qims.amegroups.com/article/view/10.21037/qims-2025-1190/coif). The authors have no conflicts of interest to declare.

Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. This study received approval from the Ethics Review Committee of Shanghai General Hospital (No. 2023364) and was exempt from the requirement for informed consent due to the retrospective nature of the study.

Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.


References

  1. Sirichand S, Killu AM, Padmanabhan D, Hodge DO, Chamberlain AM, Brady PA, Kapa S, Noseworthy PA, Packer DL, Munger TM, Gersh BJ, McLeod CJ, Shen WK, Cha YM, Asirvatham SJ, Friedman PA, Mulpuru SK. Incidence of Idiopathic Ventricular Arrhythmias: A Population-Based Study. Circ Arrhythm Electrophysiol 2017;10:e004662. [Crossref] [PubMed]
  2. Kostis JB, McCrone K, Moreyra AE, Gotzoyannis S, Aglitz NM, Natarajan N, Kuo PT. Premature ventricular complexes in the absence of identifiable heart disease. Circulation 1981;63:1351-6. [Crossref] [PubMed]
  3. Higuchi K, Bhargava M. Management of premature ventricular complexes. Heart 2022;108:565-72. [Crossref] [PubMed]
  4. Pedersen CT, Kay GN, Kalman J, Borggrefe M, Della-Bella P, Dickfeld T, Dorian P, Huikuri H, Kim YH, Knight B, Marchlinski F, Ross D, Sacher F, Sapp J, Shivkumar K, Soejima K, Tada H, Alexander ME, Triedman JK, Yamada T, Kirchhof P, Lip GY, Kuck KH, Mont L, Haines D, Indik J, Dimarco J, Exner D, Iesaka Y, Savelieva I. EHRA/HRS/APHRS expert consensus on ventricular arrhythmias. Europace 2014;16:1257-83. [Crossref] [PubMed]
  5. Ahmad MI, Soliman MZ, Soliman EZ. Relationship between premature ventricular complexes and stroke mortality in the general population. J Electrocardiol 2023;77:41-5. [Crossref] [PubMed]
  6. Ataklte F, Erqou S, Laukkanen J, Kaptoge S. Meta-analysis of ventricular premature complexes and their relation to cardiac mortality in general populations. Am J Cardiol 2013;112:1263-70. [Crossref] [PubMed]
  7. Scorza R, Jonsson M, Corander JM, Rosenqvist M, Frykman V. Prognostic impact of morphology and duration of premature ventricular contractions in a population without structural heart disease. Ann Noninvasive Electrocardiol 2023;28:e13067. [Crossref] [PubMed]
  8. Dukes JW, Dewland TA, Vittinghoff E, Mandyam MC, Heckbert SR, Siscovick DS, Stein PK, Psaty BM, Sotoodehnia N, Gottdiener JS, Marcus GM. Ventricular Ectopy as a Predictor of Heart Failure and Death. J Am Coll Cardiol 2015;66:101-9. [Crossref] [PubMed]
  9. Refaat MM, Gharios C, Moorthy MV, Abdulhai F, Blumenthal RS, Jaffa MA, Mora S. Exercise-Induced Ventricular Ectopy and Cardiovascular Mortality in Asymptomatic Individuals. J Am Coll Cardiol 2021;78:2267-77. [Crossref] [PubMed]
  10. Huizar JF, Tan AY, Kaszala K, Ellenbogen KA. Clinical and translational insights on premature ventricular contractions and PVC-induced cardiomyopathy. Prog Cardiovasc Dis 2021;66:17-27. [Crossref] [PubMed]
  11. Qin M, Song ZL, Zhu XY, Zhang Y, Jiang WF, Wu SH, Shen XY, Liu T, Liu X. Temporal and Spatial Changes of Proarrhythmic Substrate in Premature Ventricular Contraction-Induced Cardiomyopathy. JACC Clin Electrophysiol 2023;9:173-88. [Crossref] [PubMed]
  12. Ailoaei Ș, Șorodoc L, Ureche C, Sîtari N, Sandu G, Ceasovschih A, Grecu M, Sascău RA, Stătescu C. Predictors of Myocardial Fibrosis Detected by CMR in Patients With Idiopathic PVCs. Ann Noninvasive Electrocardiol 2024;29:e70039. [Crossref] [PubMed]
  13. Nikolaidou C, Kotanidis CP, Wijesurendra R, Leal-Pelado J, Kouskouras K, Vassilikos VP, Karvounis H, Ntusi N, Antoniades C, Neubauer S, Karamitsos TD. Cardiac Magnetic Resonance to Detect the Underlying Substrate in Patients with Frequent Idiopathic Ventricular Arrhythmias. Diagnostics (Basel) 2021;11:1109. [Crossref] [PubMed]
  14. Al-Khatib SM, Stevenson WG, Ackerman MJ, Bryant WJ, Callans DJ, Curtis AB, Deal BJ, Dickfeld T, Field ME, Fonarow GC, Gillis AM, Granger CB, Hammill SC, Hlatky MA, Joglar JA, Kay GN, Matlock DD, Myerburg RJ, Page RL. 2017 AHA/ACC/HRS Guideline for Management of Patients With Ventricular Arrhythmias and the Prevention of Sudden Cardiac Death: A Report of the American College of Cardiology/American Heart Association Task Force on Clinical Practice Guidelines and the Heart Rhythm Society. Circulation 2018;138:e272-e391. Erratum in: Circulation 2018;138:e419-e420. [Crossref] [PubMed]
  15. Marcus GM. Evaluation and Management of Premature Ventricular Complexes. Circulation 2020;141:1404-18. [Crossref] [PubMed]
  16. Priori SG, Blomström-Lundqvist C, Mazzanti A, Blom N, Borggrefe M, Camm J, Elliott PM, Fitzsimons D, Hatala R, Hindricks G, Kirchhof P, Kjeldsen K, Kuck KH, Hernandez-Madrid A, Nikolaou N, Norekvål TM, Spaulding C, Van Veldhuisen DJTask Force for the Management of Patients with Ventricular Arrhythmias and the Prevention of Sudden Cardiac Death of the European Society of Cardiology (ESC). 2015 ESC Guidelines for the management of patients with ventricular arrhythmias and the prevention of sudden cardiac death: The Task Force for the Management of Patients with Ventricular Arrhythmias and the Prevention of Sudden Cardiac Death of the European Society of Cardiology (ESC)Endorsed by: Association for European Paediatric and Congenital Cardiology (AEPC). Europace 2015;17:1601-87.
  17. Koester C, Ibrahim AM, Cancel M, Labedi MR. The Ubiquitous Premature Ventricular Complex. Cureus 2020;12:e6585. [Crossref] [PubMed]
  18. Arnar DO, Mairesse GH, Boriani G, Calkins H, Chin A, Coats A, et al. Management of asymptomatic arrhythmias: a European Heart Rhythm Association (EHRA) consensus document, endorsed by the Heart Failure Association (HFA), Heart Rhythm Society (HRS), Asia Pacific Heart Rhythm Society (APHRS), Cardiac Arrhythmia Society of Southern Africa (CASSA), and Latin America Heart Rhythm Society (LAHRS). Europace 2019;21:844-845. [Crossref] [PubMed]
  19. Mahida S, Sacher F, Dubois R, Sermesant M, Bogun F, Haïssaguerre M, Jaïs P, Cochet H. Cardiac Imaging in Patients With Ventricular Tachycardia. Circulation 2017;136:2491-507. [Crossref] [PubMed]
  20. Disertori M, Rigoni M, Pace N, Casolo G, Masè M, Gonzini L, Lucci D, Nollo G, Ravelli F. Myocardial Fibrosis Assessment by LGE Is a Powerful Predictor of Ventricular Tachyarrhythmias in Ischemic and Nonischemic LV Dysfunction: A Meta-Analysis. JACC Cardiovasc Imaging 2016;9:1046-55. [Crossref] [PubMed]
  21. Muser D, Nucifora G, Muser D, Nucifora G, Pieroni M, Castro SA, et al. Prognostic Value of Nonischemic Ringlike Left Ventricular Scar in Patients With Apparently Idiopathic Nonsustained Ventricular Arrhythmias. Circulation 2021;143:1359-73. [Crossref] [PubMed]
  22. Ghannam M, Siontis KC, Kim MH, Cochet H, Jais P, Eng MJ, Attili A, Sharaf-Dabbagh G, Latchamsetty R, Jongnarangsin K, Morady F, Bogun F. Risk stratification in patients with frequent premature ventricular complexes in the absence of known heart disease. Heart Rhythm 2020;17:423-30. [Crossref] [PubMed]
  23. Andreini D, Dello Russo A, Pontone G, Mushtaq S, Conte E, Perchinunno M, et al. CMR for Identifying the Substrate of Ventricular Arrhythmia in Patients With Normal Echocardiography. JACC Cardiovasc Imaging 2020;13:410-21. [Crossref] [PubMed]
  24. Muser D, Santangeli P, Castro SA, Casado Arroyo R, Maeda S, Benhayon DA, et al. Risk Stratification of Patients With Apparently Idiopathic Premature Ventricular Contractions: A Multicenter International CMR Registry. JACC Clin Electrophysiol 2020;6:722-35. [Crossref] [PubMed]
  25. Warnica W, Al-Arnawoot A, Stanimirovic A, Thavendiranathan P, Wald RM, Pakkal M, Karur GR, Wintersperger BJ, Rac V, Hanneman K. Clinical Impact of Cardiac MRI T1 and T2 Parametric Mapping in Patients with Suspected Cardiomyopathy. Radiology 2022;305:319-26. [Crossref] [PubMed]
  26. Muser D, Liang JJ, Witschey WR, Pathak RK, Castro S, Magnani S, Zado ES, Garcia FC, Desjardins B, Callans DJ, Frankel DS, Marchlinski FE, Santangeli P. Ventricular arrhythmias associated with left ventricular noncompaction: Electrophysiologic characteristics, mapping, and ablation. Heart Rhythm 2017;14:166-75. [Crossref] [PubMed]
  27. Chen Z, Sohal M, Voigt T, Sammut E, Tobon-Gomez C, Child N, Jackson T, Shetty A, Bostock J, Cooklin M, O'Neill M, Wright M, Murgatroyd F, Gill J, Carr-White G, Chiribiri A, Schaeffter T, Razavi R, Rinaldi CA. Myocardial tissue characterization by cardiac magnetic resonance imaging using T1 mapping predicts ventricular arrhythmia in ischemic and non-ischemic cardiomyopathy patients with implantable cardioverter-defibrillators. Heart Rhythm 2015;12:792-801. [Crossref] [PubMed]
  28. Baig M, Galazka P, Dakwar O, Syed SA, Sawlani R, Shahir K, Muthukumar L, Jain R, Jan MF, Port S, Khandheria BK, Jahangir A, Tajik AJ. Prevalence of myocardial edema with T2 mapping in hypertrophic cardiomyopathy. Journal of the American College of Cardiology 2021;77:1303-1303.
  29. Halliday BP, Senior R, Pennell DJ. Assessing left ventricular systolic function: from ejection fraction to strain analysis. Eur Heart J 2021;42:789-97. [Crossref] [PubMed]
  30. Gatti M, Palmisano A, Esposito A, Fiore S, Monti CB, Andreis A, Pistelli L, Vergara P, Bergamasco L, Giustetto C, De Cobelli F, Fonio P, Faletti R. Feature tracking myocardial strain analysis in patients with bileaflet mitral valve prolapse: relationship with LGE and arrhythmias. Eur Radiol 2021;31:7273-82. [Crossref] [PubMed]
  31. Pu C, Fei J, Lv S, Wu Y, He C, Guo D, Mabombo PU, Chooah O, Hu H. Global Circumferential Strain by Cardiac Magnetic Resonance Tissue Tracking Associated With Ventricular Arrhythmias in Hypertrophic Cardiomyopathy Patients. Front Cardiovasc Med 2021;8:670361. [Crossref] [PubMed]
  32. Tang X, Yu S, Yu Y, Ren H, Li S, Zhou L, Yang Z, Wu H, Zhou W, Gong L. Left ventricular myocardial strain in ventricular arrhythmia without structural heart disease using cardiac magnetic resonance. Am J Transl Res 2017;9:3006-16.
  33. Scorza R, Shahgaldi K, Rosenqvist M, Frykman V. Evaluation of patients with high burden of premature ventricular contractions by comprehensive transthoracic echocardiography. Int J Cardiol Heart Vasc 2022;42:101124. [Crossref] [PubMed]
  34. Waldman LK, Fung YC, Covell JW. Transmural myocardial deformation in the canine left ventricle. Normal in vivo three-dimensional finite strains. Circ Res 1985;57:152-63. [Crossref] [PubMed]
  35. Waldman LK, Covell JW. Effects of ventricular pacing on finite deformation in canine left ventricles. Am J Physiol 1987;252:H1023-30. [Crossref] [PubMed]
  36. Helm RH, Leclercq C, Faris OP, Ozturk C, McVeigh E, Lardo AC, Kass DA. Cardiac dyssynchrony analysis using circumferential versus longitudinal strain: implications for assessing cardiac resynchronization. Circulation 2005;111:2760-7. [Crossref] [PubMed]
  37. Yao J, Yang R, Xu D, Zhuang Y, Yong Y, Cao K. Circumferential myocardial contraction patterns in patients with idiopathic frequent premature ventricular complexes from the right ventricular outflow tract. Int J Cardiol 2013;166:166-72. [Crossref] [PubMed]
  38. Prinzen FW, Augustijn CH, Arts T, Allessie MA, Reneman RS. Redistribution of myocardial fiber strain and blood flow by asynchronous activation. Am J Physiol 1990;259:H300-8. [Crossref] [PubMed]
  39. Delhaas T, Arts T, Prinzen FW, Reneman RS. Regional fibre stress-fibre strain area as an estimate of regional blood flow and oxygen demand in the canine heart. J Physiol 1994;477:481-96. [Crossref] [PubMed]
  40. Prinzen FW, Hunter WC, Wyman BT, McVeigh ER. Mapping of regional myocardial strain and work during ventricular pacing: experimental study using magnetic resonance imaging tagging. J Am Coll Cardiol 1999;33:1735-42. [Crossref] [PubMed]
  41. Topaloglu S, Aras D, Cagli K, Yildiz A, Cagirci G, Cay S, Gunel EN, Baser K, Baysal E, Boyaci A, Korkmaz S. Evaluation of left ventricular diastolic functions in patients with frequent premature ventricular contractions from right ventricular outflow tract. Heart Vessels 2007;22:328-34. [Crossref] [PubMed]
  42. van Huls van Taxis CF, Piers SR, de Riva Silva M, Dekkers OM, Pijnappels DA, Schalij MJ, Wijnmaalen AP, Zeppenfeld K. Fatigue as Presenting Symptom and a High Burden of Premature Ventricular Contractions Are Independently Associated With Increased Ventricular Wall Stress in Patients With Normal Left Ventricular Function. Circ Arrhythm Electrophysiol 2015;8:1452-9. [Crossref] [PubMed]
  43. Patel AR, Kramer CM. Role of Cardiac Magnetic Resonance in the Diagnosis and Prognosis of Nonischemic Cardiomyopathy. JACC Cardiovasc Imaging 2017;10:1180-93. [Crossref] [PubMed]
  44. Emrich T, Hahn F, Fleischmann D, Halfmann MC, Düber C, Varga-Szemes A, Escher F, Pefani E, Münzel T, Schultheiss HP, Kreitner KF, Wenzel P. T1 and T2 mapping to detect chronic inflammation in cardiac magnetic resonance imaging in heart failure with reduced ejection fraction. ESC Heart Fail 2020;7:2544-52. [Crossref] [PubMed]
  45. Díez J, González A, Kovacic JC. Myocardial Interstitial Fibrosis in Nonischemic Heart Disease, Part 3/4: JACC Focus Seminar. J Am Coll Cardiol 2020;75:2204-18. [Crossref] [PubMed]
Cite this article as: Yang J, Zhang Y, Yu Q, Hu X, Luo Q, Zhang L. Comprehensive cardiac magnetic resonance imaging in patients with idiopathic premature ventricular contractions: role of feature-tracking strain and T1/T2 mapping in detecting subclinical myocardial dysfunction. Quant Imaging Med Surg 2026;16(2):164. doi: 10.21037/qims-2025-1190

Download Citation