Q (sonnet): In lower grade glioma, are there genes which are overexpressed in any of the molecular subtypes?

▶ list_guides
{}

▶ list_studies
{
 "search": "lower grade glioma"
}

◀ result
{"result":[{"uri":"cbioportal://mutation-frequency-guide","description":"Comprehensive guide for calculating gene mutation frequencies with gene-specific profiling denominators"},{"uri":"cbioportal://clinical-data-guide","description":"Guide for querying clinical data including patient vs sample level considerations"},{"uri":"cbioportal://sample-filtering-guide","description":"Guide for filtering samples and studies in cBioPortal queries"},{"uri":"cbioportal://common-pitfalls","description":"Guide to avoid common mistakes when querying cBioPortal data. If you already know which numbered pitfall applies, fetch just that section via read_guide(\"cbioportal://common-pitfalls#<number>\") (e.g. #16) instead of the full guide"},{"uri":"cbioportal://treatment-guide","description":"Guide for querying treatment/clinical event data including drug agents, timelines, and linking to genomic data"},{"uri":"cbioportal://faq-guide","description":"General cBioPortal FAQ: history, how to cite, data types, reference genome, abbreviations, GISTIC thresholds, API access"},{"uri":"cbioportal://statistical-tests-guide","description":"Statistical test selection guide — decision matrix for choosing Fisher's exact, Wilcoxon, chi-squared, t-test, ANOVA, etc. based on data type and group count"},{"uri":"cbioportal://gene-expression-guide","description":"Gene expression / copy-number / methylation analysis. Covers genetic_alteration_derived, profile_type discovery, and the gene_pair_coexpression view for Spearman correlation between two genes"},{"uri":"cbioportal://external-resources-guide","description":"Guide for finding external linked resources such as imaging, pathology, Minerva, HTAN, or other resource_* table links before declaring data unavailable"},{"uri":"cbioportal://gene-resolution-guide","description":"Guide for resolving ambiguous gene symbols, aliases, gene families, and shorthand such as CD3 before querying expression or alteration data"},{"uri":"cbioportal://study-resolution-guide","description":"Guide for resolving requested studies, avoiding silent substitute cohorts, and redirecting to known external cBioPortal instances when data is not in this deployment"},{"uri":"cbioportal://germline-guide","description":"Guide for querying germline variant data — storage columns, study discovery, query patterns, and somatic vs germline considerations"},{"uri":"cbioportal://study-guide/{study_id}","description":"Dynamic study-specific guide - use get_study_guide(study_id) tool to generate"}]}

◀ result
{"result":[{"cancer_study_identifier":"lgg_tcga","name":"Brain Lower Grade Glioma (TCGA, Firehose Legacy)","type_of_cancer_id":"difg","sample_count":530,"has_guide":false,"url":"https://www.cbioportal.org/study/summary?id=lgg_tcga"},{"cancer_study_identifier":"lgg_tcga_pan_can_atlas_2018","name":"Brain Lower Grade Glioma (TCGA, PanCancer Atlas)","type_of_cancer_id":"difg","sample_count":514,"has_guide":false,"url":"https://www.cbioportal.org/study/summary?id=lgg_tcga_pan_can_atlas_2018"}]}

▶ read_guide
{
 "uri": "cbioportal://gene-expression-guide"
}

▶ read_guide
{
 "uri": "cbioportal://statistical-tests-guide"
}

◀ result
# Gene Expression Analysis Guide

This guide covers continuous-value genomic data: gene **expression**, **copy number** values, **methylation**, and related profile types. Mutation/CNA/SV *frequency* analysis lives in `cbioportal://mutation-frequency-guide`.

## Where this data lives

Continuous per-sample-per-gene values are stored in `genetic_alteration_derived`:

| Column | Description |
|---|---|
| `sample_unique_id` | `<study>_<sample_stable_id>` |
| `cancer_study_identifier` | study scope |
| `hugo_gene_symbol` | gene |
| `profile_type` | which assay/normalization (see below) |
| `alteration_value` | the actual value — stored as Nullable(String); cast with `toFloat64OrNull` |

`alteration_value` is a string because the same column hosts many different value scales. The `''` and `'NA'` sentinels mean "missing"; always filter them out and use `toFloat64OrNull(alteration_value) IS NOT NULL` for downstream math.

## Discovering profile types for a study

Different studies expose different profile types depending on what assays were run and how the data was normalized. Always check what a specific study supports before picking one:

```sql
SELECT DISTINCT profile_type
FROM genetic_alteration_derived
WHERE cancer_study_identifier = 'brca_metabric'
ORDER BY profile_type;
```

Common values across the public portal:

| Family | Profile types |
|---|---|
| mRNA expression | `mrna`, `mrna_median_Zscores`, `mrna_seq_v2_rsem`, `mrna_seq_v2_rsem_Zscores`, `mrna_seq_cpm`, `mrna_seq_fpkm`, `mrna_U133`, `mrna_outliers` |
| Copy number (continuous) | `cna`, `linear_CNA`, `log2CNA`, `cna_consensus`, `cna_rae`, `gistic` |
| Methylation | `methylation_hm27`, `methylation_hm450`, `methylation_epic`, `methylation_promoters_rrbs` |
| miRNA | `mirna`, `mirna_median_Zscores` |
| Protein | `protein_quantification`, `protein_level`, `RPPA` |

**Z-score vs raw choice.** When the user asks "is X correlated with Y", either works for Spearman (rank-based) — Pearson would care. Default to the non-Z-score variant if both exist, and call out which one in the response.

## Canonical recipe — Spearman correlation between two genes

```sql
SELECT *
FROM gene_pair_coexpression(
    study        = 'brca_metabric',
    gene_a       = 'TP53',
    gene_b       = 'MYC',
    profile_type = 'mrna'
);
```

Returns one row: `(gene_a, gene_b, profile_type, spearman_correlation, num_samples)`.

- `spearman_correlation` in [−1, 1]; `NULL` when fewer than 3 valid paired samples.
- Mirrors cbioportal-backend's `ClickhouseCoExpressionMapper.getCoExpressions`, simplified to a pair lookup (the backend computes one ref gene vs ALL other genes for the coexpression page; here the agent asks about a specific pair).

### Verified examples

| Study | gene_a | gene_b | profile_type | spearman | n |
|---|---|---|---|---|---|
| `brca_metabric` | TP53 | MYC | `mrna` | 0.118 | 1980 |
| `brca_metabric` | ESR1 | PGR | `mrna` | 0.487 | 1980 |

ESR1↔PGR is the textbook breast-cancer estrogen-receptor coregula … (2785 more chars)

▶ read_guide
{
 "uri": "cbioportal://clinical-data-guide"
}

◀ result
Statistical Test Selection Guide
================================

Purpose
-------
This guide ensures the correct statistical test is selected before performing any group comparison, matching cBioPortal's own Group Comparison defaults.

HARD RULES — NEVER FABRICATE A STATISTIC
----------------------------------------
ClickHouse cannot run statistical tests. The agent therefore must NEVER produce a derived statistic that is not a literal column value from a SQL result. Specifically:

1. **Never invent a p-value.** Not "p < 0.001", not "p ≈ 0.05", not any p-value. If the user asks "what is the p-value?", the answer is *"I can't compute that — here is the 2x2 contingency table (or group statistics). Run it in cBioPortal's Group Comparison tab, in R with `fisher.test(...)` / `wilcox.test(...)`, or in Python with `scipy.stats.fisher_exact(...)` / `mannwhitneyu(...)`."*
2. **Never claim mutual exclusivity (or co-occurrence) from a contingency table alone.** A 2x2 table is not a test. The shape "altered/not altered × group A/group B" needs Fisher's exact + a defined direction (odds ratio < 1 with significant p). Without that test, the agent presents the table and stops. Descriptive phrasing is also a claim: "largely/mostly mutually exclusive", "rarely co-occur", "tend to co-occur" are forbidden without the test — point to cBioPortal's Mutual Exclusivity tab instead.
3. **Never report a "median" that came from `AVG(...)` or any non-median aggregate.** "Median" and "mean" are different statistics; for skewed clinical distributions (especially survival) they differ substantially. Use ClickHouse's `quantile(0.5)(...)` for actual median, and label arithmetic averages as "mean", never "median".
4. **Never report a hazard ratio, odds ratio, risk ratio, or relative risk** that wasn't computed by an external tool. These require regression / model fitting that ClickHouse does not do.
5. **Never report median overall survival from `AVG(OS_MONTHS)` or even `quantile(0.5)(OS_MONTHS)`.** Median OS requires Kaplan-Meier estimation, which handles censoring (`OS_STATUS = 0:LIVING` means the event hasn't happened yet). Naive medians/means over `OS_MONTHS` ignore censoring and are systematically wrong. The correct handoff: return the raw `(OS_MONTHS, OS_STATUS)` pairs (or descriptive counts: N events, N censored, follow-up range) and tell the user to run KM in R (`survival::survfit`) or Python (`lifelines.KaplanMeierFitter`), or use cBioPortal's Survival comparison.

If asked to do any of (1)–(5), respond with the appropriate handoff template from the "Approved Response Templates" section below — do not produce the number.

Ambiguous-Term Clarification (before any analysis)
--------------------------------------------------
If the user's question uses an ambiguous term, ask before computing anything. Examples:

- **"aggressive"** → could mean shorter overall survival, higher metastasis rate, higher grade/stage, higher TMB, or specific molecular features. Ask which definit … (9984 more chars)

◀ result
# Clinical Data Query Guide

## Overview
Clinical data in cBioPortal is stored at both patient and sample levels. Understanding the distinction is crucial for accurate analysis.

## Data Organization

### Patient-Level vs Sample-Level Data
- **Patient-level**: Demographics, overall survival, disease stage (stored once per patient)
- **Sample-level**: Sample type, sequencing platform, purity (can have multiple per patient)

### Key Tables
- `clinical_patient`: Patient-level clinical attributes
- `clinical_sample`: Sample-level clinical attributes
- `clinical_data_derived`: Pre-joined view combining both levels
- `clinical_attribute_meta`: Metadata about available clinical attributes

## Recommended Approach: Use clinical_data_derived

The `clinical_data_derived` table is pre-joined and optimized for most queries:

```sql
-- Get clinical data for specific attributes
SELECT
    sample_unique_id,
    patient_unique_id,
    attribute_name,
    attribute_value
FROM clinical_data_derived
WHERE
    cancer_study_identifier = 'your_study_id'
    AND attribute_name IN ('CANCER_TYPE', 'SAMPLE_TYPE', 'TMB_NONSYNONYMOUS');
```

## Clinical Attribute Discovery

### Use clinical_attribute_meta for Discovering Available Attributes
**Always start here** to see what clinical attributes are available for a specific study:

```sql
-- Discover available clinical attributes for a study
SELECT
    attr_id,
    description,
    patient_attribute,
    cancer_study_id
FROM clinical_attribute_meta
WHERE
    cancer_study_id = (
        SELECT cancer_study_id
        FROM cancer_study
        WHERE cancer_study_identifier = 'your_study_id'
    )
ORDER BY patient_attribute, attr_id;
```

**Key Fields:**
- **attr_id**: matches attr_id in clinical_sample/clinical_patient tables
- **description**: human-readable description of the attribute
- **patient_attribute**: true = patient attribute, false = sample attribute
- **cancer_study_id**: links to cancer_study table (filter by study)

## Attribute Semantics and Matching

### Case-Insensitive Matching for Attribute Values

Clinical values are free text across studies and may differ only by case. For example, a controlled-looking value such as germline mutation status may appear as `GERMLINE`, `Germline`, or another case variant.

When filtering `clinical_data_derived.attribute_value`, use case-insensitive matching unless you have already profiled the exact values in the target study:

```sql
-- Correct: case-insensitive clinical value filter
SELECT DISTINCT sample_unique_id, patient_unique_id
FROM clinical_data_derived
WHERE cancer_study_identifier = 'your_study_id'
  AND attribute_name = 'MUTATION_STATUS'
  AND upper(attribute_value) = 'GERMLINE';
```

Do not write `attribute_value = 'GERMLINE'` without first checking all distinct values for that attribute in the study.

### Query the Requested Attribute, Not a Proxy

Do not infer one clinical attribute from a related subtype or marker. Query the actual requested attribute when it … (17073 more chars)

▶ clickhouse_run_select_query
{
 "query": "SELECT DISTINCT attribute_name FROM clinical_data_derived WHERE cancer_study_identifier = 'lgg_tcga_pan_can_atlas_2018' AND (upper(attribute_name) LIKE '%SUBTYPE%' OR upper(attribute_name) LIKE '%CODEL%' OR upper(attribute_name) LIKE '%IDH%')"
}

▶ clickhouse_run_select_query
{
 "query": "SELECT DISTINCT profile_type FROM genetic_alteration_derived WHERE cancer_study_identifier = 'lgg_tcga_pan_can_atlas_2018' ORDER BY profile_type"
}

◀ result
{"rows":[{"attribute_name":"SUBTYPE"}]}

◀ result
{"rows":[{"profile_type":"gistic"},{"profile_type":"log2CNA"},{"profile_type":"rna_seq_v2_mrna"},{"profile_type":"rna_seq_v2_mrna_median_Zscores"},{"profile_type":"rna_seq_v2_mrna_median_all_sample_Zscores"},{"profile_type":"rppa"},{"profile_type":"rppa_Zscores"}]}

▶ clickhouse_run_select_query
{
 "query": "SELECT attribute_value, COUNT(DISTINCT patient_unique_id) as n
FROM clinical_data_derived
WHERE cancer_study_identifier = 'lgg_tcga_pan_can_atlas_2018' AND attribute_name = 'SUBTYPE'
GROUP BY attribute_value
ORDER BY n DESC"
}

◀ result
{"rows":[{"attribute_value":"LGG_IDHmut-non-codel","n":248},{"attribute_value":"LGG_IDHmut-codel","n":167},{"attribute_value":"LGG_IDHwt","n":92},{"n":7}]}

▶ clickhouse_run_select_query
{
 "query": "
WITH subtype AS (
    SELECT patient_unique_id, attribute_value AS subtype
    FROM clinical_data_derived
    WHERE cancer_study_identifier = 'lgg_tcga_pan_can_atlas_2018' AND attribute_name = 'SUBTYPE'
      AND attribute_value != ''
),
expr AS (
    SELECT g.sample_unique_id, g.hugo_gene_symbol, toFloat64OrNull(g.alteration_value) AS z
    FROM genetic_alteration_derived g
    WHERE g.cancer_study_identifier = 'lgg_tcga_pan_can_atlas_2018'
      AND g.profile_type = 'rna_seq_v2_mrna_median_Zscores'
      AND g.alteration_value NOT IN ('', 'NA')
      AND toFloat64OrNull(g.alteration_value) IS NOT NULL
),
joined AS (
    SELECT e.hugo_gene_symbol, s.subtype, e.z
    FROM expr e
    INNER JOIN subtype s ON s.patient_unique_id = substring(e.sample_unique_id, 1, length(e.sample_unique_id)-3)
)
SELECT hugo_gene_symbol,
       subtype,
       count() AS n,
       round(avg(z),3) AS mean_z
FROM joined
GROUP BY hugo_gene_symbol, subtype
HAVING n >= 30
ORDER BY mean_z DESC
LIMIT 5
"
}

◀ result
{"rows":[{"hugo_gene_symbol":"SEPTIN14","subtype":"LGG_IDHwt","n":92,"mean_z":235.919},{"hugo_gene_symbol":"KLK11","subtype":"LGG_IDHwt","n":92,"mean_z":31.377},{"hugo_gene_symbol":"SEC61G","subtype":"LGG_IDHwt","n":92,"mean_z":23.448},{"hugo_gene_symbol":"CYP27B1","subtype":"LGG_IDHwt","n":92,"mean_z":20.864},{"hugo_gene_symbol":"ASZ1","subtype":"LGG_IDHwt","n":92,"mean_z":18.326}]}

▶ clickhouse_run_select_query
{
 "query": "
WITH subtype AS (
    SELECT patient_unique_id, attribute_value AS subtype
    FROM clinical_data_derived
    WHERE cancer_study_identifier = 'lgg_tcga_pan_can_atlas_2018' AND attribute_name = 'SUBTYPE'
      AND attribute_value != ''
),
expr AS (
    SELECT g.sample_unique_id, g.hugo_gene_symbol, toFloat64OrNull(g.alteration_value) AS z
    FROM genetic_alteration_derived g
    WHERE g.cancer_study_identifier = 'lgg_tcga_pan_can_atlas_2018'
      AND g.profile_type = 'rna_seq_v2_mrna_median_Zscores'
      AND g.alteration_value NOT IN ('', 'NA')
      AND toFloat64OrNull(g.alteration_value) IS NOT NULL
),
joined AS (
    SELECT e.hugo_gene_symbol, s.subtype, e.z
    FROM expr e
    INNER JOIN subtype s ON s.patient_unique_id = substring(e.sample_unique_id, 1, length(e.sample_unique_id)-3)
)
SELECT hugo_gene_symbol,
       subtype,
       count() AS n,
       round(avg(z),3) AS mean_z,
       round(quantile(0.5)(z),3) AS median_z,
       round(min(z),2) AS min_z,
       round(max(z),2) AS max_z
FROM joined
GROUP BY hugo_gene_symbol, subtype
HAVING n >= 30 AND median_z > 1.5
ORDER BY median_z DESC
LIMIT 15
"
}

◀ result
{"rows":[{"hugo_gene_symbol":"MEOX2","subtype":"LGG_IDHwt","n":92,"mean_z":6.642,"median_z":4.498,"min_z":-0.22,"max_z":28.26},{"hugo_gene_symbol":"TMEM248","subtype":"LGG_IDHwt","n":92,"mean_z":4.365,"median_z":3.978,"min_z":-1.67,"max_z":40.55},{"hugo_gene_symbol":"FKBP9P1","subtype":"LGG_IDHwt","n":92,"mean_z":6.44,"median_z":3.534,"min_z":-0.52,"max_z":56.55},{"hugo_gene_symbol":"MRM2","subtype":"LGG_IDHwt","n":92,"mean_z":3.281,"median_z":3.495,"min_z":-2.92,"max_z":12.3},{"hugo_gene_symbol":"FKBP9","subtype":"LGG_IDHwt","n":92,"mean_z":4.06,"median_z":3.316,"min_z":0.36,"max_z":21.25},{"hugo_gene_symbol":"PSMC2","subtype":"LGG_IDHwt","n":92,"mean_z":3.131,"median_z":3.156,"min_z":-1.19,"max_z":11.17},{"hugo_gene_symbol":"DENND2A","subtype":"LGG_IDHwt","n":92,"mean_z":3.359,"median_z":3.153,"min_z":-1.15,"max_z":10.87},{"hugo_gene_symbol":"GSTK1","subtype":"LGG_IDHwt","n":92,"mean_z":3.151,"median_z":3.138,"min_z":-1.17,"max_z":9.73},{"hugo_gene_symbol":"PDIA4","subtype":"LGG_IDHwt","n":92,"mean_z":3.334,"median_z":3.085,"min_z":-1.74,"max_z":12.28},{"hugo_gene_symbol":"RPS2P32","subtype":"LGG_IDHwt","n":92,"mean_z":3.757,"median_z":3.068,"min_z":-0.85,"max_z":15.89},{"hugo_gene_symbol":"FZD1","subtype":"LGG_IDHwt","n":92,"mean_z":3.16,"median_z":3.025,"min_z":-1.29,"max_z":31.82},{"hugo_gene_symbol":"DNAAF5","subtype":"LGG_IDHwt","n":92,"mean_z":2.817,"median_z":2.924,"min_z":-2.27,"max_z":12.73},{"hugo_gene_symbol":"HOXA1","subtype":"LGG_IDHwt","n":92,"mean_z":5.152,"median_z":2.882,"min_z":-0.5,"max_z":70.34},{"hugo_gene_symbol":"PAXIP1-AS2","subtype":"LGG_IDHwt","n":92,"mean_z":3.24,"median_z":2.869,"min_z":-1.24,"max_z":13.31},{"hugo_gene_symbol":"TMEM209","subtype":"LGG_IDHwt","n":92,"mean_z":2.94,"median_z":2.858,"min_z":-1.54,"max_z":13.62}]}

▶ clickhouse_run_select_query
{
 "query": "
WITH subtype AS (
    SELECT patient_unique_id, attribute_value AS subtype
    FROM clinical_data_derived
    WHERE cancer_study_identifier = 'lgg_tcga_pan_can_atlas_2018' AND attribute_name = 'SUBTYPE'
      AND attribute_value != ''
),
expr AS (
    SELECT g.sample_unique_id, g.hugo_gene_symbol, toFloat64OrNull(g.alteration_value) AS z
    FROM genetic_alteration_derived g
    WHERE g.cancer_study_identifier = 'lgg_tcga_pan_can_atlas_2018'
      AND g.profile_type = 'rna_seq_v2_mrna_median_Zscores'
      AND g.alteration_value NOT IN ('', 'NA')
      AND toFloat64OrNull(g.alteration_value) IS NOT NULL
),
joined AS (
    SELECT e.hugo_gene_symbol, s.subtype, e.z
    FROM expr e
    INNER JOIN subtype s ON s.patient_unique_id = substring(e.sample_unique_id, 1, length(e.sample_unique_id)-3)
),
agg AS (
    SELECT hugo_gene_symbol, subtype, count() AS n, avg(z) AS mean_z, quantile(0.5)(z) AS median_z
    FROM joined
    GROUP BY hugo_gene_symbol, subtype
    HAVING n >= 30
)
SELECT subtype, hugo_gene_symbol, round(mean_z,3) AS mean_z, round(median_z,3) AS median_z
FROM agg
WHERE subtype = 'LGG_IDHmut-codel'
ORDER BY median_z DESC
LIMIT 8
"
}

▶ clickhouse_run_select_query
{
 "query": "
WITH subtype AS (
    SELECT patient_unique_id, attribute_value AS subtype
    FROM clinical_data_derived
    WHERE cancer_study_identifier = 'lgg_tcga_pan_can_atlas_2018' AND attribute_name = 'SUBTYPE'
      AND attribute_value != ''
),
expr AS (
    SELECT g.sample_unique_id, g.hugo_gene_symbol, toFloat64OrNull(g.alteration_value) AS z
    FROM genetic_alteration_derived g
    WHERE g.cancer_study_identifier = 'lgg_tcga_pan_can_atlas_2018'
      AND g.profile_type = 'rna_seq_v2_mrna_median_Zscores'
      AND g.alteration_value NOT IN ('', 'NA')
      AND toFloat64OrNull(g.alteration_value) IS NOT NULL
),
joined AS (
    SELECT e.hugo_gene_symbol, s.subtype, e.z
    FROM expr e
    INNER JOIN subtype s ON s.patient_unique_id = substring(e.sample_unique_id, 1, length(e.sample_unique_id)-3)
),
agg AS (
    SELECT hugo_gene_symbol, subtype, count() AS n, avg(z) AS mean_z, quantile(0.5)(z) AS median_z
    FROM joined
    GROUP BY hugo_gene_symbol, subtype
    HAVING n >= 30
)
SELECT subtype, hugo_gene_symbol, round(mean_z,3) AS mean_z, round(median_z,3) AS median_z
FROM agg
WHERE subtype = 'LGG_IDHmut-non-codel'
ORDER BY median_z DESC
LIMIT 8
"
}

◀ result
{"rows":[{"subtype":"LGG_IDHmut-codel","hugo_gene_symbol":"TNNT1","mean_z":2.567,"median_z":1.923},{"subtype":"LGG_IDHmut-codel","hugo_gene_symbol":"ATCAY","mean_z":1.27,"median_z":1.11},{"subtype":"LGG_IDHmut-codel","hugo_gene_symbol":"CACNG2","mean_z":1.097,"median_z":1.084},{"subtype":"LGG_IDHmut-codel","hugo_gene_symbol":"SYCE2","mean_z":1.149,"median_z":0.968},{"subtype":"LGG_IDHmut-codel","hugo_gene_symbol":"CHGB","mean_z":1.042,"median_z":0.964},{"subtype":"LGG_IDHmut-codel","hugo_gene_symbol":"FEM1A","mean_z":1.124,"median_z":0.943},{"subtype":"LGG_IDHmut-codel","hugo_gene_symbol":"NALF1","mean_z":0.93,"median_z":0.936},{"subtype":"LGG_IDHmut-codel","hugo_gene_symbol":"LRRTM4","mean_z":0.964,"median_z":0.903}]}

◀ result
{"rows":[{"subtype":"LGG_IDHmut-non-codel","hugo_gene_symbol":"RIPK2","mean_z":0.811,"median_z":0.704},{"subtype":"LGG_IDHmut-non-codel","hugo_gene_symbol":"ZNF212","mean_z":0.868,"median_z":0.682},{"subtype":"LGG_IDHmut-non-codel","hugo_gene_symbol":"STMP1","mean_z":0.827,"median_z":0.657},{"subtype":"LGG_IDHmut-non-codel","hugo_gene_symbol":"LUC7L2","mean_z":0.823,"median_z":0.636},{"subtype":"LGG_IDHmut-non-codel","hugo_gene_symbol":"CHCHD3","mean_z":0.808,"median_z":0.624},{"subtype":"LGG_IDHmut-non-codel","hugo_gene_symbol":"CHST9","mean_z":0.678,"median_z":0.61},{"subtype":"LGG_IDHmut-non-codel","hugo_gene_symbol":"NUDT5","mean_z":0.883,"median_z":0.608},{"subtype":"LGG_IDHmut-non-codel","hugo_gene_symbol":"PTCD1","mean_z":0.852,"median_z":0.597}]}

▶ resolve_and_route
{
 "studyIds": [
  "lgg_tcga_pan_can_atlas_2018"
 ]
}

◀ result
{"success":true,"message":"Found 1 study. Use the metadata to call the appropriate navigation tool(s).","data":{"totalCount":1,"studiesWithMetadata":[{"studyId":"lgg_tcga_pan_can_atlas_2018","name":"Brain Lower Grade Glioma (TCGA, PanCancer Atlas)","sampleCount":514,"studyViewUrl":"https://www.cbioportal.org/study?id=lgg_tcga_pan_can_atlas_2018","metadata":{"clinicalAttributeIds":["AGE","AJCC_PATHOLOGIC_TUMOR_STAGE","AJCC_STAGING_EDITION","ANEUPLOIDY_SCORE","BUFFA_HYPOXIA_SCORE","CANCER_TYPE","CANCER_TYPE_ACRONYM","CANCER_TYPE_DETAILED","DAYS_LAST_FOLLOWUP","DAYS_TO_BIRTH","DAYS_TO_INITIAL_PATHOLOGIC_DIAGNOSIS","DFS_MONTHS","DFS_STATUS","DSS_MONTHS","DSS_STATUS","ETHNICITY","FORM_COMPLETION_DATE","FRACTION_GENOME_ALTERED","GENETIC_ANCESTRY_LABEL","GRADE","HISTORY_NEOADJUVANT_TRTYN","ICD_10","ICD_O_3_HISTOLOGY","ICD_O_3_SITE","INFORMED_CONSENT_VERIFIED","IN_PANCANPATHWAYS_FREEZE","MSI_SCORE_MANTIS","MSI_SENSOR_SCORE","MUTATION_COUNT","NEW_TUMOR_EVENT_AFTER_INITIAL_TREATMENT","ONCOTREE_CODE","OS_MONTHS","OS_STATUS","OTHER_PATIENT_ID","PATH_M_STAGE","PATH_N_STAGE","PATH_T_STAGE","PERSON_NEOPLASM_CANCER_STATUS","PFS_MONTHS","PFS_STATUS","PRIMARY_LYMPH_NODE_PRESENTATION_ASSESSMENT","PRIOR_DX","RACE","RADIATION_THERAPY","RAGNUM_HYPOXIA_SCORE","SAMPLE_COUNT","SAMPLE_TYPE","SEX","SOMATIC_STATUS","SUBTYPE","TBL_SCORE","TISSUE_PROSPECTIVE_COLLECTION_INDICATOR","TISSUE_RETROSPECTIVE_COLLECTION_INDICATOR","TISSUE_SOURCE_SITE","TISSUE_SOURCE_SITE_CODE","TMB_NONSYNONYMOUS","TUMOR_TISSUE_SITE","TUMOR_TYPE","WEIGHT","WINTER_HYPOXIA_SCORE"],"molecularProfileIds":["lgg_tcga_pan_can_atlas_2018_gistic","lgg_tcga_pan_can_atlas_2018_log2CNA","lgg_tcga_pan_can_atlas_2018_mutations","lgg_tcga_pan_can_atlas_2018_rna_seq_v2_mrna","lgg_tcga_pan_can_atlas_2018_rna_seq_v2_mrna_median_Zscores","lgg_tcga_pan_can_atlas_2018_rna_seq_v2_mrna_median_all_sample_Zscores","lgg_tcga_pan_can_atlas_2018_rppa","lgg_tcga_pan_can_atlas_2018_rppa_Zscores","lgg_tcga_pan_can_atlas_2018_structural_variants"],"genericAssayProfiles":["lgg_tcga_pan_can_atlas_2018_armlevel_cna","lgg_tcga_pan_can_atlas_2018_genetic_ancestry","lgg_tcga_pan_can_atlas_2018_methylation_hm27_hm450_merge","lgg_tcga_pan_can_atlas_2018_methylation_hm450"],"heatmapProfileIds":["lgg_tcga_pan_can_atlas_2018_rna_seq_v2_mrna_median_Zscores","lgg_tcga_pan_can_atlas_2018_rna_seq_v2_mrna_median_all_sample_Zscores","lgg_tcga_pan_can_atlas_2018_rppa_Zscores"],"availableComparisonTabs":["overlap","clinical","survival","alterations","mutations","mrna","protein","generic_assay_armlevel_cna","generic_assay_genetic_ancestry","generic_assay_methylation"]}}]}}

▶ get_studyviewfilter_options
{
 "studyId": "lgg_tcga_pan_can_atlas_2018",
 "attributeIds": [
  "SUBTYPE"
 ]
}

◀ result
{"success":true,"message":"Successfully retrieved StudyView filter options","data":{"studyId":"lgg_tcga_pan_can_atlas_2018","attributes":[{"attributeId":"SUBTYPE","displayName":"Subtype","description":"Subtype","datatype":"STRING","values":["NA","LGG_IDHwt","LGG_IDHmut-non-codel","LGG_IDHmut-codel"]}]}}

▶ navigate_to_group_comparison
{
 "studyIds": [
  "lgg_tcga_pan_can_atlas_2018"
 ],
 "groups": [
  {
   "name": "LGG_IDHwt",
   "studyViewFilter": {
    "clinicalDataFilters": [
     {
      "attributeId": "SUBTYPE",
      "values": [
       {
        "value": "LGG_IDHwt"
       }
      ]
     }
    ]
   }
  },
  {
   "name": "LGG_IDHmut-codel",
   "studyViewFilter": {
    "clinicalDataFilters": [
     {
      "attributeId": "SUBTYPE",
      "values": [
       {
        "value": "LGG_IDHmut-codel"
       }
      ]
     }
    ]
   }
  },
  {
   "name": "LGG_IDHmut-non-codel",
   "studyViewFilter": {
    "clinicalDataFilters": [
     {
      "attributeId": "SUBTYPE",
      "values": [
       {
        "value": "LGG_IDHmut-non-codel"
       }
      ]
     }
    ]
   }
  }
 ],
 "tab": "mrna"
}

◀ result
{"success":true,"message":"Navigating to https://www.cbioportal.org/comparison/mrna?comparisonId=6ab5c836c2115c492d884e03","url":"https://www.cbioportal.org/comparison/mrna?comparisonId=6ab5c836c2115c492d884e03","data":{"description":"Group comparison (3 custom groups)","studies":["lgg_tcga_pan_can_atlas_2018"],"totalGroups":3,"groups":[{"name":"LGG_IDHwt","sampleCount":92},{"name":"LGG_IDHmut-codel","sampleCount":167},{"name":"LGG_IDHmut-non-codel","sampleCount":248}],"studyViewUrl":"https://www.cbioportal.org/study?id=lgg_tcga_pan_can_atlas_2018","groupUrls":[{"groupName":"LGG_IDHwt","url":"https://www.cbioportal.org/study?id=lgg_tcga_pan_can_atlas_2018#filterJson=%7B%22clinicalDataFilters%22%3A%5B%7B%22attributeId%22%3A%22SUBTYPE%22%2C%22values%22%3A%5B%7B%22value%22%3A%22LGG_IDHwt%22%7D%5D%7D%5D%7D"},{"groupName":"LGG_IDHmut-codel","url":"https://www.cbioportal.org/study?id=lgg_tcga_pan_can_atlas_2018#filterJson=%7B%22clinicalDataFilters%22%3A%5B%7B%22attributeId%22%3A%22SUBTYPE%22%2C%22values%22%3A%5B%7B%22value%22%3A%22LGG_IDHmut-codel%22%7D%5D%7D%5D%7D"},{"groupName":"LGG_IDHmut-non-codel","url":"https://www.cbioportal.org/study?id=lgg_tcga_pan_can_atlas_2018#filterJson=%7B%22clinicalDataFilters%22%3A%5B%7B%22attributeId%22%3A%22SUBTYPE%22%2C%22values%22%3A%5B%7B%22value%22%3A%22LGG_IDHmut-non-codel%22%7D%5D%7D%5D%7D"}],"pageDescription":"mRNA expression for all genes compared across LGG_IDHwt (92) vs LGG_IDHmut-codel (167) vs LGG_IDHmut-non-codel (248)."}}

═ answer (success)
## Summary

**Study:** Brain Lower Grade Glioma (TCGA, PanCancer Atlas) — `lgg_tcga_pan_can_atlas_2018`, 514 samples. The `SUBTYPE` clinical attribute defines three molecular subtypes:

| Subtype | N patients |
|---|---|
| LGG_IDHmut-non-codel | 248 |
| LGG_IDHmut-codel | 167 |
| LGG_IDHwt | 92 |

Using `rna_seq_v2_mrna_median_Zscores` (z-scored relative to diploid samples in the cohort), I screened per-gene mean/median z-score within each subtype (genes with ≥30 profiled samples in the subtype). This is a descriptive screen only — no formal test was run.

**Yes — clear group-wide overexpression signals exist, concentrated in the IDH-wildtype subtype:**

| Subtype | Top genes by median z-score | Median z (top gene) |
|---|---|---|
| **LGG_IDHwt** | MEOX2, TMEM248, FKBP9P1, MRM2, FKBP9, PSMC2, DENND2A, GSTK1, PDIA4, FZD1, HOXA1 | up to ~4.5 |
| LGG_IDHmut-codel | TNNT1, ATCAY, CACNG2, SYCE2, CHGB | up to ~1.9 |
| LGG_IDHmut-non-codel | RIPK2, ZNF212, STMP1, LUC7L2, CHCHD3 | up to ~0.7 |

The IDHwt subtype shows the strongest, most robust elevation (positive median z-score across essentially the whole 92-sample group, not driven by outliers), consistent with its known transcriptional divergence from IDH-mutant tumors. The other two subtypes show only modest shifts by comparison.

**Statistical caveat:** This is a 3-group comparison of continuous expression data — the appropriate test is **one-way ANOVA** (or Kruskal-Wallis if normality is questionable) per gene, ideally with Benjamini-Hochberg FDR correction across all genes tested. ClickHouse cannot run ANOVA; the ranking above is based on raw z-score magnitude only, not significance. Run the actual differential expression analysis with FDR-corrected p-values in the link below.

[View mRNA Expression Comparison: LGG_IDHwt vs LGG_IDHmut-codel vs LGG_IDHmut-non-codel](https://www.cbioportal.org/comparison/mrna?comparisonId=6ab5c836c2115c492d884e03)

mRNA expression for all genes compared across LGG_IDHwt (92) vs LGG_IDHmut-codel (167) vs LGG_IDHmut-non-codel (248) — this tab runs the proper per-gene test with volcano plot and q-values.