Maternal SARS-CoV-2 impacts fetal placental macrophage programs and placenta-derived microglial models of neurodevelopment

Hofbauer cells are a heterogeneous population with subclusters demonstrating both M1- and M2-like transcriptional signatures

Placental chorionic villous tissues were collected from N=12 birthing individuals: N=4 from individuals who had a positive SARS-CoV-2 nasopharyngeal PCR test during pregnancy, and N=8 from birthing individuals with a negative PCR at delivery and no history of a positive SARS-CoV-2 test during pregnancy. In SARS-CoV-2 positive maternal cases, infections occurred remote from delivery (median 19.5 weeks). No participants had received a COVID-19 vaccine prior to or during pregnancy and no placental samples were infected with SARS-CoV-2 at delivery (defined as having detectable SARS-CoV-2 viral load in a validated assay sensitive to 40 copies/mL) (24, 25). Thus, SARS-CoV-2 positive cases were defined by maternal infection in pregnancy, not placental or fetal infection. Additional participant characteristics are provided in Table 1.

To assess the HBC transcriptome, we first used a previously-described protocol to obtain primarily HBCs from placental villi; in this protocol, a Percoll-based gradient and negative bead-based selection steps are used to isolate putative HBCs from other cell types present in the chorionic villi (including trophoblasts, fibroblasts) (26). Single-cell RNA sequencing was then performed on all cell suspensions (10x Genomics). After quality control filtering to remove putative doublets and cells with less than 300 identified genes, we obtained a dataset comprised of a total of 70,817 cells. We then performed sample integration and graph-based clustering to identify broad cell types (Figure S1A). Based on analyses of marker gene expression (Figure S1B), we found that the majority of cells in our dataset had marker gene expression consistent with monocytes/macrophages, and that other cell types were represented to a lesser extent, including some fibroblasts, vascular endothelial cells (VECs), extravillous trophoblasts (EVTs), leukocytes (NK cells, CD8+ T cells, B cells), neutrophils and red blood cells (Figure S1C). From this dataset, we excluded all cell types that were not identifiable as macrophages/monocytes. After additional quality control filtering for nUMIs, gene counts, and percent mitochondrial reads (see Methods), this resulted in a dataset containing 31,719 high-quality placental macrophages/monocytes. All subsequent analyses were performed with this final dataset.

After re-processing selected cells for quality control as described, we identified 10 total subclusters of macrophages/monocytes (Figure 1A), with representation of each subcluster across donors from both SARS-CoV-2 positive cases and controls (Figure S2A). To distinguish HBCs, which are placental macrophages of fetal origin, from macrophages or monocytes of maternal origin, we used cells isolated from male placentas (N=10). Male fetal origin was confirmed by high expression of DDX3Y and low expression of XIST in 8 subclusters; these were labeled HBC 0-7 (Figure 1B). The macrophage cluster with high expression of XIST (consistent with maternal origin) was annotated as placenta-associated maternal macrophages and monocytes (PAMMs, Figure 1B) (27). A small cluster of monocytes – identified as such by high expression of monocyte marker genes S100A8, S100A9, and TIMP1 – demonstrated equal expression levels of both DDX3Y and XIST, suggesting that both fetal and maternal cells were present in this monocyte cluster. To further support HBC cluster annotation, we next compared the overall gene expression profiles of each cluster to a previously published single-cell dataset derived from human first-trimester placenta and decidua (28). In this analysis, all putative HBC subclusters showed highest correlation with HBC expression profiles from this dataset, whereas the monocyte and PAMM clusters had higher correlation with decidual macrophages than HBCs (Figure 1C).

• Download figure
• Open in new tab

Figure 1. Transcriptomic profiles of fetal and maternal macrophages and monocytes isolated from term placentas with and without SARS-CoV-2 infection during pregnancy.

(A) Uniform Manifold Approximation and Projection (UMAP) visualization of 31,719 high-quality placental macrophage and monocyte cells enriched from placentas of pregnancies with (N=4) and without (N=8) SARS-CoV-2 infection shows 10 clusters. HBC: Hofbauer cell; PAMM: placenta-associated maternal monocyte/macrophage. (B) Cluster-specific expression of DDX3Y, expressed only in fetal cells, and XIST, expressed only in maternal cells, in placentas from individuals carrying a male fetus (N=10). (C) Correlation of cluster-average gene expression with annotated cell types identified by Suryawanshi et al., Sci Adv, 2018. Each heatmap shows Spearman correlation coefficients. Highest correlation coefficient per cluster is indicated by black dots. HBC clusters were most highly correlated with Suryawanshi HBC clusters, PAMM cluster most correlated with decidual macrophages. (D) Heatmap displaying expression (log₂ fold change) of the top 5 marker genes per cluster. (E) Gene Ontology (GO) Biological Process enrichment analysis for cluster marker genes. GO terms displayed were curated from the top significant GO terms in each cluster, selecting the processes most relevant to macrophage function, and reducing redundancy. Gene Count gives the number of genes in the query set that are annotated by the relevant GO category. GO terms with an adjusted p-value < 0.05 were considered significantly enriched.

To delineate differences in the identity and functions of HBC subclusters, we next assessed top marker genes defining each subcluster, shown as an average heatmap (Figure 1D). Marker genes for HBC clusters 3, 4, 5, and 7 suggested involvement in classic M1 macrophage/pro-inflammatory activities. HBC cluster 3 demonstrated high expression of chemokine (C-X-C motif) ligand genes (CXCLs) as well as pro-inflammatory marker genes IL1B, IL1A, TNF, and NFKB1. HBC cluster 4 demonstrated high expression of multiple CC chemokine ligand genes (CCLs) in a profile similar to that observed in HBCs responding to lipopolysaccharide stimulation in vitro (29). HBC cluster 5 was characterized by high expression of genes encoding major histocompatibility complex (MHC) class II molecules (human leukocyte antigen (HLA)-DRA/B1 and -DP) and Fc-gamma binding protein (FCGBP), suggesting a role in antigen presentation to CD4+ T cells. MHC class II molecules are critical to antigen-specific responses, and upregulation of HLA complexes and antigen presentation pathways has been observed in proteomic analyses of HBCs stimulated with the viral dsRNA mimic poly(I:C) (29). HBC cluster 7 demonstrated marker genes from the interferon-induced protein with tetratricopeptide repeats (IFIT) family (IFIT2, IFIT3), CXCL10, ISG15, and MX1, associated with a pro-inflammatory type 1 interferon antiviral response. To gain further insight into the biological processes reflected in each HBC subcluster, we performed Gene Ontology (GO) enrichment analyses of cluster marker genes (Figure 1E). As would be expected given marker gene expression noted previously, pathways involved in M1-like immune/inflammatory responses were enriched in HBC 3, 4, 5, and 7, including “response to interleukin-1,” “response to interferon-gamma,” “response to tumor necrosis factor,” and “positive regulation of cytokine production.”

Gene signatures of HBC clusters 0, 1, and 2 reflected engagement in specific stress response processes, particularly response to inflammation and/or tissue damage, to ultimately support placental function. HBC 0 and HBC 1 were characterized by genes encoding heat shock proteins and other proteins involved in endoplasmic reticulum stress and the unfolded protein response, such as HSPA6, HSPA1B, DNAJB1, HSP90B1, HSPA5 and BAG3. The unfolded protein response represents a homeostatic response to restore balance when endoplasmic reticulum stress is sensed and to modulate and/or resolve inflammation (30). Additionally, HBC 0’s high expression of PDK4 and KLF2 may suggest involvement in attenuating oxidative stress responses and reducing pro-inflammatory cytokine production (31, 32), and HBC 1’s high expression of FABP5 and HMOX1 suggests engagement in anti-inflammatory responses against heme-induced toxicity and induction towards M2 polarization (33, 34). GO enrichment analysis of these clusters similarly demonstrated enrichment of pathways such as “response to unfolded protein”, “response to heat”, “response to endoplasmic reticulum stress”, and pathways related to cellular stress response and apoptosis (e.g. “extrinsic apoptotic signaling pathway”, “response to oxidative stress” and “ERK1 and ERK2 cascade”).

GO enrichment analysis also suggested both HBC 0 and HBC 1 were engaged in homeostatic functions including “receptor-mediated endocytosis”, “regulation of angiogenesis” and “vascular development”, and nutrient-sensing functions such as “response to nutrient levels” and “response to starvation.” HBC 2 was characterized by high expression of the genes encoding secreted phosphoprotein 1 (SPP1) and Nicotinamide N -methyltransferase (NNMT), both associated with M2 (anti-inflammatory) macrophage polarization in the context of tumor-associated macrophages (35, 36); SPP1, also known as osteopontin, is secreted by HBCs and plays an important role in endothelial biology and angiogenesis (37). GO analysis of HBC cluster 2 also demonstrated enrichment in “receptor-mediated endocytosis” (involved in intracellular transport of macromolecules), as well as “cellular lipid catabolic process,” and processes associated with stromal tissue development.

HBC 6 was characterized by expression of genes involved in regulation of actin polymerization and cytoskeleton organization (TMSB4X and AIF1, which encodes the canonical microglial marker Iba1 (38)) and several ribosomal proteins including RPS18, RPS23, and RPS4Y1. GO enrichment analysis of this cluster demonstrates highly specific enrichment of protein processing pathways (e.g. “protein targeting” pathways, “cytoplasmic translation,” “translational initiation”) and pathways related to RNA catabolism, oxidative phosphorylation, and ATP metabolism. Marker genes of the maternal PAMM subcluster included APOE, APOC1, VIM, LGALS1, and GPNMB among others, an expression pattern consistent with previously reported maternal placental macrophage transcriptional profiles (39, 40). High expression of LGALS1/3 and GPNMB by the PAMM cluster suggests a role in inflammation regulation (41–43), which was echoed by GO analyses identifying enrichment in immune response and immunomodulatory pathways (e.g. “antigen processing and presentation, “positive regulation of cytokine production”, and “regulation of innate immune response”) . In addition, GO enrichment analysis demonstrated PAMM were engaged in lipid metabolic processes and receptor-mediated endocytosis, consistent with their known involvement in lipid engulfment, and tissue repair/scar formation (37, 39).

Maternal SARS-CoV-2 infection drives cluster-specific differences in immune signaling and metabolic pathways

Once the baseline functions of HBC and PAMM subclusters had been established, we then sought to characterize the impact of maternal SARS-CoV-2 infection on the transcriptomic profile of HBC subclusters. To do so, we identified differentially expressed genes (DEG) by maternal SARS-CoV-2 status within each cluster. DEG were defined using log fold-change threshold of 0.2 and adjusted p-value of 0.05 (see Methods for full details). We first verified that each subcluster included representation from both SARS-CoV-2+ cases and controls (Figure 2A, top panel). The proportion of cells from cases versus controls was consistent across all subclusters, except for HBC 0, which demonstrated a relatively high contribution of control donor cells (Figure S2A). Of the 8 HBC clusters, a majority (5) were significantly impacted by maternal SARS-CoV-2 infection: HBC 0, 1, 2, 3, and 5 (Figure 2A, bottom panel). In contrast, HBC clusters 4, 6, and 7 had very few DEG in the setting of maternal SARS-CoV-2 exposure, with three, zero, and two DEG respectively. Of the 5 highly impacted clusters, HBC 1 and HBC 5 had the highest number of DEG by maternal SARS-CoV-2, with 723 and 566 DEG, respectively. PAMMs were impacted by SARS-CoV-2 to a lesser extent, with 67 DEG identified. Both up- and down-regulated DEG were identified across all impacted clusters.

• Download figure
• Open in new tab

Figure 2. Impact of maternal SARS-CoV-2 infection on Hofbauer cell subclusters.

DEG: differentially expressed genes. (A) Barplot demonstrating proportion of cells per cluster from SARS-CoV-2 positive cases (red) and negative controls (gray), top panel. Number of DEG upregulated (dark blue) and downregulated (light blue) by SARS-CoV-2 exposure per cluster, bottom panel. (B) Gene Ontology (GO) Biological Process enrichment analysis for DEG. Gene Count gives the number of genes in the query set that are annotated by the relevant GO category. GO terms with an adjusted p-value < 0.05 were considered significantly enriched. (C) Ingenuity Pathway Analysis (IPA) of DEG for HBC clusters 0 (left) and 1 (right). Canonical pathways with absolute Z-score ≥1 and adjusted p-value < 0.05 are shown. IPA analysis for remaining HBC clusters depicted in Supplement.

GO pathway enrichment analysis of DEG indicated that in the context of maternal SARS-CoV-2 infection, specific pathways involved in immune responses were enriched in all impacted HBC clusters (HBC 0, 1, 2, 3 and 5), including response to lipopolysaccharide, response to interferon-gamma, cellular response to tumor necrosis factor, and regulation of T cell activation (Figure 2B). Additionally, all subclusters were enriched for pathways related to cellular movement, such as cell chemotaxis, leukocyte cell-cell adhesion, and myeloid leukocyte migration; heat shock-related pathways (unfolded protein response); and phagocytosis pathways.

GO enrichment analysis also indicated some biological processes that were only dysregulated in specific clusters in the setting of maternal SARS-CoV-2 infection (Figure 2B). For example, regulation of vascular development was only dysregulated in HBC 0 and the PAMM cluster, regulation of lipid metabolism and transport were dysregulated in all clusters except HBC 5, more cellular energy utilization pathways (e.g. ATP metabolism, electron transport chain/oxidative phosphorylation and cellular respiration) and cellular stress/apoptosis pathways were impacted in HBC clusters relative to the PAMM cluster, and protein processing and actin cytoskeleton organization pathways were only dysregulated in HBC clusters but not PAMMs. Taken together, these functional analyses suggest that in the context of maternal SARS-CoV-2 infection, HBC subclusters and PAMMs are differentially impacted, with key dysregulated biological processes including innate immune and pro-inflammatory signaling, cell chemotaxis and migration, cellular ATP and lipid metabolism, cellular phagocytosis, and the unfolded protein response.

To better understand the impact of SARS-CoV-2 on HBC functions, we next used Ingenuity Pathway Analysis (IPA), which predicts strength and directionality (i.e. activation or suppression) of enriched canonical pathways by subcluster. In this analysis, pathways with absolute Z-score value greater than 1 (consistent with IPA being able to “call” a direction of dysregulation of the pathway) and Benjamini Hochburg-adjusted p<0.05 were included and displayed by subcluster (Figure 2C and Figure S2B-F). Z-scores ≥ 1 indicate upregulated signaling in the pathway and ≤ 1 indicate downregulated pathway signaling (44). Both HBC 1 and HBC 2 subclusters exhibited a primarily anti-inflammatory response to SARS-CoV-2, with activation of PPAR signaling and oxidative phosphorylation, a metabolic profile associated with an anti-inflammatory/pro-resolution phase macrophage signature (45, 46). In HBC 1 and 2, suppression of IL-6, IL-1, and IL-17 pathways, and activation of LXR/RXR signaling pathways in SARS-CoV-2+ cases suggests involvement in resolution of inflammation, as LXR/RXR pathway activation in macrophages is associated with inhibition of inflammatory gene expression and promotion of lipid metabolism (47). Also consistent with an anti-inflammatory role, HBC 2 showed strong suppression of the Coronavirus Pathogenesis Pathway and activation of Oxytocin Signaling Pathway, the latter of which is involved in attenuating oxidative and cellular inflammatory responses in macrophages (48).

Conversely, HBC 0 and 3 demonstrated primarily activated pro-inflammatory immune responses in SARS-CoV-2+ cases, with increases in LPS/IL-1 mediated inhibition of RXR (HBC 3), Interferon induction (HBC 3), Neuroinflammation signaling (HBC 0 and 3), T-cell signaling (HBC 0 and 3), and Production of nitric oxide and reactive oxygen species (HBC 0). Metabolic processes were suppressed in both clusters, including Oxytocin signaling pathway (HBC 0), Sirtuin signaling (HBC 3), and MSP-RON signaling (HBC 0 and 3) (48–50). In the context of maternal SARS-CoV-2 infection, subcluster HBC 5 presented a mixed picture of both pro- and anti-inflammatory signaling, with upregulation of interferon, EIF 2, neuroinflammation and T cell related signaling pathways, balanced by upregulation of anti-inflammatory pathways such as PPAR signaling and downregulation of pro-inflammatory signaling pathways such as Coronavirus Pathogenesis pathway, FAK and TNF-mediated signaling pathways.

Compared to HBC subclusters, PAMMs were less impacted overall by maternal SARS-CoV-2 at a transcriptomic level, with 67 DEG identified. In the setting of maternal SARS-CoV-2 infection, PAMMs showed activation of pathways involved in immune responses including Production of Nitric Oxide and Reactive Oxygen Species, B-cell signaling pathways, Interferon induction, and activation of the pattern recognition receptor TREM-1 signaling. Similar pro-inflammatory patterns were observed for monocytes, including activation of antiviral response pathways and Th1 signaling pathways, and suppression of MSP-RON signaling (Figure S2C). Taken together, these analyses point to transcriptional shifts in some but not all subclusters in response to SARS-CoV-2, with a greater response by HBCs compared to PAMMs, driven by a combination of immune activation/pro-inflammatory signature in subclusters HBC 0 and HBC 3 and an anti-inflammatory tissue repair signature in clusters HBC 1 and HBC 2.

Maternal SARS-CoV-2 infection impacts HBC transcriptional programs associated with phagocytosis, neuroinflammation, and neurological disorders

Tissue-resident macrophages promote resolution of inflammation through phagocytosis of apoptotic cells, invading pathogens, or cellular debris (51, 52). Phagocytosis is also a key function of microglia in early brain development (53–55). IPA functional analysis of SARS-CoV-2-specific HBC signatures demonstrated that macrophage phagocytosis (Figure 3A) and neurological disease-related pathways (Figure 3B) were key functions and pathways implicated by the cluster-specific gene expression signatures. Figure 3 summarizes the impact of maternal SARS-CoV-2 infection on placental macrophage phagocytosis (Figure 3A, 3C), illustrating the potential for altered HBC gene programs to provide insight into both fetal brain microglial function and the impact of maternal SARS-CoV-2 infection on neurodevelopment (Figure 3B, 3D). These analyses predicted SARS-CoV-2-associated suppression of phagocyte chemotaxis and cell movement pathways (e.g. reduced “activation of phagocytes”, “recruitment of phagocytes”, “cell movement of phagocytes”, “adhesion of phagocytes”) in HBC 1, 2 and 5, consistent with the suppression of synaptosome phagocytosis (a proxy for synaptic pruning) observed in subsequent experiments using in vitro Hofbauer cell induced microglial assays, detailed below. In contrast to the consistent suppression of phagocytosis in HBC clusters 1, 2 and 5, HBC clusters 0 and 3 and the PAMM cluster demonstrated activation of phagocytosis-related pathways including “Phagocytosis” (HBC 0), “Immune response of phagocytes” (HBC 0), “Phagocytosis by macrophages” (HBC 3) and “Cellular infiltration by phagocytes” (HBC 3). A representative phagocytosis pathway from IPA and expression of its constituent genes by cluster is depicted as a heatmap in Figure 3C. Cluster-specific alterations in phagocyte movement in the setting of maternal SARS-CoV-2 infection were primarily driven by expression differences in CXCL2, NFKB1A, NFKB1Z, IL1B, CXCL8, CD36, and ICAM1 by cluster (Figure 3C). Concordant with patterns observed in canonical pathways analyses, HBC 1 and 2 (and to a lesser extent HBC 5), which showed primarily immunomodulatory signatures, also show suppressed phagocytosis and phagocytic movement pathways, versus proinflammatory clusters HBC 0 and 3, which demonstrate activation of phagocytosis.

• Download figure
• Open in new tab

Figure 3. Impact of maternal SARS-CoV-2 on HBC gene programs associated with phagocytosis and neurologic disease.

(A-B) Ingenuity Pathway Analysis (IPA) phagocytosis diseases and functions pathways (A), and neurologic diseases or functions (B), enriched for ≥ 3 DEGs, with absolute Z-score ≥1 and adjusted P-value < 0.05. Activation Z-score represented by color and number of DEGs by circle size, with red color indicating pathway activation and blue color indicating suppression. (C-D). Heatmap of gene expression in “Cellular Infiltration by Phagocytes” IPA Pathway (C) and “Inflammation of Central Nervous System” (D) by cluster. Color represents gene expression level (log₂ fold change), *adjusted P-value < 0.05. (E). Module score by subcluster in comparison to cluster-specific gene expression of single-cell datasets from human brain: Microglia and Border Associated Macrophages (Askenase et al., Sci Immunol, 2021) and Yolk Sac Macrophages and Monocytes (Bian et al., Nature, 2020). Color indicates module score.

In addition to phagocytosis, pathways relevant to neurologic disease and microglial functions emerged as key dysregulated pathways in the setting of maternal SARS-CoV-2 infection. We therefore assessed whether DEG of HBC subclusters map to neuroinflammatory/neurodevelopmental pathways and functions in IPA, and plotted pathway activation Z-scores by subcluster (Figure 3B). The transcriptional signature of HBC 1 – in which Fc-gamma receptor-mediated phagocytosis (Figure 2C) and other previously-described phagocytosis pathways (Figure 3A) are suppressed in SARS-CoV-2+ cases – is also associated with increased neuroinflammation, including positive activation Z-scores for “Inflammation of central nervous system”, “Myelitis”, and “Encephalitis.” Cluster-specific expression of the genes in the “Inflammation of central nervous system pathway” is depicted in Figure 3D, with upregulation of signaling in this pathway driven by increased expression of ANXA1, FN1, CCL3, SLC1A3, NLRP3, and

MERTK, among others. Interestingly, HBC 3 – in which “Inflammation of Central Nervous System” is also predicted to be activated after maternal SARS-CoV-2 infection and whose transcriptional signature is consistent with activation in phagocytosis pathways, also implicates increased “Apoptosis of neurons” and “Neuronal cell death” (Figure 3B). In a developmental context, this pattern may represent a functional rather than pathologic gene signature in response to SARS-CoV-2, as microglia (resident brain macrophages) play a key role in neuronal cell turnover, regulation of neural progenitors, and synaptic rewiring in early neurodevelopment, all via phagocytosis (56). Thus, increased phagocytosis by tissue-resident macrophages might be an adaptive response to SARS-CoV-2-associated inflammation, while reduced phagocytosis could be a pathologic or maladaptive response to maternal immune activation (e.g., reduced microglial phagocytosis and reduced synaptic pruning associated with maternal immune activation is thought to be a key aspect of the pathogenesis of autism (57–59)). Taken together, these data support the notion that HBC transcriptional signatures provide insight into protective versus pathologic microglial programming in the setting of an immune challenge such as SARS-CoV-2.

Prior work from our group in a mouse model has shown that HBCs and fetal brain microglia share transcriptional profiles and responses to maternal obesity, an immune-activating exposure (19). To further probe the potential connection between transcriptional signatures of HBC subclusters and brain microglia in humans, we next compared marker genes from HBC subclusters with gene modules from published human single cell atlases of macrophages derived from adult and embryonic brain (Figure 3E) (60, 61). Nearly all HBC subclusters scored highly for gene signatures found in microglia, yolk sac macrophages, or CNS border associated macrophages, compared with monocytes and PAMMs. HBC 2 and 1 exhibited greatest similarity to yolk sac macrophages whereas HBC 3 and 4 were most like microglia isolated from adult brain samples. In contrast, monocytes and PAMMs were most similar to circulating monocytes, which is concordant with their shared myeloid lineage (37). This analysis supports the concept that HBCs isolated from full term human placenta share transcriptional signatures with yolk sac macrophages and fetal brain microglia, and thus may offer insights into global reprogramming of fetal macrophage populations, including those of the fetal brain, in the setting of maternal immune-activating exposures.

HBCs isolated from placentas of SARS-CoV-2 positive pregnancies can be transdifferentiated to microglia-like cells (HBC-iMGs)

To gain insight into the functional consequences of SARS-CoV-2 exposure on HBC populations, we used a previously-validated model in which HBCs isolated from SARS-CoV-2 positive cases (N=4) and a subset of SARS-CoV-2 negative controls (N=4) were cultured in media containing IL-34 and GM-CSF to obtain transdifferentiated microglia-like cells (HBC-iMGs), as we have previously described (see Methods) (62–64). Following culture, we assessed the expression of multiple markers associated with microglial identity, including IBA1, TMEM119, PU.1, and P2RY12 (65, 66), and identified expression of all markers in the majority of HBC-iMGs from both SARS-CoV-2 positive cases and negative controls (Figure 4A, Figure S3A-B).

• Download figure
• Open in new tab

Figure 4. Phenotypic characterization of HBC-iMGs by direct cytokine reprogramming.

HBC-iMGs: Hofbauer cells transdifferentiated toward microglia-like cells. (A) Images of HBC-iMGs from SARS-CoV-2 positive cases (n=4) and negative controls (n=4), immunostained for microglial markers: IBA1, PU.1, P2RY12, TMEM119. Scale bar = 100 µm. (B) Morphology-smoothed density plots (solidity vs. eccentricity) for SARS-CoV-2 negative and positive samples as indicated. Cells from negative controls exhibit a more ramified morphology than those from positive cases, suggestive of a less activated phenotype. Red color shows high density and blue is low density. Representative confocal microscopy images of amoeboid, bipolar, and ramified HBC-iMGs. Scale bar = 50 µm. (C) Violin plots represent distribution of cell solidity (left) and eccentricity (right) measurements from SARS-CoV-2 negative controls (blue, n=5223 cells) and positive cases (orange, n=237 cells). Solid lines represent median values and dashed lines interquartile range. Group differences assessed by Mann-Whitney U Test. ****P<0.0001. ns = not significant.

SARS-CoV-2-exposed HBC-iMGs demonstrate increased amoeboid morphology and impaired phagocytic behavior compared to HBC-iMGs from uninfected control placentas

To assess differences in cell phenotype by SARS-CoV-2 exposure, we first evaluated cellular morphology of IBA1-positive HBC-iMGs by quantitative assessment of two morphological characteristics: eccentricity (amoeboid vs bipolar shape) and solidity (amoeboid/bipolar vs ramified shape) (Figure 4B-C). In this analysis, HBC-iMGs from SARS-CoV-2 negative controls demonstrated a more ramified morphology than those from positive cases, indicated by their lower solidity and high eccentricity (Figure 4B-C). More ramified microglial morphology is generally typical of resting-state, tissue-surveilling microglia in vivo (67, 68). In contrast, a greater proportion of HBC-iMGs generated from SARS-CoV-2 positive cases demonstrated higher solidity and smaller cell size (Figure 4B-C, Figure S3C), consistent with a more amoeboid appearance. While this morphology is classically attributed to an immune-activated state (69), it is also typical of microglial patterns observed in fetal states (67, 70).

Transcriptional analyses of HBC clusters pointed to a cluster-specific impact of maternal SARS-CoV-2 on phagocytosis pathways. HBC clusters with the greatest similarity to embryonic/yolk sac microglia (e.g. HBC 1, 2) also exhibited cellular programs suggestive of impaired phagocytosis. Using a previously-validated model of synaptic pruning (62–64), a key physiologic function of microglia in early brain development, we next tested the functional capability of HBC-iMGs to engage in phagocytosis. In this assay, HBC-iMGs were co-cultured for 3 hours with pHrodo Red-labeled neuronal synaptosomes derived from human induced pluripotent stem cells prior to fixation and imaging. This pH-sensitive label fluoresces following intracellular engulfment (see Methods). Synaptosome engulfment by IBA1-positive cells was then measured by quantifying fluorescence using confocal microscopy images with CellProfiler software applied for segmentation and thresholding (Figure 5A). Compared to SARS-CoV-2 negative cases, HBC-iMGs from positive cases demonstrated significant impairment in synaptosome phagocytosis, reflected by a reduced phagocytic index (Figure 5B). Phagocytic index was reduced across all SARS-CoV-2 positive samples, and was driven by reduced phagocytic uptake per cell, regardless of the proportion of cells engaged in phagocytosis in any given sample. (Figure S3D-E). In conjunction with the morphologic changes, and consistent with the transcriptomic signatures observed in a subset of HBC, these functional phenotypes support a dysregulated activation state following SARS-CoV-2 infection.

• Download figure
• Open in new tab

Figure 5. Impact of maternal SARS-CoV-2 on HBC-iMG synaptosome engulfment.

HBC-iMGs: Hofbauer cells transdifferentiated toward microglia-like cells. (A) Representative image showing colocalization of pHrodo-red labeled synaptosomes (SYN) and IBA1 positive HBC-iMGs. Hoechst = nuclear stain. Scale bar = 100 µm. (B) Violin plots of phagocytic index of image fields from SARS-CoV-2 negative controls (blue, n=187 fields) and positive cases (orange, n=32 fields). Phagocytic index is calculated as synaptosome area in pixels divided by cell count per image field. Solid lines represent median values and dashed lines interquartile range. Group differences assessed by Mann-Whitney U Test. ****P<0.0001.

Study design and participant enrollment

In this study, 12 pregnant individuals with full-term, singleton pregnancies delivering at Massachusetts General Hospital (March 2021 - August 2021) were included. Participants were classified as SARS-CoV-2 positive (N=4) if they had symptomatic COVID-19 infection during pregnancy, confirmed by positive SARS-CoV-2 nasopharyngeal PCR test. Participants were classified as SARS-CoV-2 negative (n=8) if they did not have a positive SARS-CoV-2 nasopharyngeal PCR or COVID-19 symptoms during pregnancy and had a negative SARS-CoV-2 nasopharyngeal PCR at delivery upon universal COVID-19 screening on Labor and Delivery. Pregnant individuals were eligible for inclusion if they were 18 years or older and were delivering during the COVID-19 pandemic. For this study, individuals with prior COVID-19 vaccination were excluded. A study questionnaire and review of the electronic health record was used to determine key demographic variables such as maternal age, gestational age at delivery, gestational age at positive COVID-19 test, COVID-19 disease severity at diagnosis, any prior diagnoses of COVID-19 or history of COVID-19 vaccination, and infant sex and birthweight.

Placenta collection and processing

Placentas were obtained within 20 minutes after delivery and submerged in Cytowash media (Dulbecco’s Modified Eagle Medium (DMEM) containing 2.5% FBS, 1% Penicillin-Streptomycin, 0.1% Gentamicin) and stored at 4°C until cell isolation. Isolation of Hofbauer cells was performed using previously described protocols (26); reagents are listed in Supplemental Table S1 and isolation workflow and study procedures are depicted in Supplemental Figure S4. Briefly, placental chorionic villi were separated from fetal membranes and decidua, washed in DPBS wash, and mechanically homogenized. Placental tissue was then serially digested in Collagenase Digestion Buffer, Trypsin Digestion Buffer, and Collagenase Digestion Buffer 2. Undigested tissue was removed by passage through sterile gauze and 100uM filter. The cell suspension was centrifuged at 257g for 8 min at 4°C, washed, spun again, and resuspended in media. The cells were then suspended in 4mL of 20% Percoll and 5mL of 35% Percoll was underlayed through a #1 glass Pasteur pipette (92). After centrifugation for 30 minutes at 4°C without brake at 1000g, cells were isolated from the interphase layer, washed in media, and spun at 257g for 8 minutes at 4°C. Cell pellets were immunopurified by negative selection by incubation with anti-EGFR (to remove syncytiotrophoblasts) and anti-CD10 (to remove fibroblasts) conjugated to magnetic Dynabeads, prepared as previously described (76), for 20 minutes at 4°C. Tubes were placed on a DynaMagTM magnet for 5 minutes to magnetically bind cytotrophoblasts (anti-EGFR) and fibroblasts (anti-CD10) – allowing media containing unbound placental macrophages to pass through into collection tubes. Cells were cryopreserved in 90% FBS and 10% dimethyl sulfoxide (DMSO) at 1-10 million cells/vial and stored at -80°C for downstream analyses. SARS-CoV-2 viral loads were assessed in all placental tissues using qPCR as previously described, with 40 copies/mL as limit of detection (24, 25).

Single-cell RNA-sequencing (scRNA-Seq) data analysis

Derivation of Hofbauer cells transdifferentiated toward microglia-like cells (HBC-iMGs) by direct cytokine reprogramming

HBC-iMGs were derived from HBCs using previously described methods (62–64), with modifications as noted. Briefly, thawed HBCs were plated on Geltrex-coated 24-well plates (1 × 10⁶ cells in 0.5 mL per well) or 96-well plates (2 x 10⁵ cells in 0.1mL per well) depending on cell availability. After cells were incubated at 37°C for 24 h, the media was completely replaced with RPMI 1640 including GlutaMAX, 1% penicillin–streptomycin, 100 ng/mL of human recombinant IL-34 (Peprotech), and 10 ng/mL of GM-CSF (Peprotech). At day 6 of transdifferentiation, the cultures were assayed and subsequently fixed with 4% PFA to perform endpoint analysis using immunocytochemistry.

HBC-iMG Immunocytochemistry

HBC-iMGs were washed twice with PBS and blocked for 1 h with 5% FBS and 0.3% Triton-X (Sigma Aldrich) in PBS at room temperature. Next, they were washed three times with 1% FBS in PBS and incubated with primary antibodies in 5% FBS and 0.1% Triton-X overnight at 4°C (Anti-IBA1, 1:500; Abcam #ab5076; Anti-TMEM119, 1:500; Abcam #ab18537; Anti-CX3CR1, 1:100, Abcam #ab8021; Anti-PU.1, 1:1000, Abcam #ab183327, and Anti-P2RY12, 1:100, Alomone Labs). Cells were then washed three times with 1% FBS in PBS and incubated in secondary antibodies (Invitrogen Alexa Fluor, 1:500) and Hoechst 33342 (1:5000) in 5% FBS and 0.1% Triton-X in PBS for 45 min at 4°C. Cells were washed two final times and imaged using the IN Cell Analyzer 6000 (Cytiva). Marker characterization was analyzed using CellProfiler (94). Cells were segmented using one of the four microglia markers used and percent marker positive cells calculated by dividing the number of marker positive cells by the number of identified nuclei, per image. A total of 12 20x images per sample were analyzed.

Synaptosome derivation and phagocytosis assay

Image analysis

CellProfiler (v4.2.4) was used to identify cellular and subcellular structures in the confocal images (94). The module CorrectIlluminationCalculate and CorrectIlluminationApply were used in all channels to correct uneven illumination and uneven background. Nuclei and cell bodies were each identified using IdentifyPrimaryObjects. Specifically, pixel diameter ranges and the automatic thresholding method Otsu were applied. The module RelateObjects was used to drop structures incorrectly identified as nuclei by ensuring they were only accepted when they have a surrounding microglia-like cell. IdentifySecondaryObjects was used to more accurately outline cells around these nuclei and avoid debris. The module IdentifyTertiaryObjects identified cytoplasm by subtracting the area of the nucleus from the cell. IdentifySecondaryObjects was used with Otsu thresholding to identify and measure synaptosomes. Background red signal was eliminated by increasing the lower bounds on the automatic threshold, using reference images from Cytochalasin treatment as a positive control of diminished phagocytosis. MaskObjects was used with the cell and synaptosome objects to omit red signal from outside the cell. Overlays of the outlines of all generated image structures were created for quality check purposes using the OverlayOutlines module. Cell area, count, and signal intensity were created with MeasureObjectSizeShape, MeasureObjectIntensity, ExportToSpreadSheet, and ExportToDatabase. RStudio 2 (1.4) was used to organize the metadata exported from CellProfiler. Phagocytic index was calculated as area of Synaptosomes divided by cell count per image. As a confirmation, the integrated intensity of Synaptosome signal divided by cell area was also checked to make sure both measures corresponded. Images containing >80 cells were omitted due to procedure inaccuracy with dense fields. Outliers were excluded by calculating a phagocytic index threshold of 3 SD above the mean. Morphology data was produced using cell level metadata from CellProfiler followed by a cleaning process matching the field level dataset cleaning. Cells with an area or synaptosome area of greater than the mean plus 3 SD were omitted. In total, 12 20x images per sample were analyzed.

Statistics

Group differences were assessed using Mann-Whitney U tests. P values less than 0.05 were considered statistically significant. Dark lines represent median and dotted lines interquartile range, unless otherwise specified. Statistical analyses were performed in GraphPad Prism (version 9.3).