Rigorously comparing gene expression and chromatin accessibility in the same single cells could illuminate the logic of how coupling or decoupling of these mechanisms regulates fate commitment. Here we present MIRA, probabilistic multimodal models for integrated regulatory analysis, a comprehensive methodology that systematically contrasts transcription and accessibility to infer the regulatory circuitry driving cells along cell state trajectories. MIRA leverages topic modeling of cell states and regulatory potential modeling of individual gene loci. MIRA thereby represents cell states in an efficient and interpretable latent space, infers high-fidelity cell state trees, determines key regulators of fate decisions at branch points and exposes the variable influence of local accessibility on transcription at distinct loci. Applied to epidermal differentiation and embryonic brain development from two different multimodal platforms, MIRA revealed that early developmental genes were tightly regulated by local chromatin landscape whereas terminal fate genes were titrated without requiring extensive chromatin remodeling.

Pharmacologic inhibitors of cyclin-dependent kinases 4 and 6 (CDK4/6) were designed to induce cancer cell cycle arrest. Recent studies have suggested that these agents also exert other effects, influencing cancer cell immunogenicity, apoptotic responses, and differentiation. Using cell-based and mouse models of breast cancer together with clinical specimens, we show that CDK4/6 inhibitors induce remodeling of cancer cell chromatin characterized by widespread enhancer activation, and that this explains many of these effects. The newly activated enhancers include classical super-enhancers that drive luminal differentiation and apoptotic evasion, as well as a set of enhancers overlying endogenous retroviral elements that is enriched for proximity to interferon-driven genes. Mechanistically, CDK4/6 inhibition increases the level of several Activator Protein-1 (AP-1) transcription factor proteins, which are in turn implicated in the activity of many of the new enhancers. Our findings offer insights into CDK4/6 pathway biology and should inform the future development of CDK4/6 inhibitors.

Giant congenital melanocytic nevi are NRAS-driven proliferations that may cover up to 80% of the body surface. Their most dangerous consequence is progression to melanoma. This risk often triggers preemptive extensive surgical excisions in childhood, producing severe lifelong challenges. We have presented preclinical models, including multiple genetically engineered mice and xenografted human lesions, which enabled testing locally applied pharmacologic agents to avoid surgery. The murine models permitted the identification of proliferative versus senescent nevus phases and treatments targeting both. These nevi recapitulated the histologic and molecular features of human giant congenital nevi, including the risk of melanoma transformation. Cutaneously delivered MEK, PI3K, and c-KIT inhibitors or proinflammatory squaric acid dibutylester (SADBE) achieved major regressions. SADBE triggered innate immunity that ablated detectable nevocytes, fully prevented melanoma, and regressed human giant nevus xenografts. These findings reveal nevus mechanistic vulnerabilities and suggest opportunities for topical interventions that may alter the therapeutic options for children with congenital giant nevi.