High-throughput single-cell functional elucidation of neurodevelopmental disease-associated genes reveals convergent mechanisms altering neuronal differentiation.

The overwhelming success of exome- and genome-wide association studies in discovering thousands of disease-associated genes necessitates developing novel high-throughput functional genomics approaches to elucidate the molecular mechanisms of these genes. Here, we have coupled multiplexed repression of neurodevelopmental disease-associated genes to single-cell transcriptional profiling in differentiating human neurons to rapidly assay the functions of multiple genes in a disease-relevant context, assess potentially convergent mechanisms, and prioritize genes for specific functional assays. For a set of 13 autism spectrum disorder (ASD)-associated genes, we show that this approach generated important mechanistic insights, revealing two functionally convergent modules of ASD genes: one that delays neuron differentiation and one that accelerates it. Five genes that delay neuron differentiation (ADNP, ARID1B, ASH1L, CHD2, and DYRK1A) mechanistically converge, as they all dysregulate genes involved in cell-cycle control and progenitor cell proliferation. Live-cell imaging after individual ASD-gene repression validated this functional module, confirming that these genes reduce neural progenitor cell proliferation and neurite growth. Finally, these functionally convergent ASD gene modules predicted shared clinical phenotypes among individuals with mutations in these genes. Altogether, these results show the utility of a novel and simple approach for the rapid functional elucidation of neurodevelopmental disease-associated genes.

The spring-loaded genome: nucleosome redistributions are widespread, transient, and DNA-directed.

Nucleosome occupancy plays a key role in regulating access to eukaryotic genomes. Although various chromatin regulatory complexes are known to regulate nucleosome occupancy, the role of DNA sequence in this regulation remains unclear, particularly in mammals. To address this problem, we measured nucleosome distribution at high temporal resolution in human cells at hundreds of genes during the reactivation of Kaposi's sarcoma-associated herpesvirus (KSHV). We show that nucleosome redistribution peaks at 24 h post-KSHV reactivation and that the nucleosomal redistributions are widespread and transient. To clarify the role of DNA sequence in these nucleosomal redistributions, we compared the genes with altered nucleosome distribution to a sequence-based computer model and in vitro-assembled nucleosomes. We demonstrate that both the predicted model and the assembled nucleosome distributions are concordant with the majority of nucleosome redistributions at 24 h post-KSHV reactivation. We suggest a model in which loci are held in an unfavorable chromatin architecture and "spring" to a transient intermediate state directed by DNA sequence information. We propose that DNA sequence plays a more considerable role in the regulation of nucleosome positions than was previously appreciated. The surprising findings that nucleosome redistributions are widespread, transient, and DNA-directed shift the current perspective regarding regulation of nucleosome distribution in humans.

ORF33 and ORF38 of Kaposi's Sarcoma-Associated Herpesvirus Interact and Are Required for Optimal Production of Infectious Progeny Viruses.

UNLABELLED: We recently showed that the interaction between Kaposi's sarcoma-associated herpesvirus (KSHV) tegument proteins ORF33 and ORF45 is crucial for progeny virion production, but the exact functions of KSHV ORF33 during lytic replication were unknown (J. Gillen, W. Li, Q. Liang, D. Avey, J. Wu, F. Wu, J. Myoung, and F. Zhu, J Virol 89:4918-4931, 2015, http://dx.doi.org/10.1128/JVI.02925-14). Therefore, here we investigated the relationship between ORF33 and ORF38, whose counterparts in both alpha- and betaherpesviruses interact with each other. Using specific monoclonal antibodies, we found that both proteins are expressed during the late lytic cycle with similar kinetics and that both are present in mature virions as components of the tegument. Furthermore, we confirmed that ORF33 interacts with ORF38. Interestingly, we observed that ORF33 tightly associates with the capsid, whereas ORF38 associates with the envelope. We generated ORF33-null, ORF38-null, and double-null mutants and found that these mutants apparently have identical phenotypes: the mutations caused no apparent effect on viral gene expression but reduced the yield of progeny virion by about 10-fold. The progeny virions also lack certain virion component proteins, including ORF45. During viral lytic replication, the virions associate with cytoplasmic vesicles. We also observed that ORF38 associates with the membranes of vesicles and colocalizes with the Golgi membrane or early endosome membrane. Further analyses of ORF33/ORF38 mutants revealed the reduced production of virion-containing vesicles, suggesting that ORF33 and ORF38 are involved in the transport of newly assembled viral particles into cytoplasmic vesicles, a process important for viral maturation and egress. IMPORTANCE: Herpesvirus assembly is an essential step in virus propagation that leads to the generation of progeny virions. It is a complicated process that depends on the delicate regulation of interactions among virion proteins. We previously revealed an essential role of ORF45-ORF33 binding for virus assembly. Here, we report that ORF33 and its binding partner, ORF38, are required for infectious virus production due to their important role in the tegumentation process. Moreover, we found that both ORF33 and ORF38 are involved in the transportation of virions through vesicles during maturation and egress. Our results provide new insights into the important roles of ORF33 and ORF38 during viral assembly, a process critical for virus propagation that is intimately linked to KSHV pathobiology.

Discovery of a Coregulatory Interaction between Kaposi's Sarcoma-Associated Herpesvirus ORF45 and the Viral Protein Kinase ORF36.

UNLABELLED: Kaposi's sarcoma-associated herpesvirus (KSHV) is the causative agent of three human malignancies. KSHV ORF36 encodes a serine/threonine viral protein kinase, which is conserved throughout all herpesviruses. Although several studies have identified the viral and cellular substrates of conserved herpesvirus protein kinases (CHPKs), the precise functions of KSHV ORF36 during lytic replication remain elusive. Here, we report that ORF36 interacts with another lytic protein, ORF45, in a manner dependent on ORF36 kinase activity. We mapped the regions of ORF36 and ORF45 involved in the binding. Their association appears to be mediated by electrostatic interactions, since deletion of either the highly basic N terminus of ORF36 or an acidic patch of ORF45 abolished the binding. In addition, the dephosphorylation of ORF45 protein dramatically reduced its association with ORF36. Importantly, ORF45 enhances both the stability and kinase activity of ORF36. Consistent with previous studies of CHPK homologs, we detected ORF36 protein in extracellular virions. To investigate the roles of ORF36 in the context of KSHV lytic replication, we used bacterial artificial chromosome mutagenesis to engineer both ORF36-null and kinase-dead mutants. We found that ORF36-null/mutant virions are moderately defective in viral particle production and are further deficient in primary infection. In summary, our results uncover a functionally important interaction between ORF36 and ORF45 and indicate a significant role of ORF36 in the production of infectious progeny virions. IMPORTANCE: Kaposi's sarcoma-associated herpesvirus (KSHV) is a human tumor virus with a significant public health burden. KSHV ORF36 encodes a serine/threonine viral protein kinase, whose functions throughout the viral life cycle have not been elucidated. Here, we report that ORF36 interacts with another KSHV protein, ORF45. We mapped the regions of ORF36 and ORF45 involved in their association and further characterized the consequences of this interaction. We engineered ORF36 mutant viruses in order to investigate the functional roles of ORF36 in the context of KSHV lytic replication, and we confirmed that ORF36 is a component of KSHV virions. Moreover, we found that ORF36 mutants are defective in virion production and primary infection. In summary, we discovered and characterized a functionally important interaction between KSHV ORF36 and ORF45, and our results suggest a significant role of ORF36 in the production of infectious progeny virions, a process critical for KSHV pathogenesis.

Kaposi's Sarcoma-Associated Herpesvirus Inhibitor of cGAS (KicGAS), Encoded by ORF52, Is an Abundant Tegument Protein and Is Required for Production of Infectious Progeny Viruses.

UNLABELLED: Although Kaposi's sarcoma-associated herpesvirus (KSHV) ORF52 (also known as KSHV inhibitor of cGAS [KicGAS]) has been detected in purified virions, the roles of this protein during KSHV replication have not been characterized. Using specific monoclonal antibodies, we revealed that ORF52 displays true late gene expression kinetics and confirmed its cytoplasmic localization in both transfected and KSHV-infected cells. We demonstrated that ORF52 comigrates with other known virion proteins following sucrose gradient centrifugation. We also determined that ORF52 resides inside the viral envelope and remains partially associated with capsid when extracellular virions are treated with various detergents and/or salts. There results indicate that ORF52 is a tegument protein abundantly present in extracellular virions. To characterize the roles of ORF52 in the KSHV life cycle, we engineered a recombinant KSHV ORF52-null mutant virus and found that loss of ORF52 results in reduced virion production and a further defect in infectivity. Upon analysis of the virion composition of ORF52-null viral particles, we observed a decrease in the incorporation of ORF45, as well as other tegument proteins, suggesting that ORF52 is important for the packaging of other virion proteins. In summary, our results indicate that, in addition to its immune evasion function, KSHV ORF52 is required for the optimal production of infectious virions, likely due to its roles in virion assembly as a tegument protein. IMPORTANCE: The tegument proteins of herpesviruses, including Kaposi's sarcoma-associated herpesvirus (KSHV), play key roles in the viral life cycle. Each of the three subfamilies of herpesviruses (alpha, beta, and gamma) encode unique tegument proteins with specialized functions. We recently found that one such gammaherpesvirus-specific protein, ORF52, has an important role in immune evasion during KSHV primary infection, through inhibition of the host cytosolic DNA sensing pathway. In this report, we further characterize ORF52 as a tegument protein with vital roles during KSHV lytic replication. We found that ORF52 is important for the production of infectious viral particles, likely through its role in virus assembly, a critical process for KSHV replication and pathogenesis. More comprehensive investigation of the functions of tegument proteins and their roles in viral replication may reveal novel targets for therapeutic interventions against KSHV-associated diseases.

Phosphoproteomic Analysis of KSHV-Infected Cells Reveals Roles of ORF45-Activated RSK during Lytic Replication.

Kaposi's Sarcoma-Associated Herpesvirus (KSHV) is an oncogenic virus which has adapted unique mechanisms to modulate the cellular microenvironment of its human host. The pathogenesis of KSHV is intimately linked to its manipulation of cellular signaling pathways, including the extracellular signal-regulated kinase (ERK) mitogen-activated protein kinase (MAPK) pathway. We have previously shown that KSHV ORF45 contributes to the sustained activation of both ERK and p90 ribosomal S6 kinase (RSK, a major functional mediator of ERK/MAPK signaling) during KSHV lytic replication. ORF45-activated RSK is required for optimal KSHV lytic gene expression and progeny virion production, though the underlying mechanisms downstream of this activation are still unclear. We hypothesized that the activation of RSK by ORF45 causes differential phosphorylation of cellular and viral substrates, affecting biological processes essential for efficient KSHV lytic replication. Accordingly, we observed widespread and significant differences in protein phosphorylation upon induction of lytic replication. Mass-spectrometry-based phosphoproteomic screening identified putative substrates of ORF45-activated RSK in KSHV-infected cells. Bioinformatic analyses revealed that nuclear proteins, including several transcriptional regulators, were overrepresented among these candidates. We validated the ORF45/RSK-dependent phosphorylation of several putative substrates by employing KSHV BAC mutagenesis, kinase inhibitor treatments, and/or CRISPR-mediated knockout of RSK in KSHV-infected cells. Furthermore, we assessed the consequences of knocking out these substrates on ORF45/RSK-dependent regulation of gene expression and KSHV progeny virion production. Finally, we show data to support that ORF45 regulates the translational efficiency of a subset of viral/cellular genes with complex secondary structure in their 5' UTR. Altogether, these data shed light on the mechanisms by which KSHV ORF45 manipulates components of the host cell machinery via modulation of RSK activity. Thus, this study has important implications for the pathobiology of KSHV and other diseases in which RSK activity is dysregulated.

ORF45-Mediated Prolonged c-Fos Accumulation Accelerates Viral Transcription during the Late Stage of Lytic Replication of Kaposi's Sarcoma-Associated Herpesvirus.

UNLABELLED: Kaposi's sarcoma-associated herpesvirus (KSHV) encodes multiple viral proteins that activate extracellular signal-regulated kinase (ERK)-mitogen-activated protein kinase (MAPK) cascades. One of these viral proteins, ORF45, mediates sustained ERK-p90 ribosomal S6 kinase (RSK) activation during KSHV lytic replication and facilitates viral translation through the phosphorylation of a eukaryotic translation initiation factor, eIF4B. The importance of ERK-RSK activation for KSHV viral transcription has been shown; however, which transcription factor senses the sustained MAPK signaling and leads to viral transcription remains poorly understood. Here we show that the presence of ORF45 leads to the prolonged accumulation of c-Fos during the late stage of KSHV lytic replication through ERK-RSK-dependent phosphorylation and stabilization and that the depletion of c-Fos disrupts viral lytic transcription. Genome-wide screening revealed that c-Fos directly binds to multiple viral gene promoters and enhances viral transcription. Mutation of the ERK-RSK phosphorylation sites of c-Fos restrains KSHV lytic gene expression and virion production. These results indicate that the prolonged accumulation of c-Fos promotes the progression of viral transcription from early to late stages and accelerates viral lytic replication upon sustained ORF45-ERK-RSK activation during the KSHV lytic life cycle. IMPORTANCE: During KSHV lytic replication, transient activation and sustained activation of ERK-RSK induce viral immediate early (IE) transcription and late transcription, respectively. Studies have revealed that ERK-RSK activates several transcription factors involved in IE gene expression, including Ets, AP-1, CREB, and C/EBP, which lead to the transient ERK-RSK activation-dependent IE transcription. Whereas c-Fos acts as a sensor of sustained ERK-RSK activation, ORF45-ERK-RSK signaling mediates c-Fos phosphorylation and accumulation during late KSHV lytic replication, consequently promoting viral transcription through the direct binding of c-Fos to multiple KSHV promoters. This finding indicates that c-Fos mediates distinct viral transcriptional progression following sustained ERK-RSK signaling during the KSHV lytic life cycle.

A survey of the interactome of Kaposi's sarcoma-associated herpesvirus ORF45 revealed its binding to viral ORF33 and cellular USP7, resulting in stabilization of ORF33 that is required for production of progeny viruses.

UNLABELLED: The ORF45 protein of Kaposi's sarcoma-associated herpesvirus (KSHV) is a gammaherpesvirus-specific immediate-early tegument protein. Our previous studies have revealed its crucial roles in both early and late stages of KSHV infection. In this study, we surveyed the interactome of ORF45 using a panel of monoclonal antibodies. In addition to the previously identified extracellular regulated kinase (ERK) and p90 ribosomal S6 kinase (RSK) proteins, we found several other copurified proteins, including prominent ones of ∼38 kDa and ∼130 kDa. Mass spectrometry revealed that the 38-kDa protein is viral ORF33 and the 130-kDa protein is cellular USP7 (ubiquitin-specific protease 7). We mapped the ORF33-binding domain to the highly conserved carboxyl-terminal 19 amino acids (aa) of ORF45 and the USP7-binding domain to the reported consensus motif in the central region of ORF45. Using immunofluorescence staining, we observed colocalization of ORF45 with ORF33 or USP7 both under transfected conditions and in KSHV-infected cells. Moreover, we noticed ORF45-dependent relocalization of a portion of ORF33/USP7 from the nucleus to the cytoplasm. We found that ORF45 caused an increase in ORF33 protein accumulation that was abolished if either the ORF33- or USP7-binding domain in ORF45 was deleted. Furthermore, deletion of the conserved carboxyl terminus of ORF45 in the KSHV genome drastically reduced the level of ORF33 protein in KSHV-infected cells and abolished production of progeny virions. Collectively, our results not only reveal new components of the ORF45 interactome, but also demonstrate that the interactions among these proteins are crucial for KSHV lytic replication. IMPORTANCE: Kaposi's sarcoma-associated herpesvirus (KSHV) is the causative agent of several human cancers. KSHV ORF45 is a multifunctional protein that is required for KSHV lytic replication, but the exact mechanisms by which ORF45 performs its critical functions are unclear. Our previous studies revealed that all ORF45 protein in cells exists in high-molecular-weight complexes. We therefore sought to characterize the interactome of ORF45 to provide insights into its roles during lytic replication. Using a panel of monoclonal antibodies, we surveyed the ORF45 interactome in KSHV-infected cells. We identified two new binding partners of ORF45: the viral protein ORF33 and cellular ubiquitin-specific protease 7 (USP7). We further demonstrate that the interaction between ORF45 and ORF33 is crucial for the efficient production of KSHV viral particles, suggesting that the targeted interference with this interaction may represent a novel strategy to inhibit KSHV lytic replication.

Activation of p90 ribosomal S6 kinases by ORF45 of Kaposi's sarcoma-associated herpesvirus is critical for optimal production of infectious viruses.

UNLABELLED: We have previously shown that ORF45, an immediate-early and tegument protein of Kaposi's sarcoma-associated herpesvirus (KSHV), causes sustained activation of p90 ribosomal S6 kinases (RSKs) and extracellular regulated kinase (ERK) (E. Kuang, Q. Tang, G. G. Maul, and F. Zhu, J Virol 82:1838-1850, 2008, http://dx.doi.org/10.1128/JVI.02119-07). We now have identified the critical region of ORF45 that is involved in RSK interaction and activation. Alanine scanning mutagenesis of this region revealed that a single F66A point mutation abolished binding of ORF45 to RSK or ERK and, consequently, its ability to activate the kinases. We introduced the F66A mutation into BAC16 (a bacterial artificial chromosome clone containing the entire infectious KSHV genome), producing BAC16-45F66A. In parallel, we also repaired the mutation and obtained a revertant, BAC16-45A66F. The reconstitution of these mutants in iSLK cells demonstrated that the ORF45-F66A mutant failed to cause sustained ERK and RSK activation during lytic reactivation, resulting in dramatic differences in the phosphoproteomic profile between the wild-type virus-infected cells and the mutant virus-infected cells. ORF45 mutation or deletion also was accompanied by a noticeable decreased in viral gene expression during lytic reactivation. Consequently, the ORF45-F66A mutant produced significantly fewer infectious progeny virions than the wild type or the revertant. These results suggest a critical role for ORF45-mediated RSK activation in KSHV lytic replication. IMPORTANCE: KSHV is the causative agent of three human malignancies. KSHV pathogenesis is intimately linked to its ability to modulate the host cell microenvironment and to facilitate efficient production of progeny viral particles. We previously described the mechanism by which the KSHV lytic protein ORF45 activates the cellular kinases ERK and RSK. We now have mapped the critical region of ORF45 responsible for binding and activation of ERK/RSK to a single residue, F66. We mutated this amino acid of ORF45 (F66A) and introduced the mutation into a newly developed bacterial artificial chromosome containing the KSHV genome (BAC16). This system has provided us with a useful tool to characterize the functions of ORF45-activated RSK upon KSHV lytic reactivation. We show that viral gene expression and virion production are significantly reduced by F66A mutation, indicating a critical role for ORF45-activated RSK during KSHV lytic replication.

SR-TWAS: Leveraging Multiple Reference Panels to Improve TWAS Power by Ensemble Machine Learning.

Multiple reference panels of a given tissue or multiple tissues often exist, and multiple regression methods could be used for training gene expression imputation models for TWAS. To leverage expression imputation models (i.e., base models) trained with multiple reference panels, regression methods, and tissues, we develop a Stacked Regression based TWAS (SR-TWAS) tool which can obtain optimal linear combinations of base models for a given validation transcriptomic dataset. Both simulation and real studies showed that SR-TWAS improved power, due to increased effective training sample sizes and borrowed strength across multiple regression methods and tissues. Leveraging base models across multiple reference panels, tissues, and regression methods, our real application studies identified 6 independent significant risk genes for Alzheimer's disease (AD) dementia for supplementary motor area tissue and 9 independent significant risk genes for Parkinson's disease (PD) for substantia nigra tissue. Relevant biological interpretations were found for these significant risk genes.

Uncovering Plaque-Glia Niches in Human Alzheimer's Disease Brains Using Spatial Transcriptomics.

Amyloid-beta (Aß) plaques and surrounding glial activation are prominent histopathological hallmarks of Alzheimer's Disease (AD). However, it is unclear how Aß plaques interact with surrounding glial cells in the human brain. Here, we applied spatial transcriptomics (ST) and immunohistochemistry (IHC) for Aß, GFAP, and IBA1 to acquire data from 258,987 ST spots within 78 postmortem brain sections of 21 individuals. By coupling ST and adjacent-section IHC, we showed that low Aß spots exhibit transcriptomic profiles indicative of greater neuronal loss than high Aß spots, and high-glia spots present transcriptomic changes indicative of more significant inflammation and neurodegeneration. Furthermore, we observed that this ST glial response bears signatures of reported mouse gene modules of plaque-induced genes (PIG), oligodendrocyte (OLIG) response, disease-associated microglia (DAM), and disease-associated astrocytes (DAA), as well as different microglia (MG) states identified in human AD brains, indicating that multiple glial cell states arise around plaques and contribute to local immune response. We then validated the observed effects of Aß on cell apoptosis and plaque-surrounding glia on inflammation and synaptic loss using IHC. In addition, transcriptomic changes of iPSC-derived microglia-like cells upon short-interval Aß treatment mimic the ST glial response and mirror the reported activated MG states. Our results demonstrate an exacerbation of synaptic and neuronal loss in low-Aß or high-glia areas, indicating that microglia response to Aß-oligomers likely initiates glial activation in plaque-glia niches. Our study lays the groundwork for future pathology genomics studies, opening the door for investigating pathological heterogeneity and causal effects in neurodegenerative diseases.

SR-TWAS: leveraging multiple reference panels to improve transcriptome-wide association study power by ensemble machine learning.

Multiple reference panels of a given tissue or multiple tissues often exist, and multiple regression methods could be used for training gene expression imputation models for transcriptome-wide association studies (TWAS). To leverage expression imputation models (i.e., base models) trained with multiple reference panels, regression methods, and tissues, we develop a Stacked Regression based TWAS (SR-TWAS) tool which can obtain optimal linear combinations of base models for a given validation transcriptomic dataset. Both simulation and real studies show that SR-TWAS improves power, due to increased training sample sizes and borrowed strength across multiple regression methods and tissues. Leveraging base models across multiple reference panels, tissues, and regression methods, our real studies identify 6 independent significant risk genes for Alzheimer's disease (AD) dementia for supplementary motor area tissue and 9 independent significant risk genes for Parkinson's disease (PD) for substantia nigra tissue. Relevant biological interpretations are found for these significant risk genes.

Differential gene expression analysis based on linear mixed model corrects false positive inflation for studying quantitative traits.

Differential gene expression (DGE) analysis has been widely employed to identify genes expressed differentially with respect to a trait of interest using RNA sequencing (RNA-Seq) data. Recent RNA-Seq data with large samples pose challenges to existing DGE methods, which were mainly developed for dichotomous traits and small sample sizes. Especially, existing DGE methods are likely to result in inflated false positive rates. To address this gap, we employed a linear mixed model (LMM) that has been widely used in genetic association studies for DGE analysis of quantitative traits. We first applied the LMM method to the discovery RNA-Seq data of dorsolateral prefrontal cortex (DLPFC) tissue (n = 632) with four continuous measures of Alzheimer's Disease (AD) cognitive and neuropathologic traits. The quantile-quantile plots of p-values showed that false positive rates were well calibrated by LMM, whereas other methods not accounting for sample-specific mixed effects led to serious inflation. LMM identified 37 potentially significant genes with differential expression in DLPFC for at least one of the AD traits, 17 of which were replicated in the additional RNA-Seq data of DLPFC, supplemental motor area, spinal cord, and muscle tissues. This application study showed not only well calibrated DGE results by LMM, but also possibly shared gene regulatory mechanisms of AD traits across different relevant tissues.

Inferring protein expression changes from mRNA in Alzheimer's dementia using deep neural networks.

Identifying the molecular systems and proteins that modify the progression of Alzheimer's disease and related dementias (ADRD) is central to drug target selection. However, discordance between mRNA and protein abundance, and the scarcity of proteomic data, has limited our ability to advance candidate targets that are mainly based on gene expression. Therefore, by using a deep neural network that predicts protein abundance from mRNA expression, here we attempt to track the early protein drivers of ADRD. Specifically, by applying the clei2block deep learning model to 1192 brain RNA-seq samples, we identify protein modules and disease-associated expression changes that were not directly observed at the mRNA level. Moreover, pseudo-temporal trajectory inference based on the predicted proteome became more closely correlated with cognitive decline and hippocampal atrophy compared to RNA-based trajectories. This suggests that the predicted changes in protein expression could provide a better molecular representation of ADRD progression. Furthermore, overlaying clinical traits on protein pseudotime trajectory identifies protein modules altered before cognitive impairment. These results demonstrate how our method can be used to identify potential early protein drivers and possible drug targets for treating and/or preventing ADRD.

Peripheral nerve resident macrophages share tissue-specific programming and features of activated microglia.

Whereas microglia are recognized as fundamental players in central nervous system (CNS) development and function, much less is known about macrophages of the peripheral nervous system (PNS). Here, by comparing gene expression across neural and conventional tissue-resident macrophages, we identified transcripts that were shared among neural resident macrophages as well as selectively enriched in PNS macrophages. Remarkably, PNS macrophages constitutively expressed genes previously identified to be upregulated by activated microglia during aging, neurodegeneration, or loss of Sall1. Several microglial activation-associated and PNS macrophage-enriched genes were also expressed in spinal cord microglia at steady state. We further show that PNS macrophages rely on IL-34 for maintenance and arise from both embryonic and hematopoietic precursors, while their expression of activation-associated genes did not differ by ontogeny. Collectively, these data uncover shared and unique features between neural resident macrophages and emphasize the role of nerve environment for shaping PNS macrophage identity.

Single-Cell RNA-Seq Uncovers a Robust Transcriptional Response to Morphine by Glia.

Molecular and behavioral responses to opioids are thought to be primarily mediated by neurons, although there is accumulating evidence that other cell types play a prominent role in drug addiction. To investigate cell-type-specific opioid responses, we performed single-cell RNA sequencing (scRNA-seq) of the nucleus accumbens of mice following acute morphine treatment. Differential expression analysis uncovered unique morphine-dependent transcriptional responses by oligodendrocytes and astrocytes. We examined the expression of selected genes, including Cdkn1a and Sgk1, by FISH, confirming their induction by morphine in oligodendrocytes. Further analysis using RNA-seq of FACS-purified oligodendrocytes revealed a large cohort of morphine-regulated genes. The affected genes are enriched for roles in cellular pathways intimately linked to oligodendrocyte maturation and myelination, including the unfolded protein response. Altogether, our data illuminate the morphine-dependent transcriptional response by oligodendrocytes and offer mechanistic insights into myelination defects associated with opioid abuse.

Hierarchical regulation of the genome: global changes in nucleosome organization potentiate genome response.

Nucleosome occupancy is critically important in regulating access to the eukaryotic genome. Few studies in human cells have measured genome-wide nucleosome distributions at high temporal resolution during a response to a common stimulus. We measured nucleosome distributions at high temporal resolution following Kaposi's-sarcoma-associated herpesvirus (KSHV) reactivation using our newly developed mTSS-seq technology, which maps nucleosome distribution at the transcription start sites (TSS) of all human genes. Nucleosomes underwent widespread changes in organization 24 hours after KSHV reactivation and returned to their basal nucleosomal architecture 48 hours after KSHV reactivation. The widespread changes consisted of an indiscriminate remodeling event resulting in the loss of nucleosome rotational phasing signals. Additionally, one in six TSSs in the human genome possessed nucleosomes that are translationally remodeled. 72% of the loci with translationally remodeled nucleosomes have nucleosomes that moved to positions encoded by the underlying DNA sequence. Finally we demonstrated that these widespread alterations in nucleosomal architecture potentiated regulatory factor binding. These descriptions of nucleosomal architecture changes provide a new framework for understanding the role of chromatin in the genomic response, and have allowed us to propose a hierarchical model for chromatin-based regulation of genome response.

Recent advances in the study of Kaposi's sarcoma-associated herpesvirus replication and pathogenesis.

It has now been over twenty years since a novel herpesviral genome was identified in Kaposi's sarcoma biopsies. Since then, the cumulative research effort by molecular biologists, virologists, clinicians, and epidemiologists alike has led to the extensive characterization of this tumor virus, Kaposi's sarcoma-associated herpesvirus (KSHV; also known as human herpesvirus 8 (HHV-8)), and its associated diseases. Here we review the current knowledge of KSHV biology and pathogenesis, with a particular emphasis on new and exciting advances in the field of epigenetics. We also discuss the development and practicality of various cell culture and animal model systems to study KSHV replication and pathogenesis.