Glycan reachability analysis: A bottleneck-aware framework for inferring tissue-specific glycan biosynthetic potential from transcriptomics

Glycan biosynthesis requires the coordinated expression of glycosyltransferases, modifying enzymes, and nucleotide-sugar synthesis and transport machinery. Existing computational tools predict glycan structures from gene expression using binary thresholds, losing quantitative information about relative biosynthetic capacity across tissues. Here we present glycan biosynthetic reachability analysis, which integrates expression-based Z-scores across curated pathway steps using AND/OR logic and minimum aggregation to produce continuous, tissue-comparable scores and an explicit expression-limiting step. Applied to 17,382 GTEx v8 RNA-seq samples across 54 human tissue types, reachability resolved quantitative differences hidden by presence/absence calls; for example, pancreas was 96% binary-positive for the sLeX pathway but had low median reachability ( Z = − 1.86 ). Bottleneck stability was evaluated across all 19 multi-step metrics. In independent HEK293 knockout glycomics, the minimum score contained information beyond a within-knockout permutation null but did not outperform naive mean aggregation or binary topology. Nested leave-one-knockout-out selection favored a relaxed low quantile ( q = 0.2) rather than validating the strict minimum. Within-GTEx associations between reachability and signaling-response transcripts are reported only as transcriptomic coherence because predictors and readouts share the same RNA-seq source. The mouse tissue-glycome comparison remained null. Reachability is therefore a hypothesis-generating rank of transcriptomic potential, not a measure of enzymatic activity or glycan abundance.