China's prevention and control experience of echinococcosis: A 19-year retrospective

Echinococcosis poses a significant threat to public health. The Chinese government has implemented prevention and control measures to mitigate the impact of the disease. By analyzing data from the Chinese Center for Disease Control and Prevention and the State Council of the People's Republic of China, we found that implementation of these measures has reduced the infection rate by nearly 50% between 2004 to 2022 (from 0.3975 to 0.1944 per 100,000 person-years). Nonetheless, some regions still bear a significant disease burden, and lack of detailed information limites further evaluation of the effects on both alveolar and cystic echinococcosis. Our analysis supports the continuing implementation of these measures and suggests that enhanced wildlife management, case-based strategies, and surveillance systems will facilitate disease control.

Dosimetry and treatment efficiency of SBRT using TaiChiB radiotherapy system for two-lung lesions with one overlapping organs at risk

Purpose

This study aims to assess the dosimetry and treatment efficiency of TaiChiB-based Stereotactic Body Radiotherapy (SBRT) plans applying to treat two-lung lesions with one overlapping organs at risk.

Methods

For four retrospective patients diagnosed with two-lung lesions each patient, four treatment plans were designed including Plan Edge, TaiChiB linac-based, RGS-based, and a linac-RGS hybrid (Plan TCLinac, Plan TCRGS, and Plan TCHybrid). Dosimetric metrics and beam-on time were employed to evaluate and compare the TaiChiB-based plans against Plan Edge.

Results

For Conformity Index (CI), Plan TCRGS outperformed all other plans with an average CI of 1.06, as opposed to Plan Edge's 1.33. Similarly, for R50 %, Plan TCRGS was superior with an average R50 % of 3.79, better than Plan Edge's 4.28. In terms of D2 cm, Plan TCRGS also led with an average of 48.48%, compared to Plan Edge's 56.25%. For organ at risk (OAR) sparing, Plan TCRGS often displayed the lowest dosimetric values, notably for the spinal cord (Dmax 5.92 Gy) and lungs (D1500cc 1.00 Gy, D1000cc 2.61 Gy, V10 Gy 15.14%). However, its high Dmax values for the heart and great vessels sometimes exceeded safety thresholds. Plan TCHybrid presented a balanced approach, showing doses comparable to or better than Plan Edge without crossing safety limits. In terms of beam-on time, Plan TCLinac emerged as the most efficient treatment option in three out of four cases, followed closely by Plan Edge in one case. Plan TCRGS, despite its dosimetric advantages, was the least efficient, recording notably longer beam-on times, with a peak at 33.28 minutes in Case 2.

Conclusion

For patients with two-lung lesions treated by SBRT whose one lesion overlaps with OARs, the Plan TCHybrid delivered by TaiChiB digital radiotherapy system can be recommended as a clinical option.

AS-NeSt: A Novel 3D Deep Learning Model for Radiation Therapy Dose Distribution Prediction in Esophageal Cancer Treatment With Multiple Prescriptions

PURPOSE

Implementing artificial intelligence technologies allows for the accurate prediction of radiation therapy dose distributions, enhancing treatment planning efficiency. However, esophageal cancers present unique challenges because of tumor complexity and diverse prescription types. Additionally, limited data availability hampers the effectiveness of existing artificial intelligence models. This study developed a deep learning model, trained on a diverse data set of esophageal cancer prescriptions, to improve dose prediction accuracy.

METHODS AND MATERIALS

We retrospectively collected data from 530 patients with esophageal cancer, including single-target and simultaneous integrated boost prescriptions, for model building. The proposed Asymmetric ResNeSt (AS-NeSt) model features novel 3-dimensional (3D) ResNeSt blocks and an asymmetrical architecture. We constructed a loss function targeting global and local doses and validated the model's performance against existing alternatives. Model-assisted experiments were used to validate its clinical benefits.

RESULTS

The AS-NeSt model maintained an absolute prediction error below 5% for each dosimetric metric. The average Dice similarity coefficient for isodose volumes was 0.93. The model achieved an average relative prediction error of 2.02%, statistically lower than Hierarchically Densely Connected U-net (4.17%), DoseNet (2.35%), and Densely Connected Network (3.65%). It also demonstrated significantly fewer parameters and shorter prediction times. Clinically, the AS-NeSt model raised physicians' ability to accurately preassess appropriate treatment methods before planning from 95.24% to 100%, reduced planning time by over 61% for junior dosimetrists and 52% for senior dosimetrists, and decreased both inter- and intra-dosimetrist discrepancies by more than 50%.

CONCLUSIONS

The AS-NeSt model, developed with innovative 3D ResNeSt blocks and an asymmetrical encoder-decoder structure, has been validated using clinical esophageal cancer patient data. It accurately predicts 3D dose distributions for various prescriptions, including simultaneous integrated boost, showing potential to improve the management of esophageal cancer treatment in a clinical setting.

Optimization of isocenter position for multiple brain metastases single-isocenter stereotactic radiosurgery to minimize dosimetric variations due to rotational uncertainty

PURPOSE

This paper studied a novel calculation framework that can determine the optimal value isocenter position of single isocenter SRS treatment plan for multiple brain metastases, in order to minimize the dosimetric variations caused by rotational uncertainty.

MATERIALS AND METHODS

21 patients with 2-4 GTVswho received SRS treatment for multiple brain metastases in our institution were selected for the retrospective study. The PTVwas obtained by expanding GTV 1 mm isotropic margin. We studied a stochastic optimization framework, which determined the optimal value isocenter location by maximizing the average target dose coverageC_target,meanwith a rotation error of no more than 1°. We evaluated the performance of the optimal isocenter by comparing theC_target,meanand average dice similarity coefficient (DSC)with the optimal value and the center of mass (CM) respectively as the treatment isocenter. The extra PTV margin to achieve 100% target dose coverage was calculated by our framework.

RESULTS

Compared to the CM method, the optimal value isocenter method increased the averageC_target,meanof all targets from 97.0% to 97.7%and the average DSC from 0.794to 0.799. Throughout all the cases, the average extra PTV margin to obtain full target dose coverage was 0.7 mmwhen using the optimal value isocenter as the treatment isocenter.

CONCLUSION

We studied a novel computational framework using stochastic optimization to determine the optimal isocenter position of SRS treatment plan for multiple brain metastases. At the same time, our framework gave the extra PTV margin to obtain full target dose coverage.

[A preliminary study on the tear matrix metalloproteinase 9 point-of-care assay using a domestic kit]

Objective: To compare the point-of-care assays for tear matrix metalloproteinase 9 (MMP-9) using domestic and InflammaDry kits, and to evaluate the feasibility of diagnosing dry eye with the domestic kit. Methods: It was a cross-sectional study. Thirty dry eye patients and 30 age-and sex-matched normal volunteers were continuously enrolled in this cross-sectional study from June 2022 to July 2022. Both domestic and InflammaDry kits were used to detect the tear MMP-9 levels. The positive rates were recorded for qualitative analysis, and the gray ratios of bands (the gray value of detection bands to that of control bands) were collected for quantitative analysis. The correlations of MMP-9 levels with age, ocular surface disease index, fluorescence tear break-up time, tear meniscus height, Schirmer's Ⅰ test score, corneal fluorescein staining score, and meibomian gland dropout were analyzed. The Mann-Whitney U test, paired Chi-square test, Kappa test, and Spearman's correlation coefficient were used for statistical analysis. Results: There were 14 males and 16 females (30 eyes) in the control group, and their age was (39.37±19.55) years. In the dry eye group, 11 males and 19 females (30 eyes), aged (46.87±17.85) years, had moderate to severe dry eye. The positive rates of MMP-9 in tear fluid were significantly different between dry eye patients (InflammaDry: 86.67%; domestic kit: 70.00%) and controls (InflammaDry: 16.67%, P<0.001; domestic kit: 6.67%, P<0.001). Although the sensitivity of the domestic kit was lower than that of the InflammaDry kit (70.0% vs. 86.7%, P=0.001), the specificity was higher (93.3% vs. 83.3%, P=0.001). In dry eye patients, the positive coincidence rate was 80.7% (21/26), the negative coincidence rate was 100% (4/4), and the total coincidence rate was 83.3% (25/30), with no significant difference between the two kits (McNemar test: χ²=3.20, P>0.05), and the results of both kits were consistent (Kappa=0.53, P=0.001). The Spearman's correlation coefficient showed the gray ratios using both kits were positively correlated with the corneal fluorescein staining score (InflammaDry: ρ=0.48, P<0.05; domestic kit: ρ=0.52, P=0.003). Conclusion: The performances of the domestic and InflammaDry kits are consistent in the point-of-care assay for tear MMP-9, and the domestic kit has lower sensitivity but higher specificity.

Image-based features in machine learning to identify delivery errors and predict error magnitude for patient-specific IMRT quality assurance

OBJECTIVE

To identify delivery error type and predict associated error magnitude by image-based features using machine learning (ML).

METHODS

In this study, a total of 40 thoracic plans (including 208 beams) were selected, and four error types with different magnitudes were introduced into the original plans, including 1) collimator misalignment (COLL), 2) monitor unit (MU) variation, 3) systematic multileaf collimator misalignment (MLCS), and 4) random MLC misalignment (MLCR). These dose distributions of portal dose predictions for the original plans were defined as the reference dose distributions (RDD), while those for the error-introduced plans were defined as the error-introduced dose distributions (EDD). Both distributions were calculated for all beams with portal dose image prediction (PDIP). Besides, 14 image-based features were extracted from RDD and EDD of portal dose predictions to obtain the feature vectors. In addition, a random forest was adopted for the multiclass classification task, and regression prediction for error magnitude.

RESULTS

The top five features extracted with the highest weight included 1) the relative displacement in the x direction, 2) the ratio of the absolute minimum residual error to the maximal RDD value, 3) the product of the maximum and minimum residuals, 4) the ratio of the absolute maximum residual error to the maximal RDD value, and 5) the ratio of the absolute mean residual value to the maximal RDD value. The relative displacement in the x direction had the highest weight. The overall accuracy of the five-class classification model was 99.85% for the validation set and 99.30% for the testing set. This model could be applied to the classification of the error-free plan, COLL, MU, MLCS, and MLCR with an accuracy of 100%, 98.4%, 99.9%, 98.0%, and 98.3%, respectively. MLCR had the worst performance in error magnitude prediction (70.1-96.6%), while others had better performance in error magnitude prediction (higher than 93%). In the error magnitude prediction, the mean absolute error (MAE) between predicted error magnitude and actual error ranged from 0.03 to 0.33, with the root mean squared error (RMSE) varying from 0.17 to 0.56 for the validation set. The MAE and RMSE ranged from 0.03 to 0.50 and 0.44 to 0.59 for the test set, respectively.

CONCLUSION

It could be demonstrated in this study that the image-based features extracted from RDD and EDD can be employed to identify different types of delivery errors and accurately predict error magnitude with the assistance of ML techniques. They can be used to associate traditional gamma analysis with clinically based analysis for error classification and magnitude prediction in patient-specific IMRT quality assurance.

[Hsa_circ_0000670 promoted the progression of gastric cancer through the miR-515-5p/SIX1 molecular axis]

Objective: To explore whether hsa_circ_0000670 promotes the progression of gastric cancer by regulating the miR-515-5p/SIX1 molecular axis. Methods: The gastric cancer and adjacent normal tissues of 35 gastric cancer patients admitted to Rugao Hospital Affiliated to Nantong University from 2014 to 2015 were collected. The expression levels of circ_0000670, miR-515-5p and Sine oculis homeobox 1 (SIX1) in gastric cancer tissues and cells were detected by quantitative real-time polymerase chain reaction (qRT-PCR) and western blot. The correlations between circ_0000670 and miR-515-5p, miR-515-5p and SIX1, circ_0000670 and SIX1 were analyzed by the Pearson method. Patients were divided into low circ_0000670 expression group (17 cases) and high circ_0000670 expression group (18 cases) based on the median of circ_0000670 expression level, and Kaplan-Meier was used to analyze the 5-year survival of patients. Cell proliferation was assessed via clone formation assay. Cell cycle and apoptosis were detected by flow cytometry. Wound healing and Transwell assays were used to detect cell migration and invasion ability. The targeting relationship between miR-515-5p and circ_0000670 or SIX1 was confirmed by the dual luciferase reporter assay. Nude mice were injected into HGC-27 cells transfected with sh-NC or sh-circ_0000670, and the volume and weight of the transplanted tumor were measured, also, the levels of circ_0000670, miR-515-5p and SIX1 in the transplanted tumor tissue were detected. Results: The expression levels of circ_0000670 and SIX1 in gastric cancer tissues and cell lines were significantly increased (P<0.05), while the expression levels of miR-515-5p were significantly decreased (P<0.05). The survival rate of patients in the low circ_0000670 expression group (82.4%) was significantly higher than that in the high circ_0000670 expression group (28.7%, P=0.034). Circ_0000670 was negatively correlated with miR-515-5p (r=-0.846, P<0.001), and miR-515-5p was negatively correlated with SIX1 (r=-0.615, P<0.001), but circ_0000670 was positively correlated with SIX1 (r=0.814, P<0.001). Transfection of si-circ_0000670 or miR-515-5p mimic could significantly reduce the number of clone-forming cells, migration distance, migration and invasion cells (P<0.05), and increase the ratio of G(0)/G(1) phase cells, apoptosis rate and the protein level of E-cadherin (P<0.05), decreased the proportion of S-phase cells and the protein level of Vimentin (P<0.05). The dual luciferase report assay confirmed that circ_0000670 could target miR-515-5p, and miR-515-5p could bind to SIX1. Co-transfection of si-circ_0000670 and miR-515-5p inhibitor could significantly attenuate the effects of si-circ_0000670 on cell proliferation, migration, invasion, cell cycle and apoptosis (P<0.05). Co-transfection of miR-515-5p mimic and pcDNA-SIX1 could significantly reduce the effects of miR-515-5p mimic on cell proliferation, migration, invasion, cell cycle and apoptosis (P<0.05). Compared with the sh-NC group [volume=(596.20±125.46) mm(3) and weight=(538.00±114.39) g], the volume and weight of transplanted tumors in the sh-circ_0000670 group [volume=(299.20±47.58) mm 3 and weight=(289.80±48.73 g)] were significantly reduced (P<0.05), the expression levels of circ_0000670 and SIX1 were significantly reduced (P<0.05), and the expression level of miR-515-5p was significantly increased (P<0.05). Conclusion: Knockdown of circ_0000670 could inhibit cell proliferation, migration, invasion of gastric cancer cells, induce cell cycle arrest in G(0)/G(1) phase and promote cell apoptosis by regulating the miR-515-5p/SIX1 axis.

Improving CBCT image quality to the CT level using RegGAN in esophageal cancer adaptive radiotherapy

OBJECTIVE

This study aimed to improve the image quality and CT Hounsfield unit accuracy of daily cone-beam computed tomography (CBCT) using registration generative adversarial networks (RegGAN) and apply synthetic CT (sCT) images to dose calculations in radiotherapy.

METHODS

The CBCT/planning CT images of 150 esophageal cancer patients undergoing radiotherapy were used for training (120 patients) and testing (30 patients). An unsupervised deep-learning method, the 2.5D RegGAN model with an adaptively trained registration network, was proposed, through which sCT images were generated. The quality of deep-learning-generated sCT images was quantitatively compared to the reference deformed CT (dCT) image using mean absolute error (MAE), root mean square error (RMSE) of Hounsfield units (HU), and peak signal-to-noise ratio (PSNR). The dose calculation accuracy was further evaluated for esophageal cancer radiotherapy plans, and the same plans were calculated on dCT, CBCT, and sCT images.

RESULTS

The quality of sCT images produced by RegGAN was significantly improved compared to the original CBCT images. ReGAN achieved image quality in the testing patients with MAE sCT vs. CBCT: 43.7 ± 4.8 vs. 80.1 ± 9.1; RMSE sCT vs. CBCT: 67.2 ± 12.4 vs. 124.2 ± 21.8; and PSNR sCT vs. CBCT: 27.9 ± 5.6 vs. 21.3 ± 4.2. The sCT images generated by the RegGAN model showed superior accuracy on dose calculation, with higher gamma passing rates (93.3 ± 4.4, 90.4 ± 5.2, and 84.3 ± 6.6) compared to original CBCT images (89.6 ± 5.7, 85.7 ± 6.9, and 72.5 ± 12.5) under the criteria of 3 mm/3%, 2 mm/2%, and 1 mm/1%, respectively.

CONCLUSION

The proposed deep-learning RegGAN model seems promising for generation of high-quality sCT images from stand-alone thoracic CBCT images in an efficient way and thus has the potential to support CBCT-based esophageal cancer adaptive radiotherapy.

Adherence to the Dietary Approaches to Stop Hypertension (DASH) Eating Pattern Reduces the Risk of Head and Neck Cancer in American Adults Aged 55 Years and Above: A Prospective Cohort Study

OBJECTIVES

Dietary Approaches to Stop Hypertension (DASH) pattern has been found to aid in the reduction of obesity, oxidative stress, and chronic inflammation, which are all strongly linked to the development of head and neck cancer (HNC). Nevertheless, no epidemiological studies have investigated the association between this dietary pattern and HNC risk. This study was conducted with the purpose of bridging this gap in knowledge.

DESIGN

A prospective cohort study involving 98,459 American adults aged 55 years and older.

SETTING AND PARTICIPANTS

Data were drawn from the Prostate, Lung, Colorectal, and Ovarian (PLCO) Trial. In the present study, participants with dependable energy intake data who furnished baseline and dietary history information were identified as the study population.

METHODS

Diet was assessed by food frequency questionnaires and the DASH score was calculated to assess each participant's adherence to DASH eating pattern. Cox proportional hazards models were used to calculate multivariable adjusted hazard ratios (HRs) and 95% confidence intervals (CIs) for the occurrence of HNC. To visualize the variation in cancer risk for HNC and its subtypes across the entire spectrum of DASH scores, restricted cubic spline plots were utilized. Additionally, a series of predefined subgroup analyses were performed to identify potential effect modifiers, and several sensitivity analyses were conducted to assess the stability of the findings.

RESULTS

During a follow-up period of 871,879.6 person-years, 268 cases of HNC were identified, comprising 161 cases pertaining to oral cavity and pharynx cancers, as well as 96 cases of larynx cancer. In the fully adjusted model, adherence to the DASH diet was associated with a remarkable 57% reduction in the risk of HNC when comparing extreme quartiles (HR quartile 4 vs 1: 0.43; 95% CI: 0.28, 0.66; P for trend < 0.001). The restricted cubic spline plots demonstrated a linear dose-response relationship between the DASH score and the risk of HNC as well as its subtypes. Subgroup analysis revealed that the protective effect of the DASH diet against HNC was particularly pronounced in individuals with lower daily energy intake. The primary association remained robust in the sensitivity analysis.

CONCLUSIONS

In American middle-aged and older population, adherence to the DASH diet may help prevent HNC, particularly for individuals with lower daily energy intake.

Radiation pneumonitis prediction after stereotactic body radiation therapy based on 3D dose distribution: dosiomics and/or deep learning-based radiomics features

Background

This study was designed to establish radiation pneumonitis (RP) prediction models using dosiomics and/or deep learning-based radiomics (DLR) features based on 3D dose distribution.

Methods

A total of 140 patients with non-small cell lung cancer who received stereotactic body radiation therapy (SBRT) were retrospectively included in this study. These patients were randomly divided into the training (n = 112) and test (n = 28) sets. Besides, 107 dosiomics features were extracted by Pyradiomics, and 1316 DLR features were extracted by ResNet50. Feature visualization was performed based on Spearman’s correlation coefficients, and feature selection was performed based on the least absolute shrinkage and selection operator. Three different models were constructed based on random forest, including (1) a dosiomics model (a model constructed based on dosiomics features), (2) a DLR model (a model constructed based on DLR features), and (3) a hybrid model (a model constructed based on dosiomics and DLR features). Subsequently, the performance of these three models was compared with receiver operating characteristic curves. Finally, these dosiomics and DLR features were analyzed with Spearman’s correlation coefficients.

Results

In the training set, the area under the curve (AUC) of the dosiomics, DLR, and hybrid models was 0.9986, 0.9992, and 0.9993, respectively; the accuracy of these three models was 0.9643, 0.9464, and 0.9642, respectively. In the test set, the AUC of these three models was 0.8462, 0.8750, and 0.9000, respectively; the accuracy of these three models was 0.8214, 0.7857, and 0.8571, respectively. The hybrid model based on dosiomics and DLR features outperformed other two models. Correlation analysis between dosiomics features and DLR features showed weak correlations. The dosiomics features that correlated DLR features with the Spearman’s rho |ρ| ≥ 0.8 were all first-order features.

Conclusion

The hybrid features based on dosiomics and DLR features from 3D dose distribution could improve the performance of RP prediction after SBRT.

[Association of body mass index and waist circumference with frailty among people aged 80 years and older in Chinese]

Objective: To examine the association of body mass index (BMI) and waist circumference (WC) with frailty among oldest-old adults in China. Methods: A total of 7 987 people aged 80 years and older (oldest-old) who participated in the Chinese Longitudinal Healthy Longevity Survey (CLHLS) in 2017-2018 were included. Information on demographic characteristics, behavior pattern, diet, activities of daily living, cognitive function, health status, disease condition were collected by questionnaire and physical examination. Generalized linear mixed model and restricted cubic splines (RCS) were used to analyze the association of BMI and WC with frailty. Results: The mean age of all participants was 91.7 years, and their mean BMI and WC were (21.3±3.5) kg/m² and (82.9±10.5) cm, respectively. The proportion of male was 42.3% (3 377/7 987), and the proportion of people with frailty was 33.7% (2 664/7 987). After controlling confounding factors, compared with T₂ (19.1-22.1 kg/m²) of BMI, the OR (95%CI) of the female T₁ (<19.1 kg/m²) and T₃ (≥22.2 kg/m²) group was 1.39 (1.17-1.65) and 1.27 (1.07-1.52), respectively. Compared with T₂ (77-85 cm) of WC, the OR (95%CI) of female T₁ (<77 cm) and T₃ (≥86 cm) group was 1.20 (1.01-1.42) and 1.10 (0.93-1.31), respectively. The results of multiple linear regression model with restrictive cubic spline showed that there was a non-linear association of BMI and WC with frailty in female. Conclusion: There is a U-shaped association of BMI and WC with frailty in female participants.

[Impact of composite clinical worsening events on outcome of patients with pulmonary arterial hypertension associated with congenital heart disease]

Objective: To explore the impact of composite clinical worsening (cCW) events and its components on the prognosis of patients with pulmonary arterial hypertension associated with congenital heart disease (CHD-PAH). Methods: This is a retrospective study. Patients who were diagnosed with CHD-PAH in Beijing Anzhen Hospital between January 2007 and July 2018, were included, and their baseline clinical data including demographic, clinical manifestations and New York Heart Association (NYHA) classification were collected retrospectively. All-cause deaths and clinical worsening events were recorded, which included syncope, PAH related hospitalization, NYHA classification deterioration and ≥ 2 PAH related clinical symptoms (dyspnea, hemoptysis, edema, chest pain, palpitations, cyanosis) appearance/progress. Three kinds of cCW events were defined: cCW1 (included PAH related hospitalization, NYHA classification deterioration), cCW2 (increased syncope on the basis of cCW1) and cCW3 (increased ≥ 2 PAH related clinical symptoms appearance/progress on the basis of cCW2). The Kaplan-Meier survival curve was used to analyze the long-term survival of the included patients. Univariate and multivariate Cox regression models were used to evaluate the impact of cCW events and their components on the risk of all-cause mortality. Results: A total of 525 patients with CHD-PAH were included in this study. The median age at diagnosis was 20.7 (11.2, 30.3) years. There were 43.8% children (<18 years), and 68.8% female patients. There were 431 patients (82.1%) with NYHA classification II. A total of 180 patients had PAH symptoms at diagnosis. The median follow-up time was 4.5 (2.6, 6.7) years. Forty-seven patients (9.0%) died during the follow-up period. Survival rates at 1, 5 and 10 years after diagnosis of PAH were 98.0%, 89.9% and 84.4%, respectively. Cox multivariate analysis showed that NYHA classification deterioration (HR=3.901, 95%CI 1.863-8.169, P<0.001), ≥2 PAH symptoms appearance/progress (HR=4.458, 95%CI 1.870-10.625, P<0.001), PAH-related hospitalization (HR=4.058, 95%CI 1.851-8.896, P<0.001) and syncope (HR=11.313, 95%CI 4.860-26.332, P<0.001) were independent predictors of increased risk of death. All 3 kinds of cCW events were significantly associated with the significantly increased risk of death, and cCW2 was highly predictive to increased risk of death (HR=15.476, 95%CI 4.346-37.576, P<0.001). Conclusions: The overall long-term prognosis of CHD-PAH patients in this study is relatively good. cCW events and its components (NYHA classification deterioration, ≥2 PAH symptoms occurrence/worsening, PAH-related hospitalization and syncope) have adverse influence on all-cause death in this patient cohort.

Deep Learning for Patient-Specific Quality Assurance: Predicting Gamma Passing Rates for IMRT Based on Delivery Fluence Informed by log Files

Objectives: In this study, we propose a deep learning-based approach to predict Intensity-modulated radiation therapy (IMRT) quality assurance (QA) gamma passing rates using delivery fluence informed by log files. Methods: A total of 112 IMRT plans for chest cancers were planned and measured by portal dosimetry equipped on TrueBeam linac. The convolutional neural network (CNN) based learning model was trained using delivery fluence as inputs and gamma passing rates (GPRs) of 4 different criteria (3%/3 mm, 2%/3 mm, 3%/2 mm, and 2%/2 mm) as outputs. Model performance for both validation and test sets was assessed using mean absolute error (MAE), mean squared error (MSE), root MSE (RMSE), Spearman rank correlation coefficients (Sr), and Determination coefficient (R²) between the measured and predicted GPR values. Results: In the test set, the MAE of the prediction model were 0.402, 0.511, 1.724, and 2.530, the MSE were 0.640, 0.986, 6.654, and 9.508, the RMSE were 0.800, 0.993, 2.580, and 3.083, the Sr were 0.643, 0.684, 0.821, and 0.824 (P < .001) and the R² were 0.4110, 0.4666, 0.6677, and 0.6769 for 3%/3 mm, 3%/2 mm, 2%/3 mm, and 2%/2 mm, respectively. The MAE and RMSE of the prediction model decreased with stricter gamma criteria while the Sr and R² between measured and predicted GPR values increased. Conclusions: The CNN prediction model based on delivery fluence informed by log files could accurately predict IMRT QA passing rates for different gamma criteria. It could reduce QA workload and improve efficiency in pretreatment QA. Our results suggest that the CNN prediction model based on delivery fluence informed by log files may be a promising tool for the gamma evaluation of IMRT QA.

Voxel-Level BED Corrected Dosimetric and Radiobiological Assessment of 2 Kinds of Hybrid Radiotherapy Planning Methods for Stage III NSCLC

Background/purpose: To access the comparative dosimetric and radiobiological advantages of two methods of intensity-modulated radiation therapy (IMRT)-based hybrid radiotherapy planning for stage III nonsmall cell lung cancer (NSCLC). Methods: Two hybrid planning methods were respectively characterized by conventional fraction radiotherapy (CFRT) and stereotactic body radiotherapy (SBRT) and CFRT and simultaneous integrated boost (SIB) planning. All plans were retrospectively completed using the 2 methods for 20 patients with stage III NSCLC. CFRT and SBRT dose regimes 2 Gy  ×  30 f and 12.5 Gy  ×  4 f were, respectively, used for planning target volume of lymph node (PTV_LN) and planning target volume of the primary tumor (PTV_PT), while dose regimes 2 Gy  ×  26 f for PTV_LN and sequential 2 Gy  ×  4 f for PTV_LN combined with 12.5 Gy  ×  4 f for PTV_PT were adopted for CFRT and SIB plans. SBRT and SIB EQD₂ dose were calculated voxel by voxel, and then, respectively, superimposed with 30-fraction and 26-fraction CFRT plan dose to achieve biological equivalent dose (BED) dosimetric parameters of CFRT and SBRT and CFRT and SIB plans. Tumor control probability (TCP)/normal tissue complication probability (NTCP) was, respectively, calculated by equivalent uniform dose/Lyman–Kutcher–Burman models. BED plan parameters and TCP/NTCP were analyzed between 2 methods of hybrid planning. Primary tumor/lymph node (LN)/total TCP values were, respectively, evaluated as a function of the radiation dose needed to control 50% of tumor (TCD₅₀) for 20 patients. Dosimetric errors were analyzed by nontransit electronic portal imaging device dosimetry measurement during hybrid plan delivery. Results: Statistically lower BED plan parameters of PTV_LN D₂ and homogeneity index resulted in slightly lower averaged LN/total TCP curves by CFRT and SIB planning. The gaps between Max and Min LN/total TCP curves were significantly closer for CFRT and SIB planning, which indicated better robustness of LN/total TCPs. A lower esophagus dose resulted in a lower esophagus NTCP by CFRT and SIB planning, which may be compromised by 1 week shorter overall treatment time by CFRT and SIB irradiation. Spinal cord D_max was significantly reduced by CFRT and SIB plans. The dose verification results of the subplans involved in hybrid plans were acceptable, which showed that the 2 methods of hybrid planning could be delivered accurately in our center. Conclusion: CFRT and SIB plannings have more advantages on BED plan parameters and TCP/NTCP than CFRT and SBRT planning, and both methods of IMRT-based hybrid planning could be executed accurately for stage III NSCLC. The effectiveness of the results needs to be validated in the hybrid trial.

Account for the Full Extent of Esophagus Motion in Radiation Therapy Planning: A Preliminary Study of the IRV of the Esophagus

Objective

Accounting for esophagus motion in radiotherapy planning is an important basis for accurate assessment of toxicity. In this study, we calculated how much the delineations of the esophagus should be expanded based on three-dimensional (3D) computed tomography (CT), four-dimensional (4D) average projection (AVG), and maximum intensity projection (MIP) scans to account for the full extent of esophagus motion during 4D imaging acquisition.

Methods and Materials

The 3D and 4D CT scans of 20 lung cancer patients treated with conventional radiotherapy and 20 patients treated with stereotactic ablative radiation therapy (SBRT) were used. Radiation oncologists contoured the esophagus on the 3DCT, AVG, MIP and 25% exhale scans, and the combination of the esophagus in every phase of 4DCT. The union of all 4D phase delineations (U4D) represented the full extent of esophagus motion during imaging acquisition. Surface distances from U4D to 3D, AVG, and MIP volumes were calculated. Distances in the most extreme surface points (1.5 cm most superoinferior, 10% most right/left/anteroposterior) were used to derive margins accounting only for systematic (delineation) errors.

Results

Esophagus delineations on the MIP were the closest to the full extent of motion, requiring only 6.9 mm margins. Delineations on the AVG and 3D scans required margins up to 7.97 and 7.90 mm, respectively. The largest margins were for the inferior, right, and anterior aspects for the delineations on the 3D, AVG, and MIP scans, respectively.

Conclusion

Delineations on 3D, AVG, or MIP scans required extensions for representing the esophagus’s full extent of motion, with the MIP requiring the smallest margins. Research including daily imaging to determine the random components for the margins and dosimetric measurements to determine the relevance of creating a planning organ at risk volume (PRV) of the esophagus is required.

How Does the Gradient Measure of the Lung SBRT Treatment Plan Depend on the Tumor Volume and Shape?

Purpose

Gradient measure (GM) is a critical index related to normal tissue sparing in radiosurgery. This study aims to describe the dependence of GM on target volume and target shape for lung stereotactic body radiation therapy (SBRT) treatment plans.

Methods

A total of 307 peripheral and 119 central lung SBRT treatment plans were enrolled for this study. A least-squares regression was used for data analysis. First, the equations with different functional forms were established to determine the dependence of GM on a univariaty (V_P or Sp) and bivariaty (V_P and Sp), respectively. Then, the correlation coefficients and p-values of variables for all equations were compared and analyzed to determine the dependence of GM on PTV volume (VP) and sphericity (Sp).

Results

The power equations had the highest coefficient of determination (R²) in the dependence results of GM on univariate V_P. The equations were GM = 0.674 V P 0.178 and GM = 0.660 V P 0.185 for peripheral and central lesions, respectively. On the other hand, the R² of all functional forms were less than 0.25 when the relationship of GM versus univariate Sp was analyzed. Similarly, the power equation also obtained the highest R² in bivariaty V_P and Sp analysis, whether for central or peripheral. However, the R² of the bivariate equations were not improved compared with those of univariate equations. Moreover, the p-values of the variable Sp were greater than 0.05.

Conclusions

The GM of the lung SBRT plan is shape-independent and volume-dependent. The dependence of GM on PTV volume for peripheral and central lung cancer can be described by two different power equations. The results of this study can be used as a potential tool to assist dosimetric quality control during the radiosurgery process.

Influence of Clinical and Tumor Factors on Interfraction Setup Errors With Rotation Correction for Vacuum Cushion in Lung Stereotactic Body Radiation Therapy

Purpose

To explore the influence of clinical and tumor factors over interfraction setup errors with rotation correction for non-small cell lung cancer (NSCLC) stereotactic body radiation therapy (SBRT) patients immobilized in vacuum cushion (VC) to better understand whether patient re-setup could further be optimized with these parameters.

Materials and Methods

This retrospective study was conducted on 142 NSCLC patients treated with SBRT between November 2017 to July 2019 in the local institute. Translation and rotation setup errors were analyzed in 732 cone-beam computed tomography (CBCT) scans before treatment. Differences between groups were analyzed using independent sample t-test. Logistic regression test was used to analyze possible correlations between patient re-setup and clinical and tumor factors.

Results

Mean setup errors were the largest in anterior-posterior (AP) direction (3.2 ± 2.4 mm) compared with superior-inferior (SI) (2.8 ± 2.1 mm) and left-right (LR) (2.5 ± 2.0 mm) directions. The mean values were similar in pitch, roll, and rtn directions. Of the fractions, 83.7%, 90.3%, and 86.6% satisfied setup error tolerance limits in AP, SI, and LR directions, whereas 95% had rotation setup errors of <2° in the pitch, roll, or rtn directions. Setup errors were significantly different in the LR direction when age, body mass index (BMI), and "right vs. left" location parameters were divided into groups. Both univariate and multivariable model analyses showed that age (p = 0.006) and BMI (p = 0.002) were associated with patient re-setup.

Conclusions

Age and BMI, as clinical factors, significantly influenced patient re-setup in the current study, whereas all other clinical and tumor factors were not correlated with patient re-setup. The current study recommends that more attention be paid to setup for elderly patients and patients with larger BMI when immobilized using VC, especially in the left-right direction.

Investigation of Predictors to Achieve Acceptable Lung Dose in T-Shaped Upper and Middle Esophageal Cancer With IMRT and VMAT

Purpose

The purpose of this study is to investigate whether there are predictors and cutoff points that can predict the acceptable lung dose using intensity-modulated radiation therapy (IMRT) and volume-modulated arc therapy (VMAT) in radiotherapy for upper ang middle esophageal cancer.

Material and Methods

Eighty-two patients with T-shaped upper-middle esophageal cancer (UMEC) were enrolled in this retrospective study. Jaw-tracking IMRT plan (JT-IMRT), full-arc VMAT plan (F-VMAT), and pactial-arc VMAT plan (P-VMAT) were generated for each patient. Dosimetric parameters such as MLD and V20 of total lung were compared among the three plannings. Ten factors such as PCTV_inferior length and PCTV_inferior length/total lung length were calculated to find the predictors and cutoff points of the predictors. All patients were divided into two groups according to the cutoff points, and the dosimetric differences between the two groups of the three plans were compared. ANOVA, receiver operating characteristic (ROC) analysis, and Mann–Whitney U-test were performed for comparisons between datasets. A p <0.05 was considered statistically significant.

Result

The quality of the targets of the three plannings was comparable. The total lung dose in P-VMAT was significantly lower than that in JT IMRT and F-VMAT. Monitor unit (MU) of F-VMAT and P-VMAT was significantly lower than that of JT IMRT. ROC analysis showed that among JT IMRT, F-VMAT, and P-VMAT, PCTV_i-L, and PCTV_i-L/TL_L had diagnostic power to predict the suitability of RT plans according to lung dose constraints of our department. For JT IMRT, the cutoff points of PCTV_i-L and PCTV_i-L/TL_L were 16.6 and 0.59. For F-VMAT, the cutoff points of PCTV_i-L and PCTV_i-L/TL_L were 16.75 and 0.62. For P-VMAT, the cutoff points of PCTV_i-L and PCTV_i-L/TL_L were 15.15 and 0.59. After Mann–Whitney U-test analysis, it was found that among the three plannings, the group with lower PCTV_i-L and PCTV_i-L/TL_L could significantly reduce the dose of total lung and heart (p <0.05).

Conclusion

PCTV_{i-L <}16.6 and PCTV_i-L/TL_L <0.59 for JT IMRT, PCTV_{i-L <}16.75 and PCTV_i-L/TL_L <0.62 for F-VMAT and PCTV_{i-L <}15.15, and PCTV_i-L/TL_L <0.59 for P-VMAT can predict whether patients with T-shaped UMEC can meet the lung dose limits of our department.

Functional outcomes after submandibular gland-sparing neck dissection in patients with oral squamous cell carcinoma

The feasibility of submandibular gland (SMG) preservation during neck dissection has been described. The aim of this study was to analyse the functional outcomes in patients undergoing SMG preservation during neck dissection for cT1-2N0 oral squamous cell carcinoma. Consecutive patients were divided into two groups based on the management of the SMG, and underwent a saliva flow test before surgery, 7 days after surgery, and at 3, 6, 9, and 12 months after surgery. All enrolled patients completed the fourth version of the University of Washington Quality of Life (UWQOL) questionnaire at 12 months after surgery. In patients who underwent SMG preservation during neck dissection, the flow rate at 7 days after surgery was significantly lower than that preoperative; however, it gradually returned to baseline at 9 months after surgery. The saliva flow rate at 9 months after surgery was similar to that at 12 months after surgery. Further, patients with SMG preservation had higher scores for the activity, swallowing, chewing, and saliva domains than patients without SMG preservation. The results of the study suggest that saliva secretion ability can be preserved following SMG-sparing neck dissection, and that SMG preservation improves postoperative quality of life.

First-phase ejection fraction, a measure of pre-clinical heart failure, is strongly associated with increased mortality in patients with COVID-19

Introduction

Presence of heart failure is associated with a poor prognosis in patients with COVID-19. The aim of the present study was to examine whether first-phase ejection fraction (EF1), the ejection fraction measured in early systole up to the time of peak aortic velocity, a sensitive measure of pre-clinical heart failure, is associated with survival in patients hospitalised with COVID-19.

Methods

A retrospective outcome study was performed in patients hospitalised with COVID-19 who underwent echocardiography (n=380) at the West Branch of the Union Hospital, Wuhan, China and in patients admitted to King's Health Partners in South London UK. Association of EF1 with survival was performed using Cox proportional hazards regression. EF1 was compared in patients with COVID-19 and in historical controls with similar co-morbidities (n=266) who had undergone echocardiography before the COVID-19 pandemic.

Results

In patients with COVID-19, EF1 was a strong predictor of survival in each patient group (Wuhan and London). In the combined group, EF1 was a stronger predictor of survival than other clinical, laboratory and echocardiographic characteristics including age, co-morbidities and biochemical markers (figure 1). A cut-off value of 25% for EF1 gave a hazard ratio of 5.23 (95% CI: 2.85–9.60, p<0.001) unadjusted and 4.83 (95% CI: 2.35–9.95, p<0.001) when adjusted for demographics, co-morbidities, hs-cTnI and CRP (figure 2). EF1 was similar in patients with and without COVID-19 (23.2±7.3 vs 22.0±7.6%, p=0.092, adjusted for prevalence of risk factors and co-morbidities).

Conclusion

Impaired first-phase ejection fraction is strongly associated with mortality in COVID-19 and probably reflects pre-existing, pre-clinical heart failure.

Funding Acknowledgement

Type of funding sources: Public grant(s) – National budget only. Main funding source(s): National Institute for Health Research (NIHR) UKBritish Heart Foundation (BHF) UK

Automatic segmentation of lung tumors on CT images based on a 2D & 3D hybrid convolutional neural network

OBJECTIVE

A stable and accurate automatic tumor delineation method has been developed to facilitate the intelligent design of lung cancer radiotherapy process. The purpose of this paper is to introduce an automatic tumor segmentation network for lung cancer on CT images based on deep learning.

METHODS

In this paper, a hybrid convolution neural network (CNN) combining 2D CNN and 3D CNN was implemented for the automatic lung tumor delineation using CT images. 3D CNN used V-Net model for the extraction of tumor context information from CT sequence images. 2D CNN used an encoder-decoder structure based on dense connection scheme, which could expand information flow and promote feature propagation. Next, 2D features and 3D features were fused through a hybrid module. Meanwhile, the hybrid CNN was compared with the individual 3D CNN and 2D CNN, and three evaluation metrics, Dice, Jaccard and Hausdorff distance (HD), were used for quantitative evaluation. The relationship between the segmentation performance of hybrid network and the GTV volume size was also explored.

RESULTS

The newly introduced hybrid CNN was trained and tested on a dataset of 260 cases, and could achieve a median value of 0.73, with mean and stand deviation of 0.72 ± 0.10 for the Dice metric, 0.58 ± 0.13 and 21.73 ± 13.30 mm for the Jaccard and HD metrics, respectively. The hybrid network significantly outperformed the individual 3D CNN and 2D CNN in the three examined evaluation metrics (p < 0.001). A larger GTV present a higher value for the Dice metric, but its delineation at the tumor boundary is unstable.

CONCLUSIONS

The implemented hybrid CNN was able to achieve good lung tumor segmentation performance on CT images.

ADVANCES IN KNOWLEDGE

The hybrid CNN has valuable prospect with the ability to segment lung tumor.

A 2D-3D hybrid convolutional neural network for lung lobe auto-segmentation on standard slice thickness computed tomography of patients receiving radiotherapy

Background

Accurate segmentation of lung lobe on routine computed tomography (CT) images of locally advanced stage lung cancer patients undergoing radiotherapy can help radiation oncologists to implement lobar-level treatment planning, dose assessment and efficacy prediction. We aim to establish a novel 2D–3D hybrid convolutional neural network (CNN) to provide reliable lung lobe auto-segmentation results in the clinical setting.

Methods

We retrospectively collected and evaluated thorax CT scans of 105 locally advanced non-small-cell lung cancer (NSCLC) patients treated at our institution from June 2019 to August 2020. The CT images were acquired with 5 mm slice thickness. Two CNNs were used for lung lobe segmentation, a 3D CNN for extracting 3D contextual information and a 2D CNN for extracting texture information. Contouring quality was evaluated using six quantitative metrics and visual evaluation was performed to assess the clinical acceptability.

Results

For the 35 cases in the test group, Dice Similarity Coefficient (DSC) of all lung lobes contours exceeded 0.75, which met the pass criteria of the segmentation result. Our model achieved high performances with DSC as high as 0.9579, 0.9479, 0.9507, 0.9484, and 0.9003 for left upper lobe (LUL), left lower lobe (LLL), right upper lobe (RUL), right lower lobe (RLL), and right middle lobe (RML), respectively. The proposed model resulted in accuracy, sensitivity, and specificity of 99.57, 98.23, 99.65 for LUL; 99.6, 96.14, 99.76 for LLL; 99.67, 96.13, 99.81 for RUL; 99.72, 92.38, 99.83 for RML; 99.58, 96.03, 99.78 for RLL, respectively. Clinician's visual assessment showed that 164/175 lobe contours met the requirements for clinical use, only 11 contours need manual correction.

Conclusions

Our 2D–3D hybrid CNN model achieved accurate automatic segmentation of lung lobes on conventional slice-thickness CT of locally advanced lung cancer patients, and has good clinical practicability.

Dosimetric Comparison, Treatment Efficiency Estimation, and Biological Evaluation of Popular Stereotactic Radiosurgery Options in Treating Single Small Brain Metastasis

Objectives

This study aimed to show the advantages of each stereotactic radiosurgery (SRS) treatment option for single small brain metastasis among Gamma Knife (GK), Cone-based VMAT (Cone-VMAT), and MLC-based CRT (MLC-CRT) plans.

Materials and Methods

GK, Cone-VMAT, and MLC-CRT SRS plans were retrospectively generated for 11 patients with single small brain metastasis whose volume of gross tumor volume (GTV) ranged from 0.18 to 0.76 cc (median volume 0.60 cc). Dosimetric parameters, treatment efficiency, and biological parameters of the three techniques were compared and evaluated. The metric variation with the planning target volume (PTV) was also studied.

Results

The conformity index (CI) was similar in GK and MLC-CRT plans, higher than Cone-VMAT. Cone-VMAT achieved comparable volume covered by 12 Gy (V12) and gradient index (GI) as GK, lower than MLC-CRT. The heterogeneity index (HI) of GK, Cone-VMAT, and MLC-CRT decreased sequentially. GK gave the lowest volume covered by 3 Gy (V3) and 6 Gy (V6), while MLC-CRT got the highest. The beam-on time and treatment time of GK, Cone-VMAT, and MLC-CRT decreased in turn. Tumor control probability (TCP) of all three SRS plans was greater than 98%, and normal tissue complication probability (NTCP) of all organs at risk (OARs) was below 0.01%. GK and Cone-VMAT resulted in superior TCP and NTCP of the normal brain tissue than MLC-CRT. The relative value of Cone-VMAT and GK for all metrics hardly changed with the target volume. Except for the unchanged HI and TCP, the other results of MLC-CRT with respect to GK improved as the target volume increased. MLC-CRT could produce higher CI than GK and Cone-VMAT when the target volume increased above 2 and 1.44 cc, respectively.

Conclusion

For single small brain metastases, Cone-VMAT may be used as an alternative to GK-free centers. In addition to the advantage of short treatment time, MLC-CRT showed superiority in CI as the target volume increased. Treatment centers can choose appropriate SRS technique on a case-by-case basis according to institutional conditions and patients’ individual needs.

Dosimetric effect of intensity-modulated radiation therapy for postoperative non-small cell lung cancer with and without air cavity in the planning target volume

To evaluate the dosimetric effect of intensity-modulated radiation therapy (IMRT) for postoperative non-small cell lung cancer (NSCLC), with and without the air cavity in the planning target volume (PTV). Two kinds of IMRT plans were made for 21 postoperative NSCLC patients. In Plan-0: PTV included the tracheal air cavity, and in Plan-1: the air cavity was subtracted from the PTV. For PTV, the dosimetric parameters, including D_mean, D₉₈, D₉₅, D₂, D_0.2, conformity index (CI), and homogeneity index (HI) were evaluated. For organs at risk (OARs), the evaluation indexes, included the V₅, V₂₀ and the mean lung dose (MLD) of total lung, the V₃₀, V₄₀, and the mean heart dose (MHD) of heart, the spinal cord D_max, and the V₃₅ and the mean esophageal dose (MED) of esophagus. The number of segments and MUs were also recorded. Additionally, the correlation between the Plan-1 dosimetric change value relative to Plan-0, the size of air cavity, and the volume proportion of the cavity in the PTV was also analyzed. The D_mean of PTV, D₂, D_0.2, HI and CI in Plan-1 decreased compared with those in Plan-0. For OARs, the V₃₀, MHD, and MED also decreased. The CI change value of Plan-1 relative to Plan-0 had a significantly negative correlation with the size and the volume proportion of air cavity, and the MED change value also had a significantly negative correlation with the air cavity size. The IMRT plans for patients with postoperative NSCLC can achieve a better target dose distribution and offer an improved sparing of the heart and esophagus by removing the PTV air cavity, while reducing the target conformity. The change value of CI and MED had a significantly negative correlation with the air cavity size.

Dosimetric comparison and biological evaluation of fixed-jaw intensity-modulated radiation therapy for T-shaped esophageal cancer

Background

To evaluate the dosimetric and biological benefits of the fixed-jaw (FJ) intensity-modulated radiation therapy (IMRT) technique for patients with T-shaped esophageal cancer.

Methods

FJ IMRT plans were generated for thirty-five patients and compared with jaw tracking (JT) IMRT, static jaw (SJ) IMRT and JT volumetric modulated arc therapy (VMAT). Dosimetric parameters, tumor control probability (TCP) and normal tissue complication probability (NTCP), monitor units (MUs), delivery time and gamma passing rate, as a measure of dosimetric verification, were compared. The correlation between the length of PTV-C below the upper boundary of lung tissue (PTV-C_inferior) and dosimetric parameters and NTCP of the lung tissue were analyzed.

Results

The homogeneity and conformity of the target in the four plans were basically equivalent. When compared to the JT IMRT and SJ IMRT plans, FJ IMRT plan led to a statistically significant improvement in the NTCP and low-middle dosimetric parameters of the lung, and the improvement had a moderately positive correlation with the length of PTV-C_inferior, with a correlation coefficient ranging from 0.523 to 0.797; the FJ IMRT plan exhibited better lung sparing in low-dose volumes than the JT VMAT plan. The FJ IMRT plan had similar MUs (888 ± 99) and delivery times (516.1 ± 54.7 s) as the JT IMRT plan (937 ± 194, 522 ± 5.6 s) but higher than SJ IMRT (713 ± 137, 488.8 ± 45.2 s) and JT VMAT plan (517 ± 59, 263.7 ± 43.3 s).

Conclusions

The FJ IMRT technique is superior in reducing the low-dose volumes of lung tissues for patients with T-shaped esophageal cancer.

Geometric and Dosimetric Changes in Tumor and Lung Tissue During Radiotherapy for Lung Cancer With Atelectasis

Background and Purpose

This article retrospectively characterized the geometric and dosimetric changes in target and normal tissues during radiotherapy for lung cancer patients with atelectasis.

Materials and Methods

A total of 270 cone beam computed tomography (CBCT) scans of 18 lung patients with atelectasis were collected. The degree and time of resolution or expansion of the atelectasis were recorded. The geometric, dosimetric, and biological changes in the target and lung tissue were also quantified.

Results

There were two patients with expansion, four patients with complete regression, six patients with partial regression, and six patients with no change. The time of resolution or expansion varied. The tumor volume increased by 3.8% in the first seven fractions, then decreased from the 9th fraction, and by 33.4% at the last CBCT. In the LR direction, the average center of mass (COM), boundaries of the tumors gradually shifted mediastinally. In the AP direction, the COM of the tumors was shifted slightly in the posterior direction and then gradually shifted to the anterior direction; the boundaries of the tumors all moved mediastinally. In the SI direction, the COM of the tumors on the right side of the body was substantially shifted toward the head direction. The boundaries of the tumors varied greatly. D₂, D₉₈, D_mean, V₉₅, V₁₀₇, and TCP of the PTV were reduced during radiotherapy and were reduced to their lowest values during the last two fractions. The volume of the ipsilateral lung tended to increase gradually. The V₅, V₁₀, V₂₀, V₃₀, V₄₀, and NTCP of the total lung gradually increased with the fraction.

Conclusions

For most patients, regression of the atelectasis occurred, and the volume of the ipsilateral lung tended to increase while the tumor volume decreased, and the COM and boundary of the tumors shifted toward mediastinum, which caused an insufficient dose to the target and an overdose to the lungs. Regression or expansion may occur for any fraction, and it is therefore recommended that CBCT be performed at least every other day.

Virtual Patient-Specific Quality Assurance of IMRT Using UNet++: Classification, Gamma Passing Rates Prediction, and Dose Difference Prediction

The dose verification in radiotherapy quality assurance (QA) is time-consuming and places a heavy workload on medical physicists. To provide a clinical tool to perform patient specific QA accurately, the UNet++ is investigated to classify failed or pass fields (the GPR lower than 85% is considered “failed” while the GPR higher than 85% is considered “pass”), predict gamma passing rates (GPR) for different gamma criteria, and predict dose difference from virtual patient-specific quality assurance in radiotherapy. UNet++ was trained and validated with 473 fields and tested with 95 fields. All plans used Portal Dosimetry for dose verification pre-treatment. Planar dose distribution of each field was used as the input for UNet++, with QA classification results, gamma passing rates of different gamma criteria, and dose difference were used as the output. In the test set, the accuracy of the classification model was 95.79%. The mean absolute error (MAE) were 0.82, 0.88, 2.11, 2.52, and the root mean squared error (RMSE) were 1.38, 1.57, 3.33, 3.72 for 3%/3mm, 3%/2 mm, 2%/3 mm, 2%/2 mm, respectively. The trend and position of the predicted dose difference were consistent with the measured dose difference. In conclusion, the Virtual QA based on UNet++ can be used to classify the field passed or not, predict gamma pass rate for different gamma criteria, and predict dose difference. The results show that UNet++ based Virtual QA is promising in quality assurance for radiotherapy.

Regression models for predicting physical and EQD2 plan parameters of two methods of hybrid planning for stage III NSCLC

Background/purpose

To establish regression models of physical and equivalent dose in 2 Gy per fraction (EQD₂) plan parameters of two kinds of hybrid planning for stage III NSCLC.

Methods

Two kinds of hybrid plans named conventional fraction radiotherapy & stereotactic body radiotherapy (C&S) and conventional fraction radiotherapy & simultaneous integrated boost (C&SIB) were retrospectively made for 20 patients with stage III NSCLC. Prescription dose of C&S plans was 2 Gy × 30f for planning target volume of lymph node (PTV_LN) and 12.5 Gy × 4f for planning target volume of primary tumor (PTV_PT), while prescription dose of C&SIB plans was 2 Gy × 26f for PTV_LN and sequential 2 Gy × 4f for PTV_LN combined with 12.5 Gy × 4f for PTV_PT. Regression models of physical and EQD₂ plan parameters were established based on anatomical geometry features for two kinds of hybrid plans. The features were mainly characterized by volume ratio, min distance and overlapping slices thickness of two structures. The possibilities of regression models of EQD₂ plan parameters were verified by spearman’s correlation coefficients between physical and EQD₂ plan parameters, and the influence on the consistence of fitting goodness between physical and EQD₂ models was investigated by the correlations between physical and EQD₂ plan parameters. Finally, physical and EQD₂ models predictions were compared with plan parameters for two new patients.

Results

Physical and EQD₂ plan parameters of PTV_LN CI_60Gy have shown strong positive correlations with PTV_LN volume and min distance_{(PT to LN)}, and strong negative correlations with PTV_PT volume for two kinds of hybrid plans. PTV_(PT+LN) CI_60Gy is not only correlated with above three geometry features, but also negatively correlated with overlapping slices thickness_{(PT and LN)}. When neck lymph node metastasis was excluded from PTV_LN volume, physical and EQD₂ total lung V₂₀ showed a high linear correlation with corrected volume ratio_{(LN to total lung).} Meanwhile, physical total lung mean dose (MLD) had a high linear correlation with corrected volume ratio_{(LN to total lung)}, while EQD₂ total lung MLD was not only affected by corrected volume ratio_{(LN to total lung)} but also volume ratio_{(PT to total lung).} Heart D₅, D₃₀ and mean dose (MHD) would be more susceptible to overlapping structure_{(heart and LN)}. Min distance_{(PT to ESO)} may be an important feature for predicting EQD₂ esophageal max dose for hybrid plans. It’s feasible for regression models of EQD₂ plan parameters, and the consistence of the fitting goodness of physical and EQD₂ models had a positive correlation with spearman’s correlation coefficients between physical and EQD₂ plan parameters. For total lung V₂₀, ipsilateral lung V₂₀, and ipsilateral lung MLD, the models could predict that C&SIB plans were higher than C&S plans for two new patients.

Conclusion

The regression models of physical and EQD₂ plan parameters were established with at least moderate fitting goodness in this work, and the models have a potential to predict physical and EQD₂ plan parameters for two kinds of hybrid planning.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13014-021-01848-9.

A novel specific grading standard study of auto-segmentation of organs at risk in thorax: subjective-objective-combined grading standard

BACKGROUND

To develop a novel subjective-objective-combined (SOC) grading standard for auto-segmentation for each organ at risk (OAR) in the thorax.

METHODS

A radiation oncologist manually delineated 13 thoracic OARs from computed tomography (CT) images of 40 patients. OAR auto-segmentation accuracy was graded by five geometric objective indexes, including the Dice similarity coefficient (DSC), the difference of the Euclidean distance between centers of mass (ΔCMD), the difference of volume (ΔV), maximum Hausdorff distance (MHD), and average Hausdorff distance (AHD). The grading results were compared with those of the corresponding geometric indexes obtained by geometric objective methods in the other two centers. OAR auto-segmentation accuracy was also graded by our subjective evaluation standard. These grading results were compared with those of DSC. Based on the subjective evaluation standard and the five geometric indexes, the correspondence between the subjective evaluation level and the geometric index range was established for each OAR.

RESULTS

For ΔCMD, ΔV, and MHD, the grading results of the geometric objective evaluation methods at our center and the other two centers were inconsistent. For DSC and AHD, the grading results of three centers were consistent. Seven OARs' grading results in the subjective evaluation standard were inconsistent with those of DSC. Six OARs' grading results in the subjective evaluation standard were consistent with those of DSC. Finally, we proposed a new evaluation method that combined the subjective evaluation level of those OARs with the range of corresponding DSC to determine the grading standard. If the DSC ranges between the adjacent levels did not overlap, the DSC range was used as the grading standard. Otherwise, the mean value of DSC was used as the grading standard.

CONCLUSIONS

A novel OAR-specific SOC grading standard in thorax was developed. The SOC grading standard provides a possible alternative for evaluation of the auto-segmentation accuracy for thoracic OARs.

[A cohort study on plasma uric acid levels and the risk of type 2 diabetes mellitus among the oldest old in longevity areas of China]

Objective: To investigate the effect of plasma uric acid level on the incident risk of type 2 diabetes mellitus (T2DM) among the oldest old (those aged ≥80 years). Methods: Participants were recruited from the Healthy Aging and Biomarkers Cohort Study (HABCS), which conducted a baseline survey in 2008-2009 and follow-up of 3 times in 2011-2012, 2014, and 2017-2018, respectively. A total of 2 213 oldest old were enrolled in this study. The general demographic, socioeconomic, lifestyle and disease data of the oldest old were collected, and physical measurements were made for the oldest old. Fasting venous blood was collected for uric acid and blood glucose detection. Information on the incident and death of T2DM were collected through the follow-up. Cox proportional hazard regression model was used to explore the association of hyperuricemia and plasma uric acid level with the incidence of T2DM. Restricted cubic spline (RCS) function was used to explore the dose-response relationship of plasma uric acid levels with the risk of T2DM. Results: The age of participants was (93.2±7.6) years old, and 66.7% of the participants (1 475) were female. The plasma uric acid level at baseline was (289.1±88.0)μmol/L, and the prevalence of hyperuricemia was 13.3% (294 cases). During 9 years of cumulative follow-up of 7 471 person-years (average of 3.38 years for each), 122 new cases of T2DM occurred and the incidence density was 1 632.98/10⁵ person year. Cox proportional hazards regression analysis showed that per 10μmol/L increase in plasma uric acid level, the risk of T2DM increased by 1.1% [HR (95%CI): 1.011 (1.004, 1.017)]. Compared with the participants with the lowest quintile of plasma uric acid (Q₁), the risk of diabetes increased by 20.7 % among the oldest old with uric acid in the highest quintile (Q₅) [HR (95%CI):1.207 (1.029, 1.416)]. The risk of T2DM was 19.2% higher in the hyperuricemia group than that in the oldest old with normal plasma uric acid [HR (95%CI): 1.192 (1.033, 1.377)]. RCS function showed that the risk of T2DM increased with the increase in plasma uric acid levels in a nonlinear dose-response relationship (P=0.016). Conclusion: The incident risk of T2DM increases with the elevates of plasma uric acid levels in the oldest old.

Prognosis of patients in prolonged coma after severe carbon monoxide poisoning

BACKGROUND

Late recovery in patients following prolonged coma from carbon monoxide poisoning have been reported, but the probability is unclear. The purpose of this research was to assess the prognosis of patients in prolonged coma after severe carbon monoxide poisoning and related clinical and imaging features.

METHODS

There were 13 patients who had been in a state of coma for >7 days after acute carbon monoxide poisoning in the retrospective observational study, and demographic data, clinical data, laboratory data, complications, and image data were collected. Outcome was assessed by means of the Glasgow outcome scale after 1 year. The relationship between complications and imaging manifestations and prognosis was also analyzed.

RESULTS

One year after severe carbon monoxide poisoning, two patients (15.4%) had died (GOS 1), nine (69.2%) were in a persistent vegetative state (GOS 2), one (7.7%) was moderately disabled (GOS 4), and one (7.7%) achieved a good recovery (GOS 5) with minimal disability.

CONCLUSIONS

Most patients with prolonged coma after severe carbon monoxide poisoning had a poor prognosis, although the younger patients had a better prognosis. Respiratory failure, hypotension and renal failure during the course of the disease were associated with a poor prognosis. The prognosis of patients with injuries in two sites in early CT was poor. Multiple lesions (≥3) and extensive white matter damage (Fazekas grade (PVH or DWMH) = 3) on MRI of chronic phase were also associated with a poor prognosis.

[Association of blood lead level with cognition impairment among elderly aged 65 years and older in 9 longevity areas of China]

Objective: To investigate the association between blood lead concentrations and cognition impairment among Chinese older adults aged 65 or over. Method: Data was collected in 9 longevity areas from Heathy Aging and Biomarkers Cohort Study between 2017 and 2018. This study included 1 684 elderly aged 65 years and older. Information about demographic characteristics, socioeconomic factors, health status and cognitive function score of respondents were collected by questionnaire survey and physical examination. Venous blood of the subjects was collected to detect the blood lead concentration. Subjects were stratified into four groups (Q₁-Q₄) by quartile of blood lead concentration. Multivariate logistic regression model was used to analyze the association between blood lead concentration and cognitive impairment. The linear or non-linear association between blood lead concentration and cognitive impairment were described by restrictive cubic splines (RCS). Results: Among the 1 684 respondents, 843 (50.1%) were female and 191 (11.3%) suffered from cognition impairment. After adjusting for confounding factors, the OR value and 95%CI of cognition impairment was 1.05 (1.01-1.10) for every 10 μg/L increase in blood lead concentration in elderly; Compared with the elderly in Q₁, the elderly with higher blood lead concentration had an increased risk of cognitive impairment. The OR value and 95%CI of Q₂, Q₃ and Q₄ groups were 1.19 (0.69-2.05), 1.45 (0.84-2.51) and 1.92 (1.13-3.27), respectively. Conclusion: Higher blood lead concentration is associated with cognitive impairment among the elderly aged 65 years and older in 9 longevity areas in China.

[Association of blood oxidative stress level with hypertriglyceridemia in the elderly aged 65 years and older in 9 longevity areas of China]

Objective: To investigate the association of blood oxidative stress level with hypertriglyceridemia in the elderly aged 65 years and older in China. Methods: A total of 2 393 participants aged 65 years and older were recruited in 9 longevity areas from Heathy Aging and Biomarkers Cohort Study, during 2017 to 2018. Information on demographics characteristic, life style and health status were collected by questionnaire and physical examination, and venous blood was collected to detect the levels of blood oxidative stress and hypertriglyceridemia. The linear or non-linear association between oxidative stress and hypertriglyceridemia was described by restrictive cubic splines (RCS) fitting multiple linear regression model. The generalized linear mixed effect model was conducted to assess the association between oxidative stress and hypertriglyceridemia. Results: A total of 2 393 participants, mean age was 84.6 years, the youngest was 65 and the oldest was 112, the male was 47.9%(1 145/2 393), the triglyceride level was (1.4±0.8) mmol/L. The hypertriglyceridemia detection rate was 9.99%(239/2 393). The results of multiple linear regression model with restrictive cubic spline fitting showed that MDA level was linear association with triglyceride level; SOD level was nonlinear association with triglyceride level. MDA level had significantly association with hypertriglyceridemia, and the corresponding OR value was 1.063 (95%CI: 1.046,1.081) with 1 nmol/ml increment of blood MDA; SOD level had significantly association with hypertriglyceridemia, and the corresponding OR value was 0.986(95%CI: 0.983,0.989) with 1 U/ml increment of blood SOD. Conclusion: Among the elderly aged 65 and older in 9 longevity areas in China, MDA and SOD levels were associated with the risk of hypertriglyceridemia.

First-phase ejection fraction predicts response to cardiac resynchronization therapy and adverse outcomes

Background

Cardiac Resynchronization Therapy (CRT) is an important therapeutic treatment for chronic heart failure. However, even in carefully selected cases up to 40% of patient fail to respond. First-phase ejection fraction (EF1), the ejection fraction up to the time of maximal ventricular contraction, is a novel and more sensitive echocardiographic measure of early systolic function.

Purpose

We examined the value of EF1, to predict response to CRT and clinical outcomes after CRT.

Methods

We analysed echocardiograms from 197 patients who underwent CRT between 2009 and 2018 and were followed to determine clinical outcomes.

Results

Volumetric response rate (reduction in end-systolic volume ≥15%) was 92.3% vs. 12.1%, for those with EF1 in the highest vs. lowest tertiles (p&lt;0.001). A cut-off value of 11.9% for EF1 had &gt;85% sensitivity and specificity for prediction of response to CRT; on multivariate logistic regression analysis incorporating previously defined predictors, EF1 was the strongest predictor of response (OR: 1.563, 95% CI: 1.371–1.782, p&lt;0.001) (table 1). EF1 was also the strongest predictor of a clinical composite score (OR: 1.115, 95% CI: 1.044–1.191, p=0.001). Improvement in EF1 at 6 months after CRT implantation (6.5±5.8% vs 1.8±4.3% in responders vs. non-responders, p&lt;0.001) was the best predictor of heart failure re-hospitalization and death after median follow-up of 29.2 months (HR: 0.819, 95% CI:0.765–0.876, p&lt;0.001).

Conclusion

EF1 is a promising marker to identify patients likely to respond to CRT and most likely provides a measure of myocardial viability that determines response.

Funding Acknowledgement

Type of funding source: Foundation. Main funding source(s): British Heart Foundation, Wellcome/EPSRC Centre for Medical Engineering

A joint analysis strategy reveals genetic changes associated with artificial selection between egg-type and meat-type ducks

Egg-type ducks and meat-type ducks are predominantly commercial or indigenous and have been subjected to artificial directional selection. These two duck types differ substantially in body shape, production performance and reproductivity. However, the genetic changes associated with phenotypic differences remain unclear. Here, we compared the two duck types at the genomic and transcriptomic levels. We identified a large number of SNPs and genes in genomic divergent regions in terms of F_ST and θπ values. The corresponding genes were mainly enriched in embryonic development function and metabolic pathway. RNA-seq analysis also revealed differential gene expression in the liver and gonads. The differentially expressed genes were functionally associated with signal transmission and substance metabolism respectively. Furthermore, we found that seven genes were related to differentiation between the two types by both g genome and transcriptome analysis and were plausible candidate genes. These genes were annotated to GO categories of cell development and disease immunity. These findings will enable a better understanding of the artificial selection history of meat and egg ducks and provide a valuable resource for future research on the breeding of these two lineages.

LncRNA KCNQ1OT1 delayed fracture healing through the Wnt/β-catenin pathway

OBJECTIVE

The aim of this study was to explore the effect of long non-coding ribonucleic acid (lncRNA) KCNQ1 overlapping transcript 1 (KCNQ1OT1) on fracture healing and its possible mechanism.

PATIENTS AND METHODS

Abnormal lncRNAs were compared between patients with delayed fracture healing and those with normal fracture healing using gene expression profiling method. LncRNA expression in patients was verified by quantitative Reverse Transcription-Polymerase Chain Reaction (qRT-PCR). Subsequently, the model of tibial fracture was successfully established in rabbits. The effect of lncRNA KCNQ1OT1 expression on tibial fracture healing in rabbits was explored. Meanwhile, the effects of lncRNA KCNQ1OT1 on cell proliferation and apoptosis were investigated by knockdown and overexpression experiments with HC-a as a cell model. Furthermore, Western blotting was used to explore the expressions of proteins in signaling pathway affected by lncRNA KCNQ1OT1.

RESULTS

Gene expression profiling and qRT-PCR revealed that lncRNA KCNQ1OT1 was significantly down-regulated in bone tissues of patients with delayed fracture healing. Compared with the control group, knocking down lncRNA KCNQ1OT1 remarkably reduced the serum levels of alkaline phosphatase (ALP) and osteoprotegerin (OPG) in rabbits, and markedly decreased bone trabecular growth index (p<0.05). In HC-a cells, overexpression of lncRNA KCNQ1OT1 activated the Wnt/β-catenin signaling pathway, which could be suppressed by knocking down lncRNA KCNQ1OT1. Cell Counting Kit-8 (CCK-8) assay and 5(6)-carboxyfluorescein diacetate, succinimidyl ester (CFSE) results manifested that lncRNA KCNQ1OT1 remarkably promoted the proliferation and inhibited apoptosis of HC-a cells by activating the Wnt/β-catenin signaling pathway.

CONCLUSIONS

LncRNA KCNQ1OT1 plays a vital role in delayed fracture healing. Moreover, it can induce cell proliferation and inhibit cell apoptosis by activating the Wnt/β-catenin signaling pathway. Therefore, KCNQ1OT1 may be used as a biomarker to predict the occurrence of delayed fracture healing.

On the optimal number of dose-limiting shells in the SBRT auto-planning design for peripheral lung cancer

PURPOSE

The number of dose-limiting shells in the optimization process is one of the key factors determining the quality of stereotactic body radiotherapy (SBRT) auto-planning in the Pinnacle treatment planning system (TPS). This study attempted to derive the optimal number of shells by evaluating the auto-plans designed with different number of shells for peripheral lung cancer patients treated with SBRT.

METHODS

Identical treatment technique, optimization process, constraints, and dose calculation algorithm in the Pinnacle TPS were retrospectively applied to 50 peripheral lung cancer patients who underwent SBRT in our center. For each of the patients, auto-plans were optimized based on two shells, three shells, four shells, five shells, six shells, seven shells, eight shells, respectively. The optimal number of shells for the SBRT auto-planning was derived through the evaluations and comparisons of various dosimetric parameters of planning target volume (PTV) and organs at risk (OARs), monitor units (MU), and optimization time of the plans.

RESULTS

The conformity index (CI) and the gradient index (GI) of PTV, the maximum dose outside the 2 cm of PTV (D2cm ), Dmax of spinal cord (SCmax ), the percentage of volume of total lung excluding ITV receiving 20 Gy (V20) and 10 Gy (V10), and the mean lung dose (MLD) were improved when the number of shell increased, but the improvement became not significant as the number of shell reached six. The monitor units (MUs) varied little among different plans where no statistical differences were found. However, as the number of shell increased, the auto-plan optimization time increased significantly.

CONCLUSIONS

It appears that for peripheral lung SBRT plan using six shells can yield satisfactory plan quality with acceptable beam MUs and optimization time in the Pinnacle TPS.

[Risk factors for death and the clinical features of different subtypes of patients with pulmonary arterial hypertension related to congenital heart disease]

Objective: To explore the risk factors for death in patients with pulmonary arterial hypertension related to congenital heart disease (PAH-CHD) and the clinical characteristics of different subtypes in patients with PAH-CHD. Methods: It was a retrospective study. A total of 507 PAH-CHD patients, who were hospitalized in the Department of Pediatric Cardiology of Beijing Anzhen Hospital between September 2005 and May 2019, were included. Patients were divided into 4 subgroups: (1) Eisenmenger syndrome(ES) group. (2) PAH associated with prevalent systemic-to-pulmonary shunts(SP) group. (3) PAH associated with small defects(SD) group. (4) PAH after defect correction(CD) group. According to the complexity of cardiac malformation, patients were divided into simple-medium complex malformation group and complex malformation group. According to the location of shunts, patients were divided into pre-tricuspid group, post-tricuspid group, and mixed group or complex deformity group. Baseline clinical data of patients were collected from the electronic medical record system. Demographic data (age, gender, etc.), percutaneous oxygen saturation(SPO(2)), New York Heart Association(NYHA) cardiac function classification, 6 minutes walking distance(6MWD), and B type natriuretic peptide(BNP), systolic pulmonary arterial pressure(sPAP) estimated by echocardiography and mean pulmonary artery pressure (mPAP), mean right atrial pressure(mRAP), cardiac index(CI), and calculated pulmonary vascular resistance (PVR) estimated by right heart catheterization were compared among various groups. The results of regular follow-up of all enrolled patients were collected, including the status of monotherapy or combination of PAH-targeted drugs during the follow-up period, cardiac-related adverse events(hemopysis, syncope, edema, arrhythmia, etc.) and primary endpoint event(all-cause death) were obtained and analyzed. Risk factors for all-cause death were analyzed using univariate and multivariate Cox regression analysis model. Results: The median age at diagnosis was 23.1(13.9,32.1) years, 345 cases(68.0%) were female. Two hunderds and thirty-five cases(46.4%) were diagnosed with ES; 193 cases(38.1%) were diagnosed with CD, 47 cases (9.3%) were diagnosed with SD. Among them, 32 cases(6.3%) were in the SP group. All 507 patients underwent echocardiography examination, there were significant differences in sPAP among different clinical subgroups(P<0.001). A total of 289 patients(57.0%) received right heart catheterization examination, the results showed that the ES group had the highest mPAP and PVR and the lowest mRAP(all P<0.001), the CD group had the highest mRAP and CI(both P<0.001). The 6MWD in the ES group was significantly shorter than that in the SP, SD, and CD groups(all P<0.001). The proportion of patients with NYHA class Ⅲ/Ⅳ was higher in SD group than in SP group(P<0.001), which was similar between SD, ES and CD groups (P values were 0.077 and 0.072, respectively). At admission, the proportion of patients with NYHA class Ⅰ/Ⅱwas the highest in SP group(96.9% (31/32) ), followed by CD group (85.5%(165/193)) and the ES group(85.1%(200/235)), and the SD group(75.0%(35/47)). The BNP level at admission was also higher in SD group than in SP, ES and CD groups(P<0.001). Of the 507 patients, 379(74.8%) patients received PAH-targeted drug therapy at the last follow-up, and the treatment plan was mainly monotherapy(75.7%(287/379)). The median follow-up time was 3.6(2.0, 5.6) years and 37(7.3%) patients died, including 13 in the CD group, 17 in the ES group, and 7 in the SD group. No deaths occurred in the SP group. Right heart failure was the most common cause of death(11(29.7%)), followed by severe hemoptysis dyspnea(7(18.9%)), sudden cardiac death(6(16.2%)), and pulmonary hypertensive crisis(4(10.8%)). Kaplan-Meier curve showed that survival rates of end-point-free events at 1, 3, 5 and 10 years after diagnosis of PAH were 98.0%, 95.4%, 89.9%, and 84.4%, respectively; there were statistically significant differences in survival among the subgroups(P=0.026); there was no significant difference in the survival rate between the ES group and the CD group(P=0.918), and both were higher than the SD group(P values were 0.011 and 0.013, respectively). Univariate Cox regression analysis showed that NYHA class Ⅲ/Ⅳ and BNP>100 ng/L at admission were the risk factors for all-cause death in patients with PAH-CHD(HR=6.452, 95%CI 3.378-12.346, P<0.001, and HR=2.481, 95%CI 1.225-5.025, P=0.012). Multivariate Cox regression analysis showed that NYHA class Ⅲ/Ⅳ was an independent risk factor for all-cause death in patients with PAH-CHD(HR=4.998, 95%CI 1.246-20.055, P=0.023). Conclusions: PAH-CHD patients with different clinical subtypes have different clinical symptoms, cardiac functional class, hemodynamic characteristics, and mid to long-term survival rates. SP patients have the best prognosis, outcome of ES and CD patients is similar, and SD patients have the worst prognosis. NYHA class Ⅲ/Ⅳ is an independent risk factor for all-cause death in patients with PAH-CHD.