Atoshi Banerjee
Yimei Lu
Kenny Do
Travis Mize1,2
Xiaogang Wu
Jingchun Chen- 1Nevada Institute of Personalized Medicine, University of Nevada, Las Vegas, NV, United States
- 2Department of Psychology, Nevada Institute of Personalized Medicine, University of Nevada, Las Vegas, NV, United States
- 3Department of Genomic Medicine, The University of Texas MD Anderson Cancer Center, Houston, TX, United States
- 4410 AI, LLC, Germantown, MD, United States
Microglia are the primary resident immune cells of the central nervous system that maintain physiological homeostasis in the brain and contribute to the pathogenesis of many psychiatric disorders and neurodegenerative diseases. Due to the lack of appropriate human cellular models, it is difficult to study the basic pathophysiological processes linking microglia to brain diseases. In this study, we adopted a microglia-like cellular model derived from peripheral blood monocytes with granulocyte-macrophage colony-stimulating factor (GM-CSF) and interleukin-34 (IL-34). We characterized and validated this in vitro cellular model by morphology, immunocytochemistry, gene expression profiles, and functional study. Our results indicated that the iMG cells developed typical microglial ramified morphology, expressed microglial specific surface markers (P2RY12 and TMEM119), and possessed phagocytic activity. Principal component analyses and multidimensional scaling analyses of RNA-seq data showed that iMG cells were distinct from monocytes and induced macrophages (iMacs) but clustered closer to human microglia and hiPSC-induced microglia. Heatmap analyses also found that iMG cells, but not monocytes, were closely clustered with human primary microglia. Further pathway and relative expression analysis indicated that unique genes from iMG cells were involved in the regulation of the complement system, especially in the synapse and ion transport. Overall, our data demonstrated that the iMG model mimicked many features of the brain resident microglia, highlighting its utility in the study of microglial function in many brain diseases, such as schizophrenia and Alzheimer's disease (AD).
Introduction
Microglia are the primary immune cells of the central nervous system (CNS) that belong to the myeloid lineage. The origins of microglia and their functions in the CNS have been extensively debated since they were first described by Pio del Rio-Hortega in 1932 (Río-Hortega, 1932). In the mouse brain, microglia originate from myeloid progenitors in the yolk sac that migrate to the brain during early embryonic stages before the blood-brain barrier is formed (Ginhoux et al., 2013). A similar pattern of events may occur in humans, but under certain inflammatory conditions, the recruitment of bone marrow-derived progenitors may supplement the microglial population to some extent (Ginhoux and Prinz, 2015). The origins of human microglia are more heterogeneous, with one population from myeloid/mesenchymal origin in the yolk sac and a second population representing a developmental and transitory form of fetal macrophage (possibly monocytes). During normal brain development, microglia proliferate and self-renew continually throughout life to maintain cell numbers, without any contribution from bone marrow-derived macrophages (Elmore et al., 2014). Microglia serve a variety of important functions, including immune surveillance, neurogenesis, synaptic pruning, and phagocytosis (Stansley et al., 2012). Microglia-mediated synaptic loss by phagocytosis strongly correlates with cognitive decline in patients with schizophrenia or AD.
In response to environmental stimuli/stresses, microglia are highly plastic and can be activated into two distinctive phenotypes, the classical pro-inflammatory state (M1) and the alternative anti-inflammatory state (M2) (Ma et al., 2017). Microglial activation and neuroinflammation are hallmarks of many CNS pathologies. In response to lipopolysaccharide (LPS) and interferon (IFN)-γ, microglia are activated to pro-inflammatory M1 state and release cytokines, such as Interleukin (IL)-1β, IL-6, Tumor Necrosis Factor (TNF)-α, and neurotoxic factors. Prolonged hyperactivity of microglia may lead to neuronal apoptosis and brain damage, as seen in schizophrenia (Pérez-Neri et al., 2006; Monji et al., 2013) and AD (Akiyama et al., 2000; Heneka et al., 2015; Wang et al., 2015). Alternatively, in response to IL-4, microglia can be activated to the M2 state and release numerous anti-inflammatory factors, including CD206 and IL-10. M2 microglia promote healing and improve phagocytic capabilities (Cherry et al., 2014).
It is noteworthy that many studies have utilized a mouse model to investigate microglial phenotypes and functions. However, microglial transcriptomic signatures vary between mice and humans (Deczkowska et al., 2018; Friedman et al., 2018; Zhou et al., 2020); most animal-based clinical trials have failed in human diseases, including AD (Franco and Cedazo-Minguez, 2014). Other common approaches to study microglia in human diseases include post-mortem examination of the brain or brain imaging from living patients. Nevertheless, in vitro studies on human primary microglia are challenging due to technical and ethical issues. Human microglial cell lines, such as HMO6 cells (Nagai et al., 2005) and HMC3 cells (Janabi et al., 1995), can be an alternative model for human microglia. However, these cell lines have limited capacity in capturing disease- or patient-specific features of microglia. Furthermore, cell lines may change their cellular phenotypes and functions over time due to genetic manipulations and multiple passages.
Recently, many researchers are experimenting with different reprogramming methods to generate human induced pluripotent stem cell (hiPSC)-derived microglia (Muffat et al., 2016; Abud et al., 2017; Douvaras et al., 2017; Haenseler et al., 2017; Pandya et al., 2017; Takata et al., 2017). While research on hiPSC-derived microglia has been shown to be promising, the field still faces technical challenges, such as long differentiation period, which may not be suitable for high throughput clinical trials, especially for AD patients. Therefore, it is urgent to establish an efficient cellular model representing patient-specific microglia for neurological disease studies.
In this study, we generated and modified a human microglial model, with a 10–14-day conversion from human peripheral blood monocytes, according to previous studies (Ohgidani et al., 2014). The main purpose of this study was to validate the iMG cellular model with morphology, surface markers, functional assays, and RNA-seq profile. Our goals are (1) establish and validate a personalized microglia-like model from human monocytes; (2) develop an in vitro cellular system to study microglial functions, especially phagocytosis and polarization, at molecular and gene expression levels; (3) demonstrate a potential application in psychiatric and neurological diseases, such as schizophrenia and AD. Our data showed that iMG cells mimicked human primary microglia with their typical morphology, surface markers, phagocytic capacity, and plasticity for polarization. With whole-genome transcript profile, we also demonstrated that iMG cells were clustered closely with human microglia.
Materials and Methods
Materials
PBMCs were purchased from AllCells, LLC (Alameda, CA, USA, #PB003F and #PB001). Fresh blood was purchased from Gulf Coast Regional Blood Center. All other reagents were RPMI-1640 without glutamine (vWR, Radnor, PA, USA, #vWRL0106-0500), RPMI-1640 Glutamax (Life Technologies, Carlsbad, CA), Heat-inactivated fetal bovine serum (FBS, Endotoxin <1.0 EU/mL, Rockland Immunochemicals Inc, Pottstown, PA, #FBS-01-0100), Penicillin-streptomycin (P/S) (Gibco), cell culture plates (Corning, NY), Geltrex (Thermo Fisher Scientific Waltham, MA, #A1413202), Trizol (Invitrogen, Carlsbad, CA, #15596018), SuperScript™ IV VILO™ Master Mix (Invitrogen, Carlsbad, CA), High-Capacity cDNA Reverse Transcription Kit with RNase Inhibitor (Thermo Fisher Scientific Waltham, MA, #4374966), 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA), Phosphate buffer saline (PBS) (Genesee Scientific, San Diego, CA), Histopaque-1077 (Sigma Chemical Co., St. Louis, MO), SepMate™-50 tube (Stem cell technologies), IL-34 (R&D Systems Minneapolis, MN, #5265-IL-010/CF), GM-CSF (R&D Systems Minneapolis, MN, #215-GM-010/CF), Lipopolysaccharide (LPS) (Sigma, St. Louis, MO, #L6529), human IL-4 (Biolegend, San Diego, CA, #714904), DMSO (Sigma, St. Louis, MO), Paraformaldehyde (PFA) (Sigma, St. Louis, MO), Fluoromount-G™ Mounting Medium with DAPI (Invitrogen, Carlsbad, CA), human HiLyte™ Fluor 488-labeled Aβ protein (1–42) (Anaspec, #AS-60479-01), and hexafluoroisopropanol (HFIP) (#H8508, Sigma).
Antibodies
Primary antibodies were anti-TMEM119 (Proteintech, #27585-1-AP), anti-P2RY12 (Alomone Labs, #APR-020), anti-Iba1 (Abcam, #ab15690), anti-PU.1 (Abcam, #ab88082), anti-CX3CR1 (Invitrogen, #PA5-19910), and anti-TREM2 (Abcam, #ab209814). Conjugated secondary antibodies were Alexa Fluor 488 (#A11034) and Alexa Fluor 594 (#A10036) from Invitrogen Carlsbad, CA.
Generation of Induced Microglia-Like Cells (iMG cells) From Monocytes
Human peripheral blood mononuclear cells (PBMCs) were isolated from fresh blood by density gradient centrifugation with Histopaque-1077 in SepMate™-50 tubes. iMG cells were generated following the previously reported instructions with some modifications (Ohgidani et al., 2014). Briefly, 1 × 106 cells/ml fresh PBMCs were first resuspended in complete media (RPMI-1640 containing 10% FBS and 1% P/S) in plates pre-coated with 2% Geltrex and incubated overnight at 37°C with 5% CO2. For frozen PBMCs, cells were thawed and washed with complete media before plating. The next day (Day 0), non-adherent cells were carefully removed. The adherent monocytes were cultured with serum-free media (RPMI-1640 Glutamax supplemented with 0.1 μg/ml IL-34, 0.01 μg/ml GM-CSF, and 1% P/S) for 10–14 days. For comparison, induced macrophages (iMacs) were also generated from monocytes by induction with RPMI-1640 Glutamax media containing 10% FBS, 0.01 μg/ml GM-CSF, and 1% P/S for 7 days. During the induction, we replaced fresh media every 3–4 days.
Cell Morphology
During the induction, cellular morphology was monitored with phase-contrast microscopy (Accu-Scope EXI-300, ACCU-SCOPE Inc., Commack, NY, USA). Images were captured with an AU-300-HD digital camera (ACCU-SCOPE Inc.) through CaptaVision PC Imaging Software.
Immunocytochemistry
For immunocytochemistry, iMG cells or iMacs were first fixed with 4% PFA for 10 min and washed thrice with PBS. For intracellular staining, cells were permeabilized with 0.01% tween 20/PBS for 30 min. Cells were then blocked with 5% BSA/PBS for 30 min, followed by incubation with primary antibody (RT for 1 h or 4°C overnight) and then secondary antibody (RT for 1 h at dark). The primary antibodies used in this study were anti-P2RY12 (1:200), anti-TREM2 (1:500), anti-CX3CR1 (1:1000), anti-Iba1 (1:300), anti-PU.1 (1:200), and anti-TMEM119 (1:500) diluted in blocking solution. Depending on the species of primary antibodies, 1:500 Alexa Fluor 488- or Fluor 594-conjugated antibodies were used for staining. Cells were then mounted on a glass slide with Fluoromount G containing DAPI. Images were taken with a confocal laser scanning microscope (Nikon A1R) equipped with the NIS-Elements Advanced Research software. Exposure settings within each set of experiments were kept constant for comparison.
Fibrillary Amyloid-Beta (fAβ42) Preparation and Phagocytosis Assay
The preparation of fibrillary fluorescent Aβ was modified from the previously described procedure (Stine et al., 2003; Jungbauer et al., 2009). Briefly, Human HiLyte™ Fluor 488-labeled Aβ (1–42) was first resuspended as 0.1 mM HFIP and left to dry in a fume hood overnight. The resulting peptide film was stored (desiccated) at −20°C for later use. On the day before the experiment, the films were dissolved in 0.1 mM anhydrous DMSO and sonicated at 22°C for 10 min in a bath sonicator (Bioruptor Pico sonication device). To form fibrillary amyloid-beta 42 (fAβ42), the solution was diluted to 100 μM in 10 mM HCl with vigorous shaking at 37°C for 24 h. iMG cells were incubated with 100 nM fAβ42 in serum-free DMEM/F12 media at 37°C for 4 h. Cells were then quenched with 50 mM NH4Cl for 5 min and fixed with 4% PFA for 5 min at RT. Next, cells were stained with anti-Iba1 and mounted on a glass slide with Fluoromount G containing DAPI. Cells without fAβ42 treatment were used as negative controls for background setting. Fluorescent images were captured as described in the previous section.
Phagocytosis assay was also quantitated by flow cytometry using Sony SH800 (Sony Biotechnology). To this end, cells with or without fAβ42 incubation were first detached with 10 mM EDTA for 5 min and harvested by cell lifter. Cells were then fixed with 4% PFA and acquired by a flow cytometer using 488 laser and analyzed by Flowjo software (Tree Star, San Carlos, CA).
Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) for in vitro Polarization
In this study, we used LPS and IL-4 for in vitro microglial activation and polarization. On Day 11, LPS (10 μg/ml) or IL-4 (20 ng/ml) were added to the iMG cells and incubated for 24 h at 37°C with 5% CO2. The next day, cell morphology was monitored and captured for each group. The total RNA was isolated using Trizol according to the manufacturer's instruction and subjected to cDNA synthesis using SuperScript™ IV VILO™ Master Mix. With the following primers, for M1 (TNF-α) and M2 (CD206) gene-specific primers, qRT-PCR was performed using PowerUp SYBR Green PCRmix. The relative expression was normalized to the housekeeping control gene β-actin. Fold changes in mRNA levels between stimulated and non-stimulated cells were calculated using the 2−(ΔΔCT) method.
The sequences of primers were as follows:
β-actin F 5′-GATGCAGAAGGAGATCACTGC-3′
β-actin R 5′-ATACTCCTGCTTGCTGATCCA-3′
TNF-α F 5′-CCCAGGCAGTCAGATCATCTTCT-3′
TNF-α R 5′-ATGAGGTACAGGCCCTCTGAT-3′
CD206 F 5′-GTCATTCCGGGTGCTGTTCTCC-3′
CD206 R 5′-TTCGGCATCCTGGTTGCAAG-3′
RNA-seq Data Processing
For RNA-seq, total RNA was extracted from five subjects, including three monocytes, three iMacs, and four iMG cells, by RNeasy Mini kit (Qiagen) following the manufacturer's instruction. The subjects' information, including age, sex, etc., were described in Supplementary Table 1. RNA quality was assessed by Agilent 2100 Bioanalyzer. Total RNA with RIN > 8, 200 ng was subjected to cDNA library construction with RNA Library Prep Kit v2 protocol (Illumina, San Diego, CA). RNA-seq was performed on the polyadenylated fraction of RNA using Illumina NextSeq 500 platform at a depth of more than 20 million reads. FASTQ files from RNA-seq experiments were aligned to reference 1000 Genomes GRCh37 with STAR v2.5.2a (Dobin et al., 2013), and quality-controlled bam files were generated by GATK (v3.7 and v4.beta.2) RNA-seq pipeline (Al-Ars et al., 2020). Read counts from each transcript were finally collected using FeatureCounts v1.5.3 (Liao et al., 2014). The whole genome RNA-seq project PRJNA678841 generated during the current study in the format of fastq files have been deposited at the National Center for Biotechnology Information (NCBI) repository (http://www.ncbi.nlm.nih.gov/bioproject/678841).
Principal Component Analysis (PCA), Multidimensional Scaling (MDS), and Hierarchical Clustering Analysis
To determine the uniformity among the cell types, we conducted PCA using the R package pcaExplorer to perform a simple overview of the counts from their expression profiles (Marini and Binder, 2019). To further determine how the iMG cells resemble the human microglia and hiPSC-induced microglia, we performed MDS analysis with microglia expression datasets from Muffat et al. (2016), Gosselin et al. (2017), and Grubman et al. (2020). Read counts were transformed to fragments per kilobase of transcript per million mapped reads (FPKM) with the edgeR v.3.22.3 (Robinson et al., 2010). Data were then normalized with the removebatcheffect function and rescaled with R scales package. For MDS calculation, we used cmdscale function in R stats and magrittr package (Stefan Milton et al., 2020). The MDS plot was generated with the R ggplot package (Warnes et al., 2020).
For hierarchical clustering analysis, we started with a total of 914 non-overlapped microglial signature genes, including the 881 genes from Gosselin et al.'s study (Gosselin et al., 2017) and 42 microglial genes from McQuade and Blurton-Jones's study (McQuade and Blurton-Jones, 2019) (See Supplementary Table 2). Variance stabilized transformation (VST) data obtained from DESeq2 v1.28.1 were used to process the two datasets (Love et al., 2014). From the 773 overlapped genes between our data and reference data, we compared the expression consistency with the formula for consistency scores = 2xy/(x2 + y2)1/2, where x and y represented the fold change (FC) from monocytes to microglia in the two datasets, respectively. Genes with consistency scores >0.50 were considered as consistent genes (312 genes). We then conducted differential expression analysis DESeq2, as mentioned above. Raw count data were normalized by the library size and log-transformation. Differentially expressed genes were estimated using negative binomial general model statistics. Heatmap analysis was conducted with the consistent genes that overlapped with the differentially expressed genes (201 genes). We used Heatmap.2 function in the R ggplot package (Warnes et al., 2020) for the Heatmap graph.
For Venn analysis, genes with a mean of the normalized counts >10 were assigned to each cell type. The Venn diagram for monocytes, iMacs, and iMG cells was generated according to the website (http://bioinformatics.psb.ugent.be/webtools/Venn/). For the unique genes obtained from the Venn diagram in each specific cell type, we further perform the gene enrichment analysis according to the instruction from Zhou et al.'s paper (Zhou et al., 2019).
Of particular interest, we further explored the gene regulation of some microglia-specific genes reported in Gosselin et al.'s study (Gosselin et al., 2017) or Ryan et al.'s study (Ryan et al., 2017). For this purpose, we used R package pcaExplorer to compare specific gene expression between iMG cells and monocytes from our sample and the sample from (Gosselin et al., 2017).
Statistical Analyses
Statistical analyses for RT-PCR were conducted with Student's t-test or by ANNOVA for comparing more than two groups, and a p-value ≤ 0.05 were considered significant.
Results
Morphology of iMG Cells
In this study, we used fresh or frozen PBMCs to generate human induced microglia-like cells (iMG cells) using a cocktail of IL-34 (100 ng/ml) and GM-CSF (10 ng/ml) as described previously (Ohgidani et al., 2014; Sellgren et al., 2017). Untreated monocytes showed a round morphology (Figure 1A). Upon induction, the microglia-like cells started to grow branches on day 5 (Supplementary Figure 1A). On day 10–14, they displayed a typical morphology of resting ramified microglia with a small cell soma and multiple branching processes (Figure 1B). It was worth noting that the iMG cells could survive up to 36 days when culture media was replaced every 3 or 4 days (Supplementary Figure 1B). For macrophage induction, we treated monocytes with GM-CSF for 5–7 days. On day 5, the phase-contrast image showed the typical “fried-egg” shape of macrophages with enlarged circular cell body (Figure 1C), which was consistent with the previous description (Jin and Kruth, 2016).
Figure 1. iMG cells and iMacs were induced from human peripheral blood monocytes. Human peripheral blood monocytes were isolated from blood and stimulated with IL-34 and GM-CSF for microglia-like cells and GM-CSF only for macrophages. Representative images (A) monocytes after adherence without any treatment, (B) iMG cells on day 14 after monocytes treated with GM-CSF and IL-34, and (C) iMacs on day 5 after monocytes treated with GM-CSF.
Microglia-Specific Markers Expressed by iMG Cells
As a specialized macrophage, microglia share many common markers with macrophages, including Iba1 (Ionized calcium binding adapter molecule 1), transcription factor PU.1, chemokine receptor CX3CR1, and TREM2. On the other hand, transmembrane protein 119 (TMEM119) and purinergic receptor P2RY12 are considered unique markers for microglia but are absent on macrophages (Satoh et al., 2016; Zhu et al., 2017). In this study, we used both myeloid and microglia-specific markers to verify the specificity of iMG cells. As expected, TMEM119 and P2RY12 markers were expressed by iMG cells, as shown in Figure 2A(b,d). These two markers were absent on iMacs, as shown in Figure 2A(a,c). These results were consistent with recent studies from two different groups (Muffat et al., 2016; Douvaras et al., 2017), where hiPSC-derived microglia-like cells also expressed TMEM119 and P2RY12 markers. On the other hand, myeloid lineage markers, such as Iba1, PU.1, CX3CR1, and TREM2, were all expressed by both iMG cells s shown in Figure 2B(b,d,f,h) and iMacs as shown in Figure 2B(a,c,e,g). The results from immunocytochemistry indicated that monocytes-derived iMG cells and iMacs expressed common markers of myeloid lineage but only iMG cells expressed microglia-specific markers, TMEM119 and P2RY12. Therefore, TMEM119 and P2RY12 can be used as markers to distinguish iMG cells from monocyte-derived macrophages.
Figure 2. iMG cells expressed microglia-specific markers (TMEM119 and P2RY12) and other myeloid-specific markers. The expression of surface markers on the iMG cells and iMacs were performed by immunocytochemistry. Peripheral monocytes were incubated with GM-CSF (iMacs) for 5–7 days or cocktail of GM-CSF and IL-34 (iMG cells) for 10–14 days. Representative images of iMacs and iMG cells were shown. (A) iMG cells expressed unique microglial markers TMEM119 (b) and P2RY12 (d), which were absent from iMacs (a,c). (B) Myeloid lineage markers were expressed on both iMacs and iMG cells, such as Iba1 (a,b), PU.1 (c,d), CX3CR1 (e,f), and TREM2 (g,h). Images were captured at 40× with a scale representing 20 μm, where P2RY12, CX3CR1, and TREM2 were stained green, TMEM119, Iba1, and PU.1 were stained red. Cellular nuclei were stained blue with DAPI.
Phagocytosis Assay in iMG Cells
Phagocytosis is one of the most fundamental functions of microglia, and it maintains CNS homeostasis (Fu et al., 2014). Aβ42 are found abundantly in the AD brain, and fibrillary Aβ42 (fAβ42) is closely linked to AD pathogenesis (Crouch et al., 2008). To determine the phagocytic capacity of iMG cells, we used fAβ42 to mimic the physiological condition of AD patients. Considering that Iba1 is a microglia/macrophage-specific calcium-binding protein that participates in membrane ruffling and phagocytosis in activated-microglia, we co-stained the iMG cells with anti-Iba1 antibody to ensure that fAβ42 were physically engulfed inside the cells. As a negative control, no green fluorescent fAβ42 were detected inside iMG cells in Figure 3A(a1–a4), where the iMG cells were not treated with fAβ42. However, green fluorescent fAβ42 were detected inside iMG cells, along with the red fluorescent Iba1 staining in the cytoplasm, as seen in Figure 3A(b1–b4). The confocal images clearly indicated that iMG cells could phagocytose extracellular protein aggregates from their surroundings. Next, we quantitated the percentage of iMG cells with phagocytic activity by flow cytometry. Our result showed that 44.73 ± 4.20% of iMG cells had fluorescent-labeled fAβ42 (Figures 3B,C). Together, we showed that iMG cells indeed could engulf exogenous proteins (fAβ42), implicating potential application for AD study. With the flow cytometry for fluorescent fAβ42 detection, our model has the potential to measure and compare the microglial capacity between different groups, such as between AD patients and controls.
Figure 3. iMG cells phagocytosed fibrillary Aβ42. Confocal imaging and flow cytometry analysis of iMG cells were incubated with human HiLyte™ Fluor 488-labeled Aβ42 proteins. Phagocytosis activity was measured after 4 h incubation. (A) Representative confocal images were shown where controls cells were panel a (a1–a4) and iMG cells with phagocytosed fAβ42 were panel b (b1–b4). Images were captured at 100× with a scale representing 10 μm. fAβ42 and Iba1were stained green and red, respectively. Nuclei were stained blue with DAPI. Merged images of all three channels were shown in a4 and b4. Merged image (b4) indicated Aβ42 proteins were present within the cytoplasm of iMG cells. (B) The Aβ42 protein uptake was quantified by flow cytometry. Representative image of three experiments were shown, where the blue graph was for control and red for iMG cells with engulfed fAβ42. (C) Percentage of cells with engulfed fAβ42 as measured by flow cytometry. The result was from three independent experiments, and the error bars represented ± SD. *p ≤ 0.05 was considered as significant.
M1 and M2 Phenotypic Changes in iMG Cells During Polarization
Upon stimulation, microglia have the capacity to polarize toward a spectrum of activation states ranging from M1 (pro-inflammatory state) to M2 (anti-inflammatory state) (Tang and Le, 2016). Phenotypic changes, including morphology, gene expression, and cytokine release, occur during the polarization (Peng et al., 2017). LPS is a well-known agent for M1 polarization, whereas IL-4 is often used for M2 polarization. With LPS treatment, we noticed that M1 activated-iMG cells had undergone morphological changes showing super-ramified long branches (Figure 4B) compared to iMG cells without activation (Figure 4A). On the other hand, iMG cells activated by IL-4 treatment, exhibited an amoeboid shape, having a round cytoplasm with fewer short branches (Figure 4C). These changes were consistent with previous reports that microglia in M1-state encompassed long processes with ramified shapes, and microglia in M2-state possessed short and thick processes (Lecours et al., 2018; Geirsdottir et al., 2019).
Figure 4. In vitro activation and polarization of iMG cells. The morphology of iMG cells was captured after stimulation with LPS or IL-4 for 24 h. (A) Control iMG cells without any treatment, (B) iMG cells stimulated with 10 μg/ml LPS showing hyper-ramification as indicated with arrows, (C) iMG cells stimulated with 20 ng/ml IL-4 showing ameboid structure as indicated with arrows. Changes in gene expression were analyzed by qRT-PCR showed (D) increased TNF-α levels in LPS-stimulated iMG cells and (E) increased CD206 levels in IL-4-stimulated iMG cells. The experiments were all repeated at least three times. Data were represented as a mean ± SD. Significance between groups were determined by ANOVA test, where ** denoted p < 0.005.
We then tested whether M1 and M2 activated-microglia expressed specific activation markers during the polarization. TNF-α, CD86, CD16/32, and inducible nitric oxide (iNOS) are commonly used as markers for M1 activated-microglia (David and Kroner, 2011; Kobayashi et al., 2013). M2 activated-microglia, however, often shows upregulation of CD163 or CD206 (Cherry et al., 2014). In this study, we used the mRNA expression levels of TNF-α and CD206 as signals for M1 and M2 activation, respectively. As expected, TNF-α mRNA expression in iMG cells was significantly upregulated (12.20 ± 4.00 folds, p < 0.005) when treated with LPS (Figure 4D). IL-4 treatment did not change TNF-α mRNA levels in the iMG cells. Similarly, CD206 mRNA expression was increased (2.41 ± 0.90 folds, p < 0.005) in iMG cells with IL-4 treatment (Figure 4E). LPS did not change CD206 mRNA expression in iMG cells. The upregulation for M1 or M2 specific genes and the typical morphological changes demonstrated that iMG cells were highly plastic that they were able to polarize into M1 or M2 state when responding to exogenous stimulants (LPS or IL-4).
iMG Cells Closely Clustered With Brain-Resident Microglia From RNA-seq Analyses
We showed that iMG cells were similar with primary microglia in morphology, specific markers, and functions. Now we asked whether the iMG cells were similar to human brain-resident microglia or hiPSC-induced microglia from their gene expression profile. For this purpose, we did RNA-seq for monocytes, iMacs, and iMG cells. We also downloaded RNA-seq data from Muffat et al. (2016), Gosselin et al. (2017), and Grubman et al. (2020) for comparison (see Supplementary Table 3 for more details about the datasets). We first conducted PCA analyses for the whole-genome RNA-seq data for the three groups of cells (monocytes, iMacs, and iMG cells) in our induction system. Indeed, we identified three distinctive groups for monocytes, iMacs, and iMG cells from the PCA plot (Supplementary Figure 2), suggesting that the monocytes were successfully differentiated into two other cell types (iMacs or iMG cells) with different induction cytokines. To further determine whether iMG cells resembled the brain-resident microglia or hiPSC-induced microglia, we obtained gene expression data from the above mentioned studies (Muffat et al., 2016; Gosselin et al., 2017; Grubman et al., 2020). We then conducted MDS analyses with normalized data. As shown in PCA plot (Supplementary Figure 2), the MDS plot clearly showed that the three groups of cells were well-separated in our studies (Figure 5A). hiPSC-induced microglia from Muffat et al. and Grubman et al.'s groups, while clustered by each study, there were some distances between the studies. The resident microglia and monocytes from the Gosselin et al.'s study were scattered. This might be due to the different micro-environment where these cells resided. Overall, iMG cells from our study were the closest to the hiPSC-induced microglia and the resident microglia from human brain, but they were separated from monocytes. This confirmed that our iMG cells did have the features of microglia when mRNA expression data were compared.
Figure 5. Multidimensional scaling analysis (MDS) and 2D hierarchical clustering from RNA-seq profile of iMG cells. (A) MDS analysis for multiple RNA-seq datasets, including iMG cells, iMacs, and monocytes from our study, hiPSC-induced microglia (iMGL) and fetal microglia (FML) from Muffat et al.'s study, iMGL from Grubman et al.'s study, and human brain microglia (MGL) and monocytes from Gosselin et al.'s study. (B) 2D hierarchical clustering with 201 top differentially expressed genes from our study and Gosselin et al.'s study, where Ref_hMG and Ref_Mono were human brain-resident microglia and monocytes from Gosselin et al.'s study (Gosselin et al., 2017). iMG cells were induced microglia-like cells from monocytes and iMacs were induced macrophages from monocytes in our study.
We further performed 2D hierarchical clustering analysis of the RNA-seq data from our iMG cells against previously published RNA-seq data from human brain-resident microglia (Ref_hMG) (Gosselin et al., 2017). As mentioned above, a total of 914 microglial signature genes were included in this study (Supplementary Table 2). After matching with our RNA-seq profiles, we obtained an overlap of 805 genes shared between iMG cells and Ref_hMG, among which 312 were consistent between the two groups, as compared with fold changes with monocytes. Based on further differential analysis across iMG cells and monocytes with criteria |FC| (fold changes) ≥1.2, p ≤ 0.001, and average expression ≥5, we generated top 201 microglial differentially expressed genes. With RNA profiles from iMG cells and Ref_hMG, we attempted to evaluate how close the induced-microglia resembled human brain microglia at RNA expression levels. With the top 201 microglial differential genes, we conducted 2D hierarchical clustering analysis. Our data showed that iMG cells clustered closely with Ref_hMG, whereas monocytes from our experiments clustered closely with Ref_Monocytes (Figure 5B). We also noticed that iMG cells and iMacs were clustered closely to each other. These results demonstrated that monocytes treated with IL-34 and GM-CSF underwent significant changes in gene regulation and differentiated toward human microglia-like cells. The expression pattern of surface markers were consistent with a previous report (Ohgidani et al., 2014; Sellgren et al., 2017). As compared to monocytes, the iMG cells and iMacs were clustered more closely, reflecting that iMG cells were specialized macrophages in the brain.
Venn Analysis and Differential Analysis Across iMG Cells, iMacs, and Monocytes
To comprehend the differences between monocytes, iMacs, and iMG cells, we performed Venn analysis across the three groups. Unique genes to each cell type and shared genes across three or two cell types were shown in the Venn diagram (Figure 6A), where 350 genes were uniquely expressed in iMG cells. We also found that 409 genes and 859 genes were specifically expressed in iMacs and monocytes, respectively. 15,801 genes were shared across all three groups.
Figure 6. Venn and differential expression analysis and enrichment analysis. (A) Venn diagram of genes expressed by iMG cells, iMacs, and monocytes (defined by the mean of the normalized counts >10). Blue: Monocytes; Red: iMacs; Green: iMG cells. Enrichment analysis with the unique genes expressed in iMG cells (B), monocytes (C), and iMacs (D), ranked by pathway p-value. (E) Volcano plot depicting differentially expressed genes in iMG cells against monocytes. The dashed line indicated p < 1 × 10−16 (horizontal line) and |log2FC| >2 (Vertical line). Two well-known AD risk genes, APOE and TREM2, were both upregulated in iMG cells as labeled in the Volcano plot. (F) Relative expression of significant genes (C1QA, C1QB, C1QC, C2, C3, C9, CD59, and SEPRING1) from the complement system in iMG cells vs. monocytes, and iMacs and. p values were obtained from differential expression analysis DESeq2 v1.28.1, where **** denotes p < 0.001.
Gene enrichment analysis was performed with unique genes in each group. Our results showed that the unique genes from iMG cells were enriched in the pathways involving in the regulation of synaptic vesicle fusion and ion transport, and TNFs bind their physiological receptors, which aligned with microglial function in regard to their interaction with neurons and response to the environment (Figure 6B). The unique genes from monocytes were enriched in innate immune-related pathways, such as complement activation, phagocytosis, and cytokine-cytokine receptor interaction (Figure 6C). Surprisingly, 409 unique genes from iMacs were not enriched in many pathways, but only associated with the transport of small molecules pathway (Figure 6D).
The differences in gene expression between iMG cells and monocytes were also shown in the volcano plot, where significant genes were those that passed the FDR (p < 1 × 10−16) and |log2FC| > 2, as indicated in red (Figure 6E). Interestingly, APOE and TREM2, two of the most important risk genes for AD, were both significantly upregulated (Figure 6E), with p = 7.09 × 10−31 and 2.35 × 10−40, respectively. Differentially expressed genes (DEGs) between iMG cells and monocytes, or Ref_hMG and monocytes, were shown in Supplementary Tables 4, 5.
Genes in the complement system are highly expressed in microglia and play an important role in synaptic pruning that may underlie the mechanisms of brain diseases, such as schizophrenia and AD. Ryan et al.'s group recently reported a microglia-like model generated from human monocytes in the presence of polarizing cytokines (IL-34, GM-CSF, M-CSF, NGF-β, and CCL2). We wanted to compare some of the microglia-specific genes from those of previous studies. Our data showed that genes involved in the synaptic remodeling, such as complement factors C1QA, C1QB, C1QC, C2, C3, C9, C59, and SERPING1, were highly expressed in our iMG cells, consistent with Gosselin et al.'s report and Ryan et al.'s report (Figure 6F) (Gosselin et al., 2017; Ryan et al., 2017). Among those, C1QA, C1QB, C1QC, and C2 were only upregulated in iMG cells but not iMacs. We also found that more than half of genes (19) in iMG cells overlapped with 37 genes from MDMi (human monocyte-derived microglia-like cells) from Ryan et al.'s group (Ryan et al., 2017) (Supplementary Table 6). Interestingly, four genes mentioned as being microglia-specific in mice in Ryan et al.'s report, TGFβR1, PROS1, P2RX7, and C1QB, were consistently upregulated in our iMG cells (Supplementary Figure 3). As mentioned above, two of the most important risk genes for AD, APOE and TREM2, were also significantly upregulated (Supplementary Figure 3).
Discussion
Microglia are essential for the maintenance of normal brain functions. Massive evidence shows that microglial dysfunction contributes to neurodevelopmental and neurodegenerative processes. Over the past 100 years of research on microglia, advances have been made in understanding microglial plasticity, transcriptional diversity, function-related signaling pathways, and stem-cell-based microglial induction systems (Umpierre and Wu, 2020). Technical developments also have advanced scientists to explore human microglia-like cellular models from human blood monocytes (Ohgidani et al., 2014; Ryan et al., 2017; Sellgren et al., 2019), or hiPSC (Muffat et al., 2016; Abud et al., 2017; Douvaras et al., 2017; Haenseler et al., 2017; Pandya et al., 2017; Svoboda et al., 2019). In vitro differentiation from hiPSC to microglia is challenging as the reprogramming process is time-consuming (Liang and Zhang, 2013). It is essential to establish a more efficient microglial model closely resembling the human microglial phenotypic and functional features. It is also crucial to build the model from a relatively easy and non-invasive sample collection process. With modification from Ohgidani et al.'s report (Ohgidani et al., 2014), we generated and validated the microglia-like model from human peripheral monocytes. Consistent with the previous studies (Ohgidani et al., 2014; Sellgren et al., 2017), our data showed that iMG cells had the typical ramified morphology as resting microglia in the brain and expressed microglia-specific surface markers, such as TMEM119 and P2RY12. On the contrary, iMacs with GM-CSF treatment showed a typical “fried-egg” morphology without TMEM119 and P2RY12 expression. Another group has used a complex cocktail (IL-34, GM-CSF, M-CSF, NGF-β, and CCL2) for microglia induction (MDMi) (Ryan et al., 2017). While comparing microglial genes, iMG cells and MDMi shared more than half of those microglia-specific genes and 4 specific genes were significantly upregulated in iMG cells, MDMi, and human microglia. These results demonstrated that microglia induced from monocytes with different inductions shared some common microglia features, though, at this stage it is not clear what the difference is between the two inductions. We know that differentiation from hiPSC to microglia-like cells usually takes around 40 days (Muffat et al., 2016; Grubman et al., 2020). In this study, our induction from monocyte to microglia took <2 weeks. In addition, blood collection from patients is widely accepted for clinical research. Thus, there is no doubt that our model will provide a more feasible and efficient approach to study microglial functions in clinical trials.
As mentioned above, microglia are professional phagocytic cells of the brain that remove apoptotic or necrotic cells (Green et al., 2016) and eliminate unfolded proteins, such as Aβ (Daniel Lee and Landreth, 2010). To evaluate the phagocytic capacity of iMG cells, we chose fibrillary Aβ as an exogenous proteins to best mimic the physiological condition for phagocytosis in most AD patients. Nonetheless, it should be noted that different Aβ preparation methods can affect Aβ aggregation kinetics and structures (Yates and Legleiter, 2014; Nirmalraj et al., 2020). It is still not clear under what circumstances various types of Aβ are formed physiologically inside the brain (Chen et al., 2017). Microglia are believed to initially clear Aβ deposits at a preclinical stage of AD. With the progress of AD, they somehow lose their ability to clear Aβ as they are activated to a pro-inflammatory state where cytokines, chemokines, and reactive oxygen species (ROS) are released (Hale et al., 2008). To our knowledge, we are the first to show the phagocytic capacity of monocyte-derived microglia engulfing exogenous Aβ proteins. This is crucial as utilizing this microglial model will help us to understand the underlying mechanism by which microglia remove Aβ proteins under physiological condition and how they lose the ability of Aβ clearance during pathological condition. Our cellular model provides a feasible tool to measure microglial function of phagocytosis. With different compounds or with genetic engineering approaches, iMG cells can be used to evaluate their effects on removing internal or external proteins, including Aβ proteins in AD clinical trials. As a result, this model would offer a prompt monitor system for microglial functions at the cellular level, which facilitates new drug discovery or disease-modified treatment.
Microglia are highly plastic; they can be activated into the pro-inflammatory (M1) state or into the anti-inflammatory (M2) state in response to various stimuli. Such microglia remodeling is critical for microglial phagocytosis, synaptic modifications, and neuroinflammatory responses and normal brain functioning. As an example, pro-inflammatory M1 microglia are usually dominant in AD and schizophrenia pathology. With morphological changes and specific polarization markers upregulation, our data demonstrated that iMG cells could be transformed into M1 state with LPS stimulation or M2 state with IL-4 stimulation. Thus, iMG cells provide a human microglial cellular model that can be used to study the mechanisms underlying the microglial polarization (Svoboda et al., 2019). Future studies with the iMG cells model for comprehensive comparison between patients and healthy controls have the potential to identify the microglia-relevant mechanisms associated with AD, schizophrenia, and other diseases.
In this study, we also conducted a comprehensive gene expression profiling of iMG cells and compared it with expression data from brain-derived microglia. From cluster analysis, our data indicated that iMG cells were closely clustered with brain-resident microglia or hiPSC-induced microglia. Interestingly, the unique genes from iMG cells were enriched in pathways involved in the regulation of synaptic vesicle and ion transport, and synapse organization. Multiple genes in the complement system, C1QA, C1QB, C1QC, C2, C3, C9, CD59, and SEPRING1, were also highly up-regulated. This confirms that microglia-specific genes had been remodeled during the differentiation from monocytes to iMG cells and possess specialized brain functions (e.g., synaptic remodeling). Excessive evidence has shown that microglial dysfunction contributes to synaptic abnormalities and cognitive decline, which are often seen in neurodevelopmental and neurodegenerative disorders. Indeed, Sellgren et al. recently applied the model to investigate schizophrenia. They found that iMG cells derived from patients with schizophrenia eliminated more synapses than those from healthy controls (Sellgren et al., 2019). Sellgren et al. demonstrated that iMG cells would be a good model to study schizophrenia. Several other groups have successfully applied this iMG model or similar model to study the microglial functions in different disorders, including bipolar disorder, fibromyalgia, and schizophrenia (Ohgidani et al., 2017a,b; Ormel et al., 2020). In this study, we demonstrated that the iMG cellular model would also be a good model for studying cellular changes and gene regulation during the development of AD, mainly the phagocytosis of fAβ42 and the gene expression profiles. The two genes TREM2 and APOE are widely accepted as high risk genes for non-familial AD (Shi and Holtzman, 2018; Wolfe et al., 2018). Our RNA expression data indicated the two genes were significantly upregulated in iMG cells as compared to monocytes.
Our iMG cellular model has some limitations. In spite of the promising results, the model might still not fully recapitulate the microglial phenotypes from the human brain. We understand this limitation and acknowledge that this is a challenge to any surrogate cellular model in the field. Nevertheless, our study showed that the iMG cellular model had extensive human microglial features at phenotypic, functional, and gene expression levels that warrant its application in schizophrenia and AD studies. Other limitations include subject to subject variations, cell viability during induction, and high sensitivity to culture microenvironments. We acknowledge that these limitations may create some variations across batches of experiments. However, the results from our in vitro iMG cells were highly compatible and consistent with studies from other groups. Therefore, our model derived from patients with psychiatric diseases or neurodegenerative diseases would be a good alternative for the study of gene regulations and characterization of disease-related microglial phenotypes, provided that each experiment has age- and sex-matched controls. Furthermore, our study indicated that this microglial model can be extended to study the molecular mechanisms underlying synapse remodeling during brain development and neurodegeneration. Finally, this iMG model can also be used for high throughput studies for the discovery of new biomarkers or new drugs.
Overall, our findings confirmed that iMG cellular model exhibited cellular and molecular characteristics that were similar to human brain-resident microglia. Functional studies of such patient-specific microglial model will facilitate a better understanding of neuroinflammatory processes in psychiatric and neurodegenerative diseases. We were the first to validate the iMG cellular model with comprehensive RNA profile, although we were not the first study to generate this human monocyte-derived microglial model. Additionally, this study investigated the plasticity of the induced-microglia with two distinct polarization states. Therefore, the present study provides valuable data regarding the application of this model to study the pathogenesis of psychiatric and neurological disorders, including schizophrenia and AD, two diseases with cognitive deficits closely linking to microglial activation, synaptic pruning, and neuroinflammation. Further studies using this cellular model derived from patient samples are required to confirm our findings.
Data Availability Statement
The whole genome RNA-seq project PRJNA678841 generated during the current study in the format of fastq files have been deposited at the National Center for Biotechnology Information (NCBI) repository (http://www.ncbi.nlm.nih.gov/bioproject/678841).
Author Contributions
JC conceived the idea. JC and AB conducted the experiments and analyzed the data. KD performed the qRT-PCR. TM collected and processed the RNA-seq raw data in this study. JC, YL, and XW conducted the RNA-seq and bioinformatics analysis. JC, AB, and XC contributed to the writing of the manuscript. All authors reviewed the manuscript and agreed on its content.
Funding
This work was supported by grant R01 MH101054 to XC from the National Institute of Mental Health, pilot grant from grant U54GM104944 to JC from NIGMS-CTRIN, NIH grant P20 GM121325 to JC from NIGMS-COBRE, NIH P20 GM121325 03S1, NIH P20 GM121325-02S1, and GM121325_03S2 to JC from NIGMS-COBRE.
Conflict of Interest
XC is employed by 410 AI, LLC.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
RNA-seq data analyses were performed with computational resources from the Nevada Institute of Personalized Medicine (NIPM) and the National Supercomputing Institute at University of Nevada Las Vegas (UNLV). The authors would like to thank the Genome Acquisition and Analysis (GAA) Core from the NIH Center of Biomedical Research Excellence (COBRE) at NIPM. Confocal imaging was performed with the assistance of Sophie Choe at the Confocal and Biological Imaging Core at UNLV. The Nikon A1Rmp microscope system was funded by NSF MRI 1726925. We want to acknowledge Lorena Samentar and Nora Caberoy (School of Life Sciences, UNLV) for helping us with Aβ protein preparation. We also thank Bridgette Choi for helping with manuscript formatting and proofreading.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncel.2021.629279/full#supplementary-material
Supplementary Figure 1. The morphology of iMG cells on days 5 and 36. (A) Branched cells were observed on day 5 after the treatment with IL-34 and GM-CSF. (B) The iMG cells survived up to 36 days.
Supplementary Figure 2. Principal component analysis (PCA) from RNA-seq profile. PCA plot of monocytes, iMacs, and iMG cells from whole-genome RNA-seq profile.
Supplementary Figure 3. Upregulation of microglia-specific genes and AD risk genes in iMG Cells. Four genes mentioned as being mice microglia-specific in Ryan et al.'s report, TGFβR1, PROS1, P2RX7, and C1QB, were consistently upregulated in our iMG cells. Two AD-risk genes, APOE and TREM2, were significantly upregulated in iMG cells. P-values were obtained from differential expression analysis DESeq2 v1.28.1., where **** denotes p < 0.001, *** denotes p < 0.005.
Supplementary Table 1. Sample information for RNA-seq profiling.
Supplementary Table 2. Microglia signature genes in this study.
Supplementary Table 3. Datasets for MDS analysis in this study.
Supplementary Table 4. DEGs between iMG cells and monocytes.
Supplementary Table 5. DEGs between Ref_hMG and monocytes.
Supplementary Table 6. Overlapping genes between iMG cells and Ryan et al.'s data or Ref_hMG and Ryan et al.'s data.
References
Abud, E. M., Ramirez, R. N., Martinez, E. S., Healy, L. M., Nguyen, C. H. H., Newman, S. A., et al. (2017). iPSC-derived human microglia-like cells to study neurological diseases. Neuron 94, 278–293.e9. doi: 10.1016/j.neuron.2017.03.042
Akiyama, H., Barger, S., Barnum, S., Bradt, B., Bauer, J., Cole, G. M., et al. (2000). Inflammation and Alzheimer's disease. Neurobiol. Aging 21, 383–421. doi: 10.1016/S0197-4580(00)00124-X
Al-Ars, Z., Wang, S., and Mushtaq, H. (2020). SparkRA: enabling big data scalability for the GATK RNA-seq pipeline with apache spark. Genes 11:53. doi: 10.3390/genes11010053
Chen, G., Xu, T., Yan, Y., Zhou, Y., Jiang, Y., Melcher, K., et al. (2017). Amyloid beta: structure, biology and structure-based therapeutic development. Acta Pharmacol. Sin. 38, 1205–1235. doi: 10.1038/aps.2017.28
Cherry, J. D., Olschowka, J. A., and O'Banion, M. K. (2014). Neuroinflammation and M2 microglia: the good, the bad, and the inflamed. J. Neuroinflammation 11:98. doi: 10.1186/1742-2094-11-98
Crouch, P. J., Harding, S. M. E., White, A. R., Camakaris, J., Bush, A. I., and Masters, C. L. (2008). Mechanisms of A beta mediated neurodegeneration in Alzheimer's disease. Int. J. Biochem. Cell Biol. 40, 181–198. doi: 10.1016/j.biocel.2007.07.013
Daniel Lee, C. Y., and Landreth, G. E. (2010). The role of microglia in amyloid clearance from the AD brain. J. Neural Transm. 117, 949–960. doi: 10.1007/s00702-010-0433-4
David, S., and Kroner, A. (2011). Repertoire of microglial and macrophage responses after spinal cord injury. Nat. Rev. Neurosci. 12, 388–399. doi: 10.1038/nrn3053
Deczkowska, A., Keren-Shaul, H., Weiner, A., Colonna, M., Schwartz, M., and Amit, I. (2018). Disease-associated microglia: a universal immune sensor of neurodegeneration. Cell 173, 1073–1081. doi: 10.1016/j.cell.2018.05.003
Dobin, A., Davis, C. A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., et al. (2013). STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21. doi: 10.1093/bioinformatics/bts635
Douvaras, P., Sun, B., Wang, M., Kruglikov, I., Lallos, G., Zimmer, M., et al. (2017). Directed differentiation of human pluripotent stem cells to microglia. Stem Cell Rep. 8, 1516–1524. doi: 10.1016/j.stemcr.2017.04.023
Elmore, M. R. P., Najafi, A. R., Koike, M. A., Dagher, N. N., Spangenberg, E. E., Rice, R. A., et al. (2014). Colony-stimulating factor 1 receptor signaling is necessary for microglia viability, unmasking a microglia progenitor cell in the adult brain. Neuron 82, 380–397. doi: 10.1016/j.neuron.2014.02.040
Franco, R., and Cedazo-Minguez, A. (2014). Successful therapies for Alzheimer's disease: why so many in animal models and none in humans? Front. Pharmacol. 5:146. doi: 10.3389/fphar.2014.00146
Friedman, B. A., Srinivasan, K., Ayalon, G., Meilandt, W. J., Lin, H., Huntley, M. A., et al. (2018). Diverse brain myeloid expression profiles reveal distinct microglial activation states and aspects of Alzheimer's disease not evident in mouse models. Cell Rep. 22, 832–847. doi: 10.1016/j.celrep.2017.12.066
Fu, R., Shen, Q., Xu, P., Luo, J., and Tang, Y. (2014). Phagocytosis of microglia in the central nervous system diseases. Mol. Neurobiol. 49, 1422–34. doi: 10.1007/s12035-013-8620-6
Geirsdottir, L., David, E., Keren-Shaul, H., Weiner, A., Bohlen, S. C., Neuber, J., et al. (2019). Cross-species single-cell analysis reveals divergence of the primate microglia program. Cell 179, 1609–1622.e16. doi: 10.1016/j.cell.2019.11.010
Ginhoux, F., Lim, S., Hoeffel, G., Low, D., and Huber, T. (2013). Origin and differentiation of microglia. Front. Cell. Neurosci. 7:45. doi: 10.3389/fncel.2013.00045
Ginhoux, F., and Prinz, M. (2015). Origin of microglia: current concepts and past controversies. Cold Spring Harb. Perspect. Biol. 7:a020537. doi: 10.1101/cshperspect.a020537
Gosselin, D., Skola, D., Coufal, N. G., Holtman, I. R., Schlachetzki, J. C. M., Sajti, E., et al. (2017). An environment-dependent transcriptional network specifies human microglia identity. Science 356:eaal3222. doi: 10.1126/science.aal3222
Green, D. R., Oguin, T. H., and Martinez, J. (2016). The clearance of dying cells: table for two. Cell Death Differ. 23, 915–926. doi: 10.1038/cdd.2015.172
Grubman, A., Vandekolk, T. H., Schröder, J., Sun, G., Hatwell-Humble, J., Chan, J., et al. (2020). A CX3CR1 reporter hESC line facilitates integrative analysis of in-vitro-derived microglia and improved microglia identity upon neuron-glia co-culture. Stem Cell Rep. 14, 1018–1032. doi: 10.1016/j.stemcr.2020.04.007
Haenseler, W., Sansom, S. N., Buchrieser, J., Newey, S. E., Moore, C. S., Nicholls, F. J., et al. (2017). A highly efficient human pluripotent stem cell microglia model displays a neuronal-co-culture-specific expression profile and inflammatory response. Stem Cell Rep. 8, 1727–1742. doi: 10.1016/j.stemcr.2017.05.017
Hale, C., Véniant, M., Wang, Z., Chen, M., McCormick, J., Cupples, R., et al. (2008). Structural Characterization and pharmacodynamic effects of an orally active 11β-hydroxysteroid dehydrogenase type 1 inhibitor. Chem. Biol. Drug Des. 71, 36–44. doi: 10.1111/j.1747-0285.2007.00603.x
Heneka, M. T., Carson, M. J., El Khoury, J., Landreth, G. E., Brosseron, F., Feinstein, D. L., et al. (2015). Neuroinflammation in Alzheimer's disease. Lancet Neurol. 14, 388–405. doi: 10.1016/S1474-4422(15)70016-5
Janabi, N., Peudenier, S., Héron, B., Ng, K. H., and Tardieu, M. (1995). Establishment of human microglial cell lines after transfection of primary cultures of embryonic microglial cells with the SV40 large T antigen. Neurosci. Lett. 195, 105–108. doi: 10.1016/0304-3940(94)11792-h
Jin, X., and Kruth, H. S. (2016). Culture of macrophage colony-stimulating factor differentiated human monocyte-derived macrophages. J. Vis. Exp. 112:e54244. doi: 10.3791/54244
Jungbauer, L. M., Yu, C., Laxton, K. J., and LaDu, M. J. (2009). Preparation of fluorescently-labeled amyloid-beta peptide assemblies: the effect of fluorophore conjugation on structure and function. J. Mol. Recognit. 22, 403–413. doi: 10.1002/jmr.948
Kobayashi, K., Imagama, S., Ohgomori, T., Hirano, K., Uchimura, K., Sakamoto, K., et al. (2013). Minocycline selectively inhibits M1 polarization of microglia. Cell Death Dis. 4:e525. doi: 10.1038/cddis.2013.54
Lecours, C., Bordeleau, M., Cantin, L., Parent, M., Paolo, T. D., and Tremblay, M. È. (2018). Microglial implication in Parkinson's disease: loss of beneficial physiological roles or gain of inflammatory functions? Front. Cell. Neurosci. 12:282. doi: 10.3389/fncel.2018.00282
Liang, G., and Zhang, Y. (2013). Genetic and epigenetic variations in iPSCs: potential causes and implications for application. Cell Stem Cell 13, 149–159. doi: 10.1016/j.stem.2013.07.001
Liao, Y., Smyth, G. K., and Shi, W. (2014). featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930. doi: 10.1093/bioinformatics/btt656
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8
Ma, Y., Wang, J., Wang, Y., and Yang, G. Y. (2017). The biphasic function of microglia in ischemic stroke. Prog. Neurobiol. 157, 247–272. doi: 10.1016/j.pneurobio.2016.01.005
Marini, F., and Binder, H. (2019). pcaExplorer: an R/Bioconductor package for interacting with RNA-seq principal components. BMC Bioinformatics 20:331. doi: 10.1186/s12859-019-2879-1
McQuade, A., and Blurton-Jones, M. (2019). Microglia in Alzheimer's disease: exploring how genetics and phenotype influence risk. J. Mol. Biol. 431, 1805–1817. doi: 10.1016/j.jmb.2019.01.045
Monji, A., Kato, T. A., Mizoguchi, Y., Horikawa, H., Seki, Y., Kasai, M., et al. (2013). Neuroinflammation in schizophrenia especially focused on the role of microglia. Prog. Neuropsychopharmacol. Biol. Psychiatry 42, 115–121. doi: 10.1016/j.pnpbp.2011.12.002
Muffat, J., Li, Y., Yuan, B., Mitalipova, M., Omer, A., Corcoran, S., et al. (2016). Efficient derivation of microglia-like cells from human pluripotent stem cells. Nat. Med. 22, 1358–1367. doi: 10.1038/nm.4189
Nagai, A., Mishima, S., Ishida, Y., Ishikura, H., Harada, T., Kobayashi, S., et al. (2005). Immortalized human microglial cell line: phenotypic expression. J. Neurosci. Res. 81, 342–348. doi: 10.1002/jnr.20478
Nirmalraj, P. N., List, J., Battacharya, S., Howe, G., Xu, L., Thompson, D., et al. (2020). Complete aggregation pathway of amyloid β (1-40) and (1-42) resolved on an atomically clean interface. Sci. Adv. 6:eaaz6014. doi: 10.1126/sciadv.aaz6014
Ohgidani, M., Kato, T. A., Haraguchi, Y., Matsushima, T., Mizoguchi, Y., Murakawa-Hirachi, T., et al. (2017a). Microglial CD206 gene has potential as a state marker of bipolar disorder. Front. Immunol. 7:676. doi: 10.3389/fimmu.2016.00676
Ohgidani, M., Kato, T. A., Hosoi, M., Tsuda, M., Hayakawa, K., Hayaki, C., et al. (2017b). Fibromyalgia and microglial TNF-α: translational research using human blood induced microglia-like cells. Sci. Rep. 7:11882. doi: 10.1038/s41598-017-11506-4
Ohgidani, M., Kato, T. A., Setoyama, D., Sagata, N., Hashimoto, R., Shigenobu, K., et al. (2014). Direct induction of ramified microglia-like cells from human monocytes: dynamic microglial dysfunction in Nasu-Hakola disease. Sci. Rep. 4:4957. doi: 10.1038/srep04957
Ormel, P. R., Böttcher, C., Gigase, F. A. J., Missall, R. D., van Zuiden, W., Fernández Zapata, M. C., et al. (2020). A characterization of the molecular phenotype and inflammatory response of schizophrenia patient-derived microglia-like cells. Brain. Behav. Immun. 90, 196–207. doi: 10.1016/j.bbi.2020.08.012
Pandya, H., Shen, M. J., Ichikawa, D. M., Sedlock, A. B., Choi, Y., Johnson, K. R., et al. (2017). Differentiation of human and murine induced pluripotent stem cells to microglia-like cells. Nat. Neurosci. 20, 753–759. doi: 10.1038/nn.4534
Peng, H., Geil Nickell, C. R., Chen, K. Y., McClain, J. A., and Nixon, K. (2017). Increased expression of M1 and M2 phenotypic markers in isolated microglia after four-day binge alcohol exposure in male rats. Alcohol 62, 29–40. doi: 10.1016/j.alcohol.2017.02.175
Pérez-Neri, I., Ramírez-Bermúdez, J., Montes, S., and Ríos, C. (2006). Possible mechanisms of neurodegeneration in schizophrenia. Neurochem. Res. 31, 1279–1294. doi: 10.1007/s11064-006-9162-3
Río-Hortega, P. (1932). Estructuray Sistematización de los Gliomasy Paragliomas. Arch. Esp. (De Oncología), 411–677.
Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi: 10.1093/bioinformatics/btp616
Ryan, K. J., White, C. C., Patel, K., Xu, J., Olah, M., Replogle, J. M., et al. (2017). A human microglia-like cellular model for assessing the effects of neurodegenerative disease gene variants. Sci. Transl. Med. 9:eaai7635. doi: 10.1126/scitranslmed.aai7635
Satoh, J., Kino, Y., Asahina, N., Takitani, M., Miyoshi, J., Ishida, T., et al. (2016). TMEM119 marks a subset of microglia in the human brain. Neuropathology 36, 39–49. doi: 10.1111/neup.12235
Sellgren, C. M., Gracias, J., Watmuff, B., Biag, J. D., Thanos, J. M., Whittredge, P. B., et al. (2019). Increased synapse elimination by microglia in schizophrenia patient-derived models of synaptic pruning. Nat. Neurosci. 22, 374–385. doi: 10.1038/s41593-018-0334-7
Sellgren, C. M., Sheridan, S. D., Gracias, J., Xuan, D., Fu, T., and Perlis, R. H. (2017). Patient-specific models of microglia-mediated engulfment of synapses and neural progenitors. Mol. Psychiatry 22, 170–177. doi: 10.1038/mp.2016.220
Shi, Y., and Holtzman, D. M. (2018). Interplay between innate immunity and Alzheimer's disease: APOE and TREM2 in the spotlight. Nat. Rev. Immunol. 18, 759–772. doi: 10.1038/s41577-018-0051-1
Stansley, B., Post, J., and Hensley, K. (2012). A comparative review of cell culture systems for the study of microglial biology in Alzheimer's disease. J. Neuroinflammation 9:115. doi: 10.1186/1742-2094-9-115
Stefan Milton, B., Wickham, H., Henry, L., and RStudio (2020). magrittr: A Forward-Pipe Operator for R. Available online at: https://CRAN.R-project.org/package=magrittr (accessed January 10, 2021).
Stine, W. B., Dahlgren, K. N., Krafft, G. A., and LaDu, M. J. (2003). In vitro characterization of conditions for amyloid-beta peptide oligomerization and fibrillogenesis. J. Biol. Chem. 278, 11612–11622. doi: 10.1074/jbc.M210207200
Svoboda, D. S., Barrasa, M. I., Shu, J., Rietjens, R., Zhang, S., Mitalipova, M., et al. (2019). Human iPSC-derived microglia assume a primary microglia-like state after transplantation into the neonatal mouse brain. Proc. Natl. Acad. Sci. U.S.A. 116, 25293–25303. doi: 10.1073/pnas.1913541116
Takata, K., Kozaki, T., Lee, C. Z. W., Thion, M. S., Otsuka, M., Lim, S., et al. (2017). Induced-pluripotent-stem-cell-derived primitive macrophages provide a platform for modeling tissue-resident macrophage differentiation and function. Immunity 47, 183–198.e6. doi: 10.1016/j.immuni.2017.06.017
Tang, Y., and Le, W. (2016). Differential roles of M1 and M2 microglia in neurodegenerative diseases. Mol. Neurobiol. 53, 1181–1194. doi: 10.1007/s12035-014-9070-5
Umpierre, A. D., and Wu, L. J. (2020). Microglia research in the 100th year since its discovery. Neurosci. Bull. 36, 303–306. doi: 10.1007/s12264-020-00477-8
Wang, W. Y., Tan, M. S., Yu, J. T., and Tan, L. (2015). Role of pro-inflammatory cytokines released from microglia in Alzheimer's disease. Ann. Transl. Med. 3:136. doi: 10.3978/j.issn.2305-5839.2015.03.49
Warnes, G. R., Bolker, B., Bonebakker, L., Gentleman, R., Huber, W., Liaw, A., et al. (2020). gplots: Various R Programming Tools for Plotting Data. Available online at: https://CRAN.R-project.org/package=gplots (accessed March 6, 2020).
Wolfe, C. M., Fitz, N. F., Nam, K. N., Lefterov, I., and Koldamova, R. (2018). The role of APOE and TREM2 in Alzheimer's disease-current understanding and perspectives. Int. J. Mol. Sci. 20:81. doi: 10.3390/ijms20010081
Yates, E. A., and Legleiter, J. (2014). Preparation protocols of aβ(1-40) promote the formation of polymorphic aggregates and altered interactions with lipid bilayers. Biochemistry 53, 7038–7050. doi: 10.1021/bi500792f
Zhou, Y., Song, W. M., Andhey, P. S., Swain, A., Levy, T., Miller, K. R., et al. (2020). Human and mouse single-nucleus transcriptomics reveal TREM2-dependent and TREM2-independent cellular responses in Alzheimer's disease. Nat. Med. 26, 131–142. doi: 10.1038/s41591-019-0695-9
Zhou, Y., Zhou, B., Pache, L., Chang, M., Khodabakhshi, A. H., Tanaseichuk, O., et al. (2019). Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat. Commun. 10:1523. doi: 10.1038/s41467-019-09234-6
Keywords: induced microglia-like cells (iMG cells), RNA-seq profile, peripheral blood mononuclear cells, phagocytosis, schizophrenia, Alzheimer's disease, synaptic pruning
Citation: Banerjee A, Lu Y, Do K, Mize T, Wu X, Chen X and Chen J (2021) Validation of Induced Microglia-Like Cells (iMG Cells) for Future Studies of Brain Diseases. Front. Cell. Neurosci. 15:629279. doi: 10.3389/fncel.2021.629279
Received: 17 November 2020; Accepted: 03 March 2021;
Published: 09 April 2021.
Edited by:
Adam Bachstetter, University of Kentucky, United StatesReviewed by:
Wayne W. Poon, University of California, Irvine, United StatesTakahiro A. Kato, Kyushu University, Japan
Copyright © 2021 Banerjee, Lu, Do, Mize, Wu, Chen and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jingchun Chen, jingchun.chen@unlv.edu