OBJECTIVE To investigate the effects and mechanisms of Lycium ruthenicum Murr. in improving sleep. METHODS Network pharmacology was employed to identify the active components of L. ruthenicum and their associated disease targets， followed by enrichment analysis. A caffeine‑induced zebrafish model of sleep deprivation was established ， and the zebrafish were treated with L. ruthenicum Murr. extract （LRME） at concentrations of 0.1， 0.2 and 0.4 mg/mL， respectively； 24 h later， behavioral changes of zebrafish and pathological alterations in brain neurons were subsequently observed. The levels of inflammatory factors [interleukin-6 （IL-6）， IL-1β， IL-10， tumor necrosis factor-α （TNF-α）]， oxidative stress markers [superoxide dismutase （SOD）， malondialdehyde （MDA）， glutathione peroxidase （GSH-Px）， catalase （CAT）]， and neurotransmitters [5- hydroxytryptamine （5-HT）， γ-aminobutyric acid （GABA）， glutamic acid （Glu）， dopamine （DA）， and norepinephrine （NE）] were measured. The protein expression levels of protein kinase B1 （AKT1）， phosphorylated AKT1 （p-AKT1）， epidermal growth factor receptor （EGFR）， B-cell lymphoma 2 （Bcl-2）， sarcoma proto-oncogene，non-receptor tyrosine kinase （SRC）， and heat shock protein 90α family class A member 1 （HSP90AA1） in the zebrafish were also determined. RESULTS A total of 12 active components and 176 intersecting disease targets were identified through network pharmacology analysis. Among these， apigenin， naringenin and others were recognized as core active compounds， while AKT1， EGFR and others served as key targets； EGFR tyrosine kinase inhibitor resistance signaling pathway was identified as the critical pathway. The sleep improvement rates in zebrafish of LRME low-， medium-， and high-dose groups were 54.60%， 69.03% and 77.97%， 开发。E-mail：hjp_yft@163.com respectively， while the inhibition ratios of locomotor distance were 0.57， 0.83 and 0.95， respectively. Compared with the model group， the number of resting counts， resting time and resting distance were significantly increased/extended in LRME medium- and high-dose groups （P＜0.05）. Neuronal damage in the brain was alleviated. Additionally， the levels of IL-6， IL-1β， TNF-α， MDA， Glu， DA and NE， as well as the protein expression levels of AKT1， p-AKT1， EGFR， SRC and HSP90AA1， were markedly reduced （P＜0.05）， while the levels of IL-10， SOD， GSH-Px， CAT， 5-HT and GABA， as well as Bcl-2 protein expression， were significantly elevated （P＜0.05）. CONCLUSIONS L. ruthenicum Murr. demonstrates sleep-improving effects， and its specific mechanism may be related to the regulation of inflammatory responses， oxidative stress， neurotransmitter balance， and the EGFR tyrosine kinase inhibitor resistance signaling pathway.

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.

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.

Plant-derived extracellular vesicles (PDEVs) are nanoscale extracellular vesicles secreted by plant cells, characterized by a lipid bilayer structure. These vesicles carry a variety of bioactive molecules, including proteins, nucleic acids, and lipids, and play essential roles in intercellular communication and physiological regulation in plants. Compared to animal-derived extracellular vesicles, PDEVs offer several advantages, such as a broad range of sources, high biocompatibility, low immunogenicity, and low production costs. Furthermore, PDEVs have demonstrated remarkable potential as natural nanocarriers for drug delivery, due to their ability to efficiently traverse biological barriers, such as the blood-brain barrier, making them promising candidates for drug delivery systems. This review systematically elaborates on the complex composition of PDEVs, which consists of lipids, proteins, and nucleic acids, the typical structural characteristics of their lipid bilayers ranging from 30 to 150 nm, and their versatile loading capabilities as drug carriers, efficiently encapsulating various types of therapeutic agents such as hydrophilic small molecules, hydrophobic drugs, nucleic acids, and proteins. We systematically summarize the recent advancements in strategies for enhancing the loading efficiency of PDEVs, which include methods such as co-incubation, ultrasound-assisted loading, electroporation, freeze-thaw cycles, and microfluidic technology. These techniques are evaluated based on their underlying principles, suitable drug types, and their respective advantages. In addition to loading strategies, we focus on the engineered approaches to achieve targeted delivery using PDEVs, such as genetic engineering modifications, chemical ligand conjugation, membrane fusion technology, and polyethylene glycol (PEG) modification. We discuss the mechanisms of these strategies in enhancing targeting efficiency, prolonging in vivo circulation time, and improving therapeutic efficacy. Further, this review highlights the application of PDEVs in various disease models, including tumor, skin inflammation, metabolic disorders, and neurodegenerative diseases, showcasing their therapeutic potential as multifunctional delivery platforms. The ability of PDEVs to encapsulate diverse therapeutic agents and target specific tissues or cells opens up new avenues for the treatment of complex diseases, offering advantages over conventional drug delivery systems. However, despite the promising applications of PDEVs, several challenges remain in their development and clinical translation. These challenges include variability in source materials, standardization of preparation processes, quality control, scalability of production, and the need for clinical validation. To overcome these obstacles, the integration of advanced technologies such as artificial intelligence-assisted design and multi-omics analysis is proposed as a way to facilitate the precise development of PDEVs. These emerging technologies hold the potential to further enhance the precision and effectiveness of plant-based drug delivery systems, ultimately advancing the field of precision medicine. In conclusion, the use of PDEVs as a platform for drug delivery represents a promising area of research with the potential to revolutionize therapeutic strategies. Their ability to encapsulate and deliver a wide variety of bioactive molecules, along with their inherent advantages in biocompatibility and versatility, makes them a valuable tool in the development of more efficient and targeted therapeutic interventions. Continued research and innovation in this field will pave the way for the clinical implementation of PDEVs in the treatment of various diseases, offering new hope for more effective and sustainable therapeutic options.

Objective To evaluate the diagnostic efficacy of different pathogen detection methods in active tuberculosis (ATB) of kidney transplant recipients (KTRs) and analyze their diagnostic value in different types of clinical specimens. Methods A retrospective observational study was conducted, including adult KTRs with suspected ATB from the Eighth Medical Center of Chinese People's Liberation Army General Hospital from January 2020 to January 2025. Baseline data of the recipients were collected and analyzed, and the diagnostic criteria were based on a composite reference standard (including clinical diagnosis and bacteriological confirmation). The single-detection efficacy of acid-fast staining smear, mycobacterium tuberculosis culture, Xpert test for the detection of mycobacterium tuberculosis nucleic acid and resistance to rifampicin (Xpert MTB/RIF), as well as tuberculosis or non-tuberculosis mycobacterial DNA detection was compared. The differences in sensitivity among different methods were analyzed. Results A total of 95 KTRs with suspected ATB and 136 specimens were included. And 52 cases were diagnosed as ATB, among which 48 cases were confirmed by bacteriological means (92%). The sensitivity of Xpert MTB/RIF and DNA detection (0.55 and 0.52, respectively) was significantly higher than that of acid-fast staining smear (0.30) and mycobacterium tuberculosis culture (0.35) (P < 0.01). There was a statistically significant difference in the bacteriological diagnosis rate among different specimen types (P = 0.035). The detection sensitivity of each method was the highest in pus, puncture tissue and bronchoalveolar lavage fluid, followed by sputum and pleural and abdominal effusions, and the lowest in cerebrospinal fluid. Conclusions The sensitivity of nucleic acid amplification technology for immunosuppressed KTRs is slightly reduced, but it should still be recommended as the preferred initial screening tool. The comprehensive strategy of "multi-site sampling (prior to obtaining high bacterial load specimens) combined with multiple method detection" may improve the diagnostic efficacy.

Background Work-related musculoskeletal disorders (WMSDs) are common across industries. Workers in the hairdressing industry are frequently exposed to repetitive upper-limb movements and awkward working postures, which may increase their susceptibility to WMSDs. However, evidence on WMSDs among workers in this industry remains limited. Objective This study aimed to investigate the prevalence of self-reported WMSDs symptoms and associated factors among workers in the hairdressing industry in Hongkou District, Shanghai. Methods A cross-sectional survey was conducted using quota sampling in all eight subdistricts of Hongkou District, Shanghai. Hair salons in each subdistrict were ranked by the number of hairdressers and shampoo workers, and the 6 salons with the largest staff size in each subdistrict were selected, resulting in 48 survey sites. Hairdressers and shampoo workers who met the inclusion and exclusion criteria were surveyed using the Chinese version of the Musculoskeletal Disorders Questionnaire. Information on their individual characteristics, musculoskeletal symptoms, and work-related factors was collected. Chi-square tests and Firth logistic regression were used to estimate the prevalence of self-reported WMSDs symptoms and identify associated factors. Results The overall prevalence of self-reported WMSDs symptoms among workers in the hairdressing industry was 45.72%. Among hairdressers, the most commonly affected body regions were the neck (31.51%), shoulders (25.57%), lower back (15.98%), and feet (14.61%). Among shampoo workers, the most commonly affected regions were the neck (24.17%), lower back (24.17%), shoulders (20.83%), and upper back (10.83%). Firth logistic regression showed that fatigue in the corresponding body region after several hours of continuous work was positively associated with self-reported WMSDs symptoms in both types of workers (OR=6.544-52.838). Regarding individual characteristics, shorter height (female≤160 cm, male≤165 cm) was positively associated with self-reported neck and foot symptoms among hairdressers, and self-reported shoulder symptoms among shampoo workers (OR=16.324, 15.084, 9.680). Female shampoo workers had higher odds of self-reported neck symptoms than male shampoo workers (OR=4.924). Smoking was positively associated with self-reported upper-back and lower-back symptoms among shampoo workers (OR=10.544, 4.381). Among work-related factors, maintaining a fixed neck posture for long periods among hairdressers and prolonged standing among shampoo workers were positively associated with self-reported neck symptoms (OR=6.123, 3.080). Among shampoo workers, frequent overtime and maintaining a twisted posture were positively associated with self-reported upper-back symptoms (OR=11.710, 35.306). Sufficient rest time was negatively associated with self-reported foot symptoms among hairdressers (OR=0.238). Conclusion The prevalence of self-reported WMSDs symptoms was relatively high among workers in the hairdressing industry. Work-related fatigue and awkward working postures were strongly associated with WMSDs symptoms. Measures such as optimizing work-rest schedule, improving ergonomic workplace design, implementing fatigue monitoring and management, promoting work-break exercises, and providing training on appropriate working posture may help reduce the burden of WMSDs symptoms in this occupational group.

ObjectiveTo establish a model that can quickly identify the aroma components in different parts of Nardostachys jatamansi， so as to provide a quality control basis for the market circulation and clinical use of N. jatamansi. MethodsHeadspace solid-phase microextraction-gas chromatography-mass spectrometry（HS-SPME-GC-MS） combined with Smart aroma database and National Institute of Standards and Technology（NIST） database were used to characterize the aroma components in different parts of N. jatamansi， and the aroma components were quantified according to relative response factor（RRF） and three internal standards， and the markers of aroma differences in different parts of N. jatamansi were identified by orthogonal partial least squares-discriminant analysis（OPLS-DA） and cluster thermal analysis based on variable importance in the projection（VIP） value >1 and P<0.01. The odor data of different parts of N. jatamansi were collected by Heracles Ⅱ Neo ultra-fast gas phase electronic nose， and the correlation between compound types of aroma components collected by the ultra-fast gas phase electronic nose and the detection results of HS-SPME-GC-MS was investigated by drawing odor fingerprints and odor response radargrams. Chromatographic peak information with distinguishing ability≥0.700 and peak area≥200 was selected as sensor data， and the rapid identification model of different parts of N. jatamansi was established by principal component analysis（PCA）， discriminant factor alysis（DFA）， soft independent modeling of class analogies（SIMCA） and statistical quality control analysis（SQCA）. ResultsThe HS-SPME-GC-MS results showed that there were 28 common components in the underground and aboveground parts of N. jatamansi， of which 22 could be quantified and 12 significantly different components were screened out. Among these 12 components， the contents of five components（ethyl isovalerate， 2-pentylfuran， benzyl alcohol， nonanal and glacial acetic acid，） in the aboveground part of N. jatamansi were significantly higher than those in the underground part（P<0.01）， the contents of β-ionone， patchouli alcohol， α-caryophyllene， linalyl butyrate， valencene， 1，8-cineole and p-cymene in the underground part of N. jatamansi were significantly higher than those in the aboveground part（P<0.01）. Heracles Ⅱ Neo electronic nose results showed that the PCA discrimination index of the underground and aboveground parts of N. jatamansi was 82， and the contribution rates of the principal component factors were 99.94% and 99.89% when 2 and 3 principal components were extracted， respectively. The contribution rate of the discriminant factor 1 of the DFA model constructed on the basis of PCA was 100%， the validation score of the SIMCA model for discrimination of the two parts was 99， and SQCA could clearly distinguish different parts of N. jatamansi. ConclusionHS-SPME-GC-MS can clarify the differential markers of underground and aboveground parts of N. jatamansi. The four analytical models provided by Heracles Ⅱ Neo electronic nose（PCA， DFA， SIMCA and SQCA） can realize the rapid identification of different parts of N. jatamansi. Combining the two results， it is speculated that terpenes and carboxylic acids may be the main factors contributing to the difference in aroma between the underground and aboveground parts of N. jatamansi.

Patients suffering from nerve injury often experience exacerbated pain responses and complain of memory deficits. The dorsal hippocampus (dHPC), a well-defined region responsible for learning and memory, displays maladaptive plasticity upon injury, which is assumed to underlie pain hypersensitivity and cognitive deficits. However, much attention has thus far been paid to intracellular mechanisms of plasticity rather than extracellular alterations that might trigger and facilitate intracellular changes. Emerging evidence has shown that nerve injury alters the microarchitecture of the extracellular matrix (ECM) and decreases ECM rigidity in the dHPC. Despite this, it remains elusive which element of the ECM in the dHPC is affected and how it contributes to neuropathic pain and comorbid cognitive deficits. Laminin, a key element of the ECM, consists of α-, β-, and γ-chains and has been implicated in several pathophysiological processes. Here, we showed that peripheral nerve injury downregulates laminin β1 (LAMB1) in the dHPC. Silencing of hippocampal LAMB1 exacerbates pain sensitivity and induces cognitive dysfunction. Further mechanistic analysis revealed that loss of hippocampal LAMB1 causes dysregulated Src/NR2A signaling cascades via interaction with integrin β1, leading to decreased Ca²⁺ levels in pyramidal neurons, which in turn orchestrates structural and functional plasticity and eventually results in exaggerated pain responses and cognitive deficits. In this study, we shed new light on the functional capability of hippocampal ECM LAMB1 in the modulation of neuropathic pain and comorbid cognitive deficits, and reveal a mechanism that conveys extracellular alterations to intracellular plasticity. Moreover, we identified hippocampal LAMB1/integrin β1 signaling as a potential therapeutic target for the treatment of neuropathic pain and related memory loss.
Animals ; Laminin/genetics* ; Hippocampus/metabolism* ; Neuralgia/metabolism* ; Cognitive Dysfunction/etiology* ; Male ; Peripheral Nerve Injuries/metabolism* ; Extracellular Matrix/metabolism* ; Integrin beta1/metabolism* ; Pyramidal Cells/metabolism* ; Signal Transduction

OBJECTIVE: To establish and validate a novel diabetic retinopathy (DR) risk-prediction model using a whole-exome sequencing (WES)-based machine learning (ML) method. METHODS: WES was performed to identify potential single nucleotide polymorphism (SNP) or mutation sites in a DR pedigree comprising 10 members. A prediction model was established and validated in a cohort of 420 type 2 diabetic patients based on both genetic and demographic features. The contribution of each feature was assessed using Shapley Additive explanation analysis. The efficacies of the models with and without SNP were compared. RESULTS: WES revealed that seven SNPs/mutations ( rs116911833 in TRIM7, 1997T>C in LRBA, 1643T>C in PRMT10, rs117858678 in C9orf152, rs201922794 in CLDN25, rs146694895 in SH3GLB2, and rs201407189 in FANCC) were associated with DR. Notably, the model including rs146694895 and rs201407189 achieved better performance in predicting DR (accuracy: 80.2%; sensitivity: 83.3%; specificity: 76.7%; area under the receiver operating characteristic curve [AUC]: 80.0%) than the model without these SNPs (accuracy: 79.4%; sensitivity: 80.3%; specificity: 78.3%; AUC: 79.3%). CONCLUSION Novel SNP sites associated with DR were identified in the DR pedigree. Inclusion of rs146694895 and rs201407189 significantly enhanced the performance of the ML-based DR prediction model.
Diabetic Retinopathy/diagnosis* ; Humans ; Machine Learning ; Male ; Female ; Polymorphism, Single Nucleotide ; Middle Aged ; Exome Sequencing ; Aged ; Adult ; Pedigree ; Diabetes Mellitus, Type 2/complications* ; Genetic Predisposition to Disease ; Mutation

Objective:Previous studies have shown that insomnia is closely related to erectile dysfunction(ED).However,the causal relationship between them is still unclear.Mendelian randomization(MR)provides a new method for studying the relationship between the two,and the theory of heart-kidney interaction in TCM provides a new idea for exploring the causal relationship between them.Methods:Based on the statistical data collected by genome-wide association studies(GWAS),the causal relationship be-tween insomnia and ED was discussed by MR.Inverse variance weighted(IVW)is the main analysis method,and weighted median(WME),simple mode(SM),weighted mode(WM)and MR Egger method were the supplementary analysis to evaluate the causal effect.MR-Egger intercept test,Cochran Q test and leave-one-out method were used in sensitivity analysis to verify the reliability of MR results.Results:Thirty-nine SNPs significantly related to insomnia were finally included for MR analysis.The results of IVW method in MR analysis showed that insomnia had a significant causal relationship with the increased risk of ED(OR=3.111,95％CI=1.566-6.181,P=1.193 × 10-3).The results obtained by MR-Egger method,WME method,WM method and SM method were consistent with IVW method in the direction of effect.The sensitivity results suggested that the results of this study were robust.Conclusion:Our study reveals the causal relationship between insomnia and ED,which provides a new basis for future clinical practice and prevention and treatment of ED.