Subjects
In the present study, we retrospectively analyzed the medical records and clinical characteristics of CTD patients diagnosed with PAH by RHC hospitalized in the First Affiliated Hospital of Nanjing Medical University between January 1, 2010 and September 1, 2020. Inclusion criteria were listed as follows: (1) a confirmed diagnosis of PAH according to the 2015 guidelines[1]; (2) meeting the diagnostic criteria for each subcategory of CTD: Systemic lupus erythematosus (SLE) was diagnosed according to 1997 American Rheumatism Association (ACR) criteria[28], Sjogren’s syndrome (SS) was diagnosed according to 2002 international classification criteria[29], systemic sclerosis (SSc) was defined according to 1980 ACR criteria[30]; mixed CTD (MCTD) was defined by Sharp’s criteria[31]; (3) having undergone chest CT within 1 month before or after RHC and the patients’ condition was stable during the two examination intervals. The exclusion criteria included: (1) congenital heart disease; (2) other causes of precapillary PH (such as PH due to respiratory diseases, chronic thromboembolic PH, or other miscellaneous causes of PAH); (3) significant valvular heart disease of more than moderate to severe or LV ejection fraction < 50% diagnosed by echocardiography; (4) coexisting pulmonary conditions on computed tomography scan affected quantitative CT measurements: moderate or severe pulmonary interstitial fibrosis, current pneumonia, massive pleural effusion.
This study was approved by the Medical Ethics Committee of the First Affiliated Hospital of Nanjing Medical University (number: 2018-SR-333). As all variables were obtained retrospectively from available clinical data, the need for patients to sign informed consent was waived.
CT measurement of pulmonary vessels
The diameter of the MPA and AAO was measured at the level of MPA bifurcation in its maximum dimension (Figure 2A)[32].
To measure the CSA of small pulmonary vessels, three slices were selected. Three plain CT axial slices: 1 cm above the upper margin of the aortic arch (upper slice), 1 cm below the carina (middle slice), and 1cm below the right inferior pulmonary vein (lower slice).
Subsequently, images were analyzed with semiautomatic quantitative image-processing Image J software (version 1.48; National Institutes of Health, Bethesda, MD, USA). %CSA<5 we calculated according to the method reported by Matsuoka et al[23]. On each CT slice, %CSA<5 was obtained with the “Analyze Particles” function to count and measure objects on binary images, the number of vessels at a specified size and the CSA of each size range were obtained. Notably, the vessels that ran obliquely or parallel to the slice were excluded using the “Circularity” function in Image J (Figure 2B-D)[24].
Clinical outcome
The primary endpoint was a composite clinical worsening endpoint (including all-cause mortality, worsening World Health Organization functional class, ≥15% reduction in 6-MWD, all-cause hospitalization, or the introduction of parenteral prostacyclin analog therapy)[33]. The time of follow-up was calculated as the time from RHC examination to the end of the study (September 1, 2020) or to the composite endpoint of clinical deterioration, whichever came first. The follow-up data were obtained from hospital records. All patients were contacted to reconfirm survival status by telephone personal interview of the patient or family members at the end of the study.
Statistical Analyses
All statistical analyses were performed using Statistical Package for the Social Sciences (SPSS) (ver 25.0, International Business Machines, Inc. Armonk, New York, USA). Qualitative data were expressed as frequency (percentage) and compared using χ2 test or Fisher exact test. Quantitative data were expressed as mean ± standard or median (interquartile range [IQR]) and compared among groups by Student’s-test or Mann-Whitney U test. The correlation between hemodynamic parameters and quantitative CT parameters was assessed by Pearson’s correlation coefficient. Receiver-operating characteristic curves were generated to assess the effectiveness of the %CSA<5 and MPA as targeting risk factors to predict the endpoint by evaluating the sensitivity and specificity of the scales. Results were expressed in terms of area under the curve (AUC) and 95% confidence interval for this area. We define the corresponding cutoff for each variable by the difference maximization method, and the patients were further divided into 3 groups according to the cutoff of the %CSA<5 and MPA. We use “pulmonary vasculature metrics” to define classification based on the cutoff values of %CSA<5 and MPA. The prognostic value of selected baseline parameters was tested using Cox’s univariate proportional hazards regression analysis, and the variables that were significant in the univariate model were then entered into a multivariate Cox model. The results were expressed as hazard ratios (HR) with 95% CI. Kaplan-Meier method was used to calculate time-to-event function among 3 groups and differences were assessed using the log-rank statistic. All P-values were two-sided and a P-value < 0.05 was considered statistically significant.