Pericyte-derived cancer-associated fibroblasts correlate with poor survival and are enriched after chemoradiotherapy in glioblastoma.
The 5 matches
- [1] § Materials and methods › Bulk RNA-seq analysis and CAF quantification ↔ gse190504_cellanneal_deconvolution.py, lines 200–290 · score 0.98 · bulk_max, bulk_min, disp_min, IDH1 mut, IDH1 wt, CAF proportion
- [2] § Materials and methods › Bulk RNA-seq analysis and CAF quantification ↔ gse190504_CAF_rank_score.py, lines 306–437 · score 0.91 · CAF rank score, IDH1 mut, IDH1 wt, BH adjusted, Wilcoxon, CAF signature
- [3] § Materials and methods › Single cell analysis ↔ functions.R, lines 616–653 · score 0.72 · decontX, contamination, copy, Seurat, single cell, neighbor
- [4] § Materials and methods › Single cell analysis ↔ load_data.R, lines 114–147 · score 0.53 · algorithm, GSE103224, GSM3828672, GSE135045, Harmony, neighbor
- [5] § Results › CAF transcriptional signature reveals association with poor prognosis and therapy-related expansion ↔ gse190504_cellanneal_deconvolution.py, lines 200–290 · score 0.52 · CAF proportions, untreated, Spearman, histology, IDH1, bulk
Paper
Loaded from Europe PMC by your browser, not stored by OSCR: doi.org · Europe PMC
The paper is loaded when this pane is shown.
The authors' code
Python · 294 lines · 11 KB · no license · 2 matches
- from __future__ import annotations
- import re
- import sys
- import zipfile
- import xml.etree.ElementTree as ET
- from pathlib import Path
- import numpy as np
- import pandas as pd
- from scipy.stats import mannwhitneyu, spearmanr, pearsonr
- # ---------------------------------------------------------------------------
- # paths
- ROOT = Path(__file__).resolve().parent
- BASE = ROOT / "GSE190504"
- EXPR_XLSX = BASE / "GSE190504_Processed_Data_Spreadsheet_Glioma_Study.xlsx"
- META_TXT = BASE / "GSE190504_series_matrix.txt"
- REF_CSV = ROOT / "bulk_ref" / "reference_celltype_mean.csv"
- OUT_DIR = BASE / "glio_deconv_bulkref_cellanneal"
- OUT_DIR.mkdir(parents=True, exist_ok=True)
- OUT_PROP = OUT_DIR / "GSE190504_glio_cellanneal_proportions.csv"
- OUT_STATS = OUT_DIR / "GSE190504_glio_cellanneal_stats.txt"
- CELLANNEAL_REPO = ROOT / "cellanneal_repo"
- sys.path.insert(0, str(CELLANNEAL_REPO))
- from scipy.optimize import dual_annealing as _scipy_dual_annealing
- from cellanneal import general as _ca_general
- _ca_general.dual_annealing = _scipy_dual_annealing
- from cellanneal.pipelines import run_cellanneal
- # cellanneal hyper-parameter
- DISP_MIN = 0.5
- BULK_MIN = 1e-8
- BULK_MAX = 1.0
- MAXITER = 500
- # ---------------------------------------------------------------------------
- # self-contained XLSX reader (avoids the openpyxl dependency)
- # ---------------------------------------------------------------------------
- _NS = {"s": "http://schemas.openxmlformats.org/spreadsheetml/2006/main"}
- def _col_letters_to_index(letters: str) -> int:
- idx = 0
- for ch in letters:
- idx = idx * 26 + (ord(ch.upper()) - ord("A") + 1)
- return idx - 1
- def read_xlsx_sheet(xlsx_path: Path, sheet_index: int = 0) -> list[list[str | None]]:
- with zipfile.ZipFile(xlsx_path) as zf:
- shared: list[str] = []
- if "xl/sharedStrings.xml" in zf.namelist():
- tree = ET.parse(zf.open("xl/sharedStrings.xml"))
- for si in tree.getroot().findall("s:si", _NS):
- pieces = [t.text or "" for t in si.iter("{%s}t" % _NS["s"])]
- shared.append("".join(pieces))
- sheet_names = sorted(
- n for n in zf.namelist() if n.startswith("xl/worksheets/sheet")
- )
- tree = ET.parse(zf.open(sheet_names[sheet_index]))
- rows = tree.getroot().find("s:sheetData", _NS)
- out: list[list[str | None]] = []
- max_cols = 0
- for row in rows.findall("s:row", _NS):
- cells: dict[int, str | None] = {}
- for c in row.findall("s:c", _NS):
- col_letters = re.match(r"[A-Z]+", c.attrib["r"]).group(0)
- col_idx = _col_letters_to_index(col_letters)
- t = c.attrib.get("t", "n")
- v_node = c.find("s:v", _NS)
- if t == "s":
- cells[col_idx] = shared[int(v_node.text)] if v_node is not None else None
- elif t == "inlineStr":
- is_node = c.find("s:is", _NS)
- txt = "".join((tn.text or "") for tn in (is_node.iter("{%s}t" % _NS["s"]) if is_node is not None else []))
- cells[col_idx] = txt
- else:
- cells[col_idx] = v_node.text if v_node is not None else None
- if cells:
- max_cols = max(max_cols, max(cells) + 1)
- out.append([cells.get(i) for i in range(max_cols)])
- for r in out:
- if len(r) < max_cols:
- r.extend([None] * (max_cols - len(r)))
- return out
- # ---------------------------------------------------------------------------
- # helpers
- # ---------------------------------------------------------------------------
- def fmt_p(p: float) -> str:
- if p is None or (isinstance(p, float) and (np.isnan(p) or np.isinf(p))):
- return "NA"
- if p < 1e-4:
- return f"{p:.4e}"
- return f"{p:.4f}"
- def sanitize_ct(name: str) -> str:
- return re.sub(r"[^A-Za-z0-9]+", "_", str(name))
- def clean_chr(values) -> list[str]:
- out = []
- for v in values:
- s = "" if v is None else str(v)
- s = s.strip()
- s = re.sub(r'^"|"$', "", s)
- out.append(s.strip())
- return out
- def parse_series_metadata(path: Path) -> pd.DataFrame:
- lines = path.read_text(errors="replace").splitlines()
- def get_row(prefix: str) -> list[str]:
- for ln in lines:
- if ln.startswith(prefix):
- return clean_chr(ln.split("\t")[1:])
- return []
- def get_rows(prefix: str) -> list[list[str]]:
- return [clean_chr(ln.split("\t")[1:]) for ln in lines if ln.startswith(prefix)]
- titles = get_row("!Sample_title")
- chars = get_rows("!Sample_characteristics_ch1")
- if len(chars) < 5:
- raise RuntimeError("Expected >=5 !Sample_characteristics_ch1 lines.")
- histology = [re.sub(r"^histology:\s*", "", v) for v in chars[2]]
- genotype = [re.sub(r"^genotype:\s*", "", v) for v in chars[3]]
- treatment = [re.sub(r"^treatment:\s*", "", v) for v in chars[4]]
- meta = pd.DataFrame({
- "sample_id": titles,
- "histology": histology,
- "genotype": genotype,
- "treatment": treatment,
- })
- meta["idh1"] = np.where(meta["genotype"].str.contains("IDH1-mut", regex=False), "IDH1-mut", "IDH1-wt")
- meta["treatment_simple"] = np.where(meta["treatment"] == "Untreated", "Untreated", "Treated")
- return meta
- def load_expression(xlsx_path: Path) -> pd.DataFrame:
- """Return gene x sample DataFrame, gene symbols as index, linear scale."""
- sheet = read_xlsx_sheet(xlsx_path, sheet_index=0)
- if not sheet:
- raise RuntimeError("XLSX appears empty")
- ensg = [s if s is not None else "" for s in sheet[0][2:]]
- sym = [(s.upper() if s is not None else "") for s in sheet[1][2:]]
- for k in range(len(sym)):
- if not sym[k]:
- sym[k] = ensg[k].upper()
- data_rows = sheet[5:]
- sample_title = [r[1] if len(r) > 1 else None for r in data_rows]
- sample_title = [s if s is not None else "" for s in sample_title]
- values = np.full((len(data_rows), len(sym)), np.nan, dtype=float)
- for i, r in enumerate(data_rows):
- for j in range(len(sym)):
- v = r[2 + j] if 2 + j < len(r) else None
- if v is None or v == "":
- continue
- try:
- values[i, j] = float(v)
- except ValueError:
- continue
- values = np.nan_to_num(values, nan=0.0)
- expr = pd.DataFrame(values.T, index=pd.Index(sym, name="symbol"), columns=sample_title)
- expr.index = expr.index.astype(str).str.upper().str.replace(r"\..*$", "", regex=True)
- expr = expr.loc[(expr.index != "") & (~expr.index.isna())]
- if expr.index.duplicated().any():
- expr = expr.groupby(level=0, sort=False).sum()
- return expr
- def load_reference(ref_csv: Path) -> pd.DataFrame:
- ref = pd.read_csv(ref_csv, index_col=0)
- ref.index = ref.index.astype(str).str.upper().str.replace(r"\..*$", "", regex=True)
- ref = ref.loc[(ref.index != "") & (~ref.index.isna())]
- ref = ref[~ref.index.duplicated(keep="first")]
- ref.columns = [sanitize_ct(c) for c in ref.columns]
- ref = ref.apply(pd.to_numeric, errors="coerce").fillna(0.0)
- return ref
- # ---------------------------------------------------------------------------
- # main
- # ---------------------------------------------------------------------------
- def main() -> None:
- print(f"[info] reading expression: {EXPR_XLSX}")
- expr = load_expression(EXPR_XLSX)
- print(f"[info] expression matrix: {expr.shape[0]} genes x {expr.shape[1]} samples")
- print(f"[info] reading metadata: {META_TXT}")
- meta = parse_series_metadata(META_TXT)
- glio_samples = [s for s in expr.columns if s in set(meta.loc[meta["histology"] == "Glio", "sample_id"])]
- bulk = expr.loc[:, glio_samples].copy()
- print(f"[info] Glio samples: {bulk.shape[1]}")
- print(f"[info] reading reference: {REF_CSV}")
- ref = load_reference(REF_CSV)
- print(f"[info] reference: {ref.shape[0]} genes x {ref.shape[1]} cell types")
- common = bulk.index.intersection(ref.index)
- print(f"[info] common genes (raw): {len(common)}")
- bulk2 = bulk.loc[common].copy()
- ref2 = ref.loc[common].copy()
- keep = (bulk2.sum(axis=1) > 0) & (ref2.sum(axis=1) > 0)
- bulk2 = bulk2.loc[keep]
- ref2 = ref2.loc[keep]
- print(f"[info] common genes after non-zero filter: {bulk2.shape[0]}")
- if bulk2.shape[0] < 1000:
- raise RuntimeError("Too few common genes after filtering; aborting.")
- # composition normalisation per sample / per cell type, matching the
- # GSE222515 cellanneal_bulkref pipeline.
- bulk2 = bulk2.div(bulk2.sum(axis=0), axis=1).fillna(0.0)
- ref2 = ref2.div(ref2.sum(axis=0), axis=1).fillna(0.0)
- print(f"[info] running cellanneal (maxiter={MAXITER}) ...")
- mix = run_cellanneal(
- celltype_df=ref2,
- bulk_df=bulk2,
- disp_min=DISP_MIN,
- bulk_min=BULK_MIN,
- bulk_max=BULK_MAX,
- maxiter=MAXITER,
- )
- res = mix.copy().reset_index().rename(columns={"index": "sample_id"})
- # standardise CAF column name
- caf_col = next((c for c in res.columns if c.lower() == "cafs"), None)
- if caf_col is None:
- raise RuntimeError("CAFs column not found in cellanneal output")
- if caf_col != "CAFs":
- res = res.rename(columns={caf_col: "CAFs"})
- res = res.merge(meta[["sample_id", "idh1", "treatment_simple", "histology"]], on="sample_id", how="left")
- res.to_csv(OUT_PROP, index=False)
- print(f"[saved] {OUT_PROP}")
- # ---- statistics
- caf_wt = res.loc[res["idh1"] == "IDH1-wt", "CAFs"].to_numpy()
- caf_mut = res.loc[res["idh1"] == "IDH1-mut", "CAFs"].to_numpy()
- p_idh = mannwhitneyu(caf_wt, caf_mut, alternative="two-sided").pvalue if len(caf_wt) and len(caf_mut) else np.nan
- caf_unt = res.loc[res["treatment_simple"] == "Untreated", "CAFs"].to_numpy()
- caf_trt = res.loc[res["treatment_simple"] == "Treated", "CAFs"].to_numpy()
- p_trt = mannwhitneyu(caf_unt, caf_trt, alternative="two-sided").pvalue if len(caf_unt) and len(caf_trt) else np.nan
- # report cellanneal fit QC if present
- qc_cols = [c for c in res.columns if c.lower() in {"rho_spearman", "rho_pearson"}]
- qc_summary = ""
- if qc_cols:
- rho_s = res["rho_Spearman"].median() if "rho_Spearman" in res.columns else np.nan
- rho_p = res["rho_Pearson"].median() if "rho_Pearson" in res.columns else np.nan
- qc_summary = f"Median rho_Spearman: {rho_s:.4f} | Median rho_Pearson: {rho_p:.4f}"
- lines = [
- f"Run: {pd.Timestamp.now()}",
- f"Glio samples: {len(res)}",
- f"Common genes used: {len(common)} (after nonzero filter: {bulk2.shape[0]})",
- f"Reference cell types: {', '.join(ref2.columns)}",
- f"Median CAF proportion: {res['CAFs'].median():.6f}",
- f"Nonzero CAF count: {int((res['CAFs'] > 0).sum())}",
- f"IDH1 wt vs mut p (MWU): {fmt_p(p_idh)} (n_wt={len(caf_wt)}, n_mut={len(caf_mut)})",
- f"Untreated vs Treated p (MWU): {fmt_p(p_trt)} (n_unt={len(caf_unt)}, n_trt={len(caf_trt)})",
- ]
- if qc_summary:
- lines.append(qc_summary)
- OUT_STATS.write_text("\n".join(lines) + "\n")
- print(f"[saved] {OUT_STATS}")
- print("\n".join(lines))
- if __name__ == "__main__":
- main()
gse190504_cellanneal_deconvolution.py at commit d0ffbb4, no license · at the source
Overview
Abstract
The role of cancer-associated fibroblasts (CAFs) in glioblastoma remains unclear, as their existence in the brain tumor microenvironment is still debated, given that the normal brain parenchyma is devoid of fibroblasts. It is unclear whether cells described as CAFs represent a distinct stromal population or a transcriptional state of perivascular cells such as pericytes. The aim of this study was to determine the identity, origin, and functional relevance of CAFs in glioblastoma. We analyzed 54 single-cell RNA sequencing datasets together with 88 bulk RNA sequencing samples. We identified a continuous transcriptional spectrum linking endothelial cells, pericytes, and CAFs, supporting pericytes as the most likely source of CAFs in glioblastoma. We further derived and validated robust CAF- and pericyte-specific gene signatures, enabling clear separation of these populations across cohorts. Reproducible CAF-associated ligand–receptor interactions were enriched in angiogenesis and immune modulation pathways. In bulk RNA-seq data, both CAF signature scoring and deconvolution consistently demonstrated increased CAF abundance in IDH-wildtype gliomas and further enrichment after chemoradiotherapy, while selective CYP1B1 expression in CAFs suggested a potential association with therapy-induced tumor adaptation. Overall, CAFs represent a distinct, pericyte-related stromal population in glioblastoma with conserved transcriptional and signaling programs. High CAF signature scores were associated with poorer overall and progression-free survival and were enriched in IDH-wildtype and post-chemoradiotherapy gliomas, suggesting a role for CAFs in therapy-associated remodeling of the tumor microenvironment in aggressive disease.
Reproduced under the paper's license (CC BY), from the paper cited above.
Repository
Its files are read in the Code ↔ Paper reader above, with 5 matches between paragraphs and lines of code.
neuropromotion/CAFs-in-glioblastoma-microenvironment
d0ffbb4fd540b418b830368d456905d97cfc059d, 11 September 2026Availability: 1 check, the latest on 26 September 2026: the link answers
- 26 September 2026: the link answers
12 files
- CNA.R, R, 163 lines
- Diffusion_maps.R, R, 343 lines
- ESCAPE.R, R, 265 lines
- GAM.R, R, 285 lines
- MASTER_SCRIPT.R, R, 8 lines
- cellchat.R, R, 457 lines
- for_deconvolution.Rmd, R, 119 lines
- functions.R, R, 1,155 lines, 1 match
- gse190504_CAF_rank_score
.py , Python, 441 lines, 1 match - gse190504_cellanneal_dec
onvolution.py , Python, 294 lines, 2 matches - load_data.R, R, 199 lines, 1 match
- README.md, Text, 50 lines
The paper's code and data availability statement is in the Data section.
Tracing map
Proposed by the machine: these links were found in the paper and verified at the source, without human review. The map will receive a Zenodo DOI once one of the paper's authors has validated it with their ORCID.
What the map holds:
- 1 repository of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
- 11 scripts, each with its path and the digest of its content;
- 5 matches between paragraphs of the paper and lines of the code (method lexical-v1);
- neither the text of the paper nor the code itself.
Its JSON (tracing-map.json) is deposited on Zenodo with its DOI once the map is validated.
Data
Datasets cited
- geo:GSE173278, at NCBI GEO; found in the text, “Single cell analysis”
Data Availability
All code used for data processing, analysis, figure generation, and 3D diffusion map visualization is publicly available at: https://
Reproduced under the paper's license (CC BY), from the paper cited above.
Versions
The history of this record: each version stored by the harvester or made by a correction of its authors or of the maintainers of its code, and what changed in its facts. The texts of the paper (its abstract, its availability statements) are not part of it; versions that changed only those are not listed.
Version 1, 27 September 2026: the first record
Recorded: type, language, journal, volume, issue, pages, dates, 2 authors, 8 MeSH terms, 1 funder, 50 references.
Cite
This paper
Ismailov, A., & Poptsova, M. (2026). Pericyte-derived cancer-associated fibroblasts correlate with poor survival and are enriched after chemoradiotherapy in glioblastoma. PloS one, 21(9), e0355902. https://
BibTeX
@article{ismailov2026per
author = {Ismailov, Aly and Poptsova, Maria},
title = {{Pericyte-derived cancer-associated fibroblasts correlate with poor survival and are enriched after chemoradiotherapy in glioblastoma}},
journal = {PloS one},
year = {2026},
month = sep,
volume = {21},
number = {9},
pages = {e0355902},
publisher = {PLOS},
issn = {1932-6203},
doi = {10.1371/
url = {https://
pmid = {42721145},
pmcid = {PMC13561374}
}
RIS
TY - JOUR
AU - Ismailov, Aly
AU - Poptsova, Maria
TI - Pericyte-derived cancer-associated fibroblasts correlate with poor survival and are enriched after chemoradiotherapy in glioblastoma
T2 - PloS one
J2 - PLoS One
PY - 2026
DA - 2026/
VL - 21
IS - 9
SP - e0355902
SN - 1932-6203
PB - PLOS
DO - 10.1371/
UR - https://
LA - en
ER -
CSL-JSON
{
"id": "10.1371/
"type": "article-journal",
"title": "Pericyte-derived cancer-associated fibroblasts correlate with poor survival and are enriched after chemoradiotherapy in glioblastoma",
"container-title": "PloS one",
"author": [
{
"family": "Ismailov",
"given": "Aly"
},
{
"family": "Poptsova",
"given": "Maria"
}
],
"container-title-short":
"volume": "21",
"issue": "9",
"page": "e0355902",
"DOI": "10.1371/
"PMID": "42721145",
"PMCID": "PMC13561374",
"ISSN": "1932-6203",
"publisher": "PLOS",
"URL": "https://
"language": "en",
"issued": {
"date-parts": [
[
2026,
9,
10
]
]
}
}
The tracing map gets a citation of its own once an author has validated it and it has a DOI.
Similar papers
The papers with a page that share the most with this one: the tools found in their code, their categories, datasets, cited references and authors, the rarest counting most.
- [1] doi:10.1016/j.xcrm.2026.102766 [code]
- A longitudinal single-cell and spatial multiomic atlas of pediatric high-grade glioma.Journal: Cell reports. MedicineIn common: Harmony, SingleCellExperiment, circlize, 9 other tools, genetics / omics, other condition, cellular / molecular, 3 references
- [2] doi:10.1016/j.xcrm.2026.102651 [code]
- Integrative CSF profiling identifies disease-specific immune responses in leptomeningeal disease.Journal: Cell reports. MedicineIn common: Harmony, SingleCellExperiment, circlize, 9 other tools, genetics / omics, other condition, cellular / molecular
- [3] doi:10.1002/imt2.70163 [code]
- Spatial multi-omics unveils sphingolipid metabolic reprogramming within the retinal pathological niche.Journal: iMetaIn common: Harmony, SingleCellExperiment, circlize, 9 other tools, genetics / omics, cellular / molecular
- [4] doi:10.1038/s41586-026-10214-2 [code]
- Multidimensional profiling of heterogeneity in supratentorial ependymomas.Journal: NatureIn common: Harmony, SingleCellExperiment, circlize, 8 other tools, genetics / omics, other condition, 1 reference
- [5] doi:10.1038/s41586-026-10612-6 [code]
- Acquired genetic and cell-state changes in IDH-mutant glioma progression.Journal: NatureIn common: Harmony, SingleCellExperiment, ComplexHeatmap, 6 other tools, other condition, cellular / molecular, 2 references
- [6] doi:10.1038/s41593-026-02300-5 [code]
- Integrated single-cell and spatial transcriptomic profiling in ALS uncovers peripheral-to-central immune infiltration and reprogramming.Journal: Nature neuroscienceIn common: Harmony, circlize, ComplexHeatmap, 8 other tools, genetics / omics, other condition, cellular / molecular
- [7] doi:10.1038/s41514-026-00391-9 [code]
- Region-specific transcriptional signatures of brain aging in the absence of neuropathology at the single-cell level.Journal: npj agingIn common: Harmony, circlize, ComplexHeatmap, 8 other tools, genetics / omics, cellular / molecular
- [8] doi:10.1038/s41467-026-71595-6 [code]
- A single-cell and spatial atlas of early human olfactory development.Journal: Nature communicationsIn common: Harmony, SingleCellExperiment, circlize, 7 other tools, genetics / omics
- [9] doi:10.1073/pnas.2523130123 [code]
- FABP7 controls radial glial scaffold stability during human cortical development.Journal: Proceedings of the National Academy of Sciences of the United States of AmericaIn common: Harmony, circlize, ComplexHeatmap, 8 other tools
- [10] doi:10.1186/s13073-026-01704-z [code]
- Gene expression profiling enables refined parcellation of cortical layers in the heterogeneous human cerebral cortex.Journal: Genome medicineIn common: Harmony, SingleCellExperiment, ComplexHeatmap, 7 other tools, genetics / omics, cellular / molecular
Contribute
The authors of this paper can claim it, correct its record and validate its tracing map, and the maintainers of its code (its owner, or a public member of its organization) correct what it says of their repository; anyone signed in can ask for its removal. Every request goes to OSCR's own machine, which answers it; your account page follows them.
Sign in with ORCID to claim this paper as one of its authors, correct its record or validate its tracing map: when the paper's metadata lists your ORCID iD, you are recognized at once. Maintainers of its code: sign in with GitHub, then claim the repository on your account page.
Claim this paper
Correct its record
Say what each link of this record is, remove the ones that are not the paper's, add the ones that are missing. The correction becomes a new version of the record, in its Versions section.
Validate its tracing map
You validate the map as this page shows it: 1 repository of the authors' code, each at its verified commit and with its license, 11 scripts, and 5 matches between paragraphs and code (see the Code and Map sections). It then receives a DOI on Zenodo, with you (your ORCID iD) and OSCR as its creators; the code itself is not deposited.
The map's fingerprint: sha256:2ac383601fd06c30…
Add the badge to its README
The badge links the code to this page. Copy one of these into the README of the paper's code: only you decide where it goes, and nothing is changed for you.
Markdown
[, paste the snippet at the top, then “Commit changes…” and, to review it first, “Create a new branch and start a pull request”. You open the pull request; OSCR asks for no permission.
Request its removal
To ask OSCR to remove this record, the copies of its authors' scripts or its tracing map, use the removal request page: signed in, you say who you are, what to remove and why, then review and confirm the request. Published rules decide every request (how).
Discussion, reproductions, activity
Discussion: questions and error reports about this paper and its code, from signed-in readers and its authors. It opens with sign-in.
Reproductions: reports from readers who ran the authors' code: what they reproduced, with which environment, commit and data. It opens with sign-in.
Activity: what happens around this paper: new versions of its record, its map's validation, discussions and reproductions. It opens with sign-in.
