Tumor bulk transcriptomes for endometrioid (430 samples) and serous (112 samples) endometrial adenocarcinoma were downloaded from TCGA. Cellular signal analysis was then applied
61 (link) to identify the major transcriptional programs used by tumor cells based on our single-cell endometrial atlas. This method fits the raw bulk messenger RNA counts to a weighted linear combination of transcriptomic signals derived from reference single-cell data. To limit the effect of the cell cycle, we only included cells in the G1 phase and excluded a proliferative
SOX9 cluster. Proportions of epithelial-derived signals in the bulk samples were computed as the fraction of samples for which signals (exposures) derived from an epithelial cell cluster exceed the intercept term of the model. Clinical data associated with the samples that exhibited exposure for
SOX9+LGR5+ or
SOX9+LGR5− above the intercept value were further investigated; we used the Kruskal–Wallis test to confirm that
SOX9+LGR5+ signature contribution (exposure) in tumors differs in different cancer stages (as defined in TCGA clinical data), and used the Wilcoxon test and
t-test to assess the significance of the increase in exposure noted in later stages of cancer, as well as several clinical data (Supplementary Table
3).
Garcia-Alonso L., Handfield L.F., Roberts K., Nikolakopoulou K., Fernando R.C., Gardner L., Woodhams B., Arutyunyan A., Polanski K., Hoo R., Sancho-Serra C., Li T., Kwakwa K., Tuck E., Lorenzi V., Massalha H., Prete M., Kleshchevnikov V., Tarkowska A., Porter T., Mazzeo C.I., van Dongen S., Dabrowska M., Vaskivskyi V., Mahbubani K.T., Park J.E., Jimenez-Linan M., Campos L., Kiselev V.Y., Lindskog C., Ayuk P., Prigmore E., Stratton M.R., Saeb-Parsy K., Moffett A., Moore L., Bayraktar O.A., Teichmann S.A., Turco M.Y, & Vento-Tormo R. (2021). Mapping the temporal and spatial dynamics of the human endometrium in vivo and in vitro. Nature Genetics, 53(12), 1698-1711.