Objective:To explore the distribution characteristics of glioma-associated oncogene homolog 1 (Gli1) positive cells during orthodontic tooth movement process and conduct a proteomic analysis of these cells.Methods:Forty Gli1-LacZ transgenic mice were used to establish an in vivo orthodontic tooth movement (OTM) model for labeling Gli1 positive cells in Gli1-LacZ transgenic mice (OTM group) and an unforced control group, with tooth movement distance measured using micro-CT. The spatial relationship and distribution characteristics of Gli1 positive cells and H-type vessels of CD31 and endomucin (EMCN) in periodontal tissues were detected by immunofluorescence staining. Twenty Gli1-membrane-targeted tandem dimer Tomato (mT)/membrane-targeted green fluorescent protein (mG) double-genotype mice were bred and Gli1 positive cells were sorted for proteomic sequencing after tamoxifen induction. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases were used for enrichment analysis. Results:The micro-CT three-dimensional reconstruction results showed that the mesial movement of the maxillary first molar in mice after 7 days of force application was (69±15) μm, indicating the successful establishment of the Gli1-LacZ transgenic mouse OTM model. Immunofluorescence staining showed that the blood vessels in periodontal tissue were mostly H-type vessels of CD31 and EMCN. The blood vessels in the periodontal tissues are predominantly H-type vessels positive for both CD31 and EMCN. The percentage of Gli1 positive cells in the OTM group, expressed as (54.5±13.2)%, and the relative fluorescence intensity, expressed as 2.6±0.9, were both significantly greater than those in the control group, which had a Gli1 positive cell percentage of (36.3±9.1)% ( t=3.60 , P=0.002) and a relative fluorescence intensity of 1.0±0.3 ( t=5.20, P<0.001). In contrast to the control group where only a small number of Gli1 positive cells were consistent with the distribution of H-type vessels, in the OTM group the number of Gli1 positive cells increased on the tension side were closely associated with the spatial distribution of H-type vessels. GO enrichment analysis of biological processes found that a large number of proteins in Gli1 positive cells were enriched in pathways such as angiogenesis and tissue remodeling. KEGG enrichment analysis found that related proteins were mainly enriched in pathways related to angiogenesis and Gli1, such as hypoxia-inducing factor 1 signaling pathway, vascular endothelial growth factor signaling pathway and hedgehog signaling pathway. Conclusions:The number of Gli1 positive cells increased on tension side and were closely related to H-type blood vessels in response to mechanical force during orthodontic tooth movement. This may be related to profile of inducing blood vessel formation and tissue remodeling.

Objective:To investigate the association between serum IgG concentration and renal prognosis in patients with IgA nephropathy (IgAN).Methods:It was a multi-center retrospective cohort study, patients with biopsy proven primary IgAN who were recorded in the Chinese IgA Nephropathy Information Registration System between April 1996 and September 2018 were included. Exclusion criteria were: (1) age <18 years; (2) <8 glomeruli in biopsy specimens; (3) estimated glomerular filtration rate (eGFR) <15 ml·min -1·(1.73 m 2) -1 at biopsy; (4) missing baseline serum IgG values; (5) incomplete follow-up data; (6) follow-up duration <12 months. Enrolled patients were divided into 3 groups according to the baseline tertiles of serum IgG: ≤9.50 g/L (G1 group), 9.51-11.99 g/L (G2 group), and ≥12.00 g/L (G3 group). Clinical, and pathological parameters were compared across groups. The endpoint events were defined as doubled serum creatinine level from baseline, or end-stage renal disease (ESRD). Results:A total of 1 976 IgAN patients were included in this study, 631 were in G1 group, 664 in G2 group, and 681 in G3 group. The comparison of baseline clinical data showed that there were statistically significant differences among the three groups in terms of gender, age, microscopic hematuria, edema, body mass index, systolic blood pressure, diastolic blood pressure, hemoglobin, serum creatinine, eGFR, 24-hour urine protein quantity, blood uric acid, blood albumin, serum IgA, serum IgM, the proportion of using immunosuppressants, and the proportion of using glucocorticoids (all P<0.05). In terms of pathology, the higher the serum IgG concentration, the relatively less severe the overall renal pathological damage. The results of univariate Cox regression analysis showed that gender, systolic blood pressure, diastolic blood pressure, hemoglobin, serum creatinine, eGFR, 24-hour urine protein quantity, total protein, serum albumin, globulin, serum IgG, Oxford renal pathological classification, glomerular sclerosis ratio, and glomerular IgM deposition were all associated with the occurrence of renal endpoint events (all P<0.05). Based on clinical practice and previous studies, after adjusting for gender, age, systolic blood pressure, diastolic blood pressure, eGFR, 24-hour urine protein quantity, body mass index, Oxford renal pathological classification, glomerular sclerosis ratio, and the use of renin-angiotensin-aldosterone system inhibitors, glucocorticoids, and immunosuppressants, multivariate Cox regression analysis showed that as a continuous variable, the baseline serum IgG level ( HR=0.91, 95% CI 0.87-0.96) was independently associated with the risk of renal endpoint events in IgAN patients; as a categorical variable, with serum IgG ≤ 9.50 g/L as the reference, serum IgG 9.51-11.99 g/L and serum IgG ≥ 12.00 g/L were independent factors for the occurrence of renal endpoint events in IgAN patients ( HR=0.69, 95% CI 0.49-0.96, P=0.027; HR=0.50, 95% CI 0.34-0.74, P<0.001). During a median follow-up of 33(21, 53) months started from the date of renal biopsy and continued until December 31, 2019, the median follow-up duration was 33 (21, 53) months, and a total of 232 patients (11.74%) reached the composite endpoint. Kaplan-Meier survival analysis showed that the higher the serum IgG concentration in patients with IgAN, the higher their cumulative renal survival rate (Log-rank test, χ2=47.176, P<0.001). Conclusion:The higher level of serum IgG at diagnosis is associated with better clinicopathologic features and renal outcomes, and may portend better renal survival in IgAN patients.

Objective:To construct a nomogram based on preoperative MRI imaging features for the prediction of muscle-invasive upper urinary tract urothelial carcinoma（UTUC）and evaluate its performance.Methods:This retrospective cohort study analyzed the clinical data of 99 UTUC patients treated at the First Affiliated Hospital of Nanjing Medical University from April 2018 to May 2024. Among them，69（69.7%）were male and 30（30.3%）were female，with a median age of 67.0 years. All patients underwent preoperative MRI and radical nephroureterectomy. According to postoperative pathology，tumors staged ≥ T 2 were assigned to the muscle-invasive group，and those staged ≤ T 1 were assigned to the non-muscle-invasive group. Baseline data，pathological information，and imaging characteristics were collected and compared between the two groups. Logistic regression analysis was performed to identify risk factors for muscle-invasive UTUC，and a nomogram was constructed. The diagnostic performance of the model was assessed using receiver operating characteristic（ROC）curves，calibration curves，and decision curve analysis（DCA）. Results:Among the 99 patients，70（70.7%）were diagnosed with muscle-invasive UTUC，and 29（29.3%）with non-muscle-invasive UTUC. The muscle-invasive group had significantly larger tumor size［4.5（2.8，7.0）cm vs. 3.0（2.3，4.5）cm， P = 0.029］，a higher incidence of multifocal tumors［37.1%（26/70）vs. 3.5%（1/29）， P < 0.001］，patchy tumors［30.0%（21/70）vs. 6.9%（2/29）， P = 0.019］，spiculated tumor margins［52.9%（37/70）vs. 17.2%（5/29）， P = 0.001］，tumor compression on renal parenchyma or periureteral/peripelvic fat［68.6%（48/70）vs. 10.3%（3/29）， P < 0.001］，high-grade pathology［92.9%（65/70）vs. 75.9%（22/29）， P = 0.043］，lymph node metastasis［28.6%（20/70）vs. 0， P = 0.001］，and lymphovascular invasion［42.9%（30/70）vs. 10.3%（3/29）， P=0.002］. The apparent diffusion coefficient（ADC）values［0.9（0.8，1.1）× 10 -3 mm2/s vs. 1.1（1.0，1.4）× 10 -3 mm2/s， P < 0.001］and normalized ADC（NADC）values［0.8（0.7，1.0）vs. 0.9（0.8，1.1）， P = 0.002］were significantly lower in the muscle-invasive group. Univariate logistic regression identified multifocality，patchy tumor patterns，spiculated tumor margins，tumor compression on renal parenchyma or periureteral/peripelvic fat，and low NADC values as risk factors for muscle-invasive UTUC（all P < 0.05）. Multivariate analysis revealed multifocality（ OR = 17.903，95% CI 1.650 - 194.253， P = 0.018），tumor compression on renal parenchyma or perirenal / ureteral fat（ OR = 14.690，95% CI 3.069 - 70.323， P < 0.001），and low NADC value（ OR = 0.016，95% CI 0.001 - 0.471， P = 0.017）as independent risk factors. A nomogram was constructed based on these factors. The area under the ROC curve（AUC）of the model was 0.898（95% CI 0.838 - 0.957），with an optimal cutoff value of 0.639. The model showed an accuracy of 83.8%，sensitivity of 81.4%，and specificity of 89.7%. Calibration curves indicated good calibration，and DCA showed that the model provided substantial clinical net benefit. Conclusions:This study constructed a nomogram based on preoperative MRI features，including tumor multifocality，compression on renal parenchyma or periureteral/peripelvic fat and NADC value，which demonstrates good predictive performances for muscle-invasive UTUC.

Objective:To evaluate the diagnostic value of the vesical imaging-reporting and data system（VI-RADS）for determining muscle invasion in variant histology urothelial carcinoma（VUC）of the bladder.Methods:A retrospective analysis was performed on the pathological and imaging data of 518 bladder cancer patients admitted to Jiangsu Province Hospital between January 2013 and January 2023. Patients were stratified into pure urothelial carcinoma（PUC）group（ n = 457）and variant urothelial carcinoma（VUC）group（ n = 61）based on the presence of histological variants. In the PUC group，there were 390 males（85.3%）and 67 females（14.7%），with a mean age of（66.9 ± 11.2）years. Tumor characteristics included maximum diameter ≥ 30 mm in 149（32.6%），< 30 mm in 308（67.4%），multiple tumors in 147（32.2%），solitary in 310（67.8%），pedunculated morphology in 143（31.3%）and non-pedunculated in 314（68.7%）. Histological grading identified high-grade tumors in 319 patients（69.8%）and low-grade tumors in 138（30.2%）. Pathological stage distribution included 191 of T a（41.8%），127 of T 1（27.8%），76 of T 2（16.6%），47 of T 3（10.3%），and 16 of T 4（3.5%）patients. The VUC group included 61 patients，comprising 51 males（83.6%）and 10 females（16.4%），with a mean age of（65.8 ± 11.4）years. Tumor characteristics were maximum diameter ≥ 30 mm in 38（62.3%），< 30 mm in 23（37.7%），multiple tumors in 16（26.2%），solitary in 45（73.8%），pedunculated morphology in 11（18.0%）and non-pedunculated in 50（82.0%）. Histological grading identified high-grade tumors in 59 patients（96.7%）and low-grade tumors in 2（3.3%）. Pathological stage distribution included 3 of T a（4.9%），15 of T 1（24.6%），15 of T 2（24.6%），20 of T 3（32.8%），and 8 of T 4（13.1%）patients. No statistically significant differences were found between the two groups in gender，age，or tumor multiplicity（ P > 0.05）. Statistically significant differences were found in pathological grade，pathological stage，maximum tumor diameter，and pedunculated morphology（ P < 0.05）. Furthermore，an external validation cohort of 278 bladder cancer patients treated between February 2023 and February 2024 from multiple centers（Jiangsu Provincial People’s Hospital，The First Affiliated Hospital of Zhengzhou University，Union Hospital Tongji Medical College Huazhong University of Science and Technology，Jiangsu Provincial Hospital of Traditional Chinese Medicine，Suzhou Municipal Hospital，Huaian First People’s Hospital，Yixing People’s Hospital）was retrospectively analyzed to externally validate the performance of VI-RADS scoring in predicting muscle invasion of VUC. This cohort included a PUC subgroup of 241 patients，comprising 196 males（81.3%）and 45 females（18.7%），with a mean age of（68.0 ± 10.7）years. Tumor characteristics were maximum diameter ≥ 30 mm in 85（35.3%），< 30 mm in 156（64.7%），multiple tumors in 65（27.0%），solitary in 176（73.0%），pedunculated morphology in 76（31.5%）and non-pedunculated in 165（68.5%）. Histological grading identified high-grade tumors in 175 patients（72.6%）and low-grade tumors in 66（27.4%）. Pathological staging comprised 107 patients of T a（44.4%），78 of T 1（32.4%），22 of T 2（9.1%），22 of T 3（9.1%），and 12 of T 4（5.0%）. The VUC subgroup consisted of 37 patients，comprising 29 males（78.4%）and 8 females（21.6%），with a mean age of（70.5 ± 9.5）years. Tumor characteristics were maximum diameter ≥ 30 mm in 23（62.2%），< 30 mm in 14（37.8%），multiple tumors in 9（24.3%），solitary in 28（75.7%），pedunculated morphology in 7（18.9%）and non-pedunculated in 30（81.1%）. Histological grading identified high-grade tumors in 36 patients（97.3%）and low-grade tumors in 1（2.7%）. Pathological staging comprised 1 patient of T a（2.7%），9 of T 1（24.3%），7 of T 2（18.9%），19 of T 3（51.4%），and 1 of T 4（2.7%）. In this validation cohort，no significant differences were found in gender，age，tumor multiplicity，or pedunculated morphology between the PUC and VUC subgroups（ P > 0.05）. Significant differences were observed in pathological grade，pathological stage，and maximum tumor diameter（ P < 0.05）. Three radiologists independently reviewed and scored the multiparametric MRI（mp-MRI）in a blinded manner. Inter-reader agreement was assessed using the weighted kappa statistic. Differences in variables between the two groups were compared using t-tests，chi-square tests，or Fisher’s exact test. The diagnostic performance of VI-RADS for muscle invasion in VUC and PUC was comprehensively evaluated using receiver operating characteristic（ROC）curves，the area under the curve（AUC），and cut-off values determined by the Youden’s index. The DeLong test was used to assess whether the diagnostic performance of VI-RADS differed between VUC and PUC. Results:In the retrospective single-center cohort，the AUC of VI-RADS for assessing muscle invasion was 0.895（95% CI 0.864?0.922）in the PUC group，with a cut-off value of > 3，and the AUC was 0.896（95% CI 0.791-0.960）in the VUC group，with a cut-off value of > 3. The difference between the two groups was not statistically significant（ P = 0.986）. Using a VI-RADS score > 3 as the cut-off value，the accuracy，sensitivity，specificity，positive predictive value（PPV），and negative predictive value（NPV）for diagnosing muscle invasion status in the PUC group were 85.8%（392/457），70.5%（98/139），92.5%（294/318），80.3%（98/122），and 87.8%（294/335），respectively. The corresponding values for the VUC group were 82.0%（50/61），76.7%（33/43），94.4%（17/18），97.1%（33/34），and 63.0%（17/27）.In the retrospective multicenter cohort，the AUC of VI-RADS for assessing muscle invasion was 0.891（95% CI 0.845?0.927）in the PUC group，with a cut-off value of > 2，and the AUC was 0.898（95% CI 0.754?0.973）in the VUC group，with a cut-off value of > 3. The difference between the two groups was not statistically significant（ P = 0.897）. Using a VI-RADS score > 3 as the cut-off value，the accuracy，sensitivity，specificity，PPV，and NPV for diagnosing muscle invasion status in the PUC group were 85.9%（207/241），58.9%（33/56），94.1%（174/185），75.0%（33/44），and 88.3%（174/197），respectively. The corresponding values for the VUC group were 81.1%（30/37），77.8%（21/27），90.0%（9/10），95.5%（21/22），and 60.0%（9/15）.In the single-center cohort，the Kappa values for inter-reader agreement in assessing muscle invasion status using VI-RADS were 0.881（ P < 0.01）for the PUC group and 0.941（ P < 0.01）for the VUC group among the three readers. In the multicenter cohort，the Kappa values were 0.858（ P < 0.01）for the PUC group and 0.838（ P < 0.01）for the VUC group. Conclusions:VI-RADS demonstrates similarly high diagnostic performance for assessing muscle invasion in both PUC and VUC，which is applicable for diagnosing muscle invasion status in VUC，and shows good inter-reader agreement.

Objective:To investigate the association between blood selenium levels and sex hormones in Chinese men aged 18-79 years.Methods:Data were derived from the China National Human Biomonitoring survey conducted in 2017-2018, with a final sample size of 5 414 men. General demographic characteristics, behavioral habits, and dietary frequency were collected through questionnaires and physical examinations. Fasting blood samples were collected to measure blood lead, serum testosterone, and estradiol levels. Complex sampling linear regression models were used to analyze the associations between blood selenium levels and testosterone, estradiol, and the testosterone/estradiol ratio, adjusting for confounding factors including age, education level, marital status, smoking status, alcohol consumption, seafood intake, soy product intake, protein supplement intake, BMI, and diabetes status.Results:The mean age of the 5 414 participants was (46.85±27.91) years; 4 774 (91.65%) were of Han ethnicity and 4 505 (86.68%) were married. The median ( Q1, Q3) blood selenium concentration in men was 97.80 (80.64, 116.99) μg/L. After adjusting for confounding factors, the complex sampling linear regression model revealed negative associations between blood selenium levels and both testosterone levels and the testosterone/estradiol ratio, with a significant linear trend ( Ptrend<0.05). Compared with the Q1 group, the β (95% CI) values for testosterone in the Q2, Q3, and Q4 groups were -0.02 (-0.06 to 0.02), -0.03 (-0.08 to 0.01), and -0.06 (-0.09 to -0.02), respectively. Similarly, the β (95% CI) values for the testosterone/estradiol ratio in the Q2, Q3, and Q4 groups were -0.01 (-0.03 to 0.02), -0.01 (-0.04 to 0.04), and -0.03 (-0.06 to -0.01), respectively. Subgroup analysis indicated stronger associations between blood selenium levels and testosterone/estradiol levels in non-smoking and obese men (BMI≥28 kg/m2). Conclusion:Blood selenium levels are negatively associated with testosterone levels and the testosterone/estradiol ratio in Chinese adult males.

OBJECTIVE To deeply explore the in vivo pharmacodynamic substance basis of Zhigancao Decoction,a classic tradi-tional Chinese medicine formula,and provide scientific evidence for its rational application and development in modern clinical practice.METHODS Wistar rats were treated with 12.15 g·kg-1 Zhigancao Decoction by gavage.Rat plasma samples were collect-ed at 10 time points(5,15,30,60,120,240,360,480,600 and 720 min after administration)and rat heart(atrial and ventricu-lar)tissue samples were collected at 12 h after administration.Components in the plasma and heart samples were qualitatively identi-fied by ultra-performance liquid chromatography-quadrupole-time-of-flight mass spectrometry(UPLC-Q-TOF/MS),and the distri-bution characteristics of Zhigancao Decoction in vivo were analyzed.At the same time,the time-concentration curve of the prototype components and metabolites in Zhigancao Decoction was drawn to observe the changes of blood drug concentration.RESULTS A total of 11 prototype components(Ajugol,Nicotiflorin,Isoschaftoside,4-Hydroxycinnamic acid,Rehmapicrogenin,4-Hydroxybenzoic acid,4′,7-Dihydroxyflavone,Calycosin,3′,4′,7-Trihydroxyflavone,Pinellic acid,Truxillic acid)and 7 metabolites were identified from the plasma samples of Zhigancao Decoction,mainly including flavonoids(flavonoids glycosides),organic acids,and iridoid glyco-sides,etc.Additionally,6 prototype components(Ajugol,Isoschaftoside,Rehmapicrogenin,4′,7-Dihydroxyflavone,Liquiritigenin,3′,4′,7-Trihydroxyflavone)and 3 metabolites were identified from the cardiac samples(the atrium and the ventricle showed the same results).The metabolic pathways mainly involved Phase Ⅰ metabolism and glucuronidation.CONCLUSION The prototype compo-nents and metabolites in plasma and heart tissue of Zhigancao Decoction is preliminarily determined,providing a reference for analyzing the active components of Zhigancao Decoction in heart tissue.

Objective:To evaluate the performance of an artificial intelligence-based endoscopic virtual ruler (EVR) for non-invasive measurement of esophageal varices (EV) diameter.Methods:Patients with liver cirrhosis and EV hospitalized at the First Affiliated Hospital of Anhui Medical University between October 2022 and May 2023 were prospectively enrolled. EV diameter was measured using visual estimation, esophageal varix manometer (EVM), and EVR, with procedure times recorded. The intraclass correlation coefficient (ICC) was used to assess the consistency of EV diameter measurement obtained from the three methods, and repeated-measures ANOVA was used to compare differences in time measurements across three methods.Results:The study included 41 patients with liver cirrhosis and EV. Inter-observer ICC for visual estimation was 0.594, versus 0.840 for EVM and 0.884 for EVR. The ICC value between the EV diameters measured by EVR and EVM was higher than that of the visual assessment. The ICC value between EV diameter measurement by EVM and EVR was 0.991. Measurement times differed significantly across methods ( P<0.001): visual estimation 18.6±2.2 s (14.7-23.3 s), EVR 41.5±4.1 s (31.7-50.3 s), and EVM 170.8±26.4 s (129.3-229.3 s). Repeated measures analysis of variance (corrected using Greenhouse-Geisser) revealed significant differences in time across the three measurement methods [ F(1.033, 41.313)=1 233.800, P<0.001]. Subsequent Bonferroni post-hoc tests revealed significant differences in time between all method pairs ( P<0.001). Conclusion:EVR provides rapid, non-invasive EV diameter measurements with excellent agreement to EVM assessment, offering an efficient alternative to conventional techniques.

Schizophrenia is a common, severe, and complex psychiatric disorder worldwide. Genetic factors account for around 80% of the etiology of schizophrenia, yet objective diagnostic biomarkers remain lacking. This article reports two cases of female monozygotic twins diagnosed with schizophrenia, exhibiting a balanced translocation between 22q11.2 and 4p15.3. Reviewing the literature, we analyze and discuss the correlation between chromosomal balanced translocation regions and the pathogenesis of mental disorders. This aims to encourage psychiatrists to consider new perspectives on the diagnosis of schizophrenia.

Carotid artery perivascular adipose tissue (PVAT) can influence plaque formation and progression. Recently, carotid artery PVAT density has emerged as a novel imaging biomarker being capable of reflecting local metabolic and inflammatory states of adipose tissue. It is closely associated with vulnerable plaque characteristics, such as intraplaque hemorrhage, thinning or rupture of the fibrous cap, lipid-rich necrotic core, and calcification. Therefore, carotid artery PVAT density holds promise as a key parameter for early identification of vulnerable carotid plaques and stroke risk prediction. This article reviews the definition and pathophysiological mechanism of PVAT and application of imaging techniques in PVAT, as well as the association between carotid artery PVAT density and vulnerable characteristics of plaques, with the aim of providing references for early identification of asymptomatic high-risk plaques and individualized prevention strategies of ischemic stroke.

Objective:To investigate the efficacy and safety of lenvatinib plus drug-eluting bead transarterial chemoembolization (DEB-TACE) and FOLFOX-based hepatic arterial infusion chemotherapy (Len+DEB-TACE+HAIC) versus lenvatinib plus DEB-TACE (Len+DEB-TACE) for hepatocellular carcinoma (HCC) larger than 7 cm with portal vein tumor thrombosis (PVTT).Methods:The data from patients diagnosed with HCC (>7 cm) and PVTT who received either Len+DEB-TACE+HAIC ( n=99) or Len+DEB-TACE ( n=102) between July 2019 and June 2021 at six institutions in China were collected and retrospectively analyzed. Tumor responses were evaluated based on modified Response Evaluation Criteria in Solid Tumors. Objective response rate (ORR), disease control rate (DCR), time to progression (TTP), overall survival (OS), and treatment-related adverse event (TRAE) were compared between the two groups by propensity score matching. Subgroup analyses were performed for TTP and OS. Results:After propensity score matching, 83 pairs of patients were included in the study cohorts. The ORR for the Len+DEB-TACE+HAIC group and the Len+DEB-TACE group was 66.3% and 38.6% ( χ2=12.78, P<0.001), respectively. The DCR for the Len+DEB-TACE+HAIC group and the Len+DEB-TACE group was 91.6% and 79.5% ( χ2=4.87, P=0.027), respectively. The median TTP and median OS for the Len+DEB-TACE+HAIC group were significantly longer than those for the Len+DEB-TACE group (TTP, 10.1 months vs. 6.1 months, χ2=35.28, P<0.001; OS, 17.3 months vs. 12.9 months, χ2=16.84, P<0.001). The incidence of ≥grade 3 TRAEs was 38.6% in the Len+DEB-TACE+HAIC group and 33.7% in the Len+DEB-TACE group ( χ2=0.42, P=0.518). Conclusion:Compared with Len+DEB-TACE, Len+DEB-TACE+HAIC led to improved tumor response, TTP and OS with an acceptable safety profile in patients with large HCC and PVTT.