{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,1,23]],"date-time":"2026-01-23T10:27:49Z","timestamp":1769164069985,"version":"3.49.0"},"reference-count":30,"publisher":"Public Library of Science (PLoS)","issue":"1","license":[{"start":{"date-parts":[[2026,1,22]],"date-time":"2026-01-22T00:00:00Z","timestamp":1769040000000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"name":"Jiangsu Provincial Research Hospital","award":["Grant No.YJXYY202204 - YSB23"],"award-info":[{"award-number":["Grant No.YJXYY202204 - YSB23"]}]},{"DOI":"10.13039\/100017962","name":"Jiangsu Commission of Health","doi-asserted-by":"publisher","award":["Grant No. X202356"],"award-info":[{"award-number":["Grant No. X202356"]}],"id":[{"id":"10.13039\/100017962","id-type":"DOI","asserted-by":"publisher"}]},{"name":"Shanghai Haitian Pharmaceutical Technology Development Co., Ltd.","award":["Grant No. LCYJ20242030"],"award-info":[{"award-number":["Grant No. LCYJ20242030"]}]},{"name":"Scientific Research Project of Jiangsu Provincial Traditional Chinese Medicine Association","award":["Grant No. CYTF2024049"],"award-info":[{"award-number":["Grant No. CYTF2024049"]}]}],"content-domain":{"domain":["www.plosone.org"],"crossmark-restriction":false},"short-container-title":["PLoS One"],"abstract":"<jats:sec id=\"sec001\">\n                    <jats:title>Objective<\/jats:title>\n                    <jats:p>Hyperthyroidism and hypothyroidism are globally prevalent endocrine disorders, with their pathogenesis involving multifactorial mechanisms including genetics, immunity, and metabolism. Although genome-wide association studies (GWAS) have identified risk genes such as PDE8B, critical gaps remain in the annotation of causal variants in non-coding regions, characterization of tissue-specific regulatory networks, and understanding of ethnic heterogeneity. This study aimed to systematically identify genes associated with hyperthyroidism and hypothyroidism and unravel their underlying molecular mechanisms through multi-omics integration.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec id=\"sec002\">\n                    <jats:title>Methods<\/jats:title>\n                    <jats:p>We included data from the ThyroidOmics Consortium, comprising 1,840 hyperthyroidism cases (49,983 controls) and 3,340 hypothyroidism cases (49,983 controls). Core candidate genes were prioritized using a combination of SMR-HEIDI analysis, cross-tissue transcriptome-wide association study (TWAS), mBAT-combo rare variant analysis, and polygenic priority score (PoPS). GTEx colocalization (coloc) analysis was used to validate tissue-specific colocalization between these candidate genes and disease signals. Phenome-wide association study (PheWAS), KEGG pathway enrichment, and protein-protein interaction (PPI) network analyses were performed to explore gene functions, with potential targeted drugs predicted using the Drug Signatures Database (DSigDB).<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec id=\"sec003\">\n                    <jats:title>Results<\/jats:title>\n                    <jats:p>Cross-validation by four methods identified FAM227B, PDE8B, and PDE10A as key genes for hyperthyroidism, and PDE8B as the critical gene for hypothyroidism. GTEx coloc analysis (with PP4\u2009&gt;\u20090.8 as the threshold) confirmed significant colocalization: FAM227B with hyperthyroidism signals in the adrenal gland, lung, and minor salivary gland; PDE8B with both hyperthyroidism and hypothyroidism signals in thyroid tissue; and PDE10A with hyperthyroidism signals in thyroid tissue. As a core member of the phosphodiesterase family, PDE8B exhibited bidirectional regulatory characteristics in thyroid hormone synthesis via the cAMP signaling pathway and nucleotide metabolism network: its inhibition promoted hormone synthesis in hypothyroidism, while its interaction with PDE10A suppressed overactive cAMP signaling in hyperthyroidism. PheWAS linked FAM227B to cardiovascular diseases and PDE10A to neurological pathways. KEGG enrichment analysis highlighted the \u201cmorphine addiction\u201d pathway (p\u2009=\u20096.12\u2009\u00d7\u200910\u2009\u207b\u2009\u2075), suggesting potential neuroendocrine crosstalk. Notably, potential drugs targeting FAM227B, PDE8B, and PDE10A were identified.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec id=\"sec004\">\n                    <jats:title>Conclusion<\/jats:title>\n                    <jats:p>Through multi-omics integration, this study identifies PDE8B as a central gene associated with thyroid dysfunction, characterized by tissue-specific colocalization, and elucidates its critical roles in signaling pathways, comorbidity associations, and drug targeting. These findings provide insights into the bidirectional regulatory mechanisms of hyperthyroidism and hypothyroidism and a theoretical basis for developing phosphodiesterase family-based precision therapies. It should be noted that all samples in this study are of European ancestry, and the generalizability of the results in other ethnic groups remains to be verified.<\/jats:p>\n                  <\/jats:sec>","DOI":"10.1371\/journal.pone.0338805","type":"journal-article","created":{"date-parts":[[2026,1,22]],"date-time":"2026-01-22T18:33:05Z","timestamp":1769106785000},"page":"e0338805","update-policy":"https:\/\/doi.org\/10.1371\/journal.pone.corrections_policy","source":"Crossref","is-referenced-by-count":0,"title":["Integrative multi-omics analyses identify key genes and elucidate bidirectional regulatory mechanisms in thyroid dysfunction"],"prefix":"10.1371","volume":"21","author":[{"given":"Yilang","family":"Hu","sequence":"first","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Hong","family":"Hu","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Lijun","family":"Xu","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Xuzheng","family":"Tang","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Chongyu","family":"Wang","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Zhengyi","family":"Lian","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Chenyu","family":"Xie","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Zhiwan","family":"Xie","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0009-0008-8235-9346","authenticated-orcid":true,"given":"Qingqing","family":"Wang","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"340","published-online":{"date-parts":[[2026,1,22]]},"reference":[{"key":"pone.0338805.ref001","first-page":"5736394","article-title":"Transcriptomic Analysis Reveals Genetic Cross-Talk between Periodontitis and Hypothyroidism","volume":"2022","author":"B Yan","year":"2022","journal-title":"Dis Markers"},{"issue":"1","key":"pone.0338805.ref002","doi-asserted-by":"crossref","first-page":"108","DOI":"10.1186\/s12874-023-01898-5","article-title":"Utilization of five data mining algorithms combined with simplified preprocessing to establish reference intervals of thyroid-related hormones for non-elderly adults","volume":"23","author":"J Zhong","year":"2023","journal-title":"BMC Med Res Methodol"},{"issue":"10101","key":"pone.0338805.ref003","doi-asserted-by":"crossref","first-page":"1550","DOI":"10.1016\/S0140-6736(17)30703-1","article-title":"Hypothyroidism","volume":"390","author":"L Chaker","year":"2017","journal-title":"Lancet"},{"issue":"19","key":"pone.0338805.ref004","doi-asserted-by":"crossref","first-page":"10628","DOI":"10.3390\/ijms251910628","article-title":"Metabolic Syndrome, Thyroid Dysfunction, and Cardiovascular Risk: The Triptych of Evil","volume":"25","author":"A Pingitore","year":"2024","journal-title":"Int J Mol Sci"},{"issue":"1","key":"pone.0338805.ref005","doi-asserted-by":"crossref","first-page":"4455","DOI":"10.1038\/s41467-018-06356-1","article-title":"Genome-wide analyses identify a role for SLC17A4 and AADAT in thyroid hormone regulation","volume":"9","author":"A Teumer","year":"2018","journal-title":"Nat Commun"},{"key":"pone.0338805.ref006","doi-asserted-by":"crossref","first-page":"17","DOI":"10.1016\/j.molimm.2024.02.002","article-title":"Shared genetic architecture between hypothyroidism and rheumatoid arthritis: A large-scale cross-trait analysis","volume":"168","author":"R Liu","year":"2024","journal-title":"Mol Immunol"},{"issue":"18","key":"pone.0338805.ref007","doi-asserted-by":"crossref","first-page":"9957","DOI":"10.3390\/ijms25189957","article-title":"Myxedema in Both Hyperthyroidism and Hypothyroidism: A Hormetic Response?","volume":"25","author":"S Sciacchitano","year":"2024","journal-title":"Int J Mol Sci"},{"issue":"6","key":"pone.0338805.ref008","doi-asserted-by":"crossref","first-page":"580","DOI":"10.1038\/ng.2653","article-title":"The Genotype-Tissue Expression (GTEx) project","volume":"45","author":"GTEx Consortium","year":"2013","journal-title":"Nat Genet"},{"issue":"1","key":"pone.0338805.ref009","doi-asserted-by":"crossref","first-page":"605","DOI":"10.1186\/s12864-023-09710-w","article-title":"Summary-data based Mendelian randomization identifies gene expression regulatory polymorphisms associated with bovine paratuberculosis by modulation of the nuclear factor Kappa \u03b2 (NF-\u03ba\u00df)-mediated inflammatory response","volume":"24","author":"G Badia-Bringu\u00e9","year":"2023","journal-title":"BMC Genomics"},{"key":"pone.0338805.ref010","doi-asserted-by":"crossref","first-page":"745757","DOI":"10.3389\/fcvm.2021.745757","article-title":"Mendelian Randomization Integrating GWAS, eQTL, and mQTL Data Identified Genes Pleiotropically Associated With Atrial Fibrillation","volume":"8","author":"Y Liu","year":"2021","journal-title":"Front Cardiovasc Med"},{"issue":"5","key":"pone.0338805.ref011","doi-asserted-by":"crossref","first-page":"481","DOI":"10.1038\/ng.3538","article-title":"Integration of summary data from GWAS and eQTL studies predicts complex trait gene targets","volume":"48","author":"Z Zhu","year":"2016","journal-title":"Nat Genet"},{"issue":"6","key":"pone.0338805.ref012","doi-asserted-by":"crossref","first-page":"1035","DOI":"10.1016\/j.ajhg.2024.04.016","article-title":"An integrative framework to prioritize genes in more than 500 loci associated with body mass index","volume":"111","author":"D Hemerich","year":"2024","journal-title":"Am J Hum Genet"},{"issue":"3","key":"pone.0338805.ref013","doi-asserted-by":"crossref","first-page":"245","DOI":"10.1038\/ng.3506","article-title":"Integrative approaches for large-scale transcriptome-wide association studies","volume":"48","author":"A Gusev","year":"2016","journal-title":"Nat Genet"},{"issue":"1","key":"pone.0338805.ref014","doi-asserted-by":"crossref","first-page":"94","DOI":"10.1186\/s10194-024-01802-6","article-title":"A cross-tissue transcriptome-wide association study reveals novel susceptibility genes for migraine","volume":"25","author":"J Gui","year":"2024","journal-title":"J Headache Pain"},{"issue":"1","key":"pone.0338805.ref015","doi-asserted-by":"crossref","first-page":"111","DOI":"10.1186\/s10194-023-01649-3","article-title":"Identifying causal genes for migraine by integrating the proteome and transcriptome","volume":"24","author":"S-J Li","year":"2023","journal-title":"J Headache Pain"},{"issue":"3","key":"pone.0338805.ref016","doi-asserted-by":"crossref","first-page":"410","DOI":"10.1016\/j.ajhg.2019.01.002","article-title":"ACAT: A Fast and Powerful p Value Combination Method for Rare-Variant Analysis in Sequencing Studies","volume":"104","author":"Y Liu","year":"2019","journal-title":"Am J Hum Genet"},{"issue":"1","key":"pone.0338805.ref017","doi-asserted-by":"crossref","first-page":"30","DOI":"10.1016\/j.ajhg.2022.12.006","article-title":"mBAT-combo: A more powerful test to detect gene-trait associations from GWAS data","volume":"110","author":"A Li","year":"2023","journal-title":"Am J Hum Genet"},{"issue":"1","key":"pone.0338805.ref018","doi-asserted-by":"crossref","first-page":"130","DOI":"10.1186\/s40478-024-01841-9","article-title":"Macular structural integrity estimates are associated with Parkinson\u2019s disease genetic risk","volume":"12","author":"S Diaz-Torres","year":"2024","journal-title":"Acta Neuropathol Commun"},{"issue":"5","key":"pone.0338805.ref019","doi-asserted-by":"crossref","DOI":"10.1371\/journal.pgen.1004383","article-title":"Bayesian test for colocalisation between pairs of genetic association studies using summary statistics","volume":"10","author":"C Giambartolomei","year":"2014","journal-title":"PLoS Genet"},{"issue":"17","key":"pone.0338805.ref020","doi-asserted-by":"crossref","first-page":"1666","DOI":"10.1093\/hmg\/ddab119","article-title":"A cross-tissue transcriptome-wide association study identifies novel susceptibility genes for lung cancer in Chinese populations","volume":"30","author":"M Zhu","year":"2021","journal-title":"Hum Mol Genet"},{"issue":"2","key":"pone.0338805.ref021","doi-asserted-by":"crossref","first-page":"732","DOI":"10.1007\/s12020-024-03933-x","article-title":"Prioritization of drug targets for thyroid cancer: a multi-omics Mendelian randomization study","volume":"86","author":"H Sun","year":"2024","journal-title":"Endocrine"},{"issue":"7877","key":"pone.0338805.ref022","doi-asserted-by":"crossref","first-page":"527","DOI":"10.1038\/s41586-021-03855-y","article-title":"Rare variant contribution to human disease in 281,104 UK Biobank exomes","volume":"597","author":"Q Wang","year":"2021","journal-title":"Nature"},{"key":"pone.0338805.ref023","article-title":"Rejuvenation increases leaf biomass and flavonoid accumulation in Ginkgo biloba","volume":"9","author":"Z Lu","year":"2022","journal-title":"Hortic Res"},{"key":"pone.0338805.ref024","doi-asserted-by":"crossref","DOI":"10.1093\/nar\/gkac1000","article-title":"The STRING database in 2023: protein-protein association networks and functional enrichment analyses for any sequenced genome of interest","volume":"51","author":"D Szklarczyk","year":"2023","journal-title":"Nucleic Acids Res"},{"issue":"1","key":"pone.0338805.ref025","doi-asserted-by":"crossref","first-page":"185","DOI":"10.1186\/s13059-019-1758-4","article-title":"Cytoscape Automation: empowering workflow-based network analysis","volume":"20","author":"D Otasek","year":"2019","journal-title":"Genome Biol"},{"key":"pone.0338805.ref026","doi-asserted-by":"crossref","DOI":"10.1093\/nar\/gkq537","article-title":"The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function","volume":"38","author":"D Warde-Farley","year":"2010","journal-title":"Nucleic Acids Res"},{"issue":"18","key":"pone.0338805.ref027","doi-asserted-by":"crossref","first-page":"3069","DOI":"10.1093\/bioinformatics\/btv313","article-title":"DSigDB: drug signatures database for gene set analysis","volume":"31","author":"M Yoo","year":"2015","journal-title":"Bioinformatics"},{"issue":"9","key":"pone.0338805.ref028","doi-asserted-by":"crossref","first-page":"104992","DOI":"10.1016\/j.isci.2022.104992","article-title":"Genetic association and Mendelian randomization for hypothyroidism highlight immune molecular mechanisms","volume":"25","author":"S Mathieu","year":"2022","journal-title":"iScience"},{"issue":"1","key":"pone.0338805.ref029","doi-asserted-by":"crossref","first-page":"5071","DOI":"10.1038\/s41467-021-25392-y","article-title":"Multi-omics integration analysis identifies novel genes for alcoholism with potential overlap with neurodegenerative diseases","volume":"12","author":"M Kapoor","year":"2021","journal-title":"Nat Commun"},{"issue":"1","key":"pone.0338805.ref030","doi-asserted-by":"crossref","first-page":"136","DOI":"10.1186\/s12864-017-3518-8","article-title":"Melanoma genome evolution across species","volume":"18","author":"ER Kansler","year":"2017","journal-title":"BMC Genomics"}],"container-title":["PLOS One"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pone.0338805","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,1,22]],"date-time":"2026-01-22T18:33:11Z","timestamp":1769106791000},"score":1,"resource":{"primary":{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pone.0338805"}},"subtitle":[],"editor":[{"given":"Vijayalakshmi","family":"Kakulapati","sequence":"first","affiliation":[],"role":[{"role":"editor","vocabulary":"crossref"}]}],"short-title":[],"issued":{"date-parts":[[2026,1,22]]},"references-count":30,"journal-issue":{"issue":"1","published-online":{"date-parts":[[2026,1,22]]}},"URL":"https:\/\/doi.org\/10.1371\/journal.pone.0338805","relation":{},"ISSN":["1932-6203"],"issn-type":[{"value":"1932-6203","type":"electronic"}],"subject":[],"published":{"date-parts":[[2026,1,22]]}}}