Single-Cell Transcriptomic Analysis Links Nonmyelinating Schwann Cells to Proinflammatory Response in the Lung.

The lung is a barrier tissue with constant exposure to the inhaled environment. Therefore, innate immunity against particulates and pathogens is of critical importance to maintain tissue homeostasis. Although the lung harbors both myelinating and nonmyelinating Schwann cells (NMSCs), NMSCs represent the most abundant Schwann cell (SC) population in the lung. However, their contribution to lung physiology remains largely unknown. In this study, we used the human glial fibrillary acidic protein promoter driving tdTomato expression in mice to identify SCs in the peripheral nervous system and determine their location within the lung. Single-cell transcriptomic analysis revealed the existence of two NMSC populations (NMSC1 and NMSC2) that may participate in pathogen recognition. We demonstrated that these pulmonary SCs produce chemokines and cytokines upon LPS stimulation using in vitro conditions. Furthermore, we challenged mouse lungs with LPS and found that NMSC1 exhibits an enriched proinflammatory response among all SC subtypes. Collectively, these findings define the molecular profiles of lung SCs and suggest a potential role for NMSCs in lung inflammation.

MicrobiomeExplorer: an R package for the analysis and visualization of microbial communities.

SUMMARY: We developed the MicrobiomeExplorer R package to facilitate the analysis and visualization of microbial communities. The MicrobiomeExplorer R package allows a user to perform typical microbiome analytic workflows and visualize their results, either through the command line or an interactive Shiny application included with the package. In addition to applying common analytical workflows, the application enables automated analysis report generation. AVAILABILITY AND IMPLEMENTATION: Available at https://github.com/zoecastillo/microbiomeExplorer. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A resource for cell line authentication, annotation and quality control.

Cell line misidentification, contamination and poor annotation affect scientific reproducibility. Here we outline simple measures to detect or avoid cross-contamination, present a framework for cell line annotation linked to short tandem repeat and single nucleotide polymorphism profiles, and provide a catalogue of synonymous cell lines. This resource will enable our community to eradicate the use of misidentified lines and generate credible cell-based data.

Complex regulation of ADAR-mediated RNA-editing across tissues.

BACKGROUND: RNA-editing is a tightly regulated, and essential cellular process for a properly functioning brain. Dysfunction of A-to-I RNA editing can have catastrophic effects, particularly in the central nervous system. Thus, understanding how the process of RNA-editing is regulated has important implications for human health. However, at present, very little is known about the regulation of editing across tissues, and individuals. RESULTS: Here we present an analysis of RNA-editing patterns from 9 different tissues harvested from a single mouse. For comparison, we also analyzed data for 5 of these tissues harvested from 15 additional animals. We find that tissue specificity of editing largely reflects differential expression of substrate transcripts across tissues. We identified a surprising enrichment of editing in intronic regions of brain transcripts, that could account for previously reported higher levels of editing in brain. There exists a small but remarkable amount of editing which is tissue-specific, despite comparable expression levels of the edit site across multiple tissues. Expression levels of editing enzymes and their isoforms can explain some, but not all of this variation. CONCLUSIONS: Together, these data suggest a complex regulation of the RNA-editing process beyond transcript expression levels.

Diverse somatic mutation patterns and pathway alterations in human cancers.

The systematic characterization of somatic mutations in cancer genomes is essential for understanding the disease and for developing targeted therapeutics. Here we report the identification of 2,576 somatic mutations across approximately 1,800 megabases of DNA representing 1,507 coding genes from 441 tumours comprising breast, lung, ovarian and prostate cancer types and subtypes. We found that mutation rates and the sets of mutated genes varied substantially across tumour types and subtypes. Statistical analysis identified 77 significantly mutated genes including protein kinases, G-protein-coupled receptors such as GRM8, BAI3, AGTRL1 (also called APLNR) and LPHN3, and other druggable targets. Integrated analysis of somatic mutations and copy number alterations identified another 35 significantly altered genes including GNAS, indicating an expanded role for galpha subunits in multiple cancer types. Furthermore, our experimental analyses demonstrate the functional roles of mutant GNAO1 (a Galpha subunit) and mutant MAP2K4 (a member of the JNK signalling pathway) in oncogenesis. Our study provides an overview of the mutational spectra across major human cancers and identifies several potential therapeutic targets.

DLK initiates a transcriptional program that couples apoptotic and regenerative responses to axonal injury.

The cell intrinsic factors that determine whether a neuron regenerates or undergoes apoptosis in response to axonal injury are not well defined. Here we show that the mixed-lineage dual leucine zipper kinase (DLK) is an essential upstream mediator of both of these divergent outcomes in the same cell type. Optic nerve crush injury leads to rapid elevation of DLK protein, first in the axons of retinal ganglion cells (RGCs) and then in their cell bodies. DLK is required for the majority of gene expression changes in RGCs initiated by injury, including induction of both proapoptotic and regeneration-associated genes. Deletion of DLK in retina results in robust and sustained protection of RGCs from degeneration after optic nerve injury. Despite this improved survival, the number of axons that regrow beyond the injury site is substantially reduced, even when the tumor suppressor phosphatase and tensin homolog (PTEN) is deleted to enhance intrinsic growth potential. These findings demonstrate that these seemingly contradictory responses to injury are mechanistically coupled through a DLK-based damage detection mechanism.

Regulation of neuronal gene expression and survival by basal NMDA receptor activity: a role for histone deacetylase 4.

Neuronal gene expression is modulated by activity via calcium-permeable receptors such as NMDA receptors (NMDARs). While gene expression changes downstream of evoked NMDAR activity have been well studied, much less is known about gene expression changes that occur under conditions of basal neuronal activity. In mouse dissociated hippocampal neuronal cultures, we found that a broad NMDAR antagonist, AP5, induced robust gene expression changes under basal activity, but subtype-specific antagonists did not. While some of the gene expression changes are also known to be downstream of stimulated NMDAR activity, others appear specific to basal NMDAR activity. The genes altered by AP5 treatment of basal cultures were enriched for pathways related to class IIa histone deacetylases (HDACs), apoptosis, and synapse-related signaling. Specifically, AP5 altered the expression of all three class IIa HDACs that are highly expressed in the brain, HDAC4, HDAC5, and HDAC9, and also induced nuclear accumulation of HDAC4. HDAC4 knockdown abolished a subset of the gene expression changes induced by AP5, and led to neuronal death under long-term tetrodotoxin or AP5 treatment in rat hippocampal organotypic slice cultures. These data suggest that basal, but not evoked, NMDAR activity regulates gene expression in part through HDAC4, and, that HDAC4 has neuroprotective functions under conditions of low NMDAR activity.

Histone deacetylase (HDAC) inhibitor kinetic rate constants correlate with cellular histone acetylation but not transcription and cell viability.

Histone deacetylases (HDACs) are critical in the control of gene expression, and dysregulation of their activity has been implicated in a broad range of diseases, including cancer, cardiovascular, and neurological diseases. HDAC inhibitors (HDACi) employing different zinc chelating functionalities such as hydroxamic acids and benzamides have shown promising results in cancer therapy. Although it has also been suggested that HDACi with increased isozyme selectivity and potency may broaden their clinical utility and minimize side effects, the translation of this idea to the clinic remains to be investigated. Moreover, a detailed understanding of how HDACi with different pharmacological properties affect biological functions in vitro and in vivo is still missing. Here, we show that a panel of benzamide-containing HDACi are slow tight-binding inhibitors with long residence times unlike the hydroxamate-containing HDACi vorinostat and trichostatin-A. Characterization of changes in H2BK5 and H4K14 acetylation following HDACi treatment in the neuroblastoma cell line SH-SY5Y revealed that the timing and magnitude of histone acetylation mirrored both the association and dissociation kinetic rates of the inhibitors. In contrast, cell viability and microarray gene expression analysis indicated that cell death induction and changes in transcriptional regulation do not correlate with the dissociation kinetic rates of the HDACi. Therefore, our study suggests that determining how the selective and kinetic inhibition properties of HDACi affect cell function will help to evaluate their therapeutic utility.

ReportingTools: an automated result processing and presentation toolkit for high-throughput genomic analyses.

UNLABELLED: It is common for computational analyses to generate large amounts of complex data that are difficult to process and share with collaborators. Standard methods are needed to transform such data into a more useful and intuitive format. We present ReportingTools, a Bioconductor package, that automatically recognizes and transforms the output of many common Bioconductor packages into rich, interactive, HTML-based reports. Reports are not generic, but have been individually designed to reflect content specific to the result type detected. Tabular output included in reports is sortable, filterable and searchable and contains context-relevant hyperlinks to external databases. Additionally, in-line graphics have been developed for specific analysis types and are embedded by default within table rows, providing a useful visual summary of underlying raw data. ReportingTools is highly flexible and reports can be easily customized for specific applications using the well-defined API. AVAILABILITY: The ReportingTools package is implemented in R and available from Bioconductor (version ≥ 2.11) at the URL: http://bioconductor.org/packages/release/bioc/html/ReportingTools.html. Installation instructions and usage documentation can also be found at the above URL.

Spatially coordinated kinase signaling regulates local axon degeneration.

In addition to being a hallmark of neurodegenerative disease, axon degeneration is used during development of the nervous system to prune unwanted connections. In development, axon degeneration is tightly regulated both temporally and spatially. Here, we provide evidence that degeneration cues are transduced through various kinase pathways functioning in spatially distinct compartments to regulate axon degeneration. Intriguingly, glycogen synthase kinase-3 (GSK3) acts centrally, likely modulating gene expression in the cell body to regulate distally restricted axon degeneration. Through a combination of genetic and pharmacological manipulations, including the generation of an analog-sensitive kinase allele mutant mouse for GSK3ß, we show that the ß isoform of GSK3, not the α isoform, is essential for developmental axon pruning in vitro and in vivo. Additionally, we identify the dleu2/mir15a/16-1 cluster, previously characterized as a regulator of B-cell proliferation, and the transcription factor tbx6, as likely downstream effectors of GSK3ß in axon degeneration.

Granulocyte-colony stimulating factor promotes lung metastasis through mobilization of Ly6G+Ly6C+ granulocytes.

Priming of the organ-specific premetastatic sites is thought to be an important yet incompletely understood step during metastasis. In this study, we show that the metastatic tumors we examined overexpress granulocyte-colony stimulating factor (G-CSF), which expands and mobilizes Ly6G+Ly6C+ granulocytes and facilitates their subsequent homing at distant organs even before the arrival of tumor cells. Moreover, G-CSF-mobilized Ly6G+Ly6C+ cells produce the Bv8 protein, which has been implicated in angiogenesis and mobilization of myeloid cells. Anti-G-CSF or anti-Bv8 antibodies significantly reduced lung metastasis. Transplantation of Bv8 null fetal liver cells into lethally irradiated hosts also reduced metastasis. We identified an unexpected role for Bv8: the ability to stimulate tumor cell migration through activation of one of the Bv8 receptors, prokineticin receptor (PKR)-1. Finally, we show that administration of recombinant G-CSF is sufficient to increase the numbers of Ly6G+Ly6C+ cells in organ-specific metastatic sites and results in enhanced metastatic ability of several tumors.

Cross-species transcriptomic atlas of dorsal root ganglia reveals species-specific programs for sensory function.

Sensory neurons of the dorsal root ganglion (DRG) are critical for maintaining tissue homeostasis by sensing and initiating responses to stimuli. While most preclinical studies of DRGs are conducted in rodents, much less is known about the mechanisms of sensory perception in primates. We generated a transcriptome atlas of mouse, guinea pig, cynomolgus monkey, and human DRGs by implementing a common laboratory workflow and multiple data-integration approaches to generate high-resolution cross-species mappings of sensory neuron subtypes. Using our atlas, we identified conserved core modules highlighting subtype-specific biological processes related to inflammatory response. We also identified divergent expression of key genes involved in DRG function, suggesting species-specific adaptations specifically in nociceptors that likely point to divergent function of nociceptors. Among these, we validated that TAFA4, a member of the druggable genome, was expressed in distinct populations of DRG neurons across species, highlighting species-specific programs that are critical for therapeutic development.

Whole genome sequencing across clinical trials identifies rare coding variants in GPR68 associated with chemotherapy-induced peripheral neuropathy.

BACKGROUND: Dose-limiting toxicities significantly impact the benefit/risk profile of many drugs. Whole genome sequencing (WGS) in patients receiving drugs with dose-limiting toxicities can identify therapeutic hypotheses to prevent these toxicities. Chemotherapy-induced peripheral neuropathy (CIPN) is a common dose-limiting neurological toxicity of chemotherapies with no effective approach for prevention. METHODS: We conducted a genetic study of time-to-first peripheral neuropathy event using 30× germline WGS data from whole blood samples from 4900 European-ancestry cancer patients in 14 randomized controlled trials. A substantial number of patients in these trials received taxane and platinum-based chemotherapies as part of their treatment regimen, either standard of care or in combination with the PD-L1 inhibitor atezolizumab. The trials spanned several cancers including renal cell carcinoma, triple negative breast cancer, non-small cell lung cancer, small cell lung cancer, bladder cancer, ovarian cancer, and melanoma. RESULTS: We identified a locus consisting of low-frequency variants in intron 13 of GRID2 associated with time-to-onset of first peripheral neuropathy (PN) indexed by rs17020773 (p = 2.03 × 10-8, all patients, p = 6.36 × 10-9, taxane treated). Gene-level burden analysis identified rare coding variants associated with increased PN risk in the C-terminus of GPR68 (p = 1.59 × 10-6, all patients, p = 3.47 × 10-8, taxane treated), a pH-sensitive G-protein coupled receptor (GPCR). The variants driving this signal were found to alter predicted arrestin binding motifs in the C-terminus of GPR68. Analysis of snRNA-seq from human dorsal root ganglia (DRG) indicated that expression of GPR68 was highest in mechano-thermo-sensitive nociceptors. CONCLUSIONS: Our genetic study provides insight into the impact of low-frequency and rare coding genetic variation on PN risk and suggests that further study of GPR68 in sensory neurons may yield a therapeutic hypothesis for prevention of CIPN.

Nav1.7 is essential for nociceptor action potentials in the mouse in a manner independent of endogenous opioids.

Loss-of-function mutations in Nav1.7, a voltage-gated sodium channel, cause congenital insensitivity to pain (CIP) in humans, demonstrating that Nav1.7 is essential for the perception of pain. However, the mechanism by which loss of Nav1.7 results in insensitivity to pain is not entirely clear. It has been suggested that loss of Nav1.7 induces overexpression of enkephalin, an endogenous opioid receptor agonist, leading to opioid-dependent analgesia. Using behavioral pharmacology and single-cell RNA-seq analysis, we find that overexpression of enkephalin occurs only in cLTMR neurons, a subclass of sensory neurons involved in low-threshold touch detection, and that this overexpression does not play a role in the analgesia observed following genetic removal of Nav1.7. Furthermore, we demonstrate using laser speckle contrast imaging (LSCI) and in vivo electrophysiology that Nav1.7 function is required for the initiation of C-fiber action potentials (APs), which explains the observed insensitivity to pain following genetic removal or inhibition of Nav1.7.

A systems biology approach uncovers novel disease mechanisms in age-related macular degeneration.

Age-related macular degeneration (AMD) is a leading cause of blindness, affecting 200 million people worldwide. To identify genes that could be targeted for treatment, we created a molecular atlas at different stages of AMD. Our resource is comprised of RNA sequencing (RNA-seq) and DNA methylation microarrays from bulk macular retinal pigment epithelium (RPE)/choroid of clinically phenotyped normal and AMD donor eyes (n = 85), single-nucleus RNA-seq (164,399 cells), and single-nucleus assay for transposase-accessible chromatin (ATAC)-seq (125,822 cells) from the retina, RPE, and choroid of 6 AMD and 7 control donors. We identified 23 genome-wide significant loci differentially methylated in AMD, over 1,000 differentially expressed genes across different disease stages, and an AMD Müller state distinct from normal or gliosis. Chromatin accessibility peaks in genome-wide association study (GWAS) loci revealed putative causal genes for AMD, including HTRA1 and C6orf223. Our systems biology approach uncovered molecular mechanisms underlying AMD, including regulators of WNT signaling, FRZB and TLE2, as mechanistic players in disease.

Microarray analysis of retinal endothelial tip cells identifies CXCR4 as a mediator of tip cell morphology and branching.

The development of the vertebrate vascular system is mediated by both genetic patterning of vessels and by angiogenic sprouting in response to hypoxia. Both of these processes depend on the detection of environmental guidance cues by endothelial cells. A specialized subtype of endothelial cell known as the tip cell is thought to be involved in the detection and response to these cues, but the molecular signaling pathways used by tip cells to mediate tissue vascularization remain largely uncharacterized. To identify genes critical to tip cell function, we have developed a method to isolate them using laser capture microdissection, permitting comparison of RNA extracted from endothelial tip cells with that of endothelial stalk cells using microarray analysis. Genes enriched in tip cells include ESM-1, angiopoietin-2, and SLP-76. CXCR4, a receptor for the chemokine stromal-cell derived factor-1, was also identified as a tip cell-enriched gene, and we provide evidence for a novel role for this receptor in mediating tip cell morphology and vascular patterning in the neonatal retina.

Identification and functional analysis of endothelial tip cell-enriched genes.

Sprouting of developing blood vessels is mediated by specialized motile endothelial cells localized at the tips of growing capillaries. Following behind the tip cells, endothelial stalk cells form the capillary lumen and proliferate. Expression of the Notch ligand Delta-like-4 (Dll4) in tip cells suppresses tip cell fate in neighboring stalk cells via Notch signaling. In DLL4(+/-) mouse mutants, most retinal endothelial cells display morphologic features of tip cells. We hypothesized that these mouse mutants could be used to isolate tip cells and so to determine their genetic repertoire. Using transcriptome analysis of retinal endothelial cells isolated from DLL4(+/-) and wild-type mice, we identified 3 clusters of tip cell-enriched genes, encoding extracellular matrix degrading enzymes, basement membrane components, and secreted molecules. Secreted molecules endothelial-specific molecule 1, angiopoietin 2, and apelin bind to cognate receptors on endothelial stalk cells. Knockout mice and zebrafish morpholino knockdown of apelin showed delayed angiogenesis and reduced proliferation of stalk cells expressing the apelin receptor APJ. Thus, tip cells may regulate angiogenesis via matrix remodeling, production of basement membrane, and release of secreted molecules, some of which regulate stalk cell behavior.

NMDA receptor-dependent prostaglandin-endoperoxide synthase 2 induction in neurons promotes glial proliferation during brain development and injury.

Astrocytes play critical roles in brain development and disease, but the mechanisms that regulate astrocyte proliferation are poorly understood. We report that astrocyte proliferation is bi-directionally regulated by neuronal activity via NMDA receptor (NMDAR) signaling in neurons. Prolonged treatment with an NMDAR antagonist reduced expression of cell-cycle-related genes in astrocytes in hippocampal cultures and suppressed astrocyte proliferation in vitro and in vivo, whereas neuronal activation promoted astrocyte proliferation, dependent on neuronal NMDARs. Expression of prostaglandin-endoperoxide synthase 2 (Ptgs2) is induced specifically in neurons by NMDAR activation and is required for activity-dependent astrocyte proliferation through its product, prostaglandin E2 (PGE2). NMDAR inhibition or Ptgs2 genetic ablation in mice reduced the proliferation of astrocytes and microglia induced by mild traumatic brain injury in the absence of secondary excitotoxicity-induced neuronal death. Our study defines an NMDAR-mediated signaling mechanism that allows trans-cellular control of glial proliferation by neurons in brain development and injury.

Disease-associated oligodendrocyte responses across neurodegenerative diseases.

Oligodendrocyte dysfunction has been implicated in the pathogenesis of neurodegenerative diseases, so understanding oligodendrocyte activation states would shed light on disease processes. We identify three distinct activation states of oligodendrocytes from single-cell RNA sequencing (RNA-seq) of mouse models of Alzheimer's disease (AD) and multiple sclerosis (MS): DA1 (disease-associated1, associated with immunogenic genes), DA2 (disease-associated2, associated with genes influencing survival), and IFN (associated with interferon response genes). Spatial analysis of disease-associated oligodendrocytes (DAOs) in the cuprizone model reveals that DA1 and DA2 are established outside of the lesion area during demyelination and that DA1 repopulates the lesion during remyelination. Independent meta-analysis of human single-nucleus RNA-seq datasets reveals that the transcriptional responses of MS oligodendrocytes share features with mouse models. In contrast, the oligodendrocyte activation signature observed in human AD is largely distinct from those observed in mice. This catalog of oligodendrocyte activation states (http://research-pub.gene.com/OligoLandscape/) will be important to understand disease progression and develop therapeutic interventions.

Inferring the functional effects of mutation through clusters of mutations in homologous proteins.

Inferring functional consequences is a bottleneck in high-throughput cancer mutation discovery and genetic association studies. Most polymorphisms and germline mutations are unlikely to have functionally significant consequences. Most cancer somatic mutations do not contribute to tumorigenesis and are not under selective pressure. Identifying and understanding functionally important mutations can clarify disease biology and lead to new therapeutic and diagnostic opportunities. We investigated the extent to which protein mutations with functional consequences are enriched in clusters at conserved positions across related proteins. We found that disease-causing mutations form clusters more than random mutations or single nucleotide polymorphisms, confirming that mutation hotspots occur at the domain level. In addition to helping to identify functionally significant mutations, analysis of clustered mutations can indicate the mechanism and consequences for protein function. Our analysis focused on somatic cancer mutations suggests functional impact for many, including singleton mutations in FGFR1, FGFR3, GFI1B, PIK3CG, RALB, RAP2B, and STK11. This provides evidence and generates mechanistic hypotheses for the contribution of such mutations to cancer. The same approach can be applied to mutations suspected of involvement in other diseases. An interactive Web application for browsing mutation clusters is available at http://www.mcluster.org.