Faculty Bibliography

BACKGROUND:Developing genomic resources for a diverse range of species is an important step towards understanding the mechanisms underlying complex traits. Specifically, organisms that exhibit unique and accessible phenotypes-of-interest allow researchers to address questions that may be ill-suited to traditional model organisms. We sequenced the genome and transcriptome of Alston's singing mouse (Scotinomys teguina), an emerging model for social cognition and vocal communication. In addition to producing advertisement songs used for mate attraction and male-male competition, these rodents are diurnal, live at high-altitudes, and are obligate insectivores, providing opportunities to explore diverse physiological, ecological, and evolutionary questions. RESULTS:Using PromethION, Illumina, and PacBio sequencing, we produced an annotated genome and transcriptome, which were validated using gene expression and functional enrichment analyses. To assess the usefulness of our assemblies, we performed single nuclei sequencing on cells of the orofacial motor cortex, a brain region implicated in song coordination, identifying 12 cell types. CONCLUSIONS:These resources will provide the opportunity to identify the molecular basis of complex traits in singing mice as well as to contribute data that can be used for large-scale comparative analyses.

We introduce a pioneering approach that integrates pathology imaging with transcriptomics and proteomics to identify predictive histology features associated with critical clinical outcomes in cancer. We utilize 2,755 H&E-stained histopathological slides from 657 patients across 6 cancer types from CPTAC. Our models effectively recapitulate distinctions readily made by human pathologists: tumor vs. normal (AUROC = 0.995) and tissue-of-origin (AUROC = 0.979). We further investigate predictive power on tasks not normally performed from H&E alone, including TP53 prediction and pathologic stage. Importantly, we describe predictive morphologies not previously utilized in a clinical setting. The incorporation of transcriptomics and proteomics identifies pathway-level signatures and cellular processes driving predictive histology features. Model generalizability and interpretability is confirmed using TCGA. We propose a classification system for these tasks, and suggest potential clinical applications for this integrated human and machine learning approach. A publicly available web-based platform implements these models.

Immune-checkpoint blockade has revolutionized cancer treatment, but some cancers, such as acute myeloid leukemia (AML), do not respond or develop resistance. A potential mode of resistance is immune evasion of T cell immunity involving aberrant major histocompatibility complex class I (MHC-I) antigen presentation (AP). To map such mechanisms of resistance, we identified key MHC-I regulators using specific peptide-MHC-I-guided CRISPR-Cas9 screens in AML. The top-ranked negative regulators were surface protein sushi domain containing 6 (SUSD6), transmembrane protein 127 (TMEM127), and the E3 ubiquitin ligase WWP2. SUSD6 is abundantly expressed in AML and multiple solid cancers, and its ablation enhanced MHC-I AP and reduced tumor growth in a CD8⁺ T cell-dependent manner. Mechanistically, SUSD6 forms a trimolecular complex with TMEM127 and MHC-I, which recruits WWP2 for MHC-I ubiquitination and lysosomal degradation. Together with the SUSD6/TMEM127/WWP2 gene signature, which negatively correlates with cancer survival, our findings define a membrane-associated MHC-I inhibitory axis as a potential therapeutic target for both leukemia and solid cancers.

The chromosome 9p21 locus comprises several tumor suppressor genes including MTAP, CDKN2A and CDKN2B, and its homo- or heterozygous deletion is associated with reduced survival in multiple cancer types. We report that mice with germline monoallelic deletion or induced biallelic deletion of the 9p21-syntenic locus (9p21s) developed a fatal myelodysplastic syndrome/myeloproliferative neoplasm (MDS/MPN)-like disease associated with aberrant trabecular bone formation and/or fibrosis in the bone marrow (BM). Reciprocal BM transfers and conditional targeting of 9p21s suggested that the disease originates in the BM stroma. Single-cell analysis of 9p21s-deficient BM stroma revealed the expansion of chondrocyte and osteogenic precursors, reflected in increased osteogenic differentiation in vitro. It also showed reduced expression of factors maintaining hematopoietic stem/progenitor cells, including Cxcl12. Accordingly, 9p21s-deficient mice showed reduced levels of circulating Cxcl12 and concomitant upregulation of the pro-fibrotic chemokine Cxcl13 and osteogenesis- and fibrosis-related multifunctional glycoprotein Osteopontin (OPN)/Spp1. Our study highlights the potential of mutations in the BM microenvironment to drive MDS/MPN-like disease.

Investigating how chromatin organization determines cell-type-specific gene expression remains challenging. Experimental methods for measuring three-dimensional chromatin organization, such as Hi-C, are costly and have technical limitations, restricting their broad application particularly in high-throughput genetic perturbations. We present C.Origami, a multimodal deep neural network that performs de novo prediction of cell-type-specific chromatin organization using DNA sequence and two cell-type-specific genomic features-CTCF binding and chromatin accessibility. C.Origami enables in silico experiments to examine the impact of genetic changes on chromatin interactions. We further developed an in silico genetic screening approach to assess how individual DNA elements may contribute to chromatin organization and to identify putative cell-type-specific trans-acting regulators that collectively determine chromatin architecture. Applying this approach to leukemia cells and normal T cells, we demonstrate that cell-type-specific in silico genetic screening, enabled by C.Origami, can be used to systematically discover novel chromatin regulation circuits in both normal and disease-related biological systems.

BH3-mimetics are used as an efficient strategy to induce cell death in several blood malignancies, including acute myeloid leukemia (AML). Venetoclax, a potent BCL-2 antagonist, is used clinically in combination with hypomethylating agents for the treatment of AML. Moreover, MCL-1 or dual BCL-2/BCL-xL antagonists are under investigation. Yet, resistance to single or combinatorial BH3-mimetics therapies eventually ensues. Integration of multiple genome-wide CRISPR/Cas9 screens revealed that loss of mitophagy modulators sensitizes AML cells to various BH3-mimetics targeting different BCL-2 family members. One such regulator is MFN2, whose protein levels positively correlate with drug resistance in patients with AML. MFN2 overexpression is sufficient to drive resistance to BH3-mimetics in AML. Insensitivity to BH3-mimetics is accompanied by enhanced mitochondria-endoplasmic reticulum interactions and augmented mitophagy flux which acts as a pro-survival mechanism to eliminate mitochondrial damage. Genetic or pharmacologic MFN2 targeting synergizes with BH3-mimetics by impairing mitochondrial clearance and enhancing apoptosis in AML.

The incidence of melanoma, being one of the most commonly occurring cancers, has been rising since the past decade. Patients at advanced stages of the disease have very poor prognoses, as opposed to at the earlier stages. The conventional targeted therapy is well defined and effective for advanced-stage melanomas for patients not responding to the standard-of-care immunotherapy. However, targeted therapies do not prove to be as effective as patients inevitably develop V-Raf Murine Sarcoma Viral Oncogene Homolog B (BRAF)-inhibitor resistance to the respective drugs. Factors which are driving melanoma drug resistance mainly involve mutations in the mitogen-activated protein kinase (MAPK) pathway, e.g., BRAF splice variants, neuroblastoma RAS viral oncogene homolog (NRAS) amplification or parallel survival pathways. However, those mechanisms do not explain all cases of occurring resistances. Therefore, other factors accounting for BRAFi resistance must be better understood. Among them there are long non-coding RNAs (lncRNAs), but these remain functionally poorly understood. Here, we conduct a comprehensive, unbiased, and integrative study of lncRNA expression, coupled with a Clustered Regularly Interspaced Short Palindromic Repeats/Cas9-mediated activation (CRISPRa) and small molecule inhibitor screening for BRAF inhibitor resistance to expand the knowledge of potentially druggable lncRNAs, their function, and pave the way for eventual combinatorial treatment approaches targeting diverse pathways in melanoma.

Somatostatin interneurons are the earliest born population of cortical inhibitory cells. They are crucial to support normal brain development and function; however, the mechanisms underlying their integration into nascent cortical circuitry are not well understood. In this study, we begin by demonstrating that the maturation of somatostatin interneurons in mouse somatosensory cortex is activity dependent. We then investigated the relationship between activity, alternative splicing, and synapse formation within this population. Specifically, we discovered that the Nova family of RNA-binding proteins are activity-dependent and are essential for the maturation of somatostatin interneurons, as well as their afferent and efferent connectivity. Within this population, Nova2 preferentially mediates the alternative splicing of genes required for axonal formation and synaptic function independently from its effect on gene expression. Hence, our work demonstrates that the Nova family of proteins through alternative splicing are centrally involved in coupling developmental neuronal activity to cortical circuit formation.

Whereas the cellular and molecular features of human inflammatory skin diseases are well characterized, their tissue context and systemic impact remain poorly understood. We thus profiled human psoriasis (PsO) as a prototypic immune-mediated condition with a high predilection for extracutaneous involvement. Spatial transcriptomics (ST) analyses of 25 healthy, active lesion, and clinically uninvolved skin biopsies and integration with public single-cell transcriptomics data revealed marked differences in immune microniches between healthy and inflamed skin. Tissue-scale cartography further identified core disease features across all active lesions, including the emergence of an inflamed suprabasal epidermal state and the presence of B lymphocytes in lesional skin. Both lesional and distal nonlesional samples were stratified by skin disease severity and not by the presence of systemic disease. This segregation was driven by macrophage-, fibroblast-, and lymphatic-enriched spatial regions with gene signatures associated with metabolic dysfunction. Together, these findings suggest that mild and severe forms of PsO have distinct molecular features and that severe PsO may profoundly alter the cellular and metabolic composition of distal unaffected skin sites. In addition, our study provides a valuable resource for the research community to study spatial gene organization of healthy and inflamed human skin.