key: cord-0842062-ajeggpn9
authors: Matute, Juan D.; Finander, Benjamin; Pepin, David; Ai, Xingbin; Smith, Neal P.; Li, Jonathan Z.; Edlow, Andrea G.; Villani, Alexandra-Chloe; Lerou, Paul H.; Kalish, Brian T.
title: Single-cell immunophenotyping of the fetal immune response to maternal SARS-CoV-2 infection in late gestation
date: 2021-11-08
journal: Pediatr Res
DOI: 10.1038/s41390-021-01793-z
sha: e00c48af3d2d380a188a8e72ac2ad63fc2912baa
doc_id: 842062
cord_uid: ajeggpn9

BACKGROUND: During the COVID-19 pandemic, thousands of pregnant women have been infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The implications of maternal SARS-CoV-2 infection on fetal and childhood well-being need to be characterized. We aimed to characterize the fetal immune response to maternal SARS-CoV-2 infection. METHODS: We performed single-cell RNA-sequencing and T cell receptor sequencing on cord blood mononuclear cells (CBMCs) from newborns of mothers infected with SARS-CoV-2 in the third trimester (cases) or without SARS-CoV-2 infection (controls). RESULTS: We identified widespread gene expression changes in CBMCs from cases, including upregulation of interferon-stimulated genes and major histocompatibility complex genes in CD14(+) monocytes, transcriptional changes suggestive of activation of plasmacytoid dendritic cells, and activation and exhaustion of natural killer cells. Lastly, we observed fetal T cell clonal expansion in cases compared to controls. CONCLUSIONS: As none of the infants were infected with SARS-CoV-2, our results suggest that maternal SARS-CoV-2 infection might modulate the fetal immune system in the absence of vertical transmission. IMPACT: The implications of maternal SARS-CoV-2 infection in the absence of vertical transmission on fetal and childhood well-being are poorly understood. Maternal SARS-CoV-2 infection might modulate the fetal immune system in the absence of vertical transmission. This study raises important questions about the untoward effects of maternal SARS-CoV-2 on the fetus, even in the absence of vertical transmission.

Millions of people worldwide have or will become infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), causing Coronavirus Disease 2019 (COVID- 19) , and the infection of pregnant women with SARS-CoV-2 infection has been widespread. [1] [2] [3] [4] [5] [6] Despite the prevalence of antepartum infection, we have a limited understanding of the implications of SARS-CoV-2 infection on fetal and offspring health. To date, there are limited case reports of vertical, mother-to-child transmission of SARS-CoV-2, 7-12 and vertical transmission remains rare in most pregnancies complicated by maternal SARS-CoV-2 infection. 1, 6, [13] [14] [15] [16] [17] Nonetheless, in the absence of direct fetal infection and toxicity, maternal SARS-CoV-2 infection may still affect fetal development. Maternal immune activation during pregnancy after viral infection without vertical transmission can have long-term consequences for the newborn, including abnormal neurologic 18, 19 or immune system development. 20, 21 Pregnancy is a complex immunologic state, and there is no data on the effect of SARS-CoV-2-dysregulated immune state during pregnancy on the fetus. Given the number of pregnant women infected with SARS-CoV-2 worldwide, it is important to determine the potential transgenerational implications of infection with SARS-CoV-2 during pregnancy beyond vertical transmission. Maternal SARS-CoV-2 test positivity has been significantly associated with admission for neonatal care and with neonatal morbidities such as respiratory distress syndrome and hyperbilirubinemia. 22 To date, research on the implications of SARS-CoV-2 infection during pregnancy on the offspring immune system has been limited to antibody transplacental transfer 6, 23, 24 and postnatal evaluation of infants born to mothers infected with SARS-CoV-2 during pregnancy without a nonexposed control group. 25, 26 The latter may be also confounded by ex utero determinants of immune development during the first week of life. 27, 28 In the present study, we characterize the composition and cell-type-specific transcriptional landscape of umbilical cord blood mononuclear cells (CBMCs) from term gestation infants (>37 weeks) born to mothers infected with SARS-CoV-2 in the third trimester without vertical transmission. This immunogenomic investigation provides evidence of both innate and adaptive fetal immune transcriptional changes in pregnancies complicated by SARS-CoV-2 infection. Our results suggest that even in the absence of vertical transmission, SARS-CoV-2 maternal infection in the third trimester might modulate the fetal immune system.

To characterize the fetal immunologic landscape in pregnancies complicated by maternal SARS-CoV-2 infection, we performed droplet-based single-cell RNA-sequencing (scRNAseq) of CBMCs from infants born to mothers with SARS-CoV-2 infection during pregnancy (cases) and infants born to mothers without SARS-CoV-2 infection (controls). CBMCs from three cases and three controls were obtained from our biorepository. 29 None of the three infants in this study born to mothers with SARS-CoV-2 were positive for SARS-CoV-2 postnatally, had detectable SARS-CoV-2 messenger RNA (mRNA) in the placenta, or developed any neonatal morbidity. All mothers with COVID-19 in the third trimester were classified as having mild disease without respiratory support. 30 Infants born to mothers negative for SARS-CoV-2 and asymptomatic (universal screening at admission for labor) during the same epoch served as controls. Maternal comorbidities including wellcontrolled thyroid dysfunction, obesity, or gestational diabetes were matched between cases and controls as feasible. The time of maternal infection and birth in cases varied between 7 and 66 days. Table 1 displays demographic and clinical data from the cases and controls.

CBMCs were processed on the 10X Genomics Single-Cell Immune platform (see "Methods"). After quality control and doublet removal, we included 14,748 cells with high-quality single-cell transcriptomes from cases and 11,222 cells from controls in our dataset. (See quality control metrics in Supplemental Fig. 1A, B .) The cell population composition was visualized using uniform manifold approximation and projection (Fig. 1a) , and cell types were inferred by cluster-specific canonical marker genes (Fig. 1b, c) . We did not observe any differences in cell cluster composition between cases and controls (Supplemental Fig. 1c, d) .

To explore transcriptional signatures in fetal immune cells associated with maternal SARS-CoV-2 infection, we performed differential gene expression (DGE) analysis within cell types comparing cases and controls. Genes with a false discovery rate (FDR) < 5% were considered statistically significant. We identified hundreds of genes across nearly all cell types with altered expression (Fig. 2a) . We used gene ontology (GO) analysis to broadly classify genes significantly disrupted by maternal SARS-CoV-2 infection based on DGE (Supplemental Table 1 ).

CD14 + monocytes were grouped into five clusters and CD16 + monocytes were grouped into one cluster (Fig. 2b) . CD14 + subpopulations demonstrated variable expression of inflammatory genes, including ACSL1, ADGRE2, CD300E, and PADI4, which aligns with prior single-cell analysis showing monocyte diversity. 31 Consistent with data from adult COVID-19 patients, 32, 33 we found that CD14 + monocytes from cases demonstrated increased expression of interferon (IFN)-stimulated genes (ISGs) (Fig. 2c ) and concomitant IFNAR2 downregulation (Supplemental Fig. 2a) , which could reflect exposure to IFN prenatally. 34 Of note, we found that there was variable upregulation of ISGs among various CD14 + monocyte clusters in cases compared to controls (Fig. 2c) , indicating that there is no homogeneous ISG regulation among fetal monocyte clusters secondary to maternal SARS-CoV2 infection. This finding is consistent with previous reports that upon stimulation with IFN-γ, fetal bone marrow monocytes have variable upregulation of ISGs when compared to adult bone marrow monocytes. 35 GO analysis of DGE in CD14 + monocytes demonstrated enrichment of genes associated with antigen presentation and viral translational termination and reinitiation (Fig. 2d) . CB CD14 + monocytes from cases also showed upregulation of major histocompatibility class (MHC) I and II genes (Supplemental Fig. 2A ) suggesting activation in response to IFN signaling. 36 However, some HLA genes were downregulated in CD14 + monocytes (HLA-E and HLA-B) (Supplemental Fig. 2A ) consistent with decreased antigen presentation capacities in fetal bone marrow CD14 + monocytes compared to adult bone marrow counterparts upon in vitro stimulation with IFN-γ. 35 Furthermore, CD14 + monocytes from cases showed upregulation of TLR receptor transcripts (TLR2, TLR4, and TLR5) paired with upregulation of FOS and downregulation of transcriptional inhibitors of nuclear factor kappa-light-chain-enhancer of activated B cells (NFKB), including NFKBIA and NFKBIE, all of which are associated with increased NFKB activation and cytokine production 37 (Supplemental Fig. 2A ). Of note, CD14 + monocytes from cases had decreased expression of autophagy (ATG14, ATG2A, and ATG3) and endoplasmic reticulum stress (XBP1 and HSPA5) genes, which may contribute to a defect in macrophage differentiation 38 (Supplemental Fig. 2A ). Expression of S100A Alarmins in CB monocytes has been associated with chorioamnionitis and fetal inflammatory response syndrome (FIRS). 39 S100A8, S100A9, and S100A12 were found to be decreased in CD14 + monocytes in our cases compared to controls, which might suggest that the response elicited from SARS-CoV-2 maternal infection in monocytes differs from the changes elicited in FIRS secondary to chorioamnionitis (Supplemental Fig. 2A ).

Similar to CD14 + monocytes, we identified induction of ISGs in nonclassical CB monocytes (CD16 + ) (Fig. 2c) . In contrast to CD14 + monocytes, we found that there was decreased expression of cell adhesion genes (including PLAUR and THBS1), attenuation of immune activation signaling pathways genes (FOS, FOSB, MAP3K8, STAT6, and FCER1G), and decreased expression of inflammatory molecules like resistin (RETN) (Supplemental Fig. 2B ). Together, these results suggest induction of ISGs in monocytes from cases compared to controls and differences in transcriptional changes in classical and nonclassical monocytes that might indicate preferential activation of classical monocytes in cases compared to controls.

We captured the transcriptomes of both plasmacytoid and conventional dendritic cells (pDC and cDC, respectively) in CB. In adults infected with SARS-CoV-2, both types of DCs are functionally impaired, and there is an increased ratio of cDCs to pDCs in severe patients. 40 In our study, CB cDC from cases showed increased expression of ISGs like IFITM3 and APOBEC3A (Fig. 2e ). Fetal cDC from cases showed a transcriptional profile suggestive of innate immune activation including increased expression of PIK3CB, which is downstream of TLR5 and TLR7, 41 as well as increased transcription of CCL5, which can be upregulated after TLR3 stimulation (Fig. 2e ). 42 Evidence of impaired cDC maturation was suggested by upregulation of ID1, which antagonizes dendritic cell differentiation and antitumor immunity in mice, 43 as well as increased MAFB transcription, which suppresses cDC maturation. 44 cDCs from cases also demonstrated decreased expression of FOSB and many MHC II genes. 45 pDCs in cases also showed markers of immune activation, including upregulation of RELB, which promotes DC activation through RelB-p50 dimer, 46 upregulation of MHC class I and class II genes, and unfolded protein response activation, as shown by increased transcription of XBP1 47 (Supplemental Fig. 2C ). Together, these transcriptional findings could be consistent with activation of pDC over cDC in the CB of cases, potentially through activation of TLRs.

In adults, SARS-CoV-2 infection is associated with fewer blood natural killer (NK) cells, but a higher activation state in circulating NK cells. 48 We identified two clusters of CB NK cells. One population of NK cells (cluster 1) expressed higher levels of GZMB, while the second population of NK cells (cluster 2) expressed IL7R and XCL1, suggesting that cluster 1 corresponded to CD56dim and cluster 2 corresponded to CD56bright NK cells, as NCAM1 (CD56) is technically not well captured in scRNAseq. 49 Similar to adult NK cells, CB NK cells from SARS-CoV-2-positive pregnancies showed signs of exposure to IFN, including induction of ISGs like IFI6, IFIT2, and IRF9 (Fig. 2f) . 48, 50, 51 We identified increased transcription of CCL4, expression of cytotoxic genes including GNLY, GZMA, GZMB, and GZMH, and increased transcription of IFNG, paired with decreased expression of NK inhibitory molecules (Fig. 2f) . 48, 50, 51 There were transcriptional changes associated with exhaustion, such as decreased expression of KLRG1 and SIGLEC7. 52 DGE in NK cells between cases and controls were enriched for genes related to the IFN-α response, regulation of NK cell cytokine production, and viral transcription (Fig. 2g) .

In adults with acute COVID-19, there is a heterogeneous adaptive immune response in peripheral blood, including B cell receptor and T cell receptor (TCR) arrangements specific to SARS-CoV-2. 50 Given these findings, we evaluated whether maternal infection with SARS-CoV-2 had any effect on CB lymphocyte gene expression. In B cells from infants exposed to SARS-CoV-2 in utero, we identified three clusters of CB B cells corresponding to nonplasma/plasmablast (Clusters 1 and 2) and plasma/plasmablast cells (Cluster 3) based on MZB1 expression 50 (Supplemental Fig. 2D ). With the possible exception of germinal center (GC) B cells, all human B cell populations found in lymphoid tissues can also be demonstrated in peripheral blood 53 with enrichment of mature, naive B cells and regulatory B cells in CB. [54] [55] [56] Therefore, transcripts of genes related to B cell receptor activation are detected in circulating B cells [57] [58] [59] and their products might participate in tonic B cell signaling under homeostatic conditions in peripheral blood. 60 Human CB B cell transcriptional program differs from adult peripheral B cell and confers accelerated but transient responsiveness to stimulation. 56 Surprisingly, in B cells from infants born to mothers infected with SARS-CoV-2, we identified decreased transcription of genes downstream of B cell receptors in all clusters compared to controls. Specifically, we found a decreased transcription of NR4A1, CD69, and CD83 in all B cells (Fig. 3a) . NR4A1 encodes Nur77, an orphan nuclear receptor that is induced upon B cell activation in peripheral blood in humans. 61 CD83 is expressed in peripheral B cells and correlates with the engagement of BCR, TLRs, or CD40. 62, 63 Concordant with transcriptional evidence of defective B cell activation in CB lymphocytes from cases, we also found decreased expression of CD69, 64 activator protein-1, and nuclear factor of activated T cells genes, 65 as well as anti-apoptotic genes, including BCL2 and BCL2A1 66 (Fig. 3a) . These transcriptional changes suggestive of potential B cell dysfunction, combined with our prior findings of decreased transplacental transmission of IgG against SARS-CoV-2 compared to IgG against other antigens, 6,23 might translate into potential impairments in antibody-mediated immunity to SARS-CoV-2 in neonates born to mothers with COVID-19. The humoral defect might be secondary to impaired passive immunity against SARS-CoV2 combined with potential impaired B cell activation in cases compared to controls, but this phenomenon requires further functional characterization. We also identified three clusters of T cells (Fig. 3b) . Cluster 1 corresponded to cytotoxic (CD8 + ) T cells, while Clusters 2 and 3 corresponded to helper T cells. 57 Naive helper T cells predominate in CB [67] [68] [69] and retain a partial expression of a fetal-associated regulatory T cell signature. 67, [70] [71] [72] Increased expression of CCR7 in T cell Cluster 2 suggests that this cluster includes either naive T cells or central memory T cells, 73 and increased CTSW and KLRB1 in Cluster 3 suggest that this cluster includes effector and memory T cells. 74, 75 In adults with COVID-19, CD8 + T cells show decreased cytotoxic potential and exhaustion driven by interleukin-6 (IL-6). 76 Similar to adults, CB CD8 + T cells from cases demonstrated transcriptional signatures suggestive of impaired function including decreased expression of GZMA, FOS, and JUN 77-79 (Supplemental Fig. 3A) . Furthermore, there was increased expression of KLRB1 and CCR7 (Supplemental Fig. 3A ). CCR7 is expressed in naive and central memory T cells and lost in effector, effector memory, and terminally differentiated T cells. 73,80 KLRB1 encodes CD161. 75 CD161 is associated with IL-17-secreting T cells [81] [82] [83] [84] [85] and is found in CB-naive T cells that become IL-17-secreting T cells. 82, 85 CD161 is also expressed in innate-like T cells including mucosalassociated invariant T cells, 86 gamma-delta-positive T cells, 82 invariant NK T cells, 87 and promyelocytic leukemia zinc-finger protein-positive CD4 + T cells. 88 In innate-like T cells, CD161 may track with infections or inflammatory complications of pregnancy. 88 GO analysis of DGEs in CD8 + T cells demonstrated enrichment for genes associated with T cell tolerance, proliferation, and the response to IFN-γ (Supplemental Fig. 3B ). In T cell Cluster 2, we found increased expression of IL-6-IL-17 axis genes including RORA, ARID5A, RBPJ, and IL6ST in cases compared to controls (Supplemental Fig. 3C ). The IL-6-IL-17 axis has been implicated in mediating the neurodevelopmental effects of maternal immune activation in mice [89] [90] [91] [92] and coordinating tissue inflammatory responses. 93 T cell antigen receptor (TCR) repertoire in T cells reflects selection by self and foreign antigens. To investigate the repertoire of TCRs in CB from SARS-CoV-2-exposed pregnancies and controls, we performed single-cell TCR sequencing. A total of 1943 T cells were analyzed, and T cells with TCR information were well equally distributed between subject and T cell populations (Supplemental Fig. 3D, E) . In CB from pregnancies complicated by maternal SARS-CoV-2 infection, there was significantly greater number of T cells with >5 clones compared to controls (40.4% in cases vs 30.9% in controls, Kolmogorov-Smirnov test p value 2.2e − 16) (Fig. 3d) . The T cell clonal expansion in CB from cases is consistent with results of T cell repertoire analysis from adults infected with SARS-CoV-2. 50 However, an analysis of T cell clonal diversity across four metrics (Shannon, inverse Simpson, Chao1, and abundancebased coverage estimator) revealed variability between samples in clonal diversity and no clear differences between cases and controls (Supplemental Fig. 3F) .

While the application of scRNAseq analysis to CBMC in the context of maternal illness is novel, our study is exploratory and has several limitations. Importantly, the small number of samples limits the generalizability of our conclusions. However, few studies have evaluated CB immune populations by single-cell transcriptomics, 94, 95 and our results illustrate an important and potentially underrecognized effect of the COVID-19 pandemic that should be further studied. Interpretation of our findings in light of prior known literature of CB leukocytes is challenging given the lack of detailed characterization of CBMC populations using scRNAseq as has been done in peripheral blood mononuclear cells from adults. 94, 95 Future studies using simultaneous epitope and transcriptome measurement in single cells will make these correlations possible. 96 All cases included in this study were classified as mild maternal SARS-CoV-2 infection; more severe maternal infection could result in more dramatic or different fetal immune genomic signatures, as it has been demonstrated that infection during pregnancy increases the possibility of intensive care unit admission. 5,97-99 Furthermore, the time from infection to delivery and CB collection likely affects the immune phenotype observed in CB. As the time of maternal infection and birth in our cohort varied between 7 and 66 days, more pronounced findings may be found in samples with more homogeneous timing between infection and collection. Lastly, all mothers affected with SARS-CoV-2 in our cohort had comorbidities including wellcontrolled thyroid dysfunction, obesity, or gestational diabetes. Although we included mothers with similar comorbidities in the control population (except for gestational diabetes) and all these comorbidities were medically managed, it is possible that our results are influenced by the comorbidities of the mothers. However, thyroid disease, obesity, or gestational diabetes in the mother have not been reported to trigger the transcriptional response patterns we observed in cases compared to controls. [100] [101] [102] [103] Lastly, we did not functionally test the CB leukocytes to evaluate the implications of the transcriptional observations we report.

Despite these limitations, the present study identifies transcriptional changes suggestive of a fetal immune response after maternal infection with SARS-CoV-2 in the absence of vertical transmission and suggests potential transplacental immune implications of maternal SARS-CoV-2 infection in the absence of mother-to-child transmission. The source of signals promoting transcriptional changes in neonatal monocytes and other immune cells in the absence of vertical transmission is unknown. Ex vivo studies have shown that transplacental transfer of IL-1β, IL-6, and tumor necrosis factor-α is limited. 104 Type I IFNs are increased in the peripheral circulation of patients with mild COVID-19, 32 but the ability of IFN to cross the human placenta is unclear. 105 Our results raise the possibility that pro-inflammatory signaling in the mother in response to SARS-CoV-2 might promote IFN signaling at the feto-maternal interface, and placental barrier dysfunction could result in loss of selective placental permeability to circulating maternal factors. Consistent with this possibility, a recent report has identified placental inflammatory responses in SARS-CoV-2-infected mothers compared to controls in the absence of viral mRNA in the placenta. 106 It is widely accepted that different viral infections, such as severe acute respiratory syndrome by influenza H1N1 virus or human immunodeficiency virus (HIV) in pregnant women, can compromise the offspring's fetal immune system even when the infection is limited to the mother or placental bed. 19, 21, [107] [108] [109] Furthermore, even without fetal infection, maternal HIV and H1N1 infections can be associated with neurological and behavioral diseases in offspring. 19, 21, [107] [108] [109] To our knowledge, this is the first time that scRNAseq has been used to explore the effect of maternal infection on the fetal immune system. Therefore, it is possible that the effects observed in our cohort are not SARS-COV-2-specific. However, S100A genes in monocytes were downregulated in our cohort in response to SARS-CoV2, while they are upregulated in monocytes of infants exposed to chorioamnionitis, 39 which might suggest, at minimum, differences in the response of fetal monocytes to maternal infection with SARS-CoV2 compared to chorioamnionitis.

Further experimental and functional data need to be collected to clarify how maternal infection with SARS-CoV-2 influences the fetal immune system, whether similar changes occur during other maternal infections without vertical transmission and the potential effect in the newborn as reported for HIV and influenza. 19, 21, 107, 108 Given the extensive literature linking maternal immune dysregulation and abnormal fetal development in viral infections, this study raises important questions about untoward effects of maternal SARS-CoV-2 on the fetus, even in the absence of vertical transmission, and highlights the need for further studies to better characterize the fetal immune response in pregnancies affected by SARS-CoV-2 infection.

The subjects were six infants born at term to mothers with or without SARS-CoV-2 infection in the third trimester. Parents of the infants provided informed consent before sample collection and study participation. The study was approved by the Institutional Review Board of the Mass General Brigham (IRB 2020P001478 and IRB2020P000804). Cord mononuclear cells were collected using Ficoll and cryopreserved as described. 29 We used dimethyl sulfoxide as our cryopreservant agent as it adequately conserves gene expression profiles in cryopreserved cells compared to fresh cells in droplet-based single-cell RNA-sequencing. 110 We excluded preterm infants, as a strong pro-inflammatory signature in CB has been reported in infants born preterm. 111 None of the infants was exposed to prenatal steroids, was diagnosed with intrauterine growth restriction, or had any neonatal morbidities. Placental viral load was measured as previously reported. 112 Single-cell RNA-sequencing CBMC aliquots were thawed in a 37°C water bath and resuspended in RPMI-1640 with 10% fetal bovine serum (FBS) (Thermo Fisher). Samples were centrifuged at 350 × g for 7 min at 4°C. Cells were resuspended in 100 μl of 1× phosphate-buffered saline with 2.5% FBS and 2 mM EDTA.

Dead cells and red blood cells were depleted using the EasySep Dead Cell Depletion Kit and EasySep RBC Depletion Reagent (STEMCELL), according to the manufacturer's instructions. Cells were resuspended in RPMI/10% FBS and counted. Cells were loaded onto the 10X Chromium controller at a targeted recovery density of 10,000 cells per sample. Samples were processed and sequencing libraries were created using the 10X the Chromium Next GEM single-cell V(D)J Reagent Kit v1.1 with human TCR V(D)J enrichment following the manufacturer's instructions.

Sequencing data were aligned to the genome and processed using the 10X Genomics Cell Ranger software, version 4.0.0. All cells were combined into a single dataset. Doublets were removed using Scrublet version 0.2.1, and the remaining cells were reclustered. Mitochondrial genes were filtered from the dataset. Cells with fewer than 250 or more than 2500 unique genes were excluded. Cells were then clustered using the Seurat R package (version 3.2.3). Specifically, the SCT functionality of Seurat was used to identify cell types that did not depend upon unique aspects of individual samples. Clustering resolution was set to 0.8, and the first 15 principal components were used. The data were log normalized and scaled to 10,000 transcripts per cell. The expression of known marker genes was used to assign each cluster to one of the main cell types. The Seurat FindMarkers function was used to identify genetic markers of cellular subtypes.

To identify differentially expressed genes by cell type, we performed a differential gene expression analysis using Monocle2. The analysis was conducted on each cell type and also certain unions of cell types with common traits. The data were modeled and normalized using a negative binomial distribution and counts data were normalized for gene length and read depth. Genes whose FDR was <5% were considered statistically significant. GO analysis was performed using gprofiler2 version 0.2.0, and terms were selected from the Biological Process category of GO terms.

TCR sequencing data were analyzed using the R package scRepertoire (version 3.12).

Sequencing data have been deposited in the Gene Expression Omnibus under accession no. GSE165193.

Clinical characteristics of pregnant women with Covid-19 in Wuhan

Pregnancy and postpartum outcomes in a universally tested population for SARS-CoV-2 in New York City: a prospective cohort study

Universal screening for SARS-CoV-2 in women admitted for delivery

Rates of maternal and perinatal mortality and vertical transmission in pregnancies complicated by severe acute respiratory syndrome coronavirus 2 (SARS-Co-V-2) infection: a systematic review

Characteristics of women of reproductive age with laboratoryconfirmed SARS-CoV-2 infection by pregnancy status-United States

Assessment of maternal and neonatal SARS-CoV-2 viral load, transplacental antibody transfer, and placental pathology in pregnancies during the COVID-19 pandemic

Synthesis and systematic review of reported neonatal SARS-CoV-2 infections

Vertical transmission of COVID-19: a systematic review and meta-analysis

Transplacental transmission of SARS-CoV-2 infection

Clinical characteristics and intrauterine vertical transmission potential of COVID-19 infection in nine pregnant women: a retrospective review of medical records

Analysis of SARS-CoV-2 vertical transmission during pregnancy

Possible vertical transmission of SARS-CoV-2 from an infected mother to her newborn

Risk factors associated with adverse fetal outcomes in pregnancies affected by coronavirus disease 2019 (COVID-19): a secondary analysis of the WAPM study on COVID-19

Outcomes of neonates born to mothers with severe acute respiratory syndrome coronavirus 2 infection at a large medical center

Maternal transmission of SARS-COV-2 to the neonate, and possible routes for such transmission: a systematic review and critical analysis

Characteristics and outcomes of 241 births to women with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection at five New York City Medical Centers

Characteristics and outcomes of pregnant women admitted to hospital with confirmed SARS-CoV-2 infection in UK: national population based cohort study

Long-term risk of neuropsychiatric disease after exposure to infection in utero

The fetal origins of mental illness

The immune system of HIV-exposed uninfected infants

Training the fetal immune system through maternal inflammation-a layered hygiene hypothesis

Association of maternal SARS-CoV-2 infection in pregnancy with neonatal outcomes

Compromised SARS-CoV-2-specific placental antibody transfer

Efficient maternal to neonatal transfer of antibodies against SARS-CoV-2 and BNT162b2 mRNA COVID-19 vaccine

The immunologic status of newborns born to SARS-CoV2-infected mothers in Wuhan

Antibodies in infants born to mothers with COVID-19 pneumonia

Bifidobacteria-mediated immune system imprinting early in life

Stereotypic immune system development in newborn children

Rapid establishment of a COVID-19 perinatal biorepository: early lessons from the first 100 women enrolled

Single-cell RNA-seq reveals new types of human blood dendritic cells, monocytes, and progenitors

Impaired type I interferon activity and inflammatory responses in severe COVID-19 patients

Severe COVID-19 is marked by a dysregulated myeloid cell compartment

Systems biological assessment of immunity to mild versus severe COVID-19 infection in humans

Distinct functional programming of human fetal and adult monocytes

Regulation of HLA class I and II expression by interferons and influenza A virus in human peripheral blood mononuclear cells

NF-κB signaling in macrophages: dynamics, crosstalk, and signal integration

Regulation of the innate immune system by autophagy: monocytes, macrophages, dendritic cells and antigen presentation

Expression of S100A alarmins in cord blood monocytes is highly associated with chorioamnionitis and fetal inflammation in preterm infants

Acute SARS-CoV-2 infection impairs dendritic cell and T cell responses

DC-ATLAS: a systems biology resource to dissect receptor specific signal transduction in dendritic cells

Effective chemokine secretion by dendritic cells and expansion of cross-presenting CD4−/CD8+ dendritic cells define a protective phenotype in the mouse model of Coxsackievirus myocarditis

Id1 suppresses anti-tumour immune responses and promotes tumour progression by impairing myeloid cell maturation

MafB, a target of microRNA-155, regulates dendritic cell maturation

Transcriptional programming of the dendritic cell network

Control of RelB during dendritic cell activation integrates canonical and noncanonical NF-κB pathways

TLR7-mediated activation of XBP1 correlates with the IFNα production in humans

Natural killer cell immunotypes related to COVID-19 disease severity

Heterogeneity of human bone marrow and blood natural killer cells defined by single-cell transcriptome

Single-cell landscape of immunological responses in patients with COVID-19

A single-cell atlas of the peripheral immune response in patients with severe COVID-19

Unique immunological profile in patients with COVID-19

Challenges and opportunities for consistent classification of human B cell and plasma cell populations

IL-10+ regulatory B cells are enriched in cord blood and may protect against cGVHD after cord blood transplantation

Characterization of the highly prevalent regulatory CD24hiCD38hi B-cell population in human cord blood

Human cord blood B cells differ from the adult counterpart by conserved ig repertoires and accelerated response dynamics

Systematic comparison of single-cell and single-nucleus RNAsequencing methods

Integrating single-cell transcriptomic data across different conditions, technologies, and species

A genome-wide transcriptomic analysis of protein-coding genes in human blood cells

Tonic signals: why do lymphocytes bother?

Endogenous Nur77 is a specific indicator of antigen receptor signaling in human T and B cells

Activated T cells induce rapid CD83 expression on B cells by engagement of CD40

The CD83 molecule-an important immune checkpoint

Antigen extraction and B cell activation enable identification of rare membrane antigen specific human B cells

The B cell antigen receptor controls AP-1 and NFAT activity through Ras-mediated activation of Ral

Role and regulation of pro-survival BCL-2 proteins in multiple myeloma

Single-cell mapping of progressive fetal-to-adult transition in human naive T cells

Flow cytometric characterization of human umbilical cord blood lymphocytes: immunophenotypic features

CD4 T cells with effector memory phenotype and function develop in the sterile environment of the fetus

Helios enhances the preferential differentiation of human fetal CD4+ naïve T cells into regulatory T cells

Fetal and adult hematopoietic stem cells give rise to distinct T cell lineages in humans

Regulation of T cell responses in the developing human fetus

CCR7 expression and memory T cell diversity in humans

Single-cell transcriptomics identifies an effectorness gradient shaping the response of CD4+ T cells to cytokines

CD161-expressing human T cells

Impaired immune cell cytotoxicity in severe COVID-19 is IL-6 dependent

Molecular signature of CD8+ T cell exhaustion during chronic viral infection

Granzyme A from cytotoxic lymphocytes cleaves GSDMB to trigger pyroptosis in target cells

The multifaceted output of c-Jun biological activity: focus at the junction of CD8 T cell activation and exhaustion

Two subsets of memory T lymphocytes with distinct homing potentials and effector functions

Pro-inflammatory human Th17 cells selectively express P-glycoprotein and are refractory to glucocorticoids

CD161 is a marker of all human IL-17-producing T-cell subsets and is induced by RORC

Analysis of CD161 expression on human CD8+ T cells defines a distinct functional subset with tissue-homing properties

Cutting edge: phenotypic characterization and differentiation of human CD8+ T cells producing IL-17

Human interleukin 17-producing cells originate from a CD161+CD4+ T cell precursor

Human MAIT and CD8αα cells develop from a pool of type-17 precommitted CD8+ T cells

Innate-like effector differentiation of human invariant NKT cells driven by IL-7

CD161 contributes to prenatal immune suppression of IFNγ-producing PLZF+ T cells

Contributions of maternal and fetal antiviral immunity in congenital disease

The maternal interleukin-17a pathway in mice promotes autismlike phenotypes in offspring

Reversing behavioural abnormalities in mice exposed to maternal inflammation

Maternal immune activation in mice disrupts proteostasis in the fetal brain

The IL-17 family of cytokines in health and disease

Single-cell transcriptomic landscape of nucleated cells in umbilical cord blood

Characterization of dendritic cell subtypes in human cord blood by single-cell sequencing

Simultaneous epitope and transcriptome measurement in single cells

Update: characteristics of symptomatic women of reproductive age with laboratory-confirmed SARS-CoV-2 infection by pregnancy status-United States

Higher severe acute respiratory syndrome coronavirus 2 infection rate in pregnant patients

Severe acute respiratory syndrome coronavirus 2 serology levels in pregnant women and their neonates

Maternal influences on the transmission of leukocyte gene expression profiles in population samples from Brisbane

Maternal autoimmune thyroid disease and the fetal immune system

Maternal obesity alters immune cell frequencies and responses in umbilical cord blood samples

Diabetic pregnancy activates the innate immune response through TLR5 or TLR1/2 on neonatal monocyte

Transfer of proinflammatory cytokines across term

Interferons and proinflammatory cytokines in pregnancy and fetal development

Maternal respiratory SARS-CoV-2 infection in pregnancy is associated with a robust inflammatory response at the maternal-fetal interface

Maternal immune responses and obstetrical outcomes of pregnant women with COVID-19 and possible health risks of offspring

The impact of maternal HIV infection on cord blood lymphocyte subsets and cytokine profile in exposed non-infected newborns

Maternal HIV status skews transcriptomic response in infant cord blood monocytes exposed to Bacillus Calmette-Guerín

DMSO cryopreservation is the method of choice to preserve cells for droplet-based single-cell RNA sequencing

Differential patterns of 27 cord blood immune biomarkers across gestational age

SARS-CoV-2 viral load is associated with increased disease severity and mortality

We thank Dr. Pablo J. Patino (Department of Immunology, School of Medicine, Universidad de Antioquia) for helpful discussions. 

Parents of the infants provided informed consent before sample collection and study participation. The study was approved by the Institutional Review Board of the Mass General Brigham (IRB 2020P001478 and IRB2020P000804).

The authors declare no competing interests.

Supplementary information The online version contains supplementary material available at https://doi.org/10.1038/s41390-021-01793-z.Correspondence and requests for materials should be addressed to Juan D. Matute or Brian T. Kalish.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.