# The Signal in the Extreme: A Systematic Outlier Framework Identifies Discrete Immunometabolic Subtypes in Human and Cellular Models

**Authors:** Julio Jesús Garcia-Coste, Karla Aidee Aguayo-Cerón, Judith Espinosa-Raya, Alexis Alejandro García-Rivero, Carina López-Leyva, Rocío Alejandra Gutiérrez-Rojas, Cruz Vargas-De-León, Rodrigo Romero-Nava

PMC · DOI: 10.3390/medsci14010128 · 2026-03-09

## TL;DR

This study shows that analyzing outliers in metabolic-inflammation data reveals hidden subtypes that traditional methods miss, offering new insights into disease mechanisms.

## Contribution

A novel framework for systematically analyzing outliers to identify discrete immunometabolic subtypes in complex diseases.

## Key findings

- An outlier subgroup in the clinical cohort showed hyperactivation of Th1/Th17 pathways and hypertriglyceridemia.
- Outlier samples in the cellular model exhibited IL-6 overproduction and IL-10 suppression.
- Multivariate analysis confirmed distinct spatial segregation of the identified immunometabolic profiles.

## Abstract

Background: Conventional omics analysis often treats outliers as noise, yet they may harbor critical biological insights. Objetive: This study proposes a paradigm shift: actively investigating outliers to discover biologically relevant subtypes within metabolic–inflammatory syndromes. Methods: We applied a comprehensive analytical framework for outlier detection based on a multi-algorithm consensus (IQR, MAD, Isolation Forest) to a clinical cohort of diabetic neuropathy (n = 93) and an in vitro 3T3-L1 adipocyte model (n = 39). The identified outliers were characterized using robust PCA, co-expression networks, unsupervised clustering, and Random Forest predictive modeling. Results: In the clinical cohort, an outlier subgroup (47.3%) exhibited an extreme immune–metabolic phenotype characterized by hyperactivation of Th1/Th17 pathways (elevated T-bet and IL-17; p < 0.001), hypertriglyceridemia, and network reconfiguration (TGFβ and STAT4 hubs). In the cellular model, outlier samples (12.8%) showed autonomous pro-inflammatory behavior characterized by IL-6 overproduction (p = 0.002) and IL-10 suppression. Conclusions: Multivariate analysis confirmed spatial segregation of these profiles. Systematic outlier investigation revealed discrete pathophysiological subtypes invisible to mean-focused analyses, demonstrating that extreme values encapsulate potent biological signals. This framework offers a generalizable approach for uncovering clinical heterogeneity and identifying therapeutic targets in complex diseases.

## Linked entities

- **Proteins:** TBX21 (T-box transcription factor 21), IL17A (interleukin 17A), TGFB1 (transforming growth factor beta 1), STAT4 (signal transducer and activator of transcription 4), IL6 (interleukin 6), IL10 (interleukin 10)
- **Diseases:** diabetic neuropathy (MONDO:0006626)

## Full-text entities

- **Genes:** ERCC1 [NCBI Gene 101699727], TNFRSF1A (TNF receptor superfamily member 1A) [NCBI Gene 7132] {aka CD120a, FPF, TBP1, TNF-R, TNF-R-I, TNF-R55}, STAT3 (signal transducer and activator of transcription 3) [NCBI Gene 6774] {aka ADMIO, ADMIO1, APRF, HIES}, IL23A (interleukin 23 subunit alpha) [NCBI Gene 51561] {aka IL-23, IL-23A, IL23P19, P19, SGRF}, TNFRSF1B (TNF receptor superfamily member 1B) [NCBI Gene 7133] {aka CD120b, TBPII, TNF-R-II, TNF-R75, TNFBR, TNFR1B}, BLK (BLK proto-oncogene, Src family tyrosine kinase) [NCBI Gene 640] {aka MODY11}, TGFB1 (transforming growth factor beta 1) [NCBI Gene 7040] {aka CAEND1, CED, DPD1, IBDIMDE, LAP, TGF-beta1}, IL12B (interleukin 12B) [NCBI Gene 3593] {aka CLMF, CLMF2, IL-12B, IMD28, IMD29, NKSF}, STAT6 (signal transducer and activator of transcription 6) [NCBI Gene 6778] {aka D12S1644, HIES6, IL-4-STAT, STAT6B, STAT6C}, CASP3 (caspase 3) [NCBI Gene 836] {aka CPP32, CPP32B, SCA-1}, SOCS1 (suppressor of cytokine signaling 1) [NCBI Gene 8651] {aka AISIMD, CIS1, CISH1, JAB, SOCS-1, SSI-1}, STAT5A (signal transducer and activator of transcription 5A) [NCBI Gene 6776] {aka MGF, STAT5}, ADIPOQ (adiponectin, C1Q and collagen domain containing) [NCBI Gene 9370] {aka ACDC, ACRP30, ADIPQTL1, ADPN, APM-1, APM1}, IL6 (interleukin 6) [NCBI Gene 3569] {aka BSF-2, BSF2, CDF, HGF, HSF, IFN-beta-2}, KRT16 (keratin 16) [NCBI Gene 3868] {aka CK16, FNEPPK, K16, K1CP, KRT16A, NEPPK}, KRT6B (keratin 6B) [NCBI Gene 3854] {aka CK-6B, CK6B, K6B, KRTL1, PC2, PC4}, IL1B (interleukin 1 beta) [NCBI Gene 3553] {aka IL-1, IL1-BETA, IL1F2, IL1beta}, IL10 (interleukin 10) [NCBI Gene 3586] {aka CSIF, GVHDS, IL-10, IL10A, TGIF}, STAT4 (signal transducer and activator of transcription 4) [NCBI Gene 6775] {aka DPMC, SLEB11}, TBX21 (T-box transcription factor 21) [NCBI Gene 30009] {aka IMD88, T-PET, T-bet, TBET, TBLYM}, SLCO1B1 (solute carrier organic anion transporter family member 1B1) [NCBI Gene 10599] {aka HBLRR, LST-1, OATP-C, OATP1B1, OATP2, OATPC}, SOCS3 (suppressor of cytokine signaling 3) [NCBI Gene 9021] {aka ATOD4, CIS3, Cish3, SOCS-3, SSI-3, SSI3}, IL17A (interleukin 17A) [NCBI Gene 3605] {aka CTLA-8, CTLA8, IL-17, IL-17A, IL17, ILA17}, NLRP3 (NLR family pyrin domain containing 3) [NCBI Gene 114548] {aka AGTAVPRL, AII, AVP, C1orf7, CIAS1, CLR1.1}, Lep (leptin) [NCBI Gene 16846] {aka ob, obese}, TNF (tumor necrosis factor) [NCBI Gene 7124] {aka DIF, IMD127, TNF-alpha, TNFA, TNFSF2, TNLG1F}, CISH (cytokine inducible SH2 containing protein) [NCBI Gene 1154] {aka BACTS2, CIS, CIS-1, G18, SOCS}, CD40 (CD40 molecule) [NCBI Gene 958] {aka Bp50, CDW40, TNFRSF5, p50}
- **Diseases:** toxicity (MESH:D064420), Diabetes Neuropathy (MESH:D003929), inflammatory syndromes (MESH:D018746), neuroinflammation (MESH:D000090862), diabetes (MESH:D003920), neuropathy (MESH:D009422), desquamation (MESH:D017490), infertility (MESH:D007246), axonal damage (MESH:D001480), neuropatia diabetica (MESH:C565484), dyslipidemia (MESH:D050171), dysregulations (MESH:D021081), injury to (MESH:D014947), Hantavirus Pulmonary Syndrome (MESH:D018804), neuropathic pain (MESH:D009437), insulin resistance (MESH:D007333), infections (MESH:D007239), CH (MESH:C566784), oncogenesis (MESH:D063646), type 2 diabetes (MESH:D003924), hypertriglyceridemia (MESH:D015228), obesity (MESH:D009765), fever (MESH:D005334), immune dysregulation (OMIM:614878), metabolic syndrome (MESH:D024821), systemic disease (MESH:D034721), axonal degeneration (MESH:D009410), acute and fatal respiratory failure (MESH:D012131), immune-mediated disorders (MESH:C567355), edema (MESH:D004487), Kawasaki disease (MESH:D009080), endothelial dysfunction (MESH:D014652), rhabdomyolysis (MESH:D012206), morbid obesity (MESH:D009767), Inflammation (MESH:D007249), erythema (MESH:D004890), metabolic (MESH:D008659), autoimmune disorders (MESH:D001327), fibrosis (MESH:D005355), microvascular complications (OMIM:603933), Cancer (MESH:D009369), hyperphagia (MESH:D006963)
- **Chemicals:** TG (MESH:D014280), Lipid (MESH:D008055), cholesterol (MESH:D002784), Gini (-), Lipofectamine 2000 (MESH:C086724), free fatty acids (MESH:D005230), glycine (MESH:D005998), urea (MESH:D014508), glucose (MESH:D005947), uric acid (MESH:D014527), creatinine (MESH:D003404)
- **Species:** Homo sapiens (human, species) [taxon 9606], Mus musculus (house mouse, species) [taxon 10090], Heterocephalus glaber (naked mole rat, species) [taxon 10181]
- **Mutations:** A-GGG-24
- **Cell lines:** 3T3-L1 — Mus musculus (Mouse), Spontaneously immortalized cell line (CVCL_0123), PC3 — Homo sapiens (Human), Prostate carcinoma, Cancer cell line (CVCL_0035)

## Figures

16 figures with captions in the complete paper: https://tomesphere.com/paper/PMC13027428/full.md

---
Source: https://tomesphere.com/paper/PMC13027428