Esophageal carcinoma ESCA
Of 60,498 genes, 2,246 are higher in the tumor than in normal esophagus and in whole blood, and are not made by blood immune cells. 935 are protein-coding, and 924 of those have supporting plasma Cell-free RNA (cfRNA)RNA fragments that circulate in blood plasma outside cells. Most come from blood cells; a small share comes from other tissues, including tumors. evidence; the rest are mostly non-coding genes the plasma sources cannot check.
Samples compared
- Tumor
- 181
- TCGA primary tumors
- Normal tissue
- 666
- 13 TCGA tumor-adjacent, 653 GTEx esophagus
- Whole blood
- 337
- GTEx samples from 328 donors
- Survival
- 181
- patients, 85 events (Survival endpointsCurated TCGA outcomes (Liu et al. 2018): OS is overall survival, PFI the progression-free interval and DSS disease-specific survival. LB2 uses whichever has the most events in each cancer type.)
Results in brief
- 2,246
- candidate genes: 924 with plasma cfRNA evidence, 17 protein-coding without, 1,305 non-coding
- 221
- genes in the tissue panel, Held-out AUCHow well an elastic-net model separates tumor from normal tissue, measured with nested cross-validation on samples the model did not train on. It describes tissue, not the accuracy of a blood test. > 0.999
- 0
- genes associated with progression-free interval (FDR < 0.05, 150 tested)
Candidates
Filter by plasma evidence or discrimination, find a gene, and select any point or row for its full record. The table holds all 2,246 candidates.
How the candidates were selected
60,498 genes tested; 7,216 higher in tumor than in normal esophagus; 5,183 also higher than in whole blood; 2,246 remain after removing genes expressed in blood immune cells; 924 of these have plasma cfRNA evidence, 17 protein-coding candidates have none, and 1,305 are non-coding genes the plasma sources cannot list.
Each comparison uses a two-sided Wilcoxon rank-sum test. A gene passes when its q-value is below 0.05, its log2 fold change is at least 1, and the lower bound of the 95% confidence interval for its AUC is at least 0.70.
Genes made by any of 18 sorted blood immune cell types above 1 nTPM are then removed, because blood cells supply most plasma RNA. Genes missing from that reference, mostly non-coding, are kept.
Plasma evidence only reorders the list. A candidate not yet seen in plasma stays in, because every dataset misses genes. Read the full methods.
Tissue classifier panel
- STC2+0.28
- IGF2BP3+0.26
- COL11A1+0.24
- CHRNA1+0.24
- COL10A1+0.24
- ZIC5+0.21
- DMP1+0.21
- PCDHGA8+0.21
- INHBA+0.20
- CST1+0.19
- CCNE1+0.19
- CBX2+0.19
- MRGBP+0.18
- COL5A3+0.18
- HOXC6+0.17
- ITPR3−0.17
- PCGF2+0.17
- HOXC11+0.17
- HOXC13+0.16
- STEAP1−0.16
Show the other 201 genes
- ADAMTS12+0.16
- CHST1+0.15
- PGF+0.15
- C17orf53+0.15
- KIAA1549L+0.15
- ATP6V0D2+0.14
- STRIP2+0.14
- RECQL4+0.14
- SLC5A12+0.14
- WISP1+0.13
- WNT5A+0.13
- GABRD+0.13
- DMC1+0.13
- CILP2+0.13
- MFI2+0.13
- EPHB2−0.13
- FOXS1+0.13
- DNMT3B+0.13
- MAP6D1+0.13
- LRRC1−0.12
- FOXD1+0.12
- HOXC9+0.12
- C17orf96+0.12
- HOXC10+0.12
- HOXC8+0.12
- OSR2−0.12
- ADGRF4+0.11
- C11orf84+0.11
- NOTCH3+0.11
- BMP8A+0.11
- SIX1+0.11
- CHSY3−0.10
- HAP1+0.10
- COL5A2+0.10
- BAMBI+0.10
- E2F1+0.10
- GPRIN1+0.10
- LOXL2+0.10
- TMEM5−0.10
- EPCAM−0.10
- ULBP1+0.10
- MFAP2+0.10
- MET−0.10
- ANLN+0.10
- DMBX1+0.09
- HOXD9+0.09
- COL7A1+0.09
- EXOC3L1−0.09
- DCLK3+0.09
- ZIC2+0.09
- STK31+0.09
- MTFR2+0.09
- LAMC2+0.09
- CCDC73+0.09
- HSF2BP+0.09
- ZRANB3−0.09
- GRIN2D+0.09
- TICRR+0.08
- SH2D5+0.08
- RAPGEF5−0.08
- NOX4+0.08
- KIAA1524+0.08
- LRFN4+0.08
- HOXA13+0.08
- THY1−0.08
- ENPEP+0.08
- NPHP4−0.08
- NAT14+0.08
- MTA1−0.08
- HOXB9+0.07
- BTBD16+0.07
- SNAPC4−0.07
- UTP20−0.07
- SEMA5B−0.07
- TTL−0.07
- ATAD5+0.07
- DPF1+0.07
- FAM64A+0.07
- CD276−0.07
- HOXA11+0.07
- TMEM26+0.07
- COL27A1−0.07
- HGH1−0.07
- COL12A1+0.06
- MTBP+0.06
- LIPG−0.06
- HOXD3−0.06
- TIGD5+0.06
- DSCC1+0.06
- LMLN−0.06
- C19orf48−0.06
- FAP+0.06
- CLDN1+0.06
- NRARP−0.06
- BRIP1+0.06
- FAM86C1+0.05
- POMT2−0.05
- IL17D+0.05
- CDH24+0.05
- HOXD8+0.05
- DBF4+0.05
- MMP1−0.05
- ZNF519+0.05
- CCDC74A+0.05
- CYP26B1+0.05
- TFAP4−0.04
- AGRN−0.04
- TAF5−0.04
- TMEM237+0.04
- SENP3+0.04
- CENPF+0.04
- TMC7−0.04
- FAM72B+0.04
- PLEKHA8−0.04
- APOE+0.04
- POLR1A−0.04
- ZNF598−0.04
- SLC39A4+0.04
- GPR176−0.04
- ECE2+0.04
- LINGO1−0.04
- ARHGAP39+0.04
- MEX3A+0.04
- FAT1−0.04
- POLQ+0.04
- FJX1+0.03
- VARS−0.03
- POLE−0.03
- ZNF462−0.03
- KRT17+0.03
- DNAH3+0.03
- FOXL1−0.03
- P4HA3−0.03
- DHX33−0.03
- RCCD1−0.03
- WDR72+0.03
- GTF2IRD1+0.03
- TEAD4−0.03
- DGKH+0.03
- OSMR−0.03
- BMP1+0.03
- MAPK15+0.03
- SPDYE3+0.03
- ZNF367+0.03
- HCN3−0.03
- ZNF628+0.03
- RTN4RL2−0.03
- DBF4B+0.02
- FAM72D+0.02
- KIF15+0.02
- ARHGEF39+0.02
- GAD1+0.02
- TWIST1+0.02
- COL1A1+0.02
- ATAT1−0.02
- SOX12+0.02
- ARHGAP11A+0.02
- PRR19+0.02
- RHBDF1−0.02
- ADAMTS14−0.02
- CLSPN−0.02
- CLCN2+0.02
- TAF4−0.02
- CASC5−0.02
- TMEM67−0.01
- ZNF707−0.01
- KIF4A+0.01
- PIF1+0.01
- THUMPD2−0.01
- ATP13A3−0.01
- CENPJ+0.01
- PLXNA1−0.01
- IER5L−0.01
- ZNF469−0.01
- STIL+0.01
- C9orf117−0.01
- PLXNA3−0.01
- LAMA3−0.01
- FBXL19−0.01
- FOXD2−0.01
- ITIH6−0.01
- RP11-2C24.9−0.01
- BGN+0.01
- TCF20+0.01
- MMP14+0.01
- TFAP2A+0.01
- CASC10+0.01
- C8orf37−0.01
- TMEM158+0.01
- MCM8+0.01
- ITGA2+0.01
- CDCA2−0.00
- CPOX+0.00
- INPP5E−0.00
- CLDN12−0.00
- QSOX2−0.00
- PSME4+0.00
- ERVMER34-1+0.00
- DUS4L+0.00
- SKA3−0.00
- NAGS+0.00
An elastic-net logistic regression was trained on the 300 top-ranked candidates to tell tumor from normal tissue. Its Held-out AUCHow well an elastic-net model separates tumor from normal tissue, measured with nested cross-validation on samples the model did not train on. It describes tissue, not the accuracy of a blood test. is > 0.999, from nested five-fold cross-validation.
Before training, the TCGA-versus-GTEx offset was removed from the expression matrix with the tumor/normal contrast protected, as in the tissue comparison. That correction is fitted once on all samples, so the held-out AUC is not fully independent of it.
This number describes tumor and normal tissue, where separation is expected to be near perfect. It is not the accuracy of a blood test: that has to be measured in plasma from patients and controls.
Association with survival
None of the 150 top-ranked candidates is associated with progression-free interval at FDR < 0.05 (181 patients, 85 events). A null result here is informative: these genes mark the presence of the cancer, not necessarily its course.
For the 150 top-ranked candidates, a Cox model relates tumor expression, as a continuous value, to progression-free interval in 181 patients (85 events). No high/low cutpoint is searched for, since optimized cutpoints inflate false positives.
A Hazard ratio per SDFrom a Cox model with the gene’s tumor expression as a continuous variable: the change in hazard for each one-standard-deviation increase. Above 1, higher expression goes with a shorter time to the event; below 1, with a longer one. above 1 means higher expression goes with a shorter progression-free interval. These are associations in tissue, not evidence that a blood level predicts outcome.
Caveats for ESCA
- The blood comparison sets TCGA tumors against GTEx blood. Study and biology cannot be separated there, which is why genes made by blood immune cells are removed as well.