Next Article in Journal
Therapeutic Effects of Aloe saponaria against Ulcerative Colitis Induced by Dextran Sulfate Sodium
Next Article in Special Issue
Identifying the Carcinogenic Mechanism of Malignant Struma Ovarii Using Whole-Exome Sequencing and DNA Methylation Analysis
Previous Article in Journal
Multiple Sclerosis: Inflammatory and Neuroglial Aspects
Previous Article in Special Issue
Human Blood Extracellular Vesicles Activate Transcription of NF-kB-Dependent Genes in A549 Lung Adenocarcinoma Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Proof-of-Concept Analysis of B Cell Receptor Repertoire in COVID-19 Patients Undergoing ECMO by Single-Cell V(D)J and Gene Expression Sequencing

1
Department of Research, IRCCS ISMETT (Istituto Mediterraneo per i Trapianti e Terapie ad Alta Specializzazione), 90127 Palermo, Italy
2
Fondazione Ri.MED, Department of Research IRCCS ISMETT, via Ernesto Tricomi 5, 90145 Palermo, Italy
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Curr. Issues Mol. Biol. 2023, 45(2), 1471-1482; https://doi.org/10.3390/cimb45020095
Submission received: 17 January 2023 / Revised: 6 February 2023 / Accepted: 8 February 2023 / Published: 9 February 2023
(This article belongs to the Special Issue Next-Generation Sequencing (NGS) Technique and Personalized Medicine)

Abstract

:
SARS-CoV-2, which causes COVID-19, has altered human activities all over the world and has become a global hazard to public health. Despite considerable advancements in pandemic containment techniques, in which vaccination played a key role, COVID-19 remains a global threat, particularly for frail patients and unvaccinated individuals, who may be more susceptible to developing ARDS. Several studies reported that patients with COVID-19-related ARDS who were treated with ECMO had a similar survival rate to those with COVID-19-unrelated ARDS. In order to shed light on the potential mechanisms underlying the COVID-19 infection, we conducted this proof-of-concept study using single-cell V(D)J and gene expression sequencing of B cells to examine the dynamic changes in the transcriptomic BCR repertoire present in patients with COVID-19 at various stages. We compared a recovered and a deceased COVID-19 patient supported by ECMO with one COVID-19-recovered patient who did not receive ECMO treatment and one healthy subject who had never been infected previously. Our analysis revealed a downregulation of FXYD, HLA-DRB1, and RPS20 in memory B cells; MTATP8 and HLA-DQA1 in naïve cells; RPS4Y1 in activated B cells; and IGHV3-73 in plasma cells in COVID-19 patients. We further described an increased ratio of IgA + IgG to IgD + IgM, suggestive of an intensive memory antibody response, in the COVID ECMO D patient. Finally, we assessed a V(D)J rearrangement of heavy chain IgHV3, IGHJ4, and IGHD3/IGHD2 families in COVID-19 patients regardless of the severity of the disease.

1. Introduction

Since December 2019, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has posed a hazard to global public health. The coronavirus disease 2019 (COVID-19) shows a wide range of clinical manifestations, ranging from asymptomatic presentation to critical illness with severe pneumonia, acute respiratory distress syndrome (ARDS), or multiple organ failure [1]. The similarities between the worst SARS-CoV-2 consequences and seasonal influenza problems, such as ARDS or multiple organ failure, have suggested a role for extracorporeal membrane oxygenation (ECMO) implantation in patients with the most severe pulmonary decompensation [2,3]. It has recently been established that ECMO can be utilized as a rescue therapy due to the temporary replacement of lung and/or cardiac function [4]. Several studies have reported that patients with COVID-19-related ARDS who were treated with ECMO showed a survival rate comparable with those with COVID-19-unrelated ARDS [5,6,7]. Clearing the SARS-CoV-2 infection and hence influencing patients’ clinical outcomes is also mediated by humoral and adaptive immune responses [8]. A crucial role is played by the B cell antigen receptor (BCR) responsible for the recognition of pathogens. The processes of recombination and assembly of the variable and constant regions of the V, D, and J segments are crucial in generating an immense repertoire of antibodies responsible for the recognition of diverse pathogens [9]. Briefly, the antigen-binding domain of immunoglobulins is composed of two polypeptide chains. The exons that encode the antigen-binding domains are assembled from V (variable), D (diversity), and J (joining) gene segments by a process defined as “cut-and-paste” DNA rearrangements. This process, named V(D)J recombination, selects a pair of segments, introduces double-strand breaks next to each segment, deletes the intermediate DNA, and ligates the segments together. Rearrangements take place in a well-ordered way, with D-to-J joining proceeding before a V segment is joined to the rearranged D–J segments (Figure 1) [10]. Meeting a pathogen then unleashes rapid hypermutation (SHM) and class-switch recombination (CSR), thereby increasing antigen binding [11].
In this study, we used a single-cell approach to look at the dynamic changes in the transcriptomic BCR repertoire in patients with COVID-19 at various stages, comparing a recovered and a deceased COVID-19 patient who had been supported with ECMO with one COVID-19-recovered patient who had not received ECMO treatment and one healthy subject who had never previously been infected nor tested positive for SARS-CoV-2 antibodies.

2. Materials and Methods

2.1. Samples

Four subjects were enrolled in this study at IRCCS ISMETT: one COVID-19 patient on ECMO (survived), one COVID-19 patient on ECMO (died), one COVID-19 patient without ECMO (recovered), and one healthy control who tested negative for both a nasopharyngeal swab (NPS) and anti-Spike and anti-N IgG/IgM. The study was approved by the IRCCS ISMETT Institutional Research Review Board (IRRB 00/21) and the Ethics Committee of IRCCS ISMETT, with all enrolled individuals signing a written informed consent form. Fresh blood samples were acquired from patients on day three of their hospitalization.

2.2. Preparation of Single-Cell Suspensions

Venous blood was collected in K3EDTA tubes (Greiner Bio-One GmbH, Kremsmunster, Austria). Peripheral blood mononuclear cells (PBMCs) were isolated by density gradient centrifugation on Lympholyte Cell Separation Media (Cedarlane Laboratories Limited, Burlington, ON, Canada). Afterwards, CD19+ B cells were separated from PBMCs by immuno-magnetic sorting using anti-CD19 magnetic microbeads (REAlease CD19 MicroBeads Kit, Miltenyi Biotec, Auburn, CA, USA). The CD19+ B cells obtained from immuno-magnetic sorting displayed a purity yield greater than 98%, which was determined by flow cytometry analysis.

2.3. Droplet-Based Single-Cell Sequencing

Using the single-cell 3′ Library and Gel Bead Kit V3.1 (10X Genomics, 1,000,121) and the Chromium Single Cell G Chip Kit (10X Genomics, 1,000,120), the cell suspension was loaded onto the Chromium Single Cell Controller (10X Genomics) to generate single-cell gel beads in an emulsion according to the manufacturer’s protocol. In short, single cells were suspended in phosphate-buffered saline containing 0.04% bovine serum albumin. Approximately 1200 cells/µL were added to each channel, and the target number of cells to be recovered was estimated to be approximately 6000. Captured cells were lysed, with the released RNA barcoded through reverse transcription in individual GEMs. Reverse transcription was performed on a Veriti 96 Well Thermal Cycler (ThermoFisher) at 53 °C for 45 min, followed by 85 °C for 5 min, and then held constant at 4 °C. The generated cDNA was then amplified, with the quality assessed using the 4200 TapeStation System (Agilent). According to the manufacturer’s instructions, scRNA-seq libraries were constructed using the Single Cell 3′ Library and Gel Bead Kit V3.1. Finally, the libraries were sequenced using an Illumina NextSeq500/550 High Output Reagent Cartridge v2 300 cycle sequencer, with a sequencing depth of at least 100,000 reads per cell using the paired-end strategy.

2.4. Single-Cell RNA Sequencing

First, Cell Ranger v 6.0.0 (10X Genomics) was used to demultiplex the cellular bar-codes and align the reads to the human transcriptome (human reference version GRCh38) for each sample [12]. Second, each output, which was a raw unique molecular identifier (UMI) count matrix, was transformed into a Loupe Browser object. We filtered out genes that were expressed in less than five cells by using the Loupe Browser v 6.0.0 function. Several criteria were then applied to each dataset to remove cells of low quality: cells with fewer than 200 genes or more than 6000 genes. Data processing and analysis were performed with Prism GraphPad V5.0d software (GraphPad Software, San Diego, CA, USA).

2.5. Single-Cell V(D)J Sequencing and Data Processing

Single-cell V(D)J sequencing was performed following the protocol provided by the 10X Genomics Chromium Single Cell Immune Profiling Solution. The analysis pipelines in Cell Ranger (10X Genomics, version 6.0.0) were used for single-cell sequencing data processing and were loaded in the Loupe V(D)J Browser v 4.0.0. V(D)J sequence assembly and paired clonotype calling were performed using cellranger vdj with reference = refdata-cellranger-vdj-GRCh38-alts-ensembl-6.0 for each sample. Data analysis was performed using the tidyverse package v1.3.0.

3. Results

3.1. Study Design and Profiling of B Cells

In order to evaluate the different signatures of B cell receptors in the B lymphocytes of COVID-19 patients that needed extracorporeal membrane oxygenation (ECMO) in response to acute respiratory distress syndrome (ARDS), we performed scRNA-seq and single-cell BCR sequencing on the CD19+ immune cells of a COVID-19 patient under ECMO who later survived the treatment (COVID ECMO S), a COVID-19 patient under ECMO who later died (COVID ECMO D), one COVID-19-recovered patient who did not receive ECMO support (COVID R), and one healthy control who tested negative for both a nasopharyngeal swab (NPS) and anti-Spike and anti-N IgG/IgM. In Table 1, we report the characteristics of the patients included in this study. After filtering, 7884 B cells were obtained for the single-cell transcriptome data and 7733 B cells were obtained for V(D)J analysis. The scRNA-seq and single-cell paired BCR analysis were then combined, with the study limited to 5830 B lymphocytes having full-length productive paired IGH-IGK/IGL.

3.2. Features of B Cell Subsets

We started our analysis by assigning cell identities based on B cell surface marker expression as indicated below. According to the average log fold change of the canonical markers, we identified four clusters: (1) memory B cells; (2) naïve cells; (3) activated B cells, and (4) plasma cells (Figure 2A). Memory B cell subsets were identified by the presence of the typical memory marker CD27, while naïve B cells expressed the heavy chain IgD immunoglobulin (IGHD) [13]. Activated B cells expressed high levels of CD79, a typical marker of B cell activation [14], while plasma cells, together with CD79, expressed a high level of X-box binding protein 1 (XBP1), a transcriptional regulator critical for supporting the cellular reprograming activities during B-to-plasma-cell transition, which permits antibody release during terminal differentiation [15]. All B cell subpopulations analyzed expressed the membrane-spanning 4A1 (MS4A1 or CD20) marker, which unequivocally identifies B cells. In order to assess differences in the B cell populations present in the four samples, we investigated the relative frequencies of the four clusters in each sample (Figure 2B). In general, naïve cells and activated B cells accounted for the largest portion of B cells in each sample. Among COVID-19 patients, it is notable that, in the COVID ECMO D patient, the naïve cell to activated B cell ratio was different from the other samples with a shift towards the activated B cell population (Figure 2B,C). We also investigated the transcriptional signatures of the COVID-19 patients’ B cell populations (Figure 2D).

3.3. Transcriptional Signatures of Different B Cell Subpopulations of COVID-19 Patients

Once the clusters present in the B cell populations of the COVID ECMO S, COVID ECMO D, COVID R, and healthy control were established, we investigated the transcriptional profiles of the clusters in depth to find clue genes and pathways potentially involved in the worsening of clinical outcomes (Figure 3A).
Among the deregulated genes in the naïve B cells, we found a downregulation of the major histocompatibility complex, Class II, DQ Alpha 1 (HLA-DQA1), and the mitochondrial gene MT-ATP8 (Mitochondrially Encoded ATP Synthase Membrane Subunit 8) in COVID-19-infected subjects compared to the uninfected subject (Figure 3B). Interestingly, we discovered that MT-ATP8 expression levels were associated with disease severity; in particular, the two COVID-19 patients receiving ECMO had near-zero MT-ATP8 expression compared to the no-ECMO patient and the healthy control. The most notable downregulated gene in activated B cells was the ribosomal protein S4, Y-Linked 1 (RPS4Y1), which has been linked to HLA-DQA1 and has been implicated in resistance to treatment with corticosteroids and cyclosporin A (Figure 3B). The memory B cell cluster showed downregulation of three genes: FXYD, Domain-Containing Ion Transport Regulator 5; HLA-DRB1, a beta chain of antigen-presenting major histocompatibility complex class II (MHCII) molecule; and RPS20, ribosomal protein S20 (Figure 3B). The profiling of plasma cell clusters revealed a considerable downregulation of an immunoglobulin heavy chain variable region, IGHV3-73, which has been described as permitting antigen-binding activity and immunoglobulin receptor-binding activity (Figure 3B).

3.4. IgH Class Switching of B Cells and CDR3 Length and Specific Rearrangements of V(D)J Genes

To delineate the dynamic changes in IgH class switching, which is crucial for a comprehensive analysis of the BCR repertoire, we evaluated the distribution of Immunoglobulins, including IgA, IgD, IgG, and IgM in the COVID ECMO S, COVID ECMO D, COVID R, and healthy control (Figure 4A). In the healthy control, 61.3% of B cells expressed the IgM isotype, followed by IgG (22.4%), IgA (16.2%), and IgD (0.1%). Of COVID ECMO D patient B cells, 55.8% expressed the IgM isotype, followed by IgG (28.6%), IgA (18.2%), and IgD (2.3%). Of COVID ECMO S patient B cells, 70.2% expressed the IgM isotype, followed by IgG (18.5%), IgA (11.3%), and IgD (0.1%), while 65.2% of COVID R patient B cells expressed the IgM isotype, followed by IgG (13.9%), IgA (18%), and IgD (2.8%) (Figure 4B). Interestingly, compared to the other subjects involved in the study, the ratio of (IgA + IgG) to (IgD + IgM) increased in the COVID ECMO D patient, suggesting an intensive antibody response.
The distribution of clonally increased B cells in the COVID ECMO S, COVID ECMO D, COVID R, and healthy control was studied. The percentage of clonally expanded cells in COVID ECMO D patient B cells was 1.97%, which was somewhat higher than in the other three subjects, specifically: 0.56% in COVID ECMO S, 1.0206% in COVID R, and 1.14% in the healthy control (Figure 4C).
The CDR3 length of the BCR heavy chain in the COVID ECMO S, COVID ECMO D, and COVID R patients, and in the healthy control, ranged from 6 to 51 amino acids (aa), with an average of 14 aa for each sample (Figure 4D), and no significant difference between the samples.
In order to investigate the V(D)J rearrangements of the BCR heavy chain in our four samples, we detected the differences in usage frequency of the V, D, and J gene segments. A total of 53 IGHV gene segments, 27 IGHD gene segments, and 6 IGHJ gene segments were identified from all of the B cells. We generated a distribution histogram of IGHV, IGHD, and IGHJ gene usage frequency for the total number of B cells in our samples. We assessed that IGHJ4 was most frequently used in all samples, followed by IGHJ6, IGHJ5, and IGHJ3, but with no significant differences related to the disease status. IGHJ1 and IGHJ2 showed the lowest utilization (Figure S1).
The analysis of selective usage of IGHD and IGHV genes is summarized in Figures S2 and S3. In our samples, we detected an over-representation of the IGHD3 and IGHD2 families. Among the IGHD3 genes, the most expressed is IGHD3-22, which was highly expressed in both COVID-19 patients receiving ECMO (22.94% COVID ECMO D; 17.51% COVID ECMO S), compared to COVID R (15.45%) and the healthy control (13.60%), while the least expressed is IGHD3-3, in both COVID patients that underwent ECMO (8.72% COVID ECMO D; 10.31% COVID ECMO S), compared to COVID R (17.89%) and the healthy control (11.48%) (Figure S2).
Regarding the IGHV genes, we detected that the most represented is the IGHV3 family, especially IGHV3-23, which was expressed less in both COVID patients that underwent ECMO (8.85% COVID ECMO D; 10.14% COVID ECMO S), compared to COVID R (13.70%) and the healthy control (13.14%), and IGHV3-33, with no significant differences related to the disease status (Figure S3).
A total of 318 unique V–J combinations were identified in all of the whole B cells (Figure 4E). The top paired V–J frequency in the COVID patients receiving ECMO was IGHV3-33/IGHJ4 (4.68% COVID ECMO D; 5.33% COVID ECMO S), compared to COVID R (4.02%) and the healthy control (4.22%).

4. Discussion

In the present study, we aimed to perform a comprehensive analysis of the BCR repertoire, which is the genetic source of neutralizing antibodies (nAbs) [16], in COVID-19 patients that underwent ECMO using a single-cell approach. Since the quantification of gene expression is necessary to uncover disease etiology and progression [17], our goal was to shed light on the complex mechanisms that drive the immune response to the virus stimulus. The cohort included a COVID-19 patient who died after receiving ECMO (COVID ECMO D), a COVID-19 patient who survived the treatment (COVID ECMO S), a COVID-19-recovered patient without ECMO support (COVID R), and one healthy control who tested negative for both a nasopharyngeal swab (NPS) and anti-Spike and anti-N IgG/IgM. The first step was to identify four clusters—memory B cells, naïve cells, activated B cells, and plasma cells—on CD19+ immune B cells based on the average log fold change of the canonical markers, and we assessed the relative frequencies of the four clusters in each sample. We found that naïve cells and activated B cells accounted for the largest portion of B cells in each sample but, in the COVID ECMO D patient, the naïve cell to activated B cell ratio was different from the other samples, with a shift towards activated B cell populations.
The transcriptomic analysis highlights different signatures of the B cell subpopulations among the subjects studied. In particular, the naïve cells of COVID-19-infected subjects showed downregulation of MT-ATP8 and HLA-DQA1. The first is a mitochondrial gene that encodes for an ATP synthase responsible for changing cellular energy capacity, increasing mitochondrial oxidative stress and/or modulating apoptosis [18]. Interestingly, we observed that its levels of expression were related to the severity of the disease. Indeed, the two COVID-19 patients receiving ECMO showed lower MT-ATP8 expression, close to zero, compared to the no-ECMO patient and the healthy control. The second is known to bind peptides derived from antigens that access the endocytic route of antigen-presenting cells (APC) and present them on the cell surface for recognition by the CD4 T-cells. In activated B cells, the ribosomal protein S4, Y-Linked 1 (RPS4Y1), described as being involved, together with HLA-DQA1, in a resistance to the treatment with corticosteroids combined with cyclosporin A, commonly used in autoimmune diseases, auto-inflammatory diseases, and transplant rejection [14], was significantly downregulated. The memory B cell cluster showed downregulation of three genes: (1) FXYD, Domain-Containing Ion Transport Regulator 5, a glycoprotein that functions in the upregulation of chemokine production and is involved in the reduction of cell adhesion via its ability to downregulate E-cadherin; (2) HLA-DRB1, a beta chain of the antigen-presenting major histocompatibility complex class II (MHCII) molecule, known to guide antigen-specific T-helper effector functions, both in antibody-mediated immune responses and macrophage activation, to ultimately eliminate infectious agents and transformed cells [19,20,21,22,23,24,25]; and (3) RPS20, ribosomal protein S20, for which downregulation contributes to a stress phenotype by suppressing genomic and cellular stability [26]. Plasma cell cluster profiling showed significant downregulation of an immunoglobulin heavy chain variable region, IGHV3-73, described as enabling antigen-binding activity and immunoglobulin receptor-binding activity and predicted to be involved in several processes, including activation of the immune response and the defense response to other organisms, as well as phagocytosis.
We also evaluated the clonality of B cells in our cohort and found no significant differences, apart from a slight difference in clonally expanded cells in COVID ECMO D patient B cells (1.97%), versus the B cells of the other three subjects studied: 0.56% in COVID ECMO S, 1.0% in COVID R, and 1.14% in the healthy control. We also evaluated the B cell immunoglobulin distribution (IgA, IgD, IgG, and IgM) and found that the ratio of IgA + IgG to IgD + IgM was greater in the COVID ECMO D patient, suggesting an intensive memory antibody response. Lastly, we assessed the V(D)J rearrangements of the BCR heavy chain of our four samples and determined that the IGHV3, IGHJ4, and IGHD3/IGHD2 families were the most frequently used in all samples that we analyzed. These results are consistent with previously published data, in which the authors also reported the highest frequency of the IGHV3/IGHJ4 pair in symptomatic patients [27]. In addition, among the IGHD3 genes, the most expressed was IGHD3-22, which was highly expressed in both COVID-19 patients that underwent ECMO (22.94% COVID ECMO D; 17.51% COVID ECMO S), compared to COVID R (15.45%) and the healthy control (13.60%). In a recent study [28], the authors identified more than 100 anti-SARS-CoV-2 antibodies containing a conserved YYDRxG motif exclusively encoded by the IGHD3-22 gene. The presence of this motif is strictly related to a high neutralizing activity against different SARS-CoV-2 variants, as well as other SARS-related coronaviruses [28]. Conversely, IGHD3-3 is expressed less in both COVID-19 patients that underwent ECMO (8.72% COVID ECMO D; 10.31% COVID ECMO S), compared to COVID R (17.89%) and the healthy control (11.48%). Finally, a total of 318 unique V–J combinations were identified from all the B cells. The top paired V–J frequency in the COVID-19 patients, COVID ECMO D and COVID ECMO S, were 4.68% and 5.33%, respectively, compared to COVID R (4.02%) and the healthy control (4.22%). Our study, through the integration of transcriptomic data and single-cell paired BCR profiles, revealed BCR repertoire changes in COVID-19 patients with severe illness, which is consistent with other studies focusing on subjects infected and/or vaccinated against SARS-CoV-2 subjects [29].

5. Conclusions

The present study aimed to perform a comparative analysis of the BCR repertoire of COVID-19 patients with different clinical profiles and outcomes: a COVID-19 patient receiving ECMO who later died, a COVID-19 patient receiving ECMO who later survived treatment, a COVID-19 patient who recovered without ECMO support, and one healthy control. After the identification of four clusters on CD19+ immune B cells—memory B cells, naïve cells, activated B cells, and plasma cells—we evaluated the BCR repertoires from different points of view. (1) Firstly, we assessed the transcriptome highlighting different signatures of the B cell subpopulations in COVID-19-infected subjects showing a downregulation of FXYD, HLA-DRB1, and RPS20 in memory B cells; a downregulation of MT-ATP8 and HLA-DQA1 in naïve cells; a downregulation of RPS4Y1 in activated B cells; and a downregulation of IGHV3-73 in the plasma cells. (2) We then evaluated the clonality of B cells with no strong differences found. (3) We appraised the B cell immunoglobulin distribution, finding an increased ratio of IgA + IgG to IgD + IgM in the COVID ECMO D patient, suggesting an intensive memory antibody response. (4) Lastly, we assessed the V(D)J rearrangements of the BCR heavy chain in our four samples and found that the IGHV3, IGHJ4, and IGHD3/IGHD2 families were the most frequently used.
The study’s main weakness is the small sample size, which prevents the results from being meaningful for stating any of the processes involved in the immune response to SARS-CoV-2. However, given the uniqueness of the disease’s characteristics and the molecular approach used, single-cell sequencing, which is based on a large number of events/cells and analyzes each cell as an independent sample, we hope that our computational analysis will provide novel insights for large population cohort studies. Our findings uncover that the V(D)J rearrangements of the BCR heavy chain of ECMO COVID-19 patients analyzed are consistent with other studies focusing on subjects infected with and/or vaccinated against SARS-CoV-2. As a result, we plan to conduct future clinical studies on large cohorts of COVID-19 patients in order to better analyze and compare the BCR repertoire in patient subgroups.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cimb45020095/s1, Figure S1: Selective usage of J gene segments in each sample; Figure S2: Selective usage of V gene segments in each sample; Figure S3: Selective usage of D gene segments in each sample.

Author Contributions

Conceptualization, A.G. and M.B.; methodology, A.G., C.C., M.B. and R.B.; software, N.C.; validation, C.C. and A.G.; formal analysis, N.C., A.G., R.B. and M.B.; investigation, A.G. and M.B.; writing—original draft preparation, M.B., R.B. and A.G.; writing—review and editing, A.G., M.B. and P.G.C.; visualization, M.B., R.B. and A.G.; supervision, A.G., M.B. and P.G.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Ministero della Salute from current research funds in 2022 (Ricerca Corrente 2022) and by EU funding within the MUR PNRR Extended Partnership initiative on Emerging Infectious Diseases (Project no. PE00000007, INF-ACT).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and was approved by the IRCCS ISMETT Institutional Research Review Board (IRRB 00/21) and the Ethics Committee of ISMETT, with all enrolled individuals signing the written informed consent form.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.

Acknowledgments

We would like to thank Stefano Anzani for assistance in using R packages.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Tsai, P.-H.; Lai, W.-Y.; Lin, Y.-Y.; Luo, Y.-H.; Chen, H.-K.; Chen, Y.-M.; Lai, Y.-C.; Kuo, L.-C.; Chen, S.-D.; Chang, K.-J.; et al. Clinical manifestation and disease progression in COVID-19 infection. J. Chin. Med. Assoc. 2021, 84, 3–8. [Google Scholar] [CrossRef] [PubMed]
  2. Martucci, G.; Słomka, A.; Lebowitz, S.E.; Raffa, G.M.; Malvindi, P.G.; Coco, V.L.; Swol, J.; Żekanowska, E.; Lorusso, R.; Wierzba, W.; et al. COVID-19 and Extracorporeal Membrane Oxygenation. Adv. Exp. Med. Biol. 2021, 1353, 173–195. [Google Scholar]
  3. Guidelines Review Committee. Clinical Management of COVID-19: Living Guideline. In WHO Guidelines; Guidelines Review Committee: Geneva, Switzerland, 2022. [Google Scholar]
  4. Gabelloni, M.; Faggioni, L.; Cioni, D.; Mendola, V.; Falaschi, Z.; Coppola, S.; Corradi, F.; Isirdi, A.; Brandi, N.; Coppola, F.; et al. Extracorporeal membrane oxygenation (ECMO) in COVID-19 patients: A pocket guide for radiologists. Radiol. Med. 2022, 127, 369–382. [Google Scholar] [PubMed]
  5. Barbaro, R.P.; MacLaren, G.; Boonstra, P.S.; Iwashyna, T.J.; Slutsky, A.S.; Fan, E.; Bartlett, R.H.; Tonna, J.E.; Hyslop, R.; Fanning, J.J.; et al. Extracorporeal membrane oxygenation support in COVID-19: An international cohort study of the Extracorporeal Life Support Organization registry. Lancet 2020, 396, 1071–1078. [Google Scholar] [CrossRef]
  6. Schmidt, M.; Hajage, D.; Lebreton, G.; Monsel, A.; Voiriot, G.; Levy, D.; Baron, E.; Beurton, A.; Chommeloux, J.; Meng, P.; et al. Extracorporeal membrane oxygenation for severe acute respiratory distress syndrome associated with COVID-19: A retrospective cohort study. Lancet Respir. Med. 2020, 8, 1121–1131. [Google Scholar] [CrossRef] [PubMed]
  7. Ramanathan, K.; Shekar, K.; Ling, R.R.; Barbaro, R.P.; Wong, S.N.; Tan, C.S.; Rochwerg, B.; Fernando, S.M.; Takeda, S.; MacLaren, G.; et al. Correction to: Extracorporeal membrane oxygenation for COVID-19: A systematic review and meta-analysis. Crit. Care 2021, 25, 375. [Google Scholar] [PubMed]
  8. Nielsen, S.C.; Yang, F.; Jackson, K.J.; Hoh, R.A.; Röltgen, K.; Jean, G.H.; Stevens, B.A.; Lee, J.Y.; Rustagi, A.; Rogers, A.J.; et al. Human B Cell Clonal Expansion and Convergent Antibody Responses to SARS-CoV-2. Cell Host Microbe 2020, 28, 516–525. [Google Scholar]
  9. Chen, H.; Zhang, Y.; Ye, A.Y.; Du, Z.; Xu, M.; Lee, C.-S.; Hwang, J.K.; Kyritsis, N.; Ba, Z.; Neuberg, D.; et al. BCR selection and affinity maturation in Peyer’s patch germinal centres. Nature 2020, 582, 421–425. [Google Scholar]
  10. Roth, D.B. V(D)J Recombination: Mechanism, Errors, and Fidelity. Microbiol. Spectr. 2014, 2, 313–324. [Google Scholar]
  11. Petrova, V.N.; Muir, L.; McKay, P.F.; Vassiliou, G.S.; Smith, K.G.; Lyons, P.A.; Russell, C.A.; Anderson, C.A.; Kellam, P.; Bashford-Rogers, R.J. Combined Influence of B-Cell Receptor Rearrangement and Somatic Hypermutation on B-Cell Class-Switch Fate in Health and in Chronic Lymphocytic Leukemia. Front. Immunol. 2018, 9, 1784. [Google Scholar] [CrossRef]
  12. Zheng, G.X.; Terry, J.M.; Belgrader, P.; Ryvkin, P.; Bent, Z.W.; Wilson, R.; Ziraldo, S.B.; Wheeler, T.D.; McDermott, G.P.; Zhu, J.; et al. Massively parallel digital transcriptional profiling of single cells. Nat. Commun. 2017, 8, 14049. [Google Scholar] [CrossRef]
  13. Bulati, M.; Caruso, C.; Colonna-Romano, G. From lymphopoiesis to plasma cells differentiation, the age-related modifications of B cell compartment are influenced by “inflamm-ageing”. Ageing Res. Rev. 2017, 36, 125–136. [Google Scholar] [PubMed]
  14. Tjiam, M.C.; Fernandez, S.; French, M.A. Characterising the Phenotypic Diversity of Antigen-Specific Memory B Cells before and after Vaccination. Front. Immunol. 2021, 12, 738123. [Google Scholar] [PubMed]
  15. Tang, T.F.; Chan, Y.T.; Cheong, H.C.; Cheok, Y.Y.; Anuar, N.A.; Looi, C.Y.; Gan, G.G.; Wong, W.F. Regulatory network of BLIMP1, IRF4, and XBP1 triad in plasmacytic differentiation and multiple myeloma pathogenesis. Cell. Immunol. 2022, 380, 104594. [Google Scholar] [CrossRef] [PubMed]
  16. Qi, F.; Cao, Y.; Zhang, S.; Zhang, Z. Single-cell analysis of the adaptive immune response to SARS-CoV-2 infection and vaccination. Front. Immunol. 2022, 13, 964976. [Google Scholar] [PubMed]
  17. Kadam, U.S.; Lossie, A.C.; Schulz, B.; Irudayaraj, J. Gene Expression Analysis Using Conventional and Imaging Methods. In DNA and RNA Nanobiotechnologies in Medicine: Diagnosis and Treatment of Diseases; Erdmann, V.A., Barciszewski, J., Eds.; Springer: Berlin/Heidelberg, Germany, 2013; pp. 141–162. [Google Scholar]
  18. Galber, C.; Acosta, M.J.; Minervini, G.; Giorgio, V. The role of mitochondrial ATP synthase in cancer. Biol. Chem. 2020, 401, 1199–1214. [Google Scholar] [CrossRef]
  19. Galperin, M.; Farenc, C.; Mukhopadhyay, M.; Jayasinghe, D.; Decroos, A.; Benati, D.; Tan, L.L.; Ciacchi, L.; Reid, H.H.; Rossjohn, J.; et al. CD4(+) T cell-mediated HLA class II cross-restriction in HIV controllers. Sci. Immunol. 2018, 3, eaat0687. [Google Scholar] [CrossRef]
  20. Kwok, W.W.; Tan, V.; Gillette, L.; Littell, C.T.; Soltis, M.A.; LaFond, R.B.; Yang, J.; James, E.A.; DeLong, J.H. Frequency of epitope-specific naive CD4(+) T cells correlates with immunodominance in the human memory repertoire. J. Immunol. 2012, 188, 2537–2544. [Google Scholar] [CrossRef]
  21. Muehling, L.M.; Mai, D.T.; Kwok, W.W.; Heymann, P.W.; Pomes, A.; Woodfolk, J.A. Circulating Memory CD4+ T Cells Target Conserved Epitopes of Rhinovirus Capsid Proteins and Respond Rapidly to Experimental Infection in Humans. J. Immunol. 2016, 197, 3214–3224. [Google Scholar]
  22. Topalian, S.L.; Gonzales, M.I.; Parkhurst, M.; Li, Y.F.; Southwood, S.; Sette, A.; Rosenberg, S.A.; Robbins, P.F. Melanoma-specific CD4+ T cells recognize nonmutated HLA-DR-restricted tyrosinase epitopes. J. Exp. Med. 1996, 183, 1965–1971. [Google Scholar] [CrossRef]
  23. Shams, H.; Klucar, P.; Weis, S.E.; Lalvani, A.; Moonan, P.K.; Safi, H.; Wizel, B.; Ewer, K.; Nepom, G.T.; Lewinsohn, D.M.; et al. Characterization of a Mycobacterium tuberculosis peptide that is recognized by human CD4+ and CD8+ T cells in the context of multiple HLA alleles. J. Immunol. 2004, 173, 1966–1977. [Google Scholar]
  24. Abelin, J.G.; Harjanto, D.; Malloy, M.; Suri, P.; Colson, T.; Goulding, S.P.; Creech, A.L.; Serrano, L.R.; Nasir, G.; Nasrullah, Y.; et al. Defining HLA-II Ligand Processing and Binding Rules with Mass Spectrometry Enhances Cancer Epitope Prediction. Immunity 2019, 51, 766–779. [Google Scholar] [CrossRef] [PubMed]
  25. Schulze zur Wiesch, J.; Lauer, G.M.; Day, C.L.; Kim, A.Y.; Ouchi, K.; Duncan, J.E.; Wurcel, A.G.; Timm, J.; Jones, A.M.; Mothe, B.; et al. Broad repertoire of the CD4+ Th cell response in spontaneously controlled hepatitis C virus infection includes dominant and highly promiscuous epitopes. J. Immunol. 2005, 175, 3603–3613. [Google Scholar] [PubMed]
  26. Yong, W.H.; Shabihkhani, M.; Telesca, D.; Yang, S.; Tso, J.L.; Menjivar, J.C.; Wei, B.; Lucey, G.M.; Mareninov, S.; Chen, Z.; et al. Ribosomal Proteins RPS11 and RPS20, Two Stress-Response Markers of Glioblastoma Stem Cells, Are Novel Predictors of Poor Prognosis in Glioblastoma Patients. PLoS ONE 2015, 10, e0141334. [Google Scholar] [CrossRef] [PubMed]
  27. Ma, J.; Bai, H.; Gong, T.; Mao, W.; Nie, Y.; Zhang, X.; Da, Y.; Wang, X.; Qin, H.; Zeng, Q.; et al. Novel skewed usage of B-cell receptors in COVID-19 patients with various clinical presentations. Immunol. Lett. 2022, 249, 23–32. [Google Scholar] [CrossRef]
  28. Liu, H.; Kaku, C.I.; Song, G.; Yuan, M.; Andrabi, R.; Burton, D.R.; Walker, L.M.; Wilson, I.A. Human antibodies to SARS-CoV-2 with a recurring YYDRxG motif retain binding and neutralization to variants of concern including Omicron. Commun. Biol. 2022, 5, 766. [Google Scholar]
  29. He, B.; Liu, S.; Xu, M.; Hu, Y.; Lv, K.; Wang, Y.; Ma, Y.; Zhai, Y.; Yue, X.; Liu, L.; et al. Comparative global B cell receptor repertoire difference induced by SARS-CoV-2 infection or vaccination via single-cell V(D)J sequencing. Emerg. Microbes Infect. 2022, 11, 2007–2020. [Google Scholar] [PubMed]
Figure 1. The antigen-binding domain of immunoglobulins is composed of two polypeptide chains, namely heavy and light chains. The exons that encode the antigen-binding domains of the heavy chain are assembled from V (variable), D (diversity), and J (joining) gene segments by a process defined as “cut-and-paste” DNA rearrangements. This process, named V(D)J recombination, selects a pair of segments, introduces double-strand breaks next to each segment, deletes the intermediate DNA, and ligates the segments together. Rearrangements take place in a well-ordered way, with D-to-J joining proceeding before a V segment is joined to the rearranged D–J segments. The rearrangement of the light chain is identical except for the absence of D gene fragments.
Figure 1. The antigen-binding domain of immunoglobulins is composed of two polypeptide chains, namely heavy and light chains. The exons that encode the antigen-binding domains of the heavy chain are assembled from V (variable), D (diversity), and J (joining) gene segments by a process defined as “cut-and-paste” DNA rearrangements. This process, named V(D)J recombination, selects a pair of segments, introduces double-strand breaks next to each segment, deletes the intermediate DNA, and ligates the segments together. Rearrangements take place in a well-ordered way, with D-to-J joining proceeding before a V segment is joined to the rearranged D–J segments. The rearrangement of the light chain is identical except for the absence of D gene fragments.
Cimb 45 00095 g001
Figure 2. Features of B cell subsets. (A) Violin plots show the expression distribution of canonical cell markers in B cell subsets. (B) Proportion of B cell clusters in each sample. (C) t-SNE projection of B cells from all four samples. Each dot corresponds to a single cell, colored according to cell clusters. (D) Heat map representation of differentially expressed genes of B cell populations in the COVID ECMO S, COVID ECMO D, COVID R, and the healthy control. Genes with a maximum adjusted p-value of 0.01 and an absolute value of log2 (fold change) >0.5 were considered to be differentially expressed genes.
Figure 2. Features of B cell subsets. (A) Violin plots show the expression distribution of canonical cell markers in B cell subsets. (B) Proportion of B cell clusters in each sample. (C) t-SNE projection of B cells from all four samples. Each dot corresponds to a single cell, colored according to cell clusters. (D) Heat map representation of differentially expressed genes of B cell populations in the COVID ECMO S, COVID ECMO D, COVID R, and the healthy control. Genes with a maximum adjusted p-value of 0.01 and an absolute value of log2 (fold change) >0.5 were considered to be differentially expressed genes.
Cimb 45 00095 g002
Figure 3. Features of B cell subpopulations. (A) Heat map representation of differentially expressed genes of B cell subpopulations (naïve cells, activated B cells, memory B cells, and plasma cells) of the COVID ECMO S, COVID ECMO D, COVID R, and healthy control. Genes with a maximum adjusted p-value of 0.01 and an absolute value of log2 (fold change) >0.5 were considered to be differentially expressed genes. (B) The graphs show the sequencing results of highly deregulated mRNA in the COVID ECMO S, COVID ECMO D, COVID R, and healthy control.
Figure 3. Features of B cell subpopulations. (A) Heat map representation of differentially expressed genes of B cell subpopulations (naïve cells, activated B cells, memory B cells, and plasma cells) of the COVID ECMO S, COVID ECMO D, COVID R, and healthy control. Genes with a maximum adjusted p-value of 0.01 and an absolute value of log2 (fold change) >0.5 were considered to be differentially expressed genes. (B) The graphs show the sequencing results of highly deregulated mRNA in the COVID ECMO S, COVID ECMO D, COVID R, and healthy control.
Cimb 45 00095 g003
Figure 4. (A) t-SNE plot of B cells colored by immunoglobulin heavy chain expression. (B) The proportion of different immunoglobulin isotypes in each sample. (C) The pie plot shows the proportion of clonally expanded B cells in the COVID ECMO S, COVID ECMO D, and COVID R patients and in the healthy control. (D) The distribution of immunoglobulin heavy chain (IGH) CDR3 length. (E) The bubble chart shows the usage frequency of V–J gene combination in each sample.
Figure 4. (A) t-SNE plot of B cells colored by immunoglobulin heavy chain expression. (B) The proportion of different immunoglobulin isotypes in each sample. (C) The pie plot shows the proportion of clonally expanded B cells in the COVID ECMO S, COVID ECMO D, and COVID R patients and in the healthy control. (D) The distribution of immunoglobulin heavy chain (IGH) CDR3 length. (E) The bubble chart shows the usage frequency of V–J gene combination in each sample.
Cimb 45 00095 g004
Table 1. Characteristics of patients and control.
Table 1. Characteristics of patients and control.
ID PatientsAgeGenderComorbiditiesSeverity Clinical ProfileCo-InfectionsTreatmentsHospitalization (Days)
COVID ECMO S74MVasculopathy and DiabetesCARDS-Dexamethasone30 total
(10 ECMO)
COVID ECMO D67MHTNCARDS and SepsisMDR Xanthomonas (VAE), MDR Acinetobacter (UTI), C.albicans (BSI)CPAP and Dexamethasone40 total
(12 NIV + 7 ECMO)
COVID R42F-A---
HEALTHY CTRL45M-----
Abbreviations: HTN, hypertension; C, critical; A, asymptomatic; ARDS, acute respiratory distress syndrome; MDR, multidrug-resistant; VAE, ventilator-associated events; UTI, urinary tract infection; BSI, bloodstream infection; CPAP, continuous positive airways pressure; NIV, non-invasive ventilation; ECMO, extra corporeal membrane oxygenation.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gallo, A.; Cuscino, N.; Carcione, C.; Busà, R.; Conaldi, P.G.; Bulati, M. Proof-of-Concept Analysis of B Cell Receptor Repertoire in COVID-19 Patients Undergoing ECMO by Single-Cell V(D)J and Gene Expression Sequencing. Curr. Issues Mol. Biol. 2023, 45, 1471-1482. https://doi.org/10.3390/cimb45020095

AMA Style

Gallo A, Cuscino N, Carcione C, Busà R, Conaldi PG, Bulati M. Proof-of-Concept Analysis of B Cell Receptor Repertoire in COVID-19 Patients Undergoing ECMO by Single-Cell V(D)J and Gene Expression Sequencing. Current Issues in Molecular Biology. 2023; 45(2):1471-1482. https://doi.org/10.3390/cimb45020095

Chicago/Turabian Style

Gallo, Alessia, Nicola Cuscino, Claudia Carcione, Rosalia Busà, Pier Giulio Conaldi, and Matteo Bulati. 2023. "Proof-of-Concept Analysis of B Cell Receptor Repertoire in COVID-19 Patients Undergoing ECMO by Single-Cell V(D)J and Gene Expression Sequencing" Current Issues in Molecular Biology 45, no. 2: 1471-1482. https://doi.org/10.3390/cimb45020095

APA Style

Gallo, A., Cuscino, N., Carcione, C., Busà, R., Conaldi, P. G., & Bulati, M. (2023). Proof-of-Concept Analysis of B Cell Receptor Repertoire in COVID-19 Patients Undergoing ECMO by Single-Cell V(D)J and Gene Expression Sequencing. Current Issues in Molecular Biology, 45(2), 1471-1482. https://doi.org/10.3390/cimb45020095

Article Metrics

Back to TopTop