PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “Causal variant”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

BICEP: Bayesian inference for rare genomic variant causality evaluation in pedigrees.

Next-generation sequencing is widely applied to the investigation of pedigree data for gene discovery. However, identifying plausible disease-causing variants within a robust statistical framework is challenging. Here, we introduce BICEP: a Bayesian inference tool for rare variant causality evaluation in pedigree-based cohorts. BICEP calculates the posterior odds that a genomic variant is causal for a phenotype based on the variant cosegregation as well as a priori evidence such as deleteriousness and functional consequence. BICEP can correctly identify causal variants for phenotypes with both Mendelian and complex genetic architectures, outperforming existing methodologies. Additionally, BICEP can correctly down-weight common variants that are unlikely to be involved in phenotypic liability in the context of a pedigree, even if they have reasonable cosegregation patterns. The output metrics from BICEP allow for the quantitative comparison of variant causality within and across pedigrees, which is not possible with existing approaches.

Pedigree↗

Coeliac disease: investigation of proposed causal variants in the CTLA4 gene region.

Coeliac disease (CD) is an immune-mediated enteropathy triggered by gluten in genetically predisposed individuals. Patients with CD have an increased prevalence of other autoimmune disorders, including type 1 diabetes (T1D) and Graves' disease (GD). CD shares with these conditions certain HLA susceptibility alleles. A number of studies have also shown association of autoimmune diseases, including CD, with the CD28-cytotoxic T lymphocyte antigen 4 (CTLA4)-inducible costimulator (ICOS) region of chromosome 2q33, but until recently the precise causal variant has remained unknown. Recently, it was shown that, in GD, CT60 (+6230G>A), a single nucleotide polymorphism (SNP) at the end of the CTLA4 transcript, is associated with an alteration in the ratio of splice forms of the CTLA4 gene and that this ratio affects disease susceptibility. A similar but weaker association was found with T1D. There is also an independent association of GD and T1D with the SNP MH30 (-23 327G>C), which possibly affects promoter region function. Hypothesizing that CT60 and MH30 may be causal variants in other autoimmune disorders, we investigated these SNPs in CD using 149 family trios and 100 unrelated/unaffected controls. No association was detected with either SNP using both the transmission disequilibrium test (TDT) and case-control methods. Our study appears to have good power to detect moderate genetic effects, but possibly these SNPs exert too weak an effect on risk of CD to have been detected in our sample. Alternatively, the previously noted association of CD with the CTLA4 gene region may be due to different causal variants. Unlike T1D and GD, CD is not a true autoimmune disease, and CD has different associations at the CTLA4 exon 1 SNP +49G>A from all other autoimmune disorders. MH30, CT60, and other SNPs in the region may still warrant further investigation in other CD samples.

Antigens, CD↗

Identifying candidate causal variants responsible for altered activity of the ABCB1 multidrug resistance gene.

The difficulty of fine localizing the polymorphisms responsible for genotype-phenotype correlations is emerging as an important constraint in the implementation and interpretation of genetic association studies, and calls for the definition of protocols for the follow-up of associated variants. One recent example is the 3435C>T polymorphism in the multidrug transporter gene ABCB1, associated with protein expression and activity, and with several clinical conditions. Available data suggest that 3435C>T may not directly cause altered transport activity, but may be associated with one or more causal variants in the poorly characterized stretch of linkage disequilibrium (LD) surrounding it. Here we describe a strategy for the follow-up of reported associations, including a Bayesian formalization of the associated interval concept previously described by Goldstein. We focus on the region of high LD around 3435C>T to compile an exhaustive list of variants by (1) using a relatively coarse set of marker typings to assess the pattern of LD, and (2) resequencing derived and ancestral chromosomes at 3435C>T through the associated interval. We identified three intronic sites that are strongly associated with the 3435C>T polymorphism. One of them is associated with multidrug resistance in patients with epilepsy (chi2 = 3.78, P = 0.052), and sits within a stretch of significant evolutionary conservation. We argue that these variants represent additional candidates for influencing multidrug resistance due to P-glycoprotein activity, with the IVS 26+80 T>C being the best candidate among the three intronic sites. Finally, we describe a set of six haplotype tagging single-nucleotide polymorphisms that represent common ABCB1 variation surrounding 3435C>T in Europeans.

ATP Binding Cassette Transporter, Subfamily B, Mem↗

Fine mapping of the MAPT locus using quantitative trait analysis identifies possible causal variants in Alzheimer's disease.

In addition to senile plaques, neurofibrillary tangles are characteristic of Alzheimer's disease (AD) pathology, suggesting a clear involvement of the microtubule-associated protein tau (MAPT) in AD. Recent findings, suggesting that the H1c haplotype is associated with increased risk, now also implicate MAPT genetically. In this study, we aim to clarify this association by a fine mapping approach using both a traditional phenotypic association analysis and a quantitative trait (QT) analysis using cerebrospinal fluid (CSF) tau protein levels in the German population. Here, we report that both methodologies identify that the H1c haplotype may play important role in AD (AD risk, P=0.007, uncorrected; CSF tau levels, P=0.027, uncorrected). Further, the use of a sliding window approach in the QT analysis allowed for the narrowing down of the region where a probable causal variant may be located. The data suggest that this may lie at or within close proximity to the rs242557 single nucleotide polymorphism as association with CSF tau levels seems to be primarily driven by rs242557 in a gene dosage-dependent manner (trend model: P=0.002, uncorrected). These findings provide functional evidence to support the genetic association of MAPT with AD.

Aged↗

Towards improved fine-mapping of candidate causal variants.

Fine-mapping in genome-wide association studies aims to identify potentially causal genetic variants among a set of candidate variants that are often highly correlated with each other owing to linkage disequilibrium. A variety of statistical approaches are used in fine-mapping, almost all of which are based on a multiple regression framework to model the relationship between genotype and phenotype, while accommodating specific assumptions about the distribution of variant effect sizes and using different inference algorithms. Owing to their modelling flexibility and the ease of making inferential statements, these approaches are predominantly Bayesian in nature. Recently, these approaches have been improved by refining modelling assumptions, integrating additional information, accommodating summary statistics, and developing scalable computational algorithms that improve computation efficiency and fine-mapping resolution.

Humans↗

Fine-mapping the CYP2A6 regional association with nicotine metabolism among African American smokers.

The nicotine metabolite ratio (NMR; 3'hydroxycotinine/cotinine) is a stable biomarker for CYP2A6 enzyme activity and nicotine clearance, with demonstrated clinical utility in personalizing smoking cessation treatment. Common genetic variation in the CYP2A6 region is strongly associated with NMR in smokers. Here, we investigated this regional association in more detail. We evaluated the association of CYP2A6 single-nucleotide polymorphisms (SNPs) and * alleles with NMR among African American smokers (N = 953) from two clinical trials of smoking cessation. Stepwise conditional analysis and Bayesian fine-mapping were undertaken. Putative causal variants were incorporated into an existing African ancestry-specific genetic risk score (GRS) for NMR, and the performance of the updated GRS was evaluated in both African American (n = 953) and European ancestry smokers (n = 933) from these clinical trials. Five independent associations with NMR in the CYP2A6 region were identified using stepwise conditional analysis, including the deletion variant CYP2A6*4 (beta = -0.90, p = 1.55 × 10-11). Six putative causal variants were identified using Bayesian fine-mapping (posterior probability, PP = 0.67), with the top causal configuration including CYP2A6*4, rs116670633, CYP2A6*9, rs28399451, rs8192720, and rs10853742 (PP = 0.09). Incorporating these putative causal variants into an existing ancestry-specific GRS resulted in comparable prediction of NMR within African American smokers, and improved trans-ancestry portability of the GRS to European smokers. Our findings suggest that both * alleles and SNPs underlie the association of the CYP2A6 region with NMR among African American smokers, identify a shortlist of variants that may causally influence nicotine clearance, and suggest that portability of GRSs across populations can be improved through inclusion of putative causal variants.

Adult↗

Functional mapping and annotation of genetic associations with FUMA.

A main challenge in genome-wide association studies (GWAS) is to pinpoint possible causal variants. Results from GWAS typically do not directly translate into causal variants because the majority of hits are in non-coding or intergenic regions, and the presence of linkage disequilibrium leads to effects being statistically spread out across multiple variants. Post-GWAS annotation facilitates the selection of most likely causal variant(s). Multiple resources are available for post-GWAS annotation, yet these can be time consuming and do not provide integrated visual aids for data interpretation. We, therefore, develop FUMA: an integrative web-based platform using information from multiple biological resources to facilitate functional annotation of GWAS results, gene prioritization and interactive visualization. FUMA accommodates positional, expression quantitative trait loci (eQTL) and chromatin interaction mappings, and provides gene-based, pathway and tissue enrichment results. FUMA results directly aid in generating hypotheses that are testable in functional experiments aimed at proving causal relations.

Chromatin↗

Refined genomic localization and ethnic differences observed for the IBD5 association with Crohn's disease.

Although the general association of the inflammatory bowel disease (IBD) 5 region on chromosome 5q31 to Crohn's disease (CD) has been replicated repeatedly, the identity of the precise causal variant within the region remains unknown. A recent report proposed polymorphisms in solute carrier family 22, member 4 (SLC22A4) organic cation transporter 1(OCTN1) and solute carrier family 22, member 5 (SLC22A5) (OCTN2) as responsible for the IBD5 association, but definitive, large-sample comparison of those polymorphisms with others known to be in strong linkage disequilibrium was not performed. We evaluated 1879 affected offspring and parents ascertained by a North American IBD Genetics Consortium for six IBD5 tag single nucleotide polymorphisms (SNPs) to evaluate association localization and ethnic and subphenotypic specificity. We confirm association to the IBD5 region (best SNP IGR2096a_1/rs12521868, P<0.0005) and show this association to be exclusive to the non-Jewish (NJ) population (P=0.00005) (risk allele undertransmitted in Ashkenazi Jews). Using Phase II HapMap data, we demonstrate that there are a set of polymorphisms, spanning genes from prolyl 4-hydroxylase (P4HA2) through interferon regulatory factor 1 (IRF1) with equivalent statistical evidence of association to the reported SLC22A4 variant and that each, by itself, could entirely explain the IBD5 association to CD. Additionally, the previously reported SLC22A5 SNP is rejected as the potential causal variant. No specificity of association was seen with respect to disease type and location, and a modest association to ulcerative colitis is also observed. We confirm the importance of IBD5 to CD susceptibility, demonstrate that the locus may play a role in NJ individuals only, and establish that IRF1, PDLIM, and P4HA2 may be equally as likely to contain the IBD5 causal variant as the OCTN genes.

Adult↗

Association of genetically proxied cancer-targeted drugs with cardiovascular diseases through Mendelian randomization analysis.

BACKGROUND: Cancer-targeted therapies are progressively pivotal in oncological care. Observational studies underscore the emergence of cancer therapy-related cardiovascular toxicity (CTR-CVT), impacting patient outcomes. We aimed to investigate the causal relationship between different types of cancer-targeted therapies and cardiovascular disease (CVD) outcomes through a two-sample Mendelian randomization (MR) study. METHODS: This genome-wide association study was conducted using a two-sample Mendelian randomization framework. Genetic instruments for drug target gene expression were extracted from the eQTLGen consortium (31684 individuals, 37 cohorts). Genome-wide association study (GWAS) summary statistics for 19 cardiovascular diseases were derived from the FinnGen database. Primary analysis was carried out using the summary-data-based MR (SMR) method, with sensitivity analysis for validation. Colocalization analysis identifies shared causal variants between exposure eQTLs and CVD-associated single-nucleotide polymorphisms (SNPs). RESULTS: Among the 39 drug target genes, 8 were identified with detectable cis-eQTLs and were subsequently validated through positive control analysis for further investigation. In the SMR and sensitivity analyses, genetically proxied VEGFA inhibition showed significantly strong association with stroke (odds ratio [OR]&#x2009;=&#x2009;1.17, 95% confidence interval [CI]&#x2009;=&#x2009;1.09-1.26, p&#x2009;=&#x2009;1.33&#x2009;&#xd7;&#x2009;10-&#x2009;5). Additionally, the inhibition of FGFR1, FLT1, and MAP2K2 exhibited suggestive association with corresponding cardiovascular disease outcomes. Nevertheless, only VEGFA expression and stroke shared a causal variant (93.6%), whereas FGFR1, MAP2K2, and FLT1 did not share causal variants with corresponding cardiovascular diseases in the colocalization analysis. CONCLUSIONS: This genetic association study revealed evidence supporting the genetic association between the use of VEGFA inhibitors and increased stroke risk, highlighting the need for enhanced pharmacovigilance. These findings underscore the delicate balance between cardiovascular toxicity risk and the benefits of cancer-targeted therapy.

Humans↗

A flexible Bayesian framework for modeling haplotype association with disease, allowing for dominance effects of the underlying causative variants.

Multilocus analysis of single-nucleotide-polymorphism (SNP) haplotypes may provide evidence of association with disease, even when the individual loci themselves do not. Haplotype-based methods are expected to outperform single-SNP analyses because (i) common genetic variation can be structured into haplotypes within blocks of strong linkage disequilibrium and (ii) the functional properties of a protein are determined by the linear sequence of amino acids corresponding to DNA variation on a haplotype. Here, I propose a flexible Bayesian framework for modeling haplotype association with disease in population-based studies of candidate genes or small candidate regions. I employ a Bayesian partition model to describe the correlation between marker-SNP haplotypes and causal variants at the underlying functional polymorphism(s). Under this model, haplotypes are clustered according to their similarity, in terms of marker-SNP allele matches, which is used as a proxy for recent shared ancestry. Haplotypes within a cluster are then assigned the same probability of carrying a causal variant at the functional polymorphism(s). In this way, I can account for the dominance effect of causal variants, here corresponding to any deviation from a multiplicative contribution to disease risk. The results of a detailed simulation study demonstrate that there is minimal cost associated with modeling these dominance effects, with substantial gains in power over haplotype-based methods that do not incorporate clustering and that assume a multiplicative model of disease risks.

Algorithms↗

Admixture-mapping analysis reveals genetic determinants of the human plasma proteome.

Protein profiling and genetic findings can be integrated to define the genetic architecture of the circulating proteome in chronic diseases. Most self-identified African American (AA) individuals have both African and European genetic ancestry. Admixture mapping can detect genomic association regions in which causal variants exist with substantial differences in allele frequency or effect sizes between genetic ancestries. We performed admixture mapping of the circulating proteome in 1,989 participants from the Jackson Heart Study (JHS), investigating the relation of local African ancestry within genomic regions with levels of circulating proteins. We conditioned protein-local ancestry association models on variants previously found to be associated with those proteins in genome-wide association studies (GWASs). We replicated findings in 196 AA participants from the Multi-Ethnic Study of Atherosclerosis (MESA). 62 proteins were associated with local African ancestry. 21 of 62 remained statistically significant after conditioning on protein-associated variants observed in previous GWASs. 48 of 54 available protein-local ancestry associations were replicated in the MESA. Proteins associated with local African ancestry included chemokines, factors associated with vascular biology and inflammation, and other biologically interesting proteins. Admixture associations unexplained by previously reported protein-associated variants in conditional analysis suggest the existence of causal variants missed by standard GWAS techniques.

Aged↗

Genetic Architecture of Idiopathic Inflammatory Myopathies From Meta-Analyses.

OBJECTIVE: Idiopathic inflammatory myopathies (IIMs, myositis) are rare systemic autoimmune disorders that lead to muscle inflammation, weakness, and extramuscular manifestations, with a strong genetic component influencing disease development and progression. Previous genome-wide association studies identified loci associated with IIMs. In this study, we imputed data from two prior genome-wide myositis studies and analyzed the largest myositis data set to date to identify novel risk loci and susceptibility genes associated with IIMs and its clinical subtypes. METHODS: We performed association analyses on 14,903 individuals (3,206 patients and 11,697 controls) with genotypes and imputed data from the Trans-Omics for Precision Medicine reference panel. Fine-mapping and expression quantitative trait locus colocalization analyses in myositis-relevant tissues indicated potential causal variants. Functional annotation and network analyses using the random walk with restart (RWR) algorithm explored underlying genetic networks and drug repurposing opportunities. RESULTS: Our analyses identified novel risk loci and susceptibility genes, such as FCRLA, NFKB1, IRF4, DCAKD, and ATXN2 in overall IIMs; NEMP2 in polymyositis; ACBC11 in dermatomyositis; and PSD3 in myositis with anti-histidyl-transfer RNA synthetase autoantibodies (anti-Jo-1). We also characterized effects of HLA region variants and the role of C4. Colocalization analyses suggested putative causal variants in DCAKD in skin and muscle, HCP5 in lung, and IRF4 in Epstein-Barr virus (EBV)-transformed lymphocytes, lung, and whole blood. RWR further prioritized additional candidate genes, including APP, CD74, CIITA, NR1H4, and TXNIP, for future investigation. CONCLUSION: Our study uncovers novel genetic regions contributing to IIMs, advancing our understanding of myositis pathogenesis and offering new insights for future research.

Humans↗

Perils and pitfalls of permutation tests for distinguishing the effects of neighbouring polymorphisms.

In a small region several marker loci may be associated with a trait, either because they directly influence the trait or because they are in linkage disequilibrium (LD) with a causal variant. Having established a potentially causal effect at a primary variant, we may ask if any other variants in the region appear to further contribute to the trait, indicating that the additional variant is either causal or is in LD with another causal locus. Methods of approaching this problem using case-parent trio data include the stepwise conditional logistic regression approach described by Cordell and Clayton ([2002] Am. J. Hum. Genet. 70:124-141), and a constrained-permutation method recently proposed by Spijker et al. ([2005] Ann. Hum. Genet. 69:90-101). Through simulation we demonstrate that the procedure described by Spijker et al. [2005], as well as unconditional logistic regression with "affected family-based controls" (AFBACs), can lead to inflated type 1 errors in situations when haplotypes are not inferable for all trios, whereas the conditional logistic regression approach gives correct significance levels. We propose an alternative to the permutation method of Spijker et al. [2005], which does not rely on haplotyping, and results in correct type 1 errors and potentially high power when assumptions of random mating, Hardy-Weinberg Equilibrium, and multiplicative effects of disease alleles are satisfied.

Genetic Markers↗

CTLA4 +49 A/G and CT60 polymorphisms in Dutch coeliac disease patients.

Coeliac disease is an autoimmune disorder, characterised by villous atrophy of the small intestine, which results from a T-cell-mediated response to gluten-derived peptides. The cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) is involved in the regulation of T-cell activation and the CTLA4 +49 A/G polymorphism in exon 1 has been implicated in several autoimmune disorders, including coeliac disease. However, this polymorphism was recently excluded as being the causal variant in Graves' disease, autoimmune hypothyroidism and type I diabetes mellitus. This causal variant was mapped to the 3' region of CTLA4, with the CT60 polymorphism showing the strongest association. The aim of this study was to determine the role of the CTLA4 gene in coeliac disease in the Dutch population. The +49 A/G and CT60 polymorphisms were genotyped in a case-control cohort of 215 patients and controls. The frequency of the +49 G-allele was increased in cases, although not significantly. However, the frequency of the CT60 G-allele was increased with borderline significance in coeliac disease patients (P = 0.048), although the genotype distributions did not show a significant difference between cases and controls. These results indicate the involvement of the CTLA4 gene in coeliac disease development. The haplotype carrying the CT60 G-allele was shown to be associated with lower mRNA levels of the soluble CTLA-4 isoform, providing a possible mechanism for the T-cell-mediated destruction of the small intestine.

Adult↗

Identifying causal genetic variants for high-altitude adaptation through blood eQTL analysis in plateau populations.

A substantial number of genetic variants have been associated with high-altitude adaptation (HAA), yet most of them are located in non-coding genomic regions, leaving their specific functions and underlying mechanisms largely unknown. In this study, we analyze whole-genome and transcriptome sequencing data from a self-established cohort comprising 61 native highlanders (NHs) and 164 acclimatized newcomers (ANs), identifying 6,586 cis- and 34,203 trans-expression quantitative trait loci (eQTLs), along with 130 cell type-specific eQTLs. By further combining these data with a large East Asia (~30% Tibetan) genome-wide association study (GWAS) cohort, we employ colocalization and causal inference analyses to prioritize 85 cis-eQTLs associated with HAA and identify several novel candidate causal genes, including EXOC8, which is experimentally confirmed to regulate erythroid differentiation. Additionally, network analysis of these causal genes uncovers multiple regulatory pathways, mainly involving energy metabolism, autophagy, ubiquitination and inflammation. Our study offers a comprehensive eQTL map and reveals causal chains of "variant-gene-phenotype" for HAA-related traits, which provides new insights into potential regulatory mechanisms and targets for prevention and treatment of altitude sickness.

Quantitative Trait Loci↗

BTS: a scalable Bayesian Tissue Score for prioritizing GWAS variants and their functional contexts across >1000s of omics datasets.

MOTIVATION: statistics from genome-wide association studies (GWAS) are widely used in fine-mapping and colocalization analyses to identify causal variants and their enrichment in functional contexts, such as affected cell types and genomic features. With the expansion of functional genomic (FG) datasets, which now include hundreds of thousands of tracks across various cell and tissue types, it is critical to establish scalable algorithms integrating thousands of diverse FG annotations with GWAS results. RESULTS: We propose BTS (Bayesian Tissue Score), a novel, highly efficient algorithm uniquely designed for (i) identifying affected cell types and functional elements (context-mapping) and (ii) fine-mapping potentially causal variants in a context-specific manner using large collections of cell type-specific FG annotation tracks. BTS leverages GWAS summary statistics and annotation-specific Bayesian models to analyze genome-wide annotation tracks, including enhancers, open chromatin, and histone marks. We evaluated BTS on GWAS summary statistics for immune and cardiovascular traits, such as Inflammatory Bowel Disease (IBD), Rheumatoid Arthritis (RA), Systemic Lupus Erythematosus (SLE), and Coronary Artery Disease (CAD). Our results demonstrate that BTS is over 100&#xd7; more efficient in estimating functional annotation effects and context-specific variant fine-mapping compared to existing methods. Importantly, this large-scale Bayesian approach prioritizes both known and novel annotations, cell types, genomic regions, and variants and provides valuable biological insights into the functional contexts of these diseases. AVAILABILITY AND IMPLEMENTATION: Docker image is available at https://hub.docker.com/r/wanglab/bts with preinstalled BTS R package (https://bitbucket.org/wanglab-upenn/BTS-R) and BTS GWAS summary statistics analysis pipeline (https://bitbucket.org/wanglab-upenn/bts-pipeline).

Genome-Wide Association Study↗

Exploring the causal relationship between plasma proteins and postherpetic neuralgia: a Mendelian randomization study.

BACKGROUND: The proteome represents a valuable resource for identifying therapeutic targets and clarifying disease mechanisms in neurological disorders. This study investigated potential causal relationships between plasma proteins and postherpetic neuralgia (PHN). METHODS: We conducted a two-sample Mendelian randomization (MR) analysis using genome-wide association study (GWAS) summary statistics from the Decode Genetics dataset (4,907 plasma proteins) and the FinnGen database (490 PHN cases and 435,371 controls). Instrumental variables (IVs) were selected based on relevance, independence, and exclusivity. Causal associations were assessed using inverse-variance weighted (IVW), MR-Egger regression, simple mode, weighted mode, and weighted median methods. Sensitivity analyses, including leave-one-out tests, evaluated result robustness, while colocalization analysis examined shared causal variants between traits. RESULTS: Eight plasma proteins showed significant associations with PHN (PFDR < 0.05). Higher levels of ATRN, PIANP, and CD48 correlated with increased PHN risk, whereas elevated KIR2DL5A, GPI, SEMG2, EIF4B, and HFE2 levels were associated with reduced risk. Sensitivity analyses supported these findings and excluded genetic pleiotropy as a major confounding factor. Colocalization analysis did not detect shared causal variants (PPH4 < 0.8). CONCLUSION: These results suggest a potential causal role for eight plasma proteins in PHN pathogenesis. While these proteins may serve as biomarkers or therapeutic candidates, further validation is required. This study advances understanding of PHN pathophysiology and supports future investigations into diagnostic and therapeutic strategies.

Mendelian randomization↗

Single-cell expression quantitative trait locus Mendelian randomization reveals immune cell-specific causal regulatory networks and actionable targets in polycystic ovary syndrome.

ObjectiveTo systematically investigate whether the pathogenesis of polycystic ovary syndrome (PCOS) is causally related to dysregulated gene expression in specific immune cell subsets, and to evaluate the potential of these causal genes as actionable drug targets.MethodsThis study employed a two-sample Mendelian randomization (MR) framework using publicly available genome-wide association study (GWAS) summary statistics. The participant data included 797 PCOS cases and 140,558 controls (no direct patient recruitment was involved). Instrumental variables were derived from high-resolution immune cell-specific single-cell expression quantitative trait locus (sc-eQTL) data (OneK1K project) across 14 immune cell types. Primary analyses utilized the inverse-variance weighted (IVW) method. Shared causal variants were validated using Bayesian colocalization. Phenome-wide association analysis (PheWAS), external transcriptomic dataset validation (GSE8157), and DrugBank database screening were conducted for pleiotropy assessment and drug repositioning.ResultsMR analysis revealed genome-wide significant causal associations for GLIPR1 in non-classical monocytes (Mono NC) and XBP1 in CD4+ effector memory T cells (CD4 ET) with PCOS risk. Higher GLIPR1 expression was associated with a decreased PCOS risk (OR = 0.669, P = 4.34&#xd7;10-6), whereas higher XBP1 expression was associated with an increased risk (OR = 1.406, P = 9.53&#xd7;10-8). Colocalization analysis confirmed that GLIPR1 shares a causal variant with PCOS (PP.H4 = 96.73%). PheWAS and external validation confirmed the safety profile and significant upregulation (P = 0.03) of GLIPR1. Drug repositioning identified SOT-107, a Phase III protein therapy drug, as a potential interacting agent for GLIPR1.ConclusionsThis sc-eQTL MR study reveals immune cell-specific causal regulatory networks in PCOS. GLIPR1 in non-classical monocytes represents a high-confidence protective target, while XBP1 provides suggestive evidence for immune-mediated pathogenesis. The candidate drug SOT-107 highlights theoretical repositioning opportunities, though rigorous preclinical validation remains required.

Female↗