SARS-CoV-2 niches in human placenta revealed by spatial transcriptomics

Fetal MAVS and type I IFN signaling pathways control ZIKV infection in the placenta and maternal decidua

The contribution of placental immune responses to congenital Zika virus (ZIKV) syndrome remains poorly understood. Here, we leveraged a mouse model of ZIKV infection to identify mechanisms of innate immune restriction exclusively in the fetal compartment of the placenta. ZIKV principally infected mononuclear trophoblasts in the junctional zone, which was limited by mitochondrial antiviral-signaling protein (MAVS) and type I interferon (IFN) signaling mechanisms. Single nuclear RNA sequencing revealed MAVS-dependent expression of IFN-stimulated genes (ISGs) in spongiotrophoblasts but not in other placental cells that use alternate pathways to induce ISGs. ZIKV infection of Ifnar1-/- or Mavs-/- placentas was associated with greater infection of the adjacent immunocompetent decidua, and heterozygous Mavs+/- or Ifnar1+/- dams carrying immunodeficient fetuses sustained greater maternal viremia and tissue infection than dams carrying wild-type fetuses. Thus, MAVS-IFN signaling in the fetus restricts ZIKV infection in junctional zone trophoblasts, which modulates dissemination and outcome for both the fetus and the pregnant mother.

Discrete placental gene expression signatures accompany diabetic disease classifications during pregnancy

BACKGROUND

Gestational diabetes mellitus affects up to 10% of pregnancies and is classified into subtypes gestational diabetes subtype A1 (GDMA1) (managed by lifestyle modifications) and gestational diabetes subtype A2 (GDMA2) (requiring medication). However, whether these subtypes are distinct clinical entities or more reflective of an extended spectrum of normal pregnancy endocrine physiology remains unclear.

OBJECTIVE

Integrated bulk RNA-sequencing (RNA-seq), single-cell RNA-sequencing (scRNA-seq), and spatial transcriptomics harbors the potential to reveal disease gene signatures in subsets of cells and tissue microenvironments. We aimed to combine these high-resolution technologies with rigorous classification of diabetes subtypes in pregnancy. We hypothesized that differences between preexisting type 2 and gestational diabetes subtypes would be associated with altered gene expression profiles in specific placental cell populations.

STUDY DESIGN

In a large case-cohort design, we compared validated cases of GDMA1, GDMA2, and type 2 diabetes mellitus (T2DM) to healthy controls by bulk RNA-seq (n=54). Quantitative analyses with reverse transcription and quantitative PCR of presumptive genes of significant interest were undertaken in an independent and nonoverlapping validation cohort of similarly well-characterized cases and controls (n=122). Additional integrated analyses of term placental single-cell, single-nuclei, and spatial transcriptomics data enabled us to determine the cellular subpopulations and niches that aligned with the GDMA1, GDMA2, and T2DM gene expression signatures at higher resolution and with greater confidence.

RESULTS

Dimensional reduction of the bulk RNA-seq data revealed that the most common source of placental gene expression variation was the diabetic disease subtype. Relative to controls, we found 2052 unique and significantly differentially expressed genes (-22 thresholds; q<0.05 Wald Test) among GDMA1 placental specimens, 267 among GDMA2, and 1520 among T2DM. Several candidate marker genes (chorionic somatomammotropin hormone 1 [CSH1], period circadian regulator 1 [PER1], phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit beta [PIK3CB], forkhead box O1 [FOXO1], epidermal growth factor receptor [EGFR], interleukin 2 receptor subunit beta [IL2RB], superoxide dismutase 3 [SOD3], dedicator of cytokinesis 5 [DOCK5], suppressor of glucose, and autophagy associated 1 [SOGA1]) were validated in an independent and nonoverlapping validation cohort (q<0.05 Tukey). Functional enrichment revealed the pathways and genes most impacted for each diabetes subtype, and the degree of proximal similarity to other subclassifications. Surprisingly, GDMA1 and T2DM placental signatures were more alike by virtue of increased expression of chromatin remodeling and epigenetic regulation genes, while albumin was the top marker for GDMA2 with increased expression of placental genes in the wound healing pathway. Assessment of these gene signatures in single-cell, single-nuclei, and spatial transcriptomics data revealed high specificity and variability by placental cell and microarchitecture types. For example, at the cellular and spatial (eg, microarchitectural) levels, distinguishing features were observed in extravillous trophoblasts (GDMA1) and macrophages (GDMA2). Lastly, we utilized these data to train and evaluate 4 machine learning models to estimate our confidence in predicting the control or diabetes status of placental transcriptome specimens with no available clinical metadata.

CONCLUSION

Consistent with the distinct association of perinatal outcome risk, placentae from GDMA1, GDMA2, and T2DM-affected pregnancies harbor unique gene signatures that can be further distinguished by altered placental cellular subtypes and microarchitectural niches.

Zika virus co-opts microRNA networks to persist in placental niches detected by spatial transcriptomics

BACKGROUND

Zika virus congenital infection evades double-stranded RNA detection and may persist in the placenta for the duration of pregnancy without accompanying overt histopathologic inflammation. Understanding how viruses can persist and replicate in the placenta without causing overt cellular or tissue damage is fundamental to deciphering mechanisms of maternal-fetal vertical transmission.

OBJECTIVE

Placenta-specific microRNAs are believed to be a tenet of viral resistance at the maternal-fetal interface. We aimed to test the hypothesis that the Zika virus functionally disrupts placental microRNAs, enabling viral persistence and fetal pathogenesis.

STUDY DESIGN

To test this hypothesis, we used orthogonal approaches in human and murine experimental models. In primary human trophoblast cultures (n=5 donor placentae), we performed Argonaute high-throughput sequencing ultraviolet-crosslinking and immunoprecipitation to identify any significant alterations in the functional loading of microRNAs and their targets onto the RNA-induced silencing complex. Trophoblasts from same-donors were split and infected with a contemporary first-passage Zika virus strain HN16 (multiplicity of infection=1 plaque forming unit per cell) or mock infected. To functionally cross-validate microRNA-messenger RNA interactions, we compared our Argonaute high-throughput sequencing ultraviolet-crosslinking and immunoprecipitation results with an independent analysis of published bulk RNA-sequencing data from human placental disk specimens (n=3 subjects; Zika virus positive in first, second, or third trimester, CD45- cells sorted by flow cytometry) and compared it with uninfected controls (n=2 subjects). To investigate the importance of these microRNA and RNA interference networks in Zika virus pathogenesis, we used a gnotobiotic mouse model uniquely susceptible to the Zika virus. We evaluated if small-molecule enhancement of microRNA and RNA interference pathways with enoxacin influenced Zika virus pathogenesis (n=20 dams total yielding 187 fetal specimens). Lastly, placentae (n=14 total) from this mouse model were analyzed with Visium spatial transcriptomics (9743 spatial transcriptomes) to identify potential Zika virus-associated alterations in immune microenvironments.

RESULTS

We found that Zika virus infection of primary human trophoblast cells led to an unexpected disruption of placental microRNA regulation networks. When compared with uninfected controls, Zika virus-infected placentae had significantly altered SLC12A8, SDK1, and VLDLR RNA-induced silencing complex loading and transcript levels (-22; adjusted P value <.05; Wald-test with false discovery rate correction q<0.05). In silico microRNA target analyses revealed that 26 of 119 transcripts (22%) in the transforming growth factor-β signaling pathway were targeted by microRNAs that were found to be dysregulated following Zika virus infection in trophoblasts. In gnotobiotic mice, relative to mock controls, Zika virus-associated fetal pathogenesis included fetal growth restriction (P=.036) and viral persistence in placental tissue (P=.011). Moreover, spatial transcriptomics of murine placentae revealed that Zika virus-specific placental niches were defined by significant up-regulation of complement cascade components and coordinated changes in transforming growth factor-β gene expression. Finally, treatment of Zika virus-infected mice with enoxacin abolished placental Zika virus persistence, rescued the associated fetal growth restriction, and the Zika virus-associated transcriptional changes in placental immune microenvironments were no longer observed.

CONCLUSION

These results collectively suggest that (1) Zika virus infection and persistence is associated with functionally perturbed microRNA and RNA interference pathways specifically related to immune regulation in placental microenvironments and (2) enhancement of placental microRNA and RNA interference pathways in mice rescued Zika virus-associated pathogenesis, specifically persistence of viral transcripts in placental microenvironments and fetal growth restriction.

Maternal SARS-CoV-2 infection in pregnancy disrupts gene expression in Hofbauer cells with limited impact on cytotrophoblasts

BACKGROUND

Hofbauer cells (HBCs) and cytotrophoblasts (CTBs) are major cell populations in placenta. The indirect impact of maternal SARS-CoV-2 disease on these cells that are not directly infected has not been extensively studied. Herein, we profiled gene expression in HBCs and CTBs isolated from placentae of recovered pregnant subjects infected with SARS-CoV-2 during all trimesters of pregnancy, placentae from subjects with active infection, SARS-CoV-2 vaccinated subjects, and those who were unexposed to the virus.

METHODS

Placentae were collected within 4 h post-delivery and membrane-free tissues were enzymatically digested for the isolation of HBCs and CTBs. RNA extracted from HBCs and CTBs were sequenced using 150bp paired-end reads. Differentially expressed genes (DEGs) were identified by DESeq2 package in R and enriched in GO Biological Processes, KEGG Pathway, Reactome Gene Sets, Hallmark Gene Sets, and Canonical Pathways. Protein-protein interactions among the DEGs were modelled using STRING and BioGrid.

RESULTS

Pregnant subjects (n = 30) were recruited and categorized into six groups: infected with SARS-CoV-2 in i) the first (1T, n = 4), ii) second (2T, n = 5), iii) third (3T, n = 5) trimester, iv) tested positive at delivery (Delivery, n = 5), v) never infected (Control, n = 6), and vi) fully mRNA-vaccinated by delivery (Vaccinated, n = 5). Compared to the Control group, gene expression analysis showed that HBCs from infected subjects had significantly altered gene expression profiles, with the 2T group having the highest number of DEGs (1,696), followed by 3T and 1T groups (1,656 and 958 DEGs, respectively). These DEGs were enriched for pathways involved in immune regulation for host defense, including production of cytokines, chemokines, antimicrobial proteins, ribosomal assembly, neutrophil degranulation inflammation, morphogenesis, and cell migration/adhesion. Protein-protein interaction analysis mapped these DEGs with oxidative phosphorylation, translation, extracellular matrix organization, and type I interferon signaling. Only 95, 23, and 8 DEGs were identified in CTBs of 1T, 2T, and 3T groups, respectively. Similarly, 11 and 3 DEGs were identified in CTBs and HBCs of vaccinated subjects, respectively. Reassuringly, mRNA vaccination did not induce an inflammatory response in placental cells.

CONCLUSIONS

Our studies demonstrate a significant impact of indirect SARS-CoV-2 infection on gene expression of inner mesenchymal HBCs, with limited effect on lining CTB cells isolated from pregnant subjects infected and recovered from SARS-CoV-2. The pathways associated with these DEGs identify potential targets for therapeutic intervention.