1.Association of microRNA gene polymorphisms with risk, clinicopathological characteristics and therapeutical efficacy among Chinese patients with Crohn's disease.
Yanlun ZHANG ; Xiaoxiao SHAO ; Daopo LIN ; Yuan XU ; Guolong MA ; Yi JIANG
Chinese Journal of Medical Genetics 2026;43(2):111-122
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
2.Causal relationship between autoimmune diseases and aplastic anemia: A Mendelian randomization study.
Wenjie LI ; Yaonan HONG ; Rui HUANG ; Yuchen LI ; Ying ZHANG ; Yun ZHANG ; Dijiong WU
Journal of Southern Medical University 2025;45(4):871-879
OBJECTIVES:
To investigate the causal associations between autoimmune diseases and aplastic anemia (AA) using Mendelian randomization analysis.
METHODS:
Publicly available genome-wide association study (GWAS) data were utilized to obtain single nucleotide polymorphisms (SNPs) associated with autoimmune diseases and AA for analysis. The inverse variance weighted (IVW) method was employed as the primary analytical approach, with MR Egger, Weighted Mode, Weighted Median, and Simple Mode methods serving as complementary analyses. Heterogeneity and pleiotropy analyses were conducted using designated functions, and the robustness of Mendelian randomization results was assessed using leave-one-out analysis.
RESULTS:
The two-sample Mendelian randomization analysis using the IVW method revealed significant positive causal associations of rheumatoid arthritis (OR=1.094, 95% CI: 1.023-1.170, P=0.009, adjusted P=0.042), systemic lupus erythematosus (OR=1.111, 95% CI: 1.021-1.208, P=0.015, adjusted P=0.036), Hashimoto thyroiditis (OR=1.206, 95% CI: 1.049-1.387, P=0.009, adjusted P=0.029), and Sicca syndrome (OR=1.173, 95% CI: 1.054-1.306, P=0.004, adjusted P=0.035) with AA, which was supported by the results from the Weighted Median method. Sensitivity analyses indicated no evidence of pleiotropy or heterogeneity, and leave-one-out analysis confirmed the robustness of the causal relationships. No direct evidence was found linking Graves' disease, ulcerative colitis, Crohn's disease, autoimmune hepatitis, primary biliary cholangitis, or primary sclerosing cholangitis with AA (P>0.05, adjusted P>0.05), indicating a lack of causal association. Reverse Mendelian randomization results and multiple corrections indicated that AA was not an influencing factor for autoimmune diseases (adjusted P>0.05).
CONCLUSIONS
Our findings support at the genetic level that rheumatoid arthritis, systemic lupus erythematosus, Hashimoto thyroiditis, and Sicca syndrome are risk factors for AA, and confirm a causal association of the these 4 autoimmune diseases with an increased risk of AA.
Humans
;
Mendelian Randomization Analysis
;
Anemia, Aplastic/genetics*
;
Autoimmune Diseases/complications*
;
Polymorphism, Single Nucleotide
;
Genome-Wide Association Study
;
Arthritis, Rheumatoid/genetics*
;
Lupus Erythematosus, Systemic/genetics*
;
Genetic Predisposition to Disease
3.HOTAIR rs920778 single nucleotide polymorphism is associated with breast cancer susceptibility and HER2-targeted therapy resistance in Chinese population.
Mingliang ZHANG ; Feifan SUN ; Zhuoqi HAN ; Yue GAO ; Yi LUO
Journal of Southern Medical University 2025;45(10):2270-2276
OBJECTIVES:
To investigate the association of HOTAIR gene rs920778 single nucleotide polymorphism (SNP) with breast cancer susceptibility and response to HER2-targeted therapy in a Chinese population.
METHODS:
TaqMan probe-based real-time quantitative PCR was used for genotyping of the rs920778 locus (chr12:54,376,218) in peripheral blood genomic DNA from 287 breast cancer patients and 260 healthy individuals from northern Anhui Province. The genotype (GG, GT and TT) and allele (G/T) distribution frequencies were compared between the two groups to evaluate their association with breast cancer risk. Multivariate logistic regression analysis was conducted to assess the relationship between SNP at this locus and aggressive clinicopathological features (including tumor size, lymph node metastasis, ER/PR/HER2 status, and molecular subtypes) of breast cancer. For the HER2-positive subgroup, the association between rs920778 genotype and responses to dual-targeted therapy (trastuzumab [6 mg/kg q3w]+pertuzumab [420 mg q3w] + docetaxel [75 mg/m²]) was analyzed. The primary endpoints included pathological complete response rate (pCR), objective response rate (ORR), and progression-free survival (PFS).
RESULTS:
The TT genotype of rs920778 was associated with a significantly increased breast cancer susceptibility (OR=1.54, 95% CI: 1.09-2.19; P=0.017), an advanced tumor stage (P<0.001), lymph node metastasis (P<0.001), and the triple-negative subtype (P<0.001). In HER2-positive patients, TT genotype carriers had a markedly reduced objective response rate to dual HER2-targeted therapy (33.3% vs 89.3%, P=0.001) and a lower pathological complete response rate after neoadjuvant therapy (P=0.018).
CONCLUSIONS
The TT genotype of HOTAIR rs920778 serves as an independent risk factor for breast cancer susceptibility and aggressive progression in Chinese population and may predict the resistance to HER2-targeted therapies, suggesting its potential as a prognostic biomarker for precision oncology.
Adult
;
Aged
;
Female
;
Humans
;
Middle Aged
;
Breast Neoplasms/drug therapy*
;
Case-Control Studies
;
China
;
Drug Resistance, Neoplasm/genetics*
;
Genetic Predisposition to Disease
;
Genotype
;
Polymorphism, Single Nucleotide
;
Receptor, ErbB-2
;
RNA, Long Noncoding/genetics*
;
East Asian People/genetics*
4.Genetic Variation A118G in the OPRM1 Gene Underlies the Dimorphic Response to Epidural Opioid-Induced Itch.
Xiaomeng ZHOU ; Ai-Lun LI ; Wan-Jie DU ; Pengyu GAO ; Bin LAI ; Fang FANG ; Qingjian HAN ; Jing CANG
Neuroscience Bulletin 2025;41(12):2272-2284
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
5.Developing a polygenic risk score for pelvic organ prolapse: a combined risk assessment approach in Chinese women.
Xi CHENG ; Lei LI ; Xijuan LIN ; Na CHEN ; Xudong LIU ; Yaqian LI ; Zhaoai LI ; Jian GONG ; Qing LIU ; Yuling WANG ; Juntao WANG ; Zhijun XIA ; Yongxian LU ; Hangmei JIN ; Xiaowei ZHANG ; Luwen WANG ; Juan CHEN ; Guorong FAN ; Shan DENG ; Sen ZHAO ; Lan ZHU
Frontiers of Medicine 2025;19(4):665-674
Pelvic organ prolapse (POP), whose etiology is influenced by genetic and clinical risk factors, considerably impacts women's quality of life. However, the genetic underpinnings in non-European populations and comprehensive risk models integrating genetic and clinical factors remain underexplored. This study constructed the first polygenic risk score (PRS) for POP in the Chinese population by utilizing 20 disease-associated variants from the largest existing genome-wide association study. We analyzed a discovery cohort of 576 cases and 623 controls and a validation cohort of 264 cases and 200 controls. Results showed that the case group exhibited a significantly higher PRS than the control group. Moreover, the odds ratio of the top 10% risk group was 2.6 times higher than that of the bottom 10%. A high PRS was significantly correlated with POP occurrence in women older than 50 years old and in those with one or no childbirths. As far as we know, the integrated prediction model, which combined PRS and clinical risk factors, demonstrated better predictive accuracy than other existing PRS models. This combined risk assessment model serves as a robust tool for POP risk prediction and stratification, thereby offering insights into individualized preventive measures and treatment strategies in future clinical practice.
Humans
;
Female
;
Pelvic Organ Prolapse/epidemiology*
;
Middle Aged
;
Risk Assessment/methods*
;
China/epidemiology*
;
Multifactorial Inheritance
;
Aged
;
Risk Factors
;
Genome-Wide Association Study
;
Genetic Predisposition to Disease
;
Case-Control Studies
;
Adult
;
Polymorphism, Single Nucleotide
;
Genetic Risk Score
;
East Asian People
6.Predicting Diabetic Retinopathy Using a Machine Learning Approach Informed by Whole-Exome Sequencing Studies.
Chong Yang SHE ; Wen Ying FAN ; Yun Yun LI ; Yong TAO ; Zu Fei LI
Biomedical and Environmental Sciences 2025;38(1):67-78
OBJECTIVE:
To establish and validate a novel diabetic retinopathy (DR) risk-prediction model using a whole-exome sequencing (WES)-based machine learning (ML) method.
METHODS:
WES was performed to identify potential single nucleotide polymorphism (SNP) or mutation sites in a DR pedigree comprising 10 members. A prediction model was established and validated in a cohort of 420 type 2 diabetic patients based on both genetic and demographic features. The contribution of each feature was assessed using Shapley Additive explanation analysis. The efficacies of the models with and without SNP were compared.
RESULTS:
WES revealed that seven SNPs/mutations ( rs116911833 in TRIM7, 1997T>C in LRBA, 1643T>C in PRMT10, rs117858678 in C9orf152, rs201922794 in CLDN25, rs146694895 in SH3GLB2, and rs201407189 in FANCC) were associated with DR. Notably, the model including rs146694895 and rs201407189 achieved better performance in predicting DR (accuracy: 80.2%; sensitivity: 83.3%; specificity: 76.7%; area under the receiver operating characteristic curve [AUC]: 80.0%) than the model without these SNPs (accuracy: 79.4%; sensitivity: 80.3%; specificity: 78.3%; AUC: 79.3%).
CONCLUSION
Novel SNP sites associated with DR were identified in the DR pedigree. Inclusion of rs146694895 and rs201407189 significantly enhanced the performance of the ML-based DR prediction model.
Diabetic Retinopathy/diagnosis*
;
Humans
;
Machine Learning
;
Male
;
Female
;
Polymorphism, Single Nucleotide
;
Middle Aged
;
Exome Sequencing
;
Aged
;
Adult
;
Pedigree
;
Diabetes Mellitus, Type 2/complications*
;
Genetic Predisposition to Disease
;
Mutation
7.Causal Associations between Particulate Matter 2.5 (PM 2.5), PM 2.5 Absorbance, and Inflammatory Bowel Disease Risk: Evidence from a Two-Sample Mendelian Randomization Study.
Xu ZHANG ; Zhi Meng WU ; Lu ZHANG ; Bing Long XIN ; Xiang Rui WANG ; Xin Lan LU ; Gui Fang LU ; Mu Dan REN ; Shui Xiang HE ; Ya Rui LI
Biomedical and Environmental Sciences 2025;38(2):167-177
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
8.Threshold-Effect Associations of Serum 25-hydroxyvitamin D on Bone Turnover Markers and GC rs2282679 Variants in Chinese Women of Childbearing Age.
Xiao Yun SHAN ; Yu Ting LI ; Xia Yu ZHAO ; Yi Chun HU ; Si Ran LI ; Hui di ZHANG ; Yang CAO ; Rui WANG ; Li Chen YANG
Biomedical and Environmental Sciences 2025;38(4):433-446
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
9.Exploration of New Susceptible Genes associated with Non-Alcoholic Fatty Liver Disease among Children with Obesity Using Whole Exome Sequencing.
Xiong Feng PAN ; Cai Lian WEI ; Jia You LUO ; Jun Xia YAN ; Xiang XIAO ; Jie WANG ; Yan ZHONG ; Mi Yang LUO
Biomedical and Environmental Sciences 2025;38(6):727-739
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
10.Association analysis between forkhead box E1 gene and non-syndromic cleft lip with or without cleft palate in Han Chinese population.
Sixuan JIA ; Sidi ZHANG ; Yue YOU ; Jialin SUN ; Shijun DUAN ; Bing SHI ; Zhonglin JIA
West China Journal of Stomatology 2025;43(1):28-36
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*

Result Analysis
Print
Save
E-mail