OBJECTIVE: To assess the association of microribonucleic acid (miRNA) gene polymorphisms with the risk and clinicopathological characteristics of Crohn's disease (CD) and the influence of miRNA gene variants on the response to ustekinumab (UST) treatment among CD patients. METHODS: From January 2018 to February 2025, 312 patients diagnosed with CD and 527 gender- and age-matched normal controls were selected as the study subjects at the Department of Gastroenterology of the Second Affiliated Hospital of Wenzhou Medical University. Genotypes of miR-155 (rs767649), miR-21 (rs13137), miR-124 (rs531564) and miR-146a (rs57095329, rs2431697) were determined with multiplex polymerase chain reaction-ligase detection reaction (PCR-LDR) technique. The patients were divided into different subgroups according to the Montreal Classification Criteria for CD. Harvey-Bradshaw index (HBI) and simplified endoscopic score for CD were respectively applied to assess the clinical and endoscopic disease activity of CD. Unconditional logistic regression model was employed to analyze the distribution of miRNA gene polymorphisms between the two groups, as well as their influence on the clinicopathological characteristics of CD patients. Among them, 185 CD patients received first-line UST treatment, with the first sufficient dose of UST (6 mg/kg) administered intravenously. Based on the changes in HBI at week 8, the response of patients to UST treatment was evaluated. Unconditional logistic regression model was employed to analyze the distribution of miRNA gene polymorphisms between clinically responsive group (the decline of HBI ≥ 3 scores compared to week 0) and non-responsive group. All of the P values were adjusted by Bonferroni correction. This study has been approved by the Medical Ethics Committee of the Second Affiliated Hospital of Wenzhou Medical University (Ethics No.: 2025-K-12-01). RESULTS: No significant difference was found in the distribution of miRNA gene polymorphisms between the two groups (all P > 0.05). The variant genotype (TC+CC) of rs2431697 was more common among patients with terminal ileal-type and ileocolic-type CD than those with the colonic-type CD (OR = 4.98, 95%CI: 1.49~16.68, P = 0.009, adjusted P = 0.045). However, the opposite conclusion was drawn for the homozygous variant genotype (TT) of rs13137 and variant genotype (GC+CC) of rs531564 (OR = 0.37, 95%CI: 0.18~0.76, P = 0.007, adjusted P = 0.035; OR = 0.36, 95%CI: 0.18~0.73, P = 0.004, adjusted P = 0.020). Compared to patients with non-stricturing and penetrating CD, the variant genotype (AG+GG) and variant allele (G) of rs57095329 were more common in those with stricturing and penetrating CD (OR = 4.06, 95%CI: 2.46~6.71, P < 0.001, adjusted P < 0.005; OR = 3.12, 95%CI: 2.06~4.73, P < 0.001, adjusted P < 0.005). However, the frequencies of variant genotype (AT+TT) and variant allele (T) of rs13137 were lower among patients with stricturing and penetrating CD than in those without (OR = 0.25, 95%CI: 0.15~0.41, P < 0.001, adjusted P < 0.005; OR = 0.45, 95%CI: 0.33~0.63, P < 0.001, adjusted P < 0.005). Additionally, the variant genotype (AG+GG) and variant allele (G) of rs57095329 were more common among those with moderately to severely endoscopic activity than those with mildly endoscopic activity (OR = 2.01, 95%CI: 1.19~3.42, P = 0.009, adjusted P = 0.045; OR = 2.04, 95%CI: 1.28~3.25, P = 0.003, adjusted P = 0.015). In total 117 cases had shown clinical response by week 8, while 68 cases showed no response. Compared with t he clinically non-responsive group, the variant genotype (TC+CC) and variant allele (C) of rs2431697 were more common in the clinically responsive group (OR = 3.86, 95%CI: 1.80~8.32, P = 0.001, adjusted P = 0.005; OR = 2.60, 95%CI: 1.34~5.06, P = 0.005, adjusted P = 0.025). However, the variant genotype (TA+AA) of rs767649 was less frequent in the clinically responsive group than the non-responsive group (OR = 0.40, 95%CI: 0.21~0.74, P = 0.004, adjusted P = 0.020). The same conclusion was drawn for the variant genotype (AT+TT) and variant allele (T) of rs13137 when the clinically responsive group was compared with the non-responsive group (OR = 0.30, 95%CI: 0.14~0.63, P = 0.002, adjusted P = 0.010; OR = 0.54, 95%CI: 0.35~0.82, P = 0.005, adjusted P = 0.025). CONCLUSION Genetic polymorphisms of miRNAs are not associated with the risk of developing CD. The miR-146a (rs57095329) variant may increase the endoscopic activity of CD and the risk for stenosis or penetration. However, the miR-146a (rs2431697) variant may increase the risk of ileal involvement. The miR-21 (rs13137) variant may reduce the risk of ileal involvement and the risk of stenosis or penetration. The miR-124 (rs531564) variant may reduce the risk of ileal involvement. Among patients receiving UST treatment, the miR-146a (rs2431697) variant may increase the clinical response by week 8. However, both the miR-155 (rs767649) and miR-21 (rs13137) variants may decrease the clinical response by week 8.
Humans ; MicroRNAs/genetics* ; Crohn Disease/pathology* ; Male ; Female ; Adult ; Polymorphism, Single Nucleotide ; Middle Aged ; Asian People/genetics* ; Genetic Predisposition to Disease ; Genotype ; Young Adult ; Case-Control Studies ; Adolescent ; East Asian People

OBJECTIVES: To investigate the effects of methylenetetrahydrofolate reductase (MTHFR) rs1801133 and γ-glutamyl hydrolase (GGH) rs11545078 gene polymorphisms on plasma concentrations and toxicity following high-dose methotrexate (MTX) therapy in children with acute lymphoblastic leukemia (ALL). METHODS: Children with ALL treated at the Xuzhou Children's Hospital of Xuzhou Medical University from January 2021 to April 2024 were selected for this study. Genotypes of MTHFR rs1801133 and GGH rs11545078 were determined using multiplex polymerase chain reaction. MTX plasma concentrations were measured by enzyme-multiplied immunoassay technique, and toxicity was graded according to the Common Terminology Criteria for Adverse Events version 5.0. The relationships between MTHFR rs1801133 and GGH rs11545078 genotypes and both MTX plasma concentrations and associated toxicities were analyzed. RESULTS: In the low-risk ALL group, the MTHFR rs1801133 genotype was associated with increased MTX plasma concentrations at 72 hours (P<0.05). In the intermediate- to high-risk group, the MTHFR rs1801133 genotype was associated with increased MTX plasma concentrations at 48 hours (P<0.05), and the GGH rs11545078 genotype was associated with increased MTX plasma concentrations at 48 hours (P<0.05). In the intermediate- to high-risk group, the MTHFR rs1801133 genotype was associated with the occurrence of reduced hemoglobin (P<0.05), and the GGH rs11545078 genotype was associated with the occurrence of thrombocytopenia (P<0.05). CONCLUSIONS Detection of MTHFR rs1801133 and GGH rs11545078 genotypes can be used to predict increased MTX plasma concentrations and the occurrence of toxic reactions in high-dose MTX treatment of ALL, enabling timely interventions to enhance safety.
Humans ; Methotrexate/toxicity* ; Methylenetetrahydrofolate Reductase (NADPH2)/genetics* ; Precursor Cell Lymphoblastic Leukemia-Lymphoma/blood* ; Male ; Female ; Child ; Child, Preschool ; gamma-Glutamyl Hydrolase/genetics* ; Antimetabolites, Antineoplastic/adverse effects* ; Infant ; Polymorphism, Genetic ; Adolescent ; Genotype ; Polymorphism, Single Nucleotide

OBJECTIVE: To screen the genetic risk factors related to severe hematology adverse effects (AEs) in patients with chronic myeloid leukemia (CML) treated with imatinib (IM), and explore the correlation of single nucleotide polymorphisms (SNPs) in IM drug metabolism and transport pathway gene polymorphism with the risk of severe hematology AEs. METHODS: 172 newly diagnosed Chinese Han patients in CML chronic phase (CML-CP) treated with IM were included and divided into severe hematology AEs group and non-severe hematology AEs group. The demographic characteristics and laboratory test results were compared between the two groups. 11 gene SNP sites in the included subjects were genotyped using SNaPshot multiplex SNPs technique. RESULTS: Compared with non-severe hematology AEs group, the severe hematology AEs group had higher white blood cell (WBC) and EOS% (both P < 0.05), but lower hemoglobin (Hb) and hematocrit (HCT) (both P < 0.01). For rs1045642 of ABCB1 gene, there were significant differences in the distribution of allele frequency and genotype frequency of this loci between severe hematology AEs group and non-severe hematology AEs group (both P < 0.05). Carriers of rs1045642 mutation allele A had an increased risk of severe hematology AEs (OR =2.09, 95% CI : 1.24-3.55, P =0.005). There was a significant difference in the distribution of NR1I2 gene rs3814055 genotype between severe hematology AEs group and non-severe hematology AEs group (P < 0.05). The additive model and recessive model of ABCB1 gene rs1045642 and the recessive model of NR1I2 gene rs3814055 were associated with the increased risk of severe hematology AEs (OR =2.14, 3.28, 5.54, all P < 0.05). CONCLUSION Peripheral blood WBC, EOS%, Hb and HCT in patients with newly diagnosed CML-CP are all related to the risk of severe hematology AEs. ABCB1 gene rs1045642 and NR1I2 gene rs3814055 related to the metabolism and transport pathway of IM are associated with severe hematology AEs after IM treatment in CML-CP patients, and they may be potential molecular markers to predict the risk of severe hematology AEs of CML patients treated by IM.
Humans ; Imatinib Mesylate ; Leukemia, Myelogenous, Chronic, BCR-ABL Positive/genetics* ; Polymorphism, Single Nucleotide ; Genotype ; ATP Binding Cassette Transporter, Subfamily B ; Gene Frequency ; Female ; Male ; Middle Aged ; Adult ; Asian People

OBJECTIVE: To explore the correlation between single nucleotide polymorphisms (SNPs) of ARID5B gene and the risk of acute lymphoblastic leukemia (ALL) and minimal residual disease (MRD) in children of Hui and Han nationality in Ningxia. METHODS: In this case-control study, 54 ALL children and control group with matched age, sex and nationality were detected for the polymorphism of ARID5B gene using fluorescence resonance energy transfer technique, and the susceptibility of different ALL genotypes and their correlation with MRD were analyzed. RESULTS: There were no significant differences in genotype and allele frequency of rs10994982, rs7089424, rs10740055, rs7073837, rs4245595 and rs7090445 between the two groups (P >0.05). At the locus of rs10821936, the frequencies of T/T genotype and T allele in ALL group were significantly higher than those in the control group (both P < 0.05). The C/C genotype of ARID5B gene SNP rs10821936 was a risk factor for early MRD positive in ALL children ( P < 0.05). CONCLUSION ARID5B gene SNP rs10821936 is related to the development of childhood ALL and MRD.
Humans ; Precursor Cell Lymphoblastic Leukemia-Lymphoma/genetics* ; Polymorphism, Single Nucleotide ; Case-Control Studies ; Neoplasm, Residual/genetics* ; DNA-Binding Proteins/genetics* ; Transcription Factors/genetics* ; Genotype ; Genetic Predisposition to Disease ; Gene Frequency ; Child ; Male ; Female ; Alleles ; Risk Factors ; Child, Preschool

Objective： To investigate the relationship between the lipid metabolism-related gene sterol carrier protein 2(SCP2) and prostate cancer (PCa) from a multi-omics perspective using single-cell transcriptomes combined with Mendelian randomization. Methods： Single-cell transcriptome data of benign and malignant prostate tissues were obtained from GSE120716, GSE157703 and GSE141445 datasets, respectively. Integration, quality control and annotation were performed on the data to categorize the epithelial cells into high and low SCP2 expression groups, followed by further differential and trajectory analyses. Single nucleotide polymorphism (SNP) data for SCP2 expression quantitative trait loci (eQTL) were subsequently downloaded from Genotype-Tissue Expression (GTEx) and investigated from the PCa Society Cancer-Related Genomic Alteration Panel for the Investigation of Cancer-Related Alterations (PRACTICAL) to obtain PCa outcome data for Mendelian randomization analysis to validate the causal relationship between SCP2 and PCa. Results： High SCP2-expressing epithelial cells had higher energy metabolism and proliferation capacity with low immunotherapy response and metastatic tendency. Trajectory analysis showed that epithelial cells with high SCP2 expression may have a higher degree of malignancy, and SCP2 may be a key marker gene for differentiation of malignant epithelial cells in the prostate. Further Mendelian randomization results showed a significant causal relationship between SCP2 and PCa development (OR=1.045, 95% CI: 1.010 －1.083, P=0.011). Conclusion： By combining single-cell transcriptome and Mendelian randomization, the role of the lipid metabolism-related gene SCP2 in PCa development has been confirmed, and new targets and therapeutic directions for PCa treatment have been provided.
Humans ; Prostatic Neoplasms/genetics* ; Male ; Mendelian Randomization Analysis ; Polymorphism, Single Nucleotide ; Quantitative Trait Loci ; Single-Cell Analysis ; Sequence Analysis, RNA ; Carrier Proteins/genetics* ; Transcriptome ; Lipid Metabolism

BACKGROUND: Many factors are associated with the development and progression of liver fat and fibrosis; however, genetics and the gut microbiota are representative factors. Moreover, recent studies have indicated a link between host genes and the gut microbiota. This study investigated the effect of patatin-like phospholipase domain-containing 3 (PNPLA3) rs738409 (C > G), which has been reported to be most involved in the onset and progression of fatty liver, on liver fat and fibrosis in a cohort study related to gut microbiota in a non-fatty liver population. METHODS: This cohort study included 214 participants from the health check-up project in 2018 and 2022 who had non-fatty liver with controlled attenuation parameter (CAP) values <248 dB/m by FibroScan and were non-drinkers. Changes in CAP values and liver stiffness measurement (LSM), liver-related items, and gut microbiota from 2018 to 2022 were investigated separately for PNPLA3 rs738409 CC, CG, and GG genotypes. RESULTS: Baseline values showed no difference among the PNPLA3 rs738409 genotypes for any of the measurement items. From 2018 to 2022, the PNPLA3 rs738409 CG and GG genotype groups showed a significant increase in CAP and body mass index; no significant change was observed in the CC genotype group. LSM increased in all genotypes, but the rate of increase was highest in the GG genotype, followed by the CG and CC genotypes. Fasting blood glucose levels increased in all genotypes; however, HOMA-IR (Homeostasis Model Assessment of Insulin Resistance) increased significantly only in the GG genotype. HDL (high-density lipoprotein) and LDL (low-density lipoprotein) cholesterol levels significantly increased in all genotypes, whereas triglycerides did not show any significant changes in any genotype. As for the gut microbiota, the relative abundance of Feacalibacterium in the PNPLA3 rs738409 GG genotype decreased by 2% over 4 years, more than 2-fold compared to CC and GG genotypes. Blautia increased significantly in the CC group. CONCLUSION The results suggest that PNPLA3 G-allele carriers of non-fatty liver develop liver fat and fibrosis due to not only obesity and insulin resistance but also the deterioration of gut microbiota, which may require a relatively long course of time, even years.
Humans ; Gastrointestinal Microbiome ; Male ; Female ; Membrane Proteins/metabolism* ; Lipase/genetics* ; Middle Aged ; Liver Cirrhosis/epidemiology* ; Cohort Studies ; Genotype ; Adult ; Non-alcoholic Fatty Liver Disease/microbiology* ; Polymorphism, Single Nucleotide ; Acyltransferases ; Phospholipases A2, Calcium-Independent

BACKGROUND: Epstein-Barr (EB) virus infection stimulates the production of vascular endothelial growth factor (VEGF), which contributes to the progression of angiogenesis. Angiogenesis plays an important role in the development of atherosclerosis. Since serum anti-early antigen EB virus IgG (EBV EA-IgG) titer is a sign of active EB virus infection, EBV EA-IgG titer could be associated with atherosclerosis. The number of minor (T) alleles in VEGF polymorphism rs3025039 has been reported to be inversely associated with serum VEGF concentration, suggesting that rs3025039 might have a strong influence on the association between EBV EA-IgG titer and atherosclerosis. By focusing on the role of VEGF in the development of atherosclerosis, this study aimed to investigate the association between active EB virus infection and atherosclerosis. METHODS: A cross-sectional study of 2,661 older Japanese individuals aged 60-89 years who participated in annual health check-ups during 2017-2019 was conducted. Logistic regression was used to evaluate the association between EBV EA-IgG titer and atherosclerosis in relation to rs3025039 genotype. The influence of rs3025039 (T) allele carrier status on the association between EBV EA-IgG titer and atherosclerosis was also evaluated by using logistic regression. RESULTS: Among rs3025039 CC-homozygotes, with the lowest EBV EA-IgG titer tertile as the reference, the multivariable odds ratio (95% confidence interval) was 1.11 (0.82, 1.50) for the medium tertile and 1.07 (0.78, 1.47) for the high tertile. Among rs3025039 (T) allele carriers, the corresponding values were 1.44 (0.88, 2.36) and 1.88 (1.15, 3.05), respectively. There was a significant interaction between rs3025039 (T) allele carrier status and the association between EBV EA-IgG titer and atherosclerosis (adjusted p = 0.0497). CONCLUSION EBV EA-IgG titer was significantly positively associated with atherosclerosis only among participants who are genetically less likely to have progressive angiogenesis. An angiogenesis-related genetic factor was revealed as a determinant of the association between EBV EA-IgG titer and atherosclerosis. These findings introduce a novel concept that could explain the association between viral infection and atherosclerosis.
Humans ; Aged ; Male ; Middle Aged ; Female ; Japan/epidemiology* ; Atherosclerosis/virology* ; Aged, 80 and over ; Vascular Endothelial Growth Factor A/genetics* ; Epstein-Barr Virus Infections/virology* ; Polymorphism, Single Nucleotide ; Cross-Sectional Studies ; Herpesvirus 4, Human ; Antigens, Viral/immunology* ; Antibodies, Viral/blood* ; Immunoglobulin G/blood* ; Genotype ; East Asian People

OBJECTIVES: The major histocompatibility complex class I chain-related gene A (MICA), a component of the human leukocyte antigen (HLA) gene complex, is involved in the pathogenesis of various diseases including cancers and autoimmune disorders. Rosacea, a chronic inflammatory skin disease with a complex pathogenesis, potentially influenced by genetic and autoimmune factors. This study aims to investigate the relationship among MICA gene polymorphisms, single nucleotide polymorphisms (SNPs), and susceptibility to rosacea, thereby offering new insights into the disease mechanism. METHODS: Peripheral blood DNA samples were collected from 84 patients with rosacea (rosacea group) and 223 healthy volunteers (control group) who visited the Dermatology Outpatient Department of Xiangya Hospital of Central South University between November 2017 and November 2019. MICA genotyping was performed using polymerase chain reaction-sequencing-based typing (PCR-SBT) and the next-generation sequencing (NGS), and the accuracy of the 2 methods was compared. The frequency distributions of MICA alleles between the 2 groups were analyzed. Amino acid clustering and SNP site analyses were conducted to identify haplotype-linked SNPs and to classify MICA polymorphic variants. Distribution differences of these classifications between groups were also examined. RESULTS: Blood tests in rosacea patients showed mildly elevated, with no significant changes in lymphocyte counts. Both PCR-SBT and NGS accurately identified MICA alleles. The most common alleles in the rosacea group were MICA*010:01, MICA*008:04, and MICA*019:01. The frequencies of MICA*002:01 and MICA*027 were significantly lower in the rosacea group compared to controls (6.55% vs 18.16% and 1.19% vs 5.38%, respectively), while and MICA*010:01 were significantly higher (7.74% vs 3.36% and 31.55% vs 18.61%, respectively; all P<0.05). Five short tandem repeat (STR) alleles were identified. Frequencies of MICA-A4 and MICA-A9 were lower in the rosacea group than in the control group (16.07% vs 23.32% and 7.74% vs 17.26%, respectively), whereas MICA-A6 was higher (10.12% vs 4.03%; all P<0.05). Clustering and SNP analysis identified 6 linked SNP sites, classifying MICA variants into Type I (C₃₆+M₁₂₉+K₁₇₃+G₂₀₆+W₂₁₀+S₂₁₅) and Type II (Y₃₆+V₁₂₉+E₁₇₃+S₂₀₆+R₂₁₀+T₂₁₅). Type I MICA variants were significantly associated with rosacea susceptibility. CONCLUSIONS MICA gene polymorphisms are associated with susceptibility to rosacea, and there are 6 linked SNP sites within the MICA gene. Based on this, MICA polymorphic variants are classified into Type I and Type II, with Type I being more closely associated with disease development of rosacea.
Humans ; Polymorphism, Single Nucleotide ; Histocompatibility Antigens Class I/genetics* ; Rosacea/genetics* ; Genetic Predisposition to Disease/genetics* ; Female ; Male ; Adult ; Middle Aged ; Genotype ; Alleles ; Gene Frequency ; Haplotypes ; Case-Control Studies ; Aged ; High-Throughput Nucleotide Sequencing

OBJECTIVES: The gut microbiota plays a crucial role in the pathophysiology of various types of diabetes. However, the causal relationship between them has yet to be systematically elucidated. This study aims to explore the potential causal associations between gut microbiota and diabetes using a two-sample Mendelian randomization (MR) analysis, based on multiple taxonomic levels. METHODS: Eligible instrumental variables were extracted from the selected genome-wide association study (GWAS) data on gut microbiota. These were combined with GWAS datasets on type 1 diabetes (T1D), type 2 diabetes (T2D), and gestational diabetes mellitus (GDM) to conduct forward MR analysis, sensitivity analysis, reverse MR analysis, and validation of significant estimates. Microbial taxa with causal effects on T1D, T2D, and GDM were identified based on a comprehensive assessment of all analytical stages. RESULTS: A total of 2 179, 2 176, and 2 166 single nucleotide polymorphisms (SNP) were included in the MR analyses for gut microbiota with T1D, T2D, and GDM, respectively. MR results indicated causal associations between: Six microbial taxa (Eggerthella, Lachnospira, Bacillales, Desulfovibrionales, Parasutterella, and Turicibacter) and T1D; 9 microbial taxa (Verrucomicrobia, Deltaproteobacteria, Actinomycetales, Desulfovibrionale, Actinomycetaceae, Desulfovibrionaceae, Actinomyces, Alcaligenaceae, and Lachnospiraceae NC2004 group) and T2D; 10 microbial taxa (Betaproteobacteria, Coprobacter, Ruminococcus2, Tenericutes, Clostridia, Methanobacteria, Mollicutes, Methanobacteriales, Methanobacteriaceae, and Methanobrevibacter) and GDM. CONCLUSIONS This study identified specific gut microbial taxa that may significantly increase or decrease the risk of developing diabetes. Some findings were fully replicated in independent validation datasets. However, the underlying biological mechanisms of these causal relationships warrant further investigation through mechanistic studies and population-based research.
Gastrointestinal Microbiome/genetics* ; Humans ; Mendelian Randomization Analysis ; Genome-Wide Association Study ; Diabetes Mellitus, Type 2/genetics* ; Diabetes Mellitus, Type 1/genetics* ; Female ; Polymorphism, Single Nucleotide ; Diabetes, Gestational/genetics* ; Pregnancy

OBJECTIVES: Keloids are fibrotic skin disorders characterized by excessive collagen deposition and a high recurrence rate, closely associated with inflammatory mediators. However, existing epidemiological studies are limited by confounding factors and reverse causality, making it difficult to establish causation. This study aims to investigate the causal relationship between circulating cytokines and keloids using Mendelian randomization analysis. METHODS: Significant single nucleotide polymorphisms (SNPs) associated with circulating cytokines (exposures) and keloids (outcomes) were extracted from genome-wide association study (GWAS) summary datasets. Eligible SNPs were selected as instrumental variables (IVs). Exposure data were derived from a cytokine GWAS including 8 293 Finnish participants, and outcome data from a keloid GWAS based on the UK Biobank. The inverse-variance weighted (IVW) method served as the primary analytical approach to estimate causal effects, supplemented by weighted median (WME), MR-Egger regression, and other sensitivity analyses. Horizontal pleiotropy was assessed using MR-Egger regression and the MR pleiotropy residual sum and outlier (MR-PRESSO) test, while Cochran's Q test evaluated heterogeneity. Leave-one-out analysis was used to verify robustness and consistency. A reverse MR analysis was also conducted, with keloid as the exposure and cytokines as outcomes, to rule out reverse causation. RESULTS: IVW analysis identified significant positive causal associations between two cytokines and keloids-macrophage migration inhibitory factor (MIF) [odds ratio (OR)=2.081, 95% confidence interval (CI) 1.219 to 3.552, P=0.007] and monocyte chemoattractant protein-1 (MCP-1) (OR=1.673, 95% CI 1.036 to 2.701, P=0.035). Conversely, stem cell factor (SCF) showed a negative causal relationship with keloids (OR=0.518, 95% CI 0.269 to 0.998, P=0.049). Results from the MR-Egger and weighted median analyses were consistent with IVW findings. No evidence of horizontal pleiotropy was observed (P>0.05). Except for interleukin-6 (P=0.014), no heterogeneity was detected in other cytokines. Leave-one-out analysis further confirmed the robustness of the causal associations. In reverse MR analysis, keloids were causally related only to β-nerve growth factor (beta-NGF) (OR=1.048, 95% CI 1.002 to 1.095, P=0.039), with no heterogeneity or pleiotropy detected in most cytokines (P>0.05). CONCLUSIONS MIF and MCP-1 exhibit positive causal associations with keloid formation, while SCF shows a negative causal relationship. These findings provide new evidence for the causal involvement of inflammatory cytokines in keloid pathogenesis and offer potential molecular targets for developing novel keloid therapies.
Humans ; Keloid/blood* ; Mendelian Randomization Analysis ; Cytokines/genetics* ; Polymorphism, Single Nucleotide ; Genome-Wide Association Study ; Chemokine CCL2/genetics* ; Interleukin-6/genetics* ; Macrophage Migration-Inhibitory Factors/genetics* ; Male ; Stem Cell Factor/blood* ; Female ; Intramolecular Oxidoreductases