{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,8,24]],"date-time":"2026-08-24T11:25:19Z","timestamp":1787570719082,"version":"build-2736575974"},"reference-count":42,"publisher":"Oxford University Press (OUP)","issue":"Supplement_2","license":[{"start":{"date-parts":[[2026,8,24]],"date-time":"2026-08-24T00:00:00Z","timestamp":1787529600000},"content-version":"vor","delay-in-days":23,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"name":"European Union\u2019s Horizon 2020 research and innovation programme"},{"name":"EJP RD COFUND-EJP N\u00b0","award":["825575"],"award-info":[{"award-number":["825575"]}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2026,8,1]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:sec>\n                    <jats:title>Motivation<\/jats:title>\n                    <jats:p>Huntington\u2019s disease (HD) exhibits substantial variability in age of onset and disease progression that is not fully explained by CAG repeat length alone. Part of this residual variation is heritable, implicating additional genetic mechanisms. cis-regulatory variation, genetic variants that alter transcription and splicing of nearby genes, represents one such mechanism that can be quantified through allele-specific expression (ASE) analysis. However, methods for integrating ASE profiles into patient stratification frameworks remain underdeveloped, particularly for rare diseases with small cohorts and sparse data.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Results<\/jats:title>\n                    <jats:p>We adapt a network-based stratification algorithm, originally developed for somatic tumour mutations, to ASE data. By propagating gene-level ASE imbalance profiles through a protein-protein interaction network, we stratified 20 HD patients into three distinct biological patient subgroups. Differential gene expression analysis highlights neuroinflammatory pathways, including microglial activation, immune cell activation, and cytokine regulation, as key sources of inter-patient heterogeneity, while differential ASE analysis implicates proteasomal and ubiquitin-dependent protein catabolic processes, immune activation, and central nervous system development. Intersection of differentially imbalanced and expressed genes identified FAM181B as a candidate gene with potential eQTL-mediated regulation, supported by independent cis-eQTL evidence for rs3780 in the caudate and putamen, the primary HD-affected striatal regions. FAM181B encodes a nuclear protein expressed in neural tissues acting as an interactor of the Hippo pathway TEAD transcription factors, implicating transcriptional regulatory variation as a potential contributor to molecular heterogeneity between patient subgroups. Differences in cortical and striatal neuropathological scores between clusters, even when adjusted for CAG repeat length, provide clinical support for the biological relevance of the identified subgroups.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Availability<\/jats:title>\n                    <jats:p>All analysis code, Docker containers, and conda environments are available at https:\/\/github.com\/macsbio\/HD-ASE-NBS.<\/jats:p>\n                  <\/jats:sec>","DOI":"10.1093\/bioinformatics\/btag592","type":"journal-article","created":{"date-parts":[[2026,8,6]],"date-time":"2026-08-06T11:59:04Z","timestamp":1786017544000},"source":"Crossref","is-referenced-by-count":0,"title":["Network-based stratification of allele-specific expression reveals patient subgroups in Huntington\u2019s disease"],"prefix":"10.1093","volume":"42","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-1273-9756","authenticated-orcid":false,"given":"Daan","family":"van Beek","sequence":"first","affiliation":[{"name":"Maastricht Centre for Systems Biology and Bioinformatics (MaCSBio), Maastricht University , Maastricht,","place":["the Netherlands"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-2878-4049","authenticated-orcid":false,"given":"Aishwarya","family":"Iyer","sequence":"additional","affiliation":[{"name":"Maastricht Centre for Systems Biology and Bioinformatics (MaCSBio), Maastricht University , Maastricht,","place":["the Netherlands"]},{"name":"Department of Translational Genomics, GROW\/NUTRIM\/MHeNs, Maastricht University , Maastricht,","place":["the Netherlands"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-7770-620X","authenticated-orcid":false,"given":"Friederike","family":"Ehrhart","sequence":"additional","affiliation":[{"name":"Department of Translational Genomics, GROW\/NUTRIM\/MHeNs, Maastricht University , Maastricht,","place":["the Netherlands"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-5301-3142","authenticated-orcid":false,"given":"Chris T","family":"Evelo","sequence":"additional","affiliation":[{"name":"Department of Translational Genomics, GROW\/NUTRIM\/MHeNs, Maastricht University , Maastricht,","place":["the Netherlands"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-2102-3626","authenticated-orcid":false,"given":"Theo M","family":"de Kok","sequence":"additional","affiliation":[{"name":"Department of Translational Genomics, GROW\/NUTRIM\/MHeNs, Maastricht University , Maastricht,","place":["the Netherlands"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-6462-6692","authenticated-orcid":false,"given":"Ilja C W","family":"Arts","sequence":"additional","affiliation":[{"name":"Maastricht Centre for Systems Biology and Bioinformatics (MaCSBio), Maastricht University , Maastricht,","place":["the Netherlands"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-4472-7119","authenticated-orcid":false,"given":"Michiel E","family":"Adriaens","sequence":"additional","affiliation":[{"name":"Maastricht Centre for Systems Biology and Bioinformatics (MaCSBio), Maastricht University , Maastricht,","place":["the Netherlands"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-7699-8191","authenticated-orcid":false,"given":"Martina","family":"Kutmon","sequence":"additional","affiliation":[{"name":"Maastricht Centre for Systems Biology and Bioinformatics (MaCSBio), Maastricht University , Maastricht,","place":["the Netherlands"]}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2026,8,21]]},"reference":[{"key":"2026082407022146300_btag592-B1","author":"Asiimwe","year":"2024"},{"key":"2026082407022146300_btag592-B2","doi-asserted-by":"crossref","first-page":"68","DOI":"10.1038\/nature15393","article-title":"A global reference for human genetic variation","volume":"526","author":"Auton","year":"2015","journal-title":"Nature"},{"key":"2026082407022146300_btag592-B3","doi-asserted-by":"crossref","first-page":"509","DOI":"10.1002\/pro.3775","article-title":"Identification of FAM181A and FAM181B as new interactors with the TEAD transcription factors","volume":"29","author":"Bokhovchuk","year":"2020","journal-title":"Protein Sci"},{"key":"2026082407022146300_btag592-B4","author":"Carlson","year":"2025"},{"key":"2026082407022146300_btag592-B5","doi-asserted-by":"crossref","first-page":"195","DOI":"10.1186\/s13059-015-0762-6","article-title":"Tools and best practices for data processing in allelic expression analysis","volume":"16","author":"Castel","year":"2015","journal-title":"Genome Biol"},{"key":"2026082407022146300_btag592-B6","doi-asserted-by":"crossref","first-page":"485","DOI":"10.1007\/s00439-016-1654-x","article-title":"Allele-specific transcription factor binding to common and rare variants associated with disease and gene expression","volume":"135","author":"Cavalli","year":"2016","journal-title":"Hum Genet"},{"key":"2026082407022146300_btag592-B7","doi-asserted-by":"crossref","first-page":"i884","DOI":"10.1093\/bioinformatics\/bty560","article-title":"fastp: an ultra-fast all-in-one FASTQ preprocessor","volume":"34","author":"Chen","year":"2018","journal-title":"Bioinformatics"},{"key":"2026082407022146300_btag592-B8","doi-asserted-by":"crossref","first-page":"364","DOI":"10.1016\/j.it.2015.04.007","article-title":"The choreography of neuroinflammation in Huntington\u2019s disease","volume":"36","author":"Crotti","year":"2015","journal-title":"Trends Immunol"},{"key":"2026082407022146300_btag592-B9","doi-asserted-by":"crossref","first-page":"279","DOI":"10.1002\/ajmg.a.20190","article-title":"Interaction of normal and expanded CAG repeat sizes influences age at onset of Huntington disease","volume":"119A","author":"Djouss\u00e9","year":"2003","journal-title":"Am J Med Genet A"},{"key":"2026082407022146300_btag592-B10","doi-asserted-by":"crossref","first-page":"15","DOI":"10.1093\/bioinformatics\/bts635","article-title":"STAR: ultrafast universal RNA-seq aligner","volume":"29","author":"Dobin","year":"2013","journal-title":"Bioinformatics"},{"key":"2026082407022146300_btag592-B11","doi-asserted-by":"crossref","first-page":"1184","DOI":"10.1038\/nprot.2009.97","article-title":"Mapping identifiers for the integration of genomic datasets with the R\/Bioconductor package biomaRt","volume":"4","author":"Durinck","year":"2009","journal-title":"Nat Protoc"},{"key":"2026082407022146300_btag592-B12","doi-asserted-by":"crossref","first-page":"D948","DOI":"10.1093\/nar\/gkae1071","article-title":"Ensembl 2025","volume":"53","author":"Dyer","year":"2025","journal-title":"Nucleic Acids Res"},{"key":"2026082407022146300_btag592-B13","doi-asserted-by":"crossref","first-page":"18696","DOI":"10.1038\/s41598-019-55177-9","article-title":"Transcriptional correlates of the pathological phenotype in a Huntington\u2019s disease mouse model","volume":"9","author":"Gallardo-Orihuela","year":"2019","journal-title":"Sci Rep"},{"key":"2026082407022146300_btag592-B14","doi-asserted-by":"crossref","first-page":"1108","DOI":"10.1038\/nmeth.2651","article-title":"Network-based stratification of tumour mutations","volume":"10","author":"Hofree","year":"2013","journal-title":"Nat Methods"},{"key":"2026082407022146300_btag592-B1556293","doi-asserted-by":"crossref","first-page":"2859","DOI":"10.1093\/bioinformatics\/bty186","article-title":"pyNBS: a Python implementation for network-based stratification of tumor mutations","volume":"34","author":"Huang","year":"2018","journal-title":"Bioinformatics"},{"key":"2026082407022146300_btag592-B15","doi-asserted-by":"crossref","first-page":"D387","DOI":"10.1093\/nar\/gkab1053","article-title":"The Sequence Read Archive: a decade more of explosive growth","volume":"50","author":"Katz","year":"2022","journal-title":"Nucleic Acids Res"},{"key":"2026082407022146300_btag592-B16","doi-asserted-by":"crossref","first-page":"583","DOI":"10.1080\/01621459.1952.10483441","article-title":"Use of ranks in one-criterion variance analysis","volume":"47","author":"Kruskal","year":"1952","journal-title":"J Am Stat Assoc"},{"key":"2026082407022146300_btag592-B17","doi-asserted-by":"crossref","first-page":"e0143563","DOI":"10.1371\/journal.pone.0143563","article-title":"RNA sequence analysis of human Huntington disease brain reveals an extensive increase in inflammatory and developmental gene expression","volume":"10","author":"Labadorf","year":"2015","journal-title":"PLoS One"},{"key":"2026082407022146300_btag592-B18","doi-asserted-by":"crossref","first-page":"e0160295","DOI":"10.1371\/journal.pone.0160295","article-title":"Correction: RNA sequence analysis of human Huntington disease brain reveals an extensive increase in inflammatory and developmental gene expression","volume":"11","author":"Labadorf","year":"2016","journal-title":"PLoS One"},{"key":"2026082407022146300_btag592-B19","doi-asserted-by":"crossref","first-page":"887","DOI":"10.1016\/j.cell.2019.06.036","article-title":"CAG repeat not polyglutamine length determines timing of Huntington\u2019s disease onset","volume":"178","author":"Lee","year":"2019","journal-title":"Cell"},{"key":"2026082407022146300_btag592-B20","doi-asserted-by":"crossref","first-page":"882","DOI":"10.1093\/bioinformatics\/bts034","article-title":"The sva package for removing batch effects and other unwanted variation in high-throughput experiments","volume":"28","author":"Leek","year":"2012","journal-title":"Bioinformatics"},{"key":"2026082407022146300_btag592-B21","doi-asserted-by":"crossref","first-page":"923","DOI":"10.1093\/bioinformatics\/btt656","article-title":"featureCounts: an efficient general purpose program for assigning sequence reads to genomic features","volume":"30","author":"Liao","year":"2014","journal-title":"Bioinformatics"},{"key":"2026082407022146300_btag592-B22","doi-asserted-by":"crossref","first-page":"82","DOI":"10.1093\/bioinformatics\/btaa1099","article-title":"A network-based deep learning methodology for stratification of tumour mutations","volume":"37","author":"Liu","year":"2021","journal-title":"Bioinformatics"},{"key":"2026082407022146300_btag592-B23","doi-asserted-by":"crossref","first-page":"580","DOI":"10.1038\/ng.2653","article-title":"The Genotype-Tissue Expression (GTEx) project","volume":"45","author":"Lonsdale","year":"2013","journal-title":"Nat Genet"},{"key":"2026082407022146300_btag592-B24","doi-asserted-by":"crossref","first-page":"550","DOI":"10.1186\/s13059-014-0550-8","article-title":"Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2","volume":"15","author":"Love","year":"2014","journal-title":"Genome Biol"},{"key":"2026082407022146300_btag592-B25","doi-asserted-by":"crossref","first-page":"4373","DOI":"10.1093\/hmg\/ddq358","article-title":"The Huntington\u2019s disease mutation impairs Huntingtin\u2019s role in the transport of NF-\u03baB from the synapse to the nucleus","volume":"19","author":"Marcora","year":"2010","journal-title":"Hum Mol Genet"},{"key":"2026082407022146300_btag592-B26","doi-asserted-by":"crossref","first-page":"438","DOI":"10.1016\/j.gene.2015.09.035","article-title":"Analysis of the Fam181 gene family during mouse development reveals distinct strain-specific expression patterns, suggesting a role in nervous system development and function","volume":"575","author":"Marks","year":"2016","journal-title":"Gene"},{"key":"2026082407022146300_btag592-B27","doi-asserted-by":"crossref","first-page":"24","DOI":"10.1111\/ene.13413","article-title":"Huntington\u2019s disease: a clinical review","volume":"25","author":"Mccolgan","year":"2018","journal-title":"Eur J Neurol"},{"key":"2026082407022146300_btag592-B28","doi-asserted-by":"crossref","first-page":"1297","DOI":"10.1101\/gr.107524.110","article-title":"The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data","volume":"20","author":"McKenna","year":"2010","journal-title":"Genome Res"},{"key":"2026082407022146300_btag592-B29","doi-asserted-by":"crossref","first-page":"91","DOI":"10.31887\/DCNS.2016.18.1\/pnopoulos","article-title":"Huntington disease: a single-gene degenerative disorder of the striatum","volume":"18","author":"Nopoulos","year":"2016","journal-title":"Dialogues Clin Neurosci"},{"key":"2026082407022146300_btag592-B30","doi-asserted-by":"crossref","first-page":"c3109","DOI":"10.1136\/bmj.c3109","article-title":"Huntington\u2019s disease","volume":"340","author":"Novak","year":"2010","journal-title":"BMJ"},{"key":"2026082407022146300_btag592-B31","doi-asserted-by":"crossref","first-page":"e4150","DOI":"10.1371\/journal.pone.0004150","article-title":"Allele-specific gene expression is widespread across the genome and biological processes","volume":"4","author":"Palacios","year":"2009","journal-title":"PLoS One"},{"key":"2026082407022146300_btag592-B32","doi-asserted-by":"crossref","first-page":"533","DOI":"10.1038\/nrg2815","article-title":"Genome-wide allele-specific analysis: insights into regulatory variation","volume":"11","author":"Pastinen","year":"2010","journal-title":"Nat Rev Genet"},{"key":"2026082407022146300_btag592-B33","volume-title":"R: a language and environment for statistical computing","author":"R Core Team","year":"2025"},{"key":"2026082407022146300_btag592-B34","doi-asserted-by":"crossref","first-page":"111007","DOI":"10.1016\/j.isci.2024.111007","article-title":"Allele-specific DNA demethylation editing leads to stable upregulation of allele-specific gene expression","volume":"27","author":"Rajaram","year":"2024","journal-title":"iScience"},{"key":"2026082407022146300_btag592-B35","doi-asserted-by":"crossref","first-page":"1098","DOI":"10.1097\/WNR.0b013e32832e34ee","article-title":"Distinct neuroinflammatory profile in post-mortem human Huntington\u2019s disease","volume":"20","author":"Silvestroni","year":"2009","journal-title":"Neuroreport"},{"key":"2026082407022146300_btag592-B36","doi-asserted-by":"crossref","first-page":"D638","DOI":"10.1093\/nar\/gkac1000","article-title":"The STRING database in 2023: protein\u2013protein association networks and functional enrichment analyses for any sequenced genome of interest","volume":"51","author":"Szklarczyk","year":"2023","journal-title":"Nucleic Acids Res"},{"key":"2026082407022146300_btag592-B37","doi-asserted-by":"crossref","first-page":"3498","DOI":"10.1073\/pnas.0308679101","article-title":"Venezuelan kindreds reveal that genetic and environmental factors modulate Huntington's disease age of onset","volume":"101","author":"Wexler","year":"2004","journal-title":"Proc Natl Acad Sci U S A"},{"key":"2026082407022146300_btag592-B38","doi-asserted-by":"crossref","first-page":"564","DOI":"10.1038\/s41598-023-27591-7","article-title":"Allele-specific expression analysis for complex genetic phenotypes applied to a unique dilated cardiomyopathy cohort","volume":"13","author":"van Beek","year":"2023","journal-title":"Sci Rep"},{"key":"2026082407022146300_btag592-B39","doi-asserted-by":"crossref","first-page":"1061","DOI":"10.1038\/nmeth.3582","article-title":"WASP: allele-specific software for robust molecular quantitative trait locus discovery","volume":"12","author":"van de Geijn","year":"2015","journal-title":"Nat Methods"},{"key":"2026082407022146300_btag592-B40","doi-asserted-by":"crossref","first-page":"1338","DOI":"10.12688\/f1000research.15931.1","article-title":"FastQ screen: a tool for multi-genome mapping and quality control","volume":"7","author":"Wingett","year":"2018","journal-title":"F1000Res"},{"key":"2026082407022146300_btag592-B41","doi-asserted-by":"crossref","first-page":"284","DOI":"10.1089\/omi.2011.0118","article-title":"clusterProfiler: an R package for comparing biological themes among gene clusters","volume":"16","author":"Yu","year":"2012","journal-title":"OMICS"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/42\/Supplement_2\/btag592\/70726140\/btag592.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/42\/Supplement_2\/btag592\/70726140\/btag592.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,8,24]],"date-time":"2026-08-24T11:02:33Z","timestamp":1787569353000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/doi\/10.1093\/bioinformatics\/btag592\/8767273"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,8,1]]},"references-count":42,"journal-issue":{"issue":"Supplement_2","published-print":{"date-parts":[[2026,8,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btag592","relation":{},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"value":"1367-4803","type":"print"},{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2026,8]]},"published":{"date-parts":[[2026,8,1]]},"article-number":"btag592"}}