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

This study investigates the genetic diversity and evolutionary relationships of different Artemisia argyi germplasm resources to provide a basis for germplasm identification, variety selection, and resource protection. A total of 192 germplasm resources of A. argyi were studied, and EST-based simple sequence repeat(EST-SSR) primers were designed based on transcriptomic data of A. argyi. Polymerase chain reaction(PCR) amplification was performed on these resources, followed by fluorescence capillary electrophoresis to detect genetic diversity and construct DNA fingerprints. From 197 pairs of primers designed, 28 pairs with polymorphic and clear bands were selected. A total of 278 alleles were detected, with an average of 9.900 0 alleles per primer pair and an average effective number of alleles of 1.407 2. The Shannon's diversity index(I) for the A. argyi germplasm resources ranged from 0.148 1 to 0.418 0, with an average of 0.255 7. The polymorphism information content(PIC) ranged from 0.454 5 to 0.878 0, with an average of 0.766 9, showing high polymorphism. Cluster analysis divided the A. argyi germplasm resources into three major groups: Group Ⅰ contained 136 germplasm samples, Group Ⅱ contained 45, and Group Ⅲ contained 11. Principal component analysis also divided the resources into three groups, which was generally consistent with the clustering results. Mantel test results showed that the genetic variation in A. argyi populations was to some extent influenced by geographic distance, but the effect was minimal. Structure analysis showed that 190 germplasm materials had Q≥ 0.6, indicating that these germplasm materials had a relatively homogeneous genetic origin. Furthermore, 8 core primer pairs were selected from the 28 designed primers, which could distinguish various germplasm types. Using these 8 core primers, DNA fingerprints for the 192 A. argyi germplasm resources were successfully constructed. EST-SSR molecular markers can be used to study the genetic diversity and phylogenetic relationships of A. argyi, providing theoretical support for the identification and molecular-assisted breeding of A. argyi germplasm resources.
Artemisia/classification* ; Microsatellite Repeats ; Genetic Variation ; Expressed Sequence Tags ; DNA Fingerprinting ; Phylogeny ; Polymorphism, Genetic ; DNA, Plant/genetics* ; Genetic Markers

OBJECTIVES: To clarify the causal relationship between the level of cytoplasmic unactivated mineralocorticoid receptor (MR) and the development of tubulointerstitial nephritis (TIN), and to evaluate the impact of MR on dyslipidemia, particularly secondary hyperlipemia, in patients with diabetic kidney disease. METHODS: We conducted a two-sample Mendelian randomization study using genome-wide association study (GWAS) summary data. Genetic variants associated with MR levels were selected as exposures, with TIN and lipid profiles [including low-density lipoprotein cholesterol (LDL-C), triglyceride, and high-density lipoprotein cholesterol] as outcomes. A two-step Mendelian randomization approach was used to assess TIN as a mediator, employing inverse variance weighted regression as the primary analysis, supplemented by Mendelian randomization-Egger, weighted median, and sensitivity analyses. RESULTS: Cytoplasmic unactivated MR level exhibited a significant causal association with a decreased risk of TIN (OR = 0.8598, 95% CI [0.7775-0.9508], P < 0.001). Although no significant causal relationship was identified between MR level and secondary hyperlipemia, a potential association of cytoplasmic unactivated MR level with lower LDL-C levels was observed (OR = 0.9901, 95% CI [0.9821-0.9983], P = 0.018). Additionally, TIN exhibited causal links with secondary hyperlipemia (OR = 1.0016, 95% CI [1.0002-1.0029], P = 0.020) and elevated LDL-C (OR = 1.0111, 95% CI [1.0024-1.0199], P = 0.012), particularly LDL-C in European males (OR = 1.0230, 95% CI [1.0103-1.0358], P < 0.001). Inverse Mendelian randomization analysis revealed causal relationships between TIN and genetically predicted triglyceride (OR = 0.7027, 95% CI [0.6189-0.7978], P < 0.001), high-density lipoprotein cholesterol (OR = 1.1247, 95% CI [1.0019-1.2626], P = 0.046), and LDL-C (OR = 0.8423, 95% CI [0.7220-0.9827], P = 0.029). Notably, TIN mediated 16.7% of the causal association between MR and LDL-C levels. CONCLUSIONS MR plays a critical role in the development of TIN and lipid metabolism, highlighting the potential of MR-antagonists in reducing renal damage and lipid metabolism-associated complications.
Humans ; Mendelian Randomization Analysis ; Nephritis, Interstitial/metabolism* ; Receptors, Mineralocorticoid/genetics* ; Lipid Metabolism/genetics* ; Genome-Wide Association Study ; Male ; Female ; Polymorphism, Single Nucleotide ; Dyslipidemias/metabolism*

Neuraxial opioids, widely used in obstetric and perioperative pain management, often lead to unwanted itch, reducing patient satisfaction. While the μ-opioid receptor has been implicated in opioid-induced itch, the genetic basis for variable itch incidence remains unknown. This study examined 3616 patients receiving epidural opioids, revealing an itch occurrence of 26.55%, with variations among opioid types and gender. Analysis of the OPRM1 gene identified six single-nucleotide polymorphisms, notably rs1799971 (A118G), that correlated with opioid-induced itch. Mouse models with an equivalent A112G mutation showed reduced neuraxial opioid-induced itch and light touch-evoked itch, mirroring human findings. The 118G allele demonstrated an anti-itch effect without impacting analgesia, addiction, or tolerance, offering insights for risk stratification and potential anti-itch pretreatment strategies.
Receptors, Opioid, mu/genetics* ; Pruritus/chemically induced* ; Humans ; Analgesics, Opioid/administration & dosage* ; Female ; Male ; Animals ; Polymorphism, Single Nucleotide/genetics* ; Adult ; Mice ; Middle Aged

OBJECTIVE: Several epidemiological observational studies have related particulate matter (PM) exposure to Inflammatory bowel disease (IBD), but many confounding factors make it difficult to draw causal links from observational studies. The objective of this study was to explore the causal association between PM _2.5 exposure, its absorbance, and IBD. METHODS: We assessed the association of PM _2.5 and PM _2.5 absorbance with the two primary forms of IBD (Crohn's disease [CD] and ulcerative colitis [UC]) using Mendelian randomization (MR) to explore the causal relationship. We conducted two-sample MR analyses with aggregated data from the UK Biobank genome-wide association study. Single-nucleotide polymorphisms linked with PM _2.5 concentrations or their absorbance were used as instrumental variables (IVs). We used inverse variance weighting (IVW) as the primary analytical approach and four other standard methods as supplementary analyses for quality control. RESULTS: The results of MR demonstrated that PM _2.5 had an adverse influence on UC risk (odds ratio [ OR] = 1.010; 95% confidence interval [ CI] = 1.001-1.019, P = 0.020). Meanwhile, the results of IVW showed that PM _2.5 absorbance was also causally associated with UC ( OR = 1.012; 95% CI = 1.004-1.019, P = 0.002). We observed no causal relationship between PM _2.5, PM _2.5 absorbance, and CD. The results of sensitivity analysis indicated the absence of heterogeneity or pleiotropy, ensuring the reliability of MR results. CONCLUSION Based on two-sample MR analyses, there are potential positive causal relationships between PM _2.5, PM _2.5 absorbance, and UC.
Humans ; Mendelian Randomization Analysis ; Particulate Matter/analysis* ; Polymorphism, Single Nucleotide ; Inflammatory Bowel Diseases/genetics* ; Air Pollutants/analysis* ; Crohn Disease/genetics* ; Colitis, Ulcerative/genetics* ; Genome-Wide Association Study ; Risk Factors ; Environmental Exposure

OBJECTIVE: This study aimed to investigate possible serum 25-hydroxyvitamin D [25(OH)D] cutoffs for the associations between 25(OH)D and Bone turnover markers (BTMs), and how GC gene variation influences such cutoffs in Chinese women of childbearing age. METHODS: In total, 1,505 non-pregnant or non-lactating women (18-45 years) were recruited from the 2015 Chinese Adult Chronic Disease and Nutrition Surveillance. Serum 25(OH)D, osteocalcin (OC), procollagen type 1 N-terminal propeptide (P1NP), β-CrossLaps of type 1 collagen containing cross-linked C-telopeptide (β-CTX), and single nucleotide polymorphisms were determined. Locally weighted regression and smoothing scatterplot and segmented regression were performed to estimate the 25(OH)D thresholds. RESULTS: The median serum 25(OH)D was 16.63 (11.96-22.55) ng/mL and the prevalence of low serum 25(OH)D (< 12 ng/mL) was 25.2%. Women with the lowest 25(OH)D had the highest β-CTX. After adjustment for the confounders, 25(OH)D cutoffs for OC [14.04 (12.84-15.23) ng/mL], β-CTX [13.94 (12.49-15.39) ng/mL], and P1NP [13.87 (12.37-15.37) ng/mL] in the whole population, cutoffs for OC [12.30 (10.68-13.91) ng/mL], β-CTX [12.23 (10.22-14.23) ng/mL], and P1NP [11.85 (10.40-13.31) ng/mL] in women with the GC rs2282679 G allele, and cutoffs for OC [12.75 (11.81-13.68) ng/mL], β-CTX [13.05 (11.78-14.32) ng/mL], and P1NP [12.81 (11.57-14.06) ng/mL] in women with the GC rs2282679 T allele, were observed. Below these cutoffs, BTMs were negatively associated with 25(OH)D, while above these cutoffs, BTMs plateaued. CONCLUSION In Chinese women of childbearing age, there were thresholds effect of serum 25(OH)D concentrations on BTMs. The results indicated that serum 25(OH)D concentrations < 13.87 ng/mL in this population had adverse influences on maintaining bone remodeling. BTMs were suppressed at a relatively lower serum 25(OH)D in women with the GC rs2282679 G allele compared with those with the T allele.
Humans ; Female ; Vitamin D/blood* ; Adult ; Middle Aged ; Polymorphism, Single Nucleotide ; Adolescent ; Young Adult ; China ; Biomarkers/blood* ; Bone Remodeling/genetics* ; Vitamin D-Binding Protein/genetics* ; Procollagen/blood* ; Osteocalcin/blood* ; Peptide Fragments/blood* ; East Asian People

OBJECTIVE: This study aimed to evaluate the association between susceptibility genes and non-alcoholic fatty liver disease (NAFLD) in children with obesity. METHODS: We conducted a two-step case-control study. Ninety-three participants were subjected to whole-exome sequencing (exploratory set). Differential genes identified in the small sample were validated in 1,022 participants using multiplex polymerase chain reaction and high-throughput sequencing (validation set). RESULTS: In the exploratory set, 14 genes from the NAFLD-associated pathways were identified. In the validation set, after adjusting for sex, age, and body mass index, ECI2 rs2326408 (dominant model: OR = 1.33, 95% CI: 1.02-1.72; additive model: OR = 1.22, 95% CI: 1.01-1.47), C6orf201 rs659305 (dominant model: OR = 1.30, 95% CI: 1.01-1.69; additive model: OR = 1.21, 95% CI: 1.00-1.45), CALML5 rs10904516 (pre-ad dominant model: OR = 1.36, 95% CI: 1.01-1.83; adjusted dominant model: OR = 1.40, 95% CI: 1.03-1.91; and pre-ad additive model: OR = 1.26, 95% CI: 1.04-1.66) polymorphisms were significantly associated with NAFLD in children with obesity ( P < 0.05). Interaction analysis revealed that the gene-gene interaction model of CALML5 rs10904516, COX11 rs17209882, and SCD5 rs3733228 was optional ( P < 0.05), demonstrating a negative interaction between the three genes. CONCLUSION In the Chinese population, the CALML5 rs10904516, C6orf201 rs659305, and ECI2 rs2326408 variants could be genetic markers for NAFLD susceptibility.
Humans ; Non-alcoholic Fatty Liver Disease/genetics* ; Child ; Male ; Female ; Genetic Predisposition to Disease ; Case-Control Studies ; Exome Sequencing ; Adolescent ; Polymorphism, Single Nucleotide ; Obesity/complications* ; Pediatric Obesity/complications* ; China

OBJECTIVES: This study aims to explore the association between single nucleotide polymorphisms (SNPs) loci near the haplotype region hg19 chr9:100560865-100660865 of the forkhead box E1 (FOXE1) gene and the occurrence of non-syndromic cleft lip with or without cleft palate (NSCL/P) in western Han Chinese population. METHODS: In the first stage, our study recruited 159 NSCL/P patients and performed targeted region sequencing to screen SNPs loci near the haplotype region of the FOXE1 gene associated with NSCL/P. In the second stage, we selected 21 common SNPs and re-enrolled 1 000 non-syndromic cleft lip only (NSCLO) patients, 1 000 non-syndromic cleft palate only (NSCPO) patients, and 1 000 normal controls to verify the association. PLINK software was used to perform Hardy-Weinberg equilibrium (HWE) test. Association analysis for common variants, gene burden analysis for rare mutations, and function prediction of SNPs with non-synonymous mutations were performed using Mutation Taster and other software programs. RESULTS: In the first stage, 126 variants, including 76 single nucleotide variants and 50 insertion-deletions were identified. All the included SNPs confirmed to HWE, and the results of gene burden analysis and prediction of functional harmfulness for rare variants were not statistically significant. Association analysis showed that rs13292899 of the FOXE1 gene was significantly associated with NSCL/P (P=1.85E-27) and was also correlated with NSCLO (P=6.41E-23) and non-syndromic cleft lip with cleft palate (NSCLP) (P=2.36E-15) subtypes. In the validation phase, rs79268293 (P=0.013, P=0.022), rs10983951 (P=0.009 2, P=0.007 6), rs117227387 (P=0.009 2, P=0.007 6), rs3758250 (P=0.009 2, P=0.007 6), and rs116899397 (P=0.009 2, P=0.007 6) were significantly associated with NSCLO and NSCPO; rs13292899 (P=0.008 5), rs74606599 (P=0.008 3), rs143226042 (P=0.008 3), and rs117236550 (P=0.01) were associated with the occurrence of NSCLO; and rs12343182 (P=0.008 7), rs10119760 (P=0.012), rs10113907 (P=0.012), and rs13299924 (P=0.012) were associated with the occurrence of NSCPO. CONCLUSIONS This study found a new susceptible SNP rs13292899 of the FOXE1 gene that is closely associated with NSCL/P and NSCLO subtype and 13 other SNPs associated with NSCLO or NSCPO.
Female ; Humans ; Male ; China ; Cleft Lip/genetics* ; Cleft Palate/genetics* ; Forkhead Transcription Factors/genetics* ; Haplotypes ; Polymorphism, Single Nucleotide ; East Asian People/genetics*

OBJECTIVES: This study aimed to investigate the association between temporomandibular disorder (TMD) and insomnia using a two-sample Mendelian randomization (MR) approach. METHODS: Bidirectional MR analyses of two samples, TMD (n=377 277) and insomnia (n=375 359), were performed using genome-wide association study statistics published in the FinnGen database. Instrumental variables were first screened, and then inverse variance weighting (IVW) and MR-Egger were used as the main-effect assessment methods. Weighted median, weighted mode, and simple mode served as supplementary methods. We used IVW and MR-Egger to test for heterogeneity, as well as MR-Egger intercepts to assess the single nucleotide polymorphism (SNP) potential level of multiplicity effects. Sensitivity analyses were conducted based on leave-one-out to identify potentially influential SNPs. All analyses were conducted by using the two-sample MR R package and were considered statistically significant when P<0.05. RESULTS: MR analysis showed the presence of TMD on insomnia (OR=1.089, 95%CI: 1.017-1.166, P=0.014). Meanwhile, no effect of insomnia on TMD (OR=0.996, 95%CI: 0.964-1.029, P=0.816) was found. The sensitivity-analysis showed that no heterogeneity existed (P>0.05), and the presence of horizontal pleiotropy was not detected (P>0.05). Leave-one-out sensitivity analysis showed no single SNP, which may affect the causal relation. All findings indicated that the causal relationship between TMD and insomnia was not significantly affected by any individual SNP and that IV did not bias the results. CONCLUSIONS Results of MR analyses showed that TMD is a risk factor for insomnia, whereas insomnia is not a risk factor for TMD.
Humans ; Sleep Initiation and Maintenance Disorders/genetics* ; Mendelian Randomization Analysis ; Polymorphism, Single Nucleotide ; Temporomandibular Joint Disorders/complications* ; Genome-Wide Association Study

We studied the genetic diversity and established the DNA molecular identify for Prunus mume with both ornamental and edible values, aiming to collect, identify, evaluate, and breed new varities of this plant and promote the upgrading of the P. mume industry chain in northern China. We employed 13 pairs of primers with good polymorphism, clear bands, and good repeatability to analyze the genetic diversity and establish the molecular identify of 68 germplasm accessions of P. mume with both ornamental and edible values from Xingtai, Hebei Province. We then employed the unweighted pair-group method with arithmetic means (UPGMA) to perform the cluster analysis based on genetic distance. After that, we analyzed the genetic structure of the 68 germplasm accessions based on a Bayesian model. The 13 pairs of SSR primers amplified a total of 124 alleles from 68 P. mume germplasm accessions, with the mean number of alleles (Na) of 9.538 5, the minor allele frequency (MAF) of 0.369 3, the mean number of effective alleles (Ne) of 4.483 5, and the mean Shannon genetic diversity index (I) of 1.712 4. The mean Nei's gene diversity index (H) of 0.763 7, the mean observed heterozygosity (Ho) of 0.719 5, the mean expected heterozygosity (He) of 0.769 3, the mean polymorphism information content (PIC) of 0.733 6, and the mean genetic similarity (GS) of 0.772 9 suggested that there were significant genetic differences and rich genetic diversity among the studied P. mume germplasm accessions. The cluster analysis revealed that the 68 accessions were classified into three groups, with the mean genetic distance of 0.622 6. The population structure analysis classified the germplasm accessions into two populations. According to the PIC of primers, we selected primers for combination and constructed the combination with the fewest primers required for germplasm differentiation of P. mume with both ornamental and edible values. This study provides a theoretical basis for the innovation and industrial upgrading of P. mume with both ornamental and edible values in gardening and the improvement of breeding efficiency.
Prunus/classification* ; Microsatellite Repeats/genetics* ; Genetic Variation ; China ; Phylogeny ; Polymorphism, Genetic ; DNA, Plant/genetics* ; Alleles