Biological Relevance

gene networks reflect tissue and condition specific biology mavatar discovery co expression gene networks successfully capture meaningful biological relationships by quantifying pairwise correlations between gene expression profiles and integration of vast number of datasets and samples, we obtained networks that represent the coordinated activity of genes within and across tissues the resulting structures revealed clear biological organization, consistent with expected tissue specific patterns and functional groupings, supporting the utility of this approach for biological interpretation methodology to evaluate biological relevance and uncover characteristic gene interactions within dina networks, we applied a network comparison framework that identifies both conserved molecular processes and context specific functional processes across the networks this approach yields two complementary outcomes dina network core genes and network specific genes capturing the distinctive functional wiring of each network the analysis revealed both underlying common core biological functions and distinct, context specific functionalities reflecting specialized activity in the biological condition of interest for each dina context specific network, iterative subnetwork comparisons are performed against all the other networks using the deltacon method ( docid\ m3ywqj3arjsy4zaf7wthc ) biologically relevant deltacon null distributions are used as a background ensuring that the similarities and differences reflect true biological organization rather than structural noise (see docid 0zknxsqa6r3cuqdntt6w for more details) this analysis yields two complementary gene lists for each network comparison similarity genes representing conserved patterns and dissimilarity genes representing network specific wiring these lists are subsequently aggregated and filtered to produce two network specific outputs dina network core genes and dina network specific signature genes for biological validation analysis we grouped these genes into louvain communities based on the network of interest and performed over representation analysis (ora) per network community for reactome, go\ bp, go\ mf and go\ cc pathway databases through g\ profiler ( docid\ m3ywqj3arjsy4zaf7wthc ) this approach provides a statistically grounded way to pinpoint functional gene modules that are conserved or rewired between biological contexts, enabling interpretation of tissue specific or disease specific pathways of the networks within the platform, users have access to the global dina network core gene list and all the dina network specific signature gene lists results and conclusions analysis of subnetworks allowed us to identify modules (gene groups) that are either conserved across tissues or uniquely present in one tissue or condition (see docid\ m3ywqj3arjsy4zaf7wthc ) shared subnetworks are likely to correspond to fundamental cellular processes active across multiple biological contexts, whereas tissue or disease specific subnetworks represent specialized or perturbed functions here we present the biological validation results from the following networks docid\ m3ywqj3arjsy4zaf7wthc , docid\ m3ywqj3arjsy4zaf7wthc , docid\ m3ywqj3arjsy4zaf7wthc , docid\ m3ywqj3arjsy4zaf7wthc , docid\ m3ywqj3arjsy4zaf7wthc and docid\ m3ywqj3arjsy4zaf7wthc across tissues, conditions, and disease states, the modules frequently showed enrichment for broad functional categories such as extracellular matrix organization, mitochondrial and respiratory processes, translation and ribosome biogenesis, immune response, rna splicing regulation, and cell cycle or dna replication in the following sections, we focus specifically on modules exhibiting tissue specific functional enrichment, which provide insight into biological processes characteristic of each network pancreas – general based on the enrichment analysis, the pancreas network specific signature genes in module 9 showed strong association with pancreatic physiological functions (table 1) the top enriched terms within this module include protein digestion and absorption related terms, which represent core functional processes of the pancreas in contrast, the similarity tail enrichment (table 2) results highlight subnetworks related to metal ion homeostasis and detoxification, including metallothionein binding metals and response to metal ions, consistent with the pancreas’ function in metal metabolism and stress response regulation co occurrence of metal ion homeostasis within the same module as the pancreatic secretion and digestion terms may suggests that the pancreas specific functional programs are built upon a conserved regulatory backbone of metal ion metabolism and detoxification ( docid\ m3ywqj3arjsy4zaf7wthc ) table 1 top enrichment results for pancreas – general module 9 specific signature genes true 102,278 8662131519274,280 1337868480726 unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type table 2 top enrichment results for pancreas – general module 9 genes, co expressed between multiple networks true 109,275 3741496598639,276 6258503401361 unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type brain – general brain network module 6 exhibited exclusive dissimilarity tail enrichment for neural and glial developmental pathways (table 3), such as axon ensheathment, myelination, oligodendrocyte differentiation, and glial cell differentiation additional terms such as neurogenesis, gliogenesis, and midbrain/substantia nigra development further support a central nervous system–specific module the absence of similarity tail enrichment indicates that these processes are strongly tissue specific and reflect specialized subnetwork organization within the brain network module 12 showed limited similarity tail enrichment (table 4) but extensive dissimilarity tail enrichment for synaptic structure and function, such as vesicle mediated transport in synapse, synaptic/postsynaptic membranes, nmda receptor–associated pathways, neurexin/neuroligin interactions, and exocytosis (table 5) this pattern indicates a brain specific synaptic module whose subnetwork organization is distinct in brain compared with other networks table 3 top enrichment results for brain – general module 6 specific signature genes true 104,277 86848072562356,279 13151927437644 unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type table 4 enrichment results for brain – general module 12 genes, co expressed between multiple networks true 106,276 8707482993197,278 1292517006803 unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type table 5 top enrichment results for brain – general module 12 specific signature genes true 107,276 3718820861678,277 6281179138322 unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type skin – general skin network module 1 showed striking dissimilarity tail enrichment for epidermis and skin development, keratinocyte differentiation, keratinization, and cornified envelope formation, reflecting canonical skin barrier and keratinocyte functionalities (table 6) the corresponding similarity tail was enriched for broad vesicle membrane processes (table 7), consistent with conserved cellular regulation skin cells use endocytosis for basic functions such as taking in nutrients and signaling molecules which are closely connected to epidermis development these findings indicate that this module captures a skin specific functional program built on a shared signaling framework module 2 showed exclusive dissimilarity tail enrichment for melanocyte and pigmentation pathways, such as melanosome, melanosome membrane, melanin biosynthesis/metabolism, and mitf dependent programs (table 8), with no significant terms in the similarity tail we therefore interpret this as a skin specific pigmentation module skin module 16 displayed strong dissimilarity tail enrichment for keratin filament organization, and epidermal development, including intermediate filament cytoskeleton organization, cornified envelope formation, keratinocyte differentiation”, and hair follicle morphogenesis (table 9) the similarity tail contained no significant terms this pattern identifies a skin specific structural module, representing keratinocyte and epidermal differentiation programs characteristic of stratified epithelial tissue table 6 top enrichment results for skin – general module 1 specific signature genes true 115,272 38095238095235,273 61904761904765 unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type table 7 enrichment results for skin – general module 1 genes, co expressed between multiple networks true 120,269 8866213151927,271 1133786848073 unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type table 8 top enrichment results for skin – general module 2 specific signature genes true 121,269 3877551020408,270 6122448979592 unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type table 9 top enrichment results for skin – general module 16 specific signature genes true 124,267 891156462585,269 108843537415 unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type liver – general liver network module 2 showed dissimilarity tail enrichment for a combination characteristic of hepatic biology, such as complement and coagulation cascades, acute phase response, γ carboxylation of protein precursors, and peroxisomal/fatty acid and amino acid metabolism (table 10) the corresponding similarity tail enrichment (table 11) highlighted general metabolic processes, such as small molecule and amino acid catabolism, nicotinamide/nucleotide metabolism, and metal ion response, consistent with conserved cellular functions together, these results support this module as liver specific on the dissimilarity side, built on a conserved metabolic backbone table 10 top enrichment results for liver – general module 2 specific signature genes true 128,265 8956916099773,267 1043083900227 unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type table 11 top enrichment results for liver – general module 2 genes, co expressed between multiple networks true 126,266 89342403628115,268 10657596371885 unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type breast – breast cancer/tissue the breast cancer network exhibited several coherent and biologically interpretable modules the similarity tail of module 4 could be linked to endoplasmic reticulum protein folding, unfolded protein response, rrna processing, and ribosome biogenesis (table 12), marking a shared secretory, biosynthetic core typical of highly secretory epithelial tumors its dissimilarity tail emphasized transcriptional and rna processing regulation, such as major mrna splicing, rna polymerase ii elongation/pause release, and wdr5/nsl histone acetyltransferase complexes (table 13), reflecting cohort specific tuning of epigenetic and transcriptional output together these patterns indicate that breast cancer networks capture conserved cell cycle, biosynthetic, and stress response machinery, while the observed dissimilarity enrichments point to meaningful variation in dna repair, telomere biology, and transcriptional control breast cancer module 26 similarity tail was characterized by metal ion homeostasis functions, such as metallothionein and zinc/copper detoxification processes, a common oxidative stress signature (table 14) its dissimilarity tail highlighted telomere maintenance and telomerase associated repair, indicating heterogeneity in replicative lifespan control (table 15) the similarity tail of a dominant cell cycle/replication module (module 27) captured chromosome segregation, mitotic checkpoints, and dna replication, reflecting a conserved proliferative program across tumors (table 16) the associated dissimilarity signal was enriched for tp53 regulation, checkpoint control, and dna damage repair pathways, such as homologous recombination, brca/rad51, and base excision repair (table 17) this was consistent with variation in replicative stress and dna repair proficiency among cohorts table 12 top enrichment results for breast – breast cancer/tissue module 4 genes, co expressed between multiple networks true 116,271 8820861678004,273 1179138321996 unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type table 13 top enrichment results for breast – breast cancer/tissue module 4 specific signature genes true 120,269 8866213151927,271 1133786848073 unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type table 14 top enrichment results for breast – breast cancer/tissue module 26 genes, co expressed between multiple networks true 123,268 39002267573693,269 60997732426307 unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type table 15 top enrichment results for breast – breast cancer/tissue module 26 specific signature genes true 130,264 89795918367344,266 10204081632656 unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type table 16 top enrichment results for breast – breast cancer/tissue module 8 genes, co expressed between multiple networks true 132,263 9002267573696,265 0997732426304 unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type table 17 top enrichment results for breast – breast cancer/tissue module 8 specific signature genes true 128,265 8956916099773,267 1043083900227 unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type pbmc – general immune related modules appeared across multiple networks, reflecting conserved co expression structures underlying innate and adaptive immunity in pbmcs, this took the form of a stable antigen presentation/mhc ii module 1, which similarity tail showed baseline immune composition of blood mononuclear cells (table 18) while similar immune wiring may appear in other contexts, the pbmc module was distinguished by minimal activation signatures in its dissimilarity tail (table 19), suggesting maintenance of immune identity rather than context specific inflammatory activation in addition to pbmc module 1, two related immune modules were identified the first module captured innate effector and antigen processing functions, including neutrophil degranulation, phagocytosis, vesicle lumen, and lysosomal organization, while the corresponding dissimilarity enrichments highlighted toll like receptor and myd88 dependent signaling the second module represented broad cytokine signaling processes, encompassing tnf, il 17, il 10, and nf κb pathways, together with downstream transcriptional regulators these patterns reflect the mixed cellular composition and variable activation states characteristic of peripheral blood, indicating that the framework recovers coherent immune modules spanning both basal myeloid functions and cytokine mediated responses table 18 top enrichment results for pbmc – general module 1 genes, co expressed between multiple networks true 127,266 3945578231292,267 6054421768708 unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type table 19 enrichment results for pbmc – general module 1 specific signature genes true 118,270 88435374149657,272 11564625850343 unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type 1 1 unhandled content type summary our investigations confirm that mavatar discovery networks represent biologically coherent structures across all comparisons, the modules identified in both similarity and dissimilarity tails align with well established tissue and disease biology shared subnetworks capture conserved cellular and metabolic functions, while unique modules reflect context specific biology, such as immune regulation, tissue specialization, or microenvironmental interactions the consistency of these findings across multiple tissues and orientations demonstrates that the networks faithfully represent the underlying biological systems and capture both common and specialized processes relevant to each context in mavatar discovery, these results have been made accessible for graph construction the “network specific gene lists” represent genes with local connectivity patterns that differ from other networks, defining the functional rewiring unique to each context in contrast, the “dina network core genes” comprise one list of genes with connectivity profiles closely matching those in other networks, representing the shared cellular scaffold common to many systems references d koutra, j t vogelsteiny, c faloutsos, deltacon a principled massive graph similarity function proceedings of the 2013 siam international conference on data mining, sdm 2013 , 1304–4657 (2013) https //doi org/10 1145/2824443 l kolberg, u raudvere, i kuzmin, p adler, j vilo, h peterson, g\ profiler interoperable web service for functional enrichment analysis and gene identifier mapping (2023 update) nucleic acids res 51, w207–w212 (2023) https //pubmed ncbi nlm nih gov/37144459 e y chen, c m tan, y kou, q duan, z wang, g v meirelles, n r clark, a ma’ayan, enrichr interactive and collaborative html5 gene list enrichment analysis tool bmc bioinformatics 14, 128 (2013)