Recent advances in the genomics of glioblastoma (GBM) led to the introduction of molecular neuropathology but failed to translate into treatment improvement. This is largely attributed to the genetic and phenotypic heterogeneity of GBM, which are considered the major obstacle to GBM therapy. Here, we use advanced human GBM-like organoid (LEGO: Laboratory Engineered Glioblastoma-like Organoid) models and provide an unprecedented comprehensive characterization of LEGO models using single-cell transcriptome, DNA methylome, metabolome, lipidome, proteome, and phospho-proteome analysis. We discovered that genetic heterogeneity dictates functional heterogeneity across molecular layers and demonstrates that NF1 mutation drives mesenchymal signature. Most importantly, we found that glycerol lipid reprogramming is a hallmark of GBM, and several targets and drugs were discovered along this line. We also provide a genotype-based drug reference map using LEGO-based drug screen. This study provides new human GBM models and a research path toward effective GBM therapy.