Objective To investigate the impact of red blood cell suspension infusion across various perioperative periods on patients with valvular heart disease. Methods The patients with valvular heart disease admitted to Tianjin Chest Hospital from 2018 to 2020 were selected. Based on the timing of perioperative red cell suspension infusion, patients were categorized into three groups: a group 1 receiving intraoperative red cell suspension infusion, a group 2 receiving red cell suspension infusion within 24 hours after entering the ICU, and a group 3 receiving red cell suspension infusion at both time points. The laboratory results, perioperative blood component infusion volume, and other relevant parameters were retrospectively analyzed. After propensity score matching, the differences in different variables among the three groups were compared. Results After propensity score matching, 102 patients were enrolled, including 52 males and 50 females, with an average age of (61.74±10.58) years. There were 34 patients in each group. The preoperative hemoglobin (Hb) value of the group 2 was significantly higher than that of the group 1 and the group 3, and the amount of red cell suspension and autoblood transfusion was the lowest (P<0.05). Group 1 had the highest postoperative Hb, as well as the highest Hb and hematocrit (HCT) levels within 24 hours post-surgery (P<0.05). The group 1 had the lowest plasma, platelet and cryoprecipitate infusion volumes, and the shortest cardiopulmonary bypass time, aortic occlusion time, postoperative ICU stay and hospital stay, and the least blood loss and total drainage volume (P<0.05). The difference between postoperative and preoperative Hb (△Hb1) was highest in group 1 (P<0.05). Conclusion For patients with valvular heart disease, intraoperative-only infusion of red blood cell suspension is associated with a better prognosis at discharge and during follow-up.

Objective To systematically evaluate the therapeutic effects of video-assisted thoracoscopic surgery (VATS) and robot-assisted thoracic surgery (RATS) in treating mediastinal tumors. Methods A computer search was conducted on PubMed, Embase, Cochrane Library, Web of Science, Wanfang, CNKI, CBM, VIP databases for literature comparing the clinical efficacy of VATS and RATS in treating mediastinal tumors, with the search time from inception to March 31, 2024. The Newcastle-Ottawa Scale (NOS) was used to evaluate the quality of the included cohort studies, and Review Manager 5.4 software was used to perform a meta-analysis. Results A total of 32 articles were included, with 7868 patients. The NOS scores of the included cohort studies were all≥7 points. Meta-analysis results showed that compared with the VATS group, the RATS group had less intraoperative blood loss [MD=−16.71, 95%CI (−23.88, −9.54), P<0.001], lower conversion rate to open thoracotomy [OR=0.41, 95%CI (0.26, 0.67), P=0.003], lower overall postoperative complication rate [OR=0.66, 95%CI (0.48, 0.92), P=0.01], shorter postoperative drainage time [MD=−0.64, 95%CI (−0.92, −0.36), P<0.001], and shorter postoperative hospital stay [MD=−1.03, 95%CI (−1.28, −0.78), P<0.001]. There was no statistical difference between the two groups in terms of tumor size [MD=−0.06, 95%CI (−0.31, 0.19), P=0.64] or operation time [MD=5.52, 95%CI (−2.35, 13.40), P=0.17]. The RATS group had higher hospitalization costs than the VATS group [MD=1.69, 95%CI (1.26, 2.13), P<0.001]. Conclusion In mediastinal tumors resection, RATS is superior to VATS in terms of intraoperative blood loss, conversion rate to open thoracotomy, overall postoperative complication rate, postoperative drainage time, and postoperative hospital stay, but it increases hospitalization costs.

ObjectiveTo clarify whether METTL14 mediates the core role of acupuncture at Neiguan (PC6) in promoting myelination and improving behavior in young autistic rats through gene intervention technology. MethodsThe ASD model was established by intraperitoneal injection of valproic acid (VPA) in pregnant rats. Male offspring were intracerebroventricularly injected with adenovirus-packaged METTL14 shRNA (sh-METTL14) or its control (sh-NC) on postnatal day 1, with a model group set as well. Subsequently, the juvenile rats were divided into model group, acupuncture group, acupuncture+sh-NC group, and acupuncture+sh-METTL14 group. The acupuncture group received acupuncture at Neiguan (PC6) from postnatal day 7, once daily for 21 consecutive days. Neurobehavioral changes were evaluated by behavioral tests; METTL14 knockdown efficiency and the expression of METTL14, METTL3, and PTEN were detected by quantitative real-time PCR (qRT-PCR) and Western blot (WB); PTEN m⁶A levels were measured by RNA immunoprecipitation-qPCR (RIP-qPCR); myelin ultrastructure, expression of myelin basic protein (MBP) and neurofascin 155 (NF155), and dendritic spine density were observed using transmission electron microscopy (TEM), enzyme-linked immunosorbent assay (ELISA), immunofluorescence, qRT-PCR, and primary neuron culture. ResultsBehaviorally, knockdown of METTL14 significantly counteracted the beneficial effects of acupuncture in improving self-grooming, open field exploration, three-chamber social interaction, and Morris water maze learning and memory (P<0.05, P<0.01). Compared with the acupuncture+sh-NC group, the acupuncture+sh-METTL14 group showed significantly decreased mRNA and protein expression of hippocampal METTL14 (P<0.01), and the upregulating effects of acupuncture on METTL3 and PTEN expression were reversed (P<0.01). Meanwhile, knockdown of METTL14 significantly inhibited the acupuncture-induced increase in PTEN m⁶A levels (P<0.01). Morphologically, knockdown of METTL14 attenuated the improvement of myelin structure by acupuncture, reversed the downregulation of MBP and upregulation of NF155 induced by acupuncture, and blocked the increase in dendritic spine density (P<0.05, P<0.01). ConclusionMETTL14 is a key molecule mediating the therapeutic effect of acupuncture at Neiguan. Acupuncture at Neiguan upregulates METTL14, thereby enhancing m⁶A methylation modification of PTEN mRNA to stabilize its expression, ultimately promoting myelin development and improving behavioral symptoms in ASD juvenile rats. This preliminarily reveals the modern biological connotation of “opening Xuanfu and dredging myelin”.

The initiation of adaptive immune responses relies on the precise recognition and interpretation of antigenic information. In this process, the specific binding of T cell receptors (TCRs) to peptide-major histocompatibility complex (pMHC) molecules represents one of the key molecular events in the initiation of adaptive immune responses. Accordingly, the structural features of TCR-pMHC complexes provide a fundamental basis for dissecting antigen recognition mechanisms and support rational vaccine design, therapeutic target discovery in TCR-based immunotherapy, and TCR identification and optimization. However, experimental determination of TCR-pMHC structures remains costly, time-consuming, and limited in coverage, making computational approaches essential for rapidly obtaining reliable structural information. Computational methods for predicting the structures of TCR-pMHC complexes have advanced rapidly in recent years, driven by progress in deep learning-based modeling frameworks and the increasing availability of structural and sequence resources. Despite these developments, most existing tools do not adequately distinguish the key structural and biophysical differences between MHC class I (MHC-I) and MHC class II (MHC-II) complexes during model construction. As a consequence, their predictive performance differs substantially between class I and class II complexes. In general, structural predictions for class I complexes outperform those for class II complexes. This discrepancy may be related to several fundamental differences between the two systems, including the architecture of the peptide-binding groove, the distribution of peptide lengths, and the properties of peptide flanking residues (PFRs). Compared with MHC-I molecules, MHC-II molecules usually bind longer antigenic peptides, which typically range from 13 to 25 amino acids in length. PFRs at both termini of these peptides participate in regulating the overall conformation of TCR-pMHC class II complexes and exert a pronounced effect on the geometric and physicochemical characteristics of the TCR-pMHC binding interface. Furthermore, within the TCR recognition interface, the complementarity-determining regions (CDRs) consist of segments that differ markedly in conformational behavior. They commonly include regions that are relatively rigid and structurally stable, together with highly flexible segments exhibiting substantial conformational plasticity. These rigidity-flexibility features constitute an essential structural basis enabling TCRs to recognize diverse peptide-MHC ligands and to accommodate conformational heterogeneity at the interface. However, many current modeling tools, in an effort to enforce global conformational stability or reduce structural noise, tend to over-constrain intrinsically flexible regions. Such oversimplification may lead to inappropriate rigidification of flexible CDR loops, resulting in local structural distortions, compromised interface geometry, or even complete modeling failure for specific complexes. Against this background, the review approaches the field from the perspective of computational differences between MHC-I and MHC-II complexes. We first systematically organize and summarize available resources related to TCRs and pMHCs, including structural datasets, sequence databases, prediction tools, and benchmarking studies. We then focus on five representative tools capable of predicting both class I and class II complexes—AlphaFold2, AlphaFold3, TCRmodel2, tFold-TCR, and TCR-pHLA_ModellerS. After excluding structures present in the training sets of these tools, we constructed a benchmark dataset comprising 25 class I and 10 class II TCR-pMHC complexes in the bound state and conducted a systematic evaluation using this dataset. We first employ widely used general evaluation metrics, including All-Atom Root Mean Square Deviation (All-Atom RMSD), Backbone RMSD, Template Modeling score (TM-score), and DockQ, to assess the global conformational accuracy and interface modeling quality of class I and class II complexes. For class II complexes, we propose for the first time a peptide flanking residue deviation index, including the PFRs-Deviation Index (PFRs-DI), N-PFR-Deviation Index (N-PFR-DI), and C-PFR-Deviation Index (C-PFR-DI), to quantitatively characterize conformational deviations in PFRs. In addition, we propose the CDR conformational consistency index (CCC) designed to qualitatively evaluate the ability of prediction tools to capture TCR CDR conformational flexibility. These metrics collectively assess a tool’s ability to model both overall conformation and critical functional regions, thereby addressing the limitations of existing evaluation criteria that overemphasize global structure while inadequately capturing modeling quality in key functional areas. This establishes a unified analytical framework for MHC-I and MHC-II complexes to guide data resource selection, modeling strategy formulation, and evaluation system development. The framework further advances computational modeling and provides crucial support for multi-scale analysis of TCR-pMHC recognition mechanisms and their biological functions.

ObjectiveTo clarify whether METTL14 mediates the core role of acupuncture at Neiguan (PC6) in promoting myelination and improving behavior in young autistic rats through gene intervention technology. MethodsThe ASD model was established by intraperitoneal injection of valproic acid (VPA) in pregnant rats. Male offspring were intracerebroventricularly injected with adenovirus-packaged METTL14 shRNA (sh-METTL14) or its control (sh-NC) on postnatal day 1, with a model group set as well. Subsequently, the juvenile rats were divided into model group, acupuncture group, acupuncture+sh-NC group, and acupuncture+sh-METTL14 group. The acupuncture group received acupuncture at Neiguan (PC6) from postnatal day 7, once daily for 21 consecutive days. Neurobehavioral changes were evaluated by behavioral tests; METTL14 knockdown efficiency and the expression of METTL14, METTL3, and PTEN were detected by quantitative real-time PCR (qRT-PCR) and Western blot (WB); PTEN m⁶A levels were measured by RNA immunoprecipitation-qPCR (RIP-qPCR); myelin ultrastructure, expression of myelin basic protein (MBP) and neurofascin 155 (NF155), and dendritic spine density were observed using transmission electron microscopy (TEM), enzyme-linked immunosorbent assay (ELISA), immunofluorescence, qRT-PCR, and primary neuron culture. ResultsBehaviorally, knockdown of METTL14 significantly counteracted the beneficial effects of acupuncture in improving self-grooming, open field exploration, three-chamber social interaction, and Morris water maze learning and memory (P<0.05, P<0.01). Compared with the acupuncture+sh-NC group, the acupuncture+sh-METTL14 group showed significantly decreased mRNA and protein expression of hippocampal METTL14 (P<0.01), and the upregulating effects of acupuncture on METTL3 and PTEN expression were reversed (P<0.01). Meanwhile, knockdown of METTL14 significantly inhibited the acupuncture-induced increase in PTEN m⁶A levels (P<0.01). Morphologically, knockdown of METTL14 attenuated the improvement of myelin structure by acupuncture, reversed the downregulation of MBP and upregulation of NF155 induced by acupuncture, and blocked the increase in dendritic spine density (P<0.05, P<0.01). ConclusionMETTL14 is a key molecule mediating the therapeutic effect of acupuncture at Neiguan. Acupuncture at Neiguan upregulates METTL14, thereby enhancing m⁶A methylation modification of PTEN mRNA to stabilize its expression, ultimately promoting myelin development and improving behavioral symptoms in ASD juvenile rats. This preliminarily reveals the modern biological connotation of “opening Xuanfu and dredging myelin”.

The initiation of adaptive immune responses relies on the precise recognition and interpretation of antigenic information. In this process, the specific binding of T cell receptors (TCRs) to peptide-major histocompatibility complex (pMHC) molecules represents one of the key molecular events in the initiation of adaptive immune responses. Accordingly, the structural features of TCR-pMHC complexes provide a fundamental basis for dissecting antigen recognition mechanisms and support rational vaccine design, therapeutic target discovery in TCR-based immunotherapy, and TCR identification and optimization. However, experimental determination of TCR-pMHC structures remains costly, time-consuming, and limited in coverage, making computational approaches essential for rapidly obtaining reliable structural information. Computational methods for predicting the structures of TCR-pMHC complexes have advanced rapidly in recent years, driven by progress in deep learning-based modeling frameworks and the increasing availability of structural and sequence resources. Despite these developments, most existing tools do not adequately distinguish the key structural and biophysical differences between MHC class I (MHC-I) and MHC class II (MHC-II) complexes during model construction. As a consequence, their predictive performance differs substantially between class I and class II complexes. In general, structural predictions for class I complexes outperform those for class II complexes. This discrepancy may be related to several fundamental differences between the two systems, including the architecture of the peptide-binding groove, the distribution of peptide lengths, and the properties of peptide flanking residues (PFRs). Compared with MHC-I molecules, MHC-II molecules usually bind longer antigenic peptides, which typically range from 13 to 25 amino acids in length. PFRs at both termini of these peptides participate in regulating the overall conformation of TCR-pMHC class II complexes and exert a pronounced effect on the geometric and physicochemical characteristics of the TCR-pMHC binding interface. Furthermore, within the TCR recognition interface, the complementarity-determining regions (CDRs) consist of segments that differ markedly in conformational behavior. They commonly include regions that are relatively rigid and structurally stable, together with highly flexible segments exhibiting substantial conformational plasticity. These rigidity-flexibility features constitute an essential structural basis enabling TCRs to recognize diverse peptide-MHC ligands and to accommodate conformational heterogeneity at the interface. However, many current modeling tools, in an effort to enforce global conformational stability or reduce structural noise, tend to over-constrain intrinsically flexible regions. Such oversimplification may lead to inappropriate rigidification of flexible CDR loops, resulting in local structural distortions, compromised interface geometry, or even complete modeling failure for specific complexes. Against this background, the review approaches the field from the perspective of computational differences between MHC-I and MHC-II complexes. We first systematically organize and summarize available resources related to TCRs and pMHCs, including structural datasets, sequence databases, prediction tools, and benchmarking studies. We then focus on five representative tools capable of predicting both class I and class II complexes—AlphaFold2, AlphaFold3, TCRmodel2, tFold-TCR, and TCR-pHLA_ModellerS. After excluding structures present in the training sets of these tools, we constructed a benchmark dataset comprising 25 class I and 10 class II TCR-pMHC complexes in the bound state and conducted a systematic evaluation using this dataset. We first employ widely used general evaluation metrics, including All-Atom Root Mean Square Deviation (All-Atom RMSD), Backbone RMSD, Template Modeling score (TM-score), and DockQ, to assess the global conformational accuracy and interface modeling quality of class I and class II complexes. For class II complexes, we propose for the first time a peptide flanking residue deviation index, including the PFRs-Deviation Index (PFRs-DI), N-PFR-Deviation Index (N-PFR-DI), and C-PFR-Deviation Index (C-PFR-DI), to quantitatively characterize conformational deviations in PFRs. In addition, we propose the CDR conformational consistency index (CCC) designed to qualitatively evaluate the ability of prediction tools to capture TCR CDR conformational flexibility. These metrics collectively assess a tool’s ability to model both overall conformation and critical functional regions, thereby addressing the limitations of existing evaluation criteria that overemphasize global structure while inadequately capturing modeling quality in key functional areas. This establishes a unified analytical framework for MHC-I and MHC-II complexes to guide data resource selection, modeling strategy formulation, and evaluation system development. The framework further advances computational modeling and provides crucial support for multi-scale analysis of TCR-pMHC recognition mechanisms and their biological functions.

OBJECTIVE To compare the risk of pneumonitis and interstitial lung disease （ILD） associated with different antibody-drug conjugates （ADC） in the treatment of breast cancer. METHODS CNKI， VIP， PubMed， Embase and other Chinese and English databases， and ClinicalTrials.gov were searched from the inception to June 15， 2025. Randomized controlled trials （RCT） about pneumonitis and ILD associated with ADC （T-DM1， T-DXd， SG， Dato-DXd， SHR-A1811，ARX788， and T-Duo） in the treatment of breast cancer were included. After literature screening， data extraction， and quality assessment， a network meta-analysis was conducted using Stata 17.0 software， and the surface under the cumulative ranking curve（SUCRA） of all interventions were ranked. RESULTS A total of 19 RCTs involving 10 556 patients were included. The overall incidence of pneumonitis with ARX788 and T-DXd was significantly higher than that with T-DM1， T-DM1 plus TPC（T-DM1combined with pertuzumab or atezolizumab）， TPC（treatment of primary care）， and SG （ P ＜0.05）， for grade 1-2 pneumonitis， ARX788 and T-DXd showed significantly higher incidence than T-DM1， T-DM1 plus TPC， and TPC （ P ＜0.05）. For both indicators， ARX788 and T-Duo were ranked as the top two by SUCRA. For the incidence of grade ≥3 pneumonitis， T-DXd and T-DM1 were significantly higher than SG （ P ＜0.05）， T-Duo and Dato-DXd were ranked as the top two by SUCRA. For overall incidence of ILD， ARX788 was significantly higher than T-DM1， SHR-A1811， TPC， and SG （ P ＜0.05）， for the incidence of grade 1-2 ILD， ARX788 was significantly higher than T-DM1， SHR-A1811， and TPC （ P ＜0.05）， for the incidence of grade ≥3 ILD， ARX788 was significantly higher than TPC （ P ＜0.05）. For three indicators above， ARX788 and T-DXd combined with pertuzumab were ranked as the top two by SUCRA. CONCLUSIONS Compared with other ADCs， ARX788 and T-Duo are associated with a higher risk of pneumonitis and ILD in patients with breast cancer.

ObjectiveTo investigate the effect of isoacteoside on the replication process of respiratory syncytial virus （RSV） and its underlying mechanism. MethodsRSV-infected HEp-2 cells （control group） were treated with isoacteoside （isoacteoside-treated group）. Real-time quantitative PCR （qPCR） was used to analyze the mRNA expression levels of the RSV fusion protein （F） gene and matrix protein （M） gene in both groups. Viral titers were determined by plaque assay to assess the inhibitory effect of isoacteoside on viral replication. The number of RSV inclusion bodies formed was determined using cell transfection and laser confocal imaging to evaluate the role of isoacteoside in the formation of inclusion bodies required for RSV replication. 18 seven-week-old female BALB/c mice were randomly divided into three groups： uninfected control group （n=6）， RSV-infected group （n=6）， and RSV-infected group treated with isoacteoside （n=6）. 5 days after infection， the mice were euthanized and lung tissues were collected for pathological analysis to assess the therapeutic effects of isoacteoside on RSV-induced lung tissue injury. ResultsCompared to the control group， the isoacteoside-treated group exhibited the reduced mRNA levels of the RSV F gene （t=17.13， P<0.001） and M gene （t=18.22， P<0.001）， as well as a decrease in viral titer （t=15.32， P<0.001）. The number of inclusion bodies was also significantly reduced （t=16.12， P<0.001）. In the RSV-infected mouse experiment， compared to the RSV-infected group， the isoacteoside-treated group showed a decrease in the levels of inflammatory cytokines IL-4 （t=14.76， P<0.01） and IL-6 （t=21.13， P<0.001） in lung tissues， as well as decreased mRNA levels of the viral load-associated genes F （t=19.52， P<0.001） and M （t=18.76， P<0.001） in lung tissues. Pathological damage to lung tissue caused by RSV infection was alleviated. ConclusionIsoacteoside inhibits RSV replication by interfering with the formation of RSV inclusion bodies and alleviates lung tissue injury induced by RSV infection in mice.

Objective To systematically evaluate the efficacy and safety of video-assisted thoracoscopic surgery (VATS) with different numbers of ports in the treatment of spontaneous pneumothorax. Methods We conducted a comprehensive search of CNKI, PubMed, The Cochrane Library, Web of Science, EMbase, Wanfang Data, and the Chinese Medical Journal Full-text Database for clinical controlled trials on VATS with different port numbers for spontaneous pneumothorax, from their inception to March 2023. Two researchers independently screened the literature and assessed its quality.The Newcastle-Ottawa Scale (NOS) was used to assess the methodological quality of cohort and case-control studies, and the Cochrane risk-of-bias tool was used to evaluate randomized controlled trials (RCT). Meta-analysis was performed using RevMan 5.4.1 software. Results A total of 107 studies were included, comprising 35 RCT, 2 cohort studies, and 70 case-control studies. All cohort and case-control studies included in the analysis had NOS scores≥7. The meta-analysis revealed that compared to two-port VATS (2P-VATS) and three-port VATS (3P-VATS), single-port thoracoscopic surgery (SPTS) was associated with less intraoperative blood loss (SMD=–1.58, 95%CI: –1.93 to –1.22, P<0.001; and SMD=–1.59, 95%CI: –2.03 to –1.14, P<0.001, respectively), shorter postoperative hospital stay (SMD=–1.05, 95%CI: –1.29 to –0.82, P<0.001; and SMD=–1.08, 95%CI: –1.39 to –0.77, P<0.001), shorter duration of postoperative chest tube drainage (SMD=–0.75, 95%CI: –1.00 to –0.50, P<0.001; and SMD=–1.23, 95%CI: –1.72 to –0.75, P<0.001), fewer postoperative complications (OR=0.34, 95%CI: 0.26 to 0.45, P<0.001; and OR=0.47, 95%CI: 0.33 to 0.68, P<0.001), and lower pain scores at 24, 48, and 72 hours after surgery (P<0.05). The operative time for SPTS was shorter than that for 2P-VATS (SMD=–0.53, 95%CI: –0.90 to –0.16, P=0.005) but showed no significant difference compared to 3P-VATS (P=0.21). When comparing 2P-VATS with 3P-VATS, 2P-VATS demonstrated less intraoperative blood loss (SMD=–1.02, 95%CI: –1.81 to –0.22, P=0.01), shorter postoperative hospital stay (SMD=–0.59, 95%CI: –1.11 to –0.06, P=0.03), shorter duration of chest tube drainage (SMD=–0.46, 95%CI: –0.85 to –0.08, P=0.02), fewer postoperative complications (OR=0.36, 95%CI: 0.22 to 0.59, P<0.001), and lower pain scores at 24, 48, and 72 hours after surgery (P≤0.05). Conclusion Both SPTS and 2P-VATS are effective and safe surgical options for spontaneous pneumothorax, deserving further promotion and application in clinical practice. However, due to limitations in the quantity and quality of the included studies, more large-sample, high-quality research is needed to validate these findings.

Objective To systematically evaluate the efficacy of immune checkpoint inhibitors (ICIs) in treating esophageal cancer patients of different genders. Methods Computer searches were conducted on PubMed, Cochrane Library, and Embase databases to collect randomized controlled trial (RCT) on ICIs treatment for esophageal cancer patients from the establishment of the databases to January 25, 2024. Two researchers independently screened the literature and extracted data according to the inclusion and exclusion criteria. The outcome indicators were overall survival (OS) and progression-free survival (PFS), and RevMan 5.4 software was used for meta-analysis. The modified Jadad scoring scale was used to evaluate the quality of the included literature. Results A total of 10 RCT involving 5364 esophageal cancer patients were included in this study, with 2684 patients in the trial group and 2680 patients in the control group. The Jadad scores of the included literature were all ≥6 points, indicating high-quality RCT. Meta-analysis results showed that female esophageal cancer patients receiving ICIs treatment [HR=0.72, 95%CI (0.59, 0.87), P<0.001] had a more significant median OS prolongation than male patients [HR=0.73, 95%CI (0.68, 0.78), P<0.001]; while male patients [HR=0.57, 95%CI (0.52, 0.64), P<0.001] had a more significant PFS prolongation than female patients [HR=0.72, 95%CI (0.55, 0.94), P=0.01]. Female patients treated with ICIs alone [HR=0.66, 95%CI (0.50, 0.87), P=0.003] had a more significant median OS prolongation than male patients [HR=0.79, 95%CI (0.72, 0.87), P<0.001]; while male patients receiving ICIs combined with chemotherapy [HR=0.67, 95%CI (0.61, 0.74), P<0.001] had a more significant median OS prolongation than female patients [HR=0.77, 95%CI (0.59, 1.01), P=0.06]. Conclusion Female patients receiving ICIs have a slight advantage in OS compared to male patients, while male patients have an advantage in PFS. Male patients receiving ICIs combined with chemotherapy have better survival benefits than female patients, while female patients using ICIs monotherapy have better survival benefits than male patients.