For a full list of Khavari lab publications please click here


Rubin AJ, Parker KR, Satpathy AT, Qi Y, Wu B, Ong AJ, Mumbach MR, Ji AL, Kim DS, Cho SW, Zarnegar BJ, Greenleaf WJ, Chang HY, Khavari PA. CELL (2019).


Perturb-ATAC, combines multiplexed CRISPRi or knockout with genome-wide chromatin accessibility profiling in single cells via simultaneous detection of sgRNAs and open chromatin sites by assay of transposase-accessible chromatin with sequencing (ATAC-seq). Perturb-ATAC application to transcription factors (TFs), chromatin regulators, & noncoding RNAs (ncRNAs) in ∼4,300 single cells, encompassed more than 63 genotype-phenotype relationships. This uncovered regulators of chromatin accessibility, TF occupancy, and nucleosome positioning and identified a hierarchy of TFs that govern cell state, variation, and disease-associated cis-regulatory elements. Perturb-ATAC in epidermal cells revealed 3 sequential of cis-element modules encoding keratinocyte fate. Combinatorial deletion of all pairs of these TFs uncovered their epistatic relationships and highlighted genomic co-localization as a basis for synergy. Perturb-ATAC is a powerful strategy to dissect gene regulatory networks in development and disease.


Rahmanathan M,* Porter DF,* Khavari PA. NATURE METHODS (2019).


Noncoding RNA sequences, including long noncoding RNAs, small nucleolar RNAs, and untranslated mRNA regions, accomplish many of their diverse functions through direct interactions with RNA-binding proteins (RBPs). The full spectrum of potential RNA-protein interactions within living cells, however, has not yet been defined. This represents a significant gap in knowledge relevant to a host of biological process, including stem cell differentiation and cancer. Recent efforts have identified hundreds of new RBPs that lack known RNA-binding domains, thus underscoring the complexity and diversity of potential RNA-protein complexes. Recent progress has expanded the number of methods for studying RNA-protein interactions in two general categories: approaches that characterize proteins bound to an RNA of interest (RNA-centric), and those that examine RNAs bound to a protein of interest (protein-centric). Each method has unique strengths and limitations, which makes it important to select optimal approaches for the biological question being addressed. This paper reviews methods for the study of RNA-protein interactions, with a focus on their suitability for specific applications.


Kovalski JR, Bhaduri A, Zehnder AM, Neela PH, Che Y, Wozniak GG, Khavari PA. MOLECULAR CELL (2019).


Proximity-dependent biotin labeling (BioID) may identify new targets for cancers driven by difficult-to-drug oncogenes, such as Ras. BioID was therefore used with wild-type (WT) and oncogenic mutant (MT) H-, K-, and N-Ras, identifying known interactors, including Raf and PI3K, as well as a common set of 130 novel proteins proximal to all Ras isoforms. A CRISPR screen of these proteins for Ras-dependence identified mTOR, which was also found proximal to MT Ras in human tumors. Oncogenic Ras directly bound two mTOR Complex 2 (mTORC2) components, mTOR and MAPKAP1, to promote mTORC2 kinase activity at the plasma membrane. mTORC2 enabled the Ras pro-proliferative cell cycle transcriptional program and perturbing the Ras-mTORC2 interaction impaired Ras-dependent neoplasia in vivo. Combining proximity-dependent proteomics with CRISPR screening identified a new set of functional Ras-associated proteins, defined mTORC2 as a new direct Ras effector, and offers a strategy for finding new proteins that cooperate with dominant oncogenes.


Ramanathan M, Maizoub K, Rao DS, Neela PH, Zarnegar BJ, Mondal S, Roth JG, Gai J, Kovalski JR, Siprashvili Z, Palmer TD, Carette JE, Khavari PA. NATURE METHODS (2018).


RNA-protein interaction detection (RAPID) tethers enhanced proximity proteomics to RNA to rapidly identify the proteins that bind any RNA sequence of interest in living cells. RAPID defined protein binding to mutant RNA motifs in human genetic disorders, uncovered potential post-transcriptional networks in breast cancer, and discovered essential host proteins that interact with Zika virus RNA. This new methodology enables direct study of RNA-protein interactions in living cells and tissues on a timescale of minutes to accelerate the growing field of RNA proteomics.

Lineage-specific Dynamic and Pre-established Enhancer-Promoter Contacts Cooperate in Terminal Differentiation.

Rubin AJ, Barajas BC, Furlan-Magaril M, Lopez-Pajares V, Mumbach MR, Howard I, Kim DS, Boxer LD, Cairns J, Spivakov M, Wingett SW, Shi M, Zhao Z, Greenleaf WJ, Kundaje A, Snyder M, Chang HY, Fraser P, Khavari PA. NATURE GENETICS (2017).


A kinetic study of 3D chromatin dynamics across the genome of normal progenitor cells undergoing terminal differentiation identified two types of transcriptional enhancers in contact with target genes.  Constitutive enhancers are pre-looped, H3K27ac-marked, and cohesin-bound while dynamic enhancers gain H3K27ac activation marks and loop to target genes only during differentiation. Distinctive sets of transcription factor (TF) bind each enhancer class and this work discovered EHF as a new essential TF required for constitutive enhancer function in this setting. These new features of genome regulation during stem cell differentiation underscore the diversity of mechanism engaged during this process.


Bao, X, Siprashvili Z, Zarnegar BJ, Shenoy, RM, Rios, EJ, Nady N. Qu, K, Mah, A, Webster, DE, Rubin, AJ, Wozniak, GG, Tao, S, Wysocka, J, Khavari PA. DEVELOPMENTAL CELL (2017).


Screening for histone modifiers essential for progenitor maintenance demonstrated that the PRMT1 histone arginine methyltransferase is required to sustain progenitor cells in the undifferentiated state in vivo in both murine and human epidermis. LC-MS/MS of PRMT1 associated proteins demonstrated that the CSNK1a1 kinase directly binds PRMT1 and controls its genomic targeting. PRMT1 and CSNK1a1 act to cooperatively suppress the pro-differentiation transcription factor, GRHL3, thereby preventing differentiation. These data identify partners and targets of arginine methyltransferases in histone modification and genome during progenitor maintenance.

The noncoding RNAs SNORD50A and SNORD50B bind K-Ras and are recurrently deleted in human cancer.
Siprashvili Z, Webster DE, Johnston D, Shenoy R, Ungewickell A, Bhaduri A, Flockhart R, Zarnegar BJ, Che Y, Meschi F, Puglisi JD & Khavari PA. NATURE GENETICS (2016).

To help define roles for specific small noncoding snoRNAs in cancer, we analyzed 5,473 pairs of tumor and matching normal genomes in 21 human cancer types to discover recurrent deletion of SNORD50A/B in 26% of human tumors. CRISPR knockouts in tumor cells combined with new methodologies to analyze RNA-protein interactions characterized SNORD50A/B as a Ras-binding tumor suppressor RNA that inhibits Ras function by altering its post-translational modification. Small noncoding RNAs therefore can play dominant roles in cancer by binding oncogenes and altering their function.


Zarnegar, BJ, Flynn RA, Shen Y, Do BT, Chang HY & Khavari PA. NATURE METHODS (2016).

The complexity of transcriptome-wide protein-RNA interaction networks is incompletely understood, with current methods resource intensive and technically challenging, irCLIP provides an untraefficient, fast, and nonisotopic method for the detection of protein-RNA interactions.

Network Analysis Identifies Mitochondrial Regulation of Epidermal Differentiation by MPZL3 and FDXR.

Bhaduri A, Ungewickell A, Boxer LD, Lopez-Pajares V, Zarnegar BJ, Khavari PA. DEVELOPMENTAL CELL (2015).

Network reconstruction approaches can help discover new regulators of stem cell differentiation. We generated Proximity Analysis to apply topologically constrained scale-free, small-world biological network construction to discover an essential role for MPZL3, FDXR and reactive oxygen species in epidermal differentiation. 

Genomic analysis of mycosis fungoides and Sézary syndrome identifies recurrent alterations in TNFR2.

Ungewickell A, Bhaduri A, Rios E, Reuter J, Lee CS, Mah A, Zehnder A, Ohgami R, Kulkarni S, Armstrong R, Weng WK, Gratzinger D, Tavallaee M, Rook A, Snyder M, Kim Y & Khavari PA. NATURE GENETICS (2015).

Cancer is characterized by enhanced cellular proliferation in concert with increased resistance to cell death. The mechanisms whereby stem cells undergoing malignant transformation achieve these dual aberrations, however, are highly diverse, depending on tissue and tumor type. Here we use high throughput sequencing to identify recurrent induction of cell survival and activation pathways in cutaneous T cell lymphoma.

CALML5 is a ZNF750- and TINCR-induced protein that binds stratifin to regulate epidermal differentiation.

Sun BK, Boxer LD, Ransohoff JD, Siprashvili Z, Qu K, Lopez-Pajares V, Hollmig ST, Khavari PA. GENES & DEVELOPMENT (2015).


Stem cell differentiation in solid tissue requires spatially coordinated gene regulation. We used LCM RNA-seq to identify CALML5 as the most highly upregulated gene in differentiated epidermis and an essential regulator of this process. The TINCR lncRNA and the ZNF750 TF controlled CALML5 levels . LC-MS/MSidentified SFN as a key CALML5 partner essential for late differentiation. 

A LncRNA-MAF:MAFB transcription factor network regulates epidermal differentiation.

Lopez-Pajares V, Qu K, Zhang J, Webster DE, Barajas BC, Siprashvili Z, Zarnegar BJ, Boxer LD, Rios EJ, Tao S, Kretz M & Khavari PA. DEVELOPMENTAL CELL (2015).

Precisely orchestrated stem cell differentiation is critical for homeostasis in epidermis and other tissues, however, the regulatory networks involved in this process are incompletely understood. Here we identify a new network that is essential for spatially precise genome regulation during epidermal differentiation. At the center of this network are the MAF and MAFB transcription factors (TFs), which are regulated by the TINCR and ANCR lncRNAs to control key effector TFs, including GRHL3, PRDM1, ZNF750 and KLF4.

Advances in skin grafting and treatment of cutaneous wounds.

Sun BK, Siprashvili Z & Khavari PA. SCIENCE (2014).

Tissue regeneration after injury is essential for survival. Advances in stem cell biology, genome editing, and tissue regeneration have laid the foundation for new approaches to cutaneous regeneration and grafting. Cas9 gene-edited somatic cells and iPS cells may be used in this process for the treatment of a host of monogenic disorders, such as epidermolysis bullosa (EB). 

ZNF750 interacts with KLF4 and RCOR1, KDM1A, and CTBP1/2 chromatin regulators to repress epidermal progenitor genes and induce differentiation genes.

Boxer LD, Barajas B, Tao S, Zhang J & Khavari PA. GENES & DEVELOPMENT (2014).

Stem cell differentiation requires down-regulation of stem cell gene expression and induction of terminal differentiation genes. Dominant TFs, such as ZNF750, accomplish both types of impacts, but how they simultaneously effect selective gene repression and activation is unclear. Here we characterize the ZNF750 protein interactome by LC-MS/MS to discover that ZNF750 targets KDM1A to repress stem cell genes while simultaneously co-binding with KLF4 at differentiation genes to activate their expression.

Recurrent point mutations in the kinetochore gene KNSTRN in cutaneous squamous cell carcinoma.

Lee CS, Bhaduri A, Mah A, Johnson WL, Ungewickell A, Aros CJ, Nguyen CB, Rios EJ, Siprashvili Z, Straight A, Kim J, Aasi SZ & Khavari PA. NATURE GENETICS (2014).

Epithelial neoplasms, which comprise ~90% of human cancers, are characterized by widespread genome damage by mechanisms that are incompletely understood. Here we used high throughput genome sequencing to define recurrent mutations in epidermal squamous cell carcinoma (SCC), the second most common cancer in humans. SCC displays hotspot mutations in the KNSTRN gene that trigger aneuploidy and accelerate tumorigenesis, identifying a new common mechanism of genomic injury in cancer.


Khavari Lab

Khavari Lab

IQGAP1 scaffold-kinase interaction blockade selectively targets RAS-MAP kinase-driven tumors.

Jameson KL, Mazur PK, Zehnder AM, Zhang J, Zarnegar B, Sage J & Khavari PA. NATURE MEDICINE (2013).

Hyperactive signaling of the Ras-Erk1/2 MAPK pathway occurs commonly in cancer, however, complete pathway ablation is lethal, underscoring the need for new selective targeting approaches. Scaffold proteins, such as IQGAP1, specify signaling output by assembling pathway components, such as Raf, Mek, and Erk kinases. Here we demonstrate that IQAP1 scaffold-kinase inhibition (SKIB) – by selectively disrupting Erk1/2 pathway assembly in cancer - represents a promising new approach to molecular oncology therapy.

ACTL6a enforces the epidermal progenitor state by suppressing SWI/SNF-dependent induction of KLF4.

Bao X, Tang J, Lopez-Pajares V, Tao S, Qu K, Crabtree GR & Khavari PA. CELL STEM CELL (2013).

Somatic stem cells suppress differentiation to maintain tissue self-renewal. Chromatin remodelers, such as the BAF (SWI/SNF) complex, control nucleosome positioning to modulate genome expression. Here we identify a selective role for the ACTL6A BAF chromatin remodeler subunit in stem cell maintenance. ACTL6A blocks BAF-mediated activation of KLF4, a TF best known for iPS cell induction whose essential role in mammals in vivo involves induction of terminal epidermal differentiation.

Control of somatic tissue differentiation by the long non-coding RNA TINCR.

Kretz M, Siprashvili Z, Chu C, Webster DE, Zehnder A, Qu K, Lee CS, Flockhart RJ, Groff A, Chow J, Johnston D, Kim G, Spitale RC, Flynn RA, Zheng G, Aiyer S, Raj A, Rinn JL, Chang HY & Khavari PA. NATURE (2013).

Long noncoding RNAs (lncRNAs) act by a host of mechanisms to impact genomic expression. To identify lncRNAs with essential roles in differentiation, we undertook a systematic screen to identify Terminal differentiation Induced NonCoding RNA (TINCR).  TINCR controlled expression of hundreds of genes by a novel post-transcriptional mechanism that involved complementary base pairing by TINCR to target mRNAs through a 25nt “TINCR” box sequence.  TINCR stabilized bound differentiation gene mRNAs in concert with the Staufen 1 protein, demonstrating an entirely new mechanism for lncRNA action.