ObjectiveTo investigate the contamination status of ochratoxin A (OTA) in commercially available food products in Shanghai, and to assess OTA exposure levels and the associated non-carcinogenic and carcinogenic risks among pregnant women by integrating dietary consumption data of this population. MethodsThe levels of OTA contamination in 1 520 food samples collected in Shanghai from 2022 to 2023 were determined using liquid chromatography-tandem mass spectrometry. An exposure assessment model was developed based on the dietary consumption levels of pregnant women from the 2016‒2017 Shanghai Pregnant Women Dietary Monitoring Survey to calculate the estimated daily intake (EDI) of OTA, the margin of exposure for non-carcinogenic toxicity (MOE₁), and the margin of exposure for carcinogenic toxicity (MOE₂). An MOE₁ greater than 200 and an MOE₂ greater than 10 000 indicate that the non-carcinogenic toxicity and carcinogenic toxicity resulting from exposure are negligible, respectively. For samples with OTA contamination levels below the limit of detection (LOD), which accounted for more than 80% of the samples, the OTA levels were assigned values of 0 and LOD, respectively, for subsequent calculations. ResultsThe detection rates of OTA in cereals, nuts, dried fruits, and alcohol samples collected in 2022 were 2.03%, 0, 0, and 0, respectively. The OTA detection rates in cereals, nuts, dried fruits, beans, and alcohol samples collected in 2023 were 2.50%, 0.39%, 2.47%, 1.67%, and 13.33%, respectively. For pregnant women in Shanghai in 2022, simulation results indicated that when assigning a value of 0 and the LOD, theP₅₀ values of EDI for dietary OTA exposure were 0.05 and 0.72 ng·(kg·d)^-1, respectively, and the P₉₅ values of EDI for dietary OTA exposure were 0.25 and 2.40 ng·(kg·d)^-1, respectively. For pregnant women in Shanghai in 2023, the P₅₀ values of EDI for dietary OTA exposure were 0.04 and 1.00 ng·(kg·d)^-1, respectively, and the P₉₅ values of EDI for dietary OTA exposure were 0.23 and 2.67 ng·(kg·d)^-1, respectively, both substantially below the tolerable daily intake (TDI) for OTA [17 ng·(kg·d)^-1]. The EDI for dietary OTA exposure in 100.0% of Shanghai pregnant women was lower than the TDI, indicating an overall low level of dietary OTA exposure among this population. For 100.0% of pregnant women, the MOE₁ for dietary OTA exposure exceeded 200. When assigned a value of 0, the MOE₂ for 100.0% of pregnant women in both 2022 and 2023 exceeded10 000. When assigned the LOD value, 72.3% and 81.8% of pregnant women in 2022 and 2023, respectively, had an MOE₂ exceeding 10 000. ConclusionFrom 2022 to 2023, samples of cereals, nuts, dried fruits, beans, and alcohol sold in Shanghai exhibited varying degrees of OTA contamination. The overall EDI of OTA exposure among pregnant women in Shanghai remained at a low level. The non-carcinogenic and carcinogenic risks associated with OTA exposure were generally low and at controllable levels.

OBJECTIVES: To study the clinical and genetic characteristics of children with congenital adrenal hyperplasia (CAH). METHODS: Clinical data, laboratory findings, and genetic test results of 63 children diagnosed with CAH at Henan Children's Hospital from January 2017 to December 2024 were retrospectively reviewed. RESULTS: Of the 63 patients, the mean age at the first visit was (21 ± 14) days; 29 (46%) were of male sex and 34 (54%) were of female sex. The predominant clinical manifestations were poor weight gain or weight loss (92%, 58/63), poor feeding (84%, 53/63), skin hyperpigmentation (83%, 52/63), and female external genital anomalies (100%, 34/34). Laboratory abnormalities included hyponatremia (87%, 55/63), hyperkalemia (68%, 43/63), metabolic acidosis (68%, 43/63), and markedly elevated 17-hydroxyprogesterone (92%, 58/63), testosterone (89%, 56/63), and adrenocorticotropic hormone (81%, 51/63). Among 49 patients who underwent genetic testing, CYP21A2 variants were identified in 90% (44/49), with c.293-13A/C>G (33%, 30/91) and large deletions/gene conversions (29%, 26/91) being the most frequent; STAR (8%, 4/49) and HSD3B2 (2%, 1/49) variants were also detected. Following hormone replacement therapy, electrolyte disturbances were corrected in 57 cases, with significant reductions in 17-hydroxyprogesterone, adrenocorticotropic hormone, and testosterone levels (P<0.001). CONCLUSIONS CAH presenting in neonates or young infants is characterized by electrolyte imbalance, external genital anomalies, and abnormal hormone levels. Genetic testing enables definitive subtype classification; in CYP21A2-related CAH, c.293-13A/C>G is a hotspot variant. These findings underscore the clinical value of genetic testing for early diagnosis and genetic counseling in CAH. Citation:Chinese Journal of Contemporary Pediatrics, 2025, 27(11): 1367-1372.
Humans ; Adrenal Hyperplasia, Congenital/diagnosis* ; Male ; Female ; Retrospective Studies ; Infant ; Infant, Newborn

ObjectiveTo analyze the current status of per⁃ and polyfluoroalkyl substances (PFASs) contamination in animal⁃derived foods in Shanghai and to assess the risk of dietary exposure among local residents, thereby providing a scientific basis for future dietary safety and risk management. MethodsA total of 300 commonly consumed animal⁃derived food samples were collected in Shanghai in 2023 and tested for 17 types of PFASs. Based on local dietary consumption data, the weekly exposure intake (EWI) of four representative PFASs was calculated using the exposure assessment model recommended by Food and Agriculture Organization (FAO) and the World Health Organization (WHO). The calculated EWI was compared with the tolerable weekly intake (TWI) set by the European Food Safety Authority (EFSA) to assess the health risks associated with dietary intake of PFASs. ResultsThe predominant PFASs detected in aquatic products were PFUnDA (98.67%), PFTrDA (98.00%), and PFOS (95.33%), with PFUnDA having the highest median mass fraction (0.378 ng·g^-1). In meat samples, PFBA (54.17%), PFOA (20.83%), and PFOS (18.33%) were mainly detected, with PFBA showing the highest median mass fraction (0.027 ng·g^-1). In egg samples, the most frequently detected compounds were PFBA (90.00%), PFOA (63.33%), PFOS (50.00%), PFDA (50.00%), PFNA (50.00%), and PFHxDA (50.00%), with PFBA again having the highest median mass fraction (0.068 ng·g^-1). The median mass fractions of the four key PFASs (PFOA, PFOS, PFNA, and PFHxS) in aquatic products were significantly higher than those in meat products and eggs (all P<0.001), and there were statistically significant differences in the median mass fractions of the four PFASs among different types of aquatic products (all P<0.05). The EWI of the four PFASs was 2.970 ng·kg^-1 for the mean consumption group and 6.676 ng·kg^-1 for the high consumption group (P₉₅), with females having higher EWI values than males. ConclusionPFOS, PFOA, and PFBA are the major PFAS contaminants in animal⁃derived foods consumed by Shanghai residents, with aquatic products having the highest levels of contamination. The EWIs of four PFASs among Shanghai residents did not exceed EFSA’s recommended TWI at average consumption levels. However, the EWIs for high⁃consumption populations exceeded the EFSA⁃recommended TWI, posing potential health risks to these populations. Aquatic products contribute the most to dietary PFAS exposure. It is necessary to strengthen PFAS monitoring in aquatic products and conduct targeted risk assessments for high⁃consumption groups.

Background Deoxynivalenol (DON), a priority contaminant for food safety risk monitoring, is produced by Fusarium spp. infesting crops, and its common derivatives are 3-acetyl-DON (3A-DON) and 15-acetyl-DON (15A-DON), which have been shown to possess gastrointestinal toxicity, immunotoxicity, reproductive toxicity, and cytotoxicity. Due to the stable physicochemical properties of the DON family of toxins (DONs), they cannot be effectively removed during food processing, thus following the food chain, entering the human body, and posing health risks. Objective To understand the contamination status of DONs in commercial foods (cereals and bakery products) in Shanghai in 2022–2023, and to assess the exposure risk of DONs in pregnant women by combining their dietary consumption data. Methods Liquid chromatography tandem mass spectrometry (LC-MS/MS) was used to determine the contamination level of DONs in 1 100 food samples (cereals and baked goods) collected in 2022 and 944 samples collected in 2023 from Shanghai. The dietary monitoring data of pregnant women in Shanghai from 2016 to 2017 were adopted. The monitoring employed the food frequency questionnaire distributed among pregnant women through a combination of online telephone enquiry and offline on-site face-to-face survey to estimate their food consumption levels. An exposure assessment model was established to calculate the exposure level to DONs, and the probability distribution of the DONs exposure level in the pregnant women group in Shanghai was obtained by applying @Risk 7.5 software and simulating the calculation according to the Monte Carlo principle. With reference to the tolerable daily intake (TDI) of DONs [1.00 µg·(kg·d)⁻¹] proposed by the Joint FAO/WHO Expert Committee on Food Additives, the risk of exposure to DONs from commercial cereals and bakery products in pregnant women in Shanghai was assessed. Results DONs were detected in cereal and bakery samples collected in 2022 and 2023 with different levels of contamination. The level of DONs in cereal foods in 2023 (mean: 36.33 µg·kg⁻¹) decreased compared to 2022 (mean: 23.64 µg·kg⁻¹). However, the positive rate (71.67%) and level (mean: 51.22 µg·kg⁻¹) of DONs in bakery products increased significantly compared with 2022 (positive rate: 10.00%, mean: 24.39 µg·kg⁻¹). The mean consumption of cereals in 783 pregnant women was 222.48 g·d⁻¹ and the mean consumption of bakery products was 36.07 g·d⁻¹, and there was no statistically significant difference in the intake of all types of cereals and bakery products across the early, middle, and late stages of pregnancy. The modelled intakes of DONs via commercial cereals and bakery products for pregnant women in Shanghai were calculated to be 0.20 and 0.57 µg·(kg·d)⁻¹ in 2022 for the mean level and the 95th percentile level, respectively, and 0.16 µg·(kg·d)⁻¹ and 0.35 µg·(kg·d)⁻¹ in 2023, respectively. The results of the health risk assessment showed that pregnant women in Shanghai had 2.6% and 1.4% probability of exposure to DONs from cereal consumption in 2022 and 2023, respectively. Conclusion The risk of exposure of pregnant women in Shanghai to DONs via commercial cereals and bakery products is relatively low (1.4%-2.6%). However, considering the physical sensitivity of pregnant women, they should avoid consuming moldy grains and appropriately reduce intake of bakery products.

AIM To investigate the effects of Congrong San on neuronal apoptosis and the Bax/Bcl-2/Caspase3 signaling pathway in a rat model of Alzheimer's disease(AD).METHODS A total of 60 2-month-old SD male rats were randomly divided into the blank group,the model group,the memantine hydrochloride group(0.025 g/kg)and low-dose,medium-dose and high-dose Congrong San groups(4.62,9.24,18.48 g/kg).All groups except the control group received stereotactic intracerebral injection of Aβ1-42 to establish AD models.Following the successful modeling,each group received its corresponding intragastric administration once daily for 28 consecutive days.After the administration,the rats had their learning and memory ability detected by the morris water maze test;their hippocampal neuronal morphology observed with HE and Nissl staining;their hippocampal neuronal apoptosis observed with TUNEL staining;and their hippocampal expressions of amyloid precursor protein(APP),β-site APP-cleaving enzyme 1(BACE1),and apoptosis-related proteins Bax,Bcl-2 and Caspase3 detected with Western blot assay.RESULTS Compared with the model group,the groups intervened with medium-dose and high-dose Congrong San exhibited improved learning and memory performance,alleviated hippocampal neuronal damage,increased Nissl body count(P＜0.01),reduced hippocampal apoptosis rate(P＜0.05,P＜0.01),decreased protein expressions of APP,BACE1,Bax and cleaved-Caspase3/Caspase3 ratio(P＜0.05,P＜0.01),and elevated Bcl-2 expression(P＜0.01).CONCLUSION Congrong San mitigates cognitive impairment,hippocampal neuronal damage,and apoptosis in AD rats,probably through inhibition of the Bax/Bcl-2/Caspase3 signaling pathway activation.

Rapid and ultrasensitive detection of pathogen-associated biomarkers is vital for the early diagnosis and therapy of bacterial infections.Herein,we developed a close-packed and ordered Au@AgPt array coupled with a cascade triggering strategy for surface-enhanced Raman scattering(SERS)and colorimetric identification of the Staphylococcus aureus biomarker micrococcal nuclease(MNase)in serum samples.The trimetallic Au@AgPt nanozymes can catalyze the oxidation of 3,3',5,5'-tetramethylbenzidine(TMB)molecules to SERS-enhanced oxidized TMB(oxTMB),accompanied by the color change from colorless to blue.In the presence of S.aureus,the secreted MNase preferentially cut the nucleobase AT-rich regions of DNA sequences on magnetic beads(MBs)to release alkaline phosphatase(ALP),which subsequently mediated the oxTMB reduction for inducing the colorimetric/SERS signal fade away.Using this"on-to-off"triggering strategy,the target S.aureus can be recorded in a wide linear range with a limit of detection of 38 CFU/mL in the colorimetric mode and 6 CFU/mL in the SERS mode.Meanwhile,the MNase-mediated strategy characterized by high specificity and sensitivity successfully discriminated between patients with sepsis(n=7)and healthy participants(n=3),as well as monitored the prog-nostic progression of the disease(n=2).Overall,benefiting from highly active and dense"hot spot"substrate,MNase-mediated cascade response strategy,and colorimetric/SERS dual-signal output,this methodology will offer a promising avenue for the early diagnosis of S.aureus infection.

Objective To investigate the effects of Sini Powder on serum hypothalamic-pituitary-adrenal(HPA)axis-related hormones and inflammatory factors in liver cancer mice with comorbid depression,and to evaluate its effect on depressive behavior and tumor proliferation activity.Methods Forty-eight specific pathogen-free female C57BL/6 mice were randomly assigned to either a blank(n=8)or model group(n=40).The modeling group was subjected to chronic unpredictable mild stress(CUMS)for six weeks.Both groups underwent orthotopically transplanted liver tumor surgery at the end of the fourth week of CUMS treatment.At the end of the sixth week of CUMS treatment,color Doppler ultrasonography was used to observe tumor formation in the orthotopic transplantation liver tumors,and the tail suspension test was used to assess depressive behavior.Non-tumor-bearing and deceased mice were excluded.The remaining model group mice were stratified by tail suspension immobility time and randomly assigned to the following groups:model group(distilled water),Fluoxetine group(5.0 mg/kg),and Sini Powder low-dose,medium-dose,and high-dose groups(5.2,10.4,and 20.8 g/kg,respectively),with six mice per group.The treatments were administered once daily for 21 consecutive days.After treatment,depressive behaviors were assessed using the open field,tail suspension,and forced swimming tests.The proliferation status of the orthotopic liver transplantation tumor was evaluated by measuring the size of the tumor,observing pathological changes in the tumor tissue through hematoxylin and eosin staining,and detecting the positive cell rate of proliferating cell nuclear antigen(Ki-67)in the tumor tissue using immunohistochemistry.The levels of HPA axis-related hormones in serum,such as corticotropin-releasing hormone(CRH),adrenocorticotropic hormone(ACTH),corticosterone(CORT),as well as inflammatory factors such as tumor necrosis factor-α(TNF-α),interleukin-1β(IL-1β),and interleukin-6(IL-6)were measured using an enzyme-linked immunosorbent assay.Western blotting was used to assess mitogen-activated protein kinase(MAPKs)phosphorylation and the expression of nuclear factor-κB(NF-κB),NOD-like receptor family pyrin domain-containing receptor 3(NLRP3),and cysteine aspartic protease-1(Caspase-1)in orthotopic tumors.Results Compared with the blank group,the model group showed reduced total distance traveled in open field test,prolonged immobility times in the tail suspension and forced swimming tests(P<0.05,P<0.01),indicating successful establishment of the liver cancer with comorbid depression mice model.Also,the model group showed increased orthotopic tumor volume(P<0.01),and elevated serum CRH,ACTH,CORT,TNF-α,IL-1β,and IL-6 levels(P<0.01).The phosphorylation of MAPKs in tumor tissues was suppressed(P<0.01),while NF-κB,NLRP3,and Caspase-1 expression levels were downregulated(P<0.01).Compared with the model group,Sini Powder medium-and high-dose groups exhibited increased total distance traveled in the open field test(P<0.05),reduced forced swimming test and prolonged total distance traveled in open field test(P<0.01),while Sini Powder high-dose group showed reduced immobility times in the tail suspension test(P<0.05).Also,Sini Powder low-dose,medium-dose,and high-dose groups showed slower tumor growth,histological changes,including vacuolization and necrosis,decreased Ki-67 positive cell rate(P<0.01),and reduced serum CRH,ACTH,CORT,TNF-α,IL-1β,and IL-6 levels(P<0.05).Additionally,the phosphorylation of MAPKs in tumor tissues was suppressed(P<0.01),and NF-κB,NLRP3,and caspase-1 expression levels were downregulated(P<0.01).Conclusion Sini Powder may alleviate depressive behaviors and suppress tumor proliferation activity in liver cancer mice with comorbid depression by modulating MAPKs activation,inhibiting NF-κB,NLRP3,and Caspase-1 expressions,and reducing serum inflammatory factors and HPA axis-related hormones levels.

Objective To analyze the correlation between lesion volume,lesion mass,and maximum lesion diameter in the assessment of advanced hepatocarcinoma with lung metastasis,and to evaluate the application value of total volume response and total mass response of lung metastatic lesions in efficacy assessment.Methods A retrospective analysis was conducted on the CT imaging data of 20 patients clinically confirmed with hepatocarcinoma and lung metastases,followed by subsequent follow-up to monitor their survival outcomes.Volume measurement software was used to measure the volume of lesions before and after treatment.We recored lesion diameter,volume measurements and CT values,calculated the mass of the lesions.The correlation between lesion volume,mass and diameter was analyzed,as well as the correlation between the change rates of volume,mass and lesion diameter.Additionally,the total volume and total mass of all lesions were calculated.The correlation between the change rates of total volume/total mass and the change rate of pulmonary lesion diameter under the RECIST 1.1 criteria,as well as the correlation with changes in patients'tumor markers,were analyzed.Furthermore,the overall volume response and overall mass response of lesions were evaluated based on changes in total volume and total mass,and their consistencies with the RECIST 1.1 criteria for efficacy evaluation were analyzed.Finally,univariate Cox regression analysis was performed to explore the association between these variables and patient survival outcomes.Results There was strong correlation between lesion volume,mass and tumor diameter(r=0.771,0.775),between the rate of change in mass and the rate of change in lesion diameter(r=0.846),and between the rates of change in total volume/total mass and the rate of change in pulmonary lesion diameter under the RECIST 1.1 criteria(r=0.800,0.896).The correlation between the rates of change in total volume/total mass and patients'tumor markers was not statistically significant.There was moderate correlation between the rate of change in volume and the rate of change in lesion diameter(r=0.692).The evaluation results of total volume response and total mass response for pulmonary lesions in advanced hepatocarcinoma with lung metastasis were generally consistent with the RECIST 1.1 criteria(Kappa=0.486,0.426).Univariate Cox regression analysis revealed that total lesion volume(P=0.047)and total lesion mass(P=0.049)were independent prognostic factors for survival outcomes.Conclusion Lesion volume,mass,and diameter,as well as their respective change rates,were found to be interrelated.Furthermore,total lesion volume and total lesion mass were identified as independent prognostic factors for survival outcomes.The total volume response and total mass response are promising evaluation methods in evaluating the efficacy of lung metastasis of hepatocarcinoma,which are different from the RECIST 1.1 evaluation criteria.

Background and purpose:The resistance of pancreatic cancer to albumin-bound paclitaxel affects the therapeutic effect and prognosis.Signal transducer and activator of transcription 3(STAT3)is one of the important molecules regulating the chemotherapy sensitivity of cancer cells.The liposome BP1003 targeting the antisense oligonucleotide of STAT3 mRNA can inhibit the expression of STAT3 and increase the chemotherapy sensitivity.However,the effect of BP1003 on the sensitivity of pancreatic cancer cells to albumin-bound paclitaxel remains unclear.The purpose of this study was to investigate the effects of liposome binding antisense oligonucleotide BP1003 on albumin-bound paclitaxel sensitivity in pancreatic cancer cells by inhibiting STAT3.Methods:Pancreatic cancer cell lines PANC-1 and ASPC-1 were cultured.They were divided into control group(without drugs),BP1003 group(200 μg/mL BP1003 intervention),different concentrations of albumin-bound paclitaxel group(5,10,20 nmol/L albumin-bound paclitaxel intervention),BP1003+different concentrations of albumin-bound paclitaxel group(200 μg/mL BP1003 combined with 5,10,20 nmol/L albumin-bound paclitaxel intervention).The proliferation viability,apoptotic rate and the protein expression levels of STAT3,STAT4,STAT6,Bcl-2,Bax and c-Myc were detected.The transplanted tumor model was established by subcutaneous injection of PANC-1 and ASPC-1 cell suspension in nude mice,which were divided into control group(normal saline intervention),BP1003 group(25 mg/kg BP1003 intervention,once every 2 weeks)and albumin-bound paclitaxel group(10 mg/kg albumin-bound paclitaxel,once a week),BP1003+albumin-bound paclitaxel group(25 mg/kg BP1003 intervention,once every 2 weeks combined with 10 mg/kg albumin-bound paclitaxel,once a week).Four weeks later,the graft volume and mass were measured,and the protein expression levels of STAT3,Bcl-2,Bax and c-Myc were detected.Results:The apoptotic rate and the protein expression levels of Bax of PANC-1 and ASPC-1 cells in BP1003 group and albumin-bound paclitaxel group were higher than those in the control group,while the proliferation viability and protein expression levels of STAT3,Bcl-2 and c-Myc were lower than those in control group(P＜0.05).There was no significant difference in the expression levels of STAT4 and STAT6 in PANC-1 and ASPC-1 cells between BP1003 group and the control group(P＞0.05).The apoptotic rate and the protein expression levels of Bax of PANC-1 and ASPC-1 cells in BP1003+different concentrations of albumin-bound paclitaxel groups were higher than those in different concentrations of albumin-bound paclitaxel groups,and the proliferation viability and protein expression levels of STAT3,Bcl-2 and c-Myc were lower than those in different concentrations of albumin-bound paclitaxel groups(P＜0.05).The volume and mass of transplanted tumor and the protein expression levels of STAT3,Bcl-2 and c-Myc of nude mice in BP1003 group,albumin-bound paclitaxel group and BP1003+albumin-bound paclitaxel group were all lower compared with the control group,the protein expression level of Bax was higher compared with the control group(P＜0.05),and the above changes in BP1003+albumin-bound paclitaxel group were more significant compared with BP1003 and albumin-bound paclitaxel group.Conclusion:BP1003 increases the sensitivity of pancreatic cancer cells to albumin-bound paclitaxel by inhibiting the expression of STAT3.

OBJECTIVE: To identify prognostic genes associated with lysosome-dependent cell death (LDCD) in patients with gastric cancer (GC). METHODS: Differentially expressed genes (DEGs) were identified using The Cancer Genome Atlas - Stomach Adenocarcinoma. Weighted gene co-expression network analysis was performed to identify the key module genes associated with LDCD score. Candidate genes were identified by DEGs and key module genes. Univariate Cox regression analysis, and least absolute shrinkage and selection operator regression and multivariate Cox regression analyses were performed for the selection of prognostic genes, and risk module was established. Subsequently, key cells were identified in the single-cell dataset (GSE183904), and prognostic gene expression was analyzed. Cell proliferation and migration were assessed using the Cell Counting Kit-8 assay and the wound healing assay. RESULTS: A total of 4,465 DEGs, 95 candidate genes, and 4 prognostic genes, including C19orf59, BATF2, TNFAIP2, and TNFSF18, were identified in the analysis. Receiver operating characteristic curves indicated the excellent predictive power of the risk model. Three key cell types (B cells, chief cells, and endothelial/pericyte cells) were identified in the GSE183904 dataset. C19orf59 and TNFAIP2 exhibited predominant expression in macrophage species, whereas TNFAIP2 evolved over time in endothelial/pericyte cells and chief cells. Functional experiments confirmed that interfering with C19orf59 inhibited proliferation and migration in GC cells. CONCLUSION C19orf59, BATF2, TNFAIP2, and TNFSF18 are prognostic genes associated with LDCD in GC. Furthermore, the risk model established in this study showed robust predictive power.
Stomach Neoplasms/pathology* ; Humans ; Prognosis ; Lysosomes/physiology* ; RNA-Seq ; Cell Death ; Single-Cell Analysis ; Gene Expression Regulation, Neoplastic ; Cell Proliferation ; Single-Cell Gene Expression Analysis