Strange IndiaStrange India


Mice and EAE

The Cx3cr1-YFPcreERT2R26tdTomato fate mapping mouse was generated by crossing B6.129P2(Cg)-Cx3cr1tm2.1(cre/ERT2)Litt/WganJ52 with B6.Cg-Gt(ROSA)26Sortm9(CAG-tdTomato)Hze/J)53, as previously described19,20. ND6-P25L mice mice (referred to here as Nd6 mice) were generated as described previously40 and were obtained from D. Wallace. The Cx3cr1-YFPcreERT2Ndufs4flox/flox mice were generated by crossing B6.129P2(Cg)-Cx3cr1tm2.1(cre/ERT2)Litt/WganJ52 with mice with conditional alleles of the Ndufs4 gene (exon 2 flanked by loxP sites)42. Wild-type C57BL/6 female mice were purchased from Charles River. EAE induction was performed by immunization with MOG35–55, as previously described13.

To induce RFP expression in Cx3cr1-YFPcreERT2R26tdTomato mice, tamoxifen (TAM) (Sigma-Aldrich) was diluted in corn oil (Sigma-Aldrich) at a concentration of 25 mg ml−1 and administered daily through intraperitoneal (i.p.) injections for 5 days at 0.125 mg TAM per g body weight (2.5 mg per 20 g mouse). After a washout period of 28 days to allow peripheral myeloid cells to be replaced de novo (and therefore lose RFP expression)54, EAE was induced.

To knockout Ndufs4 expression in CX3CR1+ cells, TAM was diluted in corn oil and administered through daily i.p. injections for 5 days (as above) starting 1 week after EAE onset to Cx3cr1YFPcreERT2−/−Ndufs4flox/flox (Ndufs4 WT) and Cx3cr1-YFPcreERT2+/−Ndufs4flox/flox (Ndufs4 KO) mice. EAE induction was performed on female mice aged 8–20 weeks. Healthy controls included both female and male mice. For tamoxifen and small molecule treatments of EAE mice, mice were randomly assigned to each treatment group on the first day of treatment (one week after onset), so that the mean baseline EAE score of each group was not statistically different. In all other occasions, mice and samples were randomly selected by blind assessors. Investigators assessing in vivo (for example, behavioural), ex vivo (for example, pathological), and in vitro (cell cultures) outcomes were blinded to group allocation. Data from published13 and preliminary work has been used to obtain estimates of standard deviation (s.d.) and population distribution. Sample size calculations were carried out for an effect size of 0.5, 80% power, 5% level of significance, and for highest s.d. of the outcomes assessed.

Body weight and EAE clinical score were recorded daily. At 7–16 days after immunization, mice developed the first clinical signs of disease (EAE onset), and at 3 days after disease onset, they reached the acute phase of disease (A-EAE). EAE scoring (0, healthy; 1, limp tail; 2, ataxia and/or paresis of hindlimbs; 3, paralysis of hindlimbs and/or paresis of forelimbs; 4, quadriplegia; 5, found dead) was adapted to include 0.25 increments in case of mice exceeding the conventional score, and a score of 4.5 was given to mice that reached the severity limit of our licence and had to be culled humanely. During the experiments, mice were housed in ventilated cages, softly lit and subject to a light dark cycle with a relative humidity of 45–65% and at a temperature range of 20–24 °C.

scRNA-seq

For Cx3cr1-YFPcreERT2R26tdTomato scRNA-seq experiments, mice were deeply anaesthetized with an i.p. injection of ketamine (10 mg ml−1, Boehringer Ingelheim) and xylazine (1.17 mg ml−1, Bayer) and perfused with ice-cold artificial CSF (aCSF). After perfusion, spinal cords were extracted from the spinal columns using a 5 ml syringe filled with ice-cold aCSF. Spinal cords were mechanically dissociated in a glass Dounce tissue homogenizer with 6 ml of homogenization buffer (aCSF plus 10 mM HEPES (Sigma-Aldrich), 1% bovine serum albumin (BSA) (Sigma-Aldrich), 1 mM EDTA (Thermo Fisher Scientific), 10 mg ml−1 of DNase (3000U, Roche) and 40 U μl−1 of RNase inhibitor (Invitrogen)). After tissue dissociation, the suspension was filtered through a pre-wet 40 μM strainer and the homogenizer rinsed with 2 ml of homogenization buffer. The samples were then transferred to 15 ml Falcon tubes and 2.7 ml of 90% Percoll (GE Healthcare) in 10× PBS (Thermo Fisher Scientific) was added to each sample to remove myelin and debris. The 15 ml Falcon tubes were then inverted ten times gently and the samples were centrifuged at 800g for 20 min at 4 °C with a brake speed of 0. Myelin debris visibly layered at the surface was carefully removed with a pipette. Ice-cold buffer (95% autoMACS rinsing solution (Miltenyi Biotec) and 5% MACS BSA (Miltenyi Biotec)) was added to the samples (to fill the 15 ml tubes) and the samples were centrifuged at 800g for 5 min at 4 °C to remove the remaining Percoll. The pellets were resuspended in ice-cold aCSF (1 ml + 7 ml) and then further centrifuged at 800g for 5 min at 4 °C. Pelleted cells were resuspended in 200 μl of FACS buffer (Cell Staining Buffer (BioLegend) plus 10 mM HEPES, 1% BSA, 1 mM EDTA, 10 mg ml−1 of DNase and 40 U μl−1 of RNase inhibitor) and 7-AAD live/dead stain (Thermo Fisher Scientific) added at a concentration of 1:50.

The samples were sorted using a BD FACS Aria III cell sorter set to 3-way purity with a 100 μm nozzle at 20 psi using the FACS gating strategy summarized in Extended Data Fig. 1b. In brief, live cells gates were set based on the unstained WT sample, 7-AAD-stained WT sample treated with DMSO (positive control for cell death) and non-TAM-treated Cx3cr1-YFPcreERT2R26tdTomato samples (to control for bleed-through). FACS-isolated cells were sequenced at a single-cell resolution using the v3 10x Genomics Chromium Single Cell 3′ Solution. Up to 18,000 cells per sample were loaded into each well and the resultant libraries were sequenced on the NovaSeq 6000 (Illumina) system to depth of at least 50,000 reads per cell as calculated by Cell Ranger.

For Nd6 and Ndufs4 scRNA-seq experiments, mice were perfused with actinomycin D (5 mg ml−1, Sigma-Aldrich) and triptolide (10 mg ml−1, Sigma-Aldrich) in ice-cold HBSS (Gibco). Next, spinal cords were extracted, mechanically disassociated and incubated in digestion buffer (collagenase (2 mg ml−1, Sigma-Aldrich), dispase (0.2 mg ml−1 MP Bio), DNase (0.1 mg ml−1, New England), actinomycin D (5 mg ml−1, Sigma-Aldrich), triptolide (10 mg ml−1, Sigma-Aldrich) and anisomycin (27.1 ml ml−1, Sigma-Aldrich)) for 15 min at 0 rpm, and additional 15 min at 300 rpm at 37 °C. After digestion, the tissue suspension was filtered through 100 μm strainer. Next, 9 ml of HBSS were added, and debris was removed from each sample using debris removal solution (Miltenyi Biotec) according to the manufacturer’s protocol. After debris removal, cells were counted, and 16,000 cells per sample were loaded for library preparation.

Before mapping and quantification, the quality of the samples was assessed using fastQC v.0.12.11, applied on raw files; the outputs were summarized using multiQC (v.1.14)55. All of the samples were aligned to the GRCm38 (Ensembl 93 from pre-built Cell Ranger reference (v.3.1.0)) reference genome using Cell Ranger (v.7.1.0)56, including intronic reads. The resulting feature–barcode raw matrices were loaded in Seurat v.4.3.0.1 and consolidated into one Seurat object. Next, cells were filtered using the following thresholds: 400 < number of genes/features < 8000; 800 < number of UMIs < 35,000; 3% < proportion of reads incident to ribosomal genes < 30%; proportion of reads incident to mitochondrial genes < 5%. All of the thresholds were determined on the basis of distributions of respective values, assessed on violin plots, and default UMAPs. The expression levels were normalized using the LogNormalize method from Seurat with the default parameters. Variable features were found using vst57; the top 5,000 most variable features were retained for subsequent steps of the analysis. The expression levels were then scaled and centred. Deterministic PCA (prcomp) was performed on the variable features. The first 50 principal components were retained for subsequent analyses. Harmony58 was used to correct batch effects assessed on the Illumina flow cell batches (θ = 2). UMAPs were then calculated on the batch-corrected, normalized expression matrix using 50 Harmony dimensions. The optimal parameters for determining partitions/clusters on the given dataset were determined using the ClustAssess pipeline59. The shared nearest-neighbour graph was created using FindNeighbors. Clusters were determined using the FindClusters function with the Louvain clustering algorithm, with a resolution parameter of 0.3. The cluster-specific marker genes were determined both using cluster versus complement and all pairwise cluster versus cluster approaches, using the FindMarkers function. The differential expression analysis was performed using the Wilcoxon test, with additional criteria on log2[FC] and Benjamini–Hochberg-adjusted P value. The enrichment analysis of the marker genes was performed using g:profiler.

For the RNA velocity analysis24, loom files were generated using velocyto (v.0.17.17) on all quantified genes. scvelo (v.0.2.5) was used to perform the analysis (that is, filtering, normalization, computing first and second-order moments, and subsequently estimating velocity stochastically). All cluster identification and velocity analyses were performed in R (v.4.2.3) on high-memory Linux servers.

Pyscenic60 was used with the metadata made available for the Mus musculus (mm10) reference genome, making use of the docker container to generate a loom object from our Seurat object. We used loompy61 to create a SCope object, subsequently explored using the SCope web application.

Data-derived filters used for the Nd6 dataset comprise number of genes/features < 2,500; number of UMIs < 5,000; 5% < proportion of reads incident to ribosomal genes < 20%; proportion of reads incident to mitochondrial genes < 15%. For the Ndufs4 dataset, the filtering thresholds were as follows: number of genes/features < 3,000; number of UMIs < 3,000; 5% < proportion of reads incident to ribosomal genes < 40%; proportion of reads incident to mitochondrial genes < 30%. The Ndufs4 data were batch-corrected using Harmony58 for different pools. For both datasets, ClustAssess59 was used to establish the features used and a stable cluster configuration.

For assigning cell types, the outputs of several methods were corroborated.

The unsupervised angle is based on a cluster-driven method relying on stable ClustAssess partitions at a high resolution (that is, a high number of clusters). The smaller clusters were assigned to cell types based on expressed marker genes, which allowed a more granular assessment of transcriptomic signatures. This approach facilitated the reconstruction of assignments for the original, larger clusters. We used this method to assign cell types to the human single-cell data.

The secondary approach was based on a strict cell-centric method for which we ranked each cell based on its median expression of a manually curated, literature-based, gene set associated with one cell type. The cell type was driven by the highest rank; cells with median expression of 0 for all predefined cell types were left unassigned.

The third angle was based on a lenient voting scheme based on a minimum number of genes (all k combinations out of the predefined n genes were permitted) expressed in any one cell; a minimum average expression value was also used. The implementation of this voting scheme is presented in ClustAssess59. For assigning cell identities, we required an average expression higher than 2, and at least 2 marker genes present.

The following marker genes were used: macrophages: Ms4a7, Ccl7, Ecm1, Mgst, Arg1; microglia: Serpine2, P2ry12, Siglech, Slc2a5, Plxdc2, Sall1, Sparc; Monocytes: Vim, Chil3, Plac8, Ccr2; dendritic cells: Flt3, H2-DMb2, Itgax, Ccl17, Ccr7, Ifitm1; Neutrophils: S100a8, S100a9, Csf3r, Cxcr2, Mmp9, Csf1, Il1r2; CNS-associated macrophages: Mrc1, Lyve1; T cells: Cd3e, Bcl11b, Cd4, Cd8a, Il2ra; NK cells: Ncam1, Gzmb, Prf1, Ifng.

Bulk RNA-seq

For comparative gene expression analysis, data obtained from our hiMG differentiation protocol were compared with previously available datasets of primary human microglia isolated from post mortem tissue and microglia derived from human iPS cells62,63. The quality checking was performed using fastQC v.0.12.12 on all raw files. The outputs were summarized using multiQC v.1.14. Initial sequencing depths ranged from 21.2 million to 25.7 million reads, aside from MG 3, which had 58.5 million reads. seqtk (v.1.3-r106)64 was used for subsampling reads, without replacement to 22 million reads, for the MG 3 sample65. Subsequently, all of the samples were aligned to the GRCh38.p13 genome using STAR v.2.7.10a (paired-end mode)66. Expression quantification was performed using featureCounts (v.1.6.3)67 against the gtf matching the reference Homo sapiens genome. Post-quantification quality assessment identified one replicate per sample as an outlier; the subsequent analysis was based on three replicates per sample. Next, noisyR (v.1.0.0)68 was used to estimate and remove noise from the count matrix. The filtered raw expression matrix was normalized using quantile normalization69, within the bulkAnalyseR framework. DEGs were identified using edgeR70 and DESeq271. The P values were adjusted using Benjamini–Hochberg multiple-testing correction. The enrichment analysis was performed using all expressed genes, above noise level, as the background set using g:profiler72. A Shiny app was generated using bulkAnalyseR (v.1.1.0)73 to provide a user-friendly interface for the analysis and visualization of the data.

Reanalysis of available RNA-seq human datasets

Raw fastQ files from previous studies27,28 were downloaded from the European Nucleotide Archive using fasterq-dump. The quality checking and mapping leading to the filtered feature–barcode matrices were performed as described above. For ref. 28, cells with a number of genes/feature between 200 and 5,000 were retained; an upper threshold of maximum 20,000 UMI counts was used, and the maximum mt% was set to 5%. For the ref. 27 dataset, cells with less than 4,000 maximum number of genes were retained; an upper threshold of 15,000 was used for the maximum UMIs per cell; 5% maximum proportion of fragments incident to mitochondrial genes; 10% maximum proportion of reads incident to ribosomal genes.

Both datasets were normalized using SCTransform. The ref. 28 dataset was batch-corrected using Harmony on the patient variable, with θ = 2. To find a reliable clustering, on each separate dataset, ClustAssess was used with 20–50 iterations, testing resolution parameters between 0.1 and 1.5 (0.1 increment steps). For the ref. 27 dataset, the top 4,500 highly variable features were selected for subsequent analyses. For ref. 28 dataset, the top 3,500 highly variable features yielded the most stable results. For both, the resolution parameter was set to 0.6.

For assessing homology, we considered the log-transformed mean expression for homologous genes shared between cluster 4 in our mouse dataset and the MAMS clusters identified using the cluster-based cell-type assignment method. Adjusted R2 values were computed by fitting a linear model, and we obtained the values 0.40 for ref. 28 and 0.38 for ref. 27. For identifying enriched gene sets, we used gprofiler74 on markers found using a receiver operator characteristic (ROC) test with a log-transformed fold change threshold of 0.25 and a threshold for the minimum percentage of cells expressed of 0.1. All analyses were performed on R v.4.2.3, on high-memory servers.

Steady-state metabolomics using LC–MS

For ex vivo LC–MS, data analysis was performed as previously described75,76. In brief, mice were deeply anaesthetized with i.p. injection of ketamine–xylazine and culled by cervical dislocation. Spinal cords were extracted from the spinal columns using a 5 ml syringe filled with ice-cold sorting buffer (DMEM no phenol red (Thermo Fisher Scientific) + 25 mM HEPES (Sigma-Aldrich) + 5% dialysed FBS (Thermo Fisher Scientific) + 1 mM EDTA (Honeywell Fluka) + 1× GlutaMAX (Thermo Fisher Scientific)). Spinal cords were mechanically dissociated in a glass Dounce tissue homogenizer with 7 ml of cold homogenization buffer (sorting buffer + 10 mg ml−1 of DNase). After tissue dissociation, the suspension was filtered through a pre-wet 40 μM strainer and the homogenizer was rinsed with 2 ml of homogenization buffer. The samples were then transferred to 15 ml Falcon tubes and 2.7 ml of 90% Percoll (GE Healthcare) in 10× PBS (Thermo Fisher Scientific) was added to each sample to remove myelin and debris. The 15 ml Falcon tubes were then inverted 10 times gently and the samples were centrifuged at 800g for 20 min at 4 °C with a brake speed of 0. Myelin debris visibly layered at the surface was carefully removed with a pipette. Ice-cold buffer (95% autoMACS rinsing solution (Miltenyi Biotec) and 5% MACS BSA (Miltenyi Biotec)) was added to the samples (to fill the 15 ml tubes) and the samples were centrifuged at 800g for 5 min at 4 °C to remove the remaining Percoll. The pellets were resuspended in sorting buffer (1 ml + 7 ml) and then further centrifuged at 800g for 5 min at 4 °C. Pelleted cells were resuspended in 100 μl of sorting buffer and 0.5% SYTOX blue dead cell stain (Thermo Fisher Scientific).

The samples were sorted using a BD FACS Aria III cell sorter set to yield with an 85 μm nozzle using the FACS gating strategy summarized in the Extended Data Fig. 1b. In brief, the live-cell gates were set based on the unstained WT sample, SYTOX-stained WT sample treated with DMSO (positive control for cell death) and non-TAM-treated Cx3cr1-YFPcreERT2R26tdTomato samples (to control for bleed-through). FACS-sorted cells were collected in polypropylene tubes, centrifuged at 800g for 3 min at 4 °C and resuspended in ice-cold PBS 1× (at a concentration of 1 million cells per ml).

An additional 1 ml of PBS 1× was added to each sample and cells were further centrifuged at 800g for 3 min at 4 °C. Pelleted cells were then resuspended in cold metabolite extraction buffer (at a concentration of 1 million cells per ml), moved to an autosampler vial and stored at −80 °C. When all of the samples were collected (before LC–MS analysis), autosampler vials were immerged in a methanol dry ice bath twice (10 min in and 5 min out). A total of 10 μl per sample was pulled together to generate a pulled sample control. All of the samples were then stored at −80 °C for LC–MS analysis.

Metabolite extraction was performed at 1 million cells per ml on frozen cell pellets in ice-cold methanol, acetonitrile and water (5:3:1) and vortexed for 30 min at 4 °C. After centrifugation (18,000g for 10 min at 4 °C), the supernatants were transferred to prechilled autosampler vials (Phenomenex) and stored at −80 °C for subsequent LC–MS analysis. LC–MS analysis was performed on randomized samples using the Vanquish UHPLC system coupled online to a Q Exactive Orbitrap mass spectrometer (Thermo Fisher Scientific) in positive- and negative-ion modes (separate runs). The LC system was fitted with a Kinetex C18 column (150 mm × 2.1 mm, 1.7 μm; Phenomenex); mobile phases and 5 min gradients were used as previously described in detail75. The acquired spectra (.raw) were converted to.mzXML format using RawConverter and the signals were subsequently annotated and integrated using Maven (Princeton University) in conjunction with the KEGG database. Quality-control samples and instrument stability were assessed as described previously76.

Analysis of LC–MS data was performed using MetaboAnalyst v.5.0 (https://www.metaboanalyst.ca/). Features with >50% missing values were removed (ADP-ribose, IMP, UDP-glucose and UDP-GlcNAc were also removed as there were too many missing values). The remaining missing values were estimated using k-nearest neighbours with no data filtering. Data were normalized to the median, followed by log-transformation and auto scaling to obtain normalized arbitrary units. Partial least squares–discriminant analysis (PLS-DA) was performed on the whole LC–MS dataset to generate a three-dimensional score plot among selected components. A hierarchical clustering heat map was generated on normalized data using autoscale features and the Ward clustering method. Statistical analysis of metabolites found to be significantly altered among all conditions was performed using one-way ANOVA (Graph Pad Prism 9 for macOS, GraphPad) followed by an unpaired t-test for multiple comparisons (Supplementary Data 3). Direct comparisons between cells belonging to two different disease stages were performed to generate volcano plots (fold change threshold: 2.0; P-value threshold: 0.1). Correlation analysis against a given feature or stage was performed with the PatternHunter tool of MetaboAnalyst using Pearson r as distance measure.

LD-REIMS analysis

Data were first preprocessed using a workflow developed in-house and described in detail previously77. Concisely, a signal-to-noise-ratio-based peak detection algorithm was used to reduce the three-dimensional raw data into lists of detected peaks at each pixel. After the effective data resampling, further reduction in data size was achieved through a semi-supervised region of interest (ROI) selection procedure, which performed automatic annotation of pixels to be either sample or background related according to some preliminary manual selection. This ROI mask file generated was then loaded in the nominal next step of the workflow, which aligned and mass corrected all of the detected peaks through a two-step algorithm using user-defined references masses (m/z = 174.0408, 255.2324, 766.5392 in this case) as tissue-unrelated peaks were filtered by the custom-built SPUTNIK package. As a result, a data cube of dimension M × N was produced for each imaging run, where M is the total number of pixels and N is the length of the common mass axis, which was also shared between all runs by considering only the masses that are common in all datasets.

To probe the metabolic differences between varying regions and between images in an untargeted fashion, unsupervised colocalization analysis based on ranked intensity correlations was performed on the preprocessed data. This effectively paired up spectral features that produced statistically significant spatial correlation, of which the relative intensity distributions can be overlayed for visual inspection. All paired features extracted this way were then manually identified through literature search and online databases (HMDB)78, then visualized on a circle plot.

To facilitate multivariate supervised analysis, further refinement of the preprocessed data was performed by computing a weighted average of structural similarity index (SSIM) and multi-SSIM for images of every spectral feature with a representative reference image; those that produced a weighted average below a perceptually defined threshold were deemed to be belonging to either background or isotopes and hence removed. The remaining data cube was then visualized by a spectral unmixing algorithm79 into (k = 2) components (corresponding to tissue and background in an image) of which the abundances can be spatially mapped to give high-contrast images. Spectra from certain anatomically important regions (such as infiltrated white matter, as seen on haematoxylin and eosin staining) were then labelled and extracted by means of manual ROI selection on the abundance image. Weighted logistic regression classification models were then built on the labelled spectra and their performance was assessed using leave-one-out cross-validation between groups. Model refinement through feature selection was then performed based on recursive feature elimination. The validity of the selected features was verified by using them to train another logistic regression model, and then comparing the ROC curves of the new and original models and their respective areas under the curve. These features were considered to be significant in driving the separation between anatomical and pathological features within the imaging data and were identified according to the same procedure as before.

In vitro microglial cultures

All cell lines were tested for mycoplasma contamination on a regular basis using a PCR-based test and were found to be negative for mycoplasma. BV2 microglia were cultured as previously described13. In brief, BV2 microglial cells (gift from A. Tolkovsky) were kept in culture and expanded in expansion medium (DMEM high glucose (Gibco), 2% FBS (Thermo Fisher Scientific), 1% penicillin–streptomycin (Thermo Fisher Scientific)). After reaching 70% confluency, BV2 cells were washed with warm PBS 1×, detached with trypsin-EDTA 0.05% (Thermo Fisher Scientific) for 3 min followed by trypsin neutralization in expansion medium.

Primary microglia were obtained from mixed glial cultures from cerebral cortices from postnatal day 1–3 (P1–P3) WT and Nd6 unsexed pups. In brief, primary microglia were obtained by peeling off the meninges and cortices were dissected, tissue mechanically dissociated and transferred to the enzyme solution (5 ml Neurobasal-A, Gibco; 6.19 mg papain, Worthington Biochemical) at 37 °C for 15 min. Next, enzyme digestion was stopped by adding microglia medium (DMEM high glucose (Gibco), 10% fetal bovine serum (FBS, Gibco), 1% penicillin–streptomycin (Gibco)), the tissue suspension was then centrifuged at 300g for 5 min and the supernatant was removed and incubated in disassociation solution (1 ml Neurobasal-A, Gibco; 5 μl 1 mg ml−1 DNase). After 2 min, the reaction was stopped by adding microglial medium, the cell suspension was centrifuged at 300g for 5 min, resuspended and plated in a 10 cm Petri dish in microglial medium. The next day, cellular debris was removed, and the mixed-glia culture was incubated in microglial medium with 20% L-929 conditioned medium obtained as previously described13. After 2 weeks, microglia were attached to the bottom of the Petri dish, whereas astrocytes formed an upper monolayer. Primary microglial cells were then separated from astrocytes by mild trypsinization, as described previously80. After astrocytes were removed, cells were scraped in ice-cold PBS, and centrifuged at 300g for 5 min at room temperature. Pelleted cells were resuspended, counted and plated in microglia medium with 20% L-929. Cells were allowed to rest for 48 h before experiments were performed.

hiMGs were derived from human iPS cells, line HPSI0214i-wibj_2. Human iPS cells were cultured in mTeSR1 (Stem Cell Technologies) on hESC-qualified Matrigel-coated plates (Corning). Human iPS eclls were routinely passaged using Accutase cell detachment solution (Stem Cell Technologies). Using a protocol with minor modifications, microglia were differentiated from iPS cells through embryoid bodies (EBs) and primitive macrophage precursors (PMPs)81,82. In brief, iPS cells were dissociated to single cells with Accutase cell detachment solution and plated at 10,000 cells per well in 96-well ultra-low attachment plates (Corning) in 100 μl of microglia embryoid body induction media (MEBIM) supplemented with the ROCK inhibitor StemMACS Y27632 (Miltenyi Biotec) before centrifugation at 300g for 3 min at room temperature. Plated cells were then stored in a 5% CO2 incubator set to 37 °C. EBs were cultured for 4 days with a half medium (that is, 50 μl) exchange with MEBIM occurring after 2 days. When performing the half-medium change, small-molecule concentrations were doubled to account for the medium already present in the wells. After 4 days, 16 EBs were plated per well of tissue-culture treated six-well plates in 3 ml of microglial differentiation medium A (MDM-A). From this point forward, 2 ml medium was exchanged every 5–6 days. When performing the two-third medium change, small-molecule concentrations were doubled to account for the medium already present in the wells. After 14 days, PMPs were collected from suspension during a medium exchange. Collected medium was first filtered through a 40 μm nylon membrane then centrifuged at 500g for 5 min at room temperature. Pelleted cells were resuspended in an appropriate volume of microglial differentiation medium B (MDM-B) and a cell count was performed using Trypan Blue (Thermo Fisher Scientific) exclusion. Cells were plated at a density of 180,000–200,000 cells per cm2 in six-well plates in 1.5 ml of MDM-B. PMPs adhered to the uncoated tissue culture-treated plastic within 24 h. The maturation of the PMPs into microglia occurred over 14+ days, with half of the medium exchanged every 3 days. When performing the half-medium change, small-molecule concentrations were doubled to account for the medium already present in the wells. All cytokines and growth factors were obtained from Peprotech. After a minimum of 14 days of maturation, microglia-like cells were incubated in 2 ml of Accutase cell detachment solution for 20–30 min in a 5% CO2 incubator set to 37 °C. After 20–30 min, any remaining attached cells were detached with gentle pipetting. Microglial-like cells were centrifuged at 500g for 5 min at room temperature. Pelleted cells were resuspended in an appropriate volume of MDM-B and a cell count was performed using Trypan Blue exclusion. Cells were allowed to rest for 48 h before experiments were performed.

In vitro cytotoxicity assay

The cytotoxicity of the small molecules on microglia was measured using the CellTox Green Cytotoxicity Assay (Promega) according to the manufacturer’s instructions. BV2 microglia and primary microglia were seeded at 100,000 cells per well and hiMG at 75,000 cells per well in black 96-well clear-bottom plates in their respective experimental media (BV2: DMEM and 1% penicillin–streptomycin; primary microglia: DMEM, 10% FBS, 20% L-929 medium and 1% penicillin–streptomycin; and hiMG: advanced DMEM/F12, 1% N2 supplement, 1% GlutaMAX, 0.1% 2-mercaptoethanol, 10 ng ml−1 GM-CSF, 50 ng ml−1 TGFβ1 and 100 ng ml−1 IL-34). For BV2 microglia, 24 h after plating the cells were stimulated with LPS (100 ng ml−1, serotype EH100, Enzo) and IFNγ (10 ng ml−1, Peprotech) for 24 h in experimental medium. For primary microglia and hiMG, the cells were stimulated for 24 h with LPS + IFNγ in experimental medium 48 h after replating. Cells were then treated with either rotenone (1 µM, Sigma-Aldrich), metformin (10 mM, Sigma-Aldrich), dimethyl malonate-DMM (10 mM, Sigma-Aldrich), disodium malonate-DSM (10 mM, Sigma-Aldrich), 4-octyl itaconate-4-OI (125 μM, StemCell Technologies) or S1QEL1.1 (0.1 μM, Caymen Chemical) individually or in combination for 30 min.

After treatment, the medium was removed, and the cells were incubated with Hoechst (1 μg ml−1) at room temperature for 10 min followed by fluorescence measurement at 358 nm excitation and 461 nm emission on the Tecan Infinite M200 Pro plate reader. Next, the cell lysis solution was added to the designated wells for 30 min at room temperature as a positive control for cell death. After 30 min, the cells were incubated with the CellTox reagent for 15 min at room temperature and shielded from light. The fluorescence was then measured at 485–500 nm excitation and 520–530 nm emission on the Tecan Infinite M200 Pro plate reader. Cell toxicity was determined by first subtracting the medium-only background fluorescence from the Hoechst and CellTox stains, and the ratio of CellTox fluorescence over Hoechst Fluorescence for each well was calculated. The percentage death was then calculated by taking the ratio of CellTox/Hoechst over the average cell death positive control values and multiplying by 100.

Assessment of ROS production in vitro and ex vivo

For the in vitro analysis of mitochondrial superoxide production after RET induction, the MitoSOX Red Mitochondrial Superoxide Indicator (Thermo Fisher Scientific) probe was used according to the manufacturer’s instructions. BV2 microglia and primary microglia were replated at 100,000 cells per well and hiMGs at 75,000 cells per well in 96-well black clear-bottom plates in their respective experimental media (as described above). For BV2 microglia, 24 h after plating the cells were stimulated with LPS + IFNγ for 24 h in experimental medium. For primary microglia and hiMG, the cells were stimulated for 24 h with LPS + IFNγ in experimental medium 48 h after replating. For all cell types, 24 h after stimulation the medium was removed, and RET was induced by adding 100 μl of the RET solution (in the cell types’ experimental medium) containing oligomycin (Sigma-Aldrich) at a final concentration of 2 μM and dimethyl succinate at a final concentration of 20 mM (Sigma-Aldrich). After 15 min at 37 °C (5% CO2 incubator), 100 μl of MitoSOX (2 μM final concentration) and 20 μl of either rotenone (1 μM), metformin (10 mM), dimethyl malonate-DMM (10 mM), disodium malonate-DSM (10 mM), 4-octyl itaconate-4-OI (125 μM) or S1QEL1.1 (0.1 μM) individually or in combination. After 30 min, the medium was removed from the wells and washed once in 150 μl of 1× HBSS (without Ca2+ and Mg2+) (Gibco) and then replaced with 150 μl of 1× HBSS (without Ca2+ and Mg2+). The MitoSOX signal was acquired at 510 nm excitation and 580 nm emission with a Tecan Infinite M200 Pro plate reader.

For the analysis of MitoSOX using flow cytometry, BV2 cells were replated at 0.5 × 106 cells per well in a 6-well plate in 1 ml of experimental media (DMEM (Thermo Fisher Scientific) with 1% penicillin–streptomycin (Thermo Fisher Scientific)). At 24 h after replating, the cells were stimulated for 24 h with LPS + IFNγ in experimental medium. At the end of the 24 h LPS + IFNγ treatment timepoint, RET was induced (as described above). After 15 min, rotenone (1 µM) was added to the well. After 30 min, the medium was collected from the well into a 15 ml conical tube. Then, 500 μl of Accutase cell-dissociation reagent (Thermo Fisher Scientific) was added to the well and the plate was incubated at 37 °C for 5 min. After 5 min, 1 ml of 1× PBS was added to the well and the cells were transferred to the 15 ml conical tube with the conditioned medium. The cells were then pelleted at 300g for 5 min. The supernatant was removed, and the cell pellet was resuspended in 100 μl of MitoSOX Red working solution (1 μl of a 5 mM stock into 1 ml of 1× PBS). The cell suspension was moved to a 1.5 ml microcentrifuge tube and incubated at 37 °C in an Eppendorf Thermomixer comfort machine set to 500 rpm for 30 min protected from light. After 30 min, labelled cells were pelleted at 300g for 5 min. The supernatant was removed, and the cell pellet was resuspended in 300 μl of 1× PBS and transferred to a polystyrene flow cytometry tube (Corning). The samples were kept on ice until acquisition. Then, 10 min before analysis, a working solution of DAPI (1:1,000 in 1× PBS) (Millipore Sigma) was added to the sample (1:100 of working solution).

Flow cytometry gates were set using the FACS-gating strategy summarized in Extended Data Fig. 7b. In brief, gates were set on DAPI- and MitoSOX-Red-positive controls and >200,000 total events were captured using the BD FACS Aria III cell analyser. The FCS files were exported, and analysis was performed using FlowJo v.10 (BD Biosciences).

Ex vivo ROS production of Cx3cr1-YFPcreERT2R26tdTomato mice was obtained after exposing CNS homogenates to 1 μM rotenone or vehicle solution for 30 min, at 37 °C with constant shaking at 500 rpm. We then stained the samples using 5 μM of the CellROX Deep Red Flow Cytometry Assay Kit (Thermo Fisher Scientific) in sorting buffer for 30 min, at 37 °C with constant shaking at 500 rpm. The samples were immediately analysed by flow cytometry on the Aria Fusion (BD Biosciences) system using the FACS gating strategy summarized in the Extended Data Fig. 1b with Sytox Blue (Thermo Fisher Scientific, S34857) to identify live cells. Analysis was performed using FlowJo v.10 (BD Biosciences) in the Sytox-Blue-negative, RFP+YFP+ population (microglia) and RFPYFP+ population (predominantly consisting of infiltrating myeloid cells). Histograms of far-red fluorescence corresponding to the CellROX signal were assessed, and geometric means were extracted.

Mitochondrial membrane potential in vitro and ex vivo

Cells were replated in their respective experimental medium (described above) in a 96-well black clear-bottom plate (Thermo Fisher Scientific) at a density of 0.5 × 105 cells per well. Then, 24 h after replating BV2 microglia and 48 h after replating hiMG, cells were stimulated with LPS + IFNγ in their respective experimental medium. After 24 h of stimulation, RET was induced as described above. After 15 min at 37 °C, cells were stained with either tetramethylrhodamine methyl ester perchlorate-TMRM (Cambridge Bioscience, final concentration 50 nM, for BV2 microglia) or tetramethylrhodamine ethyl ester-TMRE (Abcam, final concentration 200 nM, for hiMG) and rotenone (1 μM) was added to the wells and the plate was incubated for 30 min at 37 °C. For the TMRM probe, the signal was acquired at 548 nm excitation and 575 nm emission on a Tecan Infinite M200 Pro plate reader. For the TMRE probe, the medium was removed from the wells and the cells were washed twice with 100 μl of 1× PBS/0.2% BSA. After washing, 100 μl of 1× PBS/0.2% BSA was added to the well and the signal was acquired at 549 nm excitation and 575 nm emission on the Tecan Infinite M200 Pro plate reader.

The ex vivo mitochondrial membrane potential of Cx3cr1-YFPcreERT2R26tdTomato mice was obtained by staining brain and spinal cord homogenates using 200 nM of the mitochondrial membrane potential indicator MitoView633 (Biotium) in sorting buffer for 30 min, at 37 °C with constant shaking at 500 rpm. The samples were immediately analysed using flow cytometry on the Aria Fusion (BD Biosciences) system using the FACS gating strategy summarized in the Extended Data Fig. 1b with Sytox Blue (Thermo Fisher Scientific, S34857) to identify live cells. Analysis was performed using FlowJo v.10 (BD Biosciences) in the Sytox-Blue-negative, YFP+ and RFP+ population (microglia). Histograms of far-red fluorescence corresponding to the MitoView633 signal were assessed, and geometric means were extracted.

The ex vivo mitochondrial membrane potential of CD45+Cd11b+ myeloid cells from WT and Nd6 mice was obtained by staining brain and spinal cord homogenates on ice in sorting buffer using Brilliant-Violet-510-conjugated antibodies against CD45 (BD Horizon, 1 μl per 1 million cells), FITC-conjugated antibodies against CD11b (BD Pharmingen, 1 μl per 1 million cells) and Zombie Red (BioLegend, 1:500). The samples were then centrifuged at 500g for 5 min at 4 °C. Sorting medium was aspirated, and cells pellets were resuspended and stained in 200 nM of the mitochondrial membrane potential indicator MitoView633 (Biotium) in sorting buffer for 30 min, at 37 °C with constant shaking at 500 rpm. The samples were immediately analysed by flow cytometry on Aria Fusion (BD Biosciences) using the FACS gating strategy summarized in the Extended Data Fig. 9g with Zombie Red (Thermo Fisher Scientific, S34857) to identify live cells. Analysis was performed using FlowJo v.10 (BD Biosciences) in the Zombie-Red-negative, FITC-positive and Brilliant-Violet-510-positive population. Histograms of far-red fluorescence corresponding to the MitoView633 signal were assessed, and geometric means were extracted.

Microglial cells and SH-SY5Y in vitro co-cultures

SH-SY5Y cells (gift from M. Whitehead) were kept in culture and expanded in growth medium (DMEM-F12 (Gibco), 10% FBS, 1% penicillin–streptomycin). The growth medium was refreshed every 4 to 7 days. After reaching 80–90% confluency, SH-SY5Y cells were washed with warm PBS (Gibco), detached with trypsin-EDTA 0.05% (Thermo Fisher Scientific) for 3 min, trypsin neutralized in growth medium and centrifuged at 400g for 5 min. Next, SH-SY5Y cells were counted and seeded at a density of 5 × 104 cells per well in 12-well plates with differentiation medium (Neurobasal medium (Thermo Fisher Scientific) 2% B27 supplement (Thermo Fisher Scientific), 1% GlutaMAX (Thermo Fisher Scientific) and 10 μM all-trans retinoic acid (StemCell Technologies)). Differentiation medium was replaced every other day until day 9 after seeding for co-culture experiments with BV2 microglia, primary microglia and hiMGs.

BV2 microglia, primary microglia and hiMGs were maintained and cultured as described above. To time the co-cultures, on day 8 (24 h before stimulation) of the SH-SY5Y differentiation, BV2 microglia were replated inside a 0.4-μm-pore size Transwell insert (12 well-size, Millipore) at a density of 1 × 105 cells per Transwell with 1 ml of experimental medium. Transwell inserts were placed in a 12-well plate prefilled with 1 ml of experimental media per well. For primary microglia and hiMG, the cells were replated inside the 0.4-μm-pore size Transwell insert on day 7 (48 h before stimulation) of the SH-SY5Y differentiation. On day 9 of the SH-SY5Y differentiation, the BV2 microglia, primary microglia, and hiMG were stimulated with LPS/IFNγ in their respective experimental media for 12 h.

For all microglial cells, 12 h after the LPS/IFNγ stimulation (day 10 of the SH-SY5Y cell differentiation) RET was induced (as described above) in the Transwell for 15 min at 37 °C. After 15 min rotenone (1 μM) was added to the Transwell. The cells were further incubated for 30 min at 37 °C. After 30 min, the medium was removed from the Transwells with the microglia and subsequently transferred to a 12-well plate prefilled with 1 ml of 1× PBS. The Transwell inserts were washed once with prewarmed 1× PBS and put into co-culture with the SH-SY5Y cells. Then, 1 ml of experimental medium was added to the Transwell. After 6 h of co-culture, the Transwells were removed and 10× phase-contrast images of 3–5 ROIs per well were acquired using an Echo Rebel microscope (Discover Echo). After image acquisition, the cells were washed 1× in PBS and 350 μl of RLT lysis buffer (QIAGEN) was added to each well. The plates were stored at −70 °C until RNA isolation.

For the mitoTEMPO experiments, the SH-SY5Y cells were pretreated with 50 μM of mitoTEMPO 12 h before being placed into co-culture with BV2 microglia. For the S1QEL1.1 experiments, BV2 microglia were treated for 30 min with 0.1 μM of S1QEL1.1 after 15 min of RET induction before being placed into co-culture with the SH-SY5Y cells.

For neurite length analysis, the images were converted to 8-bit and contrast was adjusted so the neurites were easily visible using Fiji v.2.0.0. software. Data were quantified by semi-automatic tracking using the NeuronJ plugin from >15 neurite lengths per ROI.

For the caspase-3 experiments, SH-SY5Y cells were seeded at a density of 5 × 104 cells per well on coverslips in a 24-well plate with differentiation medium. Differentiation medium was replaced every other day until day 9 after seeding for co-culture experiments. On day 8 (24 h before stimulation) of the SH-SY5Y differentiation BV2 microglia were replated inside a 0.4-μm-pore size Transwell insert (24-well size, Millipore) at a density of 1 × 105 cells per Transwell with 500 μl of experimental medium. Transwell inserts were placed into a 24-well plate prefilled with 500 μl of experimental medium per well. On day 9 of the SH-SY5Y differentiation, the BV2 microglia were stimulated with LPS + IFNγ in experimental medium for 12 h. Then, 12 h after the LPS/IFNγ stimulation (day 10 of the SH-SY5Y cell differentiation), RET was induced (as described above) in the Transwell for 15 min at 37 °C. After 15 min, rotenone (1 μM) was added to the Transwell. The cells were further incubated for 30 min at 37 °C. After 30 min, the medium was removed from the Transwells and subsequently transferred to a 24-well plate prefilled with 500 μl of 1× PBS. The Transwell inserts were washed once with prewarmed 1× PBS and put into co-culture with the SH-SY5Y cells. Then, 500 μl of experimental medium was added to the Transwell. After 6 h of co-culture, the Transwells were removed and the SH-SY5Y cells were washed once with 1× PBS and subsequently fixed with 4% paraformaldehyde for 10 min, washed twice with 1× PBS and blocked for 1 h at room temperature in blocking buffer (0.1% Triton X-100 in 1× PBS (PBS-T) with 10% normal goat serum (NGS)). Coverslips were then incubated at 4 °C overnight with antibodies against cleaved caspase-3 (Cell Signaling, 1:1,000) in antibody blocking buffer (PBS-T with 1% NGS). Next, the coverslips were washed three times with PBS-T for 5 min, incubated with the appropriate secondary antibody (1:1,000) in antibody blocking buffer for 1 h at room temperature and then washed with PBS three times for 5 min. Lastly, the coverslips were counterstained with DAPI for 5 min then mounted onto glass slides with ProLong diamond antifade mountant (Thermo Fisher Scientific) and sealed with nail polish. Images were taken on the Leica DMI400B microscope of three ROIs per coverslip using a ×63 objective. Images were batch analysed using Fiji; the pipeline included fluorescent signal identification by minimum error thresholding followed by the measurement of the mean grey value. The values from three ROIs were averaged for each coverslip, and the resulting values were again averaged to generate the mean fluorescence intensity of each condition.

Catalase gene expression analysis was performed using RT–qPCR. At the end of co-culture experiment, total RNA from SH-SY5Y cells was collected by washing cells with ice-cold 1× PBS and adding 350 μl of RLT lysis buffer (QIAGEN). The samples were then stored at −80 °C until extraction. Total RNA from SH-SY5Y cells was extracted using the RNeasy Mini Kit (QIAGEN) according to the manufacturer’s instructions. For RT–qPCR analysis, equal amounts of RNA were reversed transcribed using the high-capacity cDNA reverse transcription kit (Thermo Fisher Scientific) according to the manufacturer’s instructions. cDNA was then quantified using the NanoDrop 2000c instrument (Thermo Fisher Scientific) and RT–qPCR was performed using the TaqMan Fast Universal PCR Master Mix (2×) (Thermo Fisher Scientific) and TaqMan Gene Expression Assays for catalase (Thermo Fisher Scientific, Hs00156308_m1), and 18S (Thermo Fisher Scientific) was used for normalization. The samples were run in triplicates using the 7500 Fast Real-Time PCR System (Applied Biosystems) and analysed using the \({2}^{-\Delta \Delta {C}_{{\rm{t}}}}\) method.

Western blot analysis of OXPHOS complexes

Ex vivo isolated microglia were homogenized in 1× RIPA buffer (Abcam) supplemented with 1× protease and phosphatase inhibitors (Thermo Fisher Scientific). Equal protein amounts (5 µg) were resolved by SDS–PAGE on 4–15% Mini-PROTEAN TGX Stain Free gels (30 μl, 10-well comb, Bio-Rad, 4568083) and transferred to polyvinylidene fluoride membranes using the Trans-Blot Turbo Transfer Pack (Bio-Rad, 1704156). The membranes were blocked with 5% milk in 0.1% TBST for 1 h at room temperature. The membranes were then immunoblotted over night at 4 °C with the indicated primary antibodies diluted in blocking buffer (SDHA (Abcam, 1:1,000), total OXPHOS Rodent WB Antibody Cocktail (Abcam, 1:500), NDUFAF1 (Abcam, 1:10,000) and TOMM20 (Santa Cruz, 1:2,000)). The membranes were washed three times with TBST for 5 min and immunoblotted with the appropriate HRP-conjugated secondary antibodies: anti-rabbit (Cell Signaling, 7074, 1:2,000) or anti-mouse (Cell Signaling, 7076, 1:2,000). Membranes were washed three times with TBST for 5 min, then imaged on the ChemiDoc MP Imaging System (Bio-Rad), soaked in Clarity Western ECL Substrate (1705061). Uncropped and unprocessed scans of the blots are provided in the Extended Data Fig. 6.

Protein band densitometry was measured in ImageJ. An identical selection frame was used to determine the mean grey value (MGV) measurement of protein bands and background areas above and below each band. Band intensity (BI) was determined using the formula: BI = (255 − MGV protein band) − [(255 − MGV background above) + (255 − MGV background below)]/2, and values were normalized to TOMM20.

Mitochondrial DNA copy number

We isolated 30,000 cells per sample using the FACS gating strategy summarized in the Extended Data Fig. 1b and centrifuged them at 500g for 5 min at 4 °C. The sorting medium was aspirated, and cells pellets were used for total DNA isolation using the DNeasy Blood & Tissue Kit (Qiagen) according to manufacturer’s instructions. Total DNA was then quantified using the NanoDrop 2000c instrument (Thermo Fisher Scientific). For each sample, 10 ng of DNA was loaded into qPCR reactions using the TaqMan Fast Universal PCR Master Mix (Thermo Fisher Scientific, 4352042) and primers for mitochondrial-genome-encoded MT-ND1 (Thermo Fisher Scientific, 4331182, Mm04225274_s1). TaqMan Copy Number Reference Assay Tfrc (Thermo Fisher Scientific, 4458366) primers were loaded into the same reactions for normalization. qPCR was run on the QuantStudio 7 Flex (Thermo Fisher Scientific, 4485701) system, and the resulting files were analysed using the Quant Studio Design & Analysis Software v.2.5.1 using the \({2}^{-\Delta \Delta {C}_{{\rm{t}}}}\) method.

Metabolic flux analysis with Seahorse

Metabolic flux measurements were performed on myeloid cells isolated ex vivo from Cx3cr1-YFPcreERT2R26tdTomato mice and on primary WT and Nd6 microglia cultured in vitro.

For the ex vivo metabolic flux analysis, a Seahorse calibration plate was rehydrated with calibration buffer (Agilent) overnight in a non-CO2 incubator at 37 °C the day before FACS isolation of cells. On the day of FACS, the Seahorse cell culture plate was coated with Cell-Tak (Corning) and washed with distilled H2O. After FACS (using the FACS gating strategy summarized in the Extended Data Fig. 1b), cells were collected in polypropylene tubes, centrifuged at 800g for 3 min at 4 °C and resuspended in 1× PBS (at a concentration of 1 million cells per ml). Cells were then centrifuged at 800g for 3 min at 4 °C and resuspended at 1 million cells per ml in sorting buffer. Cells were seeded on the Seahorse XFe96 plate at 40,000 cells per well (up to six replicates per mouse per cell type) and centrifuged at 200g for 1 min with 0 brake. The culture plate was placed in an incubator at 37 °C for 30 min and then transferred into a non-CO2 incubator for 30 min. The plate was calibrated, and the cells were washed once with mitochondrial assay solution (MAS, made of sucrose 70 mM + mannitol 220 mM + KH2PO4 10 mM + HEPES 2 mM + EGTA 1 mM, in distilled H2O at pH 7.2). Just before loading the plate into the Seahorse XFe96 analyser, cells were put in MAS + 0.2% BSA + 1 nM plasma membrane permeabilizer (Agilent) + 4 mM ADP (pH 7.4). The cell plate was equilibrated, and the baseline was measured. Glutamate (17.5 mM) and malate (17.5 mM) were used to activate complex I. Rotenone (3 µM, well final concentration) was used to inhibit complex I and disodium succinate (10 mM, well final concentration) (Sigma-Aldrich) was used to activate complex II. Antimycin A (4 µM, well final concentration) was used to inhibit complex III and stop mitochondrial respiration. Five reads were obtained (every 4.63 s) at the baseline and after each injection. If cells did not respond to the addition of antimycin A, they were excluded from further analysis. After the Seahorse run, the cells were fixed with 4% PFA for 10 min and stained with DAPI for 2 min. Using the Leica microsystems CTR4000 system to take images at ×10, the OCR values of every plate were normalized based on the mean area of DAPI of four regions of interest. Seahorse analysis was performed using Wave v.2.6.1. Complex I activity was calculated as the percentage change decrease in the mean OCR obtained from two measurements before and after rotenone injection. Complex II activity was calculated as the percentage change decrease in the mean OCR obtained from two measurements before and after succinate injection.

For the in vitro metabolic flux analysis, primary microglia were replated at 100,000 cells per well of a Seahorse XFe24 Microplate in microglia medium and incubated for 48 h to allow the cells to adhere and rest. After 48 h, the cells were stimulated with LPS + IFNγ in experimental medium for 24 h. The day before the assay, a Seahorse sensor cartridge was rehydrated with calibration buffer (Agilent) overnight in a non-CO2 incubator at 37 °C. In each microplate, four wells did not contain cells to serve as background wells, according to the manufacturer’s guidelines. On the day of the Seahorse assay, assay medium (XF DMEM (Agilent) supplemented with XF glucose (10 mM, Agilent), XF glutamine (2 mM, Agilent) and XF pyruvate (1 mM, Agilent)) was prepared and incubated in a bead bath set to 37 °C 1–2 h before the assay. After 24 h of stimulation with LPS + IFNγ, the medium was removed from the microplate containing the primary microglia and replaced with 500 μl of assay medium then placed into a non-CO2 incubator for 45 min. The drugs were then loaded into the ports of the sensor cartridge: port A contained oligomycin to inhibit mitochondrial ATP production (1 µM final concentration in well); port B contained FCCP to induce mitochondrial uncoupling and maximal mitochondrial respiration (1 µM final concentration in well); and port C contained rotenone/antimycin A to inhibit mitochondrial complex I and complex III, respectively, and therefore stop mitochondrial respiration (1 µM final concentration in well). The plate was then placed into the Seahorse XFe24 Analyzer to calibrate. After 45 min, the cells were then subjected to the XFe Cell Mito stress test (Agilent) according to the manufacturer’s instructions. The cell plate was equilibrated, and the baseline was measured. Three reads were obtained (every 8.75 min) at the baseline and after each drug injection. If cells did not respond to the addition of drugs, they were excluded from further analysis. After completion of the assay, the medium was removed and 50 μl of RIPA buffer (10× RIPA (Abcam), 100× protease/phosphatase inhibitors (Sigma-Aldrich), distilled H2O (Corning)) was added to each well. The protein amount of each well was determined using the Pierce BCA Protein Assay for normalization. Analysis was performed using the Wave v.2.6.1. software (Agilent).

Genomic PCR and RT–qPCR

The genotype of the mice was verified using genomic PCR by collecting the ear notch and extracting DNA by adding DNA Releasy (15 µl, Anachem) to each sample. The samples were lysed in a thermal cycler: 65 °C for 15 min, 96 °C for 2 min, 65 °C for 4 min, 96 °C for 1 min, 65 °C for 1 min, 96 °C for 30 s, 4 °C on hold and molecular-biology-grade water (120 µl, Thermo Fisher Scientific) was added to each sample. A total of 20–100 ng of DNA from each sample was added to a master mix containing molecular-biology-grade water (7.5 µl) and BioMix Red (12.5 µl, Bioline) with the appropriate amount of primers. The primer sequences for Cx3cr1cre mice were as follows: C: 5′-AAGACTCACGTGGACCTGCT-3′; mt, R: 5′-CGGTTATTCAACTTGCACCA-3′; WT, R: 5′-AGGATGTTGACTTCCGAGTTG-3′; expected product: mutant, 300 bp; heterozygote, 300 bp and 695 bp; WT 695 bp. The primer sequences for R26tdTomato mice were as follows: mt F: 5′-CTGTTCCTGTACGGCATGG-3′; mt R: 5′-GGCATTAAAGCAGCGTATCC-3′; WT F: 5′-AAGGGAGCTGCAGTGGAGTA-3′; WT R: 5′-GGCATTAAAGCAGCGTATCC-3′; expected products: mutant, 196 bp; heterozygote, 297 bp and 196 bp; WT, 297 bp. The primer sequences for Nd6 mice were as follows: F: 5′-TACCCGCAAACAAAGATCACCCAG-3′; R: 5′-TTAGGCAGACTCCTAGAAGG-3′. Primer sequences for ndufs4 mice were as follows: C loxP_A: 5′-AGCCTGTTCTCATACCTCGG-3′; mt R loxP_B: 5′-GCTCTCTATGAGGGTACAGAG-3′: WT R loxP_C: 5′-GGTGCATACTTATACTACTAGTAG-3′; expected products: mutant, 200 bp; heterozygote 200 bp and 150 bp; WT, 150 bp. Amplification reaction parameters were as follows: initial denaturation at 94 °C for 4 min (1×); denaturation 94 °C for 30 s, annealing temperature 65 °C for 30 s, elongation time 72 °C for 40 s (repeated 35×); final extension 72 °C for 8 min, and 4 °C on hold. The same PCR protocol was applied for the genomic analysis of cells isolated by FACS from Cx3cr1-YFPcreERT2Ndufs4flox/flox mice, where DNA was extracted using the DNeasy Blood & Tissue Kit (Qiagen) according to the manufacturer’s instruction. Next, amplified nucleic acids (1 μl) were mixed with master mix (24 μl, see above), and loaded onto a 3% agarose gel (agarose (Bioline), Tris-buffered saline TAE (Sigma-Aldrich), red nucleic acid gel stain (Biotium)) with a TrackIt 100 bp DNA adder (Thermo Fisher Scientific), and run at 90 V for 60 min in 1× TAE electrophoresis buffer. After the run was completed, DNA bands were visualized with ultraviolet light (312 nm) using the ChemiDoc Imaging System (Bio-Rad) and directly photographed.

For gene expression analysis using RT–qPCR, total RNA from primary microglia and SH-SY5Y samples was collected at a given timepoint. Next, cells were washed with 1× PBS, 350 μl of RLT lysis buffer (Qiagen) was added and the samples were stored at −80 °C until extraction. All in vitro samples for RT–qPCR experiments were processed using the RNeasy Mini Kit (Qiagen) according to the manufacturer’s protocol. For RT–qPCR analysis, equal amounts of RNA were reversed-transcribed using the high-capacity cDNA reverse transcription kit (Thermo Fisher Scientific) according to the manufacturer’s instructions. TaqMan Fast Universal PCR Master Mix (2×) (Thermo Fisher Scientific) and TaqMan Gene Expression Assays was used for Il1b (Thermo Fisher Scientific, Mm00434228_m1), Il6 (Thermo Fisher Scientific, Mm00446190_m1), Tnf (Thermo Fisher Scientific, Mm00443258_m1), Inos (Thermo Fisher Scientific, Mm00440502_m1) and CAT (Thermo Fisher Scientific, Hs00156308_m1). For normalization, mouse Actb or 18S were used (Thermo Fisher Scientific). The samples were run in technical duplicates or triplicates using the 7500 Fast Real-Time PCR System (Applied Biosystems) and analysed using the \({2}^{-\Delta \Delta {C}_{{\rm{t}}}}\) method.

Ex vivo NDUFS4 protein analysis

To determine the ex vivo NDUFS4 protein expression in CX3CR1+ cells, brain and spinal cord tissue homogenates were prepared for FACS staining as described in the ‘Steady-state metabolomics using LC–MS’ section (up to obtaining unlabelled cells). The unlabelled cell suspension was then stained with a CD16/CD32 Fc blocking antibody (BD Biosciences, 567021, 1:50) for 10 min on ice. After incubation, 1 ml of cell staining buffer (BioLegend, 420201) was added and the cells were centrifuged for 5 min at 400g. The supernatant was then discarded, and the cell pellet was resuspended in 100 μl of cell staining buffer containing antibodies against the cell surface markers CD45 (BD Biosciences, 563891, 1:100, BV510) and CX3CR1 (BioLegend, 149005, 1:500, PE) and the Zombie violet live/dead stain (BioLegend, 423113, 1:100). The cells were incubated for 40 min at 4 °C with gentle agitation of the cells occurring every 10 min. After 40 min, the cells were washed once with 1 ml of cell staining buffer and centrifuged for 5 min at 400g. The supernatant was discarded and 500 μl of 1× fix/perm buffer (BD Biosciences, 562574) was added and the cells were incubated for 40 min at 4 °C with gentle agitation of the cells occurring every 10 min. After the cells were fixed, they were washed twice with 1× perm/wash buffer (BD Biosciences, 562574) and centrifuged for 5 min at 400g. After the last wash, the supernatant was discarded, and the cell pellet was resuspended in 100 μl of 1× perm/wash buffer containing the antibody against the intracellular marker NDUFS4 (Abcam, ab137064, 1:62.5). The cells were incubated for 40 min at 4 °C with gentle agitation of the cells occurring every 10 min. After 40 min, the cells were washed once with 1 ml of 1× perm/wash buffer and centrifuged for 5 min at 400g. The supernatant was discarded, and the cell pellet was resuspended in 100 μl of 1× perm/wash buffer containing the Alexa Fluor 647 secondary antibody (Thermo Fisher Scientific, A21244, 1:100). The cells were incubated for 40 min at 4 °C with gentle agitation of the cells occurring every 10 min. After 40 min, the cells were washed once with 1 ml of 1× perm/wash buffer and centrifuged for 5 min at 350g. The supernatant was discarded, and the cell pellet was resuspended in 500 μl of cell staining buffer and immediately analysed on the Aria Fusion flow cytometer (BD Biosciences). Gates were set using the FACS gating strategy that is summarized in the Extended Data Fig. 10c. In brief, unstained single cells were used to gate the negative population; single stains against CD45, CX3CR1 and NDUFS4 were used to gate the positive populations; a single stain against Zombie violet was used to gate the live cell population (negative). Samples were processed and 20,000–50,000 total events were captured. Analysis was performed using FlowJo v.10 (BD Biosciences).

Small molecules for in vivo studies

The optimal route of administration of DMM (Sigma-Aldrich) and disodium malonate (DSM) (Sigma-Aldrich) that would lead to the lowest blood concentration and highest CNS tissue concentration was determined using a pharmacokinetics study. EAE mice at 1 week (7 days) from clinical disease onset were separated into six treatment groups (n = 2 per group): (1) i.p. DMM (160 mg per kg); (2) i.p. DSM (160 mg per kg); (3) oral (drinking water) DMM (1.5 w/v%); (4) oral DSM (1.5 w/v%); (5) i.p. saline; and (6) oral saline. i.p. injections: mice received 5 daily i.p. injections of DMM, DSM or saline, were allowed to rest for 30 min, then the blood was collected through the tail vein. For oral administration, mice were given 5 days of ad libitum access to DMM-, DSM- or saline-treated water; blood was collected through the tail vein daily for 4 days. Obtained blood was centrifuged at 10,000g for 5 min and the purified plasma was collected and stored at −80 °C until analysis. At the end of the treatment period, the brain and spinal cord was extracted from the mice and flash-frozen in liquid nitrogen. The frozen tissue was then stored at −80 °C until analysis.

Malonate from blood and tissue was quantified as previously described83. In brief, 20–30 mg of tissue was homogenized in 25 μl mg−1 extraction buffer (50% methanol, 30% acetonitrile and 20% water) and 1 nmol 13C3-malonate as an internal standard, using prefilled bead mill tubes (Thermo Fisher Scientific) in a Precellys 24 homogenizer (Bertin Instruments) (6,500 rpm; twice for 15 s, cooling on ice in between). Then, 10 μl blood was added to 450 μl extraction buffer containing the internal standard in a microcentrifuge tube and vortexed. All of the samples were subsequently incubated at −20 °C for 1 h before centrifugation at 17,000g for 10 min at 4 °C. The supernatant was transferred to a microcentrifuge tube and centrifuged under the same conditions. The resulting supernatant was transferred to precooled MS vials and stored at −80 °C until analysis by LC–MS/MS.

At 1 week from the onset of clinical signs of EAE, mice were randomized to receive treatment or control substances. These included: saline 0.9%, metformin (Sigma-Aldrich, 100 mg per kg), 4-OI (StemCell Technologies, 50 mg per kg), DMM (Sigma-Aldrich, 160 mg per kg) or DMM (160 mg per kg) + metformin (100 mg per kg). Mice were given daily i.p. injections with the treatment or control substances until the end of the study (30 days after immunization).

Mouse and human tissue pathology

At the end of the study period, EAE mice were deeply anesthetized with an i.p. injection of ketamine (10 mg ml−1, Boehringer Ingelheim) plus xylazine (1.17 mg ml−1, Bayer) and transcardially perfused with saline (0.9% NaCl + 5 M EDTA) for 5 min followed by 4% PFA for an additional 5 min. The spinal columns were isolated and post-fixed in 4% PFA at 4 °C overnight. The next day, the spinal cords were extracted, washed in 1× PBS and immersed in a 30% sucrose (in PBS) cryoprotectant solution for at least 72 h at 4 °C. The spinal cords were then embedded in optimum cutting temperature medium in a plastic mould and frozen on 2-methylbutane (Thermo Fisher Scientific) on dry ice. The tissue blocks were cryo-sectioned (30 μm for Cx3cr1-YFPcreERT2R26tdTomato EAE mice, 10 μm for all of the other EAE experiments) using a cryostat (Leica, CM1850) with a microtome blade (Feather) onto Superfrost Plus slides (Thermo Fisher Scientific). The sections were stored at −80 °C until use.

Human tissue was obtained from the Multiple Sclerosis and Parkinson’s Tissue Bank (Imperial College London) under ethical approval (IRAS reference: 279989). Autopsy flash-frozen brain tissue from one case with secondary progressive MS (44 years old, male) was cryo-sectioned with a 10 μm thickness using a cryostat (CM1850, Leica) with a microtome blade (A35, Feather). The sections were stored at −80 °C until use.

IHC and Bielschowsky staining

For the RFP and YFP double immunohistochemistry (IHC), sections were washed with 1× PBS + 0.3% Triton X-100 (Sigma-Aldrich) and incubated in Bloxall (Vector Laboratories) for 10 min, washed three times for 5 min in 1× PBS + 0.3% Triton X-100, then blocked in 1× PBS + 0.3% Triton X-100 + 10% NGS (Abcam) for 1 h. Antibody incubation was performed overnight at 4 °C with anti-RFP antibodies (rabbit, 1:400, Abcam) in PBS + 0.3% Triton X-100 + 10% NGS. The next day, the sections were washed three times for 5 min in 1× PBS + 0.3% Triton X-100 and incubated with biotinylated goat anti-rabbit (1:1,000, Abcam) for 10 min, washed three times for 5 min in 1× PBS + 0.3% Triton X-100, and incubated with avidin biotin complex (ABC, Vector Laboratories) for 30 min according to the manufacturer’s instructions. After washing three times for 5 min in 1× PBS + 0.3% Triton X-100, the RFP staining was visualized using brown 3,3′-diaminobenzidine (DAB) (Vector Laboratories). The sections were then washed in H2O for 5 min, and then once for 5 min in 1× PBS. The sections were put in Bloxall for 10 min, washed three times for 5 min in 1× PBS + 0.3% Triton X-100, blocked with avidin and then biotin, washed three times for 5 min in 1× PBS + 0.3% Triton X-100, and then blocked with 1× PBS + 0.3% Triton X-100 with 10% NGS for 1 h. Antibody incubation was performed overnight at 4 °C with an anti-YFP antibody (chicken, 1:1,000, Abcam). The next day, the sections were incubated with biotinylated goat anti-chicken (1:1,000, Vector Laboratories) antibodies for 1 h, washed three times for 5 min in 1× PBS + 0.3% Triton X-100 and incubated with ABC for 30 min. After three 5 min washes in 1× PBS + 0.3% Triton X-100, the YFP staining was visualized with blue/grey DAB (Vector Laboratories). The sections were counterstained with haematoxylin (Sigma-Aldrich) for 30 s and washed with distilled H2O. The sections were covered with mounting medium (Ibidi) and a cover slip placed over top (Thermo Fisher Scientific). Images were taken using the Olympus BX53 system and post-processed using Fiji v.2.0.0.

To quantify the amount of oxidative stress in EAE mice using IHC, we analysed GP91-PHOX expression in the spinal cord of EAE mice at 50 days after immunization26. In brief, the sections were air dried at room temperature for 25 min, rinsed with 1× PBS and blocked with Bloxall for 10 min at room temperature, washed three times for 5 min with 1× PBS and incubated using the M.O.M. Immunodetection Kit (Vector Laboratories) according to the manufacturer’s instructions. Next, the sections were incubated overnight at 4 °C with anti-GP91-PHOX antibodies (mouse, 1:200, BD Biosciences) diluted in 1× PBS + 0.3% Triton X-100 with 1% NGS. The next day, the sections were washed three times for 5 min in 1× PBS and incubated with a biotinylated goat anti-mouse secondary antibodies (1:1,000, Thermo Fisher Scientific) for 1 h at room temperature, washed again three times for 5 min and incubated with the ABC kit (Vector Laboratories) for 1 h at room temperature according to the manufacturer’s instructions. After washing three times for 5 min in 1× PBS + 0.3% Triton X-100, staining was visualized by applying DAB (Vector Laboratories). The sections were next washed twice in distilled H2O for 5 min. Finally, the sections were dehydrated using an alcohol gradient: 70% (twice for 2 min), 90% (twice for 2 min), 100% (twice for 2 min) and xylene (twice for 2 min), mounted (Ibidi) and coverslipped (Thermo Fisher Scientific).

For IHC analysis of human post mortem brain tissue, the sections were washed with cold methanol and H2O2 to eliminate endogenous peroxidase activity. After blocking with 10% NGS in PBS solution, the sections were incubated overnight at 4 °C with primary antibodies against PLP (mouse, 1:200, BioRad) and MHC-II (mouse, 1:1,000, BioLegend) in 1.25% NGS in PBS. The next day, the sections were incubated for 1 h and 30 min at room temperature with polymer-based secondary antibodies (Vector Labs). Subsequent reaction product was developed using 3,3′-diaminobenzidine (Vector Labs). The sections were counterstained with haematoxylin, sealed with Pertex Mounting Medium (CellPath) and imaged using the Olympus BX53 microscope with a motorized stage and the Neurolucida software (v.11.07 64-bit, Microbrightfield). White-matter lesions were staged using a robust, previously validated lesion classification system84,85,86. In brief, the CAL was staged according to its neuroanatomical location, cellular density, degree and pattern of demyelination and lesion immunological status. Using these criteria, CAL was defined as a sharply demarcated demyelinating plaque with high MHC-II reactivity specific to the lesion rim, accompanied by sparse myeloid activation within the periplaque.

To assess in vivo axonal loss, spinal cord sections were analysed after Bielschowsky staining, as previously performed13. In brief, the sections were rehydrated in distilled H2O for 10 min at room temperature and placed into the silver nitrate solution (AgNO3) (10%, Sigma-Aldrich) for 20 min. Next, the sections were washed in distilled H2O and incubated with ammonium silver solution (AgNO3 (10%) with NH3OH (32%, Sigma-Aldrich) added until the precipitate became clear) for 20 min at room temperature. After the sections turned brown, the tissue was incubated in NH3OH solution (0.05%, Sigma-Aldrich) for 5 min at room temperature and developed in the ammonium silver solution supplemented with 100 μl of developer solution (formaldehyde, 37%, (Thermo Fisher Scientific), citric acid (Sigma-Aldrich) and nitric acid (70%, Thermo Fisher Scientific) in distilled H2O). After 2 min of staining development, the sections were quickly incubated in the 0.05% NH3OH solution for 5 min and washed three times for 5 min in 1× PBS. Lastly, the sections were dehydrated using an alcohol gradient: 70% (twice for 2 min), 90% (twice for 2 min), 100% (twice for 2 min) and xylene (twice for 2 min), mounted (Ibidi) and coverslipped (Thermo Fisher Scientific).

Quantification of oxidative stress via IHC for GP91-PHOX and axonal loss via Bielschowsky staining in EAE was obtained from analysing n ≥ 12 from spinal cord axial sections per n ≥ 4 biological replicates positive DAB area and axonal loss area respectively of equally spaced axial spinal cord sections outlined using the Olympus BX53 microscope with a motorized stage and Neurolucida software (v.11.07 64-bit, Microbrightfield). Data are expressed as either the percentage of damaged tissue or the area of positive (%) DAB staining per section (±s.e.m.).

Immunofluorescence analysis

For the assessment of cell markers and axonal degeneration using immunofluorescence, slides were dried and washed with distilled H2O and 1× PBS for 5 min each. Tissue was blocked with 1× PBS + 0.3% Triton X-100 + 10% NGS or NDS for 1 h at room temperature. For those sections stained with primary antibodies made in mouse, slides were blocked with M.O.M mouse IgG blocking reagent (Vector Laboratories) in TBS + 0.1% Tween-20 + two drops of IgG diluent for 1 h at room temperature, followed by a block with M.O.M. protein diluent solution for 5 min at room temperature. The following primary antibodies were used (diluted in 1× PBS + 1% NGS or NDS and kept overnight at 4 °C): anti-GFP (chicken, 1:1,000, Abcam), anti-IBA1 (rabbit, 1:500, FUJIFILM Wako; or goat, 1:500, Novus biologicals), anti-CD3 (rat, 1:500, BD bioscience), anti-CD19 (rabbit, 1:500, Abcam), anti-GFAP (mouse, 1:400, Sigma-Aldrich) anti-CX3CR1 (rabbit, 1:300, Novus biologicals), anti-NEUN (mouse, 1:400, Sigma-Aldrich), anti-OLIG2 (goat, 1:500, Novus biologicals), anti-SPP1 (goat, 1:500, R&D Systems), anti-GP91-PHOX (rabbit, 1:500, Proteintech; or mouse, 1:200, BD bioscience) anti-CASPASE3 (mouse, 1:600, Novus biologicals), anti-NDUFS4 (mouse, 1:100, Abcam; or rabbit, 1:500, Novus biologicals), anti-amyloid precursor protein-APP (mouse, 1:200, Sigma-Aldrich), anti-neurofilament heavy polypeptide-NHP (rabbit, 1:500, Abcam). Notably, mouse primary antibodies were also diluted in M.O.M. protein diluent and citric buffer (10 mM) and microwave (three times for 5 min) unmasking was used for the APP staining, as previously described47. The next day, tissue was washed and incubated with the following secondary antibodies (1:1,000, in PBS + 0.1% Triton X-100 + 1% NGS or NDS for 1 h at room temperature): anti-mouse Alexa Fluor 405 (donkey, Thermo Fisher Scientific), anti-chicken Alexa Fluor 488 (goat or donkey, Thermo Fisher Scientific), anti-rabbit Alexa Fluor 488 (donkey, Thermo Fisher Scientific), anti-mouse Alexa Fluor 546 (donkey, Thermo Fisher Scientific), anti-goat Alexa Fluor 647 (donkey, Thermo Fisher Scientific), anti-mouse Alexa Fluor 647 (donkey, Thermo Fisher Scientific), anti-rabbit Alexa Fluor 647 (goat, Thermo Fisher Scientific) and anti-rat Alexa Fluor 647 (goat, Thermo Fisher Scientific). The tissue was counterstained with 1:10,000 DAPI in 1× PBS for 5 min, then washed and embedded with mounting medium for subsequent analysis.

Quantification of Cx3cr1-YFPcreERT2R26tdTomato specificity was performed by acquiring images on the Andor Dragonfly confocal microscope from control (78 days from TAM administration), A-EAE and C-EAE axial spinal cord sections. Data are expressed as the number (±s.e.m.) of IBA1+, CD3+, CD19+, GFAP+ being RFP+YFP+ or RFPYFP+, or RFPYFP cells from three spinal cord sections per mouse from ≥2 independent biological replicates.

SPP1+ microglia (RFP+YFP+) and macrophages (RFPYFP+); NDUFS4+SPP1+ microglia (RFP+YFP+); and SPP1+ microglia (CX3CR1+) were quantified by acquiring images on the Leica microsystems CTR4000 confocal microscope (Leica Biosystems) and expressed as number of cells (±s.e.m.) from three spinal cord sections per mouse from ≥4 independent biological replicates.

Caspase-3+ microglia (IBA1+) from Cx3cr1-YFPcreERT2Ndufs4flox/flox mice and GP91-PHOX+ microglia (IBA1+), astrocytes (GFAP+), oligodendrocytes (OLIG2+) and neurons (NEUN+) from EAE mice that received the small-molecule treatments were quantified from n = 3 images of the posterior columns of spinal cord axial sections acquired through the Leica microsystems DMI4000 microscope (Leica Biosystems) and expressed as the number of cells (±s.e.m.) from ≥4 independent biological replicates.

Quantification of GP91-PHOX+SPP1+ microglia (RFP+YFP+) and macrophages (RFPYFP+); GP91-PHOX+SPP1+ microglia (IBA1+); and microglial arborizations was performed by acquiring images on the Andor Dragonfly confocal microscope. Data are expressed as the number of cells (±s.e.m.) from three spinal cord sections per mouse from ≥4 independent biological replicates. Skeleton analysis was performed as described previously41 using Fiji v.2.0.0. software plugins AnalyzeSkeleton (2D/3D) after applying cut-off criteria = 4.0, pixel radius = 3 and mask weight = 0.6. Data are expressed as the average number of endpoints per cell, and average process length per μm2 from n ≥ 3 ROIs per mouse from n ≥ 3 independent biological replicates (±s.e.m.).

Human tissue immunofluorescence was performed on post mortem brain tissue sections after fixing them in 4% PFA for 10 min. After blocking with 10% NDS in 0.3% Triton X-100 PBS solution (PBS-T), consecutive sections were incubated overnight at 4 °C with primary antibodies against MHC-II (mouse, 1:500, BioLegend) and SPP1 (rabbit, 1:400, Abcam), or antibodies against MHC-II (mouse, 1:500, BioLegend) and NDUFS4 (rabbit, 1:300, Atlas Antibodies) in 1.25% NDS PBS-T. The next day, autofluorescence was quenched using TrueBlack Lipofuscin Quencher (23007, Biotium, 1:20 dilution) according to the manufacture’s recommendations. The sections were then incubated for 1 h at room temperature in 1.5% NDS PBS with the following secondary antibodies (1:250): AlexaFluor 488 donkey anti-rabbit (Thermo Fisher Scientific) and AlexaFluor 546 donkey anti-mouse (Thermo Fisher Scientific). After three PBS washes, the sections were incubated in DAPI (1:1,000) for 5 min, mounted using ProLong Diamond Antifade Mountant (P36961, Invitrogen), then imaged using a whole-slide scanner at ×20 resolution (Zeiss Axioscan Z1).

For in vitro immunocytochemistry, hiMGs were plated onto a glass coverslip (VWR) fixed with 4% PFA (Sigma-Aldrich) for 10 min and washed twice with 1× PBS. Next, cells were blocked with 10% NGS 0.1% Triton X-100 in 1× PBS for 30 min, and incubated overnight with anti-IBA1 (1:500, WAKO) and anti-CSF1R (1:100, Invitrogen) antibodies at 4 °C. The next day, cells were washed three times in 1× PBS, and incubated for 1 h with anti-rabbit Alexa 488 and anti-rat Alexa 546 secondary antibodies (1:1,000, Thermo Fisher Scientific) at room temperature. Cells were counterstained with DAPI (1:10,000, Abcam) in 1× PBS for 5 min and embedded in mounting medium (Fluoromount-G, Invitrogen). The Leica microsystems CTR4000 confocal microscope (Leica Biosystems) was used to take confocal images of the staining.

CyTOF analysis

The first groups of Ndufs4 WT and KO EAE mice (n = 4) and of small-molecule treated EAE mice (n = 6) for CyTOF were deeply anaesthetized with i.p. injections of ketamine (10 mg ml−1, Boehringer Ingelheim) and xylazine (1.17 mg ml−1, Bayer) and then transcardially perfused with ice-cold 1× PBS + 0.1% EDTA. The spinal cords were extracted from the spinal columns using a 5 ml syringe filled with ice-cold sorting buffer (DMEM no phenol red (Thermo Fisher Scientific) + 25 mM HEPES (Sigma-Aldrich) + 5% dialysed FBS (Thermo Fisher Scientific) + 1 mM EDTA (Honeywell Fluka) + 1× GlutaMAX (Thermo Fisher Scientific)). The spinal cords were mechanically dissociated in a glass Dounce tissue homogenizer with 7 ml of cold homogenization buffer (sorting buffer + 10 mg ml−1 of DNase). After tissue dissociation, the suspension was filtered through a pre-wet 40 μM strainer and the homogenizer rinsed with 2 ml of homogenization buffer. The samples were then transferred to 15 ml Falcon tubes and 2.7 ml of 90% Percoll (GE Healthcare) in 10× PBS (Thermo Fisher Scientific) was added to each sample to remove myelin and debris. The 15 ml Falcon tubes were then inverted ten times gently and the samples were centrifuged (Heraeus Megafuge 40R, Thermo Fisher Scientific) at 800g for 20 min at 4 °C with a brake speed of 0. Myelin debris visibly layered at the surface was carefully removed with a pipette. Ice-cold buffer (95% autoMACS rinsing solution (Miltenyi Biotec) and 5% MACS BSA (Miltenyi Biotec)) was added to the samples (to fill the 15 ml tubes) and the samples were centrifuged (Heraeus Megafuge 40R, Thermo Scientific Scientific) at 800g for 5 min at 4 °C to remove the remaining Percoll. Pellets were resuspended in sorting buffer (1 ml + 7 ml) and then further centrifuged at 800g for 5 min at 4 °C.

The cell pellet was resuspended in 100 μl of cold buffer and 10 μl of CD45 positive-selection MicroBeads (Miltenyi Biotec) were added to each sample and mixed well. The samples were incubated at 20 min on ice and gently mixed every 5 min. After incubation, 2 ml of cold buffer was added to wash the samples and the samples were then centrifuged at 300g for 10 min at 4 °C. After centrifugation, the supernatant was removed with a pipette, and the pellet was resuspended in 500 μl of cold buffer. The LS MACS columns (Miltenyi Biotec) were placed in the QuadroMACS Separator (Miltenyi Biotec) on the MACS MultiStand (Miltenyi Biotec) and rinsed three times with 3 ml of cold buffer. After the columns were rinsed, the 500 μl cell suspension was added onto the column. Unlabelled (that is, CD45) cells from the filtrate were collected and the column was washed three times with 3 ml of cold buffer. Finally, 5 ml of cold buffer was added to the column and subsequently flushed with a plunger to isolate the CD45+ cells. The samples were centrifuged at 300g for 5 min at 4 °C. The supernatant was removed, and the cells were resuspended in 50 μl of Maxpar PBS (Fluidigm). A working solution (1:500) of 5 μM Cell-ID Cisplatin (Fluidigm) was prepared and 50 μl was added to the cells (1:1 mix) for a final concentration of 2.5 μM. The cells were incubated for 90 s and the cisplatin was quenched with 5 ml of Maxpar cell staining buffer (Fluidigm). The samples were centrifuged at 300g for 5 min at 4 °C, the supernatant was removed and the cells were resuspended in 1 ml of 1× Fix I buffer (Fluidigm) diluted in Maxpar PBS and left to incubate for 10 min at room temperature. The samples were centrifuged at 1,000g for 7 min at room temperature, the supernatant was removed and the cells were resuspended in 1 ml of 1× barcode perm buffer (Fluidigm) and then recentrifuged at 1,000g for 7 min at room temperature. The supernatant was removed, and the cells were resuspended in 800 μl of 1× barcode perm buffer. Then, 100 μl of a Cell-ID palladium (Pd) unique barcode (Fluidigm) was added to each sample and left to incubate for 30 min at room temperature. The samples were gently mixed after 15 min. After barcoding, the samples were centrifuged at 1,000g for 7 min at room temperature, the supernatant was removed, and the cells were resuspended in 2 ml of Maxpar cell staining buffer. The samples were recentrifuged at 1,000g for 7 min at room temperature, the supernatant was removed and the cells were resuspended in 100 μl of Maxpar cell staining buffer. All of the barcoded samples were combined into one single 5 ml polypropylene tube (BD Biosciences). Each sample tube was washed with 500 μl of Maxpar cell staining buffer, transferred to the single tube (containing the barcoded samples) and centrifuged at 1,000g for 7 min at room temperature. The supernatant was removed, the cells were resuspended in 50 μl of Maxpar cell staining buffer and CD16/CD32 Fc receptor block (BD Biosciences) at 1:25 was added. The sample was left to incubate for 10 min at room temperature. After 10 min, 50 μl of the surface antibody cocktail (Supplementary Data 6) was added to the sample. The sample was mixed with gentle pipetting and incubated for 15 min at room temperature. At 15 min, the sample was gently mixed and left to incubate for an additional 15 min. After surface antibody incubation, 2 ml of Maxpar cell staining buffer was added to the sample and it was centrifuged at 1,000g for 7 min at room temperature. The supernatant was removed, 2 ml of Maxpar cell staining buffer was added to the sample and it was recentrifuged at 1,000g for 7 min at room temperature. The supernatant was removed, and the cells were fixed in 1 ml of 1× Maxpar Fix I buffer for 15 min at room temperature. Then, 2 ml of Maxpar Perm-S buffer (Fluidigm) was added and the sample was centrifuged at 1,000g for 7 min at room temperature. The supernatant was removed, 2 ml of Maxpar Perm-S buffer was added and the sample was recentrifuged at 1,000g for 7 min at room temperature. The supernatant was removed (leaving behind ~50 μl) and 50 μl of the cytoplasmic/secreted protein antibody cocktail (Supplementary Data 6) was added. The sample was gently mixed with a pipette and left to incubate for 30 min at room temperature (at 15 min the sample was gently mixed). After incubation, the cells were washed with 2 ml of Maxpar cell staining buffer, and centrifuged at 1,000g for 7 min at room temperature. The supernatant was removed, 2 ml of Maxpar cell staining buffer was added to the sample, and it was recentrifuged at 1,000g for 7 min at room temperature. The supernatant was removed, and the cells were fixed in 1 ml of 1.6% PFA for 10 min at room temperature. After 10 min, the cells were centrifuged at 1,000g for 7 min at room temperature. The supernatant was removed, and the cells were incubated with 1 ml of Cell-ID Intercalator-IR diluted in Maxpar fix and perm buffer (1:1,000 dilution of 125 μM stock solution) (Fluidigm) for 1 h at room temperature. After incubation, the cells were centrifuged at 1,000g for 7 min at room temperature. The supernatant was removed, 2 ml of Maxpar cell staining buffer was added to the sample, and it was recentrifuged at 1,000g for 7 min at room temperature. The supernatant was removed, 1 ml of Maxpar cell staining buffer was added to the sample, and it was recentrifuged at 1,000g for 10 min at room temperature. At this stage, the supernatant was removed leaving between 50–100 μl of solution covering the cell pellet. The sample was stored for 48 h at 4 °C until acquisition.

The next day, the other groups of Ndufs4 WT and KO EAE mice (n = 4) and of small-molecule-treated EAE mice (n = 6) were processed for CyTOF according to the same protocol. This set of barcoded samples was stored for 24 h at 4 °C until acquisition. On the day of acquisition, the samples were combined into a single tube of eight barcoded samples. The cells were washed in 1 ml of Maxpar Cell Acquisition Solution (Fluidigm) and centrifuged at 1,000g for 7 min at room temperature. The supernatant was removed, and the samples were left on ice until acquisition. Directly before acquisition, the cells were resuspended in 1 ml of 0.1× EQ Four Element Calibration Beads (Fluidigm) diluted in Maxpar Cell Acquisition Solution, passed through a 40 μM filter and then analysed on the Helios mass cytometer using a WB injector. The data files were debarcoded using Fluidigm software (v.7.0.8493) with a Mahalanobis distance of 2.50 (Ndufs4 WT and KO EAE mice) and 8.0 (small-molecule-treated EAE mice) then concatenated and exported as FCS files.

The individual FCS files were imported into FlowJo v.10 (BD Biosciences) and live cells were selected using Boolean gating (the gating strategy is reported in Extended Data Fig. 10j). We first gated time by 140Ce beads to capture all cells. These events were then gated on time by event length followed by gating on viable cells (194 Pt). We next further gated viable cells positively stained with the nucleic acid intercalator iridium (191Ir+193Ir+). Viable cells were then gated to select for myeloid cells (CD45+CD11b+) and lymphoid cells (CD45+CD11b) and exported as individual FCS files to be further analysed in R v.4.1.3. Fcs files were loaded using flowCore (2.6.0). For the small-molecule treated EAE samples, we used the package GdClean_0.0.0.9000 (https://github.com/JunweiLiu0208/GdClean) to estimate and remove gadolinium contamination from the samples. Flow sets were then converted to single-cell experiment objects using CATALYST (v.1.18.1). Using CATALYST, cells were assigned to 1 of 100 grid points on a self-organizing map and metaclustered into 20 groups using the consensus clustering method considering cell-surface features only87. UMAP dimensionality reduction was performed using the CATALYST wrapper around scater runUMAP with the default parameters and data from cell-surface features only. Each of the clusters was assigned to cell types on the basis of the expression of cell-type-specific protein markers (Supplementary Data 7).

All packages are publicly available at the Comprehensive R Archive Network (https://cran.r-project.org), the Bioconductor project (http://bioconductor.org) or their respective GitHub pages. Complete R analysis workflow for CyTOF is available at GitHub (https://github.com/regan-hamel/ComplexI), and CyTOF code and data has been published on Zenodo (https://doi.org/10.5281/zenodo.10510046).

Statistical analyses

Unless otherwise stated, statistical analyses were performed using Graph Pad Prism 9 for macOS (GraphPad Software). EAE scores were compared using two-way ANOVA followed by Bonferroni multiple-comparisons test. Differences between two groups were analysed using unpaired t-tests unless otherwise stated. Differences among more than two groups were analysed using one-way ANOVA followed by Tukey’s multiple-comparison test unless otherwise stated. Data are shown in the text and figures as the mean ± s.e.m. unless otherwise stated P < 0.05 was considered to be significant in all analyses unless otherwise stated. Methods of the other statistical analyses have been described in the relevant sections reported above.

Ethics statement

Animal research has been regulated under the Animals (Scientific Procedures) Act 1986 Amendment Regulations 2012 following ethical review by the University of Cambridge Animal Welfare and Ethical Review Body (AWERB). Animal work was covered by the PPL 80/2457 and PPL PP2135981 (to S.P.). Research on human tissue has been approved by the London—Queen Square Research Ethics Committee (IRAS reference: 279989) and the Cambridge University Hospitals NHS Foundation Trust Research and Development Department.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.



Source link

By AUTHOR

Leave a Reply

Your email address will not be published. Required fields are marked *