STING orchestrates the neuronal inflammatory stress response in multiple sclerosis.

Inflammation-induced neurodegeneration is a defining feature of multiple sclerosis (MS), yet the underlying mechanisms remain unclear. By dissecting the neuronal inflammatory stress response, we discovered that neurons in MS and its mouse model induce the stimulator of interferon genes (STING). However, activation of neuronal STING requires its detachment from the stromal interaction molecule 1 (STIM1), a process triggered by glutamate excitotoxicity. This detachment initiates non-canonical STING signaling, which leads to autophagic degradation of glutathione peroxidase 4 (GPX4), essential for neuronal redox homeostasis and thereby inducing ferroptosis. Both genetic and pharmacological interventions that target STING in neurons protect against inflammation-induced neurodegeneration. Our findings position STING as a central regulator of the detrimental neuronal inflammatory stress response, integrating inflammation with glutamate signaling to cause neuronal cell death, and present it as a tractable target for treating neurodegeneration in MS.

Small-molecule-induced epigenetic rejuvenation promotes SREBP condensation and overcomes barriers to CNS myelin regeneration.

Remyelination failure in diseases like multiple sclerosis (MS) was thought to involve suppressed maturation of oligodendrocyte precursors; however, oligodendrocytes are present in MS lesions yet lack myelin production. We found that oligodendrocytes in the lesions are epigenetically silenced. Developing a transgenic reporter labeling differentiated oligodendrocytes for phenotypic screening, we identified a small-molecule epigenetic-silencing-inhibitor (ESI1) that enhances myelin production and ensheathment. ESI1 promotes remyelination in animal models of demyelination and enables de novo myelinogenesis on regenerated CNS axons. ESI1 treatment lengthened myelin sheaths in human iPSC-derived organoids and augmented (re)myelination in aged mice while reversing age-related cognitive decline. Multi-omics revealed that ESI1 induces an active chromatin landscape that activates myelinogenic pathways and reprograms metabolism. Notably, ESI1 triggered nuclear condensate formation of master lipid-metabolic regulators SREBP1/2, concentrating transcriptional co-activators to drive lipid/cholesterol biosynthesis. Our study highlights the potential of targeting epigenetic silencing to enable CNS myelin regeneration in demyelinating diseases and aging.

Cellular architecture of evolving neuroinflammatory lesions and multiple sclerosis pathology.

Multiple sclerosis (MS) is a neurological disease characterized by multifocal lesions and smoldering pathology. Although single-cell analyses provided insights into cytopathology, evolving cellular processes underlying MS remain poorly understood. We investigated the cellular dynamics of MS by modeling temporal and regional rates of disease progression in mouse experimental autoimmune encephalomyelitis (EAE). By performing single-cell spatial expression profiling using in situ sequencing (ISS), we annotated disease neighborhoods and found centrifugal evolution of active lesions. We demonstrated that disease-associated (DA)-glia arise independently of lesions and are dynamically induced and resolved over the disease course. Single-cell spatial mapping of human archival MS spinal cords confirmed the differential distribution of homeostatic and DA-glia, enabled deconvolution of active and inactive lesions into sub-compartments, and identified new lesion areas. By establishing a spatial resource of mouse and human MS neuropathology at a single-cell resolution, our study unveils the intricate cellular dynamics underlying MS.

An inducible genetic tool to track and manipulate specific microglial states reveals their plasticity and roles in remyelination.

Recent single-cell RNA sequencing studies have revealed distinct microglial states in development and disease. These include proliferative-region-associated microglia (PAMs) in developing white matter and disease-associated microglia (DAMs) prevalent in various neurodegenerative conditions. PAMs and DAMs share a similar core gene signature. However, the extent of the dynamism and plasticity of these microglial states, as well as their functional significance, remains elusive, partly due to the lack of specific tools. Here, we generated an inducible Cre driver line, Clec7a-CreERT2, that targets PAMs and DAMs in the brain parenchyma. Utilizing this tool, we profiled labeled cells during development and in several disease models, uncovering convergence and context-dependent differences in PAM and DAM gene expression. Through long-term tracking, we demonstrated microglial state plasticity. Lastly, we specifically depleted DAMs in demyelination, revealing their roles in disease recovery. Together, we provide a versatile genetic tool to characterize microglial states in CNS development and disease.

Modulation of multiple sclerosis risk and pathogenesis by the gut microbiota: Complex interactions between host genetics, bacterial metabolism, and diet.

Multiple sclerosis (MS) is an autoimmune disease of the central nervous system, affecting nearly 2 million people worldwide. The etiology of MS is multifactorial: Approximately 30% of the MS risk is genetic, which implies that the remaining ~70% is environmental, with a number of factors proposed. One recently implicated risk factor for MS is the composition of the gut microbiome. Numerous case-control studies have identified changes in gut microbiota composition of people with MS (pwMS) compared with healthy control individuals, and more recent studies in animal models have begun to identify the causative microbes and underlying mechanisms. Here, we review some of these mechanisms, with a specific focus on the role of host genetic variation, dietary inputs, and gut microbial metabolism, with a particular emphasis on short-chain fatty acid and tryptophan metabolism. We put forward a model where, in an individual genetically susceptible to MS, the gut microbiota and diet can synergize as potent environmental modifiers of disease risk and possibly progression, with diet-dependent gut microbial metabolites serving as a key mechanism. We also propose that specific microbial taxa may have divergent effects in individuals carrying distinct variants of MS risk alleles or other polymorphisms, as a consequence of host gene-by-gut microbiota interactions. Finally, we also propose that the effects of specific microbial taxa, especially those that exert their effects through metabolites, are highly dependent on the host dietary intake. What emerges is a complex multifaceted interaction that has been challenging to disentangle in human studies, contributing to the divergence of findings across heterogeneous cohorts with differing geography, dietary preferences, and genetics. Nonetheless, this provides a complex and individualized, yet tractable, model of how the gut microbiota regulate susceptibility to MS, and potentially progression of this disease. Thus, we conclude that prophylactic or therapeutic modulation of the gut microbiome to prevent or treat MS will require a careful and personalized consideration of host genetics, baseline gut microbiota composition, and dietary inputs.

The influence of aging and the microbiome in multiple sclerosis and other neurologic diseases.

The human gut microbiome is well-recognized as a key player in maintaining health. However, it is a dynamic entity that changes across the lifespan. How the microbial changes that occur in later decades of life shape host health or impact age-associated inflammatory neurological diseases such as multiple sclerosis (MS) is still unclear. Current understanding of the aging gut microbiome is largely limited to cross-sectional observational studies. Moreover, studies in humans are limited by confounding host-intrinsic and extrinsic factors that are not easily disentangled from aging. This review provides a comprehensive summary of existing literature on the aging gut microbiome and its known relationships with neurological diseases, with a specific focus on MS. We will also discuss preclinical animal models and human studies that shed light on the complex microbiota-host interactions that have the potential to influence disease pathology and progression in aging individuals. Lastly, we propose potential avenues of investigation to deconvolute features of an aging microbiota that contribute to disease, or alternatively promote health in advanced age.

Iron scavenging and myeloid cell polarization.

Myeloid cells that populate all human organs and blood are a versatile class of innate immune cells. They are crucial for sensing and regulating processes as diverse as tissue homeostasis and inflammation and are frequently characterized by their roles in either regulating or promoting inflammation. Recent studies in cultured cells and mouse models highlight the role of iron in skewing the functional properties of myeloid cells in tissue damage and repair. Here, we review certain emerging concepts on how iron influences and determines myeloid cell polarization in the context of its uptake, storage, and metabolism, including in conditions such as multiple sclerosis (MS), sickle cell disease, and tumors.

Inflammasome diversity: exploring novel frontiers in the innate immune response.

Pathogens elicit complex mammalian immune responses by activating multiple sensors within inflammasomes, which recognize diverse pathogen-associated molecular patterns (PAMPs) and damage-associated molecular patterns (DAMPs). This simultaneous activation induces the formation of protein complexes referred to as multiple inflammasomes, that orchestrate a spectrum of programmed cell death pathways, including pyroptosis, apoptosis, and necroptosis. This concept is crucial for comprehending the complexity of the innate immune system's response to diverse pathogens and its implications for various diseases. Novel contributions here include emphasizing simultaneous sensor activation by pathogens, proposing the existence of multiple inflammasome complexes, and advocating for further exploration of their structural basis. Understanding these mechanisms may offer insights into disease pathogenesis, paving the way for potential therapeutic interventions targeting inflammasome-mediated immune responses.

Evidence supports a causal association between allele-specific vitamin D receptor binding and multiple sclerosis among Europeans.

Although evidence exists for a causal association between 25-hydroxyvitamin D (25(OH)D) serum levels, and multiple sclerosis (MS), the role of variation in vitamin D receptor (VDR) binding in MS is unknown. Here, we leveraged previously identified variants associated with allele imbalance in VDR binding (VDR-binding variant; VDR-BV) in ChIP-exo data from calcitriol-stimulated lymphoblastoid cell lines and 25(OH)D serum levels from genome-wide association studies to construct genetic instrumental variables (GIVs). GIVs are composed of one or more genetic variants that serve as proxies for exposures of interest. Here, GIVs for both VDR-BVs and 25(OH)D were used in a two-sample Mendelian Randomization study to investigate the relationship between VDR binding at a locus, 25(OH)D serum levels, and MS risk. Data for 13,598 MS cases and 38,887 controls of European ancestry from Kaiser Permanente Northern California, Swedish MS studies, and the UK Biobank were included. We estimated the association between each VDR-BV GIV and MS. Significant interaction between a VDR-BV GIV and a GIV for serum 25OH(D) was evidence for a causal association between VDR-BVs and MS unbiased by pleiotropy. We observed evidence for associations between two VDR-BVs (rs2881514, rs2531804) and MS after correction for multiple tests. There was evidence of interaction between rs2881514 and a 25(OH)D GIV, providing evidence of a causal association between rs2881514 and MS. This study is the first to demonstrate evidence that variation in VDR binding at a locus contributes to MS risk. Our results are relevant to other autoimmune diseases in which vitamin D plays a role.

Model selection over partially ordered sets.

In problems such as variable selection and graph estimation, models are characterized by Boolean logical structure such as the presence or absence of a variable or an edge. Consequently, false-positive error or false-negative error can be specified as the number of variables/edges that are incorrectly included or excluded in an estimated model. However, there are several other problems such as ranking, clustering, and causal inference in which the associated model classes do not admit transparent notions of false-positive and false-negative errors due to the lack of an underlying Boolean logical structure. In this paper, we present a generic approach to endow a collection of models with partial order structure, which leads to a hierarchical organization of model classes as well as natural analogs of false-positive and false-negative errors. We describe model selection procedures that provide false-positive error control in our general setting, and we illustrate their utility with numerical experiments.

Targeting the muscarinic M1 receptor with a selective, brain-penetrant antagonist to promote remyelination in multiple sclerosis.

Multiple sclerosis (MS) is a chronic and debilitating neurological disease that results in inflammatory demyelination. While endogenous remyelination helps to recover function, this restorative process tends to become less efficient over time. Currently, intense efforts aimed at the mechanisms that promote remyelination are being considered promising therapeutic approaches. The M1 muscarinic acetylcholine receptor (M1R) was previously identified as a negative regulator of oligodendrocyte differentiation and myelination. Here, we validate M1R as a target for remyelination by characterizing expression in human and rodent oligodendroglial cells (including those in human MS tissue) using a highly selective M1R probe. As a breakthrough to conventional methodology, we conjugated a fluorophore to a highly M1R selective peptide (MT7) which targets the M1R in the subnanomolar range. This allows for exceptional detection of M1R protein expression in the human CNS. More importantly, we introduce PIPE-307, a brain-penetrant, small-molecule antagonist with favorable drug-like properties that selectively targets M1R. We evaluate PIPE-307 in a series of in vitro and in vivo studies to characterize potency and selectivity for M1R over M2-5R and confirm the sufficiency of blocking this receptor to promote differentiation and remyelination. Further, PIPE-307 displays significant efficacy in the mouse experimental autoimmune encephalomyelitis model of MS as evaluated by quantifying disability, histology, electron microscopy, and visual evoked potentials. Together, these findings support targeting M1R for remyelination and support further development of PIPE-307 for clinical studies.

Electrical stimulation of the vagus nerve ameliorates inflammation and disease activity in a rat EAE model of multiple sclerosis.

Multiple sclerosis (MS) is a demyelinating central nervous system (CNS) disorder that is associated with functional impairment and accruing disability. There are multiple U.S. Food and Drug Administration (FDA)-approved drugs that effectively dampen inflammation and slow disability progression. However, these agents do not work well for all patients and are associated with side effects that may limit their use. The vagus nerve (VN) provides a direct communication conduit between the CNS and the periphery, and modulation of the inflammatory reflex via electrical stimulation of the VN (VNS) shows efficacy in ameliorating pathology in several CNS and autoimmune disorders. We therefore investigated the impact of VNS in a rat experimental autoimmune encephalomyelitis (EAE) model of MS. In this study, VNS-mediated neuroimmune modulation is demonstrated to effectively decrease EAE disease severity and duration, infiltration of neutrophils and pathogenic lymphocytes, myelin damage, blood-brain barrier disruption, fibrinogen deposition, and proinflammatory microglial activation. VNS modulates expression of genes that are implicated in MS pathogenesis, as well as those encoding myelin proteins and transcription factors regulating new myelin synthesis. Together, these data indicate that neuroimmune modulation via VNS may be a promising approach to treat MS, that not only ameliorates symptoms but potentially also promotes myelin repair (remyelination).

Astrocytic TIMP-1 regulates production of Anastellin, an inhibitor of oligodendrocyte differentiation and FTY720 responses.

Astrocyte activation is associated with neuropathology and the production of tissue inhibitor of metalloproteinase-1 (TIMP1). TIMP1 is a pleiotropic extracellular protein that functions both as a protease inhibitor and as a growth factor. Astrocytes that lack expression of Timp1 do not support rat oligodendrocyte progenitor cell (rOPC) differentiation, and adult global Timp1 knockout (Timp1KO) mice do not efficiently remyelinate following a demyelinating injury. Here, we performed an unbiased proteomic analysis and identified a fibronectin-derived peptide called Anastellin (Ana) that was unique to the Timp1KO astrocyte secretome. Ana was found to block rOPC differentiation in vitro and enhanced the inhibitory influence of fibronectin on rOPC differentiation. Ana is known to act upon the sphingosine-1-phosphate receptor 1, and we determined that Ana also blocked the pro-myelinating effect of FTY720 (or fingolimod) on rOPC differentiation in vitro. Administration of FTY720 to wild-type C57BL/6 mice during MOG35-55-experimental autoimmune encephalomyelitis ameliorated clinical disability while FTY720 administered to mice lacking expression of Timp1 (Timp1KO) had no effect. Analysis of Timp1 and fibronectin (FN1) transcripts from primary human astrocytes from healthy and multiple sclerosis (MS) donors revealed lower TIMP1 expression was coincident with elevated FN1 in MS astrocytes. Last, analyses of proteomic databases of MS samples identified Ana peptides to be more abundant in the cerebrospinal fluid (CSF) of human MS patients with high disease activity. A role for Ana in MS as a consequence of a lack of astrocytic TIMP-1 production could influence both the efficacy of fingolimod responses and innate remyelination potential in the MS brain.

Predicting protein conformational motions using energetic frustration analysis and AlphaFold2.

Proteins perform their biological functions through motion. Although high throughput prediction of the three-dimensional static structures of proteins has proved feasible using deep-learning-based methods, predicting the conformational motions remains a challenge. Purely data-driven machine learning methods encounter difficulty for addressing such motions because available laboratory data on conformational motions are still limited. In this work, we develop a method for generating protein allosteric motions by integrating physical energy landscape information into deep-learning-based methods. We show that local energetic frustration, which represents a quantification of the local features of the energy landscape governing protein allosteric dynamics, can be utilized to empower AlphaFold2 (AF2) to predict protein conformational motions. Starting from ground state static structures, this integrative method generates alternative structures as well as pathways of protein conformational motions, using a progressive enhancement of the energetic frustration features in the input multiple sequence alignment sequences. For a model protein adenylate kinase, we show that the generated conformational motions are consistent with available experimental and molecular dynamics simulation data. Applying the method to another two proteins KaiB and ribose-binding protein, which involve large-amplitude conformational changes, can also successfully generate the alternative conformations. We also show how to extract overall features of the AF2 energy landscape topography, which has been considered by many to be black box. Incorporating physical knowledge into deep-learning-based structure prediction algorithms provides a useful strategy to address the challenges of dynamic structure prediction of allosteric proteins.

Functional effects of mutations in proteins can be predicted and interpreted by guided selection of sequence covariation information.

Predicting the effects of one or more mutations to the in vivo or in vitro properties of a wild-type protein is a major computational challenge, due to the presence of epistasis, that is, of interactions between amino acids in the sequence. We introduce a computationally efficient procedure to build minimal epistatic models to predict mutational effects by combining evolutionary (homologous sequence) and few mutational-scan data. Mutagenesis measurements guide the selection of links in a sparse graphical model, while the parameters on the nodes and the edges are inferred from sequence data. We show, on 10 mutational scans, that our pipeline exhibits performances comparable to state-of-the-art deep networks trained on many more data, while requiring much less parameters and being hence more interpretable. In particular, the identified interactions adapt to the wild-type protein and to the fitness or biochemical property experimentally measured, mostly focus on key functional sites, and are not necessarily related to structural contacts. Therefore, our method is able to extract information relevant for one mutational experiment from homologous sequence data reflecting the multitude of structural and functional constraints acting on proteins throughout evolution.

Expanded T lymphocytes in the cerebrospinal fluid of multiple sclerosis patients are specific for Epstein-Barr-virus-infected B cells.

Epstein-Barr virus (EBV) infection has long been associated with multiple sclerosis (MS), but the role of EBV in the pathogenesis of MS is not clear. Our hypothesis is that a major fraction of the expanded clones of T lymphocytes in the cerebrospinal fluid (CSF) are specific for autologous EBV-infected B cells. We obtained blood and CSF samples from eight relapsing-remitting patients in the process of diagnosis. We stimulated cells from the blood with autologous EBV-infected lymphoblastoid cell lines (LCL), EBV, varicella zoster virus, influenza, and candida and sorted the responding cells with flow cytometry after 6 d. We sequenced the RNA for T cell receptors (TCR) from CSF, unselected blood cells, and the antigen-specific cells. We used the TCR Vß CDR3 sequences from the antigen-specific cells to assign antigen specificity to the sequences from the CSF and blood. LCL-specific cells comprised 13.0 ± 4.3% (mean ± SD) of the total reads present in CSF and 13.3 ± 7.5% of the reads present in blood. The next most abundant antigen specificity was flu, which was 4.7 ± 1.7% of the reads in the CSF and 9.3 ± 6.6% in the blood. The prominence of LCL-specific reads was even more marked in the top 1% most abundant CSF clones with statistically significant 47% mean overlap with LCL. We conclude that LCL-specific sequences form a major portion of the TCR repertoire in both CSF and blood and that expanded clones specific for LCL are present in MS CSF. This has important implications for the pathogenesis of MS.

Multiple retinal isomerizations during the early phase of the bestrhodopsin photoreaction.

Bestrhodopsins constitute a class of light-regulated pentameric ion channels that consist of one or two rhodopsins in tandem fused with bestrophin ion channel domains. Here, we report on the isomerization dynamics in the rhodopsin tandem domains of Phaeocystis antarctica bestrhodopsin, which binds all-trans retinal Schiff-base (RSB) absorbing at 661 nm and, upon illumination, converts to the meta-stable P540 state with an unusual 11-cis RSB. The primary photoproduct P682 corresponds to a mixture of highly distorted 11-cis and 13-cis RSB directly formed from the excited state in 1.4 ps. P673 evolves from P682 in 500 ps and contains highly distorted 13-cis RSB, indicating that the 11-cis fraction in P682 converts to 13-cis. Next, P673 establishes an equilibrium with P595 in 1.2 µs, during which RSB converts to 11-cis and then further proceeds to P560 in 48 µs and P540 in 1.0 ms while remaining 11-cis. Hence, extensive isomeric switching occurs on the early ground state potential energy surface (PES) on the hundreds of ps to µs timescale before finally settling on a metastable 11-cis photoproduct. We propose that P682 and P673 are trapped high up on the ground-state PES after passing through either of two closely located conical intersections that result in 11-cis and 13-cis RSB. Co-rotation of C11=C12 and C13=C14 bonds results in a constricted conformational landscape that allows thermal switching between 11-cis and 13-cis species of highly strained RSB chromophores. Protein relaxation may release RSB strain, allowing it to evolve to a stable 11-cis isomeric configuration in microseconds.

Multiscale embedded printing of engineered human tissue and organ equivalents.

Creating tissue and organ equivalents with intricate architectures and multiscale functional feature sizes is the first step toward the reconstruction of transplantable human tissues and organs. Existing embedded ink writing approaches are limited by achievable feature sizes ranging from hundreds of microns to tens of millimeters, which hinders their ability to accurately duplicate structures found in various human tissues and organs. In this study, a multiscale embedded printing (MSEP) strategy is developed, in which a stimuli-responsive yield-stress fluid is applied to facilitate the printing process. A dynamic layer height control method is developed to print the cornea with a smooth surface on the order of microns, which can effectively overcome the layered morphology in conventional extrusion-based three-dimensional bioprinting methods. Since the support bath is sensitive to temperature change, it can be easily removed after printing by tuning the ambient temperature, which facilitates the fabrication of human eyeballs with optic nerves and aortic heart valves with overhanging leaflets on the order of a few millimeters. The thermosensitivity of the support bath also enables the reconstruction of the full-scale human heart on the order of tens of centimeters by on-demand adding support bath materials during printing. The proposed MSEP demonstrates broader printable functional feature sizes ranging from microns to centimeters, providing a viable and reliable technical solution for tissue and organ printing in the future.

betAS: intuitive analysis and visualization of differential alternative splicing using beta distributions.

Next-generation RNA sequencing allows alternative splicing (AS) quantification with unprecedented resolution, with the relative inclusion of an alternative sequence in transcripts being commonly quantified by the proportion of reads supporting it as percent spliced-in (PSI). However, PSI values do not incorporate information about precision, proportional to the respective AS events' read coverage. Beta distributions are suitable to quantify inclusion levels of alternative sequences, using reads supporting their inclusion and exclusion as surrogates for the two distribution shape parameters. Each such beta distribution has the PSI as its mean value and is narrower when the read coverage is higher, facilitating the interpretability of its precision when plotted. We herein introduce a computational pipeline, based on beta distributions accurately modeling PSI values and their precision, to quantitatively and visually compare AS between groups of samples. Our methodology includes a differential splicing significance metric that compromises the magnitude of intergroup differences, the estimation uncertainty in individual samples, and the intragroup variability, being therefore suitable for multiple-group comparisons. To make our approach accessible and clear to both noncomputational and computational biologists, we developed betAS, an interactive web app and user-friendly R package for visual and intuitive differential splicing analysis from read count data.

Inefficient recruitment of DDX39B impedes pre-spliceosome assembly on FOXP3 introns.

Forkhead box P3 (FOXP3) is the master fate-determining transcription factor in regulatory T (Treg) cells and is essential for their development, function, and homeostasis. Mutations in FOXP3 cause immunodysregulation polyendocrinopathy enteropathy X-linked (IPEX) syndrome, and aberrant expression of FOXP3 has been implicated in other diseases such as multiple sclerosis and cancer. We previously demonstrated that pre-mRNA splicing of FOXP3 RNAs is highly sensitive to levels of DExD-box polypeptide 39B (DDX39B), and here we investigate the mechanism of this sensitivity. FOXP3 introns have cytidine (C)-rich/uridine (U)-poor polypyrimidine (py) tracts that are responsible for their inefficient splicing and confer sensitivity to DDX39B. We show that there is a deficiency in the assembly of commitment complexes (CCs) on FOXP3 introns, which is consistent with the lower affinity of U2AF2 for C-rich/U-poor py tracts. Our data indicate an even stronger effect on the conversion of CCs to pre-spliceosomes. We propose that this is due to an altered conformation that U2AF2 adopts when it binds to C-rich/U-poor py tracts and that this conformation has a lower affinity for DDX39B. As a consequence, CCs assembled on FOXP3 introns are defective in recruiting DDX39B, and this leads to the inefficient assembly of pre-spliceosome complexes.