1. Study design
This is a case-control, cross sectional, single-centre study, approved by the local Ethical Committee (Comitato Etico dell’Area Vasta Emilia Nord, protocol number 177/2020, March 11th, 2020) and by the University Hospital Committee (Direzione Sanitaria dell’Azienda Ospedaliero-Universitaria di Modena, protocol number 7531, March 11th, 2020). Each participant, including healthy controls, provided informed consent according to Helsinki Declaration, and all uses of human material have been approved by the same Committees. A total of 14 pregnant women with SARS-CoV2 infection was included in the study; they had a median age of 33.8 years (range 19-39). Patients were matched for age and gender with 28 pregnant women negative to nasopharyngeal swab (median 33.9 years, range 18-42) and a total of 15 non pregnant healthy women (CTR), median age 39 years (range 25-50 years). We recorded demographic data, medical history, symptoms, signs, temperature, and main laboratory findings from each patient. For details, see supplementary Table 1.
Pregnant women were eligible for inclusion if they were aged 18 years or older, able to provide informed consent and diagnosed with SARS- CoV-2 infection. Confirmed SARS-CoV-2 infection was defined as nasopharyngeal swab reverse transcription–polymerase chain reaction (RT-PCR) test results positive for SARS-CoV-2. According to routine methods, we could measure anti-SARS-CoV-2 IgM and IgG in 12 women positive to the swab, and 9 had detectable plasma levels of IgM and IgG, while three were negative to both antibodies. All the 15 control women that we could test were negative. In all cases, paucisymptomatic symptoms were fever <37.5°C, cough, loss of smell or taste, and no drugs were administered to treat them.
The total number and type of leukocytes in peripheral blood was analyzed by hemocytometer in the Clinical Laboratory of the University Hospital, that also analyze all the biochemical parameters quoted in the manuscript, according to routine methods.
2. Blood collection
Blood samples (up to 20 mL) were obtained after informed consent. In some donors, blood was obtained after diagnosis of SARS-CoV-2 infection. Peripheral blood mononuclear cells (PBMC) were isolated according to standard procedures and stored in liquid nitrogen until use. Plasma was collected and stored at -80°C until use. Measurements were taken from individual patients; in the case of plasma, each measurement was performed in duplicate and only the mean was considered and shown.
3. Mass cytometry analysis
Thawed PBMC were washed twice with PBS and stained with Maxpar® Direct™ Immune Profiling Assay™ (Fluidigm), a dry 30-marker antibody panel (viability marker Cell-ID™ Intercalator-103Rh included) plus the addition of 6 drop-in catalog antibodies (Fluidigm) and 2 custom-conjugated mAbs, for a total of 38 markers. The markers were the following: CD3, CD19, CD45, CD4, CD20, CD45RA, CD8, CD25, CD45RO, CD11c, CD27, CD56, CD14, CD28, CD57, CD16, CD38, CD66b, CCR7, CXCR3, CXCR5, HLA-DR,IgD, TCRγδ, CD123, CD127, CD161, CD294, CCR4, CCR6, CXCR1, PDL1, CD80, CD40, CD24, PD1-1, CD11b/MAC, CD21, IgM. See Supplementary Table 5 for the complete list of mAbs used. At least 300,000 events were acquired per sample. Data in FCS file format were normalized for intra-file and inter-file signal drift using the FCS Processing tab in the CyTOF Software 6.7. The method is a two-step algorithm that first identifies the EQ Four Element Calibration Beads and then applies the dual count values registered by the beads to calculate the normalization factor to be applied to the data.
4. Representation of high parameter cytometry
Compensated and normalized Flow Cytometry Standard (FCS) 3.0 files were imported into FlowJo software version 10 (Becton Dickinson, San Josè, CA) and pre-processed excluding EQ Four Element Calibration Beads and doublet using Gaussian Discrimination parameters. Then, were selected live undamaged CD45+ and excluded artifact cells (CD3+CD19+ or CD3+CD14+). All living CD45+ were exported for further analysis in R using Bioconductor libraries CATALYST (version 1.12.2) 21 and diffcyt (version 1.8.8) 22. The data were transformed using arcsinh with cofactor = 5 to make the distributions more symmetric and to map them to a comparable range of expression. The main cell population identification was performed through unsupervised clustering using the FlowSOM algorithm (K= 30). 2D visual representation was performed applying Uniform Manifold Approximation and Projection (UMAP). Then, the clusters identified as CD4+, CD8+, or CD19+ lymphocytes, were selected and re-clusterized separately to describe more in-depth the cellular distribution of each sub-population. We used K= 15 for CD4+ and CD8+ T cells, while K= 9 for CD19+ cells. Clusters with similar markers distribution were merged. Then we re-applied UMAP for dimensionality reduction and visualization purposes. Statistical analysis was performed using generalized linear mixed models (GLMM) applying as FDR cutoff = 0.05.
5. Polychromatic flow cytometry
5.1. T cell characterization
For the analysis of T cell skewing toward Th1, Th2, or Th17, and of chemokine receptor expression, thawed PBMC were washed twice with PBS and stained with the viability marker AQUA LIVE DEAD (ThermoFisher). Then, up to 1 million cells were washed and stained at 37 °C with the following mAbs: anti-CXCR3- AF488, -CXCR4-PE. Cells were washed again and stained at room temperature with anti-CD161-PC7, -CCR6-BV605, -CCR4-PE-CF594, -CD4-AF700, -CD8- APC-Cy7. Cells were washed, fixed, and permeabilized using Foxp3/Transcription Factor Staining Buffer Set (ThermoFisher). Finally, cells were stained with anti- GATA3-BV421 and anti-TBET-APC and washed. A minimum of 500,000 PBMC were acquired by using Attune NxT acoustic focusing flow cytometer (ThermoFisher). mAbs used are listed in Supplementary Table 6.
5.2. Proliferation assay
Cells were stimulated for 6 days in resting conditions, or after stimulation with anti-CD3 plus anti-CD28 mAbs (1 μg/mL each, Miltenyi Biotech, Bergisch Gladbach, Germany) and with 20 ng/mL IL-2. The fluorescent dye 5,6-carboxyfluorescein diacetate succinimidyl ester (CSFE) was used at a concentration of 1 μg/mL (ThermoFisher) according to standard procedures 23. Flow cytometric analyses for the identification of cycling cells belonging to different T cell populations were performed by gating CD4+, CD8+ T cells and CD19+ B cells. mAbs used are listed in Supplementary Table 2.
5.3 In vitro stimulation and intracellular cytokine staining
For functional assays on cytokine production by T cells, thawed isolated PBMCs were stimulated for 16 h at 37 °C in a 5% CO2 atmosphere with anti-CD3/CD28 (1 μg/mL) in complete culture medium (RPMI 1640 supplemented with 10% fetal bovine serum and 1% each of l- glutamine, sodium pyruvate, nonessential amino acids, antibiotics, 0.1 M HEPES, 55 μM β-mercaptoethanol). For each sample, at least 2 million cells were left unstimulated as negative control, and 2 million cells were stimulated. All samples were incubated with a protein transport inhibitor containing brefeldin A (Golgi Plug, Becton Dickinson) and previously titrated concentration of CD107a-PE. After stimulation, cells were stained with LIVE-DEAD Aqua (ThermoFisher Scientific) and surface mAbs recognizing CD4 AF700, and CD8 APC- Cy7 (Biolegend, San Diego, CA, USA). Cells were washed with stain buffer and fixed and permeabilized with the cytofix/cytoperm buffer set (Becton Dickinson) for cytokine detection. Cells were next stained with previously titrated mAbs recognizing CD3 PE- Cy5, IL-17 BV421, TNF BV605, IFN-γ FITC, IL-4 APC, or granzyme-B BV421 (all mAbs from Biolegend). Then, a minimum of 100,000 cells per sample were acquired on a Attune NxT acoustic cytometer (ThermoFisher) 24. mAbs used are listed in Supplementary Table 2.
5. Quantification of cytokine plasma levels
The plasma levels of 62 molecular species were quantified using a Luminex platform (Human Cytokine Discovery, R&D System, Minneapolis, MN) for the simultaneous detection of the following molecules: G-CSF, PDGF-AA, EGF, PDGF-AB/BB, VEGF, GM-CSF, FGF, GRZB, IL-1A, IL-1RA, IL-2, IL-27, IL-4, IL-6, IL-10, IL-13, TNF, IL-17C, IL-11, IL-18, IL-23, IL-6RA, IL-19, IFN-B, IL-3, IL-5, IL-7, IL-12p70, IL-15, IL-33, TGF-B, IFN-G, IL-1B, IL-17, IL-17E, CCL3, CCL11, CCL20, CXCCL1, CXCL2, CCL5, CCL2, CCL4, CCL19, CXCL1, CXCL10, PD-L1, FLT-3, TACI, FAS, LEPTIN R, APRIL, OPN, BAFF, LEPTIN, BMP4, CD40 LIGAND, FAS LIGAND, BMP7, BMP2, TRAIL, according to the manufacturer’s instruction. Data in the scatter plots represent the mean of two technical replicates.
6. Analysis of the correlations among all parameters
To identify possible correlations among the parameters we have studied, we have designed a table containing: i) all 27 clusters percentages obtained using the unsupervised analysis on living CD45+ cells; ii) 12 biochemical parameters reported in Supplementary Table 1, and iii) 61 plasma cytokines out of 62 because all values of IL-17p70 were identical. Pairwise correlations between variables were calculated and visualized as a correlogram using R function corrplot (Figure 5, and Supplementary Tables 2-4). Spearman’s rank correlation coefficient (ρ) was indicated by color scale; significance was indicated by *P < 0.05, **P < 0.01, and ***P < 0.001. All variables were displayed using original order without applying any hierarchical clustering.
7. Statistical analysis
High-dimensional cytometric analysis was performed by using differential discovery in high-dimensional cytometry via high-resolution clustering. Quantitative variables were compared using Mann-Whitney test. Statistical analyses were performed using Prism 6.0 (GraphPad Software Inc., La Jolla, USA).