Super-enhancer interactomes from single cells link clustering and transcription.

Regulation of gene expression hinges on the interplay between enhancers and promoters, traditionally explored through pairwise analyses. Recent advancements in mapping genome folding, like GAM, SPRITE, and multi-contact Hi-C, have uncovered multi-way interactions among super-enhancers (SEs), spanning megabases, yet have not measured their frequency in single cells or the relationship between clustering and transcription. To close this gap, here we used multiplexed imaging to map the 3D positions of 376 SEs across thousands of mammalian nuclei. Notably, our single-cell images reveal that while SE-SE contacts are rare, SEs often form looser associations we termed "communities". These communities, averaging 4-5 SEs, assemble cooperatively under the combined effects of genomic tethers, Pol2 clustering, and nuclear compartmentalization. Larger communities are associated with more frequent and larger transcriptional bursts. Our work provides insights about the SE interactome in single cells that challenge existing hypotheses on SE clustering in the context of transcriptional regulation.

Acquisition of Drug Resistance in Basal Cell Nevus Syndrome Tumors through Basal to Squamous Cell Carcinoma Transition.

Although basal cell carcinomas arise from ectopic Hedgehog pathway activation and can be treated with pathway inhibitors, sporadic basal cell carcinomas display high resistance rates, whereas tumors arising in patients with Gorlin syndrome with germline Patched (PTCH1) alterations are uniformly suppressed by inhibitor therapy. In rare cases, patients with Gorlin syndrome on long-term inhibitor therapy will develop individual resistant tumor clones that rapidly progress, but the basis of this resistance remains unstudied. In this study, we report a case of an SMO inhibitor-resistant tumor arising in a patient with Gorlin syndrome on suppressive SMO inhibitor for nearly a decade. Using a combination of multiomics and spatial transcriptomics, we define the tumor populations at the cellular and tissue level to conclude that Gorlin tumors can develop resistance to SMO inhibitors through the previously described basal to squamous cell carcinoma transition. Intriguingly, through spatial whole-exome genomic analysis, we nominate PCYT2, ETNK1, and the phosphatidylethanolamine biosynthetic pathway as genetic suppressors of basal to squamous cell carcinoma transition resistance. These observations provide a general framework for studying tumor evolution and provide important clinical insight into mechanisms of resistance to SMO inhibitors for not only Gorlin syndrome but also sporadic basal cell carcinomas.

A Spatiotemporal and Machine-Learning Platform Accelerates the Manufacturing of hPSC-derived Esophageal Mucosa.

Human pluripotent stem cell-derived tissue engineering offers great promise in designer cell-based personalized therapeutics. To harness such potential, a broader approach requires a deeper understanding of tissue-level interactions. We previously developed a manufacturing system for the ectoderm-derived skin epithelium for cell replacement therapy. However, it remains challenging to manufacture the endoderm-derived esophageal epithelium, despite both possessing similar stratified structure. Here we employ single cell and spatial technologies to generate a spatiotemporal multi-omics cell atlas for human esophageal development. We illuminate the cellular diversity, dynamics and signal communications for the developing esophageal epithelium and stroma. Using the machine-learning based Manatee, we prioritize the combinations of candidate human developmental signals for in vitro derivation of esophageal basal cells. Functional validation of the Manatee predictions leads to a clinically-compatible system for manufacturing human esophageal mucosa. Our approach creates a versatile platform to accelerate human tissue manufacturing for future cell replacement therapies to treat human genetic defects and wounds.

Acquisition of drug resistance in basal cell nevus syndrome tumors through basal to squamous cell carcinoma transition.

While basal cell carcinomas (BCCs) arise from ectopic hedgehog pathway activation and can be treated with pathway inhibitors, sporadic BCCs display high resistance rates while tumors arising in Gorlin syndrome patients with germline Patched ( PTCH1 ) mutations are uniformly suppressed by inhibitor therapy. In rare cases, Gorlin syndrome patients on long-term inhibitor therapy will develop individual resistant tumor clones that rapidly progress, but the basis of this resistance remains unstudied. Here we report a case of an SMO i -resistant tumor arising in a Gorlin patient on suppressive SMO i for nearly a decade. Using a combination of multi-omics and spatial transcriptomics, we define the tumor populations at the cellular and tissue level to conclude that Gorlin tumors can develop resistance to SMO i through the previously described basal to squamous cell carcinoma transition (BST). Intriguingly, through spatial whole exome genomic analysis, we nominate PCYT2, ETNK1, and the phosphatidylethanolamine biosynthetic pathway as novel genetic suppressors of BST resistance. These observations provide a general framework for studying tumor evolution and provide important clinical insight into mechanisms of resistance to SMO i for not only Gorlin syndrome but sporadic BCCs as well.

Skin basal cell carcinomas assemble a pro-tumorigenic spatially organized and self-propagating Trem2+ myeloid niche.

Cancer immunotherapies have revolutionized treatment but have shown limited success as single-agent therapies highlighting the need to understand the origin, assembly, and dynamics of heterogeneous tumor immune niches. Here, we use single-cell and imaging-based spatial analysis to elucidate three microenvironmental neighborhoods surrounding the heterogeneous basal cell carcinoma tumor epithelia. Within the highly proliferative neighborhood, we find that TREM2+ skin cancer-associated macrophages (SCAMs) support the proliferation of a distinct tumor epithelial population through an immunosuppression-independent manner via oncostatin-M/JAK-STAT3 signaling. SCAMs represent a unique tumor-specific TREM2+ population defined by VCAM1 surface expression that is not found in normal homeostatic skin or during wound healing. Furthermore, SCAMs actively proliferate and self-propagate through multiple serial tumor passages, indicating long-term potential. The tumor rapidly drives SCAM differentiation, with intratumoral injections sufficient to instruct naive bone marrow-derived monocytes to polarize within days. This work provides mechanistic insights into direct tumor-immune niche dynamics independent of immunosuppression, providing the basis for potential combination tumor therapies.

A scalable, GMP-compatible, autologous organotypic cell therapy for Dystrophic Epidermolysis Bullosa.

Background: Gene editing in induced pluripotent stem (iPS) cells has been hailed to enable new cell therapies for various monogenetic diseases including dystrophic epidermolysis bullosa (DEB). However, manufacturing, efficacy and safety roadblocks have limited the development of genetically corrected, autologous iPS cell-based therapies. Methods: We developed Dystrophic Epidermolysis Bullosa Cell Therapy (DEBCT), a new generation GMP-compatible (cGMP), reproducible, and scalable platform to produce autologous clinical-grade iPS cell-derived organotypic induced skin composite (iSC) grafts to treat incurable wounds of patients lacking type VII collagen (C7). DEBCT uses a combined high-efficiency reprogramming and CRISPR-based genetic correction single step to generate genome scar-free, COL7A1 corrected clonal iPS cells from primary patient fibroblasts. Validated iPS cells are converted into epidermal, dermal and melanocyte progenitors with a novel 2D organoid differentiation protocol, followed by CD49f enrichment and expansion to minimize maturation heterogeneity. iSC product characterization by single cell transcriptomics was followed by mouse xenografting for disease correcting activity at 1 month and toxicology analysis at 1-6 months. Culture-acquired mutations, potential CRISPR-off targets, and cancer-driver variants were evaluated by targeted and whole genome sequencing. Findings: iPS cell-derived iSC grafts were reproducibly generated from four recessive DEB patients with different pathogenic mutations. Organotypic iSC grafts onto immune-compromised mice developed into stable stratified skin with functional C7 restoration. Single cell transcriptomic characterization of iSCs revealed prominent holoclone stem cell signatures in keratinocytes and the recently described Gibbin-dependent signature in dermal fibroblasts. The latter correlated with enhanced graftability. Multiple orthogonal sequencing and subsequent computational approaches identified random and non-oncogenic mutations introduced by the manufacturing process. Toxicology revealed no detectable tumors after 3-6 months in DEBCT-treated mice. Interpretation: DEBCT successfully overcomes previous roadblocks and represents a robust, scalable, and safe cGMP manufacturing platform for production of a CRISPR-corrected autologous organotypic skin graft to heal DEB patient wounds.

GRHL2 and AP2a coordinate early surface ectoderm lineage commitment during development.

Ectodermal dysplasias including skin abnormalities and cleft lip/palate result from improper surface ectoderm (SE) patterning. However, the connection between SE gene regulatory networks and disease remains poorly understood. Here, we dissect human SE differentiation with multiomics and establish GRHL2 as a key mediator of early SE commitment, which acts by skewing cell fate away from the neural lineage. GRHL2 and master SE regulator AP2a balance early cell fate output, with GRHL2 facilitating AP2a binding to SE loci. In turn, AP2a restricts GRHL2 DNA binding away from de novo chromatin contacts. Integration of these regulatory sites with ectodermal dysplasia-associated genomic variants annotated within the Biomedical Data Commons identifies 55 loci previously implicated in craniofacial disorders. These include ABCA4/ARHGAP29 and NOG regulatory regions where disease-linked variants directly affect GRHL2/AP2a binding and gene transcription. These studies elucidate the logic underlying SE commitment and deepen our understanding of human oligogenic disease pathogenesis.

LY6D marks pre-existing resistant basosquamous tumor subpopulations.

Improved response to canonical therapies requires a mechanistic understanding of dynamic tumor heterogeneity by identifying discrete cellular populations with enhanced cellular plasticity. We have previously demonstrated distinct resistance mechanisms in skin basal cell carcinomas, but a comprehensive understanding of the cellular states and markers associated with these populations remains poorly understood. Here we identify a pre-existing resistant cellular population in naive basal cell carcinoma tumors marked by the surface marker LY6D. LY6D+ tumor cells are spatially localized and possess basal cell carcinoma and squamous cell carcinoma-like features. Using computational tools, organoids, and spatial tools, we show that LY6D+ basosquamous cells represent a persister population lying on a central node along the skin lineage-associated spectrum of epithelial states with local environmental and applied therapies determining the kinetics of accumulation. Surprisingly, LY6D+ basosquamous populations exist in many epithelial tumors, such as pancreatic adenocarcinomas, which have poor outcomes. Overall, our results identify the resistant LY6D+ basosquamous population as an important clinical target and suggest strategies for future therapeutic approaches to target them.

De novo and inherited variants in coding and regulatory regions in genetic cardiomyopathies.

BACKGROUND: Cardiomyopathies are a leading cause of progressive heart failure and sudden cardiac death; however, their genetic aetiology remains poorly understood. We hypothesised that variants in noncoding regulatory regions and oligogenic inheritance mechanisms may help close the diagnostic gap. METHODS: We first analysed whole-genome sequencing data of 143 parent-offspring trios from Genomics England 100,000 Genomes Project. We used gene panel testing and a phenotype-based, variant prioritisation framework called Exomiser to identify candidate genes in trios. To assess the contribution of noncoding DNVs to cardiomyopathies, we intersected DNVs with open chromatin sequences from single-cell ATAC-seq data of cardiomyocytes. We also performed a case-control analysis in an exome-negative cohort, including 843 probands and 19,467 controls, to assess the association between noncoding variants in known cardiomyopathy genes and disease. RESULTS: In the trio analysis, a definite or probable genetic diagnosis was identified in 21 probands according to the American College of Medical Genetics guidelines. We identified novel DNVs in diagnostic-grade genes (RYR2, TNNT2, PTPN11, MYH7, LZR1, NKX2-5), and five cases harbouring a combination of prioritised variants, suggesting that oligogenic inheritance and genetic modifiers contribute to cardiomyopathies. Phenotype-based ranking of candidate genes identified in noncoding DNV analysis revealed JPH2 as the top candidate. Moreover, a case-control analysis revealed an enrichment of rare noncoding variants in regulatory elements of cardiomyopathy genes (p = .035, OR = 1.43, 95% Cl = 1.095-1.767) versus controls. Of the 25 variants associated with disease  (p< 0.5), 23 are novel and nine are predicted to disrupt transcription factor binding motifs. CONCLUSION: Our results highlight complex genetic mechanisms in cardiomyopathies and reveal novel genes for future investigations.

Gibbin mesodermal regulation patterns epithelial development.

Proper ectodermal patterning during human development requires previously identified transcription factors such as GATA3 and p63, as well as positional signalling from regional mesoderm1-6. However, the mechanism by which ectoderm and mesoderm factors act to stably pattern gene expression and lineage commitment remains unclear. Here we identify the protein Gibbin, encoded by the Xia-Gibbs AT-hook DNA-binding-motif-containing 1 (AHDC1) disease gene7-9, as a key regulator of early epithelial morphogenesis. We find that enhancer- or promoter-bound Gibbin interacts with dozens of sequence-specific zinc-finger transcription factors and methyl-CpG-binding proteins to regulate the expression of mesoderm genes. The loss of Gibbin causes an increase in DNA methylation at GATA3-dependent mesodermal genes, resulting in a loss of signalling between developing dermal and epidermal cell types. Notably, Gibbin-mutant human embryonic stem-cell-derived skin organoids lack dermal maturation, resulting in p63-expressing basal cells that possess defective keratinocyte stratification. In vivo chimeric CRISPR mouse mutants reveal a spectrum of Gibbin-dependent developmental patterning defects affecting craniofacial structure, abdominal wall closure and epidermal stratification that mirror patient phenotypes. Our results indicate that the patterning phenotypes seen in Xia-Gibbs and related syndromes derive from abnormal mesoderm maturation as a result of gene-specific DNA methylation decisions.

The cell-surface 5'-nucleotidase CD73 defines a functional T memory cell subset that declines with age.

Memory T cells exhibit considerable diversity that determines their ability to be protective. Here, we examine whether changes in T cell heterogeneity contribute to the age-associated failure of immune memory. By screening for age-dependent T cell-surface markers, we identify CD4 and CD8 memory T cell subsets that are unrelated to previously defined subsets of central and effector memory cells. Memory T cells expressing the ecto-5'-nucleotidase CD73 constitute a functionally distinct subset of memory T cells that declines with age. They resemble long-lived, polyfunctional memory cells but are also poised to display effector functions and to develop into cells resembling tissue-resident memory T cells (TRMs). Upstream regulators of differential chromatin accessibility and transcriptomes include transcription factors that facilitate CD73 expression and regulate TRM differentiation. CD73 is not just a surrogate marker of these regulatory networks but is directly involved in T cell survival.

c-FOS drives reversible basal to squamous cell carcinoma transition.

While squamous transdifferentiation within subpopulations of adenocarcinomas represents an important drug resistance problem, its underlying mechanism remains poorly understood. Here, using surface markers of resistant basal cell carcinomas (BCCs) and patient single-cell and bulk transcriptomic data, we uncover the dynamic roadmap of basal to squamous cell carcinoma transition (BST). Experimentally induced BST identifies activator protein 1 (AP-1) family members in regulating tumor plasticity, and we show that c-FOS plays a central role in BST by regulating the accessibility of distinct AP-1 regulatory elements. Remarkably, despite prominent changes in cell morphology and BST marker expression, we show using inducible model systems that c-FOS-mediated BST demonstrates reversibility. Blocking EGFR pathway activation after c-FOS induction partially reverts BST in vitro and prevents BST features in both mouse models and human tumors. Thus, by identifying the molecular basis of BST, our work reveals a therapeutic opportunity targeting plasticity as a mechanism of tumor resistance.

Biomedical Data Commons (BMDC) prioritizes B-lymphocyte non-coding genetic variants in Type 1 Diabetes.

The repurposing of biomedical data is inhibited by its fragmented and multi-formatted nature that requires redundant investment of time and resources by data scientists. This is particularly true for Type 1 Diabetes (T1D), one of the most intensely studied common childhood diseases. Intense investigation of the contribution of pancreatic ß-islet and T-lymphocytes in T1D has been made. However, genetic contributions from B-lymphocytes, which are known to play a role in a subset of T1D patients, remain relatively understudied. We have addressed this issue through the creation of Biomedical Data Commons (BMDC), a knowledge graph that integrates data from multiple sources into a single queryable format. This increases the speed of analysis by multiple orders of magnitude. We develop a pipeline using B-lymphocyte multi-dimensional epigenome and connectome data and deploy BMDC to assess genetic variants in the context of Type 1 Diabetes (T1D). Pipeline-identified variants are primarily common, non-coding, poorly conserved, and are of unknown clinical significance. While variants and their chromatin connectivity are cell-type specific, they are associated with well-studied disease genes in T-lymphocytes. Candidates include established variants in the HLA-DQB1 and HLA-DRB1 and IL2RA loci that have previously been demonstrated to protect against T1D in humans and mice providing validation for this method. Others are included in the well-established T1D GRS2 genetic risk scoring method. More intriguingly, other prioritized variants are completely novel and form the basis for future mechanistic and clinical validation studies The BMDC community-based platform can be expanded and repurposed to increase the accessibility, reproducibility, and productivity of biomedical information for diverse applications including the prioritization of cell type-specific disease alleles from complex phenotypes.

AP-1 and TGFß cooperativity drives non-canonical Hedgehog signaling in resistant basal cell carcinoma.

Tumor heterogeneity and lack of knowledge about resistant cell states remain a barrier to targeted cancer therapies. Basal cell carcinomas (BCCs) depend on Hedgehog (Hh)/Gli signaling, but can develop mechanisms of Smoothened (SMO) inhibitor resistance. We previously identified a nuclear myocardin-related transcription factor (nMRTF) resistance pathway that amplifies noncanonical Gli1 activity, but characteristics and drivers of the nMRTF cell state remain unknown. Here, we use single cell RNA-sequencing of patient tumors to identify three prognostic surface markers (LYPD3, TACSTD2, and LY6D) which correlate with nMRTF and resistance to SMO inhibitors. The nMRTF cell state resembles transit-amplifying cells of the hair follicle matrix, with AP-1 and TGFß cooperativity driving nMRTF activation. JNK/AP-1 signaling commissions chromatin accessibility and Smad3 DNA binding leading to a transcriptional program of RhoGEFs that facilitate nMRTF activity. Importantly, small molecule AP-1 inhibitors selectively target LYPD3+/TACSTD2+/LY6D+ nMRTF human BCCs ex vivo, opening an avenue for improving combinatorial therapies.

Endogenous Retrovirus-Derived lncRNA BANCR Promotes Cardiomyocyte Migration in Humans and Non-human Primates.

Transposable elements (TEs) comprise nearly half of the human genome and are often transcribed or exhibit cis-regulatory properties with unknown function in specific processes such as heart development. In the case of endogenous retroviruses (ERVs), a TE subclass, experimental interrogation is constrained as many are primate-specific or human-specific. Here, we use primate pluripotent stem-cell-derived cardiomyocytes that mimic fetal cardiomyocytes in vitro to discover hundreds of ERV transcripts from the primate-specific MER41 family, some of which are regulated by the cardiogenic transcription factor TBX5. The most significant of these are located within BANCR, a long non-coding RNA (lncRNA) exclusively expressed in primate fetal cardiomyocytes. Functional studies reveal that BANCR promotes cardiomyocyte migration in vitro and ventricular enlargement in vivo. We conclude that recently evolved TE loci such as BANCR may represent potent de novo developmental regulatory elements that can be interrogated with species-matching pluripotent stem cell models.

Clinical trial in a dish using iPSCs shows lovastatin improves endothelial dysfunction and cellular cross-talk in LMNA cardiomyopathy.

Mutations in LMNA, the gene that encodes lamin A and C, causes LMNA-related dilated cardiomyopathy (DCM) or cardiolaminopathy. LMNA is expressed in endothelial cells (ECs); however, little is known about the EC-specific phenotype of LMNA-related DCM. Here, we studied a family affected by DCM due to a frameshift variant in LMNA Human induced pluripotent stem cell (iPSC)-derived ECs were generated from patients with LMNA-related DCM and phenotypically characterized. Patients with LMNA-related DCM exhibited clinical endothelial dysfunction, and their iPSC-ECs showed decreased functionality as seen by impaired angiogenesis and nitric oxide (NO) production. Moreover, genome-edited isogenic iPSC lines recapitulated the EC disease phenotype in which LMNA-corrected iPSC-ECs showed restoration of EC function. Simultaneous profiling of chromatin accessibility and gene expression dynamics by combining assay for transposase-accessible chromatin using sequencing (ATAC-seq) and RNA sequencing (RNA-seq) as well as loss-of-function studies identified Krüppel-like factor 2 (KLF2) as a potential transcription factor responsible for the EC dysfunction. Gain-of-function studies showed that treatment of LMNA iPSC-ECs with KLF2 agonists, including lovastatin, rescued the EC dysfunction. Patients with LMNA-related DCM treated with lovastatin showed improvements in clinical endothelial dysfunction as indicated by increased reactive hyperemia index. Furthermore, iPSC-derived cardiomyocytes (iPSC-CMs) from patients exhibiting the DCM phenotype showed improvement in CM function when cocultured with iPSC-ECs and lovastatin. These results suggest that impaired cross-talk between ECs and CMs can contribute to the pathogenesis of LMNA-related DCM, and statin may be an effective therapy for vascular dysfunction in patients with cardiolaminopathy.

TFAP2C- and p63-Dependent Networks Sequentially Rearrange Chromatin Landscapes to Drive Human Epidermal Lineage Commitment.

Tissue development results from lineage-specific transcription factors (TFs) programming a dynamic chromatin landscape through progressive cell fate transitions. Here, we define epigenomic landscape during epidermal differentiation of human pluripotent stem cells (PSCs) and create inference networks that integrate gene expression, chromatin accessibility, and TF binding to define regulatory mechanisms during keratinocyte specification. We found two critical chromatin networks during surface ectoderm initiation and keratinocyte maturation, which are driven by TFAP2C and p63, respectively. Consistently, TFAP2C, but not p63, is sufficient to initiate surface ectoderm differentiation, and TFAP2C-initiated progenitor cells are capable of maturing into functional keratinocytes. Mechanistically, TFAP2C primes the surface ectoderm chromatin landscape and induces p63 expression and binding sites, thus allowing maturation factor p63 to positively autoregulate its own expression and close a subset of the TFAP2C-initiated surface ectoderm program. Our work provides a general framework to infer TF networks controlling chromatin transitions that will facilitate future regenerative medicine advances.