Single-cell multiomic analysis identifies a HOX-PBX gene network regulating the survival of lymphangioleiomyomatosis cells

Lymphangioleiomyomatosis (LAM) is a rare, progressive lung disease that predominantly affects women. LAM cells carry TSC1/TSC2 mutations, causing mTORC1 hyperactivation and uncontrolled cell growth. mTORC1 inhibitors stabilize lung function; however, sustained efficacy requires long-term administration, and some patients fail to tolerate or respond to therapy. Although the genetic basis of LAM is known, mechanisms underlying LAM pathogenesis remain elusive. We integrated single-cell RNA sequencing and single-nuclei ATAC-seq of LAM lungs to construct a gene regulatory network controlling the transcriptional program of LAM cells. We identified activation of uterine-specific HOX-PBX transcriptional programs in pulmonary LAMCORE cells as regulators of cell survival depending upon HOXD11-PBX1 dimerization. Accordingly, blockage of HOXD11-PBX1 dimerization by HXR9 suppressed LAM cell survival in vitro and in vivo. PBX1 regulated STAT1/3, increased the expression of antiapoptotic genes, and promoted LAM cell survival in vitro. The HOX-PBX gene network provides promising targets for treatment of LAM/TSC mTORC1-hyperactive cancers.

[1]  Timothy L. Tickle,et al.  Guided construction of single cell reference for human and mouse lung , 2022, bioRxiv.

[2]  Suzanne Miller,et al.  Cross-talk between LAM cells and fibroblasts may influence alveolar epithelial cell behavior in Lymphangioleiomyomatosis. , 2021, American journal of physiology. Lung cellular and molecular physiology.

[3]  C. Lareau,et al.  Single-cell chromatin state analysis with Signac , 2021, Nature Methods.

[4]  A. R. Costa-Pinto,et al.  Roles of the HOX Proteins in Cancer Invasion and Metastasis , 2020, Cancers.

[5]  Avi Srivastava,et al.  Multimodal single-cell chromatin analysis with Signac , 2020, bioRxiv.

[6]  Kyle J. Gaulton,et al.  Single-cell multiomic profiling of human lungs reveals cell-type-specific and age-dynamic control of SARS-CoV2 host genes , 2020, eLife.

[7]  Henry W. Long,et al.  MITF is a driver oncogene and potential therapeutic target in kidney angiomyolipoma tumors through transcriptional regulation of CYR61 , 2020, Oncogene.

[8]  D. Kwiatkowski,et al.  Single Cell Transcriptomic Analysis Identifies a Unique Pulmonary Lymphangioleiomyomatosis Cell. , 2020, American journal of respiratory and critical care medicine.

[9]  Yan Lu,et al.  Single-Cell Transcriptomic Analysis. , 2020, Comprehensive Physiology.

[10]  Tudor I. Oprea,et al.  Therapies for rare diseases: therapeutic modalities, progress and challenges ahead , 2019, Nature Reviews Drug Discovery.

[11]  A. Shilatifard,et al.  Single-Cell Transcriptomic Analysis of Human Lung Provides Insights into the Pathobiology of Pulmonary Fibrosis , 2019, American journal of respiratory and critical care medicine.

[12]  E. Silverman,et al.  A genome-wide association study implicates NR2F2 in lymphangioleiomyomatosis pathogenesis , 2019, European Respiratory Journal.

[13]  Paul J. Hoffman,et al.  Comprehensive Integration of Single-Cell Data , 2018, Cell.

[14]  Allon M Klein,et al.  Scrublet: Computational Identification of Cell Doublets in Single-Cell Transcriptomic Data. , 2019, Cell systems.

[15]  Andrew J. Hill,et al.  The single cell transcriptional landscape of mammalian organogenesis , 2019, Nature.

[16]  Samriddha Ray,et al.  Single cell RNA analysis identifies cellular heterogeneity and adaptive responses of the lung at birth , 2019, Nature Communications.

[17]  Fan Zhang,et al.  Fast, sensitive, and accurate integration of single cell data with Harmony , 2018, bioRxiv.

[18]  Christoph Hafemeister,et al.  Comprehensive integration of single cell data , 2018, bioRxiv.

[19]  R. Morgan,et al.  Targeting HOX-PBX interactions causes death in oral potentially malignant and squamous carcinoma cells but not normal oral keratinocytes , 2018, BMC Cancer.

[20]  Matthew D. Young,et al.  SoupX removes ambient RNA contamination from droplet-based single-cell RNA sequencing data , 2018, bioRxiv.

[21]  Laleh Haghverdi,et al.  Batch effects in single-cell RNA-sequencing data are corrected by matching mutual nearest neighbors , 2018, Nature Biotechnology.

[22]  Ruikang K. Wang,et al.  Robust principal component analysis in optical micro-angiography. , 2017, Quantitative imaging in medicine and surgery.

[23]  W. Wong,et al.  Modeling gene regulation from paired expression and chromatin accessibility data , 2017, Proceedings of the National Academy of Sciences.

[24]  R. Morgan,et al.  Targeting HOX/PBX dimers in cancer , 2017, Oncotarget.

[25]  R. Morgan,et al.  HOX transcription factors are potential targets and markers in malignant mesothelioma , 2016, BMC Cancer.

[26]  S. Langley,et al.  Importance of HOX genes in normal prostate gland formation, prostate cancer development and its early detection , 2014, BJU international.

[27]  S. Hammes,et al.  Uterine-specific loss of Tsc2 leads to myometrial tumors in both the uterus and lungs. , 2013, Molecular endocrinology.

[28]  W. Travis,et al.  Lymphangioleiomyomatosis: calling it what it is: a low-grade, destructive, metastasizing neoplasm. , 2012, American journal of respiratory and critical care medicine.

[29]  Kazuhisa Takahashi,et al.  Prevalence of Uterine and Adnexal Involvement in Pulmonary Lymphangioleiomyomatosis: A Clinicopathologic Study of 10 Patients , 2011, The American journal of surgical pathology.

[30]  Mathieu Lupien,et al.  PBX1 Genomic Pioneer Function Drives ERα Signaling Underlying Progression in Breast Cancer , 2011, PLoS genetics.

[31]  Y. Kida Efficacy and safety of sirolimus in lymphangioleiomyomatosis. , 2011, The New England journal of medicine.

[32]  Philip Machanick,et al.  MEME-ChIP: motif analysis of large DNA datasets , 2011, Bioinform..

[33]  Chuhsing Kate Hsiao,et al.  Identification of a Novel Biomarker, SEMA5A, for Non–Small Cell Lung Carcinoma in Nonsmoking Women , 2010, Cancer Epidemiology, Biomarkers & Prevention.

[34]  Jilly F. Evans,et al.  Effects of prolactin on TSC2-null Eker rat cells and in pulmonary lymphangioleiomyomatosis. , 2010, American journal of respiratory and critical care medicine.

[35]  Kevin J Harrington,et al.  Targeting HOX and PBX transcription factors in ovarian cancer , 2010, BMC Cancer.

[36]  E. Goncharova,et al.  Signal Transducer and Activator of Transcription 3 Is Required for Abnormal Proliferation and Survival of TSC2-Deficient Cells: Relevance to Pulmonary Lymphangioleiomyomatosis , 2009, Molecular Pharmacology.

[37]  Jing Chen,et al.  ToppGene Suite for gene list enrichment analysis and candidate gene prioritization , 2009, Nucleic Acids Res..

[38]  E. Henske,et al.  Estrogen promotes the survival and pulmonary metastasis of tuberin-null cells , 2009, Proceedings of the National Academy of Sciences.

[39]  J. Moss,et al.  Chemokine-Enhanced Chemotaxis of Lymphangioleiomyomatosis Cells with Mutations in the Tumor Suppressor TSC2 Gene1 , 2009, The Journal of Immunology.

[40]  R. Morgan,et al.  HOX transcription factors are potential therapeutic targets in non-small-cell lung cancer (targeting HOX genes in lung cancer) , 2009, British Journal of Cancer.

[41]  D. Kwiatkowski,et al.  mTOR-raptor binds and activates SGK1 to regulate p27 phosphorylation. , 2008, Molecular cell.

[42]  J. Moss,et al.  Lymphatic Involvement in Lymphangioleiomyomatosis , 2008, Annals of the New York Academy of Sciences.

[43]  S. Wacholder,et al.  Gene Expression Signature of Cigarette Smoking and Its Role in Lung Adenocarcinoma Development and Survival , 2008, PloS one.

[44]  R. Morgan,et al.  Antagonism of HOX/PBX dimer formation blocks the in vivo proliferation of melanoma. , 2007, Cancer research.

[45]  O. Pourquié,et al.  Hox genes in time and space during vertebrate body formation , 2007, Development, growth & differentiation.

[46]  G. Russ,et al.  Superior Outcomes in Renal Transplantation after Early Cyclosporine Withdrawal and Sirolimus Maintenance Therapy, Regardless of Baseline Renal Function , 2005, Transplantation.

[47]  Hongbing Zhang,et al.  Estrogen-induced Smooth Muscle Cell Growth Is Regulated by Tuberin and Associated with Altered Activation of Platelet-derived Growth Factor Receptor-β and ERK-1/2* , 2004, Journal of Biological Chemistry.

[48]  Hongbing Zhang,et al.  Perturbed IFN-γ-Jak-Signal Transducers and Activators of Transcription Signaling in Tuberous Sclerosis Mouse Models , 2004, Cancer Research.

[49]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[50]  Joel Moss,et al.  Stimulation of Signal Transducer and Activator of Transcription-1 (STAT1)-dependent Gene Transcription by Lipopolysaccharide and Interferon-γ Is Regulated by Mammalian Target of Rapamycin* , 2003, Journal of Biological Chemistry.

[51]  K. Nilsson,et al.  Ser727/Tyr701-phosphorylated Stat1 is required for the regulation of c-Myc, cyclins, and p27Kip1 associated with ATRA-induced G0/G1 arrest of U-937 cells. , 2003, Blood.

[52]  J. Testa,et al.  Recurrent lymphangiomyomatosis after transplantation: genetic analyses reveal a metastatic mechanism. , 2003, American journal of respiratory and critical care medicine.

[53]  Rossella Fattori,et al.  Drug-eluting stents in vascular intervention , 2003, The Lancet.

[54]  E. Henske,et al.  Tuberin, the tuberous sclerosis complex 2 tumor suppressor gene product, regulates Rho activation, cell adhesion and migration , 2002, Oncogene.

[55]  E. Henske,et al.  Mutational analysis of the tuberous sclerosis geneTSC2 in patients with pulmonary lymphangioleiomyomatosis , 2000, Journal of medical genetics.

[56]  J. Ritz,et al.  Fludarabine-induced immunosuppression is associated with inhibition of STAT1 signaling , 1999, Nature Medicine.

[57]  C. Walker,et al.  Estrogen stimulation and tamoxifen inhibition of leiomyoma cell growth in vitro and in vivo. , 1995, Endocrinology.

[58]  L. Rassenti,et al.  Oral fludarabine in untreated patients with B-cell chronic lymphocytic leukemia , 2010, Leukemia.

[59]  D. Kwiatkowski,et al.  Interferon-gamma-Jak-Stat signaling in pulmonary lymphangioleiomyomatosis and renal angiomyolipoma: a potential therapeutic target. , 2005, American journal of respiratory cell and molecular biology.

[60]  W. Plunkett,et al.  Cellular and Clinical Pharmacology of Fludarabine , 2002, Clinical pharmacokinetics.