Language-encoded network topology enables large language models to reason about complex networks
Summary
BioGlyph is a method that translates network topology into interpretable language roles to enhance large language models' ability to reason about complex networks, demonstrating significant performance improvements across various domains.
View Cached Full Text
Cached at: 09/04/26, 06:22 AM
# Language-encoded network topology enables large language models to reason about complex networks
Source: [https://arxiv.org/html/2609.03229](https://arxiv.org/html/2609.03229)
Ucchwas Talukder UtshaSakib MostafaAffiliation:Department of Radiation Oncology, Stanford University, Stanford, California, USAJames ZouAffiliation:Department of Electrical Engineering, Stanford University, Stanford, California, USAAffiliation:Department of Biomedical Data Science, Stanford University School of Medicine, Stanford, California, USAAffiliation:Department of Computer Science, Stanford University, Stanford, California, USAMd Tauhidul IslamAffiliation:Department of Radiation Oncology, Stanford University, Stanford, California, USAAffiliation:correspondence: tauhid@stanford\.edu
###### Abstract
Networks describe diverse systems in biology and beyond, from protein interactions and social relationships to power grids and citation records\. Reasoning about such systems requires understanding their structure: which elements are central, which connections bridge otherwise separate communities, how the network is organized, and how its structure changes when elements are removed\. Although large language models \(LLMs\) excel at reasoning over natural language, they often struggle with such structural questions when networks are represented as edge lists, sentences or tables of numerical measurements, because the structural meaning of these representations must be inferred\. Here we introduce BioGlyph, a method that compiles network topology into an interpretable and transferable language of structural roles\. BioGlyph combines graph\-partitioning algorithms and structural measurements to identify roles such as hubs, community cores and cross\-community connectors, and uses fixed rules to translate their algorithmic signatures into a universal vocabulary\. The resulting representation describes each network element through its structural role, supporting evidence and semantic consequences, while leaving both the original network and the LLM unchanged\. Across twenty networks spanning five different domains including biological, social and information systems, BioGlyph substantially improves the ability of open LLMs to answer structural reasoning questions, outperforming edge\-based, numerical and learned representations by up to 26 percentage points in system accuracy\. Ablation experiments show that the improvement comes from explicitly encoding structural roles in semantically interpretable terms by Bioglyph\. The performance improvement becomes even more prominent in dense, community\-structured networks and diminishes in sparse networks whose topology can be more readily inferred from direct textual representations\. Applied to a budding\-yeast protein\-interaction network, BioGlyph also exposes biologically meaningful organization: cross\-community connectors are enriched for essential genes, whereas peripheral proteins are depleted\. BioGlyph thus provides a common, interpretable representation through which both language models and scientists can reason about complex network structure\.
## Introduction
Figure 1:BioGlyph translates network structure into interpretable roles\.a, An illustrative network of 31 nodes organized into five communities\. BioGlyph identifies complementary structural roles, including hubs, cut nodes, bridges, community cores, cross\-community connectors and bottleneck links\. For example, node A is a highly connected hub, whereas removal of node B separates the network into three components\. The dashed outline indicates a fragile community connected to the rest of the network through a single link\.b, BioGlyph applies classical network\-analysis algorithms and fixed rules to identify eleven structural roles\. The rules use either exact perturbation tests or thresholds derived from the network structure\. The procedure is deterministic for a given random seed and uses network topology alone, without node labels or biological annotations\.c, BioGlyph converts the resulting structural information into a human\- and language\-model\-readable description\. Each entry specifies the structural role, the evidence supporting it and, where relevant, the consequence of removing the node or link\. The resulting description provides an interpretable summary of network organization that can be used by a frozen language model to answer questions about connectivity, importance and perturbation\. Answers can then be checked against exact graph algorithms\.Networks provide a natural way to describe how parts of a system are connected\. Protein interactions inside a cell, gene\-regulatory circuits, and connections in the brain can all be represented as graphs\. The same representation is used for power grids, transport systems, citation records, the World Wide Web, and social networks\. In a graph, nodes represent the parts of a system and edges represent the relationships between them\. Decades of network science have shown that the behavior of a system depends not only on the individual nodes, but also on the pattern of connections among them\[[45](https://arxiv.org/html/2609.03229#bib.bib1),[4](https://arxiv.org/html/2609.03229#bib.bib2)\]\. Several structural patterns appear across very different networks\. Hubs have unusually many connections\. Cut nodes, or articulation points, separate a connected network when removed\. Bridges provide the only link between otherwise disconnected regions\. Many networks also form dense communities connected by a small number of nodes that link one community to another\[[16](https://arxiv.org/html/2609.03229#bib.bib3),[15](https://arxiv.org/html/2609.03229#bib.bib4)\]\. Such structural features can have direct scientific meaning\. Network position has been linked to gene essentiality\[[25](https://arxiv.org/html/2609.03229#bib.bib5),[51](https://arxiv.org/html/2609.03229#bib.bib6)\], infrastructure vulnerability\[[2](https://arxiv.org/html/2609.03229#bib.bib7)\], and the flow of information through social systems\. Understanding a network therefore requires more than knowing which nodes are connected\. The important question is often what those connections imply for the system, especially when a node or edge is removed\.
Large language models are strong general\-purpose systems for reasoning over text\. Models trained on large text collections can learn new tasks from examples in a prompt\[[7](https://arxiv.org/html/2609.03229#bib.bib8)\]and can solve multi\-step problems through explicit reasoning\[[46](https://arxiv.org/html/2609.03229#bib.bib9)\]\. Network reasoning presents a different challenge\. When a graph is written as an edge list or as sentences describing which nodes are connected, a language model must recover the relevant structure from the text before answering the question\. The model must determine, for example, whether a node connects two regions, whether removing an edge separates the graph, or which candidate is most important for connectivity\. Language models often make errors on such problems even when standard graph algorithms can compute the answers exactly\[[44](https://arxiv.org/html/2609.03229#bib.bib10),[13](https://arxiv.org/html/2609.03229#bib.bib11)\]\.
Perturbation\-based questions make this limitation particularly clear\. Can two nodes still reach each other after a a third node is removed? Does deleting a single edge split the network? Which node causes the greatest disruption when removed? Related questions arise when a router fails in a communication system or a connection is disrupted in the brain\. They are especially important in biomedicine, where many questions can be framed as network perturbations\. What happens when a gene is upregulated or downregulated? Which protein is most important for maintaining a pathway? Which component connects two biological processes? Which proteins are most vulnerable to disruption? Answering these questions requires understanding how each component is positioned within the wider biological system\. A list of interactions alone does not explain why a protein is important or predict what will happen when it is perturbed\. In each setting, the effect of a perturbation depends on the structure of the network\.
Current approaches generally provide network information to LLMs either as connections or as learned numerical representations\. Connections can be written as edge lists, structured data or natural\-language descriptions\[[19](https://arxiv.org/html/2609.03229#bib.bib12),[44](https://arxiv.org/html/2609.03229#bib.bib10),[13](https://arxiv.org/html/2609.03229#bib.bib11)\]\. These formats preserve the underlying network, but the model still has to infer important structural properties from the connections\. Other approaches learn numerical representations of nodes or network regions and provide these representations to language models\[[32](https://arxiv.org/html/2609.03229#bib.bib13),[37](https://arxiv.org/html/2609.03229#bib.bib14),[8](https://arxiv.org/html/2609.03229#bib.bib15)\]\. Although such representations can capture complex patterns, they are difficult for researchers to interpret\. A numerical representation does not readily explain why a protein is important, which other proteins it connects, or what might happen when it is removed\. These limitations raise a simple question: can network structure be expressed directly in terms that both language models and researchers can understand?
Here we introduce BioGlyph, a method that translates network structure into a simple language of structural roles \(Fig\.[1](https://arxiv.org/html/2609.03229#Sx1.F1)\)\. BioGlyph first uses established network\-analysis algorithms to identify properties such as hubs, community cores, cut nodes, bridges and proteins that connect different communities\. It then uses fixed rules to describe these properties in words\. Each protein or interaction can therefore be described not only by its structural role, but also by the measurements supporting that role and by what its position means for the surrounding network\. For example, rather than providing only the number of interactions of a protein, BioGlyph can describe it as a hub or as a connector between two communities and explain how the network would change if that protein were removed\. The method uses only the network itself and does not require biological labels or task\-specific training\. The LLM is also kept unchanged\. BioGlyph therefore changes how network information is presented rather than changing the underlying model\. By expressing network properties in words, BioGlyph makes the results of network analysis easier to interpret and incorporate into a broader reasoning process\.
We evaluated BioGlyph across eight networks representing biological, social and information systems\. We asked LLMs questions about network connectivity, importance and the consequences of removing nodes or interactions, while keeping the underlying network and model fixed and changing only how the network was represented\. We compared BioGlyph with edge lists, natural\-language descriptions, numerical network measurements and learned network representations\. BioGlyph consistently improved performance on questions that required understanding network structure, and approached the performance of a graph neural network trained specifically on the task\. The improvement was particularly pronounced in dense networks with many interconnected communities, where directly reading the individual connections is difficult\. In sparse networks, where the connections are already relatively easy to read, the advantage was smaller\. These results indicate that explicitly describing meaningful structural properties can make complex network information easier for LLMs to use\.
We next asked whether the same structural descriptions could reveal biologically meaningful patterns independently of their use for language\-model reasoning\. In the budding\-yeast protein\-interaction network, BioGlyph identified structural roles using network connectivity alone, without using gene\-essentiality annotations\. Proteins that connected different communities were substantially more likely to be essential, whereas peripheral proteins were less likely to be essential\. Cross\-community connectors were also more strongly enriched for essential genes than proteins selected solely by degree at the same set size\. These findings suggest that the structural vocabulary captures biologically meaningful organization rather than simply providing a convenient format for prompting an LLM\.
Together, these results show that network structure can be translated into a common language that is accessible to both computational models and researchers\. BioGlyph preserves information obtained from established network\-analysis methods while making important structural relationships easier to interpret, communicate and use in biomedical reasoning\. This provides a general framework for using LLMs with biological networks and, more broadly, for connecting quantitative network analysis with human\-readable scientific interpretation\.
## Results
### BioGlyph improves reasoning over the yeast interactome and identifies biologically meaningful roles
We first evaluated BioGlyph on STRING\-Yeast, a budding\-yeast protein\-interaction network comprising 3,384 proteins and 43,030 physical associations\[[36](https://arxiv.org/html/2609.03229#bib.bib22)\]\. Independent gene\-essentiality annotations enabled us to test whether roles inferred solely from network topology correspond to biological function\. The benchmark comprised 156 questions, 10\.3% of which referred to a protein absent from the network and therefore required the model to state that the question could not be answered\. For each question, we retrieved the network region surrounding the named proteins and represented it in six forms: BioGlyph; no network information \(R0\); an edge list \(R1\); the same edges expressed as sentences \(R2\); a table of raw graph measurements \(R3\); and structural roles learned using GraphSAGE \(R4\)\. GraphSAGE was the strongest of the four graph encoders evaluated; results for all four are reported in Section[S1](https://arxiv.org/html/2609.03229#S1)\. Each representation was provided to the same two frozen 8\-billion\-parameter language models, Qwen3\-8B and Llama\-3\.1\-8B\. Answers were graded against exact graph\-algorithm outputs rather than by another language model\.
Figure[2](https://arxiv.org/html/2609.03229#Sx2.F2)illustrates how a single BioGlyph description can support sustained reasoning about one network region\. We compiled a region containing 120 proteins and provided the description once, with the first question\. Qwen3\-8B identified ACT1 as the protein whose removal would cause the greatest disruption, citing the stated consequence that its removal leaves eight connected components\. The model subsequently used the same description to answer four further graded questions, recognized ACT1 as a cross\-community connector and retained its original answer when challenged with the alternative suggestion that CBF5 would be more disruptive\. All graded responses agreed with the corresponding exact graph\-algorithm outputs\. ACT1 is independently annotated by SGD as essential, although this annotation was not available to the model\.
Figure 2:A single BioGlyph description supports multi\-turn reasoning about the yeast interactome\.aA STRING\-Yeast region containing 120 proteins and 781 links was compiled once and provided to a frozen Qwen3\-8B\. Nodes are coloured according to the eight connected components that remain after removal of ACT1\. Three of eight conversation turns are shown\. Turn 1 asks which of six proteins would cause the greatest disruption if removed; turn 5 asks whether ACT1 also connects distinct communities; and turn 8 challenges the initial answer by proposing CBF5 instead\. Gene names are shown for readability: node 1027 corresponds to YFL039C \(ACT1\), and node 2164 to YLR175W \(CBF5\)\. The model received only node identifiers and did not receive the SGD annotation\. Responses are taken from stored model outputs and shortened only at sentence boundaries\. Each graded response was evaluated against an exact graph\-algorithm output\.bThe same BioGlyph description was used throughout all eight turns\. Exact graph algorithms supplied the grading reference; no language model graded another model\. Further multi\-turn examples are provided in Section[S2](https://arxiv.org/html/2609.03229#S2)\.Across both models, BioGlyph yielded an accuracy of 75\.6% on the yeast benchmark \(95% confidence interval \(CI\), 72\.3–78\.9%\), compared with 47\.3% without network information \(Fig\.[3](https://arxiv.org/html/2609.03229#Sx2.F3)a\)\. Direct textual or numerical representations of the retrieved graph performed less well: accuracy was 40\.1% for edge lists, 37\.5% for sentences and 36\.1% for raw graph measurements\. GraphSAGE\-derived roles yielded 49\.2% accuracy\. These representations also frequently exceeded the 24,576\-token context limit: 56\.6% of raw\-measurement prompts, 38\.5% of sentence prompts and 30\.6% of edge\-list prompts were over the limit\. By contrast, every BioGlyph description fitted within the context window\.
Figure 3:BioGlyph improves network reasoning on STRING\-Yeast, and its structural roles are associated with gene essentiality\.Results are pooled across Qwen3\-8B and Llama\-3\.1\-8B unless otherwise indicated\.aSystem accuracy for each representation, with 95% bootstrap CIs\. Percentages within selected bars denote the fraction of prompts that exceeded the context window; the dotted line indicates accuracy without network information\.bPairwise comparisons of BioGlyph with each alternative, restricted to questions for which both representations fitted within the context window and produced an answer\. The number of paired answers is shown at right\.cAblation analysis of the BioGlyph description\. Solid bars show system accuracy; hatched bars show accuracy on the common subset of questions readable under every ablation\.dAccuracy across the seven benchmark question families\.eProportion of proteins annotated by SGD as essential within each structural role\. The dashed line indicates the network\-wide essentiality rate of 30\.9%\.fEssentiality among equally sized sets of 85 proteins selected as cross\-community connectors or ranked using a single centrality measure\.To distinguish representational quality from context\-window truncation, we performed pairwise analyses restricted to questions for which both representations fitted within the context window and produced an answer \(Fig\.[3](https://arxiv.org/html/2609.03229#Sx2.F3)b\)\. On the subset shared with the raw\-measurement table, BioGlyph achieved 79\.5% accuracy, compared with 73\.7% for the table\. The corresponding comparisons were 74\.6% versus 57\.8% for sentences, 75\.0% versus 53\.1% for edge lists and 78\.1% versus 47\.1% for GraphSAGE\-derived roles\. Thus, BioGlyph retained an accuracy advantage when both representations were available in full\.
Ablation experiments identified the components of the representation that contributed to this advantage \(Fig\.[3](https://arxiv.org/html/2609.03229#Sx2.F3)c\)\. On the common subset of questions readable under every ablation, role names alone yielded 69\.9% accuracy\. Adding the graph measurements supporting each assignment yielded 68\.8%, whereas additionally stating the structural consequence of each role increased accuracy to 82\.7%\. The semantic content of the role names also contributed\. In an opaque control, each role name was replaced by a meaningless token while the role assignments and all other prompt content were held fixed\. This substitution reduced overall system accuracy from 66\.8% to 51\.1%\.
The gains extended across multiple forms of structural reasoning \(Fig\.[3](https://arxiv.org/html/2609.03229#Sx2.F3)d\)\. For cut\-node questions, BioGlyph achieved 90\.2% accuracy, compared with 48\.9% for sentences and 31\.5% for GraphSAGE\. For questions comparing the degrees of two proteins, BioGlyph and GraphSAGE both achieved 95\.7%, whereas sentences achieved 37\.0%\. BioGlyph also outperformed the alternative representations across the perturbation\-oriented question families\.
We next asked whether the structural roles identified by BioGlyph were associated with an independent biological property\. STRING\-Yeast proteins were matched to SGD gene\-essentiality annotations\[[9](https://arxiv.org/html/2609.03229#bib.bib44)\], which were withheld from both the BioGlyph compiler and the language models\. Of all proteins in the network, 30\.9% were essential\. This proportion increased to 57\.6% among the 85 proteins classified as cross\-community connectors, 46\.9% among community cores and 42\.1% among hubs \(Fig\.[3](https://arxiv.org/html/2609.03229#Sx2.F3)e\)\. By contrast, only 12\.9% of peripheral proteins were essential\. The proportions among cut nodes and boundary nodes were close to the network\-wide rate\.
Finally, we examined whether the cross\-community connector role merely recapitulated a standard centrality ranking\. We selected equally sized sets of 85 proteins using degree, betweenness or PageRank \(Fig\.[3](https://arxiv.org/html/2609.03229#Sx2.F3)f\)\. Essential proteins constituted 56\.5% of the highest\-betweenness set, 48\.2% of the highest\-PageRank set and 24\.7% of the highest\-degree set, compared with 57\.6% of cross\-community connectors\. The connector role therefore identifies a group enriched for essential proteins by combining high betweenness with participation across multiple communities, rather than by applying any one of these centrality measures alone\. Results for the other benchmark networks, additional model families and further controls are reported in Section[S1](https://arxiv.org/html/2609.03229#S1)\.
### BioGlyph remains effective on a dense drug\-interaction network
We next evaluated BioGlyph on ChCh\-Miner, a BioSNAP network in which 1,514 drugs are connected by 48,514 reported drug–drug interactions\[[52](https://arxiv.org/html/2609.03229#bib.bib23)\]\. With a mean degree of approximately 64, ChCh\-Miner is the densest network in our benchmark\. We evaluated 156 questions using the same representations and two frozen 8\-billion\-parameter language models as in the yeast analysis\. Of these questions, 10\.3% referred to a drug absent from the network and therefore required the model to state that they could not be answered\. Across both models, BioGlyph achieved 77\.7% system accuracy \(95% confidence interval \(CI\), 74\.5–80\.9%; Fig\.[4](https://arxiv.org/html/2609.03229#Sx2.F4)a\), compared with 50\.8% without network information\. GraphSAGE\-derived roles yielded 51\.4% accuracy\. Direct representations of the retrieved graph performed less well: accuracy was 41\.7% for edge lists, 38\.3% for sentences and 25\.3% for raw graph measurements\.
The differences in accuracy were accompanied by substantial differences in prompt length\. The median raw\-measurement prompt contained 35,601 tokens, exceeding the 24,576\-token context limit; consequently, 67\.1% of these prompts could not be passed to the model\. The corresponding proportions were 43\.3% for sentence prompts and 30\.1% for edge\-list prompts\. By contrast, the median BioGlyph description contained 6,353 tokens, and every BioGlyph description fitted within the context window \(Fig\.[4](https://arxiv.org/html/2609.03229#Sx2.F4)e\)\. To distinguish representational quality from context\-window truncation, we again restricted pairwise analyses to questions for which both representations fitted within the context window and produced an answer \(Fig\.[4](https://arxiv.org/html/2609.03229#Sx2.F4)b\)\. On the subset shared with the raw\-measurement table, BioGlyph achieved 72\.4% accuracy, compared with 59\.2% for the table\. The corresponding comparisons were 75\.9% versus 60\.5% for sentences, 72\.8% versus 55\.0% for edge lists and 79\.4% versus 49\.7% for GraphSAGE\-derived roles\. BioGlyph therefore retained an accuracy advantage when both representations were available in full\.
Figure 4:BioGlyph supports network reasoning on ChCh\-Miner, a dense network of 1,514 drugs and 48,514 reported interactions\.Results are pooled across Qwen3\-8B and Llama\-3\.1\-8B\.aSystem accuracy for each representation, with 95% bootstrap CIs\. Percentages within selected bars denote the fraction of prompts that exceeded the context window; the dotted line indicates accuracy without network information\.bPairwise comparisons of BioGlyph with each alternative, restricted to questions for which both representations fitted within the context window and produced an answer\. The number of paired answers is shown at right\.cAblation analysis of the BioGlyph description\. Solid bars show system accuracy; hatched bars show accuracy on the common subset of questions readable under every ablation\.dAccuracy across the seven benchmark question families\.eMedian prompt length plotted against system accuracy\. The vertical dashed line indicates the 24,576\-token context limit\. Connected points show the progression from raw measurements to role names, role names with supporting evidence and the complete BioGlyph description\.Ablation experiments further identified the components underlying this advantage \(Fig\.[4](https://arxiv.org/html/2609.03229#Sx2.F4)c\)\. On the common subset of questions readable under every ablation, role names alone yielded 65\.4% accuracy, and adding the measurements supporting each assignment yielded 66\.2%\. Additionally stating the structural consequence of each role increased accuracy to 74\.6%\. The semantic content of the role names also contributed: replacing each name with a meaningless token while preserving the role assignments and all other prompt content reduced overall system accuracy from 66\.8% to 55\.1%\.
The improvement extended across multiple forms of structural reasoning \(Fig\.[4](https://arxiv.org/html/2609.03229#Sx2.F4)d\)\. For cut\-node questions, BioGlyph achieved 89\.1% accuracy, compared with 46\.7% for sentences and 43\.5% for GraphSAGE\. For questions comparing the degrees of two drugs, BioGlyph achieved 94\.6% and GraphSAGE 95\.7%, whereas sentences achieved 54\.3%\. BioGlyph also outperformed the alternative representations across the principal perturbation and connectivity question families\. Two additional model families showed the same overall pattern on ChCh\-Miner \(Section[S1](https://arxiv.org/html/2609.03229#S1)\)\. A comparison of ChCh\-Miner with the other principal benchmark networks is provided in Fig\.[S1](https://arxiv.org/html/2609.03229#S1.F1)\.
### BioGlyph identifies structurally important proteins in a human pathway network
We next applied BioGlyph to the Reactome functional interaction network\[[48](https://arxiv.org/html/2609.03229#bib.bib45)\], which comprises 10,022 proteins and 194,494 functional interactions, with a mean degree of approximately 39\. Its largest connected component contains 9,823 proteins, including 402 cut nodes whose removal disconnects at least one protein from the remaining main component\. To quantify the structural consequence of each cut node, we removed each one in turn and counted the number of proteins separated from the largest remaining component \(Fig\.[5](https://arxiv.org/html/2609.03229#Sx2.F5)a\)\. Removal of EP300 produced the largest separation, detaching 76 proteins, followed by GPLD1 with 53 and CTCF with 51\. This effect was not determined by the number of interaction partners alone: GPLD1 has 71 partners and detaches 53 proteins, whereas MYC has 382 partners but detaches only 20\.
We compared these structural predictions with gene\-dependency measurements from DepMap 24Q4\[[39](https://arxiv.org/html/2609.03229#bib.bib46),[11](https://arxiv.org/html/2609.03229#bib.bib47)\]\. DepMap provided annotations for 9,387 proteins in the largest connected component across 1,178 cancer cell lines\. Of these proteins, 392 were cut nodes and 8,995 were not\. Selective dependencies constituted 13\.3% of the cut nodes, compared with 8\.8% of the other proteins \(Fig\.[5](https://arxiv.org/html/2609.03229#Sx2.F5)c\)\. The proportions classified as common essentials were similar, at 11\.0% and 10\.8%, respectively\. At the protein level, EP300 was a dependency in 28\.4% of cell lines; CTCF and YY1 were common essentials, whereas GPLD1 was rarely a dependency\. Neither the BioGlyph compiler nor the structural screen used these dependency annotations\.
Figure 5:BioGlyph identifies structurally important proteins in the Reactome functional interaction network\.Panelsa–ccombine compiler outputs with DepMap annotations and do not involve a language model\. Paneldpools results from Qwen3\-8B and Llama\-3\.1\-8B\.aThe 12 cut nodes whose removal detaches the largest number of proteins from the main component\. Bar colours indicate DepMap dependency categories; columns at right report the fraction of cell lines dependent on each protein and its number of network partners\.bProportion of selective dependencies among the highest\-ranked proteins according to the number detached after removal, betweenness or degree\. The dotted line indicates the proportion among all screened proteins\.cProportions of selective dependencies and common essentials among screened cut nodes and all other screened proteins\.dLanguage\-model reasoning across 60 Reactome modules\. Left, the proportion of responses within two proteins of the exact region\-specific detachment count; median absolute errors are shown beneath the bars\. Right, the proportion of checkable structural claims that agree with exact region\-level graph analysis\.We next compared the detachment\-based ranking with rankings based on betweenness and degree \(Fig\.[5](https://arxiv.org/html/2609.03229#Sx2.F5)b\)\. Among the ten highest\-ranked proteins, selective dependencies constituted 20\.0% of the detachment\-ranked set, 50\.0% of the betweenness\-ranked set and 40\.0% of the degree\-ranked set\. Among the top 300 proteins, the corresponding proportions were 15\.0%, 15\.7% and 13\.0%\. Thus, the three criteria prioritized partially distinct sets of proteins: detachment quantifies the size of the network region disrupted by removal, whereas betweenness and degree capture different aspects of structural prominence\. We then asked whether language models could use BioGlyph descriptions to reason about protein\-removal effects\. We compiled 60 Reactome modules and asked each model how many proteins would become detached after removal of a named protein \(Fig\.[5](https://arxiv.org/html/2609.03229#Sx2.F5)d\)\. A response was considered correct when its estimate was within two proteins of the exact count for that region\. Neither the BioGlyph descriptions nor the raw\-measurement tables exceeded the context limit in this analysis, allowing every prompt to reach the models\.
Across 120 responses, BioGlyph yielded 42\.5% accuracy and a median absolute error of one protein\. By comparison, sentences yielded 24\.2% accuracy with a median error of seven proteins, edge lists yielded 11\.7% with a median error of 13 and raw graph measurements yielded 1\.7% with a median error of 27\. The models supplied a numerical count in 74\.2% of BioGlyph responses, compared with 32\.5% of responses based on raw measurements\. When analysis was restricted to responses containing a count, accuracy was 57\.3% with BioGlyph and 5\.1% with raw measurements\. We also compared the models’ qualitative structural claims with exact graph analysis\. Among checkable claims, 86\.5% of those generated from BioGlyph descriptions were correct, compared with 79\.0% from raw measurements, 72\.8% from edge lists and 69\.1% from sentences\. The models recovered several exact perturbation effects, including the detachment of 54 proteins after removal of RAD21 and the exact counts for E2F1 and GATA1\. Further protein\-level examples and controls are reported in Fig\.[S3f](https://arxiv.org/html/2609.03229#S2.F3e)and Section[S1](https://arxiv.org/html/2609.03229#S1)\.
### BioGlyph improves structural reasoning in a dense social network
We next evaluated BioGlyph on ego\-Facebook, a SNAP network comprising 4,039 individuals and 88,234 friendship links\[[29](https://arxiv.org/html/2609.03229#bib.bib50)\]\. The network has a mean degree of approximately 44; all nodes belong to a single connected component, and only 11 are cut nodes\. We evaluated 147 questions using the same representations and two frozen 8\-billion\-parameter language models as in the preceding experiments\. Of these questions, 10\.2% referred to an individual absent from the network and therefore required the model to state that they could not be answered\. Across both models, BioGlyph achieved 79\.3% system accuracy \(95% confidence interval \(CI\), 75\.9–82\.5%; Fig\.[6](https://arxiv.org/html/2609.03229#Sx2.F6)a\), compared with 49\.1% without network information\. GraphSAGE\-derived roles yielded 54\.6% accuracy\. Direct representations of the retrieved graph performed less well: accuracy was 39\.5% for edge lists, 41\.3% for sentences and 29\.3% for raw graph measurements\.
Prompt lengths again differed substantially among representations\. The median raw\-measurement prompt contained 35,161 tokens, exceeding the 24,576\-token context limit; consequently, 63\.9% of these prompts could not be passed to the model\. The corresponding proportions were 35\.5% for sentence prompts and 24\.0% for edge\-list prompts\. By contrast, the median BioGlyph and GraphSAGE descriptions contained 6,012 and 6,547 tokens, respectively, and both representations remained within the context window\. To distinguish representational quality from context\-window truncation, we restricted pairwise analyses to questions for which both representations fitted within the context window and produced an answer \(Fig\.[6](https://arxiv.org/html/2609.03229#Sx2.F6)b\)\. On the subset shared with the raw\-measurement table, BioGlyph achieved 78\.6% accuracy, compared with 65\.5% for the table\. The corresponding comparisons were 73\.1% versus 58\.5% for sentences, 76\.2% versus 47\.9% for edge lists and 79\.5% versus 52\.7% for GraphSAGE\-derived roles\. BioGlyph therefore retained an accuracy advantage when both representations were available in full\.
Figure 6:BioGlyph improves structural reasoning on ego\-Facebook, a network of 4,039 individuals and 88,234 friendship links\.Results are pooled across Qwen3\-8B and Llama\-3\.1\-8B\.aSystem accuracy at two retrieval sizes, with regions containing at most 120 or 300 individuals\. Bars show results for each retrieval size; diamonds show pooled accuracy with 95% bootstrap CIs\. Percentages within selected bars denote the fraction of prompts that exceeded the context window\. The dotted line indicates accuracy without network information\.bPairwise comparisons of BioGlyph with each alternative, restricted to questions for which both representations fitted within the context window and produced an answer\. The number of paired answers is shown at right\.cAblation analysis of the BioGlyph description\. Solid bars show system accuracy; hatched bars show accuracy on the common subset of questions readable under every ablation\.dAccuracy across the seven benchmark question families\.Ablation experiments identified the components contributing to this advantage \(Fig\.[6](https://arxiv.org/html/2609.03229#Sx2.F6)c\)\. On the common subset of questions readable under every ablation, role names alone yielded 71\.5% accuracy, and adding the measurements supporting each assignment yielded the same accuracy\. Additionally stating the structural consequence of each role increased accuracy to 80\.8%\. The semantic content of the role names also contributed: replacing each name with a meaningless token while preserving the role assignments and all other prompt content reduced overall system accuracy from 71\.3% to 57\.7%\. The largest gains occurred on questions concerning network connectivity \(Fig\.[6](https://arxiv.org/html/2609.03229#Sx2.F6)d\)\. For cut\-node questions, BioGlyph achieved 93\.5% accuracy, compared with 45\.7% for sentences and 37\.0% for GraphSAGE\. For questions comparing the degrees of two individuals, BioGlyph and GraphSAGE both achieved 95\.7%\. BioGlyph further achieved 74% accuracy on candidate\-cut questions, 75% on fragmentation questions, 72% on edge\-removal questions and 84% on reachability questions\.
Performance at the two retrieval sizes showed how the representations responded as the retrieved region expanded \(Fig\.[6](https://arxiv.org/html/2609.03229#Sx2.F6)a\)\. When the maximum region size increased from 120 to 300 individuals, sentence accuracy decreased from 55\.8% to 26\.9%, while the proportion of sentence prompts exceeding the context window increased from 12\.2% to 58\.8%\. Raw\-measurement accuracy similarly decreased from 43\.5% to 15\.0%\. By contrast, GraphSAGE accuracy was 52\.7% and 56\.5% at the two retrieval sizes, and BioGlyph accuracy remained stable at 78\.6% and 79\.9%, respectively\. Results for the remaining benchmark networks and additional controls are reported in Section[S1](https://arxiv.org/html/2609.03229#S1)\.
### BioGlyph supports multi\-turn reasoning on an email network
We finally evaluated BioGlyph on email\-Eu\-core, a communication network from a European research institution\[[29](https://arxiv.org/html/2609.03229#bib.bib50)\]\. The network comprises 986 individuals and 16,064 email links, with a mean degree of approximately 33\. All nodes belong to a single connected component; 73 nodes are cut nodes and 95 links are bridges\. We evaluated 156 questions using the same six representations and two frozen 8\-billion\-parameter language models as in the preceding experiments\. Of these questions, 10\.3% referred to an individual absent from the network and therefore required the model to state that they could not be answered\. Across both models, BioGlyph achieved 78\.8% system accuracy \(95% confidence interval \(CI\), 75\.6–81\.9%; Fig\.[7](https://arxiv.org/html/2609.03229#Sx2.F7)a\), compared with 49\.2% without network information and 48\.4% with GraphSAGE\-derived roles\. Direct representations of the retrieved graph performed less well: accuracy was 41\.8% for edge lists, 39\.9% for sentences and 26\.0% for raw graph measurements\.
Prompt lengths again differed substantially among representations\. The raw\-measurement table exceeded the 24,576\-token context limit in 66\.3% of cases, compared with 41\.3% for sentence prompts and 22\.8% for edge\-list prompts\. By contrast, the median BioGlyph description contained 5,825 tokens, and every BioGlyph prompt fitted within the context window\. To distinguish representational quality from context\-window truncation, we restricted pairwise analyses to questions for which both representations fitted within the context window and produced an answer \(Fig\.[7](https://arxiv.org/html/2609.03229#Sx2.F7)b\)\. Among the 147 paired answers shared with the raw\-measurement table, BioGlyph achieved 73\.5% accuracy, compared with 62\.6% for the table\. The corresponding comparisons were 74\.4% versus 63\.0% for sentences, 74\.0% versus 51\.4% for edge lists and 79\.1% versus 46\.2% for GraphSAGE\-derived roles\. BioGlyph therefore retained an accuracy advantage when both representations were available in full\.
Figure 7:BioGlyph supports structural and multi\-turn reasoning on email\-Eu\-core\.Results are pooled across Qwen3\-8B and Llama\-3\.1\-8B\.aSystem accuracy for each representation, with 95% bootstrap CIs\. Percentages within selected bars denote the fraction of prompts that exceeded the context window; the dotted line indicates accuracy without network information\.bPairwise comparisons of BioGlyph with each alternative, restricted to questions for which both representations fitted within the context window and produced an answer\. The number of paired answers is shown at right\.cAblation analysis of the BioGlyph description\. Solid bars show system accuracy; hatched bars show accuracy on the common subset of questions readable under every ablation\.dAccuracy across the seven benchmark question families\.eTurn\-level accuracy in multi\-turn conversations\. The left panel shows a three\-turn study using either BioGlyph or raw graph measurements\. The right panel shows eight\-turn conversations supported by a single BioGlyph description provided at the first turn; turn 8 challenges an earlier answer\.fOutcomes after the model was told that its preceding answer was incorrect, distinguishing cases in which the model retained a correct answer, corrected an incorrect answer, abandoned a correct answer or remained incorrect\.Ablation experiments identified the components contributing to this advantage \(Fig\.[7](https://arxiv.org/html/2609.03229#Sx2.F7)c\)\. On the common subset of questions readable under every ablation, role names alone yielded 62\.3% accuracy\. Adding the measurements supporting each assignment increased accuracy to 64\.8%, whereas additionally stating the structural consequence of each role increased it to 74\.6%\. The semantic content of the role names also contributed: replacing each name with a meaningless token while preserving the role assignments and all other prompt content reduced overall system accuracy from 68\.6% to 54\.8%\. The largest gains occurred on questions concerning network connectivity \(Fig\.[7](https://arxiv.org/html/2609.03229#Sx2.F7)d\)\. For cut\-node questions, BioGlyph achieved 93\.5% accuracy, compared with 41\.3% for sentences and 29\.3% for GraphSAGE\. BioGlyph further achieved 85\.9% accuracy on fragmentation questions and 83\.0% on edge\-removal questions\. Degree comparisons were comparatively simple: BioGlyph, GraphSAGE and both shorter BioGlyph variants achieved 100% accuracy on this question family\.
We next asked whether a single BioGlyph description could support reasoning across successive conversation turns \(Fig\.[7](https://arxiv.org/html/2609.03229#Sx2.F7)e,f\)\. In a three\-turn study, the network region was provided with the initial question, followed by a related question and then a challenge to the preceding answer\. BioGlyph yielded accuracies of 59\.9%, 71\.7% and 52\.5% across the three turns, respectively\. Raw graph measurements yielded 62\.2%, 55\.6% and 37\.6%\. Thus, although the two representations performed similarly on the initial turn, BioGlyph better supported the follow\-up question and the subsequent challenge\.
We then conducted 219 eight\-turn conversations for each model, providing a single BioGlyph description only at the first turn\. Accuracy across all graded turns was 75\.2%\. When an earlier answer was challenged, the models retained a correct answer in 39\.2% of cases and corrected an incorrect answer in 9\.7% \(Fig\.[7](https://arxiv.org/html/2609.03229#Sx2.F7)f\)\. In the three\-turn study, BioGlyph preserved 49\.2% of initially correct answers after a challenge, compared with 33\.5% for raw measurements\. Performance was also insensitive to question phrasing: rewriting the benchmark questions in ordinary prose changed BioGlyph accuracy from 68\.6% to 67\.9%, while raw\-measurement accuracy changed from 27\.5% to 29\.3%\. Complete conversation examples are provided in Section[S2](https://arxiv.org/html/2609.03229#S2)\.
## Discussion
Networks remain challenging objects for language models\. Although these models process natural language fluently, they often fail on questions that require structural reasoning, such as predicting the consequences of removing a node or edge, identifying elements that connect otherwise distinct regions, or determining which elements are important for network integrity\. Here we show that these limitations can be reduced by changing the representation supplied to the model\. Rather than presenting the network directly or requiring the model to interpret a table of graph measurements, BioGlyph compiles topology into descriptions that identify structural roles, provide the evidence supporting each assignment and state its structural consequence\. These descriptions enabled small, frozen language models to reason across 20 biological, social, information and infrastructure networks\. The largest performance improvement obtained by BioGlyph over existing methods occurred in dense or community\-structured networks, for which direct representations were often both difficult to interpret and too long for the context window\.
The performance improvement of BioGlyph were due to two unique features of the compiled representation \(Section[S4](https://arxiv.org/html/2609.03229#S4)\)\. The first was selective compression\. By determining which structural information was relevant to report, the compiler produced substantially shorter descriptions than edge lists, adjacency sentences or raw\-measurement tables\. The compression allowed nearly all BioGlyph prompts to reach the model\. Direct graph representations, by contrast, frequently exceeded the context limit, thus language models failed to produce coherent response\. The second feature was semantic interpretation\. Instead of adding the measurements supporting a role assignment, producing little improvement over role names alone, we explicitly stated the structural consequence of the role consistently which resulted in substantial accuracy improvement\. Replacing meaningful role names with arbitrary tokens also reduced performance, indicating that both the role vocabulary and the accompanying consequence statements help the model use the encoded structure\.
This representation also supported forms of interaction that are difficult to obtain from a table\. A single BioGlyph description could be supplied once and used throughout an eight\-turn conversation, including follow\-up questions and challenges to earlier answers \(Section[S7](https://arxiv.org/html/2609.03229#S7)\)\. Performance changed little when benchmark questions were rewritten in ordinary prose, suggesting that the descriptions were not narrowly coupled to the templates used to generate the benchmark\. Raw measurements supported these interactions less effectively and, when the initial prompt exceeded the context window, could not support a conversation at all\. The models nevertheless sometimes abandoned correct answers when challenged, showing that access to an informative representation does not eliminate the conversational instabilities of the underlying language model\.
Comparisons with learned graph representations further distinguish BioGlyph from approaches that encode topology without assigning explicit structural meaning\. Unsupervised roles produced by four graph encoders performed close to the no\-network baseline, and incorporating these roles into the compiled description did not improve its accuracy\. Supplying continuous node embeddings as soft prompts\[[32](https://arxiv.org/html/2609.03229#bib.bib13)\]also performed no better than supplying embeddings from mismatched nodes \(Section[S9](https://arxiv.org/html/2609.03229#S9)\)\. These findings complement previous work that represents graphs as edge lists, sentences or templates\[[13](https://arxiv.org/html/2609.03229#bib.bib11),[44](https://arxiv.org/html/2609.03229#bib.bib10),[19](https://arxiv.org/html/2609.03229#bib.bib12)\], or trains graph\-aware adapters\[[8](https://arxiv.org/html/2609.03229#bib.bib15),[37](https://arxiv.org/html/2609.03229#bib.bib14)\]\. In our experiments, direct textual renderings were difficult for small frozen models to use, whereas a semantically explicit intermediate representation made the same topology more accessible without modifying either the graph or the language model\.
The structural roles in BioGlyph also had biological relevance independently of language\-model performance\. In the yeast interactome, cross\-community connectors were essential more frequently than the network average and at a rate comparable to, or slightly higher than, equally sized groups selected using individual centrality measures\. Community cores and hubs were also enriched in essential proteins, whereas peripheral proteins were depleted\. These findings extend longstanding observations that network position is associated with gene essentiality\[[25](https://arxiv.org/html/2609.03229#bib.bib5),[22](https://arxiv.org/html/2609.03229#bib.bib18),[51](https://arxiv.org/html/2609.03229#bib.bib6)\]by expressing combinations of structural properties as named, inspectable roles\. Importantly, not every role was associated with essentiality, indicating that the vocabulary distinguishes structural properties rather than imposing a uniformly positive biological interpretation\.
In the Reactome functional interaction network, the same compiler provided an annotation\-independent perturbation screen\. Cut nodes contained a higher proportion of selective cancer dependencies than other screened proteins, although individual proteins varied substantially: some structurally disruptive proteins were broadly or selectively essential, whereas others showed little dependency\. Two additional human interactomes similarly showed associations between selected structural roles and gene essentiality or broad expression, but these associations were not uniform across roles or networks \(Section[S5](https://arxiv.org/html/2609.03229#S5)\)\. These results suggest that compiled structural roles can complement biological annotations by identifying proteins with distinct forms of network importance, rather than serving as direct proxies for biological essentiality\. Because each assignment is linked to explicit measurements and exact graph operations, both the role and any subsequent language\-model claim can be audited against the underlying network analysis\.
A compiled description is not a substitute for a graph algorithm\. When a question can be expressed as a single graph operation, the appropriate algorithm can answer it directly; indeed, exact algorithms supplied the reference answers throughout our benchmark\. The value of the compiled representation lies instead at the interface between algorithmic analysis and natural\-language interaction\. A user may pose a question without specifying the required graph operation, request an explanation that integrates several structural properties or challenge a previous answer\. BioGlyph acts as a front end to classical graph analysis: it executes predefined algorithms, records their outputs and translates them into a closed\-world description that a language model can interpret and explain\. This separation preserves exact computation while enabling a more flexible mode of interaction\.
BioGlyph has several limitations\. First, its performance depends on network density\. In sparse networks, retrieved subgraphs may be small enough to fit directly within the model’s context window; in such cases, adjacency sentences or raw graph measurements can perform comparably to BioGlyph\. Conversely, BioGlyph descriptions occasionally exceeded the prompt budget, reducing accuracy on two of the 20 networks\. Increasing the context window recovered performance on one of these networks \(Section[S3](https://arxiv.org/html/2609.03229#S3)\)\. BioGlyph therefore reduces, but does not always eliminate, the trade\-off between preserving sufficient structural information and limiting prompt length\. Second, language models remain susceptible to misleading feedback\. Although BioGlyph increased the likelihood that models retained a correct answer when challenged, the models still failed to reach the correct answer in some conversations\. Structural claims should therefore be verified using the underlying graph algorithms before they are revised or accepted\. Third, the advantage of BioGlyph varies based on the used model: on four hosted frontier models with sufficient context, the compiled description led the raw\-measurement table in one of twelve comparisons, matched it in four and trailed it in seven \(Section[S6](https://arxiv.org/html/2609.03229#S6)\), and Qwen3\-32B closes the controlled gap to parity \(Section[S3](https://arxiv.org/html/2609.03229#S3)\)\. BioGlyph therefore substitutes for a selection that larger models increasingly perform unaided, and is most useful for the open models that can be run locally\.
Our benchmarking approach also has limitations\. First, the scored benchmark was restricted to questions with answers that could be determined exactly\. We additionally examined more open\-ended biological questions, but these often require the integration of multiple structural criteria with experimental or clinical evidence and could not be assessed using the same grading procedure\. Second, retrieval for one family of benchmark questions used the predicate that was subsequently evaluated during grading\. An oracle\-free retrieval analysis preserved the ordering of the representations, indicating that this information was unlikely to account for the main comparison \(Section[S10](https://arxiv.org/html/2609.03229#S10)\)\. Nevertheless, future benchmarks should fully separate retrieval from evaluation\. Finally, all language models were kept frozen\. This design isolated the effect of network representation and avoided task\-specific adaptation, but it may also have limited the performance attainable from BioGlyph descriptions\. In a preregistered screen evaluated against biological outcomes rather than graph\-derived answers, a frozen BioGlyph reader provided competitive results to the best deterministic decision rule \(Section[S12](https://arxiv.org/html/2609.03229#S12)\)\.
These limitations motivate three extensions of our study\. First, training models on compiled descriptions could determine how much of the remaining error arises from the representation and how much from the frozen reader\. Second, using the compiler to generate structurally justified candidates for subsequent evaluation by a language model could support biological investigation without asking the model to replace deterministic analysis\. Third, an agent could invoke exact graph algorithms whenever a question reduces to a defined operation and use the compiled description for interpretation, explanation and multi\-step reasoning\. Such a system would combine the reliability of graph algorithms with the flexibility of language\-based interaction\. We release the compiler, benchmark questions and exact answers, stored model responses and scored results to enable independent verification and further development\.
## Conclusion
BioGlyph converts network topology into concise, interpretable descriptions that small, frozen language models can use without modification to either the model or the underlying network\. Across networks from multiple domains, these descriptions improved structural reasoning relative to edge lists, adjacency sentences, raw graph measurements and unsupervised graph representations, particularly for dense and community\-structured networks\. Ablation experiments showed that the improvement reflects both selective compression and the explicit communication of structural consequences\. The same descriptions also supported multi\-turn reasoning and produced claims that could be verified against exact graph analysis\. In biological networks, BioGlyph identified structural roles associated with independent measures of gene essentiality and dependency while remaining blind to those annotations\. The framework is therefore best viewed not as a replacement for graph algorithms, but as an interpretable interface between exact network analysis and natural\-language reasoning\. By making structural evidence accessible, concise and auditable, BioGlyph provides a practical foundation for language\-model\-assisted investigation of complex biomedical networks\.
## Methods
BioGlyph converts a network region into a short description of its structural roles\. For each benchmark question, we retrieve the region around the named nodes, compute classical graph quantities on that region, and apply one fixed rule for each of eleven roles\. Each assigned role becomes a line of readable text that is given to a frozen language model\. Exact graph algorithms provide the ground\-truth answers for all graded benchmark questions; no language model grades another language model\. We make stochastic steps reproducible with fixed seeds and deterministic tie\-breaking\.
### Structural signals
LetG=\(V,E\)G=\(V,E\)denote a network withn=\|V\|n=\|V\|nodes, and letG~\\widetilde\{G\}denote its simple undirected view, obtained by removing self\-loops and ignoring edge direction\. We compute connected components, cut nodes, bridges, thekk\-core numberc\(v\)c\(v\), and the community partition onG~\\widetilde\{G\}\. Degreed\(v\)d\(v\), PageRankπ\(v\)\\pi\(v\), node betweennessb\(v\)b\(v\), and edge betweennessb\(e\)b\(e\)are computed onGG, preserving direction when the input graph is directed\. The networks used in the experiments reported here are treated as undirected, soG=G~G=\\widetilde\{G\}\.
We writeN\(v\)N\(v\)for the neighbors of nodevv\. Graph quantities are computed with NetworkX\[[20](https://arxiv.org/html/2609.03229#bib.bib29)\]\. We use Brandes’ algorithm\[[6](https://arxiv.org/html/2609.03229#bib.bib30)\]for node and edge betweenness and the standard core decomposition\[[33](https://arxiv.org/html/2609.03229#bib.bib31)\]forkk\-core number\. Betweenness is computed exactly whenn≤5,000n\\leq 5\{,\}000\. For larger graphs, we estimate node and edge betweenness from1,0001\{,\}000pivot nodes selected with seed00\. The signal record stores whether exact or sampled betweenness was used\.
We partitionG~\\widetilde\{G\}into communitiesℳ\\mathcal\{M\}using the Leiden algorithm\[[38](https://arxiv.org/html/2609.03229#bib.bib32)\]with the modularity objective\[[31](https://arxiv.org/html/2609.03229#bib.bib33)\]and seed00\. We relabel communities by decreasing size and break ties using the smallest node identifier\. This ordering keeps community identifiers stable across repeated runs\. LetC\(v\)C\(v\)denote the community containing nodevv, and let
kv,m=\|\{u∈N\(v\):C\(u\)=m\}\|k\_\{v,m\}=\\left\|\\left\\\{u\\in N\(v\):C\(u\)=m\\right\\\}\\right\|\(1\)be the number of neighbors ofvvthat belong to communitymm\. We compute the participation coefficient of Guimerà and Amaral\[[18](https://arxiv.org/html/2609.03229#bib.bib34)\]as
P\(v\)=1−∑m∈ℳ\(kv,md\(v\)\)2\.P\(v\)=1\-\\sum\_\{m\\in\\mathcal\{M\}\}\\left\(\\frac\{k\_\{v,m\}\}\{d\(v\)\}\\right\)^\{2\}\.\(2\)
We setP\(v\)=0P\(v\)=0whend\(v\)=0d\(v\)=0\. A value near zero indicates that most neighbors belong to one community; larger values indicate neighbors spread across several\. We also compute the intra\-community degree
kint\(v\)=kv,C\(v\),k^\{\\mathrm\{int\}\}\(v\)=k\_\{v,C\(v\)\},\(3\)which counts the neighbors ofvvthat remain inside its own community\.
### The BioGlyph compiler
The compiler assigns eleven structural roles at three levels\. Eight roles describe nodes, two describe edges, and one describes a community\. Each emitted role contains its name, the measurements that triggered the assignment, and a short statement describing the corresponding structural consequence\. Two roles are assigned directly from exact connectivity decompositions\.CUT\_NODEmarks an articulation point ofG~\\widetilde\{G\}, andBRIDGE\_EDGEmarks a bridge\. Both assignments are exact and use no numerical threshold\.
Six roles identify unusually large structural values\. For a structural quantityXX, letμX\\mu\_\{X\}andsXs\_\{X\}denote its mean and standard deviation over the relevant comparison set\. An elementuureceives the role when
X\(u\)\>μX\+σXsX,X\(u\)\>\\mu\_\{X\}\+\\sigma\_\{X\}s\_\{X\},\(4\)whereσX\\sigma\_\{X\}is fixed for that role\. We do not assign a threshold\-based role whensX=0s\_\{X\}=0\. Unless stated otherwise, the comparison distribution is computed within the retrieved region rather than over the full network\.
The six threshold\-based roles are:
- •HUB: degreed\(v\)d\(v\), withσ=2\.0\\sigma=2\.0, compared across all nodes in the region\.
- •AUTHORITY: PageRankπ\(v\)\\pi\(v\), withσ=2\.0\\sigma=2\.0\. We define this role for directed networks; it is not used in the undirected networks reported here\.
- •BOTTLENECK\_LINK: edge betweennessb\(e\)b\(e\), withσ=2\.0\\sigma=2\.0, compared across all edges in the region\.
- •CROSS\_COMMUNITY\_CONNECTOR: node betweennessb\(v\)b\(v\), withσ=2\.0\\sigma=2\.0, together with the requirement thatvvhas neighbors in at least two communities\.
- •BOUNDARY\_NODE: participation coefficientP\(v\)P\(v\), withσ=1\.5\\sigma=1\.5\.
- •COMMUNITY\_CORE: intra\-community degreekint\(v\)k^\{\\mathrm\{int\}\}\(v\), withσ=1\.0\\sigma=1\.0\. The comparison is performed separately within each community containing at least three nodes\.
Three additional roles use fixed structural criteria:
- •ISOLATE: a node whose connected component contains one node\.
- •PERIPHERAL: a node withc\(v\)≤1c\(v\)\\leq 1that belongs to a component containing more than one node\.
- •FRAGILE\_REGION: a community whose induced subgraph has edge connectivityλ≤1\\lambda\\leq 1\. We evaluate edge connectivity for communities containing between33and1,5001\{,\}500nodes\. Communities outside this range receive no edge\-connectivity value\.
For aCUT\_NODE, the evidence field records the number of connected components in the whole region before and after removal, which is not the component\-local count that Eq\. \([9](https://arxiv.org/html/2609.03229#Sx5.E9)\) grades\. For small regions we recompute that count after removal exactly; for larger graphs the field can hold a lower bound instead\. The distinction affects the reported evidence, not the role assignment, which is exact at every network size\. Algorithm[1](https://arxiv.org/html/2609.03229#alg1)summarizes the compiler\.
Algorithm 1BioGlyph compilation of one network region1:Input:region
GG, seed
ss, fixed role parameters
\{σX\}\\\{\\sigma\_\{X\}\\\}
2:Construct the simple undirected view
G~\\widetilde\{G\}
3:Compute structural signals on
GGand
G~\\widetilde\{G\}
4:Partition
G~\\widetilde\{G\}into communities using Leiden with seed
ss
5:Compute
P\(v\)P\(v\)and
kint\(v\)k^\{\\mathrm\{int\}\}\(v\)
6:Initialize the role set
𝒢←∅\\mathcal\{G\}\\leftarrow\\emptyset
7:AssignCUT\_NODEto every articulation point
8:AssignBRIDGE\_EDGEto every bridge
9:AssignISOLATEandPERIPHERALusing their fixed structural criteria
10:foreach threshold\-based roledo
11:Assign the role to elements satisfying Eq\. \([4](https://arxiv.org/html/2609.03229#Sx5.E4)\)
12:endfor
13:foreach community containing
33to
1,5001\{,\}500nodesdo
14:AssignFRAGILE\_REGIONif its edge connectivity is at most
11
15:endfor
16:Attach the triggering measurements and structural consequence to every assigned role
17:Return:
𝒢\\mathcal\{G\}, sorted by level, role name, and target
We render the compiler output at several levels of detail\.*Names*lists each target and its assigned structural roles\.*Names with evidence*adds the measurements that triggered each role\. The*full BioGlyph description*adds the short structural consequence associated with each role\. It also states a closed\-world guarantee: every cut node and every bridge in the region is listed, so a node without aCUT\_NODEentry is not an articulation point of that region\.
The*opaque*control keeps the same assigned roles as the*names*rendering but replaces each semantic role name with a meaningless token\. The control carries no evidence and no structural consequence, which makes it the*names*rendering with the vocabulary withheld\. All BioGlyph renderings also declare the nodes present in the retrieved region, including nodes that receive no structural role\.
### Retrieval and prompt construction
Each benchmark question names one or more nodes of a full network\. We retrieve a region around those targets and render it under each representation\. Named targets are admitted first and remain in the region even above the node cap\. If a target identifier does not exist in the full network, the pipeline marks the question as unanswerable before language\-model inference\. We use two retrieval settings\. The first expands two breadth\-first\-search rings from the named targets and caps the retrieved region at120120nodes\. Within each expansion step, we admit nodes first by increasing graph distance and then by node identifier\. The second setting allows up to300300nodes and routes on the question family\. Six families again begin with two breadth\-first rings\. Counterfactual reachability questions use paths instead of rings: retrieval admits the named targets, the shortest paths between them, the shortest path that survives the removal, and the immediate target neighborhoods, each in deterministic order\.
After initial retrieval, a family\-specific sufficiency check can trigger at most two additional expansion rings while capacity remains\. For question families whose answer depends on the complete local neighborhood of a target, the check asks whether that neighborhood is present\. For counterfactual reachability, the check asks whether the two endpoints remain connected inside the retrieved region once the named node is removed\. The check therefore evaluates the graded predicate itself, so the stopping rule for this one family uses privileged information\. Section[S10](https://arxiv.org/html/2609.03229#S10)reports an oracle\-free retrieval arm that routes by personalized PageRank and never runs the loop\. The oracle described below, rather than the retrieval heuristic, determines whether the retrieved region actually supports the benchmark answer\.
Every representation receives the same block of full\-network target facts\. For each node named by the question, the block states whether the node exists and gives its degree in the full network\. These facts are identical across representations\. The final prompt contains the rendered network region, the shared target\-fact block, and the question\. A fixed system instruction asks the model to reason through the problem and finish with
> ANSWER: <answer\>
for the main benchmark\.
We treat context length as part of system performance\. Each language model is served with a40,96040\{,\}960\-token context window\. We reserve16,38416\{,\}384tokens for generation, leaving a maximum prompt length of24,57624\{,\}576tokens\. Prompt length is measured with the model’s own tokenizer and chat template\. A prompt that exceeds this budget is not sent to the model\.
### Network representations
We hold the retrieved region, question, and language model fixed and vary only how the region is represented\.
- •R0, no graph: the model receives only the shared target facts and the question\.
- •R1, edge list: the region is represented as a list of edges\.
- •R2, sentences: the same adjacency information is expressed as natural\-language sentences\.
- •R3, raw measurements: the region is represented as a table of graph quantities\.
- •R4, learned roles: each node is represented by a role learned with an unsupervised graph encoder\.
- •R5, BioGlyph: the deterministic compiler converts the region into named structural roles with evidence and structural consequences\.
The R3 table is designed to expose the numerical inputs used by the BioGlyph compiler\. For each node it contains degree, normalized degree, betweenness, PageRank, core number, and community identifier\. The remaining node columns are intra\-community degree, cross\-community degree, participation coefficient, number of communities touched, component size, articulation\-point status, and the component count after removal when available\. For each edge it contains edge betweenness and bridge status\. For each community it contains size and edge connectivity when computed\. R3 and R5 are generated from the same structural\-signal record\. Because R3 prints edge\-level quantities for individual edges, those rows also reveal which edges are present\. We retain these columns because edge betweenness and bridge status are inputs to BioGlyph edge roles\.
### Learned\-role baselines
For the learned\-role baseline, each node receives an eleven\-dimensional input vector\. Five graph quantities are also printed with the learned role for interpretation: degree, betweenness, PageRank, core number, and clustering coefficient\. Six additional quantities are used by the graph encoder but are not printed: eigenvector centrality, average neighbor degree, maximum neighbor degree, standard deviation of neighbor degree, log triangle count, and two\-step random\-walk return probability\. The feature set contains no articulation, bridge, community, or removal information, because those quantities are what the BioGlyph compiler derives\. We standardize these features and pass them through two message\-passing layers with hidden width6464, output dimension3232, and a PReLU activation\. We evaluate four architectures: GCN\[[26](https://arxiv.org/html/2609.03229#bib.bib35)\], GraphSAGE\[[21](https://arxiv.org/html/2609.03229#bib.bib36)\], GAT\[[41](https://arxiv.org/html/2609.03229#bib.bib37)\], and GIN\[[49](https://arxiv.org/html/2609.03229#bib.bib16)\], implemented with PyTorch Geometric\[[14](https://arxiv.org/html/2609.03229#bib.bib49)\]\.
We train each encoder without task labels using Deep Graph Infomax\[[42](https://arxiv.org/html/2609.03229#bib.bib38)\]\. Letzvz\_\{v\}denote the embedding of nodevv, and letz~v\\widetilde\{z\}\_\{v\}denote its embedding after corruption by row\-wise permutation of the input feature matrix\. We define the graph summary
s=ς\(1n∑v∈Vzv\),s=\\varsigma\\left\(\\frac\{1\}\{n\}\\sum\_\{v\\in V\}z\_\{v\}\\right\),\(5\)whereς\\varsigmais the logistic function\. With bilinear discriminator
D\(z,s\)=z⊤Ws,D\(z,s\)=z^\{\\top\}Ws,\(6\)the training objective is
ℒ=−1n∑v∈Vlogς\(D\(zv,s\)\)−1n∑v∈Vlog\[1−ς\(D\(z~v,s\)\)\]\.\\mathcal\{L\}=\-\\frac\{1\}\{n\}\\sum\_\{v\\in V\}\\log\\varsigma\\\!\\left\(D\(z\_\{v\},s\)\\right\)\-\\frac\{1\}\{n\}\\sum\_\{v\\in V\}\\log\\left\[1\-\\varsigma\\\!\\left\(D\(\\widetilde\{z\}\_\{v\},s\)\\right\)\\right\]\.\(7\)We train for200200epochs with Adam, a learning rate of0\.010\.01, and seed00\. We then cluster the learned embeddings intoK=8K=8groups using seededkk\-means\+\+\[[3](https://arxiv.org/html/2609.03229#bib.bib39)\]\. Cluster identifiers are relabeled by decreasing cluster size\.
The encoder runs once on the full network\. We standardize features using full\-network means and standard deviations and assign each node to the nearest learned centroid,
r\(v\)=argmink∈\{1,…,K\}‖zv−μk‖22,r\(v\)=\\arg\\min\_\{k\\in\\\{1,\\ldots,K\\\}\}\\left\\\|z\_\{v\}\-\\mu\_\{k\}\\right\\\|\_\{2\}^\{2\},\(8\)whereμk\\mu\_\{k\}denotes centroidkk\. A retrieved region displays each node with its full\-network learned role and the five printed graph quantities, above a legend that names every learned role from the feature profile of its cluster\. The main text reports GraphSAGE; results for all four architectures are provided in the Supplementary Information\.
### Benchmark construction and exact oracle
The benchmark contains seven question families\.*Cut node*asks whether a named node is an articulation point\.*Fragmentation*asks how many pieces remain after removing a named node\.*Edge removal*asks how many pieces remain after removing an edge\.*Reachability*asks whether two nodes remain connected after a third node is removed\.*Compare degree*asks which candidate has the highest degree\.*Candidate cut*asks which candidate causes the largest fragmentation when removed\.*Compositional*asks which candidates satisfy both an articulation condition and a degree condition\.
We target twenty questions per family and network using seed00\. All candidate pools are drawn from the largest connected component\. The cut\-node, fragmentation, and edge\-removal families draw half of their questions from positive cases and half from negative cases\. For reachability we construct both answers directly rather than filtering sampled pairs\. Half of the questions remove an articulation point and take the two endpoints from different components of the remaining graph, so the answer is no\. The remaining questions draw the removed node and the first endpoint from the largest component\. We then take the second endpoint from the nodes lying at distance two or three from the first after the removal, so the answer is yes\. Candidate\-cut and compositional questions include at least one articulation point among the candidates\. Degree\-comparison questions draw one candidate from the highest\-degree decile and the rest uniformly, resampling when necessary so that the maximum\-degree candidate is unique\.
For fragmentation, candidate\-cut, and compositional questions, we restrict sampled targets to degree at most6060\. This constraint allows the target neighborhood to fit inside the largest retrieval region\. Where no candidate meets the bound on a given network, we fall back to the unbounded pool rather than leave the family empty\. We also create unanswerable examples by replacing a referenced node identifier with an identifier that does not occur in the network\. When every quota is available, a network contributes140140answerable questions and1616unanswerable questions\. Some networks contribute fewer because the required structural cases do not exist\. Across eight benchmark networks from five domains, the final benchmark contains1,2391\{,\}239questions\.
We grade every benchmark family with exact graph algorithms\. Letβ\(v\)\\beta\(v\)denote the number of biconnected blocks containing nodevv\. The number of pieces left in the original component after removingvvis
pieces\(v\)=\{0,d\(v\)=0,β\(v\),vis an articulation point,1,otherwise\.\\mathrm\{pieces\}\(v\)=\\begin\{cases\}0,&d\(v\)=0,\\\\ \\beta\(v\),&v\\text\{ is an articulation point\},\\\\ 1,&\\text\{otherwise\}\.\\end\{cases\}\(9\)For an edgeee,
pieces\(e\)=1\+\[eis a bridge\],\\mathrm\{pieces\}\(e\)=1\+\\mathbf\{1\}\\\!\\left\[e\\text\{ is a bridge\}\\right\],\(10\)where𝟏\[⋅\]\\mathbf\{1\}\[\\cdot\]is the indicator function\. These counts refer to the connected component containing the target rather than to unrelated components elsewhere in the network\.
We evaluate each question with two oracle calls\. The full\-network oracle provides the ground\-truth answer used for grading\. The region oracle determines what can be inferred from the retrieved region\. We call a retrieved region sufficient when the answer obtained from that region agrees with the full\-network answer\.
### Language models, inference and grading
We evaluate every representation using the same two frozen 8B language models, Qwen3\-8B\[[50](https://arxiv.org/html/2609.03229#bib.bib41)\]and Llama\-3\.1\-8B\-Instruct\[[17](https://arxiv.org/html/2609.03229#bib.bib42)\], served with vLLM\[[28](https://arxiv.org/html/2609.03229#bib.bib40)\]and seed00\. Qwen3\-8B runs in thinking mode with temperature0\.60\.6, following its recommended configuration, whereas Llama\-3\.1\-8B\-Instruct uses greedy decoding\. For the main benchmark, we read the text following the finalANSWER:marker\. If the marker is absent, we parse the complete reply\. A deterministic rule\-based parser first checks for an explicit statement that the question is unanswerable\. Otherwise, it extracts the answer type required by the question family\. Exact\-match grading is then performed against the graph oracle\. When several answers are equally valid, any oracle\-approved answer is accepted\. A reply with no recoverable answer is scored as incorrect, and so is an over\-budget prompt, which reaches the model as an empty reply\. A generation that reaches the output limit is parsed like any other reply and is incorrect only when no answer can be recovered from it\.
The basic evaluation unit is one prediction for one question, model, retrieval setting, and retrieval cap\. LetUUdenote the set of evaluation units andyu∈\{0,1\}y\_\{u\}\\in\\\{0,1\\\}the correctness of unituu\. System accuracy is
Asys=1\|U\|∑u∈Uyu\.A\_\{\\mathrm\{sys\}\}=\\frac\{1\}\{\|U\|\}\\sum\_\{u\\in U\}y\_\{u\}\.\(11\)A unit is correct when the parsed model answer matches the oracle\. Questions whose missing identifier is detected before model inference are counted as correct for every representation, because the pipeline reports them as unanswerable without calling a model\.
For pairwise comparisons, we also evaluate only units for which both representations fit within the context window and produce a completed reply\. For representationsaaandbb, we define
Uab=\{u∈U:uis answerable, andaandbboth fit and complete\}\.U\_\{ab\}=\\left\\\{u\\in U:\\begin\{array\}\[\]\{l\}u\\text\{ is answerable, and\}\\\\ a\\text\{ and \}b\\text\{ both fit and complete\}\\end\{array\}\\right\\\}\.\(12\)We score the model replies onUabU\_\{ab\}directly\. We recompute this common set for every pairwise comparison\. For the rendering ablation, raw measurements, names, names with evidence, and the full BioGlyph description share one common feasible set so that all steps are evaluated on the same questions\.
### Conversation experiments
#### Three\-turn conversations\.
We send the retrieved region once with the first question, using either BioGlyph or the raw\-measurement representation\. Turn 2 asks a follow\-up about the same region without resending the network representation\. Turn 3 challenges the previous answer with a scripted incorrect alternative derived from the oracle answer\. For a Boolean answer we negate the correct value; for a count we use a neighboring incorrect value; and for a candidate question we name a candidate excluded by the oracle\. The challenge template is fixed within each question family and does not depend on the model’s previous response\. We run the conversation turn by turn and append the visible model response to the conversation history\. Hidden reasoning traces are not added to later prompts\. These experiments reserve8,1928\{,\}192tokens for each model response, leaving a prompt budget of32,76832\{,\}768tokens\.
#### Eight\-turn conversations\.
Each eight\-turn thread starts from one retrieved region and one BioGlyph description\. We use four turn types\. Graded turns have an exact graph\-oracle answer\. Lookup turns ask for information stated directly in the BioGlyph description\. Opinion turns ask for an explanation and are not scored\. Push\-back turns challenge a previous answer with an incorrect alternative and are paired with the earlier turn for before\-and\-after analysis\. Gold answers are computed on the region shown to the model, except for full\-network degree facts that are explicitly supplied in the shared fact block\. The region description is sent only at the first turn\. These experiments reserve4,0964\{,\}096tokens for each response, leaving a prompt budget of36,86436\{,\}864tokens\.
#### Natural\-language questions\.
We asked the same questions about the same retrieved regions, but replaced the benchmark answer format with a request for a short answer in ordinary language\. We extracted the final answer with a deterministic rule\-based parser and graded it with the same exact graph oracle used in the main benchmark\. The parser handles negation and ignores node identifiers or intermediate numbers that are not part of the final answer\.
#### Push\-back outcomes\.
For every completed challenge pair, we record correctness before and after the challenge\. We report four outcomes: the model held a correct answer, corrected a wrong answer, abandoned a correct answer, or remained wrong\. Both turns must complete for the pair to enter this analysis\.
### Biological validation
#### Yeast gene essentiality\.
We match STRING\-Yeast proteins to SGD gene\-essentiality annotations\[[9](https://arxiv.org/html/2609.03229#bib.bib44)\]\. Neither the BioGlyph compiler nor the language models receive these annotations\. For this analysis the compiler runs once on the whole STRING\-Yeast network, so the threshold roles use full\-network comparison distributions rather than regional ones\. For each node\-level structural role, we compare the fraction of essential proteins among nodes carrying that role with the essential fraction across the full labeled network\. We use two\-sided Fisher exact tests and adjust across roles with the Benjamini–Hochberg procedure\[[5](https://arxiv.org/html/2609.03229#bib.bib17)\]\. AdjustedPPvalues are included with the released results\. To compare the cross\-community connector role with individual centrality measures, we construct equally sized groups from the highest\-ranked proteins by degree, betweenness, and PageRank\.
#### Reactome perturbation screen\.
We apply the structural screen to the largest connected component of the Reactome functional interaction network\[[48](https://arxiv.org/html/2609.03229#bib.bib45)\]\. For each proteinvv, we measure how many proteins separate from the largest remaining component aftervvis removed\. LetNNdenote the number of proteins in the original component and𝒞\(G−v\)\\mathcal\{C\}\(G\-v\)the connected components after removal\. We define
s\(v\)=\(N−1\)−maxC∈𝒞\(G−v\)\|C\|\.s\(v\)=\(N\-1\)\-\\max\_\{C\\in\\mathcal\{C\}\(G\-v\)\}\|C\|\.\(13\)For a non\-articulation point,s\(v\)=0s\(v\)=0\. We compare this structural ranking with DepMap 24Q4 gene\-effect data\[[39](https://arxiv.org/html/2609.03229#bib.bib46),[11](https://arxiv.org/html/2609.03229#bib.bib47)\], which reports17,91617\{,\}916screened genes across1,1781\{,\}178cancer cell lines\. The comparisons below are restricted to the proteins of the Reactome main component that DepMap screened\. For each protein, the dependency fraction is the fraction of screened cell lines with Chronos score below−0\.5\-0\.5\.
We use four descriptive dependency categories in Fig\.[5](https://arxiv.org/html/2609.03229#Sx2.F5)\. A protein on the DepMap common\-essential list is labeled*common essential*\. Among proteins not on that list, a dependency fraction below5%5\\%is labeled*rarely a dependency*\. A fraction from5%5\\%to50%50\\%is labeled*selective dependency*, and a fraction above50%50\\%is labeled*a dependency in most lines*\. We evaluate the structural screen in two ways\. First, fork∈\{10,25,50,100,200,300\}k\\in\\\{10,25,50,100,200,300\\\}, we compute the fraction of selective dependencies among the topkkproteins ranked bys\(v\)s\(v\)\. We compare this curve with rankings by betweenness and degree and with the overall selective\-dependency rate\. Second, we compare screened cut nodes with all other screened proteins for selective dependency and common essentiality using one\-sided Fisher exact tests\.
#### Counting perturbation size\.
We select the sixty proteins with the highests\(v\)s\(v\)values and retrieve a module containing between2020and120120proteins around each target\. Each module is rendered as an edge list, adjacency sentences, raw measurements, and BioGlyph\. The prompt names the candidate protein, states its degree within the retrieved module, and asks how many proteins detach after the candidate is removed\.
The exact gold count is computed on the retrieved module rather than on the full Reactome network\. Lety^\\widehat\{y\}be the predicted count andyythe exact module\-level count\. We score a prediction as correct when
\|y^−y\|≤2\.\|\\widehat\{y\}\-y\|\\leq 2\.\(14\)A truncated response or a response from which no count can be recovered is incorrect\. Among replies containing a recoverable count, we also report the median absolute error\|y^−y\|\|\\widehat\{y\}\-y\|\.
We separately evaluate structural statements in the model’s reasoning\. A sentence claiming that removal of one named node disconnects the module can be checked directly against the exact cut nodes of that module\. We evaluate only sentences that refer to exactly one node in the retrieved region\. Correct denials and correct assertions receive equal credit\. IfTTofMMcheckable claims agree with the exact graph analysis, we report
### Statistical analysis
All contrasts are paired on the evaluation unit by an inner join, so a unit attempted by only one representation contributes nothing\. For confidence intervals within one network, we use percentile bootstrap intervals\[[12](https://arxiv.org/html/2609.03229#bib.bib43)\]from3,0003\{,\}000resamples with seed00, resampling evaluation units with replacement\. For effects pooled across networks, we use a two\-stage bootstrap with5,0005\{,\}000resamples and seed00\. We first sample networks with replacement, and we then resample evaluation units with replacement within each sampled network\. The point estimate remains the plain paired mean, and only the interval changes\. The procedure accounts for between\-network variation, because questions drawn from one network come from the same graph and often from the same neighborhoods\. Units that share a question stay correlated inside the second stage, so the interval is not corrected at that finer level of clustering\. For pairwise representation contrasts, we construct the common feasible set of evaluation units before resampling\. We also report the mean per\-network effect with att\-interval across the eight benchmark networks, which leaves seven degrees of freedom\.
### Datasets and implementation
#### STRING\-Yeast\.
We use the budding\-yeast physical interaction network from STRING v12\[[36](https://arxiv.org/html/2609.03229#bib.bib22)\]and retain interactions with combined score at least700700\. Nodes represent proteins and edges represent physical associations\.
#### ChCh\-Miner\.
We use the BioSNAP ChCh\-Miner drug interaction network\[[52](https://arxiv.org/html/2609.03229#bib.bib23)\]\. Nodes represent drugs and edges represent reported drug–drug interactions\.
#### ego\-Facebook\.
We use the ego\-Facebook friendship network from SNAP\[[29](https://arxiv.org/html/2609.03229#bib.bib50)\]\. Nodes represent people and edges represent friendships\.
#### email\-Eu\-core\.
We use the email\-Eu\-core network from SNAP\[[29](https://arxiv.org/html/2609.03229#bib.bib50)\]\. Nodes represent anonymized members of a European research institution and edges represent email interactions\. Department annotations are not used\.
The eight\-network benchmark that yields the1,2391\{,\}239questions comprises STRING\-Yeast, ego\-Facebook, email\-Eu\-core, Wiki\-Vote, Amazon\-Photo, ogbn\-arxiv, Cora, and Coauthor\-CS, drawn from five domains\. ChCh\-Miner is evaluated under the same protocol but sits outside that pool, so its questions are additional to the1,2391\{,\}239\. Section[S1](https://arxiv.org/html/2609.03229#S1)describes the remaining networks and their preprocessing\.
#### Reactome FI\.
We use the Reactome functional interaction network\[[48](https://arxiv.org/html/2609.03229#bib.bib45)\]for the perturbation screen rather than for the main seven\-family benchmark\. We map proteins to DepMap using approved gene symbols\[[40](https://arxiv.org/html/2609.03229#bib.bib48)\]\.
#### Implementation\.
The BioGlyph compiler and graph\-analysis pipeline run on CPU using NetworkX\[[20](https://arxiv.org/html/2609.03229#bib.bib29)\]and igraph\. The benchmark oracle computes cut nodes, bridges, connected components, and perturbation answers exactly\. Only betweenness on graphs larger than5,0005\{,\}000nodes is sampled, as described above\. We use seed00for stochastic procedures, including community detection, sampled betweenness, graph\-encoder initialization,kk\-means clustering, question sampling, language\-model decoding where applicable, and bootstrap resampling\. The Supplementary Information reports the additional networks and graph encoders \(Sections[S1](https://arxiv.org/html/2609.03229#S1)and[S9](https://arxiv.org/html/2609.03229#S9)\), the further model families \(Section[S6](https://arxiv.org/html/2609.03229#S6)\), the retrieval controls \(Section[S10](https://arxiv.org/html/2609.03229#S10)\), the trained\-model reference \(Section[S9](https://arxiv.org/html/2609.03229#S9)\), and the pre\-registered candidate\-knockout screen \(Section[S12](https://arxiv.org/html/2609.03229#S12)\)\.
## Data availability
Every network we analyzed is public\. Section[S1](https://arxiv.org/html/2609.03229#S1)of the Supplementary Information describes the four networks of the main text, and Table[S1](https://arxiv.org/html/2609.03229#S1.T1)gives their sizes\. The protein interaction networks are STRING v12\[[36](https://arxiv.org/html/2609.03229#bib.bib22)\]\(the budding yeast, human and*E\. coli*physical subnetworks\), HuRI\[[30](https://arxiv.org/html/2609.03229#bib.bib19)\]and PP\-Pathways\[[1](https://arxiv.org/html/2609.03229#bib.bib20)\]\. ego\-Facebook, email\-Eu\-core, Wiki\-Vote and the Enron email graph\[[27](https://arxiv.org/html/2609.03229#bib.bib24)\]come from the SNAP collection\[[29](https://arxiv.org/html/2609.03229#bib.bib50)\], and ChCh\-Miner and DG\-AssocMiner come from BioSNAP\[[52](https://arxiv.org/html/2609.03229#bib.bib23)\]\. The western US power grid analyzed by Watts and Strogatz\[[45](https://arxiv.org/html/2609.03229#bib.bib1)\]comes from Newman’s network data collection\. Cora and CiteSeer are the standard citation networks\[[34](https://arxiv.org/html/2609.03229#bib.bib25)\], and Amazon\-Photo and Coauthor\-CS follow Shchur et al\.\[[35](https://arxiv.org/html/2609.03229#bib.bib52)\]; we load those four through PyTorch Geometric\[[14](https://arxiv.org/html/2609.03229#bib.bib49)\]and ogbn\-arxiv through the Open Graph Benchmark package\[[24](https://arxiv.org/html/2609.03229#bib.bib51)\]\. The disease comorbidity network is HuDiNe\[[23](https://arxiv.org/html/2609.03229#bib.bib26)\]\. The glioblastoma patient\-similarity network is built by similarity network fusion\[[43](https://arxiv.org/html/2609.03229#bib.bib28)\]from the public TCGA cohort\[[47](https://arxiv.org/html/2609.03229#bib.bib27)\], and the glioblastoma single\-cell network is built from the Smart\-seq2 expression matrix deposited under GEO accession GSE131928\. The knockout screen runs on the Reactome functional interaction network\[[48](https://arxiv.org/html/2609.03229#bib.bib45)\]\. We validated the roles against external labels, none of which we showed to any model: SGD gene essentiality\[[9](https://arxiv.org/html/2609.03229#bib.bib44)\], DepMap 24Q4 gene effect\[[39](https://arxiv.org/html/2609.03229#bib.bib46),[11](https://arxiv.org/html/2609.03229#bib.bib47)\], GTEx tissue expression\[[10](https://arxiv.org/html/2609.03229#bib.bib21)\]and approved gene symbols from HGNC\[[40](https://arxiv.org/html/2609.03229#bib.bib48)\]\. We release the compiled descriptions, the question sets with their exact answers, every stored model reply and the scored result tables behind every number in this paper together with the analysis code\.
## Code availability
The code used for all analyses is available to editors and reviewers through a private Code Ocean capsule\. The capsule contains the complete BioGlyph implementation, including data preprocessing, representation construction, model training and evaluation pipelines\. An interactive browser\-based demonstration is available at[https://islamlab\.org/bioglyph](https://islamlab.org/bioglyph), allowing users to reproduce the worked examples and apply BioGlyph to their own networks without local installation\. Upon publication, the Code Ocean capsule will be made publicly available under a persistent digital object identifier\.
## References
- \[1\]M\. Agrawal, M\. Zitnik, and J\. Leskovec\(2018\)Large\-scale analysis of disease pathways in the human interactome\.InPacific Symposium on Biocomputing\. Pacific Symposium on Biocomputing,Vol\.23,pp\. 111\.Cited by:[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
- \[2\]R\. Albert, H\. Jeong, and A\. Barabási\(2000\)Error and attack tolerance of complex networks\.nature406\(6794\),pp\. 378–382\.Cited by:[Introduction](https://arxiv.org/html/2609.03229#Sx1.p1.1)\.
- \[3\]D\. Arthur S\. Vassilvitskiiet al\.\(2007\)K\-means\+\+: the advantages of careful seeding\.InSoda,Vol\.7,pp\. 1027–1035\.Cited by:[Learned\-role baselines](https://arxiv.org/html/2609.03229#Sx5.SSx5.p2.4)\.
- \[4\]A\. Barabási and R\. Albert\(1999\)Emergence of scaling in random networks\.science286\(5439\),pp\. 509–512\.Cited by:[Introduction](https://arxiv.org/html/2609.03229#Sx1.p1.1)\.
- \[5\]Y\. Benjamini and Y\. Hochberg\(1995\)Controlling the false discovery rate: a practical and powerful approach to multiple testing\.Journal of the Royal statistical society: series B \(Methodological\)57\(1\),pp\. 289–300\.Cited by:[Yeast gene essentiality\.](https://arxiv.org/html/2609.03229#Sx5.SSx9.SSS0.Px1.p1.1)\.
- \[6\]U\. Brandes\(2001\)A faster algorithm for betweenness centrality\.Journal of mathematical sociology25\(2\),pp\. 163–177\.Cited by:[Structural signals](https://arxiv.org/html/2609.03229#Sx5.SSx1.p2.1)\.
- \[7\]T\. Brown, B\. Mann, N\. Ryder, M\. Subbiah, J\. D\. Kaplan, P\. Dhariwal, A\. Neelakantan, P\. Shyam, G\. Sastry, A\. Askell,et al\.\(2020\)Language models are few\-shot learners\.Advances in neural information processing systems33,pp\. 1877–1901\.Cited by:[Introduction](https://arxiv.org/html/2609.03229#Sx1.p2.1)\.
- \[8\]R\. Chen, T\. Zhao, A\. Jaiswal, N\. Shah, and Z\. Wang\(2024\)Llaga: large language and graph assistant\.arXiv preprint arXiv:2402\.08170\.Cited by:[Introduction](https://arxiv.org/html/2609.03229#Sx1.p4.1),[Discussion](https://arxiv.org/html/2609.03229#Sx3.p4.1)\.
- \[9\]J\. M\. Cherry, E\. L\. Hong, C\. Amundsen, R\. Balakrishnan, G\. Binkley, E\. T\. Chan, K\. R\. Christie, M\. C\. Costanzo, S\. S\. Dwight, S\. R\. Engel,et al\.\(2012\)Saccharomyces genome database: the genomics resource of budding yeast\.Nucleic acids research40\(D1\),pp\. D700–D705\.Cited by:[BioGlyph improves reasoning over the yeast interactome and identifies biologically meaningful roles](https://arxiv.org/html/2609.03229#Sx2.SSx1.p7.1),[Yeast gene essentiality\.](https://arxiv.org/html/2609.03229#Sx5.SSx9.SSS0.Px1.p1.1),[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
- \[10\]G\. Consortium\(2020\)The gtex consortium atlas of genetic regulatory effects across human tissues\.Science369\(6509\),pp\. 1318–1330\.Cited by:[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
- \[11\]J\. M\. Dempster, I\. Boyle, F\. Vazquez, D\. E\. Root, J\. S\. Boehm, W\. C\. Hahn, A\. Tsherniak, and J\. M\. McFarland\(2021\)Chronos: a cell population dynamics model of crispr experiments that improves inference of gene fitness effects\.Genome biology22\(1\),pp\. 343\.Cited by:[BioGlyph identifies structurally important proteins in a human pathway network](https://arxiv.org/html/2609.03229#Sx2.SSx3.p2.1),[Reactome perturbation screen\.](https://arxiv.org/html/2609.03229#Sx5.SSx9.SSS0.Px2.p1.2),[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
- \[12\]B\. Efron\(1992\)Bootstrap methods: another look at the jackknife\.InBreakthroughs in statistics: Methodology and distribution,pp\. 569–593\.Cited by:[Statistical analysis](https://arxiv.org/html/2609.03229#Sx5.SSx10.p1.1)\.
- \[13\]B\. Fatemi, J\. Halcrow, and B\. Perozzi\(2024\)Talk like a graph: encoding graphs for large language models\.InInternational conference on learning representations,Vol\.2024,pp\. 43909–43934\.Cited by:[Introduction](https://arxiv.org/html/2609.03229#Sx1.p2.1),[Introduction](https://arxiv.org/html/2609.03229#Sx1.p4.1),[Discussion](https://arxiv.org/html/2609.03229#Sx3.p4.1)\.
- \[14\]M\. Fey and J\. E\. Lenssen\(2019\)Fast graph representation learning with pytorch geometric\.arXiv preprint arXiv:1903\.02428\.Cited by:[Learned\-role baselines](https://arxiv.org/html/2609.03229#Sx5.SSx5.p1.1),[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
- \[15\]S\. Fortunato\(2010\)Community detection in graphs\.Physics reports486\(3\-5\),pp\. 75–174\.Cited by:[Introduction](https://arxiv.org/html/2609.03229#Sx1.p1.1)\.
- \[16\]M\. Girvan and M\. E\. Newman\(2002\)Community structure in social and biological networks\.Proceedings of the national academy of sciences99\(12\),pp\. 7821–7826\.Cited by:[Introduction](https://arxiv.org/html/2609.03229#Sx1.p1.1)\.
- \[17\]A\. Grattafiori, A\. Dubey, A\. Jauhri, A\. Pandey, A\. Kadian, A\. Al\-Dahle, A\. Letman, A\. Mathur, A\. Schelten, A\. Vaughan,et al\.\(2024\)The llama 3 herd of models\.arXiv preprint arXiv:2407\.21783\.Cited by:[Language models, inference and grading](https://arxiv.org/html/2609.03229#Sx5.SSx7.p1.1)\.
- \[18\]R\. Guimera and L\. A\. Nunes Amaral\(2005\)Functional cartography of complex metabolic networks\.nature433\(7028\),pp\. 895–900\.Cited by:[Structural signals](https://arxiv.org/html/2609.03229#Sx5.SSx1.p3.2)\.
- \[19\]J\. Guo, L\. Du, H\. Liu, M\. Zhou, X\. He, and S\. Han\(2023\)Gpt4graph: can large language models understand graph structured data? an empirical evaluation and benchmarking\.arXiv preprint arXiv:2305\.15066\.Cited by:[Introduction](https://arxiv.org/html/2609.03229#Sx1.p4.1),[Discussion](https://arxiv.org/html/2609.03229#Sx3.p4.1)\.
- \[20\]A\. A\. Hagberg, D\. A\. Schult, and P\. J\. Swart\(2008\)Exploring network structure, dynamics, and function using networkx\.InProceedings of the python in science conference,pp\. 11–15\.Cited by:[Structural signals](https://arxiv.org/html/2609.03229#Sx5.SSx1.p2.1),[Implementation\.](https://arxiv.org/html/2609.03229#Sx5.SSx11.SSS0.Px6.p1.1)\.
- \[21\]W\. Hamilton, Z\. Ying, and J\. Leskovec\(2017\)Inductive representation learning on large graphs\.Advances in neural information processing systems30\.Cited by:[Learned\-role baselines](https://arxiv.org/html/2609.03229#Sx5.SSx5.p1.1)\.
- \[22\]X\. He and J\. Zhang\(2006\)Why do hubs tend to be essential in protein networks?\.PLoS genetics2\(6\),pp\. e88\.Cited by:[Discussion](https://arxiv.org/html/2609.03229#Sx3.p5.1)\.
- \[23\]C\. A\. Hidalgo, N\. Blumm, A\. Barabási, and N\. A\. Christakis\(2009\)A dynamic network approach for the study of human phenotypes\.PLoS computational biology5\(4\),pp\. e1000353\.Cited by:[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
- \[24\]W\. Hu, M\. Fey, M\. Zitnik, Y\. Dong, H\. Ren, B\. Liu, M\. Catasta, and J\. Leskovec\(2020\)Open graph benchmark: datasets for machine learning on graphs\.Advances in Neural Information Processing Systems33,pp\. 22118–22133\.Cited by:[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
- \[25\]H\. Jeong, S\. P\. Mason, A\. Barabási, and Z\. N\. Oltvai\(2001\)Lethality and centrality in protein networks\.Nature411\(6833\),pp\. 41–42\.Cited by:[Introduction](https://arxiv.org/html/2609.03229#Sx1.p1.1),[Discussion](https://arxiv.org/html/2609.03229#Sx3.p5.1)\.
- \[26\]T\. N\. Kipf and M\. Welling\(2016\)Semi\-supervised classification with graph convolutional networks\.arXiv preprint arXiv:1609\.02907\.Cited by:[Learned\-role baselines](https://arxiv.org/html/2609.03229#Sx5.SSx5.p1.1)\.
- \[27\]B\. Klimt and Y\. Yang\(2004\)The enron corpus: a new dataset for email classification research\.InEuropean conference on machine learning,pp\. 217–226\.Cited by:[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
- \[28\]W\. Kwon, Z\. Li, S\. Zhuang, Y\. Sheng, L\. Zheng, C\. H\. Yu, J\. Gonzalez, H\. Zhang, and I\. Stoica\(2023\)Efficient memory management for large language model serving with pagedattention\.InProceedings of the 29th symposium on operating systems principles,pp\. 611–626\.Cited by:[Language models, inference and grading](https://arxiv.org/html/2609.03229#Sx5.SSx7.p1.1)\.
- \[29\]J\. Leskovec and A\. Krevl\(2014\)SNAP Datasets: stanford large network dataset collection\.Note:[http://snap\.stanford\.edu/data](http://snap.stanford.edu/data)Cited by:[BioGlyph improves structural reasoning in a dense social network](https://arxiv.org/html/2609.03229#Sx2.SSx4.p1.1),[BioGlyph supports multi\-turn reasoning on an email network](https://arxiv.org/html/2609.03229#Sx2.SSx5.p1.1),[ego\-Facebook\.](https://arxiv.org/html/2609.03229#Sx5.SSx11.SSS0.Px3.p1.1),[email\-Eu\-core\.](https://arxiv.org/html/2609.03229#Sx5.SSx11.SSS0.Px4.p1.1),[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
- \[30\]K\. Luck, D\. Kim, L\. Lambourne, K\. Spirohn, B\. E\. Begg, W\. Bian, R\. Brignall, T\. Cafarelli, F\. J\. Campos\-Laborie, B\. Charloteaux,et al\.\(2020\)A reference map of the human binary protein interactome\.Nature580\(7803\),pp\. 402–408\.Cited by:[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
- \[31\]M\. E\. Newman\(2006\)Modularity and community structure in networks\.Proceedings of the national academy of sciences103\(23\),pp\. 8577–8582\.Cited by:[Structural signals](https://arxiv.org/html/2609.03229#Sx5.SSx1.p3.1)\.
- \[32\]B\. Perozzi, B\. Fatemi, D\. Zelle, A\. Tsitsulin, M\. Kazemi, R\. Al\-Rfou, and J\. Halcrow\(2024\)Let your graph do the talking: encoding structured data for llms\.arXiv preprint arXiv:2402\.05862\.Cited by:[Introduction](https://arxiv.org/html/2609.03229#Sx1.p4.1),[Discussion](https://arxiv.org/html/2609.03229#Sx3.p4.1)\.
- \[33\]S\. B\. Seidman\(1983\)Network structure and minimum degree\.Social networks5\(3\),pp\. 269–287\.Cited by:[Structural signals](https://arxiv.org/html/2609.03229#Sx5.SSx1.p2.1)\.
- \[34\]P\. Sen, G\. Namata, M\. Bilgic, L\. Getoor, B\. Gallagher, and T\. Eliassi\-Rad\(2008\)Collective classification in network data\.AI magazine29\(3\),pp\. 93–106\.Cited by:[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
- \[35\]O\. Shchur, M\. Mumme, A\. Bojchevski, and S\. Günnemann\(2018\)Pitfalls of graph neural network evaluation\.arXiv preprint arXiv:1811\.05868\.Cited by:[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
- \[36\]D\. Szklarczyk, R\. Kirsch, M\. Koutrouli, K\. Nastou, F\. Mehryary, R\. Hachilif, A\. L\. Gable, T\. Fang, N\. T\. Doncheva, S\. Pyysalo,et al\.\(2023\)The string database in 2023: protein–protein association networks and functional enrichment analyses for any sequenced genome of interest\.Nucleic acids research51\(D1\),pp\. D638–D646\.Cited by:[BioGlyph improves reasoning over the yeast interactome and identifies biologically meaningful roles](https://arxiv.org/html/2609.03229#Sx2.SSx1.p1.1),[STRING\-Yeast\.](https://arxiv.org/html/2609.03229#Sx5.SSx11.SSS0.Px1.p1.1),[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
- \[37\]J\. Tang, Y\. Yang, W\. Wei, L\. Shi, L\. Su, S\. Cheng, D\. Yin, and C\. Huang\(2024\)Graphgpt: graph instruction tuning for large language models\.InProceedings of the 47th International ACM SIGIR Conference on Research and Development in Information Retrieval,pp\. 491–500\.Cited by:[Introduction](https://arxiv.org/html/2609.03229#Sx1.p4.1),[Discussion](https://arxiv.org/html/2609.03229#Sx3.p4.1)\.
- \[38\]V\. A\. Traag, L\. Waltman, and N\. J\. Van Eck\(2019\)From louvain to leiden: guaranteeing well\-connected communities\.Scientific reports9\(1\),pp\. 5233\.Cited by:[Structural signals](https://arxiv.org/html/2609.03229#Sx5.SSx1.p3.1)\.
- \[39\]A\. Tsherniak, F\. Vazquez, P\. G\. Montgomery, B\. A\. Weir, G\. Kryukov, G\. S\. Cowley, S\. Gill, W\. F\. Harrington, S\. Pantel, J\. M\. Krill\-Burger,et al\.\(2017\)Defining a cancer dependency map\.Cell170\(3\),pp\. 564–576\.Cited by:[BioGlyph identifies structurally important proteins in a human pathway network](https://arxiv.org/html/2609.03229#Sx2.SSx3.p2.1),[Reactome perturbation screen\.](https://arxiv.org/html/2609.03229#Sx5.SSx9.SSS0.Px2.p1.2),[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
- \[40\]S\. Tweedie, B\. Braschi, K\. Gray, T\. E\. Jones, R\. L\. Seal, B\. Yates, and E\. A\. Bruford\(2021\)Genenames\. org: the hgnc and vgnc resources in 2021\.Nucleic acids research49\(D1\),pp\. D939–D946\.Cited by:[Reactome FI\.](https://arxiv.org/html/2609.03229#Sx5.SSx11.SSS0.Px5.p1.1),[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
- \[41\]P\. Veličković, G\. Cucurull, A\. Casanova, A\. Romero, P\. Lio, and Y\. Bengio\(2017\)Graph attention networks\.arXiv preprint arXiv:1710\.10903\.Cited by:[Learned\-role baselines](https://arxiv.org/html/2609.03229#Sx5.SSx5.p1.1)\.
- \[42\]P\. Veličković, W\. Fedus, W\. L\. Hamilton, P\. Liò, Y\. Bengio, and R\. D\. Hjelm\(2018\)Deep graph infomax\.arXiv preprint arXiv:1809\.10341\.Cited by:[Learned\-role baselines](https://arxiv.org/html/2609.03229#Sx5.SSx5.p2.1)\.
- \[43\]B\. Wang, A\. M\. Mezlini, F\. Demir, M\. Fiume, Z\. Tu, M\. Brudno, B\. Haibe\-Kains, and A\. Goldenberg\(2014\)Similarity network fusion for aggregating data types on a genomic scale\.Nature methods11\(3\),pp\. 333–337\.Cited by:[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
- \[44\]H\. Wang, S\. Feng, T\. He, Z\. Tan, X\. Han, and Y\. Tsvetkov\(2023\)Can language models solve graph problems in natural language?\.Advances in Neural Information Processing Systems36,pp\. 30840–30861\.Cited by:[Introduction](https://arxiv.org/html/2609.03229#Sx1.p2.1),[Introduction](https://arxiv.org/html/2609.03229#Sx1.p4.1),[Discussion](https://arxiv.org/html/2609.03229#Sx3.p4.1)\.
- \[45\]D\. J\. Watts and S\. H\. Strogatz\(1998\)Collective dynamics of ‘small\-world’networks\.nature393\(6684\),pp\. 440–442\.Cited by:[Introduction](https://arxiv.org/html/2609.03229#Sx1.p1.1),[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
- \[46\]J\. Wei, X\. Wang, D\. Schuurmans, M\. Bosma, F\. Xia, E\. Chi, Q\. V\. Le, D\. Zhou,et al\.\(2022\)Chain\-of\-thought prompting elicits reasoning in large language models\.Advances in neural information processing systems35,pp\. 24824–24837\.Cited by:[Introduction](https://arxiv.org/html/2609.03229#Sx1.p2.1)\.
- \[47\]J\. N\. Weinstein, E\. A\. Collisson, G\. B\. Mills, K\. R\. Shaw, B\. A\. Ozenberger, K\. Ellrott, I\. Shmulevich, C\. Sander, and J\. M\. Stuart\(2013\)The cancer genome atlas pan\-cancer analysis project\.Nature genetics45\(10\),pp\. 1113–1120\.Cited by:[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
- \[48\]G\. Wu, X\. Feng, and L\. Stein\(2010\)A human functional protein interaction network and its application to cancer data analysis\.Genome biology11\(5\),pp\. R53\.Cited by:[BioGlyph identifies structurally important proteins in a human pathway network](https://arxiv.org/html/2609.03229#Sx2.SSx3.p1.1),[Reactome FI\.](https://arxiv.org/html/2609.03229#Sx5.SSx11.SSS0.Px5.p1.1),[Reactome perturbation screen\.](https://arxiv.org/html/2609.03229#Sx5.SSx9.SSS0.Px2.p1.1),[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
- \[49\]K\. Xu, W\. Hu, J\. Leskovec, and S\. Jegelka\(2018\)How powerful are graph neural networks?\.arXiv preprint arXiv:1810\.00826\.Cited by:[Learned\-role baselines](https://arxiv.org/html/2609.03229#Sx5.SSx5.p1.1)\.
- \[50\]A\. Yang, A\. Li, B\. Yang, B\. Zhang, B\. Hui, B\. Zheng, B\. Yu, C\. Gao, C\. Huang, C\. Lv,et al\.\(2025\)Qwen3 technical report\.arXiv preprint arXiv:2505\.09388\.Cited by:[Language models, inference and grading](https://arxiv.org/html/2609.03229#Sx5.SSx7.p1.1)\.
- \[51\]H\. Yu, P\. M\. Kim, E\. Sprecher, V\. Trifonov, and M\. Gerstein\(2007\)The importance of bottlenecks in protein networks: correlation with gene essentiality and expression dynamics\.PLoS computational biology3\(4\),pp\. e59\.Cited by:[Introduction](https://arxiv.org/html/2609.03229#Sx1.p1.1),[Discussion](https://arxiv.org/html/2609.03229#Sx3.p5.1)\.
- \[52\]M\. Zitnik, R\. Sosič, S\. Maheshwari, and J\. Leskovec\(2018\)BioSNAP Datasets: Stanford biomedical network dataset collection\.Note:[http://snap\.stanford\.edu/biodata](http://snap.stanford.edu/biodata)Cited by:[BioGlyph remains effective on a dense drug\-interaction network](https://arxiv.org/html/2609.03229#Sx2.SSx2.p1.1),[ChCh\-Miner\.](https://arxiv.org/html/2609.03229#Sx5.SSx11.SSS0.Px2.p1.1),[Data availability](https://arxiv.org/html/2609.03229#Sx6.p1.1)\.
## Supplementary Information
## S1The networks of the main text, and the result pooled over the whole benchmark
Table[S1](https://arxiv.org/html/2609.03229#S1.T1)describes the four networks that the main text takes one at a time\. For each network it states what a node is, what an edge means, and how large the network is\. Two of the columns need a word of explanation\. A cut node is a node whose removal breaks its part of the network into separate pieces, and a bridge is an edge that provides the only route between two sides\. We computed both exactly rather than estimating them, so the two counts together show how much of a network rests on a single point\. The yeast interactome contains 363 cut nodes and 553 bridges among its 3,384 proteins\. Ego\-Facebook contains only 11 cut nodes among 4,039 people, so almost nothing in that network depends on one person\. The final column names the label that comes from outside the network itself\. We never showed any of these labels to a model, which is what allows us to use one of them, SGD gene essentiality, as an independent check on the compiled roles\. Table[S2](https://arxiv.org/html/2609.03229#S1.T2)lists all twenty benchmark networks with their sizes\.
Figure[S1](https://arxiv.org/html/2609.03229#S1.F1)shows every representation on all four networks under both ways of scoring, andBioGlyphleads on each of them\. We then pooled the whole benchmark, which asks 1,239 questions over eight networks from five domains and answers each of them with two frozen 8B models, Qwen3\-8B and Llama\-3\.1\-8B\. Across that pool the models answered 70\.6% of the questions correctly fromBioGlyph\(95% CI 69\.4 to 71\.9\)\. When we showed them no network at all they answered 50\.7%, and this blind floor is the number that any representation has to beat\. Most representations do not beat it\. The edge list reached 48\.7%, adjacency text 49\.4%, and the four learned\-role encoders 49\.1 to 50\.0%, so close to the floor as to be indistinguishable from it\. The raw\-metric table fell below the floor at 39\.5%, and length explains why: 41\.7% of its prompts exceeded the 24,576\-token budget and never reached the model at all\.BioGlyphbeat every encoder and the raw\-metric table in all 24 of the network\-by\-model combinations we ran, and it beat adjacency text in 18 of them\. It came within one point of a graph neural network that we trained directly on the benchmark questions \(74\.1 against 74\.8,n=381n=381\), and it stayed below a retrieval ceiling of 82\.6%\.
We also scored the formats on only those questions that both of them could deliver\. Across the pool the models answered 66\.8% of those questions fromBioGlyphand 61\.4% from the raw\-measurement table\. We then unfolded the description into the same three rungs that we used on the yeast interactome and scored all of them on one common set of questions\. The role names alone brought the models to 61\.9%, adding the measurements behind each name left them at 62\.0%, and adding the line that says what removal would cost lifted them to 68\.4%\. The last rung raises the score on 19 of the 20 networks that carry the ladder\. Twelve further networks carry the supporting studies\.
Figure[S2](https://arxiv.org/html/2609.03229#S1.F2)and Table[S3](https://arxiv.org/html/2609.03229#S1.T3)give the pooled result in full: every arm under both ways of scoring, with each arm’s overflow share and prompt length\. Every arm here includes the four learned\-role encoders separately, the five fusion arms and the opaque control\. In each of the three models taken alone,BioGlyphleads and the raw\-measurement table trails every other arm\. The environment sweep locates the six losses to adjacency text: all of them sit on the three sparse citation and co\-authorship graphs, where a retrieved region is small enough to read directly\.
Three model families beyond the two 8B readers took the same benchmark, and the main text keeps them out of its figures so that every panel there rests on one pool\. On the drug interaction network all four families answered far more from aBioGlyphdescription than from any conventional format\. Qwen3\-8B reached 87\.5%, Gemma\-3\-12B 71\.8%, Llama\-3\.1\-8B 67\.9% and Mistral\-Nemo\-12B 57\.7%\. Given a table of raw measurements instead, the same four managed 18\.6%, 18\.6%, 32\.1% and 15\.1%\. Scale helps as well as format\. On the yeast interactome Qwen3\-32B answered 90\.7% of the questions from aBioGlyphdescription and 38\.8% from the measurement table\.
Figure S1:Every representation on the four main\-text networks\.Results pool Qwen3\-8B and Llama\-3\.1\-8B throughout\.aSystem accuracy of the six competing representations; the R4 bar is the mean of the four encoders\. Where more than 5% of an arm’s prompts exceeded the token budget, the share printed inside the bar never reached the model\.bThe same four networks for the threeBioGlyphrendering steps and for the opaque\-name control, which is the role\-name rendering with every role name replaced by a meaningless token\.cThe accuracy difference betweenBioGlyphand each arm ofaandb, on the questions both arms fitted into the window and finished; the R4 row takes the best of the four encoders\.Figure S2:The pooled result over the eight benchmark networks\.Results pool Qwen3\-8B and Llama\-3\.1\-8B unless labeled\.aSystem accuracy of every arm with 95% bootstrap confidence intervals; reference lines mark the no\-network floor, the trained graph neural network and the retrieval ceiling\.bBioGlyphagainst every baseline, on all questions and on the questions both arms fit and finished, with two\-stage intervals over networks\.cSystem accuracy per network\.dThe environment sweep: which arm wins in each network\-by\-model cell\. Table[S3](https://arxiv.org/html/2609.03229#S1.T3)gives every arm for each of the three models separately\.Table S1:The four networks of the main text\.Each row says what a node is and what an edge means, then gives the size of the network\. Nodes and edges are those of the simple undirected graph, and edges per node is the average over all nodes\. A cut node is a node whose removal breaks its part of the network into separate pieces\. A bridge is an edge that is the only route between two sides\. Both are counted exactly, so those two columns say how much of each network rests on a single point\. The last column names the annotation that comes from outside the network itself, and no model ever saw one\. SGD gene essentiality is the only outside label this paper uses as a check, and it marks 30\.9% of the yeast proteins as essential\. The DrugBank identities give the drugs their names, ego\-Facebook carries no annotation at all, and the department labels of email\-Eu\-core play no part in any result reported here\.Table S2:All twenty networks\.Statistics of the simple undirected graph\. Questions are the benchmark items on that network \(a share of them is unanswerable from any retrieved region and is scored on the decline\)\.Table S3:The pooled headline over the eight main networks\.System accuracy over all questions with its 95% bootstrap interval, the share of prompts that never reached the model \(overflow\) and the median prompt length, for the two 8B models pooled and for each model alone; the last two columns are the reasoning\-controlled advantage of R5 BioGlyph over that arm on the answerable questions both arms fitted and finished, with the two\-stage \(network, question\) bootstrap interval, and the number of paired units\.
## S2Talking toBioGlyph: the exchanges behind the numbers
The main text reports accuracies\. This section shows the exchanges that produced them\. Figure[2](https://arxiv.org/html/2609.03229#Sx2.F2)gives one exchange, and Figs\.[S3a](https://arxiv.org/html/2609.03229#S2.F3)to[S3f](https://arxiv.org/html/2609.03229#S2.F3e)give many more, each printed on a page of its own\. Together they show what a graded turn, a lookup, an opinion and a push\-back actually look like, and what a model writes when it is right and when it is wrong\.
Every card in these figures is a stored reply from the runs that this paper scores\. We rewrote nothing\. We removed the markdown, collapsed the bullet lists, and where a reply ran long we cut it at the end of a sentence and marked the cut\. Every verdict is the one the paper’s evaluator gave, never our own reading of the reply\. We also selected the cards deterministically: the releasedqa\_items\.jsonrecords the rule that picked each card, together with the file, the identifier and the model behind it\.
The models saw less than the cards show\. Node identities, such as yeast ORF and gene names or ICD\-9 codes, sit outside the exchange, and so do the outside labels such as SGD essentiality and DepMap dependency\. No model ever received any of them\. Email\-Eu\-core is anonymized at source, so it has no identities to show in the first place\. We sent the compiled description once, with the first question, and never sent it again\.
Turns come in four kinds, and we keep them apart both on the cards and in every number we report\. A graded turn has an exact answer from the oracle, so we score it\. A lookup turn asks for something that the description already states, which tests whether the model can read rather than whether it can reason\. An opinion turn asks why a node matters; it has no right answer, and we never score it\. A push\-back turn asserts a wrong answer, and we score the model’s response against the turn that it disputes\.
Figure[S3a](https://arxiv.org/html/2609.03229#S2.F3)takes the widest view, with six single questions on three of the main networks, each one shown in all four formats side by side\. Every prompt carries its length in tokens, so a reader can see what we gave the model as well as what it answered\. The next three figures each follow one conversation from beginning to end: a yeast link \(Fig\.[S3b](https://arxiv.org/html/2609.03229#S2.F3a)\), a disease\-network node \(Fig\.[S3c](https://arxiv.org/html/2609.03229#S2.F3b)\) and an email\-network link \(Fig\.[S3d](https://arxiv.org/html/2609.03229#S2.F3c)\)\. We picked the first two mechanically, taking the first thread of its kind, in identifier order, that met the stated rule\. We picked the email thread for a different reason: it contains both a misread lookup and a capitulation under push\-back, and the failures belong in the record as much as the successes do\.
Figure[S3e](https://arxiv.org/html/2609.03229#S2.F3d)sets the raw\-metric table besideBioGlyphon the same three\-turn conversation, twice\. In the first conversation the table is already too long at turn 1 and never reaches the model, which is the common case on the yeast interactome\. In the second both formats fit and both answer\. The same figure then asks one question twice, once in benchmark style and once as ordinary prose, because a description that worked only in benchmark phrasing would be of little practical use\. Figure[S3f](https://arxiv.org/html/2609.03229#S2.F3e)closes the section with a knockout table that a biologist can check\. It lists named proteins in Reactome\-FI modules, the oracle’s count of proteins that detach when we remove each one, each model’s answer, and DepMap’s verdict beside it\. Two honest failures follow, and we give the pooled push\-back rates beneath them, so that the failures read as examples of a known frequency rather than as anecdotes\.
Figure S3a:One question, one verdict: what asking BioGlyph looks like\.Six single\-turn questions on three of the main networks \(Qwen3\-8B, khop2:120, benchmark\-style replies\)\. Each card gives the region the model saw, the identities of the named nodes where the network has them, the reply as written, the exact answer, and the same question in four formats with each prompt’s length\. In every drawing the node in question is orange and the colors are the pieces left after the removal\. Selection rules are inqa\_items\.json\.Figure S3b:Eight turns on one yeast link, and the answers hold up under push\-back\.STRING\-Yeast, Qwen3\-8B, BioGlyph\. The four key turns are shown in full and the other four as one\-line chips with their verdicts; the map is sent once with turn 1\. ORF names and SGD essentiality of the named proteins are shown outside the loop; the model never sees them\. Beneath the rail: how Llama\-3\.1\-8B fared on the same thread\.Figure S3c:A disease\-network thread with ICD\-9 identities\.HuDiNe comorbidity network, Qwen3\-8B, BioGlyph; the named nodes are ICD\-9 codes\.Figure S3d:An email\-network thread with a reading error and a capitulation\.email\-Eu\-core, Qwen3\-8B, BioGlyph\. The lookup turn is answered wrongly \(the link is a bottleneck in the map, the model says it is not\) and at turn 8 the model gives up a correct answer under push\-back; Llama\-3\.1\-8B corrected itself on the same thread\. People are anonymized at source, so there are no identity chips\.Figure S3e:The same conversation with the raw\-metric table beside it, and the same question asked as prose\.Top: two three\-turn conversations \(Qwen3\-8B\), BioGlyph on the left and the information\-matched raw\-metric table on the right; the bar beside each arm is the turn\-1 prompt against the 32,768\-token budget of the three\-turn study \(40,960 less 8,192 reserved for each reply\)\. In the first conversation the table never starts; in the second both arms fit and both answer, so only the verdicts are shown\. Bottom: one question asked in benchmark style and as prose\.Figure S3f:A knockout table a biologist can check, and two honest failures\.Top: named proteins in Reactome\-FI modules, the oracle count of proteins that detach when the protein is removed, each model’s answer under BioGlyph \(within two counts as right\), and DepMap’s dependency share, outside the loop\. Bottom: a lookup turn answered wrongly \(the map lists node 28 as a community core, not a hub\), and a capitulation under push\-back\.
## S3The raw\-measurement table: where it fails, why, and what recovers it
The raw\-measurement table is the control a careful reader asks about first, because it carries every fact the compiled description carries and more \(Methods\)\. Figure[S4](https://arxiv.org/html/2609.03229#S3.F4)and Table[S4](https://arxiv.org/html/2609.03229#S3.T4)give its anatomy\. The table exceeds the 24,576\-token prompt budget on 41\.7% of all prompts: 21\.8% at the tighter retrieval setting and 61\.6% at the wider one\. The same regions rendered asBioGlyphexceed the budget on 0\.3%\. The table is the longest rendering on every network, with a median of 19,723 tokens against 6,229 and a 90th percentile of 79,682 against 15,224\. Its median length on the same regions runs from 1\.1 times the description’s on the sparsest network to 7\.0 times on the yeast interactome at the wider setting\.
Restricting the score to the questions an arm fit and finished lifts every long arm, and the table most of all, from 39\.5% to 60\.7%;BioGlyphreaches 69\.5% on its own feasible set\. On the questions the two arms share, the table reaches 61\.4% against 66\.8% forBioGlyph, and the pooled two\-stage interval for the difference excludes zero\. The gap is wider at the tighter retrieval setting than at the wider one \(Table[S4](https://arxiv.org/html/2609.03229#S3.T4)\)\. Among the answerable replies that finished, the table mistakes its reasoning on 18\.0%, against 20\.2% for adjacency text, 24\.4% for the edge list and 13\.6% forBioGlyph\. The table is also the arm most willing to declare a question unanswerable when it is, declining 60\.2% of the finished unanswerable questions against 32\.1% forBioGlyph\. A table that prints a row for every node makes a missing node conspicuous by its absence\. A description lists only the nodes that earn a role, and beyond the region list at its head it offers no such cue\.
Two interventions recover the table, and both recover it to parity rather than beyond\. Qwen3\-32B reads the table at 41\.7% andBioGlyphat 82\.3% over all questions, yet on the questions both arms fit the two stand level, at 77\.0% and 76\.9%\. A 131,072\-token window on the two human interactomes removes the overflow: on PP\-Pathways the table’s overflow falls from 40\.7% to 1\.6% and its accuracy rises from 32\.1% to 52\.2%, level withBioGlyphat 52\.6%\. On HuRI the description itself exceeds the standard budget on 22\.8% of prompts, because the interactome’s hubs earn many roles, and the wide window brings both arms level there as well \(Fig\.[S4](https://arxiv.org/html/2609.03229#S3.F4)\)\. The honest reading is the one the main text gives\. An 8B model reasons somewhat better over named roles than over the same facts as numbers, and far better over a description that fits than over one that does not\. A larger model or a wider window removes the second effect and most of the first\.
Table S4:The raw\-metric baseline, network by network and setting by setting\.Overflow is the share of R3 prompts longer than the 24,576\-token budget \(never sent, scored wrong\); tokens are median prompt lengths; the ratio is R3 over R5 BioGlyph\. System accuracy counts every question; the controlled columns are R3 and R5 on the answerable questions both arms fitted and finished, with the paired difference, its bootstrap interval and the number of units\. Two 8B models pooled\. Below the rule: the same for Qwen3\-32B and for the wider window \(40k = 40,960 tokens, 128k = 131,072\)\.Median tokensSystem accuracy \(%\)Reasoning\-controlled \(%\)NetworkSettingOverflow \(%\)R3R5RatioR3R5R3R5Δ\\Delta\[95% CI\]nnego\-Facebookkhop2:12037\.820,6404,2624\.8043\.578\.663\.878\.3\+14\.5 \[\+6\.6, \+22\.4\]152auto:30090\.156,7108,7476\.5015\.079\.981\.281\.2\+0\.0 \[−\-18\.8, \+18\.8\]16email\-Eu\-corekhop2:12047\.423,8345,0864\.7035\.972\.162\.270\.6\+8\.4 \[\+0\.0, \+16\.8\]119auto:30085\.369,89414,1854\.9016\.085\.664\.385\.7\+21\.4 \[\+0\.0, \+42\.9\]28Wiki\-Votekhop2:12032\.420,7595,1454\.0042\.970\.858\.071\.0\+13\.0 \[\+5\.6, \+20\.4\]162auto:30087\.865,90612,5255\.3014\.775\.661\.961\.9\+0\.0 \[−\-19\.0, \+19\.0\]21Amazon\-Photokhop2:12013\.515,2235,0233\.0047\.171\.551\.267\.6\+16\.4 \[\+10\.3, \+22\.5\]213auto:30080\.838,11410,6713\.6019\.272\.164\.184\.6\+20\.5 \[\+5\.1, \+35\.9\]39STRING\-Yeastkhop2:12043\.920,7134,2454\.9041\.076\.368\.478\.2\+9\.8 \[\+2\.3, \+17\.3\]133auto:30069\.263,9599,0767\.0031\.175\.083\.381\.9−\-1\.4 \[−\-9\.8, \+6\.9\]72ogbn\-arxivkhop2:1200\.39,7957,1341\.4060\.360\.360\.258\.7−\-1\.5 \[−\-6\.2, \+3\.5\]259auto:30036\.920,42514,3691\.4038\.158\.055\.060\.3\+5\.3 \[−\-1\.3, \+11\.9\]151Corakhop2:1200\.07,0076,5921\.1067\.666\.766\.965\.8−\-1\.2 \[−\-6\.2, \+4\.2\]260auto:30010\.66,7056,2131\.1061\.566\.367\.864\.8−\-3\.0 \[−\-9\.6, \+3\.5\]230Coauthor\-CSkhop2:1200\.08,6626,6971\.3059\.060\.358\.058\.4\+0\.4 \[−\-4\.7, \+5\.4\]257auto:30033\.721,08310,3062\.0038\.161\.552\.561\.1\+8\.6 \[\+1\.9, \+15\.4\]162all eightkhop2:12021\.813,7795,19349\.769\.560\.766\.8\+6\.1 \[\+3\.9, \+8\.2\]1,555all eightauto:30061\.633,93611,15529\.371\.763\.066\.9\+3\.9 \[\+0\.6, \+7\.2\]719all eightboth41\.719,7236,22939\.570\.661\.466\.8\+5\.4 \[\+3\.6, \+7\.3\]2,274Qwen3\-32B, all eightboth53\.026,5876,59541\.782\.377\.076\.9−\-0\.1 \[−\-2\.0, \+1\.8\]1,011Llama\-3\.1\-8B, HuRIauto:300, 40k window9\.051\.936\.564\.062\.0−\-2\.0 \[−\-16\.0, \+12\.0\]50Llama\-3\.1\-8B, HuRIauto:300, 128k window0\.055\.859\.058\.156\.4−\-1\.7 \[−\-10\.3, \+6\.0\]117Llama\-3\.1\-8B, PP\-Pathwaysauto:300, 40k window78\.816\.053\.233\.350\.0\+16\.7 \[−\-8\.3, \+41\.7\]24Llama\-3\.1\-8B, PP\-Pathwaysauto:300, 128k window3\.254\.553\.855\.250\.0−\-5\.2 \[−\-12\.9, \+3\.4\]116Figure S4:The raw\-measurement table in full\.Results pool Qwen3\-8B and Llama\-3\.1\-8B unless labeled\.aShare of raw\-measurement prompts over the 24,576\-token budget per network and retrieval setting, withBioGlyphat the same setting\.bMedian prompt length of every arm per network against the budget; Table[S4](https://arxiv.org/html/2609.03229#S3.T4)gives the lengths themselves\.cSystem accuracy and reasoning\-controlled accuracy of every arm, with the share of prompts never sent\.dWhat the finished answerable answers got wrong, per arm, and the share of finished unanswerable questions the model itself declined\.eThe 131,072\-token window on the two human interactomes \(Llama\-3\.1\-8B\): prompts never sent above, system accuracy below, at both windows and both retrieval settings\.
## S4Which layer of the description does the work: the ladder on every network
The ladder unfolds the compiled description into four rungs that carry the same computed information: the raw\-measurement table, role names alone, names with evidence, and the full description with its stated consequences \(Methods\)\. Figure[S5](https://arxiv.org/html/2609.03229#S4.F5)and Tables[S5](https://arxiv.org/html/2609.03229#S4.T5)and[S6](https://arxiv.org/html/2609.03229#S4.T6)give the rungs on all twenty networks, on all questions and on one common feasible set\. On the eight main networks that set holds 2,333 units, and accuracy moves from 63\.6% at the table to 61\.9% with names, 62\.0% with evidence and 68\.4% with the consequence\. Renaming the table as bare roles costs a small model a little, and adding the measurements returns almost nothing\. Stating what each role implies earns the rung its lead, and only that last step’s interval excludes zero\.
The pattern holds beyond the main pool\. The consequence step is positive on nineteen of the twenty networks and is the largest of the three steps on seventeen\. The exceptions are PP\-Pathways, HuDiNe and the single\-cell graph, where the naming or the evidence step is as large or larger\. Across the eight studies of Table[S5](https://arxiv.org/html/2609.03229#S4.T5)the consequence step is positive in seven, and its interval excludes zero in five\. The naming step is positive in three, and its interval excludes zero on the negative side in three\.
The opaque control comes at the same question from the other side\. The control is the names rung line for line, with each role name replaced by a meaningless token, so the contrast that isolates the words is a contrast with that rung\. Eleven networks carry the control, and the named rung leads on every one\. At the extremes, the names lift system accuracy from 55\.1% to 58\.0% on PP\-Pathways and from 51\.1% to 66\.8% on the yeast interactome\. The named prompt is the longer of the two throughout, so the gain is not a saving of room\. Against the full description the control trails on ten of the eleven networks; the exception is HuRI, where the full description exceeds the budget on 25\.6% of prompts and the shorter control does not\. The two readings agree\. Turning a table into bare role names gains nothing by itself, but once the roles are there, whether they arrive as words or as symbols decides most of what the rung is worth\.
Table S5:Which part of the description does the work: the three rendering steps in every study\.Each step is the reasoning\-controlled difference between two rungs of the R5 ladder \(names−\-R3 raw metrics; names \+ evidence−\-names; BioGlyph−\-names \+ evidence\) on one common feasible set of all four rungs, with the two\-stage \(network, question\) bootstrap interval; the total is BioGlyph−\-R3 on the same units\. Asterisks mark intervals that exclude zero\. Fromladder\.parquet\(scripts/ladder\_report\.py\)\.Table S6:The four rungs on every network, two 8B models pooled\.System accuracy of R3 raw metrics, R5 names only, R5 names \+ evidence and R5 BioGlyph, then the same four on the common feasible units of all four rungs \(reasoning\-controlled\), the consequence step \(BioGlyph−\-names \+ evidence\) with its bootstrap interval, and the number of feasible units\. Rows marked s7ext are the harder question set on the two patient graphs\.Figure S5:The ladder on every network\.Results pool Qwen3\-8B and Llama\-3\.1\-8B\. The four rungs are the raw\-measurement table, role names alone, names with the measurements behind them, and the full description with its stated consequences\.a, bSystem accuracy and accuracy on the common feasible set of the four rungs, on all twenty networks, grouped by domain\.cThe consequence step per network with 95% intervals; a filled marker means the interval excludes zero\.dThe full description against the opaque\-name control on every network where the control ran; the control is the names rung with its names replaced, so this gap carries the evidence and the stated consequences as well as the words\. Table[S5](https://arxiv.org/html/2609.03229#S4.T5)gives the three steps for every study\.
## S5The sixteen further networks
Five biological networks beyond the yeast interactome and the drug network carry the same benchmark \(Fig\.[S6](https://arxiv.org/html/2609.03229#S5.F6), Table[S7](https://arxiv.org/html/2609.03229#S5.T7)\)\. The two further STRING interactomes repeat the yeast picture\. On STRING\-EcoliBioGlyphreaches 79\.3% against 39\.3% for the raw table, which loses 54\.2% of its prompts to the budget; the controlled comparison stands at 84\.7% against 76\.6%, with an interval that excludes zero\. On STRING\-Human the system accuracies are 69\.7% against 41\.2%, and the controlled comparison is level at 63\.3% against 62\.3%\.
The two networks assembled from binary or curated human interactions behave differently\. On PP\-PathwaysBioGlyphreaches 57\.4% against 32\.4%, because the table overflows on 45\.8% of prompts, and once both arms fit the two stand level at 51\.1% against 49\.3%\. On HuRI the description itself overflows on 25\.6% of prompts, because the interactome’s hubs earn many roles\. Its system accuracy of 48\.4% falls below the table at 49\.7% and below its own names\-only rung at 60\.7%, and the controlled comparison is level\. The disease–gene network DG\-AssocMiner is the one network where the description trails on all questions by more than an overflow story\. It stands at 48\.9% against 59\.5% for the table and 64\.3% for adjacency text, and level once both fit at 63\.6% against 65\.3%\. The bipartite structure gives a region few of the roles the vocabulary was built for, and we report the miss as it is\. The biological check on HuRI runs in the yeast direction\. Cut nodes are enriched for common\-essential genes, at 10\.4% against a background of 7\.6%, with an odds ratio of 1\.5 \(95% CI 1\.21 to 1\.86\); Fig\.[S7](https://arxiv.org/html/2609.03229#S5.F7)gives the enrichment by role for both interactomes\. A 131,072\-token window brings both arms level on these two interactomes \(Fig\.[S4](https://arxiv.org/html/2609.03229#S3.F4)e\)\.
Figure S6:Where the description helps, by network and by question family\.Results pool Qwen3\-8B and Llama\-3\.1\-8B\. Row colors give the domain throughout\.aSystem accuracy ofBioGlyphminus the better of adjacency text and the raw\-measurement table, for each of the nineteen networks that carry all seven question families\. Red is an advantage to the description, blue to the baseline\.bThe same comparison as one paired contrast per network, on the questions both arms fitted and finished, with 95% intervals, sorted\. Table[S7](https://arxiv.org/html/2609.03229#S5.T7)lists every arm on every network, including the single\-cell graph and the two harder co\-community sets, which carry one family each and so do not appear ina\.Figure S7:The compiled roles of the two human interactomes, against DepMap and GTEx\.Odds ratio per role with its 95% interval, one marker per outside label\. Amber marks enrichment, blue depletion, grey a test that does not separate atp<0\.05p<0\.05\.aHuRI\.bPP\-Pathways, which DepMap does not screen for selective essentiality\. No language model reads anything in this figure\.Table S7:The sixteen further networks: system accuracy of every arm and the controlled advantage of R5 BioGlyph\.Two 8B models pooled \(Qwen3\-8B and Llama\-3\.1\-8B are present in every study\)\. The R4 column is the mean of the four encoders\. The last three columns are R5 BioGlyph minus R3 raw metrics and minus R2 adjacency text on the answerable questions both arms fitted and finished, and R3’s overflow share\. A dash marks an arm not run on that network\.Three graphs built from patients rather than molecules carry the same benchmark \(Table[S7](https://arxiv.org/html/2609.03229#S5.T7)\)\. They are the HuDiNe comorbidity network of ICD\-9 diseases, a glioblastoma patient\-similarity network fused from three omics views of the TCGA cohort, and a single\-cell graph built from glioblastoma cells\. On the two patient\-similarity graphs the retrieved regions are small and regular, every arm fits, and the raw table runs close to the compiled description under both views\. On HuDiNe the system accuracies are 75\.6% against 73\.1% and the controlled ones 76\.5% against 75\.1%\. On TCGA\-GBM\-PSN they are 77\.4% against 75\.2% and 77\.6% against 75\.2%\. The ladder is unusual on both: the naming step is strongly negative, and the evidence step, negligible elsewhere, is one of the largest gains\. The single\-cell graph is hard for every arm, at 34\.1% forBioGlyphagainst 15\.9% for a table that overflows on 61\.4% of prompts, and its controlled contrast rests on 27 units\. Two harder question sets built on the same patient graphs ask about co\-community membership and are answerable from the region far less often\. On those the description leads once both arms fit, at 17\.1% against 7\.1% on the disease network and 35\.2% against 18\.3% on the patient\-similarity network\.
The remaining eight networks split into three further dense graphs and five sparse ones, and the split is the result \(Fig\.[S6](https://arxiv.org/html/2609.03229#S5.F6), Table[S7](https://arxiv.org/html/2609.03229#S5.T7)\)\. On Amazon\-Photo, Wiki\-Vote and email\-Enron the main\-text picture repeats\.BioGlyphreaches 71\.8%, 73\.2% and 79\.6% against a raw table at 33\.2%, 28\.8% and 36\.5% that loses about half of its prompts\. It stays ahead once both arms fit, at 70\.2% against 53\.2%, 69\.9% against 58\.5% and 77\.1% against 71\.7%, and the consequence step runs from 6\.8 to 10\.7 across the three\.
On the five sparse graphs the retrieved region is small, almost nothing overflows, and the raw view is readable\. There the table stands ahead of the description once both fit on Cora, CiteSeer and Power\-Grid, at 67\.3% against 65\.3%, 78\.7% against 75\.4% and 67\.4% against 62\.7%, and only the Power\-Grid interval excludes zero\. Adjacency text is ahead of the description on ogbn\-arxiv on all questions\. The consequence step stays positive on all five sparse graphs, but the naming step is negative on four of them, so the ladder nets to little\. A description helps where the raw view is unmanageable, not everywhere, and we report the absence of an effect as plainly as the effect\.
## S6Model families and scale, and the frontier ceiling
Two further open model families read the eight main networks on a reduced set of arms: adjacency text, the raw table, the four encoders and the three compiled rungs \(Fig\.[S8](https://arxiv.org/html/2609.03229#S6.F8), Table[S8](https://arxiv.org/html/2609.03229#S6.T8)\)\. Table[S8](https://arxiv.org/html/2609.03229#S6.T8)also carries the controls of Sections[S9](https://arxiv.org/html/2609.03229#S9)and[S10](https://arxiv.org/html/2609.03229#S10), so it collects every arm that no figure of its own reports\. Both keep the order\. Gemma\-3\-12B goes from 35\.2% on the raw table and 49\.2% on adjacency text to 68\.0% onBioGlyph, and Mistral\-Nemo\-12B from 30\.1% and 41\.5% to 60\.3%\. Once both arms fit, Mistral keeps its lead over the table, at 56\.3% against 49\.8% with an interval that excludes zero, and over adjacency text at 55\.1% against 46\.6%\. Gemma stands at parity under the same control, at 60\.2% against 62\.6% and 61\.3% against 58\.5%, with both intervals crossing zero\. Qwen3\-32B lifts every arm, leaves the ordering unchanged, and closes the controlled gap between the description and the table to level \(Section[S3](https://arxiv.org/html/2609.03229#S3)\)\.
Four hosted frontier readers, gpt\-5\.5, gpt\-5\.6, Claude Opus 5 and Claude Fable 5, read three of the networks with their own context windows, so nothing overflows\. There the effect is absent\. Every arm scores between 73\.9% and 87\.8%\. The description leads the table in one of the twelve network\-by\-reader cells, matches it in four and trails it in seven, and the gap runs from 6\.6 below to 0\.6 above\. The absence of a gap marks the honest ceiling of the claim\. The compiled description is for the models most people can run, and it is those models that cannot make the compiler’s choices for themselves\.
Figure S8:Model families and scale\.aGemma\-3\-12B and Mistral\-Nemo\-12B on the eight main networks at both retrieval settings; the share printed inside a bar never reached the model\.bQwen3\-32B against the two 8B models, per arm\.cThe four hosted frontier readers on three networks with their own context windows, where nothing overflows and the description and the table read equally well\. Table[S8](https://arxiv.org/html/2609.03229#S6.T8)gives the controlled comparisons for both further families\.Table S8:Further controls\.Top: hosted frontier readers on three networks \(khop2:120; nothing overflows\), system accuracy per arm\. Middle: two further open model families on the eight main networks \(system accuracy pooled over the eight, and the controlled advantage of R5 BioGlyph over R3 and R2\)\. Bottom: fusion of learned roles with BioGlyph \(R6\) against R5 alone, and the oracle\-free PPR retrieval\.
## S7Conversation statistics: threads, the three\-turn study and prose
The exchanges of Section[S2](https://arxiv.org/html/2609.03229#S2)come from three studies whose statistics this section reports in full \(Fig\.[S9](https://arxiv.org/html/2609.03229#S7.F9), Table[S9](https://arxiv.org/html/2609.03229#S7.T9)\)\. The eight\-turn threads run on email\-Eu\-core, STRING\-Yeast and HuDiNe: 719 threads and 11,504 turns over the two models, with no prompt over budget\. Turn classes are never pooled\. Qwen3\-8B answers 87\.8% of the graded turns, 76\.9% of the lookups and 65\.6% of the push\-backs; Llama\-3\.1\-8B answers 61\.4%, 60\.2% and 43\.7%\. Graded accuracy does not fall along a thread\. On turns five to eight against turns one to four, Qwen stands at 92\.3% against 84\.8% and Llama at 65\.2% against 58\.9%, so the dips by position are hard turns rather than late turns\. Under push\-back, scored against the turn each push\-back disputes, Qwen holds a correct answer 59\.1% of the time, gives one up 27\.7% and corrects a wrong one 6\.8%\. For Llama the three shares are 25\.1%, 30\.6% and 18\.4%\.
In the three\-turn study, 368 conversations run once per model and arm, giving 736 model\-conversations per arm\. The raw\-measurement table never starts 131 of its conversations, 17\.8%, because the region does not fit at turn 1;BioGlyphstarts all of them\. The second turn is the load\-bearing one, because the region is not resent\. Asked as a follow\-up,BioGlyphanswers 73\.5% of the questions it answers 75\.5% of when they stand alone\. The table answers 58\.7% as a follow\-up against 72\.6% alone\. A model reading the description still has the map two turns later, and a model reading the table largely does not\.
The style study asks the same questions as prose, with no labeled answer line\. The change movesBioGlyphfrom 71\.4% to 68\.8% and the table from 43\.2% to 41\.1%\. The advantage of the description over the table survives the phrasing\. On the questions both arms fit and finished, the description stands at 75\.3% against 69\.8% in benchmark style and 65\.9% against 60\.1% in prose, with both intervals excluding zero\. The prose extractor recovers an answer from 98\.6% of replies\. Scoring the canonicalized answer rather than the visible text raises measured accuracy from 42\.6% to 65\.6%, which is why the extractor exists \(Table[S9](https://arxiv.org/html/2609.03229#S7.T9)\)\.
Table S9:The conversation studies\.Top: the eight\-turn threads \(built for the compiled description alone, because their lookup turns ask about roles the description names\), accuracy by turn class, which are never pooled, and what the push\-back did against the turn it disputes\. Middle: the three\-turn study, R5 BioGlyph against R3 raw metrics on the same conversations, accuracy at each turn on the conversations that fitted and finished, and the number of R3 conversations that never started\. Bottom: the same questions asked in benchmark style and as prose\.Eight\-turn threads \(719 threads per model\)NetworkModelThreadsGradedLookupPush\-backHeldGave upCorrectedWrong twiceHuDiNeLlama\-3\.1\-8B25064\.671\.650\.630\.928\.519\.720\.9HuDiNeQwen3\-8B25089\.582\.065\.658\.728\.97\.05\.4STRING\-YeastLlama\-3\.1\-8B25059\.452\.043\.421\.727\.721\.728\.9STRING\-YeastQwen3\-8B25083\.878\.469\.262\.323\.17\.37\.3email\-Eu\-coreLlama\-3\.1\-8B21960\.056\.636\.122\.536\.213\.328\.0email\-Eu\-coreQwen3\-8B21990\.469\.461\.656\.031\.56\.06\.5Three\-turn study \(map at turn 1, a facts\-only follow\-up, a push\-back\)NetworkModelArmTurn 1Turn 2Turn 3nnHeldGave upWrong twiceHuDiNeLlama\-3\.1\-8BR3 raw metrics42\.346\.617\.01378\.837\.445\.6HuDiNeLlama\-3\.1\-8BR5 BioGlyph44\.262\.819\.613810\.152\.727\.7HuDiNeQwen3\-8BR3 raw metrics75\.489\.284\.513883\.85\.410\.1HuDiNeQwen3\-8BR5 BioGlyph76\.195\.389\.213885\.89\.51\.4STRING\-YeastLlama\-3\.1\-8BR3 raw metrics50\.034\.36\.1604\.530\.363\.6STRING\-YeastLlama\-3\.1\-8BR5 BioGlyph47\.241\.89\.2894\.137\.853\.1STRING\-YeastQwen3\-8BR3 raw metrics62\.252\.150\.04545\.86\.243\.8STRING\-YeastQwen3\-8BR5 BioGlyph75\.376\.574\.58968\.48\.217\.3email\-Eu\-coreLlama\-3\.1\-8BR3 raw metrics54\.035\.89\.91133\.332\.557\.5email\-Eu\-coreLlama\-3\.1\-8BR5 BioGlyph52\.249\.216\.411310\.738\.545\.1email\-Eu\-coreQwen3\-8BR3 raw metrics76\.189\.986\.86786\.82\.910\.3email\-Eu\-coreQwen3\-8BR5 BioGlyph67\.594\.388\.511487\.76\.64\.9R3 conversations that never started \(turn 1 over budget\): 131 of 736 \(17\.8%\); R5: 0Benchmark style against prose \(same questions; all answerable with overflow scored wrong, and finished replies only\)Networkmodelarmbench\. allprose allbench\. finishedprose finishedHuDiNeLlama\-3\.1\-8BR3 raw metrics60\.050\.068\.949\.2HuDiNeLlama\-3\.1\-8BR5 BioGlyph65\.759\.368\.559\.9HuDiNeQwen3\-8BR3 raw metrics74\.371\.474\.371\.4HuDiNeQwen3\-8BR5 BioGlyph81\.477\.981\.477\.9STRING\-YeastLlama\-3\.1\-8BR3 raw metrics34\.331\.463\.246\.7STRING\-YeastLlama\-3\.1\-8BR5 BioGlyph56\.459\.362\.660\.7STRING\-YeastQwen3\-8BR3 raw metrics35\.735\.082\.080\.3STRING\-YeastQwen3\-8BR5 BioGlyph87\.980\.787\.980\.7email\-Eu\-coreLlama\-3\.1\-8BR3 raw metrics46\.450\.056\.655\.2email\-Eu\-coreLlama\-3\.1\-8BR5 BioGlyph60\.761\.465\.462\.2email\-Eu\-coreQwen3\-8BR3 raw metrics8\.68\.680\.080\.0email\-Eu\-coreQwen3\-8BR5 BioGlyph76\.474\.376\.474\.3Figure S9:The conversation studies\.aGraded accuracy against turn index over the eight\-turn threads, per model and pooled; the dotted lines mark the pooled early and late halves\. Turn 8 is the push\-back and is scored against the turn it disputes, so it is not plotted here\.bThe three\-turn study:BioGlyphagainst the raw\-measurement table at each turn, on the conversations that fitted and finished, with the paired difference above each turn\.cThe same questions in benchmark phrasing and as ordinary prose\. Table[S9](https://arxiv.org/html/2609.03229#S7.T9)gives accuracy by turn class, the push\-back outcomes and the per\-network breakdown\.
## S8The Reactome\-FI knockout screen, in full
The main component of the Reactome functional interaction network holds 9,823 proteins, 402 of them cut nodes, and DepMap screened 97\.5% of them\. Each protein is scored by how much of the network detaches when it is removed \(Methods\), and the ranking is validated against DepMap with no model in the loop \(Fig\.[5](https://arxiv.org/html/2609.03229#Sx2.F5), Table[S10](https://arxiv.org/html/2609.03229#S8.T10)\)\. The top of the ranking is a list of known regulators\. EP300 detaches 76 proteins, then GPLD1 with 53, CTCF with 51, YY1 with 29 and TP53 with 26; Table[S10](https://arxiv.org/html/2609.03229#S8.T10)gives each one’s DepMap class\. Cut nodes are enriched for selective essentiality relative to the rest of the screened proteins, at 13\.3% against 8\.8% \(one\-sided Fisherp=2\.4×10−3p=2\.4\\times 10^\{\-3\}\)\.
The language\-model leg poses the detachment count for the 60 modules to three models under each rendering; the main text reports the two 8B readers, and Table[S10](https://arxiv.org/html/2609.03229#S8.T10)adds Qwen3\-32B beside them\. Pooled over the three models, the description is within two of the oracle on 47\.8% of cases with a median error of 2\.0\. Its node\-level disconnection claims agree with the exact analysis 87\.1% of the time\. Adjacency text stands at 36\.7%, the edge list at 22\.2% and the raw table at 10\.6%\. The RAD21 case quoted in the Results is one of these modules\. All three models return the oracle’s 54 for the retrieved module, whose whole\-network separation is 19, and DepMap records RAD21 as a dependency in 99\.8% of 1,178 lines\.
Table S10:Biological validation\.Top: SGD essentiality of yeast proteins by BioGlyph role, with the odds ratio against the background rate, and the essential share of the same number of proteins picked by a single centrality\. Middle: the Reactome\-FI knockout screen, deterministic leg \(top\-ranked cut nodes and how many proteins they detach\) and the language\-model leg \(each arm asked how many proteins detach when the named protein is removed from the retrieved module; a count within±\\pm2 of the oracle is scored right\)\.
## S9Learned roles, fusion and the trained yardstick
The four learned\-role encoders sit near the blind floor on every network \(Fig\.[S10](https://arxiv.org/html/2609.03229#S9.F10)\), and fusing their roles into the compiled description changes nothing\. The five fusion arms carry one encoder each, plus one with the compiler’s own whole\-graph roles\. They score between 69\.6% and 70\.0% pooled against 70\.6% forBioGlyphalone, and none of them is ahead under the controlled view \(Table[S8](https://arxiv.org/html/2609.03229#S6.T8)\)\. A further control hands the model the encoder’s embedding directly, as one continuous token through a trained projector, with the language model frozen\. The continuous token is no better than a shuffled token that belongs to another node, at 43\.9% against 44\.2% on the same 312 questions\. On those questions the discrete role reaches 47\.4% and the compiled description 81\.4%, so a frozen model reads neither form of the embedding\.
The trained yardstick is a graph neural network trained on the questions themselves and scored by the same evaluator on its held\-out split\. It reaches 74\.8% on 381 questions where the frozen readers ofBioGlyphreach 74\.1% pooled, with Qwen3\-8B above it at 81\.9% and Llama\-3\.1\-8B below it at 66\.4%\. The two trade question families\. The trained model is ahead on edge removal at 85\.7% against 67\.9% and on counterfactual reachability at 96\.4% against 83\.0%\. The frozen reader is ahead on degree comparison at 96\.4% against 80\.4% and on the compositional family at 42\.7% against 22\.9%\. The retrieval ceiling, the share of answerable questions a perfect reader of the retrieved region could answer, is 82\.6% pooled, 88\.0% on the dense networks and 72\.6% on the sparse ones\.
Figure S10:Learned roles and the trained yardstick\.aThe four learned\-role encoders against the compiler on every main network; the marker is the encoder, the diamond the compiled description with its 95% interval\.bA continuous graph token against a shuffled one and against the discrete bucket\-id encoder, on the same 312 questions\.cA graph neural network trained on the benchmark questions against a frozen reader of the compiled description, on the trained model’s own held\-out split, by question family and by model\. Table[S8](https://arxiv.org/html/2609.03229#S6.T8)gives the five fusion arms against the description alone\.
## S10Retrieval
Retrieval is a shared component rather than part of the comparison, but it sets what any reader can reach\. Figure[S11](https://arxiv.org/html/2609.03229#S10.F11)gives the study: what retrieval alone leaves answerable, the two settings the paper pools for every arm, and the oracle\-free alternative\. The wide\-window run on the human interactomes is in Fig\.[S4](https://arxiv.org/html/2609.03229#S3.F4)e\. The wider setting raises the ceiling, at 86\.6% against 78\.7%, and raises the overflow of every long arm at once\. The compiled description is the only rendering whose accuracy is essentially the same at both settings, because it is the only one that never leaves the budget\.
The oracle\-free alternative retrieves by personalized PageRank from the named targets\. It takes no route from the question family, admits no counterfactual path, and never runs the verification loop, so no part of it touches privileged information \(Methods\)\. The study covers the three families this retrieval suits, with 1,471 questions read by the two 8B models, and the ordering of the arms holds\.BioGlyphreaches 49\.9% against 33\.0% for the raw table, which loses 27\.2% of its prompts\. On the questions both arms fit, the two stand at 44\.9% against 39\.1% with an interval that excludes zero\. Adjacency text catches up under this retrieval, level with the description under the controlled view\. On the two families shared with the main benchmark at the same cap, the description scores 52\.5% under oracle\-free retrieval and 51\.5% under the standard rings\. The main results do not rest on the verification loop\.
Figure S11:Retrieval\.aWhat retrieval alone leaves answerable, by setting, by network and by question family; the dotted line is the pooled share\.bThe two retrieval settings for each competing arm: system accuracy and the share of prompts that never reached the model\.cThe oracle\-free alternative, which retrieves by personalized PageRank from the named targets and never consults the graded predicate; the dotted line is the no\-network floor on the same questions\. Table[S8](https://arxiv.org/html/2609.03229#S6.T8)gives the paired contrasts\.
## S11One region in every representation
Figure[S12](https://arxiv.org/html/2609.03229#S11.F12)shows one retrieved region of the yeast interactome, chosen deterministically\. The region is rendered as the edge list, as adjacency text, as the raw\-measurement table, as the three rungs of the compiled description and as the opaque control, each with its token count\. A reader can see what the model reads under each arm and why the raw table is the longest of them\.
Figure S12:One region, every representation\.The same retrieved region of STRING\-Yeast rendered as the edge list, adjacency text, the raw\-measurement table, the three rungs of the compiled description and the opaque control, each with its prompt length in tokens; the median length of each arm on this network is given for scale\. The role\-name rung and the opaque control differ only in the vocabulary their role lines use, which the two narrow panels of the top row show side by side\.
## S12The candidate\-knockout screen: a biology\-graded pick, pre\-registered
Every other experiment in this paper grades a model against the graph oracle\. The candidate\-knockout screen closes the remaining quadrant: the model reads compiled structure and makes a pick that biology grades\. One question shows one retrieved module and five of its proteins as anonymous node identifiers, exactly one of which carries the biological label\. A model therefore cannot answer from memorized gene biology, and chance is 20%\. Three experiments of 300 questions each are analyzed separately and never pooled: SGD essentiality on the yeast interactome, and DepMap selective and common essentiality on Reactome\-FI\. The design, the arms, the reference policies and five predictions were pre\-registered before any model ran \(docs/P3\_CANDIDATE\_SCREEN\.md\); the registered Qwen3\-32B leg never ran, so two 8B models report\. Because the biological signal lives in a protein’s whole\-network position, the structural arms render the compiler’s whole\-network output for the region’s nodes, and a degree\-only arm carries the strongest single deterministic policy into the prompt\.
Table[S11](https://arxiv.org/html/2609.03229#S12.T11)reports every arm beside every policy, and the registered predictions largely hold\. The no\-graph control answers fewer than 42% of its questions and lands below the chance line on all three experiments, so the identifiers do not leak\. The degree\-only arm lands beside the degree policy on all three, at 24\.5% against 25\.0%, 39\.3% against 42\.0% and 46\.2% against 47\.3%, so the models read a stated maximum\. As registered, the description and the raw measurements stay close on both Reactome experiments, and on yeast, where every policy sits within a few percentage units of chance, no arm separates from another reliably\. One registered prediction failed, and we report it as written\. On the common\-essential experiment, adjacency text at 46\.5% leads both structural arms\. There the description at 39\.3% trails the bare degree heuristic at 46\.2%, with an interval that excludes zero: shown the full structural story, the models are distracted from the one quantity that predicts the label\. The registered headline stands\. No arm beats the degree policy by a meaningful margin, so a frozen model turns compiled whole\-network structure into a biology\-graded pick at roughly the level of the best deterministic policy, and not beyond it\.
Table S11:The candidate\-knockout screen: picking the biologically labeled protein out of five anonymous candidates\.One question shows one retrieved module and five candidate node ids, exactly one of which carries the label; the pick is graded by SGD or DepMap, never by the graph oracle\. Model rows pool Qwen3\-8B and Llama\-3\.1\-8B \(Wilson 95% intervals\); policy rows are the deterministic reference picks stored per question at build time\. The three experiments are analyzed separately, as pre\-registered\.Similar Articles
Structural Interpretations of Protein Language Model Representations via Differentiable Graph Partitioning
This paper proposes SoftBlobGIN, a framework that enhances the interpretability of protein language model representations by projecting them onto contact graphs for structure-aware message passing. It demonstrates improved performance on enzyme classification and binding-site detection while providing auditable structural explanations.
Are Large Language Models Suitable for Graph Computation? Progress and Prospects
This survey reviews the use of large language models for graph computation, categorizing them into two paradigms: LLMs as executors and LLMs as planners. It finds LLMs promising for simple tasks but unreliable for large-scale exact computations, and suggests future directions.
Speaking the Language of Science: Toward a General-Purpose Generative Foundation Model for the Natural Sciences
LOGOS is a scientific generative language model that encodes diverse scientific objects and spatial interactions as token sequences, enabling a unified autoregressive framework for tasks across natural sciences. Models at 1B, 3B, and 8B parameters show consistent performance scaling and are released to facilitate research.
Scaling Laws for Hypernetwork-Based Knowledge Injection in Large Language Models
The paper investigates scaling laws for hypernetwork-based knowledge injection into LLMs, finding predictive power law scaling and reliable out-of-distribution generalization, establishing hypernetworks as a scalable alternative to LoRA and full fine-tuning.
Graph-Native Reinforcement Learning Enables Traceable Scientific Hypothesis Generation through Conceptual Recombination
This paper introduces Graph-PRefLexOR, a family of graph-native reasoning models fine-tuned with group relative policy optimization (GRPO) to generate traceable scientific hypotheses through explicit reasoning phases. The method achieves 40-65% improvements over base models in reasoning traceability and demonstrates enhanced semantic diversity and conceptual recombination.