SOX10 regulates expression of the SH3-domain kinase binding protein 1 (Sh3kbp1) locus in Schwann cells via an alternative promoter

Single-cell sequencing reveals glial cell involvement in development of neuropathic pain via myelin sheath lesion formation in the spinal cord

BACKGROUND

Neuropathic pain (NP), which results from injury or lesion of the somatosensory nervous system, is intimately associated with glial cells. The roles of microglia and astrocytes in NP have been broadly described, while studies on oligodendrocytes have largely focused on axonal myelination. The mechanisms of oligodendrocytes and their interactions with other glial cells in NP development remain uncertain.

METHODS

To explore the function of the interaction of the three glial cells and their interactions on myelin development in NP, we evaluated changes in NP and myelin morphology after a chronic constriction injury (CCI) model in mice, and used single-cell sequencing to reveal the subpopulations characteristics of oligodendrocytes, microglia, and astrocytes in the spinal cord tissues, as well as their relationship with myelin lesions; the proliferation and differentiation trajectories of oligodendrocyte subpopulations were also revealed using pseudotime cell trajectory and RNA velocity analysis. In addition, we identified chemokine ligand-receptor pairs between glial cells by cellular communication and verified them using immunofluorescence.

RESULTS

Our study showed that NP peaked on day 7 after CCI in mice, a time at which myelin lesions were present in both the spinal cord and sciatic nerve. Oligodendrocytes, microglia, and astrocytes subpopulations in spinal cord tissue were heterogeneous after CCI and all were involved in suppressing the process of immune defense and myelin production. In addition, the differentiation trajectory of oligodendrocytes involved a unidirectional lattice process of OPC-1-Oligo-9, which was arrested at the Oligo-2 stage under the influence of microglia and astrocytes. And the CADM1-CADM1, NRP1-VEGFA interactions between glial cells are enhanced after CCI and they had a key role in myelin lesions and demyelination.

CONCLUSIONS

Our study reveals the close relationship between the differentiation block of oligodendrocytes after CCI and their interaction with microglia and astrocytes-mediated myelin lesions and NP. CADM1/CADM1 and NRP-1/VEGFA may serve as potential therapeutic targets for use in the treatment of NP.

SH3KBP1 Promotes Glioblastoma Tumorigenesis by Activating EGFR Signaling

Glioblastoma (GBM) is the most common and aggressive brain tumor in adults. Overexpression or activation of epidermal growth factor receptor (EGFR) occurs commonly in multiple human cancers and promotes tumorigenesis. However, the underlying molecular mechanism of EGFR aberrant activation and the downstream signaling pathways remains largely unknown. In this study, we report that both SH3-domain kinase binding protein 1 (SH3KBP1) mRNA and protein levels are highly expressed in GBM and its high expression is associated with worse survival of glioma patients. In addition, we provide evidence that SH3KBP1 is prominently expressed in GBM stem cells (GSCs) and have potential to serve as a novel GSCs marker. Moreover, silencing SH3KBP1 dramatically impairs GBM cell proliferation, migration and GSCs self-renewal ability in vitro and xenograft tumors growth in vivo. Most importantly, we found that SH3KBP1 directly interacts with EGFR and may act as an adaptor protein to transduce EGFR signaling. Together, our work uncovers SH3KBP1 as a novel regulator of oncogenic EGFR signaling and also as a potential therapeutic target for GBM patients with EGFR activation.

SOX10-regulated promoter use defines isoform-specific gene expression in Schwann cells

Background

Multicellular organisms adopt various strategies to tailor gene expression to cellular contexts including the employment of multiple promoters (and the associated transcription start sites (TSSs)) at a single locus that encodes distinct gene isoforms. Schwann cells—the myelinating cells of the peripheral nervous system (PNS)—exhibit a specialized gene expression profile directed by the transcription factor SOX10, which is essential for PNS myelination. SOX10 regulates promoter elements associated with unique TSSs and gene isoforms at several target loci, implicating SOX10-mediated, isoform-specific gene expression in Schwann cell function. Here, we report on genome-wide efforts to identify SOX10-regulated promoters and TSSs in Schwann cells to prioritize genes and isoforms for further study.

Results

We performed global TSS analyses and mined previously reported ChIP-seq datasets to assess the activity of SOX10-bound promoters in three models: (i) an adult mammalian nerve; (ii) differentiating primary Schwann cells, and (iii) cultured Schwann cells with ablated SOX10 function. We explored specific characteristics of SOX10-dependent TSSs, which provides confidence in defining them as SOX10 targets. Finally, we performed functional studies to validate our findings at four previously unreported SOX10 target loci: ARPC1A, CHN2, DDR1, and GAS7. These findings suggest roles for the associated SOX10-regulated gene products in PNS myelination.

Conclusions

In sum, we provide comprehensive computational and functional assessments of SOX10-regulated TSS use in Schwann cells. The data presented in this study will stimulate functional studies on the specific mRNA and protein isoforms that SOX10 regulates, which will improve our understanding of myelination in the peripheral nerve.

A genome-wide assessment of conserved SNP alleles reveals a panel of regulatory SNPs relevant to the peripheral nerve

Background

Identifying functional non-coding variation is critical for defining the genetic contributions to human disease. While single-nucleotide polymorphisms (SNPs) within cis-acting transcriptional regulatory elements have been implicated in disease pathogenesis, not all cell types have been assessed and functional validations have been limited. In particular, the cells of the peripheral nervous system have been excluded from genome-wide efforts to link non-coding SNPs to altered gene function. Addressing this gap is essential for defining the genetic architecture of diseases that affect the peripheral nerve. We developed a computational pipeline to identify SNPs that affect regulatory function (rSNPs) and evaluated our predictions on a set of 144 regions in Schwann cells, motor neurons, and muscle cells.

Results

We identified 28 regions that display regulatory activity in at least one cell type and 13 SNPs that affect regulatory function. We then tailored our pipeline to one peripheral nerve cell type by incorporating SOX10 ChIP-Seq data; SOX10 is essential for Schwann cells. We prioritized 22 putative SOX10 response elements harboring a SNP and rapidly validated two rSNPs. We then selected one of these elements for further characterization to assess the biological relevance of our approach. Deletion of the element from the genome of cultured Schwann cells—followed by differential gene expression studies—revealed Tubb2b as a candidate target gene. Studying the enhancer in developing mouse embryos revealed activity in SOX10-positive cells including the dorsal root ganglia and melanoblasts.

Conclusions

Our efforts provide insight into the utility of employing strict conservation for rSNP discovery. This strategy, combined with functional analyses, can yield candidate target genes. In support of this, our efforts suggest that investigating the role of Tubb2b in SOX10-positive cells may reveal novel biology within these cell populations.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4692-z) contains supplementary material, which is available to authorized users.

Stringent comparative sequence analysis reveals SOX10 as a putative inhibitor of glial cell differentiation

Background

The transcription factor SOX10 is essential for all stages of Schwann cell development including myelination. SOX10 cooperates with other transcription factors to activate the expression of key myelin genes in Schwann cells and is therefore a context-dependent, pro-myelination transcription factor. As such, the identification of genes regulated by SOX10 will provide insight into Schwann cell biology and related diseases. While genome-wide studies have successfully revealed SOX10 target genes, these efforts mainly focused on myelinating stages of Schwann cell development. We propose that less-biased approaches will reveal novel functions of SOX10 outside of myelination.

Results

We developed a stringent, computational-based screen for genome-wide identification of SOX10 response elements. Experimental validation of a pilot set of predicted binding sites in multiple systems revealed that SOX10 directly regulates a previously unreported alternative promoter at SOX6, which encodes a transcription factor that inhibits glial cell differentiation. We further explored the utility of our computational approach by combining it with DNase-seq analysis in cultured Schwann cells and previously published SOX10 ChIP-seq data from rat sciatic nerve. Remarkably, this analysis enriched for genomic segments that map to loci involved in the negative regulation of gliogenesis including SOX5, SOX6, NOTCH1, HMGA2, HES1, MYCN, ID4, and ID2. Functional studies in Schwann cells revealed that: (1) all eight loci are expressed prior to myelination and down-regulated subsequent to myelination; (2) seven of the eight loci harbor validated SOX10 binding sites; and (3) seven of the eight loci are down-regulated upon repressing SOX10 function.

Conclusions

Our computational strategy revealed a putative novel function for SOX10 in Schwann cells, which suggests a model where SOX10 activates the expression of genes that inhibit myelination during non-myelinating stages of Schwann cell development. Importantly, the computational and functional datasets we present here will be valuable for the study of transcriptional regulation, SOX protein function, and glial cell biology.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3167-3) contains supplementary material, which is available to authorized users.

SOX10 regulates an alternative promoter at the Charcot-Marie-Tooth disease locus MTMR2

Schwann cells are the myelinating glia of the peripheral nervous system and dysfunction of these cells causes motor and sensory peripheral neuropathy. The transcription factor SOX10 is critical for Schwann cell development and maintenance, and many SOX10 target genes encode proteins required for Schwann cell function. Loss-of-function mutations in the gene encoding myotubularin-related protein 2 (MTMR2) cause Charcot-Marie-Tooth disease type 4B1 (CMT4B1), a severe demyelinating peripheral neuropathy characterized by myelin outfoldings along peripheral nerves. Previous reports indicate that MTMR2 is ubiquitously expressed making it unclear how loss of this gene causes a Schwann cell-specific phenotype. To address this, we performed computational and functional analyses at MTMR2 to identify transcriptional regulatory elements important for Schwann cell expression. Through these efforts, we identified an alternative, SOX10-responsive promoter at MTMR2 that displays strong regulatory activity in immortalized rat Schwann (S16) cells. This promoter directs transcription of a previously unidentified MTMR2 transcript that is enriched in mouse Schwann cells compared to immortalized mouse motor neurons (MN-1), and is predicted to encode an N-terminally truncated protein isoform. The expression of the endogenous transcript is induced in a heterologous cell line by ectopically expressing SOX10, and is nearly ablated in Schwann cells by impairing SOX10 function. Intriguingly, overexpressing the two MTMR2 protein isoforms in HeLa cells revealed that both localize to nuclear puncta and the shorter isoform displays higher nuclear localization compared to the longer isoform. Combined, our data warrant further investigation of the truncated MTMR2 protein isoform in Schwann cells and in CMT4B1 pathogenesis.

Haplotype-specific modulation of a SOX10/CREB response element at the Charcot-Marie-Tooth disease type 4C locus SH3TC2

Loss-of-function mutations in the Src homology 3 (SH3) domain and tetratricopeptide repeats 2 (SH3TC2) gene cause autosomal recessive demyelinating Charcot-Marie-Tooth neuropathy. The SH3TC2 protein has been implicated in promyelination signaling through axonal neuregulin-1 and the ERBB2 Schwann cell receptor. However, little is known about the transcriptional regulation of the SH3TC2 gene. We performed computational and functional analyses that revealed two cis-acting regulatory elements at SH3TC2-one at the promoter and one ∼150 kb downstream of the transcription start site. Both elements direct reporter gene expression in Schwann cells and are responsive to the transcription factor SOX10, which is essential for peripheral nervous system myelination. The downstream enhancer harbors a single-nucleotide polymorphism (SNP) that causes an ∼80% reduction in enhancer activity. The SNP resides directly within a predicted binding site for the transcription factor cAMP response element binding protein (CREB), and we demonstrate that this regulatory element binds to CREB and is activated by CREB expression. Finally, forskolin induces Sh3tc2 expression in rat primary Schwann cells, indicating that SH3TC2 is a CREB target gene. These findings prompted us to determine if SNP genotypes at SH3TC2 are associated with differential phenotypes in the most common demyelinating peripheral neuropathy, CMT1A. Interestingly, this revealed several associations between SNP alleles and disease severity. In summary, our data indicate that SH3TC2 is regulated by the transcription factors CREB and SOX10, define a regulatory SNP at this disease-associated locus and reveal SH3TC2 as a candidate modifier locus of CMT disease phenotypes.

Identification of HMX1 target genes: a predictive promoter model approach

PURPOSE

A homozygous mutation in the H6 family homeobox 1 (HMX1) gene is responsible for a new oculoauricular defect leading to eye and auricular developmental abnormalities as well as early retinal degeneration (MIM 612109). However, the HMX1 pathway remains poorly understood, and in the first approach to better understand the pathway's function, we sought to identify the target genes.

METHODS

We developed a predictive promoter model (PPM) approach using a comparative transcriptomic analysis in the retina at P15 of a mouse model lacking functional Hmx1 (dmbo mouse) and its respective wild-type. This PPM was based on the hypothesis that HMX1 binding site (HMX1-BS) clusters should be more represented in promoters of HMX1 target genes. The most differentially expressed genes in the microarray experiment that contained HMX1-BS clusters were used to generate the PPM, which was then statistically validated. Finally, we developed two genome-wide target prediction methods: one that focused on conserving PPM features in human and mouse and one that was based on the co-occurrence of HMX1-BS pairs fitting the PPM, in human or in mouse, independently.

RESULTS

The PPM construction revealed that sarcoglycan, gamma (35kDa dystrophin-associated glycoprotein) (Sgcg), teashirt zinc finger homeobox 2 (Tshz2), and solute carrier family 6 (neurotransmitter transporter, glycine) (Slc6a9) genes represented Hmx1 targets in the mouse retina at P15. Moreover, the genome-wide target prediction revealed that mouse genes belonging to the retinal axon guidance pathway were targeted by Hmx1. Expression of these three genes was experimentally validated using a quantitative reverse transcription PCR approach. The inhibitory activity of Hmx1 on Sgcg, as well as protein tyrosine phosphatase, receptor type, O (Ptpro) and Sema3f, two targets identified by the PPM, were validated with luciferase assay.

CONCLUSIONS

Gene expression analysis between wild-type and dmbo mice allowed us to develop a PPM that identified the first target genes of Hmx1.

Proline: the distribution, frequency, positioning, and common functional roles of proline and polyproline sequences in the human proteome

Proline is an anomalous amino acid. Its nitrogen atom is covalently locked within a ring, thus it is the only proteinogenic amino acid with a constrained phi angle. Sequences of three consecutive prolines can fold into polyproline helices, structures that join alpha helices and beta pleats as architectural motifs in protein configuration. Triproline helices are participants in protein-protein signaling interactions. Longer spans of repeat prolines also occur, containing as many as 27 consecutive proline residues. Little is known about the frequency, positioning, and functional significance of these proline sequences. Therefore we have undertaken a systematic bioinformatics study of proline residues in proteins. We analyzed the distribution and frequency of 687,434 proline residues among 18,666 human proteins, identifying single residues, dimers, trimers, and longer repeats. Proline accounts for 6.3% of the 10,882,808 protein amino acids. Of all proline residues, 4.4% are in trimers or longer spans. We detected patterns that influence function based on proline location, spacing, and concentration. We propose a classification based on proline-rich, polyproline-rich, and proline-poor status. Whereas singlet proline residues are often found in proteins that display recurring architectural patterns, trimers or longer proline sequences tend be associated with the absence of repetitive structural motifs. Spans of 6 or more are associated with DNA/RNA processing, actin, and developmental processes. We also suggest a role for proline in Kruppel-type zinc finger protein control of DNA expression, and in the nucleation and translocation of actin by the formin complex.

Genome-wide analysis of EGR2/SOX10 binding in myelinating peripheral nerve

Myelin is essential for the rapidity of saltatory nerve conduction, and also provides trophic support for axons to prevent axonal degeneration. Two critical determinants of myelination are SOX10 and EGR2/KROX20. SOX10 is required for specification of Schwann cells from neural crest, and is required at every stage of Schwann cell development. Egr2/Krox20 expression is activated by axonal signals in myelinating Schwann cells, and is required for cell cycle arrest and myelin formation. To elucidate the integrated function of these two transcription factors during peripheral nerve myelination, we performed in vivo ChIP-Seq analysis of myelinating peripheral nerve. Integration of these binding data with loss-of-function array data identified a range of genes regulated by these factors. In addition, although SOX10 itself regulates Egr2/Krox20 expression, leading to coordinate activation of several major myelin genes by the two factors, there is a large subset of genes that are activated independent of EGR2. Finally, the results identify a set of SOX10-dependent genes that are expressed in early Schwann cell development, but become subsequently repressed by EGR2/KROX20.

Distal enhancers upstream of the Charcot-Marie-Tooth type 1A disease gene PMP22

Myelin insulates axons in the peripheral nervous system to allow rapid propagation of action potentials, and proper myelination requires the precise regulation of genes encoding myelin proteins, including PMP22. The correct gene dosage of PMP22 is critical; a duplication of PMP22 is the most common cause of the peripheral neuropathy Charcot-Marie-Tooth Disease (CMT) (classified as type 1A), while a deletion of PMP22 leads to another peripheral neuropathy, hereditary neuropathy with liability to pressure palsies. Recently, duplications upstream of PMP22, but not containing the gene itself, were reported in patients with CMT1A like symptoms, suggesting that this region contains regulators of PMP22. Using chromatin immunoprecipitation analysis of two transcription factors known to upregulate PMP22-EGR2 and SOX10-we found several enhancers in this upstream region that contain open chromatin and direct reporter gene expression in tissue culture and in vivo in zebrafish. These studies provide a novel means to identify critical regulatory elements in genes that are required for myelination, and elucidate the functional significance of non-coding genomic rearrangements.