Dictys: dynamic gene regulatory network dissects developmental continuum with single-cell multiomics.

Gene regulatory networks (GRNs) are key determinants of cell function and identity and are dynamically rewired during development and disease. Despite decades of advancement, challenges remain in GRN inference, including dynamic rewiring, causal inference, feedback loop modeling and context specificity. To address these challenges, we develop Dictys, a dynamic GRN inference and analysis method that leverages multiomic single-cell assays of chromatin accessibility and gene expression, context-specific transcription factor footprinting, stochastic process network and efficient probabilistic modeling of single-cell RNA-sequencing read counts. Dictys improves GRN reconstruction accuracy and reproducibility and enables the inference and comparative analysis of context-specific and dynamic GRNs across developmental contexts. Dictys' network analyses recover unique insights in human blood and mouse skin development with cell-type-specific and dynamic GRNs. Its dynamic network visualizations enable time-resolved discovery and investigation of developmental driver transcription factors and their regulated targets. Dictys is available as a free, open-source and user-friendly Python package.

MAF functions as a pioneer transcription factor that initiates and sustains myelomagenesis.

Deregulated expression of lineage-affiliated transcription factors (TFs) is a major mechanism of oncogenesis. However, how the deregulation of nonlineage affiliated TF affects chromatin to initiate oncogenic transcriptional programs is not well-known. To address this, we studied the chromatin effects imposed by oncogenic MAF as the cancer-initiating driver in the plasma cell cancer multiple myeloma. We found that the ectopically expressed MAF endows myeloma plasma cells with migratory and proliferative transcriptional potential. This potential is regulated by the activation of enhancers and superenhancers, previously inactive in healthy B cells and plasma cells, and the cooperation of MAF with the plasma cell-defining TF IRF4. Forced ectopic MAF expression confirms the de novo ability of oncogenic MAF to convert transcriptionally inert chromatin to active chromatin with the features of superenhancers, leading to the activation of the MAF-specific oncogenic transcriptome and the acquisition of cancer-related cellular phenotypes such as CCR1-dependent cell migration. These findings establish oncogenic MAF as a pioneer transcription factor that can initiate as well as sustain oncogenic transcriptomes and cancer phenotypes. However, despite its pioneer function, myeloma cells remain MAF-dependent, thus validating oncogenic MAF as a therapeutic target that would be able to circumvent the challenges of subsequent genetic diversification driving disease relapse and drug resistance.

Systems medicine dissection of chr1q-amp reveals a novel PBX1-FOXM1 axis for targeted therapy in multiple myeloma.

Understanding the biological and clinical impact of copy number aberrations (CNAs) on the development of precision therapies in cancer remains an unmet challenge. Genetic amplification of chromosome 1q (chr1q-amp) is a major CNA conferring an adverse prognosis in several types of cancer, including in the blood cancer multiple myeloma (MM). Although several genes across chromosome 1 (chr1q) portend high-risk MM disease, the underpinning molecular etiology remains elusive. Here, with reference to the 3-dimensional (3D) chromatin structure, we integrate multi-omics data sets from patients with MM with genetic variables to obtain an associated clinical risk map across chr1q and to identify 103 adverse prognosis genes in chr1q-amp MM. Prominent among these genes, the transcription factor PBX1 is ectopically expressed by genetic amplification and epigenetic activation of its own preserved 3D regulatory domain. By binding to reprogrammed superenhancers, PBX1 directly regulates critical oncogenic pathways and a FOXM1-dependent transcriptional program. Together, PBX1 and FOXM1 activate a proliferative gene signature that predicts adverse prognosis across multiple types of cancer. Notably, pharmacological disruption of the PBX1-FOXM1 axis with existing agents (thiostrepton) and a novel PBX1 small molecule inhibitor (T417) is selectively toxic against chr1q-amp myeloma and solid tumor cells. Overall, our systems medicine approach successfully identifies CNA-driven oncogenic circuitries, links them to clinical phenotypes, and proposes novel CNA-targeted therapy strategies in MM and other types of cancer.

The innate sensor ZBP1-IRF3 axis regulates cell proliferation in multiple myeloma.

Multiple myeloma is a malignancy of plasma cells initiated and driven by primary and secondary genetic events. However, myeloma plasma cell survival and proliferation might be sustained by non-genetic drivers. Z-DNA-binding protein 1 (ZBP1; also known as DAI) is an interferon-inducible, Z-nucleic acid sensor that triggers RIPK3-MLKL-mediated necroptosis in mice. ZBP1 also interacts with TBK1 and the transcription factor IRF3 but the function of this interaction is unclear, and the role of the ZBP1-IRF3 axis in cancer is not known. Here we show that ZBP1 is selectively expressed in late B-cell development in both human and murine cells and it is required for optimal T-cell-dependent humoral immune responses. In myeloma plasma cells, the interaction of constitutively expressed ZBP1 with TBK1 and IRF3 results in IRF3 phosphorylation. IRF3 directly binds and activates cell cycle genes, in part through co-operation with the plasma cell lineage-defining transcription factor IRF4, thereby promoting myeloma cell proliferation. This generates a novel, potentially therapeutically targetable and relatively selective myeloma cell addiction to the ZBP1-IRF3 axis. Our data also show a noncanonical function of constitutive ZBP1 in human cells and expand our knowledge of the role of cellular immune sensors in cancer biology.

Chromatin-based, in cis and in trans regulatory rewiring underpins distinct oncogenic transcriptomes in multiple myeloma.

Multiple myeloma is a genetically heterogeneous cancer of the bone marrow plasma cells (PC). Distinct myeloma transcriptome profiles are primarily driven by myeloma initiating events (MIE) and converge into a mutually exclusive overexpression of the CCND1 and CCND2 oncogenes. Here, with reference to their normal counterparts, we find that myeloma PC enhanced chromatin accessibility combined with paired transcriptome profiling can classify MIE-defined genetic subgroups. Across and within different MM genetic subgroups, we ascribe regulation of genes and pathways critical for myeloma biology to unique or shared, developmentally activated or de novo formed candidate enhancers. Such enhancers co-opt recruitment of existing transcription factors, which although not transcriptionally deregulated per se, organise aberrant gene regulatory networks that help identify myeloma cell dependencies with prognostic impact. Finally, we identify and validate the critical super-enhancer that regulates ectopic expression of CCND2 in a subset of patients with MM and in chronic lymphocytic leukemia.

Single-cell profiling of human bone marrow progenitors reveals mechanisms of failing erythropoiesis in Diamond-Blackfan anemia.

Ribosome dysfunction underlies the pathogenesis of many cancers and heritable ribosomopathies. Here, we investigate how mutations in either ribosomal protein large (RPL) or ribosomal protein small (RPS) subunit genes selectively affect erythroid progenitor development and clinical phenotypes in Diamond-Blackfan anemia (DBA), a rare ribosomopathy with limited therapeutic options. Using single-cell assays of patient-derived bone marrow, we delineated two distinct cellular trajectories segregating with ribosomal protein genotypes. Almost complete loss of erythroid specification was observed in RPS-DBA. In contrast, we observed relative preservation of qualitatively abnormal erythroid progenitors and precursors in RPL-DBA. Although both DBA genotypes exhibited a proinflammatory bone marrow milieu, RPS-DBA was characterized by erythroid differentiation arrest, whereas RPL-DBA was characterized by preserved GATA1 expression and activity. Compensatory stress erythropoiesis in RPL-DBA exhibited disordered differentiation underpinned by an altered glucocorticoid molecular signature, including reduced ZFP36L2 expression, leading to milder anemia and improved corticosteroid response. This integrative analysis approach identified distinct pathways of erythroid failure and defined genotype-phenotype correlations in DBA. These findings may help facilitate therapeutic target discovery.

Anopheles coluzzii stearoyl-CoA desaturase is essential for adult female survival and reproduction upon blood feeding.

Vitellogenesis and oocyte maturation require anautogenous female Anopheles mosquitoes to obtain a bloodmeal from a vertebrate host. The bloodmeal is rich in proteins that are readily broken down into amino acids in the midgut lumen and absorbed by the midgut epithelial cells where they are converted into lipids and then transported to other tissues including ovaries. The stearoyl-CoA desaturase (SCD) plays a pivotal role in this process by converting saturated (SFAs) to unsaturated (UFAs) fatty acids; the latter being essential for maintaining cell membrane fluidity amongst other housekeeping functions. Here, we report the functional and phenotypic characterization of SCD1 in the malaria vector mosquito Anopheles coluzzii. We show that RNA interference (RNAi) silencing of SCD1 and administration of sterculic acid (SA), a small molecule inhibitor of SCD1, significantly impact on the survival and reproduction of female mosquitoes following blood feeding. Microscopic observations reveal that the mosquito thorax is quickly filled with blood, a phenomenon likely caused by the collapse of midgut epithelial cell membranes, and that epithelial cells are depleted of lipid droplets and oocytes fail to mature. Transcriptional profiling shows that genes involved in protein, lipid and carbohydrate metabolism and immunity-related genes are the most affected by SCD1 knock down (KD) in blood-fed mosquitoes. Metabolic profiling reveals that these mosquitoes exhibit increased amounts of saturated fatty acids and TCA cycle intermediates, highlighting the biochemical framework by which the SCD1 KD phenotype manifests as a result of a detrimental metabolic syndrome. Accumulation of SFAs is also the likely cause of the potent immune response observed in the absence of infection, which resembles an auto-inflammatory condition. These data provide insights into mosquito bloodmeal metabolism and lipid homeostasis and could inform efforts to develop novel interventions against mosquito-borne diseases.

Current progress and potential opportunities to infer single-cell developmental trajectory and cell fate.

Rapid technological advances in transcriptomics and lineage tracing technologies provide new opportunities to understand organismal development at the single-cell level. Building on these advances, various computational methods have been proposed to infer developmental trajectories and to predict cell fate. These methods have unveiled previously uncharacterized transitional cell types and differentiation processes. Importantly, the ability to recover cell states and trajectories has been evolving hand-in-hand with new technologies and diverse experimental designs; more recent methods can capture complex trajectory topologies and infer short- and long-term cell fate dynamics. Here, we summarize and categorize the most recent and popular computational approaches for trajectory inference based on the information they leverage and describe future challenges and opportunities for the development of new methods for reconstructing differentiation trajectories and inferring cell fates.

Brd2/4 and Myc regulate alternative cell lineage programmes during early osteoclast differentiation in vitro.

Osteoclast (OC) development in response to nuclear factor kappa-Β ligand (RANKL) is critical for bone homeostasis in health and in disease. The early and direct chromatin regulatory changes imparted by the BET chromatin readers Brd2-4 and OC-affiliated transcription factors (TFs) during osteoclastogenesis are not known. Here, we demonstrate that in response to RANKL, early OC development entails regulation of two alternative cell fate transcriptional programmes, OC vs macrophage, with repression of the latter following activation of the former. Both programmes are regulated in a non-redundant manner by increased chromatin binding of Brd2 at promoters and of Brd4 at enhancers/super-enhancers. Myc, the top RANKL-induced TF, regulates OC development in co-operation with Brd2/4 and Max and by establishing negative and positive regulatory loops with other lineage-affiliated TFs. These insights into the transcriptional regulation of osteoclastogenesis suggest the clinical potential of selective targeting of Brd2/4 to abrogate pathological OC activation.