Search PubMedSearch

SEARCH · Search PubMed

Results for “Colocalization analysis”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

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

At least 55 records · Page 3Linked to original sources

Proteins containing an uncleaved signal for glycophosphatidylinositol membrane anchor attachment are retained in a post-ER compartment.

Glycophosphatidylinositol (GPI)-anchored membrane proteins are initially synthesized with a cleavable COOH-terminal extension that signals anchor attachment. Overexpression in COS cells of hGH-DAF fusion proteins containing the GPI signal of decay accelerating factor (DAF) fused to the COOH-terminus of human growth hormone (hGH), produces both GPI-anchored hGH-DAF and uncleaved precursors that retain the GPI signal. Using hGH-DAF fusion proteins containing a mutated, noncleavable GPI signal, we show that uncleaved polypeptides are retained inside the cell and accumulate in a brefeldin A-sensitive, Golgi-like juxtanuclear structure. Retention requires the presence of either a functional or a noncleavable GPI signal; hGH-DAF fusion proteins containing only the COOH-terminal hydrophobic domain (a component of the GPI signal) are secreted. Immunofluorescence analysis shows colocalization of the retained, uncleaved fusion proteins with both a Golgi marker and with p53, a marker of the ER-Golgi intermediate compartment. Since N-linked glycosylation is postulated to facilitate the transport of proteins to the cell surface, we engineered a glycosylation site into hGH-DAF. Glycosylation failed to completely override the transport block, but allowed some uncleaved hGH-DAF to pass through the secretory pathway and acquire endoglycosidase H resistance. The retained molecules remained endoglycosidase H sensitive. We suggest that the uncleaved fusion protein is retained in a sorting compartment between the ER and the medial Golgi complex. We speculate that a mechanism exists to retain proteins containing an uncleaved GPI signal as part of a system for quality control.

Amino Acid Sequence

Dissecting the shared genetic architecture between migraine subtypes and cardiovascular diseases: a multi-layered genomic analysis.

BACKGROUND: Epidemiological studies have linked migraine to an increased risk of cardiovascular disease (CVD); however, the shared genetic basis and putative causal relationships between migraine subtypes and cardiovascular traits remain poorly understood. METHODS: Leveraging large-scale GWAS summary statistics for migraine phenotypes (overall migraine, migraine with aura [MA], and migraine without aura [MO]) from FinnGen R12, along with seven cardiovascular diseases from publicly available consortia, we conducted a multi-layered genetic analysis. This integrative framework encompassed genetic correlation [linkage disequilibrium score regression (LDSC) and high-definition likelihood (HDL)], cross-trait meta-analysis (CPASSOC and PLACO), Bayesian colocalization, summary-data-based Mendelian randomization (SMR) using GTEx v8 eQTL data, and bidirectional two-sample Mendelian randomization (MR). RESULTS: Significant genetic correlations were identified between migraine and multiple cardiovascular traits, with hypertension and coronary artery disease (CAD) showing the most robust associations. MA exhibited broader genetic overlap with cardiovascular diseases than MO, including a notably stronger correlation with ischemic stroke, whereas MO demonstrated a stronger correlation with hypertension. Cross-trait meta-analysis identified 160 pleiotropic loci across 17 of 21 trait pairs. Colocalization analysis confirmed 32 loci harboring shared causal variants, mapped to 13 candidate genes, of which 7 (PHACTR1, LRP1, SOX7, ABO, FHOD3, MEI1, XKR6) were further validated by SMR as exhibiting tissue-specific regulatory effects. Among these, PHACTR1 displayed the broadest pleiotropic profile across migraine phenotypes and vascular diseases. After MR-PRESSO outlier removal, bidirectional MR identified 10 MR-supported associations, two of which (genetic liability to hypertension on overall migraine, and CAD on MA) survived Bonferroni correction, all free of detectable horizontal pleiotropy. Genetic liability to hypertension was associated with increased migraine risk (OR = 1.90, 95% CI 1.25-2.90, P = 2.64 × 10⁻³), atherosclerotic diseases showed subtype-specific effects (inverse for MO, positive for MA), and, in the reverse direction, migraine was associated with increased ischemic stroke risk. CONCLUSIONS: This study provides a comprehensive and systematic characterization of the shared genetic architecture between migraine subtypes and cardiovascular diseases. By identifying pleiotropic genes and bidirectional putative causal relationships with subtype-specific patterns, our findings carry implications for the development of targeted therapeutics and subtype-specific cardiovascular risk stratification.

Humans

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×10-6), whereas higher XBP1 expression was associated with an increased risk (OR = 1.406, P = 9.53×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

Immunohistochemical and biochemical analysis of serotonin and substance P colocalization in the nucleus tractus solitarii and associated afferent ganglia of the rat.

In a previous study of afferent projections to the nucleus tractus solitarii (NTS), it was shown that over half of the retrogradely-labelled neurons in the nucleus raphe pallidus contained serotonin-immunoreactivity and over half of these neurons contained substance P-immunoreactivity, suggesting that these two putative neurotransmitters are colocalized in NTS-afferent neurons. The objectives of the present study were to 1) directly determine if varicosities in the NTS, the area postrema (AP), and the dorsal motor nucleus of the vagus nerve (DMN) do contain both transmitters, 2) determine if primary afferent neurons in the nodose and pretrosal ganglia might also colocalize serotonin and substance P, and 3) quantify the amount of substance P that is contained in serotonergic varicosities in the NTS. Distributions and colocalization of substance P and serotonin in the NTS were studied using dual-color immunohistochemistry, while the quantity of substance P in serotonergic varicosities was assessed by radioimmunoassay (RIA) using micropunches from the NTS of 5,7-dihydroxytryptamine-(5,7 DHT-) and vehicle-treated rats. Varicosities that contained both serotonin- and substance P-immunoreactivity were found in the NTS, the DMN, and the AP. Double-labelled varicosities were common in the NTS and DMN (i.e., qualitatively similar to the density seen in the hypoglossal nucleus and in the ventral horn of the cervical spinal cord); however, the vast majority of the varicosities in these autonomic areas only displayed immunoreactivity for one or the other of these transmitters. This paucity of doubly-labelled varicosities, in comparison to the number of singly-labelled varicosities, was reflected in the lack of a significant decrease in substance P levels as determined by RIA of micropunches taken from caudal and intermediate levels of the NTS in 5,7 DHT- and vehicle-treated rats.(ABSTRACT TRUNCATED AT 250 WORDS)

Afferent Pathways

From genetic causality to druggable targets: A multiomics framework identifies ZSCAN16 in gout pathogenesis.

ObjectiveGout is a prevalent form of inflammatory arthritis in which many patients respond suboptimally to current therapies. Drug development is hampered by a lack of genetically validated targets, leading to high clinical trial attrition. This study aimed to systematically identify and prioritize novel, druggable targets for gout via a multilayered genetic and functional genomics approach.MethodsWe performed two-sample Mendelian randomization (MR) using cis-expression quantitative trait locus (cis-eQTL) data and dual independent gout genome-wide association study (GWAS) cohorts (openGWAS and FinnGen). The candidate genes were subjected to a rigorous validation pipeline including Bayesian colocalization, phenome-wide association studies (PheWASs) to assess pleiotropy and on-target safety, and single-cell RNA sequencing (scRNA-seq) to delineate the cellular context. Molecular docking was used to evaluate the structural druggability of prioritized targets.ResultsMR analysis revealed 15 genes causally associated with gout. Colocalization analysis (PPH4 > 0.8) prioritized two targets: ZSCAN16 (risk-increasing, OR = 1.04, 95% CI [1.02-1.06]) and TRIM10 (protective, OR = 0.96, 95% CI [0.94-0.98]). Crucially, PheWAS revealed that ZSCAN16 is highly specific to gout, whereas TRIM10 exhibited extensive pleiotropy with hematological and cardiometabolic traits, indicating significant safety risks. Single-cell analysis provided orthogonal validation, demonstrating flare-specific upregulation of ZSCAN16 in cytotoxic T/NK cells. Molecular docking confirmed ZSCAN16 as a structurally druggable target, showing high-affinity binding with known compounds (e.g. digoxin, binding energy = -9.6 kcal/mol).ConclusionsOur study identifies ZSCAN16 as a high-potential, druggable therapeutic target for gout, highlighting its genetic influence on specific immune cell activities during acute flares. Conversely, TRIM10 was deprioritized owing to substantial pleiotropic liabilities and poor chemical tractability. These findings suggest that ZSCAN16 could play a crucial role in the pathogenesis of gout and may provide a valuable lead for future drug discovery efforts.

Humans

Detection of adrenocorticotropin-related and alpha-melanocyte-stimulating hormone-related substances in the anterior pituitary of larval and adult Ambystoma tigrinum (class: Amphibia).

Immunohistochemical analysis revealed the colocalization of ACTH-related immunoreactive forms and alpha-MSH-related immunoreactive forms in corticotropic cells of the anterior pituitary of larval Ambystoma tigrinum. Subsequent analysis of acid extracts of anterior pituitaries obtained from larval, neotene, and postmetamorphic adults indicated that the molar ratio of ACTH(1-39) to alpha-MSH was, respectively, 2.5:1, 1:1.5, and 1:1.3. However, in sexually mature adults the molar ratio of ACTH(1-39) to alpha-MSH was 8:1. These results indicate that before and immediately after metamorphosis, both ACTH(1-39) and alpha-MSH are major end products of corticotropic cells in this species. As postmetamorphic animals reach sexual maturity, ACTH(1-39) is a major end product of corticotropic cells and only trace amounts of alpha-MSH are produced. Thus in A. tigrinum the proteolytic processing of ACTH in corticotropic cells is developmentally regulated. The alpha-MSH-related peptide produced by the corticotropic cells appears to be an ACTH(1-13)amide-like form.

Adrenocorticotropic Hormone

Single-Cell Transcriptome-Wide Mendelian Randomization and Colocalization Uncover Potential Immunocytes-Related Therapeutic Targets for Obesity.

Weight-loss treatment is crucial for individuals with obesity to prevent various complications. The role of Immune cells in obesity has been recently recognized, whereas its translation into therapy requires identifying key target genes. We performed Mendelian randomization (MR) analysis to assess causal relationships between expression quantitative trait loci (eQTL) of 14 immune cells and obesity-related traits (obesity, body mass index and body fat percentage), and validated the results in colocalization analysis. For the putative causal genes identified by the MR and colocalization analyses, we conducted pathway enrichment, differential expressed gene (DEG) analysis and search of druggable evidence, and utilized a Tier system to prioritize drug targets for obesity. MR and colocalization evidence was observed for 1630 genes associated with one or more obesity-related traits, mainly expressed in CD4+ naive/central memory T cells and enriched in antigen processing and presentation pathways. Forty-one genes showed causal relationship with all three outcomes, among which 19 genes have not been reported for obesity previously. DEG analysis using single-cell RNA sequencing data of blood or adipose tissue indicated that the differential expression of UBE2Z in monocytes, ZCCHC7 in T cells, and FNBP4 in B cells between lean and obese individuals were consistent with the MR results. By searching drug-gene interaction databases, we found targeted drugs for PYGB and PRUNE1, and PYGB was the top gene ranked in the Tier system. This study provides evidence for the involvement of immune cells in obesity, and the potential cell-specific, immune-related targets for obesity treatment.

Obesity

Integrated multi-omics strategies for identifying novel therapies in psoriasis.

MOTIVATION: Psoriasis is a chronic, immune-mediated disorder with an unmet need for effective treatments. To systematically prioritize therapeutic targets, we integrated proteome-wide Mendelian randomization (MR) with expression validation in blood/skin, genetic susceptibility analysis, differential gene expression (DGE) from bulk and single-cell RNA sequencing (scRNA-seq), colocalization, pathway enrichment, and protein-protein interaction analyses. RESULTS: Proteome-wide MR identified 29 candidate protein targets (Bonferroni-corrected), all replicated in independent datasets. Fifteen targets showed significant expression associations in blood or skin. Eleven proteins-UBLCP1, IL23A, ASF1A, RARRES2, ICAM1, PRSS53, ICAM5, GCA, IL2RA, DBI, and NFKB1-exhibited consistent directional effects with their genes. Genetic susceptibility analysis confirmed 20 target-specific polygenic scores for psoriasis and five for psoriatic arthritis. DGE analysis identified 13 targets in bulk and 13 in scRNA-seq-primarily in keratinocytes and immune cells-with IL2RA, COMP, and A2ML1 dysregulated across both. Colocalization analysis implicated shared causal variants for psoriasis in ASF1A, CD8A, CTF1, IL7R, MMP12, RARRES2, XCL2, DBI, IL23A, IL2RA, SGSH, and TIMD4. Enrichment analyses highlighted involvement in cytotoxicity, immune regulation, and JAK-STAT signaling. Eighteen targets interacted with approved anti-psoriasis drugs. Notably, drugs targeting IL2RA, IL7R, CTF1, ICAM1, MMP12, NFKB1, CD8A, DDX58, IL12A, SGSH, and FAP are approved or in trials for other diseases, suggesting repurposing potential. Our integrative multi-omics approach prioritized 29 high-confidence targets, including 13 novel candidates (RARRES2, ASF1A, CTF1, DBI, B3GNT2, CD8A, TIMD4, CRTAM, SGSH, XCL2, DAPK2, A2ML1, and FAP). Several high-priority targets-such as IL2RA, IL23, MMP12, RARRES2, IL7R, and ICAM1-were supported across analytical layers. These findings provide a robust foundation for psoriasis drug development. AVAILABILITY AND IMPLEMENTATION: The code used for the analyses in this manuscript has been archived in Zenodo at [DOI: 10.5281/zenodo.19692128].

Psoriasis

Neuropathologic findings in surgically treated hemimegalencephaly: immunohistochemical, morphometric, and ultrastructural study.

Tissues from three cases of hemimegalencephaly (HME) causing intractable seizures treated by cortical resection were studied using immunohistochemical, ultrastructural, and morphometric techniques. Severe cortical dysplasia was seen in all cases and included lesions best characterized as hemilissencephaly and polymicrogyria. Blurring of the cortex-white matter junction, the presence of large neuronal heterotopias, and neuronal cytomegaly were frequent observations. Immunohistochemical analysis demonstrated cellular colocalization of astrocytic markers glial fibrillary acidic protein and vimentin in one case of hemilissencephaly. Morphometric data showed significant increases over controls in neuronal profile area in all cases of HME. Neuronal cell density was increased significantly above controls in one of the cases. The study shows that HME results from severe cortical dysplasia which may be caused by multiple insults, manifest in one of several ways, and reflects abnormal or altered signals that regulate cortical morphogenesis.

Brain

Genetically Predicted Gene Expression and Circulating Metabolites Associated with Cervical High-Grade Squamous Intraepithelial Lesion: A Mendelian Randomization Study.

BACKGROUND: High-grade squamous intraepithelial lesion (HSIL) is a precancerous condition of the cervix. Identifying risk factors associated with HSIL and understanding their potential mechanisms may inform prevention strategies. This study aimed to investigate the associations of genetically predicted gene expression and circulating metabolites with HSIL risk using Mendelian randomization (MR). METHODS: We performed two-sample MR analysis to evaluate the associations of genetically predicted gene expression (eQTLGen consortium, N=31,684) and circulating metabolites (genome-wide association study [GWAS], N=8,299) with HSIL risk (FinnGen R12, N=293,218; 8,291 cases). Mediation analysis was conducted to explore whether metabolites might mediate the associations between genes and HSIL. Sensitivity analyses, including Mendelian randomization pleiotropy residual sum and outlier (MR-PRESSO), leave-one-out, and colocalization, were performed to assess the robustness of the findings. All GWAS data used in this study were derived from European-ancestry populations. RESULTS: Eleven genes showed significant associations with HSIL after false discovery rate (FDR) correction (q<0.05), including VWA7, PAX8, GUSBP1, IKZF3, PAX8-AS1, NFKBIL1 (interpret with caution due to an influential single nucleotide polymorphism [SNP]), ERBB2, COL11A2, SKIV2L, TCF19, and PGAP3. Eleven circulating metabolites were also significantly associated with HSIL. Mediation analysis suggested that two phospholipid metabolites (GCST90200685 and GCST90200692) might mediate a small proportion of the total protective association of COL11A2 with HSIL (1.46% and 1.45%, respectively), indicating that the protective association of COL11A2 is largely independent of these circulating metabolites. Colocalization analysis showed strong evidence of shared causal variants for eight genes (PP.H4>0.98), while COL11A2 showed weak evidence of colocalization (PP.H4=1.58&#xd7;10-15). Functional enrichment analysis indicated that COL11A2-related genes were enriched in extracellular matrix (ECM)-receptor interaction and PI3K-Akt signaling pathways. CONCLUSION: This MR study identified 11 genes and 11 circulating metabolites associated with HSIL risk. Among these, COL11A2 showed a protective association that appeared to be largely independent of circulating phospholipid metabolites, suggesting potential local mechanisms. These findings provide genetic and metabolic clues for future studies on HSIL etiology.

COL11A2

Prioritizing Parkinson's disease risk-associated mitochondrial candidate genes via multi-omics integrative analysis.

BACKGROUND: Mitochondrial dysfunction has been implicated in Parkinson's disease (PD), but the genetically regulated mitochondrial genes associated with PD risk remain incompletely defined. METHODS: We conducted a summary-data-based genetic epidemiology study integrating summary-based Mendelian randomization (SMR), Heterogeneity in dependent instruments (HEIDI) filtering, and Bayesian colocalization to prioritize mitochondrial-related molecular features associated with PD risk. Mitochondrial-related genes were defined using MitoCarta3.0. Genetically predicted gene expression and plasma protein abundance were evaluated using expression quantitative trait loci (eQTL) data from eQTLGen and GTEx v8, and protein quantitative trait loci (pQTL) data was assessed using International Parkinson's Disease Genomics Consortium (IPDGC) as the discovery genome-wide association study (GWAS) and FinnGen as the replication dataset. Prespecified QTL analyses were interpreted using FDR correction, HEIDI filtering, and colocalization support. DNA methylation QTL analysis, mitochondrial phenotype MR, and single-nucleus RNA-seq analysis were performed as complementary analyses. RESULTS: In the primary eQTL analysis, higher genetically predicted TTC19 expression was associated with lower PD risk (OR = 0.80, 95% CI: 0.74-0.87, PPH4&#x202f;= 0.80), whereas higher MALSU1 expression was associated with increased PD risk (OR = 2.21, 95% CI: 1.59-3.06, PPH4&#x202f;= 0.96). Both associations survived FDR correction, passed HEIDI filtering, and showed colocalization support. GTEx whole-blood data supported the direction of the TTC19 association. No mitochondrial protein reached significance after FDR correction and colocalization filtering in the primary pQTL analysis. Complementary methylation analysis highlighted cg06270993 as an exploratory regulatory signal for MALSU1. CONCLUSIONS: This MR-colocalization study prioritizes TTC19 and MALSU1 as genetically supported mitochondrial-related candidate genes associated with PD risk. Further validation is required to define their functional roles in PD pathogenesis.

Humans

Multi-omics integration and colocalization analyses prioritize candidate molecular loci associated with hypothermia.

BACKGROUND: Hypothermia is a life-threatening condition lacking specific pharmacological treatments. This study aimed to prioritize genetically supported molecular loci associated with hypothermia and to explore their pharmacological tractability using multi-omics data. METHODS: Initially, 2532 druggable genes were curated from the Drug-Gene Interaction Database and established literature. These were cross-referenced with cis-eQTL and cis-pQTL datasets, encompassing 870,655 and 114,281 SNPs for blood, respectively, alongside 2379 shared SNPs across adipose, skeletal muscle, and heart tissues. Matched instrumental variables were integrated with hypothermia GWAS summary statistics for two-sample Mendelian randomization (MR) and Bayesian colocalization. Transcriptomic differential expression analysis (DEA) was subsequently conducted as an exploratory analysis of cold-exposure-associated expression changes. Database-derived compound annotations were systematically re-evaluated according to target specificity, established pharmacological mechanism, and concordance with the direction of the MR estimates. RESULTS: Among 671 gene-level MR tests, 36 genes reached nominal significance, whereas only ABCC8 remained significant after FDR correction. Colocalization was evaluable for 8 of these 36 genes, and 4 loci (COL18A1, SLC1A7, ADIPOQ, and MERTK) met the prespecified PP.H4>0.90 threshold. The remaining 28 loci were not evaluable because sufficient overlapping regional variants were unavailable after harmonization. Transcriptomic analysis identified altered expression of SLC1A3 and SLCO4A1 under cold exposure, although these findings did not directly validate the colocalization-supported loci. Re-evaluation of database-derived compound annotations did not identify any direct, selective, and directionally concordant drug-repurposing candidate for hypothermia. CONCLUSIONS: COL18A1, SLC1A7, ADIPOQ, and MERTK showed colocalization support among the 8 evaluable nominal MR-associated loci. Because colocalization coverage was limited, these genes should be regarded as preliminary candidate loci rather than established therapeutic targets. The pharmacological annotations were indirect, non-selective, unsupported, or directionally inconsistent and should be interpreted solely as hypothesis-generating information.

Bayesian colocalization

Shared genetic architecture of obesity and gastroesophageal reflux disease.

Obesity is identified as a risk factor of gastroesophageal reflux disease (GERD). This study aims to elucidate the shared genetic architecture of obesity-related phenotypes and GERD. Based on the publicly available genome-wide association studies' datasets, this genome-wide pleiotropic association study was conducted with various genetic approaches (including linkage disequilibrium score regression, high-definition likelihood inference for genetic correlations, pleiotropic analysis under composite null hypothesis, Functional Mapping and Annotation, Bayesian colocalization, summary-based Mendelian randomization, and multi-marker analysis of genomic annotation analysis) sequentially to unravel the genetic associations from single-nucleotide polymorphism to gene levels, and to reveal the underlying shared genetic architecture between obesity-related phenotypes and GERD. This study discovered shared genetic mechanisms between GERD and several obesity-related phenotypes, including arm fat percentage (left), arm fat percentage (right), leg fat percentage (left), leg fat percentage (right), trunk fat percentage, waist-to-hip ratio, and body mass index. Significant genetic correlations were observed by linkage disequilibrium score regression and high-definition likelihood inference for genetic correlations, with multiple associated pleiotropic loci and their mapped genes identified by pleiotropic analysis under composite null hypothesis, Functional Mapping and Annotation, Bayesian colocalization, summary-based Mendelian randomization, and multi-marker analysis of genomic annotation analysis. Additionally, several brain tissues were identified to be linked to both obesity and GERD by multi-marker analysis of genomic annotation. This research provided strong evidence of genetic correlations and brought novel insights into the underlying genetic connections and shared genetic architectures of obesity and GERD.

Humans

Cell Type-Resolved Causal Inference and Spatial Transcriptomic Integration Reveal Immune-Specific Genetic Drivers of Autoimmune and Malignant Thyroid Disease.

BACKGROUND: Thyroid diseases, including autoimmune thyroid disease (AITD) and thyroid cancer, are characterized by immune dysregulation, yet the cell type-specific genetic mechanisms underlying these conditions remain poorly understood. Most genome-wide association studies (GWAS) have relied on bulk tissue expression quantitative trait loci (eQTL), which cannot resolve the heterogeneity of immune cell populations. METHODS: We performed two-sample Mendelian randomization (MR) analyses using single-cell cis-eQTLs from 14 immune cell subtypes (OneK1K cohort) as instrumental variables against GWAS summary statistics for four thyroid outcomes: autoimmune hyperthyroidism, autoimmune hypothyroidism, thyroid cancer and autoimmune thyroiditis. Causal associations were validated through Bayesian colocalization, phenome-wide association analysis (PheWAS) and multi-layered transcriptomic validation encompassing spatial transcriptomics of AITD tissue (GSE248205), bulk RNA-seq of thyroid cancer (GSE3678) and single-cell RNA-seq of thyroid tumours (GSE250521). gsMap spatial LD score regression was applied to map disease heritability onto spatial tissue architecture. RESULTS: We identified six Bonferroni-significant causal gene-cell type pairs for autoimmune hyperthyroidism, including protective effects of ABHD16A in na&#xef;ve/immature B cells (OR&#xa0;=&#xa0;0.440), HIST1H3H in CD8 NC T cells (OR&#xa0;=&#xa0;0.324), HMGN4 in NK recruiting cells (OR&#xa0;=&#xa0;0.556) and ZKSCAN4 in CD8 S100B T cells (OR&#xa0;=&#xa0;0.427), with five pairs showing strong colocalization (PP.H4 &#x2265; 86%). Three pairs reached significance for autoimmune hypothyroidism, including a risk association of HLA-F in CD4 NC T cells (OR&#xa0;=&#xa0;1.139). For autoimmune thyroiditis, FAM134B/RETREG1 showed consistent suggestive protective associations across both CD4 and CD8 NC T cells (PP.H4 &#x2265; 90% for both), suggesting a possible involvement of ER phagy regulation in thyroiditis susceptibility. Thyroid cancer showed a suggestive association with HLA-G in classical monocytes (OR&#xa0;=&#xa0;1.899, PP.H4&#xa0;=&#xa0;53%). Spatial transcriptomic validation demonstrated progressive immune infiltration from control tissue to Graves' disease to Hashimoto's thyroiditis (7.7%-15.7%, 46.1%-54.1%, respectively) and strong spatial correlation between target gene expression and corresponding cell type enrichment (e.g., plasma cell-HLA-DQB1: r&#xa0;=&#xa0;0.491, p < 10-300). HLA-G was independently validated in thyroid cancer bulk (log2fc&#xa0;=&#xa0;0.542, p&#xa0;=&#xa0;9.51&#xa0;&#xd7;&#xa0;10-3, AUC&#xa0;=&#xa0;0.857) and single-cell datasets. PheWAS revealed no significant associations detected for the core candidates. gsMap identified significant enrichment of autoimmune hypothyroidism heritability in gastrointestinal tract, adrenal gland and adipose tissue (all Bonferroni p < 0.002). CONCLUSIONS: This study establishes a multi-scale analytical framework integrating cell type-resolved genetic inference with spatial tissue validation, revealing distinct immunogenetic architectures underlying autoimmune versus malignant thyroid disease. Protective genetic programs in autoimmune hyperthyroidism converge on chromatin remodelling (HIST1H3H, HMGN4, ZKSCAN4) and lipid metabolism (ABHD16A) across lymphocyte subsets, whereas thyroid cancer risk involves immune escape mediated by HLA-G in myeloid cells. The ER-phagy receptor RETREG1 represents a candidate pathway warranting further investigation in autoimmune thyroiditis. These findings provide genetically supported, cell type-specific therapeutic targets and demonstrate a generalizable strategy for dissecting the immune-mediated mechanisms of complex thyroid diseases.

Mendelian randomization

Biochemical analysis of the TPS-b subfamily reveals a cineole-centered monoterpene biosynthetic module in Medicago truncatula.

Terpenoids constitute one of the largest and most structurally diverse classes of plant specialized metabolites, with diversity generated by terpene synthases (TPSs) and downstream tailoring enzymes. In Medicago truncatula, the TPS-b subfamily comprises five putative synthases, two of which are embedded within a previously uncharacterized genomic locus containing a cytochrome P450 (CYP) and a BAHD-type acyltransferase. Here, we present a comprehensive biochemical analysis of the M. truncatula TPS-b subfamily and define a cineole-centered monoterpene biosynthetic module. Heterologous expression and in vitro assays with multiple prenyl diphosphate substrates revealed three catalytically active TPS-b enzymes with distinct substrate preferences and product profiles. MtTPS4 functions as a dedicated (E)-&#x3b2;-ocimene synthase, whereas MtTPS15 exhibits substrate-dependent bifunctionality, producing (E)-&#x3b2;-ocimene from geranyl diphosphate and &#x3b1;-farnesene from farnesyl diphosphate. MtTPS36 generates 1,8-cineole as the predominant product alongside &#x3b1;-terpineol from geranyl and neryl diphosphate. Genome analysis revealed that MtTPS36 is colocalized with a cytochrome P450 belonging to the CYP736 family. Biochemical characterization of this CYP identified a previously undescribed plant cineole hydroxylase that catalyzes oxidation of 1,8-cineole to yield 2&#x3b1;-hydroxy-1,8-cineole (also known as 2-exo-hydroxy-1,8-cineole), establishing a TPS-CYP biosynthetic module. These results define the gene-to-metabolite relationships within the TPS-b subfamily in M. truncatula and expand the known enzymatic biosynthetic capacity underlying oxygenated monoterpene biosynthesis in plants.

1,8-Cineole

Modulation of hexokinase association with mitochondria analyzed with quantitative three-dimensional confocal microscopy.

Hexokinase isozyme I is proposed to be associated with mitochondria in vivo. Moreover, it has been suggested that this association is modulated in coordination with changes in cell metabolic state. To test these hypotheses, we analyzed the subcellular distribution of hexokinase relative to mitochondria in paraformaldehyde-fixed astrocytes using immunocytochemistry and quantitative three-dimensional confocal microscopy. Analysis of the extent of colocalization between hexokinase and mitochondria revealed that approximately 70% of cellular hexokinase is associated with mitochondria under basal metabolic conditions. In contrast to the immunocytochemical studies, between 15 to 40% of cellular hexokinase was found to be associated with mitochondria after fractionation of astrocyte cultures depending on the exact fractionation conditions. The discrepancy between fractionation studies and those based on imaging of distributions in fixed cells indicates the usefulness of using techniques that can evaluate the distributions of "cytosolic" enzymes in cells whose subcellular ultrastructure is not severely disrupted. To determine if hexokinase distribution is modulated in concert with changes in cell metabolism, the localization of hexokinase with mitochondria was evaluated after inhibition of glucose metabolism with 2-deoxyglucose. After incubation with 2-deoxyglucose there was an approximate 35% decrease in the amount of hexokinase associated with mitochondria. These findings support the hypothesis that hexokinase is bound to mitochondria in rat brain astrocytes in vivo, and that this association is sensitive to cell metabolic state.

Animals

Adenovirus L1 52- and 55-kilodalton proteins are present within assembling virions and colocalize with nuclear structures distinct from replication centers.

Analysis of a temperature-sensitive mutant, Ad5ts369, had indicated that the adenovirus L1 52- and 55-kDa proteins (52/55-kDa proteins) are required for the assembly of infectious virions. By using monoclonal antibodies directed against bacterially produced L1 52-kDa protein, the L1 52/55-kDa proteins were found to be differentially phosphorylated forms of a single 48-kDa polypeptide. Both phosphoforms were shown to be present within all suspected virus assembly intermediates (empty capsids, 50 to 100 molecules; young virions, 1 to 2 molecules) but not within mature virions. The mobilities of these proteins in polyacrylamide gels were affected by reducing agents, indicating that the 52/55-kDa proteins may exist as homodimers within the cell and within assembling particles. Immunofluorescence analysis revealed that the 52/55-kDa proteins localize to regions within the infected nucleus that are distinct from viral DNA replication centers, indicating that replication and assembly of viral components likely occur in separate nuclear compartments. Immunoelectron microscopic studies determined that the 52/55-kDa proteins are found in close association with structures that appear to contain assembling virions. These results are consistent with an active but transient role for the L1 products in assembly of the adenovirus particle, perhaps as scaffolding proteins.

Adenoviridae

Monoclonal antibody MT2 identifies an extracellular matrix glycoprotein that is co-localized with tenascin during adult newt limb regeneration.

Using immunohistochemical techniques and mAb MT2, we describe here a novel extracellular matrix (ECM) molecule that is developmentally regulated during limb regeneration in adult newts. The MT2 antigen appears during preblastema stages, is most abundant during blastema stages, and persists, near undifferentiated cells, until digit stages. The MT2 antigen is located in an acellular layer under the wound epithelium and throughout the ECM of the undifferentiated mesenchyme as a thick, cord-like component. In unamputated limbs mAb MT2 reactivity is restricted to tendons, myotendinous junctions, periosteum and to a layer of material beneath the epidermis. In both unamputated limbs and regenerating limbs, the reactivity to mAb MT2 colocalizes closely with urodele tenascin. Immunoblot analysis of blastema extracts showed that the unreduced form of the MT2 antigen is a large, polydispersed protein of approximately the same size as tenascin. However, based upon (a) molecular weights of reduced subunits, (b) competition experiments on tissue sections, and (c) analysis of molecules immunoprecipitated by mAb MT2, we conclude that the MT2 substance is unrelated biochemically to tenascin. The results from immunoblots, enzyme digestions and DEAE-Sephacell binding studies suggest that the unreduced MT2 antigen is a large protein composed of subunits which are connected by disulfide bonds. Reduction of the MT2 antigen results in three components recognized by mAb MT2. The largest of these reduced components is a chondroitin sulfate-like glycoprotein with a molecular weight (Mr) of 310-325 x 10(3). A second component (Mr, 285-300 x 10(3)) is the core protein of the 310-325 x 10(3) glycoprotein.(ABSTRACT TRUNCATED AT 250 WORDS)

Animals