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.

Objective To investigate the mechanism by which Senegenin(SEN)alleviates microglial inflammatory response through the nuclear factor erythroid 2-related factor 2(Nrf2)/NOD-like receptor protein 3(NLRP3)pathway.Methods BV2 mouse microglia cells were randomly divided into control group,model group,SEN group and MCC950 group.Cells in control group were not treated,and cells in model group were added with 1 μg/ml lipopolysaccharide(LPS);Cells in SEN group were added with 1 μg/ml LPS+4 μmol/L SEN,and cells in MCC950 group were added with 1 μg/ml LPS+10 μmol/L MCC950 for 24 hours.CCK-8 method was used to detect the effect of different concentrations of SEN on the viability of BV2 cells.Griess method was used to determine the release amount of nitric oxide(NO)in the supernatant.Real-time fluorescent quantitative PCR was used to determine the mRNA expression levels of NLRP3,lymphocyte apoptosis-associated spect-like protein containing a CARD(ASC),caspase-1,interleukin(IL)-1β and IL-18 mRNA.Immunofluorescence staining was used to detect the expression levels of ASC,IL-1β,Nrf2 and heme oxygenase-1(HO-1).Western blotting was used to detect the expression levels of NLRP3,caspase-1,ASC,IL-1β,IL-18,Nrf2,HO-1,nuclear factor kappa B(NF-κB)and inducible nitric oxide synthase(iNOS).Results The results of CCK-8 method showed that there was no significant difference in the viability of BV2 cells treated with 2～20 μmol/L SEN compared with control group(P>0.05).Compared with control group,the viability of BV2 cells in model group decreased significantly(P<0.05).Compared with model group,the viability of BV2 cells in 4 μmol/L SEN group was significantly restored(P<0.05).Compared with control group,the results of Griess method showed that the release amount of NO in cells of model group increased significantly(P<0.05);the results of real-time PCR showed that the expression levels of NLRP3,ASC,caspase-1,IL-1β and IL-18 mRNA in cells of model group increased significantly(P<0.05);the results of Western blotting showed that the protein expression levels of NLRP3,ASC,caspase-1,IL-1β and IL-18 proteins in cells of model group increased significantly(P<0.05),and the immunofluorescence staining results showed that the expression levels of iNOS and NF-κB protein in cells of model group increased,and the expression levels of Nrf2 and HO-1 decreased,with statistically significant differences(P<0.05).Compared with model group,the release amount of NO in cells of SEN group and MCC950 group decreased,and the expression levels of NLRP3,ASC,caspase-1,IL-1β and IL-18 mRNA and proteins decreased,with statistically significant differences(P<0.05);in the SEN group,the expression levels of iNOS and NF-κB decreased,and immunofluorescence staining showed that Nrf2 was translocated into the nucleus,and the expression levels of Nrf2 and HO-1 proteins increased significantly,with statistically significant differences(P<0.05).Conclusions SEN could alleviate the inflammatory response of mouse microglia cells induced by LPS and inhibit the activation and expression of NLRP3 inflammasome,with an effect comparable to that of the inflammasome inhibitor MCC950.The mechanism may be related to the regulation of the expression of upstream factors Nrf2 and HO-1.

Cellular death in the body can occur through different processes,including apoptosis,pyroptosis and necrotic apoptosis,etc.PANoptosis is a newly discovered form of inflammatory cell death in recent years.It can be triggered by various stimulating factors and integrates multiple components that can induce cell death to assemble into various types of macromolecular complexes-PANoptosome,which then mediates cell death.Given the impact of PANoptosis on the entire disease spectrum,promoting or inhibiting its occurrence process may prevent the development of various diseases.The review summarizes the research progress on the occurrence mechanism of PANoptosis and its role in some diseases,and explores the crosstalk among multiple programmed cell death pathways,aiming to provide new ideas for the treatment of related diseases.

Cosmetics may be an important source of human exposure to per-and polyfluoroalkyl substances(PFASs),posing risks to human health.In this study,a nontarget screening method for PFASs in cosmetics was developed using ultra-high performance liquid chromatography-high-resolution mass spectrometry(UHPLC-Q-Orbitrap HRMS)based on the Kendrick mass defect(KMD).The sample was extracted by ultrasonic assisted extraction prior to being analyzed by UHPLC-Q-Orbitrap HRMS.Acquisition of HRMS data was achieved in both full scan and data-dependent(Full MS/dd MS2)mode.The data collected by HRMS were imported into an in-lab built R script for processing.Samples retained the mass spectra peaks with KMD values in the range of 0.85-1 or 0-0.15 for in-and out-of-library matching;when KMD deviation(δKMD)<0.001 and CF 2 mass error(δMS)<15 ppm,it was considered as a potential PFASs homologues.According to matches of parent ions(MS),fragment ions(MS2)and retention time(RT)with the in-house built PFASs database,the screened and identified potential PFASs were categorized to 5 confidence levels(CL1-CL5).A total of 15 kinds of PFASs homologues with confidence level of CL3 and above were screened from 13 cosmetics products and 8 cosmetic raw materials,including perfluoroalkyl alcohol,hydroperfluoroalkyl sulfonic acid,chloroperfluoroalkyl sulfonic acid,etc.with concentrations ranging from 1.9 ng/g to 98.1 ng/g.The nontarget screening method could be used to screen and identify PFASs homologues feasibly and therefore provided data basis for management and control of PFASs addition in cosmetics.