Cellular landscape of adrenocortical carcinoma at single-nuclei resolution

Adrenocortical carcinoma (ACC) is a rare yet devastating tumour of the adrenal gland with a molecular pathology that remains incompletely understood. To gain novel insights into the cellular landscape of ACC, we generated single-nuclei RNA sequencing (snRNA-seq) data sets from twelve ACC tumour samples and analysed these alongside snRNA-seq data sets from normal adrenal glands (NAGs). We find the ACC tumour microenvironment to be relatively devoid of immune cells compared to NAG tissues, consistent with known high tumour purity values for ACC as an immunologically "cold" tumour. Our analysis identifies three separate groups of ACC samples that are characterised by different relative compositions of adrenocortical cell types. These include cell populations that are specifically enriched in the most clinically aggressive and hormonally active tumours, displaying hallmarks of reorganised cell mechanobiology and dysregulated steroidogenesis, respectively. We also identified and validated a population of mitotically active adrenocortical cells that strongly overexpress genes POLQ, DIAPH3 and EZH2 to support tumour expansion alongside an LGR4+ progenitor-like or cell-of-origin candidate for adrenocortical carcinogenesis. Trajectory inference suggests the fate adopted by malignant adrenocortical cells upon differentiation is associated with the copy number or allelic balance state of the imprinted DLK1/MEG3 genomic locus, which we verified by assessing bulk tumour DNA methylation status. In conclusion, our results therefore provide new insights into the clinical and cellular heterogeneity of ACC, revealing how genetic perturbations to healthy adrenocortical renewal and zonation provide a molecular basis for disease pathogenesis.