[Esophageal dysmotility detection in patients with gastroesophageal reflux-related cough and its influence on the efficacy of anti-reflux therapy]

Objective: To investigate the characteristics of esophageal dysmotility in patients with an initial diagnosis of acid/non-acid gastroesophageal reflux-related cough (GERC), and its correlation with the therapeutic response to anti-reflux treatments to search for the useful indicators to screen patients with chronic cough suitable for anti-reflux therapy. Methods: A total of 173 patients with suspicious GERC who attended the Chronic Cough Specialist Clinic of Department of Pulmonary and Critical Care Medicine, Tongji Hospital, School of Medicine, Tongji University between June 2020 and December 2022 were retrospectively selected for the study. The age of the patients was (45.1±14.6) years old, including 87 males and 86 females. Their demographic characteristics, clinical manifestations, and the results of high-resolution manometry (HRM) and multichannel intraluminal impedance-pH monitoring (MII-pH) were collected. Information on the etiological identification process and final diagnosis was also recorded. The recruited cases were grouped according to therapeutic outcomes and divided into cases with a favourable response to conventional anti-reflux treatment, cases with a favourable response to intensified anti-reflux treatment, and cases with no response to anti-reflux treatment. Factors influencing the efficacy of anti-reflux treatment were investigated. Differences between groups were compared using the χ2 test, Student-Newman-Keuls test, and Kruskal-Wallis H(K) test, where applicable. Logistic regression analysis using forward stepwise regression based on maximum likelihood estimation was used to screen for influence factors. Results: The 175 patients with suspicious GERC included 45 (26.0%) patients who responded to conventional anti-reflux treatment, 54 (31.2%) who responded to intensified anti-reflux treatment and 74 (42.8%) who did not respond to anti-reflux therapies. Esophageal dysmotility was present in 52.0% of patients (90/173), but was less common in patients who responded to conventional anti-reflux treatment (χ2=8.09, P=0.018). Although the majority of reflux episodes were non-acid (136/173, 78.6%), the proportion of acid reflux (χ2=19.49, P<0.001) and acid exposure time (H=11.04, P=0.004) were significantly higher in patients who responded to conventional anti-reflux treatment. The patients with acid and non-acid GERC had comparable proportion of esophageal dysmotility (64.9% vs. 48.5%, χ2=3.11, P=0.078), with a shorter break [2.4 (0.7, 5.6) cm vs. 6.1 (1.4, 10.0) cm, Z=-2.39, P=0.017], longer upper esophageal sphincter [(4.1±0.9) cm vs. (3.7±1.3) cm, t=-2.09, P=0.038], higher percentage of normal esophageal contractions [60.0% (17.8%, 90.0%) vs. 30.0% (0, 80.0%), Z=-2.14, P=0.032], and lower percentage of large break [10.0% (0, 40.0%) vs. 50.0% (0, 100.0%), Z=-2.92, P=0.004] in the patients with non-acid GERC. The mean resting pressure of the lower esophageal sphincter was significantly lower (H=7.49, P=0.024), while the percentage of ineffective esophageal contractions was markedly higher (H=8.60, P=0.014) in the patients who responded to intensified anti-reflux treatment and in the patients who did not respond to the anti-reflux therapies. Multifactorial logistic regression analysis identified the percentage of ineffective contraction as an independent factor predicting the efficacy of conventional anti-reflux treatment, with a cut-off value of≤45% and a moderate predictive value (AUC=0.67, P=0.004). Conclusions: Esophageal dysmotility is common in GERC patients with different characteristics in acid and non-acid GERC. The percentage of ineffective esophageal contraction may be a useful indicator for selecting an anti-reflux strategy and predicting treatment outcomes.

Measurement of oesophageal hiatus surface area by multiplanar reconstruction of MDCT: relationship with lower oesophageal sphincter pressure and acid reflux

AIM

To evaluate the relationship between oesophageal hiatus surface area (OHSA) and gastro-oesophageal reflux disease (GERD).

MATERIALS AND METHODS

Patients who underwent 24-h pH monitoring, oesophageal high-resolution manometry, and upper abdominal contrast-enhanced multidetector computed tomography (MDCT) during 2014-2021 were enrolled. Patients with a hiatus hernia (HH) on MDCT or who had a history of gastro-oesophageal surgery were excluded. Multiplanar reconstruction (MPR) of the MDCT image was used for the measurement of OHSA. Correlations of OHSA with acid exposure time (AET) and lower oesophageal sphincter (LOS) pressure of all patients were analysed.

RESULTS

Seventy-eight patients were included in the study. OHSA was much less in the AET <4% group than in the AET >6% group (1.61 ± 0.42 versus 2.09 ± 0.55 cm2, p<0.001). Correlation analysis reveals that OHSA correlated positively with AET (correlation coefficient = 0.47, p<0.001). Receiver operating characteristic (ROC) curve analysis reveals that OHSA can significantly distinguish patients in different groups divided by AET (area under the ROC curve [AUC] = 0.76, 95% confidence interval [CI]: 0.63-0.90). OHSA was not related to LOS pressure (correlation coefficient = -0.268, p=0.051). There was no difference in OHSA between the low LOS pressure group and the normal LOS pressure group (1.84 ± 0.61 versus 1.74 ± 0.50 cm2, p=0.52).

CONCLUSIONS

OHSA significantly correlated with AET but has no relationship with LOS pressure. It may be an independent risk factor of GERD.

Evaluation of a Deep-Learning Auto-Segmentation Model of Cardiac Substructures

PURPOSE/OBJECTIVE(S)

Increasing evidence has suggested that limiting dose not only to the whole heart but also to cardiac substructures can potentially reduce cardiac toxicities. Manual contouring of cardiac substructure can be challenging and time-consuming. To address this concern, we developed a deep learning (DL) model, trained on convolutional neural network algorithms in large external datasets, for auto-segmentation of cardiac substructures. This study aimed to evaluate the quality of the cardiac substructure contours generated by the DL algorithm.

MATERIALS/METHODS

We identified 28 patients with esophagus or gastroesophageal junction cancer from a single institution who received radiation to the esophagus between January 2017 and December 2022. For each case, the DL-generated cardiac substructures (4 heart chambers - left/right atrium [L/RA] and L/R ventricle [L/RV], 4 coronary arteries - L common [LCA], L anterior descending [LAD], L circumflex [LCx], and R common [RCA], and great vessels - ascending aorta [AA], pulmonary artery [PA], and superior vena cava [SVC]) were modified by two radiation oncologists (RO) using the contouring atlas developed by Duane et al. Spatial overlapping of the contours were then assessed using the Dice similarity coefficient (DSC), 95th percentile Hausdorff distance (HD-95), and normalized surface dice at 2 mm tolerance (NSD-2).

RESULTS

The mean values of DSC, HD-95, and NSD-2 are shown in Table 1. Overall, the mean DSC, HD-95, and NSD-2 for the heart chambers ranged from 0.82 to 0.92, 0.40 cm to 1.52 cm, and 0.68 to 0.85, respectively. Ranges of the mean DSC, HD-95, and NSD-2 for the coronary vessels were 0.41 to 0.74, 0.18 cm to 0.98 cm, and 0.66 to 0.77, respectively. Lastly, comparison of the great vessel contours yielded the following ranges for mean DSC, HD-95, and NSD-2 respectively: 0.72 to 0.92, 0.30 cm to 1.64 cm, and 0.65 to 0.83.

CONCLUSION

Our study demonstrates that auto-segmentation of cardiac substructures by DL-powered models can be comparable to manual contours for certain cardiac substructures, namely the four heart chambers and great vessels. Further improvement of the DL on contouring of coronary vessels would be needed prior to the autosegmentation model being widely adopted.

Development and Validation of a Deep Learning-Based Auto-Delineation of Target Volume and Organs at Risk in Pancreatic Cancer Radiotherapy

PURPOSE/OBJECTIVE(S)

The delineation of the clinical target volume (CTV), gross target volume (GTV) and organs at risk (OARs) is a crucial and laborious in pancreatic cancer radiotherapy. In this work, we propose and evaluate a three-dimensional (3D) novel convolutional neural network (CNN) for automatic and accurate CTV, GTV and OARs in pancreatic cancer.

MATERIALS/METHODS

A total of 120 computed tomography (CT) scans patients with pancreatic cancer were collected. A novel 3D CNN network, called ResUNet3D, was developed to achieve auto-delineation. 96 patients chosen randomly were used for training, 12 patients for validation, and 12 patients for testing. Meanwhile, the Dice similarity coefficient (DSC) and 95th percentile Hausdorff distance (HD95%) were used to assess the performance.

RESULTS

The DSC values for the test data were 80.9±8.6%, 77.5±5.6%, 94.5±1.3%, 66.2±13.4%, 73.6±7.6%, 79.0±8.7%, 94.1±1.9%, 94.6±1.4%, 87.3±5.8% for CTV, GTV, liver, duodenum, spinal cord, bowel, kidney left, kidney right, stomach. The corresponding HD95% values were 10.7±6.9mm, 7.8±5.7mm, 11.6±5.6mm, 18.6±5.6mm, 2.7±0.7mm, 17.7±8.6mm, 3.9±1.4mm, 3.7±1.9mm, 13.4±5.7mm, respectively. The average delineation time for one patient's CT images was within 5 seconds.

CONCLUSION

The experimental results demonstrate that the CTV, GTV and OARs delineated for pancreatic cancer by ResUNet3D achieved a close agreement with the ground truth. ResUNet3D could significantly reduce the radiation oncologists' contouring time.

Pilot Study of a Novel Ring Gantry-Based PET/CT Linear Accelerator in Patients with Prostate Cancer Receiving [18F]-DCFPyL for PSMA PET Imaging

PURPOSE/OBJECTIVE(S)

The RefleXion X1® system is a hybrid PET imaging-radiotherapy system that uses real-time positron emissions from a PET tracer to deliver biologically guided radiotherapy (BgRT). This study (NCT05470699) evaluated the hypothesis that the X1 PET imaging subsystem would be able to detect [18F]-DCFPyL PSMA PET signal sufficient to generate a deliverable BgRT plan in patients with prostate cancer.

MATERIALS/METHODS

Patients with prostate cancer scheduled for a diagnostic [18F]-DCFPyL PSMA PET scan as part of standard of care were eligible. Upon completion of the diagnostic PSMA PET scan, images were transferred to the radiotherapy planning system for target identification and contouring. If at least one PET avid tumor lesion was identified, the patient was then scanned on the X1 unit. BgRT planning was performed on each X1 scanned patient. The target lesion volume, activity concentration (AC) and normalized target signal (NTS) were acquired. Successful and deliverable BgRT plans required that the target AC was ≥ 5 kBq/ml and NTS ≥ 2.7.

RESULTS

Twenty-six patients underwent [18F]-DCFPyL PET scans (13 with rising PSA after surgery or radiotherapy, 6 with known metastases and 7 with newly diagnosed high-risk prostate cancer). Median (range) PSA was 3.40 (0.04-122). In 16 patients a PET avid tumor was identified and contoured for planning (4 lymph nodes, 5 bone, 6 prostate gland, and 1 prostate bed). In 13 patients the target lesion was visualized on the X1 PET scan, while in 3 patients the target lesion was too close to the bladder to be clearly visualized. BgRT planning was feasible and met standard of care published SBRT organ dose constraints in 8 patients (3 prostate gland, 3 bone, 2 lymph nodes). BgRT planning was not feasible in 8 patients due to insufficient AC, low NTS or proximity of the target lesion to the PET avid bladder. The accompanying table compares median (range) target volume, AC and NTS for feasible versus not feasible plans.

CONCLUSION

This is the first study to investigate the feasibility of using [18F]-DCFPyL PET imaging for BgRT plan generation on the X1 system in patients with prostate cancer. Lesions that are relevant to radiotherapy of prostate cancer can be visualized including lymph node and bone metastases. A dedicated BgRT workflow with PSMA PET imaging on the X1 at 60 minutes post injection will result in higher target AC and will optimize BgRT planning. PET avid lesions < 1 cm or close to the bladder may make BgRT planning challenging. [18F]-DCFPyL-guided BgRT is technically feasible using the RefleXion X1. BgRT using targeted PET radiopharmaceuticals to biologically guide external beam radiotherapy represents a promising new dimension in radiation oncology and warrants further investigation.

Radiotherapy for Patients with Locally Advanced Esophageal Squamous Cell Carcinoma Receiving Neoadjuvant Immunotherapy Combined with Chemotherapy

PURPOSE/OBJECTIVE(S)

With the success of immunotherapy in advanced esophageal cancer, neoadjuvant chemo-immunotherapy (CIT) is being increasingly used for local staged esophageal cancer, especially in the context of clinical trials, which brings similar pCR with neoadjuvant chemoradiotherapy and shows promising results. However, there is still a part of potentially operable patients can't undergo surgery after neoadjuvant chemo-immunotherapy. The follow-up treatment and prognosis of this population remain unclear.

MATERIALS/METHODS

Patients pathologically diagnosed with ESCC, clinical stage T1-3N+M0 or T3-4aNanyM0(AJCC 8th), PS 0-1 were retrospectively enrolled from 1/2020 to 6/2021 in Zhejiang Cancer Hospital. All patients firstly received PD-1 inhibitors (Camrelizumab, Sintilimab or Tislelizumab) plus chemotherapy (albumin paclitaxel,260 mg/m²on day 1 plus carboplatin AUC = 5 on day 1) every 3 weeks for 2-4 cycles. For those patients who did not receive surgery, definitive radiotherapy with 50.4Gy/28F or 50Gy/25F was adopted using VMAT, concurrent with chemotherapy or alone. The concurrent chemotherapy regimens included weekly TC (paclitaxel 50 mg/m 2, d1, carboplatin AUC = 2, d1) or S1 (60mg bid d1-14,29-42). The survival outcomes and treatment toxicity were recorded and analyzed.

RESULTS

A total of 56 eligible patients were finally identified from 558 patients who were treated in department of thoracic surgery, 31 patients showed no response to neoadjuvant CIT (6 with PD and 25 with SD), 25 patients achieved PR but did not receive surgery due to poor performance status or refuse to operation. Median age was 66(IQR 56-72) and 55(98.2%) were males. 12(19.6%) were stage II and 44(80.4%) were stage III. Among all the patients,25 (44.6%) received radiotherapy alone, and 31 (55.4%) received chemoradiotherapy after neoadjuvant CIT. The median follow-up was 11.8 months (IQR 8.6-20.1). The median PFS and OS were 16.5 months (95CI 12.9-21.4) and 18.6 months (95CI 11.2-NA), respectively. In the subgroup analysis, the median PFS for patients with PR to CIT was 20.2 moths (95CI:17.23-NA), and 12.9 moths (95CI: 0.68-20.4) for patients with SD or PD, HR was 0.45 (95CI:0.22- 0.93, P = 0.027). No significant difference was observed for patients received radiotherapy alone or chemoradiotherapy with HR = 1.36(95CI:0.69-2.71, P = 0.37). The most common AEs observed during this study were anemia (98.2%), Leukopenia (83.9%), Thrombocytopenia (53.6%). Adverse events of grade≥3 radiation-induced pneumonitis and esophagitis were 12.5% and 32.1%, especially, 6 patients (10.7%) died from esophageal fistula and 2 patients (3.6%) died from grade 5 pneumonitis.

CONCLUSION

For local advanced ESCC patients after neoadjuvant CIT who did not receive surgery, definitive radiotherapy was an optional treatment strategy. However, those patients with no response to CIT also showed poor response to radiotherapy, and particular attention should be paid to treatment related toxicity, especially esophageal fistula.

Dosimetric Plan Evaluation of Biology Guided Radiotherapy Using [18F]-DCFPyL PSMA Radiotracer in Patients with Prostate Cancer

PURPOSE/OBJECTIVE(S)

The X1 system represents a cutting-edge solution in radiotherapy delivery, with its capability to perform Biology Guided Radiotherapy (BgRT). The system utilizes real-time positron emission tomography (PET) signal as biological fiducials to provide tracked dose delivery and is initially available for use with [18F]-Fluorodeoxyglucose (FDG). The aim of this research study is to assess the quality of BgRT treatment plans for prostate cancer using patients' PSMA PET images obtained on the X1 system.

MATERIALS/METHODS

Sixteen patients with at least one PET-avid tumor identified on their whole-body diagnostic PSMA PET scan were selected. These patients were scanned on X1 following their diagnostic scan without additional radiotracer administration. Based on the X1 PET images, a BgRT plan was created for each patient, with the prescription dose determined by the location of treatment sites. The planning objectives of organs-at-risk (OARs) were established in accordance with the 2018 Timmerman guidelines. Target coverage objective was the dose covering 95% (D95%) of the planning target volume (PTV) to be higher than 100%. The following parameters were analyzed: PTV D95%, the minimal dose (Dmin) of gross tumor volume (GTV), plan maximum dose (Dmax), conformity index (CI), gradient index (GI), and maximum point dose (D0.03cc) to the nearest OARs. The X1 BgRT planning system also generated dose volume histogram (DVH) bounds, which model variations in BgRT delivery. The low boundary of GTV Dmin, representing the minimum GTV dose in the worst-case scenario, was recorded.

RESULTS

BgRT plans were created for all patients, except for one where the target signal was indistinguishable from the bladder. The prescription dose was 2700 cGy or 3000 cGy in 3 fractions for lymph node lesions, 2400 cGy to 3000 cGy in 3 fractions for bone metastasis, and 4500 cGy in 5 fractions for lesions in prostate. All plans met the dose constraints for OARs as per the Timmerman guidelines. The Dmax of all plans was 129.9% ± 6.9% (mean ± standard deviation). The PTV D95% and GTV Dmin were 101.7% ± 1.0% and 111.0% ± 7.6%, respectively. The low boundary of GTV Dmin was 95.9% ± 5.8%. The CI and GI were 1.22 ± 0.11 and 9.40 ± 2.12, respectively. The D0.03cc to nearest OARs was 84.6% ± 25.4%. The estimated treatment time was 699 ± 228 seconds.

CONCLUSION

This study is a pioneering effort to evaluate the quality of BgRT plans for prostate cancer patients using the [18F]-DCFPyL PSMA radiotracer. Our results showed that all BgRT plans met the planning objectives defined in the Timmerman protocol. BgRT with [18F]-DCFPyL represents a promising treatment modality for patients with prostate cancer. Further research is needed to validate this approach, including a comprehensive assessment of the dosimetric and tracking accuracy through physical measurements.

Different Radiation Dose of Neoadjuvant Chemoradiation for Resectable Thoracic Esophageal Squamous Carcinoma: A Randomized Phase II Clinical Trial

PURPOSE/OBJECTIVE(S)

Radiation dose used in the neoadjuvant chemoradiotherapy (nCRT) for patients with locally advanced esophageal squamous cell carcinoma (ESCC) varies in clinical practice, no prospective studies have been conducted on this issue. The hypothesis of this study is that higher radiation dose in nCRT was associated with improved progression free survival in resectable ESCC.

MATERIALS/METHODS

A phase II prospective randomized controlled trial was conducted from February 22, 2018 to February 22, 2021 in Zhejiang Cancer Hospital (NCT03381651). Patient with locally resectable advanced thoracic ESCC (cT2-T4aNxM0, T1-3N+M0) were randomized 1:1 to receive nCRT with different radiation doses of 50.4 Gy/28F and 41.4 Gy/23F, weekly TC regimen (paclitaxel 50 mg/m2, d1, carboplatin AUC = 2, d1) was administered as concurrent chemotherapy. The primary endpoint was 2-year progression free survival (PFS), and the secondary endpoints included complete pathological response (pCR) rate, overall survival (OS), local recurrence free survival (LRFS), and distant metastasis free survival (DMFS). The last follow-up date was February 22, 2022 and data analysis was performed from May 1, 2022, to October 31, 2022.

RESULTS

A total of 147 patients were enrolled in the prospective study, including 72 patients in 50.4 Gy/28F group and 75 patients in 41.4 Gy/23F group. Among all the patients, 101 patients finally underwent surgical resection, 23/46(50%) achieved pCR in the high-dose group, while 18/55(32.7%) achieved pCR in the low-dose group (P = 0.10). In intention-to-treat population analysis, 2-year PFS in the high-dose and low-dose groups were 53.5% (95% CI: 42.9-66.7%) and 45.5% (95% CI: 35.4%-58.5%), respectively, and hazard ratio (HR) was 0.72 (95% CI: 0.46-1.12, P = 0.14). Median PFS were 44.4 months (95% CI: 16-NA) versus 20.8 months (95% CI: 15-NA) in high and low dose groups. 2-year OS in the high and low-dose groups were 61.8% (95% CI: 51.2 -74.5%) and 60.5% (95% CI: 50.3 -72.8%), respectively, with a HR = 0.78(95% CI: 0.48-1.27, P = 0.32). Median OS were not arrived in both groups up to the last follow-up. In addition, 2-year LRFS were 87.9% (95% CI: 78.4-98.5%) and 74.8% (95% CI: 63.3-88.3%) in high and low-dose groups, while 2-year DMFS were 70.6% (95% CI: 58.3-85.5%) and 73.6% (95% CI: 61.7-87.8%), respectively. No significant statistical difference was observed between the two groups in terms of side effects during nCRT and postoperative complications.

CONCLUSION

Our results showed that neoadjuvant high-dose radiotherapy failed to improve overall survival in ESCC. However, it was found that pCR and local tumor control rates tended to be increased at 50.4 Gy/28F, compared to 41.4 Gy/23F, although the trend was not statistically significant. There was no significant difference in treatment side effects between the two groups. Longer-term follow-up is required to verify our findings.

The radiomic biomarker in non-small cell lung cancer: 18F-FDG PET/CT characterisation of programmed death-ligand 1 status

AIM

To present an integrated 2-[18F]-fluoro-2-deoxy-d-glucose (18F-FDG) positron-emission tomography (PET)/computed tomography (CT) radiomic characterisation of programmed death-ligand 1 (PD-L1) status in non-small-cell lung cancer (NSCLC).

MATERIALS AND METHODS

In this retrospective study, 18F-FDG PET/CT images and clinical data of 394 eligible patients were divided into training (n=275) and test sets (n=119). Next, the corresponding nodule of interest was segmented manually on the axial CT images by radiologists. After which, the spatial position matching method was used to match the image positions of CT and PET, and radiomic features of the CT and PET images were extracted. Radiomic models were built using five different machine-learning classifiers and the performance of the radiomic models were further evaluated. Finally, a radiomic signature was established to predict the PD-L1 status in patients with NSCLC using the features in the best performing radiomic model.

RESULTS

The radiomic model based on the PET intranodular region determined using the logistic regression classifier preformed best, yielding an area under the receiver operating characteristics curve (AUC) of 0.813 (95% CI: 0.812, 0.821) on the test set. The clinical features did not improve the test set AUC (0.806, 95% CI: 0.801, 0.810). The final radiomic signature for PD-L1 status was consisted of three PET radiomic features.

CONCLUSION

This study showed that an 18F-FDG PET/CT-based radiomic signature could be used as a non-invasive biomarker to discriminate PD-L1-positive from PD-L1-negative in patients with NSCLC.

[The official-private game of ginseng trade between China and Korea in modern times]

As culturally friendly countries, China and the Korean Peninsula have had active medical exchanges throughout history, which had a significant impact on the traditional medicine industry and daily life of both sides. Ginseng is an important commodity in the drug trade between China and Korea, and its profits are lucrative and an important source of fiscal revenue for the Korean government.In modern times, lured by huge profits, some Korean and Chinese businessmen have engaged in ginseng smuggling, even government officials.In order to safeguard the interests of the official camp, the two governments have always severely punished smuggling, but ginseng smuggling has not disappeared, and the game between official and private ginseng has become a unique phenomenon in the medical exchanges between China and Korea.

Andrographolide loaded montmorillonite attenuated enterotoxigenic Escherichia coli induced intestinal barrier injury and inflammation in a mouse model

Montmorillonite (MMT), a natural absorbent agent, has widely been accepted for its antidiarrhea function in human and farm animals; however, its specific physicochemical property limits its biological function in practical use. In the current study, raw MMT was loaded by andrographolide, namely andrographolide loaded montmorillonite (AGP-MMT). The microstructure of AGP-MMT was observed by scanning electron microscope (SEM) and X-ray diffraction (XRD). The effect of AGP-MMT on the growth performance, intestinal barrier and inflammation was investigated in an enterotoxigenic Escherichia coli (ETEC) challenged mice model. The results show that the microstructure of MMT was obviously changed after andrographolide modification: AGP-MMT exhibited a large number of spheroid particles, and floccule aggregates, but lower interplanar spacing compared with MMT. ETEC infection induced body weight losses and intestinal barrier function injury, as indicated by a lower villus height and ratio of villus height/crypt depth, whereas the serum levels of diamine oxidase (DAO), D-xylose and ETEC shedding were higher in the ETEC group compared with the CON group. Mice pretreated with AGP-MMT showed alleviated body weight losses and the intestinal barrier function injury induced by ETEC challenge. The villus height and the ratio of villus height/crypt depth, were higher in mice pretreated with AGP-MMT than those pretreated with equal levels of MMT. Pretreatment with AGP-MMT also alleviated the increased concentration of serum tumor necrosis factor-α (TNF-α) and interleukin-1β (IL-1β), and the corresponding genes in the jejunum induced by ETEC infection in mice. The protein and mRNA levels of IL-1β were lower in mice pretreated with AGP-MMT than those with equal levels of MMT. The results indicate that AGP-MMT was more effective in alleviating intestinal barrier injury and inflammation in mice with ETEC challenge than MMT.

[Traditional Chinese medicine for treatment of echinococcosis: a review]

Echinococcosis is a zoonotic parasitic disease caused by infection with Echinococcus species. As the drug of first choice for treatment of echinococcosis, albendazole suffers from problems of large doses and remarkable adverse reactions in clinical therapy. Development of novel drugs against echinococcosis is of urgent need. Recently, great advances have been achieved in the research on traditional Chinese medicine for treatment of echinococcosis. This review summarizes the progress of researches on traditional Chinese medicine for treatment of echinococcosis, aiming to provide insights into development of anti-echinococcosis drugs.

The efficacy of hydrogel containing zinc oxide-loaded and minocycline serum albumin nanopartical in the treatment of peri-implantitis

BACKGROUND

We conducted this animal study to assess the efficacy of the novel hydrogel containing zinc oxide-loaded and minocycline serum albumin nanoparticals (Mino-ZnO@Alb NPs) on peri-implantitis in an experimental mouse model.

MATERIAL AND METHODS

Mino-ZnO@Alb NPs was prepared as previously reported. The peri-implantitis model was successfully established in rats, and the rats were divided into three groups randomly: Mino-ZnO@Alb NPs (Mino-ZnO) group, minocycline group, and untreated group. Four weeks later, clinical and radiographic assessments were performed to evaluate soft tissue inflammation and bone resorption level. Histologic analysis was performed to estimate the amount of remaining supporting bone tissue (SBT) around implants. ELISA tests were used to determine the concentration of inflammation factor interleukin-1-beta (IL-1β) and anti-inflammation factor tumor necrosis factor-alpha (TNF-α) around implants.

RESULTS

After one month, the Mino-ZnO group showed better results than the other two groups in regards to the results of bleeding on probing, probing pocket depth, bleeding index and gingival index. X-ray showed that SBT at mesial and distal sites around implants in the other two groups was significantly lower compared with that of Mino-ZnO group. The quantity of osteoclasts in peri-implant tissues of the Mino-ZnO group was less than that in the minocycline and untreated groups. IL-1β in the Mino-ZnO group was lower than that in the other two groups. TNF-α level was the opposite.

CONCLUSIONS

Mino-ZnO@Alb NPs can effectively treat peri-implantitis and promote soft tissue healing, and may act as a promising product.

Venous and arterial thromboembolism after colorectal cancer in the Netherlands: Incidence, predictors, and prognosis

BACKGROUND

Colorectal cancer (CRC) is the third most prevalent cancer type. CRC-patients are at increased risk of venous and arterial thromboembolism (TE), but the magnitude of the risks, their predictors and consequences are not exactly known.

OBJECTIVES

We aimed to determine incidence, predictors and prognosis of TE after incident CRC in a large, unselected population.

METHODS

Using data from Statistics Netherlands and the Netherlands Comprehensive Cancer Organization, all incident CRC-patients were identified between 2013 and 2018 plus a sample of 1:2 age- and sex-matched control subjects. Incidence rates and cumulative incidences for TE were estimated. Predictor variables for TE were explored by univariable Cox regression. The association between TE and all-cause mortality was evaluated by multivariable time-dependent Cox regression.

RESULTS

68,238 incident CRC-patients were matched to 136,476 controls. CRC-patients had a 1-year cumulative venous TE (VTE) incidence of 1.93 % (95%CI 1.83-2.04), versus 0.24 % (95%CI 0.21-0.27) in controls (HR 8.85; 95%CI 7.83-9.99). For arterial TE (ATE), this was 2.74 % (95%CI 2.62-2.87) in CRC versus 1.88 % (95%CI 1.81-1.95) in controls (HR 1.57; 95%CI 1.47-1.66). Cancer stage, surgery, chemotherapy and asthma were predictors for VTE, whereas age, prior ATE and Parkinson's disease were predictors for ATE. CRC patients with TE had an increased risk of all-cause mortality (VTE HR; 3.68 (95%CI 3.30-4.10, ATE HR; 3.05 (95%CI 2.75-3.39)) compared with CRC-patients without TE.

CONCLUSIONS

This Dutch nationwide cohort study adds detailed knowledge on the risk of VTE and ATE, their predictors and prognosis in CRC-patients. These findings may drive TE prophylactic management decisions.

[Establishment and validation of a risk prediction model for disseminated intravascular coagulation patients with electrical burns]

Objective: To establish and validate a risk prediction model of disseminated intravascular coagulation (DIC) by the screening independent risk factors for the occurrence of DIC in patients with electrical burns. Methods: The retrospective case series study was conducted. The clinical data of 218 electrical burn patients admitted to Baogang Hospital of Inner Mongolia from January 2015 to January 2023 who met the inclusion criteria were collected, including 198 males and 20 females, with the age of (38±14) years. The patients were divided into DIC group and non DIC group based on whether they were diagnosed with DIC during the treatment period. The following data of patients of two groups were collected and compared, including age, gender, total burn area, full-thickness burn area, injury voltage, whether osteofascial compartment syndrome occurred within 1 day after injury, duration of stay in burn intensive care unit, total length of hospital stay, whether combined with inhalation injury and multiple injuries, whether shock occurred upon admission, the abbreviated burn severity index score, and the acute physiology and chronic health evaluation Ⅱ score. The laboratory examination data of the patients within 24 hours after admission were also collected, including blood routine indexes: white blood cell count (WBC), hemoglobin level, platelet count (PLT), and neutrophil count; coagulation indexes: activated partial thromboplastin time (APTT), prothrombin time, thrombin time, and levels of D-dimer and fibrinogen (FIB); blood biochemistry indexes: aspartic transaminase, alanine transaminase, direct bilirubin, total bilirubin, total protein, albumin, blood glucose, creatinine, and urea nitrogen; blood gas analysis indexes: blood pH value, arterial partial pressure of oxygen, arterial partial pressure of carbon dioxide, bicarbonate, and base excess; and cardiac zymogram indexes: levels of myoglobin, troponin, lactate dehydrogenase, creatine kinase (CK), and α-hydroxybutyrate dehydrogenase. Data were statistically analyzed with chi-square test, Fisher's exact probability test, independent sample t test, and Mann-Whitney U test. For the variables with statistically significant differences in single factor analysis, the least absolute value selection and shrinkage operator (LASSO) regression was used to reduce the dimension, and the predictive factors for DIC in 218 patients with electrical burns were screened. The above-mentioned predictors were included in multivariate logistic regression analysis to find out the independent risk factors for DIC in 218 patients with electrical burns, and to draw the prediction model nomograms. The performance of the prediction model was evaluated by the receiver operating characteristic (ROC) curve and the area under the ROC curve, and the prediction model was validated by the calibration curve and clinical decision curve analysis (DCA). Results: Compared with those in non DIC group, the total burn area, full-thickness burn area, total length of hospital stay, and the proportions of high voltage caused injury, occurrence of osteofascial compartment syndrome within 1 day after injury, combination of inhalation injury, and occurrence of shock upon admission of patients in DIC group were significantly increased/prolonged (with Z values of -2.53, -4.65, and -2.10, respectively, with χ2 values of 11.46, 16.00, 7.98, and 18.93, respectively, P<0.05). Compared with those in non DIC group, the APTT, level of D-dimer, myoglobin, WBC, PLT, and levels of FIB, total bilirubin, and CK of patients within 24 hours after admission in DIC group were significantly prolonged/increased (with Z values of -2.02, -4.51, and -3.82, respectively, with t values of -3.84, -2.34, -2.77, -2.70, and -2.61, respectively), and the level of total protein and blood pH value were significantly reduced (t=-2.85, Z=-2.03), P<0.05. LASSO regression analysis was carried out for the above 17 indicators with statistically significant differences. The results showed that injury voltage, the occurrence of shock upon admission, the occurrence of osteofascial compartment syndrome within 1 day after injury, and levels of D-dimer and total protein within 24 hours after admission were predictive factors for the occurrence of DIC in 218 patients with electrical burns (with regression coefficients of 0.24, 0.52, 0.35, 0.13, and -0.001, respectively). Multivariate logistic regression analysis showed that injury voltage, the occurrence of shock upon admission, the occurrence of osteofascial compartment syndrome within 1 day after injury, and D-dimer level within 24 hours after admission were independent risk factors for DIC in 218 patients with electrical burns (with odds ratios of 3.33, 4.24, 2.68, and 1.38, respectively, with 95% confidence intervals of 1.43-7.79, 1.78-10.07, 1.17-6.13, and 1.19-1.61, respectively, P<0.05). Based on the aforementioned four independent risk factors, the nomogram of prediction model for evaluating the probability of DIC in patients was drawn. The area under the ROC curve of prediction model was 0.88, and the 95% confidence interval was 0.82-0.95, indicating that the model had good predictive ability; the curve of prediction model tended to be near the ideal curve, indicating that the model had a high calibration degree; the clinical DCA of prediction model showed that the threshold probability of patients ranged from 4% to 97%, indicating that the model had good predictive ability. Conclusions: The injury voltage, the occurrence of shock upon admission, the occurrence of osteofascial compartment syndrome within 1 day after injury, and D-dimer level within 24 hours after admission are independent risk factors for the occurrence of DIC in patients with electrical burns. The prediction model established based on the above indicators can provide early warning for the occurrence of DIC in these patients.

[A case of intellectual developmental disorder with behavioral abnormalities and craniofacial dysmorphism with or without seizures caused by PHF21A gene variation and review of literature]

Objective: To discuss the clinical and genetic features of intellectual developmental disorder with behavioral abnormalities and craniofacial dysmorphism with or without seizures (IDDBCS). Methods: The clinical and genetic records of a patient who was diagnosed with IDDBCS caused by PHF21A gene variation at Children's Hospital Capital Institute of Pediatrics in 2021 were collected retrospectively. Using " PHF21A gene" as the keyword, relevant articles were searched at CNKI, Wanfang Data and PubMed from establishment of databases to February 2023. Clinical and genetic features of IDDBCS were summarized in the combination of this case. Results: An 8 months of age boy showed overgrowth (height, weight and head circumference were all higher than the 97^th percentile of children of the same age and sex) and language and motor developmental delay after birth, and gradually showed autism-like symptoms like stereotyped behavior and poor eye contact. At 8 months of age, he began to show epileptic seizures, which were in the form of a series of spastic seizures with no reaction to adrenocorticotropic hormone but a good response to vigabatrin. Physical examination showed special craniofacial appearances including a prominent high forehead, sparse eyebrows, broad nasal bridge, and downturned mouth with a tent-shaped upper lip. The patient also manifested hypotonia. Whole exome sequencing showed a de novo heterogeneous variant, PHF21A (NM_001101802.1): c.54+1G>A, and IDDBCS was diagnosed. A total of 6 articles (all English articles) were collected, involving this case and other 14 patients of IDDBCS caused by PHF21A gene variation. Clinical manifestations were intellectual disability or developmental delay (15 patients), craniofacial anomalies (15 patients), behavioral abnormalities (12 patients), seizures (9 patients), and overgrowth (8 patients). The main pathogenic variations were frameshift variations (8 patients). Conclusions: IDDBCS should be considered when patients show nervous developmental abnormalities, craniofacial anomalies, seizures and overgrowth. PHF21A gene variation detection helps to make a definite diagnosis.

[Risk factors for pulmonary atelectasis in adults with tracheobronchial tuberculosis]

Objective: To investigate the risk factors for pulmonary atelectasis in adults with tracheobronchial tuberculosis(TBTB). Methods: Clinical data of adult patients (≥18 years old) with TBTB from February 2018 to December 2021 in Public Health Clinical Center of Chengdu were retrospectively analyzed. A total of 258 patients were included, with a male to female ratio of 1∶1.43. The median age was 31(24, 48) years. Clinical data including clinical characteristics, previous misdiagnoses/missed diagnoses before admission, pulmonary atelectasis, the time from symptom onset to atelectasis and bronchoscopy, bronchoscopy and interventional treatment were collected according to the inclusion and exclusion criteria. Patients were divided into two groups according to whether they had pulmonary atelectasis. Differences between the two groups were compared. Binary logistic regression was used to analyze the risk factors for pulmonary atelectasis. Results: The prevalence of pulmonary atelectasis was 14.7%, which was most common in the left upper lobe (26.3%). The median time from symptom onset to atelectasis was 130.50(29.75,358.50)d, and the median time from atelectasis to bronchoscopy was 5(3,7)d. The median age, the proportion of misdiagnosis of TBTB before admission, and the time from symptom onset to bronchoscopy in the atelectasis group were higher than those without atelectasis, and the proportion of receiving bronchoscopy examination and interventional therapy previously, and the proportion of pulmonary cavities were lower than those without atelectasis (all P<0.05). The proportions of cicatrices stricture type and lumen occlusion type in the atelectasis group were higher than those without atelectasis, while the proportions of inflammatory infiltration type and ulceration necrosis type were lower than those without atelectasis (all P<0.05). Older age (OR=1.036, 95%CI: 1.012-1.061), previous misdiagnosis(OR=2.759, 95%CI: 1.100-6.922), longer time from symptom onset to bronchoscopy examination (OR=1.002, 95%CI: 1.000-1.005) and cicatrices stricture type (OR=2.989, 95%CI: 1.279-6.985) were independent risk factors for pulmonary atelectasis in adults with TBTB (all P<0.05). Of the patients with atelectasis who underwent bronchoscopy interventional therapy, 86.7% had lung reexpansion or partial reexpansion. Conclusions: The prevalence of pulmonary atelectasis is 14.7% in adult patients with TBTB. The most common site of atelectasis is left upper lobe. The TBTB type of lumen occlusion is complicated by pulmonary atelectasis in 100% of cases. Being older, misdiagnosed as other diseases, longer time from onset of symptoms to bronchoscopy examination, and being the cicatrices stricture type are factors for developing pulmonary atelectasis. Early diagnosis and treatment are needed to reduce the incidence of pulmonary atelectasis and increase the rate of pulmonary reexpansion.

[Epidemiological characteristics of tick-borne infectious diseases in Hubei Province, 2016-2021]

Objective: To understand the epidemiological characteristics of tick-borne infectious diseases (TBID) and the risk factors of severe illness and death in Hubei Province from 2016 to 2021. Methods: Based on the incidence data of fever with thrombocytopenia syndrome (SFTS), tsutsugamushi disease, typhus and other TBID reported during 2016-2021, the epidemiological analysis was conducted. Field investigation results of TBID in areas with high incidence in 2021, logistic regression analysis of population characteristics, epidemiological history and other factors were used to explore the risk factors of severe and fatal cases. In the field vector investigation, free ticks and surface ticks of the host animals in the cases' home and surrounding grassland were monitored and detected. Results: A total of 3 826 TBID cases were reported in Hubei from 2016 to 2021, of which 71.30% (2 728/3 826) were SFTS, 13.04% (499/3 826) were tsutsugamushi disease and 15.66% (599/3 826) were typhus. A total of 44 cases died in 6 years; the fatality rate was 1.15% (44/3 826). In the peak seasons of incidence from May to July, the cases in people engaged in agriculture related work accounted for 84.61% (3 237/3 826). The incidence rate in women was higher than that in men, and the cases aged ≥50 years accounted for 81.02% of the total (3 100/3 826), and the incidence rate increased with age (P<0.001). The TBID cases were distributed in 86 counties and districts in 16 prefectures (municipality). The incidence rates of different areas had significant differences (P<0.05), and there was a certain spatial-temporal clustering and expasion. Bovis microplus and Haemaphysalis longicornis were captured in the field, and the positive rates in host animals and grassland ticks were 10.94% (7/64) and 40.00% (2/5), respectively. Univariate logistic regression analysis results showed that age ≥50 years and leukocyte <2.0×10⁹/L were risk factors for severe illness and death. Conclusions: The TBID reported in Hubei were mainly SFTS, tsutsugamushi disease and typhus. In order to reduce the incidence of TBID, it is necessary to strengthen the prevention and control in women aged ≥50 years and reduce field exposure and tick bites during the epidemic period.

[A case of combined oxidative phosphorylation deficiency 32 caused by MRPS34 gene variation and literature review]

Objective: To investigate the clinical features and genetic features of combined oxidative phosphorylation deficiency 32 (COXPD32) caused by MRPS34 gene variation. Methods: The clinical data and genetic test of a child with COXPD32 hospitalized in the Department of Neurology, Children's Hospital, Capital Institute of Pediatrics in March 2021 were extracted and analyzed. A literature search was implemented using Wanfang, China biology medicine disc, China national knowledge infrastructure, ClinVar, human gene mutation database (HGMD) and Pubmed databases with the key words "MRPS34" "MRPS34 gene" and "combined oxidative phosphorylation deficiency 32" (up to February 2023). Clinical and genetic features of COXPD32 were summarized. Results: A boy aged 1 year and 9 months was admitted due to developmental delay. He showed mental and motor retardation, and was below the 3^rd percentile for height, weight, and head circumference of children of the same age and gender. He had poor eye contact, esotropia, flat nasal bridge, limbs hypotonia, holding instability and tremors. In addition, Grade Ⅲ/6 systolic murmur were heard at left sternal border. Arterial blood gases suggested that severe metabolic acidosis with lactic acidosis. Brain magnetic resonance imaging (MRI) showed multiple symmetrical abnormal signals in the bilateral thalamus, midbrain, pons and medulla oblongata. Echocardiography showed atrial septal defect. Genetic testing identified the patient as a compound heterozygous variation of MRPS34 gene, c.580C>T (p.Gln194Ter) and c.94C>T (p.Gln32Ter), with c.580C>T being the first report and a diagnosis of COXPD32. His parents carried a heterozygous variant, respectively. The child improved after treatment with energy support, acidosis correction, and "cocktail" therapy (vitaminB₁, vitaminB₂, vitaminB₆, vitaminC and coenzyme Q10). A total of 8 cases with COXPD32 were collected through 2 English literature reviews and this study. Among the 8 patients, 7 cases had onset during infancy and 1 was unknown, all had developmental delay or regression, 7 cases had feeding difficulty or dysphagia, followed by dystonia, lactic acidosis, ocular symptoms, microcephaly, constipation and dysmorphic facies(mild coarsening of facial features, small forehead, anterior hairline extending onto forehead,high and narrow palate, thick gums, short columella, and synophrys), 2 cases died of respiratory and circulatory failure, and 6 were still alive at the time of reporting, with an age range of 2 to 34 years. Blood and (or) cerebrospinal fluid lactate were elevated in all 8 patients. MRI in 7 cases manifested symmetrical abnormal signals in the brainstem, thalamus, and (or) basal ganglia. Urine organic acid test were all normal but 1 patient had alanine elevation. Five patients underwent respiratory chain enzyme activity testing, and all had varying degrees of enzyme activity reduction. Six variants were identified, 6 patients were homozygous variants, with c.322-10G>A was present in 4 patients from 2 families and 2 compound heterozygous variants. Conclusions: The clinical phenotype of COXPD32 is highly heterogenous and the severity of the disease varies from development delay, feeding difficulty, dystonia, high lactic acid, ocular symptoms and reduced mitochondrial respiratory chain enzyme activity in mild cases, which may survive into adulthood, to rapid death due to respiratory and circulatory failure in severe cases. COXPD32 needs to be considered in cases of unexplained acidosis, hyperlactatemia, feeding difficulties, development delay or regression, ocular symptoms, respiratory and circulatory failure, and symmetrical abnormal signals in the brainstem, thalamus, and (or) basal ganglia, and genetic testing can clarify the diagnosis.

[Early identification and influencing factors of post-traumatic stress disorder in high-stress rescue workers]

Objective: To explore the risk intensity and related influencing factors of post-traumatic stress disorder (PTSD) among high-stress rescue workers, and to provide effective tools for the risk assessment of PTSD in military rescue workers. Method: From June to August 2022, cluster sampling was used to select the high-stress rescue personnel of an Army department as the survey subjects. The acute Stress reaction (ASD) scale and PTSD checklist were used to evaluate the risk of PTSD in military rescue personnel. Multivariate logistic regression were used to analyze the influencing factors of PTSD. Results: The age of 4 460 subjects was (24.38±4.072) years old, including 4 396 males (98.6%). The positive rate of initial screening for ASD was 2.85% (127/4 460). The positive rate of PTSD was 0.67% (30/4 460). Multivariate logistic regression model analysis showed that female, older age, recent trauma exposure history, passive smoking and alcohol consumption were at higher risk of ASD, the values of OR (95%CI) were 4.183 (1.819-9.618), 6.278 (1.363-28.912), 3.094 (1.500-6.379), 2.059 (1.298-3.267) and 2.607 (1.614-4.211), respectively; Lower education level was associated with lower risk of ASD, OR (95%CI) was 0.593 (0.359-0.978); People who are older, thinner, have a history of mental illness, and drink alcohol were at higher risk for PTSD, the values of OR (95%CI) were 20.144 (2.459-165.043), 10.287 (2.218-47.700), 91.104 (8.592-965.980) and 2.866 (1.144-7.180), respectively. Conclusion: Gender, age, education level, passive smoking, alcohol consumption, past history of mental illness and body mass index may be related to the potential risk of PTSD in rescue workers,passive smoking, alcohol consumption, and weight controlling should be focused on to reduce potential risks of PTSD.

[Clinical efficacy of entecavir combined with Biejiajian pills and its influence on TCM syndrome scores during the treatment of chronic hepatitis B with hepatic fibrosis and blood stasis syndrome]

Objective: To investigate the clinical efficacy of entecavir combined with Biejiajian pills and its influence on TCM syndrome scores during the treatment of chronic hepatitis B with hepatic fibrosis and blood stasis syndrome by prospective, randomized and controlled study. Methods: Patients with chronic hepatitis B with hepatic fibrosis and blood stasis syndrome were selected as the research subjects and randomly divided into a treatment group and a control group. Entecavir plus Biejiajian pills or entecavir plus a simulant of Biejiajian pills were given for 48 weeks. The changes in liver stiffness measurement (LSM) and TCM syndrome scores before and after treatment were compared between the two groups to analyze the correlation. The data between groups were analyzed by t-test/Wilcoxon rank sum test or χ(2) test. Pearson correlation coefficient was used to analyze the correlation between TCM syndrome scores and LSM values. Results: After 48 weeks of treatment, the LSM values of the two groups were significantly lower than those of the baseline (P < 0.001), liver fibrosis was significantly improved, and the LSM values of the treatment group were lower than those of the control group [(8.67 ± 4.60) kPa and (10.13 ± 4.43) kPa, t = -2.011, P = 0.049]. After 48 weeks of treatment, the TCM syndrome scores of the two groups were significantly reduced compared with the baseline (P < 0.001), and the clinical symptoms were significantly relieved, and the total effective rates of the improvement of the TCM syndrome scores in the two groups were 74.19% and 72.97%, respectively, but the differences between the groups were not statistically significant (χ(2) = 0.013, P = 0.910). Correlation analysis showed that there was no obvious trend between TCM syndrome scores and LSM values. There were no serious adverse reactions associated with the drug during the observation period of this study. Conclusion: Based on antiviral treatment with entecavir, regardless of whether it is combined with the Biejiajian pill, it can effectively reduce the LSM value, improve liver fibrosis, reduce TCM syndrome scores, and alleviate symptoms in patients with chronic hepatitis B with liver fibrosis and blood stasis syndrome. Compared with entecavir alone, the combined Biejia pill has greater efficacy in improving liver fibrosis and a favorable safety profile, meriting its implementation and widespread application.

[Progress on health-related quality of life and its influencing factors in patients with tuberculosis sequelae]

With the emergence of new tuberculosis patients, the number of patients with tuberculosis sequelae is increasing, which not only increases the medical burden of tuberculosis sequelae year by year, but also affects the health-related quality of life (HRQOL) of patients. The HRQOL of patients with tuberculosis sequelae has gradually received attention, but there are few relevant studies. Studies have shown that HRQOL is related to various factors such as post-tuberculosis lung disease, adverse reaction to anti-tuberculosis drugs, decreased physical activity, psychological barriers, low economic status and marital status. This article reviewed the current situation of HRQOL in patients with sequelae of tuberculosis and its influencing factors, in order to provide a reference for improving the quality of life of patients with sequelae of tuberculosis.

[Prediction of potential suitable habitats of Haemphysalis concinna in Heilongjiang Province based on the maximum entropy model]

OBJECTIVE

To predict the potential suitable habitat of Haemaphysalis concinna in Heilongjiang Province under different climatic scenarios.

METHODS

The geographic locations of ticks in Heilongjiang Province from 1980 to 2022 were captured from literature review and field ticks monitoring data from Harbin Center for Disease Control and Prevention in Heilongjiang Province, and the tick distribution sites with spatial correlations were removed using the software ArcGIS 10.2. The environment data under historical climatic scenarios from 1970 to 2000 and the climatic shared socioeconomic pathways (SSP) 126 scenario model from 2021 to 2040 and from 2041 to 2060 were downloaded from the WorldClim website, and the elevation (1 km, 2010), population (1 km grid population dataset of China, 2010) and annual vegetation index (1 km, 2010) data were downloaded from the Resource and Environmental Science and Data Center, Institute of Geographical Sciences and Natural Resources, Chinese Academy of Sciences. The contribution of environmental factors to H. concinna distribution was evaluated and environmental variables were screened using the software MaxEnt 3.4.1 and R package 4.1.0, and the areas of suitable habitats of H. concinna and changes in center of gravity were analyzed using the maximum entropy model in Heilongjiang Province under different climatic scenarios. In addition, the accuracy of the maximum entropy model for prediction of H. concinna distribution was assessed using the area under curve (AUC) of the receiver operating characteristic curve.

RESULTS

A total of 79 H. concinna distribution sites and 24 environmental variables were collected, and 70 H. concinna distribution sites and 9 environmental factors that contributed to distribution of the potential suitable habitats of H. concinna in Heilongjiang Province were screened. The three most significant contributing factors included precipitation seasonality, annual precipitation, and mean temperature of the driest quarter, with cumulative contributions of 60.7%. The total area of suitable habitats of H. concinna was 29.05 × 10⁴ km² in Heilongjiang Province under historical climatic scenarios, with the center of gravity of suitable habitats located at (47.31° N, 129.16° E), while the total area of suitable habitats of H. concinna reduced by 0.97 × 10⁴ km² in Heilongjiang Province under the climatic SSP126 scenario from 2041 to 2060, with the center of gravity shifting to (47.70° N, 129.28° E).

CONCLUSIONS

The distribution of suitable habitats of H. concinna strongly correlates with temperature and humidity in Heilongjiang Province. The total area of potential suitable habitats of H. concinna may appear a tendency towards a decline with climatic changes in Heilongjiang Province, and high-, medium- and low-suitable habitats may shift.