Daan van Beek, Aishwarya Iyer, Friederike Ehrhart, Chris T Evelo, Theo M de Kok, Ilja C W Arts, Michiel E Adriaens, Martina Kutmon
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.
MOTIVATION: Huntington's 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.
RESULTS: 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.
AVAILABILITY: All analysis code, Docker containers, and conda environments are available at https://github.com/macsbio/HD-ASE-NBS.