Next Article in Journal
Phenotypes of Floral Nectaries in Developmental Mutants of Legumes and What They May Tell about Genetic Control of Nectary Formation
Previous Article in Journal
Prevalence of Periodontal Pathogens in Slovak Patients with Periodontitis and Their Possible Aspect of Transmission from Companion Animals to Humans
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Machine-Learning-Assisted Analysis of TCR Profiling Data Unveils Cross-Reactivity between SARS-CoV-2 and a Wide Spectrum of Pathogens and Other Diseases

1
Laboratory of Hygiene and Epidemiology, Faculty of Medicine, University of Thessaly, 41222 Larisa, Greece
2
Laboratory of Genetics, Department of Biology, University of Patras, 26500 Patras, Greece
3
Department of Immunology & Histocompatibility, Faculty of Medicine, University of Thessaly, 41500 Larisa, Greece
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
These authors contributed equally to this work.
Biology 2022, 11(10), 1531; https://doi.org/10.3390/biology11101531
Submission received: 24 September 2022 / Revised: 13 October 2022 / Accepted: 17 October 2022 / Published: 19 October 2022
(This article belongs to the Section Bioinformatics)

Abstract

:

Simple Summary

For the last two years, COVID-19 has been rigorously studied aiming to identify novel prognostic and therapeutic avenues. Recently, T cell receptor profiling has emerged as a method to associate adaptive immunity with COVID-19 progression and severity. Such data are typically analyzed to explore T cell receptor properties and characteristics in the context of SARS-CoV-2 infection. However, the equally informative alternative analytic strategy of identifying any preferential recognition of viral antigens by the T-cell-mediated immune response is mostly overlooked. In this study, we propose a novel Machine-Learning-oriented approach for analyzing T cell receptor repertoires that is based on the concept of utilising the level at which each SARS-CoV-2 antigen is recognised by the available T cell receptors in each sample from COVID-19-convalescent and healthy cohorts. This approach also allowed us to observe a group of T cell receptors capable of recognising SARS-CoV-2 antigens that were already established in samples from the healthy cohort, leading us to the cross-reactivity phenomenon hypothesis. To explore this, all T cell receptors were examined for being able to recognise antigens from other pathogens and diseases, unveiling evidence of putative cross-reactivity with M. tuberculosis and Influenza virus, among others.

Abstract

During the last two years, the emergence of SARS-CoV-2 has led to millions of deaths worldwide, with a devastating socio-economic impact on a global scale. The scientific community’s focus has recently shifted towards the association of the T cell immunological repertoire with COVID-19 progression and severity, by utilising T cell receptor sequencing (TCR-Seq) assays. The Multiplexed Identification of T cell Receptor Antigen (MIRA) dataset, which is a subset of the immunoACCESS study, provides thousands of TCRs that can specifically recognise SARS-CoV-2 epitopes. Our study proposes a novel Machine Learning (ML)-assisted approach for analysing TCR-Seq data from the antigens’ point of view, with the ability to unveil key antigens that can accurately distinguish between MIRA COVID-19-convalescent and healthy individuals based on differences in the triggered immune response. Some SARS-CoV-2 antigens were found to exhibit equal levels of recognition by MIRA TCRs in both convalescent and healthy cohorts, leading to the assumption of putative cross-reactivity between SARS-CoV-2 and other infectious agents. This hypothesis was tested by combining MIRA with other public TCR profiling repositories that host assays and sequencing data concerning a plethora of pathogens. Our study provides evidence regarding putative cross-reactivity between SARS-CoV-2 and a wide spectrum of pathogens and diseases, with M. tuberculosis and Influenza virus exhibiting the highest levels of cross-reactivity. These results can potentially shift the emphasis of immunological studies towards an increased application of TCR profiling assays that have the potential to uncover key mechanisms of cell-mediated immune response against pathogens and diseases.

1. Introduction

Since SARS-CoV-2 was initially reported in Wuhan, China, there have been 515 million confirmed cases and 6.2 million deaths worldwide as of May 2022, according to the Johns Hopkins Coronavirus Resource Centre [1]. Individuals infected with SARS-CoV-2 exhibit a wide spectrum of responses, from asymptomatic to requiring admission to an intensive care unit [2]. Both the research community and pharmaceutical industry have been rigorously studying COVID-19 [3] and implications of SARS-CoV-2 infection [4], aiming to identify novel prognostic and therapeutic avenues [5,6,7], while also exerting a massive effort to bring a plethora of vaccination schemes to the public within a very limited timeframe [8].
During the past year, there has been a shift in published research highlighting the need to better understand the T cell immunological profile association with COVID-19 progression and severity [9,10,11,12,13,14,15]. T cell immunity appears to be a much more sensitive indicator of past infections in comparison with antibody response. High-throughput methods approaching T cell response can be informative by correlating concepts of clonal depth, breadth and dynamics with symptoms and disease severity [16]. Furthermore, previous studies associated with other viruses such as Middle East Respiratory Syndrome (MERS) and SARS-CoV-1 indicate that coronavirus-specific T cells appear to have long term persistence [17,18]. The same phenomenon also seems to take place in SARS-CoV-2 biology [19,20]. These observations could shape the hypothesis of cross-reactivity between different pathogens, where past infection or vaccination could be protective through long-lived T cell clones. The cross-reactivity phenomenon between SARS-CoV-2 and other coronaviruses has been reported in the literature. Neutralising antibodies isolated from memory-B cells of a SARS-CoV-1 infected individual have been described to react with SARS-CoV-2 surface glycoprotein [21]. In addition, cross-reactive T cells recognising SARS-CoV-2 seem to be acquired during previous infections by other human coronaviruses in 20% to 50% of unexposed individuals around the world [22,23,24,25,26,27]. Those preexisting cells may affect the clinical manifestations of COVID-19 infection.
The adaptive immune response to pathogenic infections is largely dependent on the CD4+ and CD8+ T cell subfamilies [28]. Upon activation, CD8+ T cells can exterminate infected cells and form the long-term memory T cell subpopulation. Conversely, CD4+ T cells control the function of myeloid cells, support CD8+ response and play a key role in the selection of antigen-specific B cells which contribute to the host organism’s neutralising antibody arsenal. T cell receptors (TCRs) are proteins localised on the surface of T cells that are products of recombined genomic sequences during the T cell developmental process. The uniqueness of each TCR sequence essentially controls the T cells’ specificity. TCRs recognise peptides presented by the major histocompatibility complex (MHC) on the surface of most cell types (MHC class-I recognised by CD8-expressing T cytotoxic cells), or on the surface of antigen-presenting cells (APC) (MHC class-II recognised by CD4-expressing T helper cells). The ability of TCRs to recognise more than one peptide-MHC structure defines cross-reactivity [29]. The cross-reactivity of T cells is considered as one arm of the well-described heterologous immunity, namely the immunity that can develop towards one pathogen after exposure to non-identical pathogens [30]. The other arm concerns the bystander activation of T cells, that can be caused by independent activation via released cytokines, or by low-affinity recognition of pMHC [29,31]. Heterologous immunity has been well established in viral infections in several animal models, but also in human viral infections where cross-reactivity of T cells could potentially influence the protection or severity of virus-associated immunopathology [32,33,34].
T cell cross-reactivity, however, comes at a cost. Pathogen-induced autoimmune disorders may also result from cross-reactive T cells, following the immune system’s initial reaction to the pathogen. This phenomenon has been termed molecular mimicry, where peptides derived from pathogens can activate autoreactive T cells due to structural similarity between pathogenic and self-peptides, causing autoimmune diseases or accelerating a previously initiated autoimmune process [35,36,37]. The mechanics of T cells cross-reactivity involves changes in complementarity-determining region (CDR) loop conformation, altered TCR docking on the pMHC, flexible changes in pMHC and structural degeneracy [29,38,39,40]. The binding of TCR with peptides differs concerning their affinity, since the bound peptides could be considerably different in their chemistry [38,39]. In this context, several previous studies have reported the presence of cross-reactivity among different viruses, or epitopes derived from the same pathogen [41,42,43].
Recent evidence in the literature highlights [11,12,13,15] an effort of the research community to explore the TCR repertoire in the context of several levels of COVID-19 infection severity, by utilising data from high-throughput TCR-Sequencing (TCR-Seq) assays. The majority of the work focuses on developing computational methods to unveil differences and similarities between healthy and infected subjects related to the TCR repertoire diversity, CDR3 length distribution and the V and J gene segment preference [11,12,13,14]. Other studies have attempted to combine the aforementioned TCR-related statistics with Machine Learning (ML), aiming to provide predictive tools to distinguish healthy and infected subjects [44,45]. To the best of our knowledge, however, no studies exist in the literature which attempt to approach this field from the viewpoint of the extent each SARS-CoV-2 antigen is recognised by the TCR repertoire of COVID-19 infected and non-infected individuals.
In this study, we propose a novel ML-based approach for analysing TCR repertoires derived from the Multiplexed Identification of T cell Receptor Antigen (MIRA) specificity assay [46] (Figure 1A,B). Such data have been recently generated through academic partnerships with the industrial sector, and were released in the form of freely accessible databases such as immunoACCESS© [15]. This resource includes the immunoSEQ dataset of sequenced TCR beta chain (TCRb) repertoires from COVID-19-exposed, -infected or -recovered individuals who have participated in the Immune Response Action to COVID-19 Events study, as well as thousands of patients’ blood samples collected by international institutions globally. The immunoACCESS© repository also includes the MIRA dataset which is complementary to immunoSEQ; the repository catalogues TCRb sequences and TCR specific information about the peptide’s interaction in amino acid level, as well as the targeting epitope molecule it comes from. Our approach (Figure 1B) focused on the MIRA dataset and utilised, for the first time, the level at which each SARS-CoV-2 antigen is recognised by the available TCRs in each sample, to train several ML algorithms that can distinguish samples from COVID-19-convalescent and healthy (no known exposure) cohorts with over 85% accuracy. The module for highlighting the importance of each antigen revealed that the TCR clones recognising ORF7b, nucleocapsid phosphoprotein, as well as ORF1ab, ORF10, ORF3a and membrane glycoprotein to a lesser degree, play a key role for the classification task. Additionally, all MIRA TCRs, regardless of their assigned cohort, were further analysed to determine their potential for recognising epitopes originating from pathogens other than SARS-CoV-2 (Figure 1C). To this end, data from public TCR databases were processed to unveil evidence of putative cross-reactivity between SARS-CoV-2 and multiple pathogens and other diseases, with M. tuberculosis and Influenza virus being the most cross-reactive.

2. Materials and Methods

2.1. Data Collection and Pre-Processing

TCRs that are able to bind to SARS-CoV-2 epitopes were retrieved from the immuneACCESS© database [15] (Figure 1A). These SARS-CoV-2 specific TCRs are part of the MIRA dataset which is based on 144 samples (experiments) obtained from cohorts with exposed subjects and healthy controls (Figure 1B and Figure 2A). It should be noted that for certain subjects, more than one sample is available in MIRA. Specifically, 90 samples originate from COVID-19-convalescent, 39 from healthy (no known exposure), 4 from COVID-19-acute, 8 from COVID-19-non-acute and 3 from COVID-19-exposed subjects. Samples from the COVID-19-acute, COVID-19-non-acute and COVID-19-exposed cohorts were excluded from the analyses presented herein, due to their limited number. TCR sequences were initially filtered to keep CDR3 regions delimited by a conserved cysteine at the start and a conserved phenylalanine or tryptophan at the end (anchors of CDR3 region). Unproductive CDR3 segments and sequences containing special characters not corresponding to amino acids (X, #, *, etc.) were also excluded. MIRA was further filtered to keep information associated only with functional V genes, removing information related to pseudogenes and ORFs according to the immunogenetics information system (IMGT) [47]. Clonotypes with ambiguous V gene family members (denoted with X) were also removed. The analysis was focused on the remaining 130,072 (120,128 unique) TCRs detected in the studied cohorts (28 healthy and 85 convalescent subjects). A total of 1006 unique TCRs originate from the minigene-detail file, 114,411 from peptide-detail-ci (CD8+) and 4651 from peptide-detail-cii (CD4+). In addition, 34 unique TCRs were found in both minigene-detail and peptide-detail-ci files, and 26 TCRs were found in both peptide-detail-ci and peptide-detail-cii files. The unique number of TCRs per 1000 TCRs in each sample is depicted in Figure 2B.
All remaining TCRs were further analysed from the SARS-CoV-2 antigens’ point of view. For every sample, each TCR was assigned to an antigen category (N = 11) based on its epitope recognition ability (Figure 1B and Figure 2C). Some TCRs were able to recognise epitopes from different antigens. For these cases, the assignment to the corresponding antigens was weighted based on the number of antigens. The number of TCRs per antigen was normalised based on the total number of TCRs per sample. This approach enabled the representation of each sample by an 11-dimensional vector and facilitated the aggregation of a dataset used to build several ML models that can classify samples into the convalescent and healthy categories and explore the underlying biology (Figure 1B).
Data catalogues derived from three public TCR databases with immunogenetic information were downloaded to investigate the MIRA dataset’s TCR involvement in immune response during other infections. Pathology-associated TCR database McPAS [48] is a manually curated dataset of TCR sequences associated with various pathological manifestations, containing information about the T cell type, organ or tissue antigen target and related MHC molecules (version 4 January 2022). TCR3d [49] is a structural repertoire database including experimentally determined TCR structures and complexes. It also contains TCR sequences and related antigenic peptides and MHC molecules. We used a data frame derived from TCR3d focused on TCRb CDR3 sequence level and the association with viruses (version 13 January 2022). VDJdb [50] is a database containing TCR sequences, their cognate antigens and related MHC molecules (version 22 March 2022). All three aforementioned databases were filtered to keep information about immune response in human species and CDR3 sequences associated with TCRb; the MIRA filtering strategy described earlier in this section was also applied here. Additionally, in the case of VDJdb, CDR3 sequences with zero confidence score were removed from downstream analyses. As stated in VDJdb’s documentation, the higher the score the more confidence there is in the antigen specificity annotation of a given TCR clonotype. Some VDJdb sequences were processed during fixing steps according to IMGT nomenclature and were included in our analysis. Furthermore, the McPAS sequences associated with antigen identification method id “3” were removed in accordance with the database’s recommendation for confidence in the accuracy of the data.

2.2. Machine Learning Model Training and Feature Importance Estimation

All 90 COVID-19-convalescent samples were labelled as positives, and the 39 healthy (no known exposure) samples as negatives (Figure 1B and Figure 2A). Each sample consists of an 11-dimensional vector of normalised values, for every SARS-CoV-2 antigen, that represents the percentage of TCRs recognising the antigen’s epitopes from the total amount of SARS-CoV-2-specific TCRs in the sample (Figure 1B and Figure 2C). Visualisation of the data in the principal component analysis (PCA) space is depicted in Figure 2C. Principal component loadings can be found in Table S1. Both positive and negative samples were randomly divided into training and test sets based on a 7:3 ratio. This process was repeated 20 times to generate an equal amount of training/test set combinations and control for any bias that could result from the splitting process (Figure 1B).
These sets were used to train and evaluate a total of 100 models based on popular ML algorithmic families (Figure 1B): Gaussian Naive Bayes (GaussianNB), Decision Trees (DT), K-Nearest Neighbours (KNN), Random Forests (RF) and Support Vector Machines (SVM). The hyperparameters of each model were tuned based on a grid search approach, and balanced accuracy was the target metric for choosing the best performing model. For GaussianNB and the variance smoothing parameter, the grid search was run on the values 1 × 10−9, 1 × 10−8 and 1 × 10−7. For DT, the benchmarked values for maximum depth were 10, 30 and 90 with the maximum features parameter set to None. In the case of KNN, the k-neighbours parameter values were 2, 5 and 10. The algorithmic options for selecting nearest neighbours were ball_tree, kde_tree and brute, and the distance metric parameter values were 1 (manhattan) and 2 (euclidean). For RF, the maximum depth values were 10, 30 and 90; the number of estimators were 10, 50 and 100. The maximum features parameter was set to None. The SVM kernel was set to radial basis function and the different ‘C’ parameter values were 0.1, 1, 10 and 100. The gamma values were 1, 0.1, 0.01 and 0.001. All models were trained based on a 10-fold cross validation scheme repeated 10 times. Subsequently, the best performing model was evaluated on its designated test set.
This approach resulted in the calculation of performance metrics such as balanced accuracy, precision, sensitivity, specificity and NPV (Figure 3A,B). Since 20 models were trained and benchmarked for each ML algorithm, all performance plots depict the metrics’ score distributions from all test sets, providing hints of putative training/test split bias and data heterogeneity. Additionally, the importance of each feature was estimated by repeated (N = 50) random shuffles of single feature values to assess the decrease or increase of the models’ performance (Figure 3C). For each ML algorithm, features with a median score above 0.01 were selected for a second round of training and evaluation, following the previously aforementioned strategy (Figure 3D).

2.3. Identification of TCRs That Recognise Epitopes from Antigens of SARS-CoV-2 and Other Pathogens and Diseases

Several definitions of the clonotype concept exist as various studies approach it with different immunogenetic characteristics. The MIRA dataset contains TCRb sequences targeting specific SARS-CoV-2 epitopes and the TCR bioidentity is described by CDR3 amino acid sequence, V and J genes. In this analysis, every clonotype consists of sequences characterised by the same V gene family member, the same J gene family member and the same CDR3 sequence in amino acid level. Clonal expansion of every clonotype in each sample (experiment) was calculated by counting the times it appears divided by the total MIRA clonotypes detected. Hence, clonal expansion was defined as a measure of the T cell proportion expressing a specific TCRb sequence. The mean clonal expansion of each clonotype was assessed by calculating the mean of expansion values from all experiments where this clonotype was detected (Figure 4A).
The six most frequent clonotypes were further analysed from the viewpoint of antigen recognition, clonal expansion range and clinical impact. The level of clonal expansion was also calculated separately for the convalescent and healthy cohorts (Figure 4B). Additionally, each clonotype was characterised for the clinical cohort distribution and statistical significance was also determined with Fisher’s exact test (Figure 4C).
The cross-reactivity phenomenon between SARS-CoV-2 and other pathogens was also examined. We used the public TCR databases McPAS [48], TCR3d [49] and VDJdb [50] to confirm the existence of different pathogens and other diseases associated with clonotypes targeting SARS-CoV-2 epitopes in the MIRA dataset. Quantitative analysis was conducted calculating the number of unique CDR3-mediated connections associated with each pathology to capture the number of putative cross-reactive TCRs (Figure 5A). In certain cases, a single CDR3 sequence can recognise epitopes from multiple antigens. Thus, the number of SARS-CoV-2 antigen connections with other pathogens’ and diseases’ antigens is not equal to the number of unique MIRA CDR3 sequences (Figure 5 and Figure S1, Tables S2–S5). This approach was applied twice, once by screening all CDR3 sequences (Figure 5, Tables S2 and S5) from the aforementioned databases (derived from CD8+ and CD4+ T cells) and once by targeting only CD8+ or CD4+ T cells, to examine any potential functionality bias (Figure S1, Tables S3 and S4). To identify locations crucial for the interaction with cross-reactive CDR3 regions, the epitopes were aligned to the pathogenic proteins they derive from (Figure 4D,E and Figure 5B). This was the first step to characterise the specific domains and functionality of antigens recognised by the same CDR3 sequences.
It should be noted that these public databases include information for both CD4+ and CD8+ T cells and a data bias exists due to the numerous study results associated with specific pathogenic cases. To briefly mention the most notable cases, there are 16,162 unique CDR3 sequences associated with M. tuberculosis (derived from 14,992 CD4+ and 1183 CD8+ T cells), 3639 with Influenza virus (derived from 155 CD4+ and 3479 CD8+ T cells), 2663 with CMV (derived from 2658 CD8+ T cells), 1583 with HIV (derived from 350 CD4+ and 1514 CD8+ T cells) and 1437 with EBV (derived from 1334 CD8+ T cells). The total number of unique CDR3 sequences does not coincide with the CD4+ and CD8+ subsets, as some CDR3 sequences seem to be associated with both subpopulations. In addition, for some CDR3 sequences there is no information about T cell type in the public TCR repositories. However, separating T cell populations reflects important information about the immune response against the studied pathogens.

3. Results

3.1. MIRA Dataset Exploration Unveils Differentially Recognised SARS-CoV-2 Antigens between Convalescent and Healthy Samples

The MIRA dataset includes 144 samples (experiments) from the immunoACCESS study [15] that are divided into five cohorts: (1) COVID-19-convalescent, (2) healthy (no known exposure), (3) COVID-19-acute, (4) COVID-19-non-acute and (5) COVID-19-exposed (Figure 2A). Most cohorts have a balanced male-to-female ratio; however, there are 26 samples from unknown gender, of which the majority (N = 21) are reported as COVID-19-convalescent cohort. The COVID-19-acute, -non-acute and -exposed cohorts were removed from all subsequent analyses due to low sample numbers. Initial analysis regarding the normalised number of unique TCRs in each sample revealed there is no statistically significant difference between the healthy and convalescent cohorts (Figure 2B).
Most existing studies focus on TCR properties such as the underlying VJ rearranged sequences, CDR3 sequences, CDR3 size and clonal depth/diversity [11,12,13,14]. This information has been frequently used to characterise the TCR repertoire of immune response against COVID-19 and to train ML algorithms that are able to distinguish between samples of distinct cohorts [44,45]. In this study, a different approach was adopted. Rather than using the aforementioned TCR-related data, we exploited the MIRA dataset [15], and the connection between TCRs and the SARS-CoV-2 epitopes to generate an 11-dimensional vector representing the level at which each SARS-CoV-2 antigen is recognised by the TCRs in each sample (Figure 1B and Figure 2C).
The results of this approach revealed that even in samples from the healthy cohort, all SARS-CoV-2 antigens are recognised by TCRs to some extent, suggesting either previous unreported COVID-19 infection of subjects in the healthy cohort, or putative cross-reactivity between SARS-CoV-2 and other pathogens. Interestingly, the ORF1ab and surface glycoprotein are the two antigens recognised by TCRs with the highest overall clonal depth, although the difference in clonal expansion between the two cohorts is not statistically significant. More importantly, the nucleocapsid phosphoprotein, ORF7b, ORF10, ORF8 and envelope exhibit statistically significant differences in the number of TCRs recognising their epitopes between the two cohorts. However, the ORF10, ORF8 and envelope proteins are recognised by TCRs with low clonal depth. The projection of these samples on PCA space (Figure 2D, Table S1) has led to the assumption that this approach could be used to develop a ML-based framework for specifically distinguishing between samples from the two MIRA cohorts (Figure 1B), aiming to extract additional information from MIRA that possibly cannot be derived from the statistical analysis described above.

3.2. Explainable ML Highlights Key SARS-CoV-2 Antigens for Classifying Samples into the Convalescent and Healthy MIRA Cohorts

The strategy of modelling the MIRA dataset involved a repeated process (N = 20) of separately splitting samples from healthy and convalescent cohorts into training and test sets. At each split, five models based on GaussianNB, DT, KNN, RF and SVM were trained and evaluated (Figure 1B).
Using a prediction score cut-off of 0.5 enabled the extraction of performance metrics on each test set (Figure 3A). SVM was the overall best performing algorithm with a median performance of at least 0.75 in all metrics. Notably, the SVM algorithm exhibits 0.856 median balanced accuracy, 0.896 precision, 0.962 sensitivity, 0.75 specificity and 0.9 negative predictive value (NPV). To observe SVM’s performance across the whole spectrum of prediction score thresholds, an incremental cut-off was applied and all metrics were calculated at each step (Figure 3B).
To assess the importance of each feature, a repeated (N = 50) feature value perturbation process was applied for all ML algorithms (Figure 3C). Overall, the most important features are ORF7b, nucleocapsid phosphoprotein as well as ORF1ab, ORF10, ORF3a as well as membrane glycoprotein to a lesser extent. After selecting only the most important features for each algorithm, the training and evaluation process was repeated. This resulted in slightly improved performance for GaussianNB, KNN and SVM algorithms, but not for DT and RF, as expected considering their innate ability to readily rely only on important features (Figure 3D).

3.3. Exploratory Analysis of MIRA TCRs Unveils Evidence of Putative Cross-Reactivity between SARS-CoV-2 and Other Pathogens and Diseases

MIRA TCRs exhibit diverse frequency of occurrence and clonal expansion levels (Figure 4A). In general, the most frequent clonotypes in the dataset present low mean expansion after being triggered with the MIRA assay, while the least frequent clonotypes are associated with the highest expansion levels.
The analysis focused on the six most common clonotypes with frequencies of greater than 11% of the total number of subjects; CASSIRSSYEQYF+V19-01+J02-07 (CD8+ T cell TCR in MIRA), CASSLAGAYEQYF+TCRBV05-01+TCRBJ02-07 (CD8+), CASSLSAPQETQYF+TCRBV27-01+TCRBJ02-05 (CD8+), CASSLSSPQETQYF+TCRBV27-01+ TCRBJ02-05 (CD8+), CASSDRGPNQPQHF+TCRBV27-01+TCRBJ01-05 (CD8+) and CASSDRGPTDTQYF+ TCRBV27-01+TCRBJ02-03 (CD8+), that were found in 23%, 15.04%, 13.27%, 12.38%, 11.5% and 11.5% of total MIRA subjects, respectively (Figure 4A, marked with arrows). The distribution of the clonal expansion level in samples belonging to the two cohorts was also calculated, not unveiling any statistically significant differential expansion between the two cohorts (Figure 4B, based on Mann-Whitney; the statistical test could not be performed for some of the TCRs, denoted with p-val N/A). For every related sample, the number of times each clonotype appears in the corresponding sample(s) was divided by the total number of the sample’s clonotypes.
Calculation of cohort distribution took place in each clonotype’s subgroup compared to the whole sample with Fisher’s exact test (Figure 4C). We observed a significant difference in the case of the most frequent clonotype CASSIRSSYEQYF+V19-01+J02-07 (p-value = 0.0019) and the fourth most common clonotype CASSLSSPQETQYF+V27-01+J02-05 (p-value = 0.038). Specifically, CASSIRSSYEQYF+V19-01+J02-07 clonotype is detected in a sample’s subpopulation where most subjects are characterised as healthy. In contrast, CASSLSSPQETQYF+V27-01+J02-05 clonotype is detected only in convalescent subjects.
Additionally, the SARS-CoV-2 antigen targets of the six most common clonotypes were identified. CASSIRSSYEQYF+TCRBV19-01+TCRBJ02-07 interacts with surface glycoprotein and ORF1ab, CASSLAGAYEQYF+TCRBV05-01+TCRBJ02-07 recognises the nucleocapsid phosphoprotein, CASSLSAPQETQYF+TCRBV27-01+TCRBJ02-05 recognises ORF1ab and envelope, and CASSLSSPQETQYF+TCRBV27-01+TCRBJ02-05, CASSDRGPNQPQHF+TCRBV27-01+TCRBJ01-05 and CASSDRGPTDTQYF+TCRBV27-01+TCRBJ02-03 interact with ORF1ab. Public databases such as McPAS, TCR3d and VDJdb were used to query the putatively cross-reactive components of the aforementioned TCRs. CASSIRSSYEQYF+TCRBV19-01+TCRBJ02-07 (catalogued as CD8+ in the MIRA dataset and the other TCR catalogues as well) was also found to interact with Influenza virus and Epstein-Barr Virus (EBV). The description of this specific clonotype as part of the immune response against Influenza virus has been previously described [51]. The most common clonotype, which seems to be cross-reactive according to our results (CASSIRSSYEQYF+V19-01+J02-07), is significantly associated with healthy individuals (p-value = 0.0019). The exact opposite phenomenon takes place in CASSLSSPQETQYF+V27-01+J02-05 analysis, where we did not detect any cross-reactivity and a significant correlation was observed with the COVID-19-convalescent cohort (p-value = 0.038), shaping the hypothesis of diversification in SARS-CoV-2 immunity dependent on past infections and/or vaccination.
These results verified the suspicions of cross-reactivity that surfaced through observations based on results of the initial analysis presented in Figure 2C. Figure 4D highlights the most common TCR’s recognition sites on surface glycoprotein from SARS-CoV-2 and matrix protein 1 (M1) from Influenza A virus. The same information is also depicted in the form of a circular plot in Figure 4E. The remaining five most common clonotypes (CASSLAGAYEQYF+TCRBV05-01+TCRBJ02-07, CASSLSAPQETQYF+TCRBV27-01+TCRBJ02-05, CASSLSSPQETQYF+TCRBV27-01+ TCRBJ02-05, CASSDRGPNQPQHF+TCRBV27-01+TCRBJ01-05 and CASSDRGPTDTQYF+ TCRBV27-01+TCRBJ02-03) were not found in VDJdb, McPAS or TCR3d to interact with other pathogens.
To further assess the cross-reactive properties of all TCRs in the MIRA dataset, the epitopes in MIRA as well as epitopes in the three aforementioned public TCR databases were used to match MIRA CDR3 sequences with antigens from SARS-CoV-2 and other pathogens and diseases (Figure 5, Tables S2–S5). In general, some CDR3 sequences can recognise epitopes from multiple antigens. Therefore, the number of SARS-CoV-2 antigen connections with other pathogens’ and diseases’ antigens might not be equal to the number of unique MIRA CDR3 sequences. In Figure 5, the connections between antigens from SARS-CoV-2 and other pathogens (Table S2) and diseases have not been divided according to the CDR3 sequence origin (CD8+ or CD4+ T cell), since we have observed that some sequences are associated with CD8+ T cells in MIRA and CD4+ T cells in McPAS, VDJdb and TCR3d databases, or vice versa. The same analysis was repeated after grouping CDR3 sequences based on their CD8+ or CD4+ T cell origin in both MIRA and external databases (Figure S1, Tables S3 and S4).
As shown in Figure 5A and Figure S1, evidence regarding the cross-reactivity phenomenon is widespread and links SARS-CoV-2 to a plethora of pathogens and other diseases that can be grouped into three major categories (Table S6). The first category includes M. tuberculosis and viruses such as Influenza virus, Cytomegalovirus (CMV), EBV, Human Immunodeficiency Virus (HIV), Hepatitis C Virus (HCV), Yellow Fever Virus (YFV), Dengue Virus (DENV) and Human T-lymphotropic Virus type 1 (HTLV-1). From 2136 CDR3 sequence connections between the MIRA and other repositories, 1792 (83.9%) can recognise epitopes from SARS-CoV-2 and members of the first category (Tables S2 and S6). The second category consists of malignancies and malignancy-related agents such as Melanoma, Breast Cancer and Neoantigens. Roughly 9.9% (211 out of 2136) of the CDR3-mediated connections between MIRA and the three public TCR profiling databases recognise epitopes from the second category (Tables S2 and S6). On the other hand, 6.2% (133 out of 2136) of the connections recognise epitopes from the third category, which reflects auto-immune states and disorders arising from external stimuli such as Celiac Disease, Inflammatory Bowel Disease (IBD), Diabetes Type 1, Psoriatic Arthritis, Allergy and Toxic Epidermal Necrolysis (Tables S2 and S6).
The most cross-reactive partner of SARS-CoV-2 was found to be M. tuberculosis, with 747 CDR3-mediated connections recognising epitopes from both pathogens (Figure 5A, Table S2). These CDR3 sequences were isolated mostly from CD8+ T cell TCRs in MIRA and CD4+ in the other TCR repositories (Table S5). Specifically, 271 connections are associated with ORF1ab and 127 with surface glycoprotein. Influenza virus exhibits the second highest number of CDR3-mediated connections (498 out of 2136) with SARS-CoV-2 (Table S2). The majority of CDR3 sequences are associated with CD8+ T cell TCRs in both MIRA and the public TCR profiling databases (Table S5), in contrast to the case of M. tuberculosis. Notably, 144 connections are related to surface glycoprotein and 123 to ORF1ab.
To have a complete view of the cross-reactivity phenomenon, we further generated a circular plot depicting the “cross-talk” between distinct antigen regions through the recognition by MIRA CDR3 sequences (Figure 5B, Table S5). To ensure concise visualisation, a subset with the most cross-reactive pathogens, as depicted in Figure 5A, was selected for generating the plot, including M. tuberculosis, Influenza virus (A subtype), CMV, EBV and HIV. Each scaled segment in the circle represents an antigen and every antigen is coloured based on the pathogen it belongs to. The links connecting antigen pairs correspond to the “cross-talk” between the antigens’ segments through their ability to be recognised by a specific CDR3 in MIRA. Although M. tuberculosis is the most cross-reactive partner of SARS-CoV-2, the matching epitope information is missing from McPAS, TCR3d and VDJdb for the majority of CDR3 sequences. Thus, the M. tuberculosis connections in Figure 5B are severely limited. The same phenomenon was also observed for the other pathogens, but to a lesser degree.

4. Discussion

Over the last two years, the global impact of COVID-19 on healthcare [52] and the socio-economic [53] field has been devastating. The response of both the scientific community and pharmaceutical industry was swift and decisive in exploring the biological aspect of SARS-CoV-2 and its pathological implications, as well as delivering pharmaceutical products that could assist in restraining COVID-19. During the past year, an observed shift in the literature was evident, highlighting the T cell immunological profile characterisation in the framework of COVID-19 progression and severity [9,10,11,12,13,14,15]. Collaborations between academia and industry resulted in the publication of immunological datasets from studies with thousands of subjects, such as the immunoACCESS© resource [15]. Specifically, the MIRA dataset provides access to thousands of TCR clonotypes that can specifically recognise SARS-CoV-2 epitopes (Figure 1A).
Our strategy is based on a novel ML-oriented TCR profiling assay analytic approach, which can highlight the targeted antigens of the immune response against pathogens and diseases (Figure 1B and Figure 3). During the last two decades we have experienced an abundance of breakthroughs in biotechnology that facilitated the dawn of the big data era for biology. We believe the scientific community should emphasise the development of efficient and accurate computational approaches, to exploit the wealth of information embedded in the ever-increasing volume of biomedical data. ML can be the ideal substrate for combining data from heterogeneous sources of information, while unveiling higher-order and more abstract connections between the underlying mechanisms of biological phenomena and the environment. In the context of COVID-19 related research and other infectious diseases, ML could be used for combining epidemiological surveillance with data from immunoassays (TCR-Seq and others), genomics, transcriptomics and even metagenomics. Such approaches can provide a solid foundation for understanding the entanglement of genetic factors and the environment, as well as their implications for the progression of pandemics, for example, within different populations.
One key observation in our study is that T cell TCRs in the MIRA dataset with the ability to recognise epitopes from ORF1ab and surface glycoprotein exhibit similar levels of clonal expansion between the two cohorts and present the highest clonal expansion levels in the healthy cohort (Figure 2C); this suggests either previously unreported COVID-19 exposure or putative cross-reactivity between SARS-CoV-2 and other pathogens. The former hypothesis could not be verified by any means, since the samples’ collection date in MIRA is not available. Therefore, we proceeded exploring whether MIRA TCRs exhibit cross-reactive properties that can be a product of an immune response against both SARS-CoV-2 and other pathogens and diseases.
Our analysis was based on the CDR3 sequences that are common between MIRA and McPAS, TCR3d and VDJdb repositories. The results unveiled widespread evidence of putative cross-reactivity that link SARS-CoV-2 to a plethora of pathogens and other diseases (Figure 5), which can be grouped into three major categories: (a) M. tuberculosis and viruses including Influenza virus, CMV, EBV and HIV among others, (b) malignancies and malignancy-related agents, and (c) auto-immune states and disorders. Interestingly, the majority of CDR3 sequences that target pathogens in the first category originate from CD8+ T Cell TCRs, according to McPAS, TCR3d and VDJdb, with M. tuberculosis being the exception. The association of BCG vaccine with CD4+ T cell response against M. tuberculosis has been previously reported in the literature [54]. In contrast to the first category, we observed the exact opposite pattern for CDR3 sequences that are associated with pathological states from the second and third categories, since they are derived mostly from CD4+ T cell TCRs.
M. tuberculosis was found to exhibit the highest levels of T cell cross-reactivity with SARS-CoV-2. These results are in accordance with published epidemiological studies conducted prior to SARS-CoV-2 vaccine implementation, which have suggested a negative association between incidence, morbidity and mortality of COVID-19 and national Bacille Calmette–Guérin (BCG) vaccination programs. Specifically, in countries where national BCG vaccination has been implemented, lower numbers COVID-19 cases and related deaths have been recorded [55,56,57]. A study conducted by Escobar et al. [58] found a 10.4% reduction in mortality from COVID-19 for every 10% increase in a country’s BCG index. Discovered in 1921, even today the function of BCG vaccine remains obscure to a certain extent. The BCG vaccine contains attenuated Mycobacterium bovis, and induces humoral and adaptive immunity, activating both non-specific and cross-reactive immune responses in the host [59,60] against a variety of infectious (viruses, bacteria, fungi and parasites) and non-infectious agents. Epigenetic and metabolic reprogramming of innate immune cells, known as trained immunity, is considered responsible for these protective effects [61,62,63,64].
Influenza virus exhibits the second highest number of CDR3 sequences that are putatively cross-reactive with SARS-CoV-2. This type of cross-reactivity could be attributed to the seasonal vaccination against Influenza virus and relevant exposure of a large portion of the population to this virus. The role of protective immunity induced by the polyvalent influenza virus vaccine (against Influenza A virus and/or Influenza B virus subtypes) and the likelihood of COVID-19 has been previously examined [65,66,67]; meanwhile, others explored this association from a clinical manifestation and disease outcome perspective [68]. Although the precise pathophysiological mechanisms underlying this association require further investigation, three main theories have been put forward. The first theory relates to antigenic mimicry which results in clonal activation and lymphocytes proliferation [69]. Depending on each individual’s HLAs, only a limited number of epitopes can be recognised, and those are the immunodominant ones. A second theory of trained immunity has also been proposed as the mechanism behind these beneficial heterologous effects of vaccines [70]. Influenza virus vaccination acts as a non-specific exciter of our immune response [71]. Debisarun et al. found that rather than binding, cytokines broaden T cell responses against SARS-CoV-2 [72]. Salem et al. suggested flu-induced bystander immune response as a probable protective mechanism [73].
Information related to T cell cross-reactivity between SARS-CoV-2 and the remaining viruses from the first category of cross-reactive pathogens is scarce in the literature. Cellular cross-reactivity against EBV and SARS-CoV-2 has not been studied to date. However, there have been reports associating the clinical manifestation of COVID-19 with reactivation of EBV infection and correlating it with severe disease progression, thus underpinning a possible entanglement [74,75,76]. Cross-reactivity between SARS-CoV-2 and HIV has also been reported in studies describing false positive HIV results in COVID-19 patients. Such cases were associated with antibody cross-reaction during immunoassay screening tests, while cross-reactive CDR3 regions were detected between the two viruses [77,78]. Sequence analysis has shown that HIV and SARS-CoV proteins share common motifs that shape a degree of homology [79]. In addition, studying similar antigenic features could help in engineering super antibodies that neutralise different pathogens [80].
Evidence in the literature regarding the connection between SARS-CoV-2 and cancer at various levels remains extremely limited. Most related studies examine the immune response to SARS-CoV-2 in patients with cancer [81]. However, the cross-reactivity phenomenon between SARS-CoV-2 and malignancy-related antigens has not been previously reported. Conversely, since the COVID-19 pandemic was declared in 2020, numerous reports in the literature have linked the immune response against SARS-CoV-2 proteins with self-antigens, thus unveiling putative COVID-19 implications for autoimmune disorders including immune thrombocytopenic purpura, Guillian–Barrė syndrome and subvariants, antiphospholipid antibodies and lupus anticoagulant, Kawasaki and multisystem inflammatory syndrome in children [82,83,84,85]. In general, the emergence of autoimmunity after viral infections involving EBV, CMV, HTLV-1, herpes and hepatitis virus among others has been thoroughly described in the literature [86].
Interestingly, in the MIRA dataset there are some CD4+ and CD8+ T cell TCRs that have common CDR3 sequences. The same observation was made when comparing CDR3 sequences between MIRA and other TCR repositories. There are numerous TCRs that share the same CDR3 sequence, which originate from CD8+ T cells in MIRA and CD4+ T cells in other databases. However, the functional relevance of CD4+ and CD8+ T cells with common CDR3 sequences remains obscure, since to the best of our knowledge this phenomenon has not been debated in the literature.
The observations in this study are of multilevel significance, ranging from putative indirect protection from severe COVID-19 infection based on vaccination against and/or previous exposure to Influenza viruses, M. tuberculosis and other pathogens, to the association of SARS-CoV-2 related TCRs with malignancies and autoimmune disorders. However, there are several limitations related to this study. The MIRA dataset includes a very limited number of samples associated with COVID-19-acute, COVID-19-non-acute and COVID-19-exposed cohorts, thus prohibiting any statistical or ML analysis that potentially connects TCR profile irregularities with disease severity. We did however manage to make statistically significant assumptions using samples derived from COVID-19-convalescent and healthy cohorts, although ideally the number of these samples should be higher. Due to the limited number of MIRA samples, the generalisation capacity of the ML models presented here would possibly be limited on other SARS-CoV-2 TCR-Seq datasets. However, these models are serving a specific purpose in this study, which is to use them for extracting meaningful biological information related to the antigens that can distinguish MIRA cohorts. Another limitation relates to the availability of HLA allele information connecting TCR clonotypes and epitopes in McPAS, TCR3d and VDJdb. The cross-reactivity analysis focused only on the CDR3 amino acid sequence comparison between the MIRA dataset and other databases, due to the limited availability of the relevant HLA information for most TCR clonotypes in the external repositories. Therefore, it should be noted that the cross-reactions described here could only take place in individuals with specific HLA alleles, enabling the presentation of relative epitopes to potential long-lived memory T cells developed during previous infection and/or vaccination. Additionally, the extent of cross-reactivity is delimited by the inherent data bias in McPAS, TCR3d and VDJdb, stemming from the scientific community’s focus on specific pathogenic cases.

5. Conclusions

We believe that our study provides a novel ML-based computational framework for analysing TCR-Seq datasets and systematically highlights the breadth and depth of “cross-talk” between antigens from different pathogens, a phenomenon that may also exhibit therapeutic implications, especially in the COVID-19 context. It is our view that the scientific community should accelerate the effort of generating TCR profiling data without omitting the relevant HLA information, and include assays based on as many human pathologies as possible. ML can play a pivotal role in combining such data with multiomics and epidemiological surveillance, to build intelligent infrastructures that can be an invaluable asset in fully understanding the underlying immune response complexity.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biology11101531/s1, Figure S1: Cross-reactivity analysis of all Multiplexed Identification of T cell Receptor Antigen (MIRA) T cell receptors (TCRs) only from CD8+ T cells; Table S1: PCA loadings; Table S2: Unique cross-reactive counts of all MIRA T cell CDR3 sequences; Table S3: Unique cross-reactive counts of MIRA CD8+ T cell CDR3 sequences; Table S4: Unique cross-reactive counts of MIRA CD4+ T cell CDR3 sequences; Table S5: All MIRA cross-reactive CDR3 sequences and the corresponding recognised epitopes from MIRA and external databases; Table S6: Unique cross-reactive counts of all MIRA T cell CDR3 sequences and the corresponding categories of cross-reactive pathogens.

Author Contributions

G.K.G. and A.P.G. designed the study under M.S.’s and C.H.’s supervision. G.K.G. performed the Machine Learning analysis, cross-reactivity plots and prepared the figures. A.P.G. performed the MIRA and external databases TCR/CDR3 data filtering and the cross-reactivity analysis. G.K.G. and A.P.G. wrote the paper with the assistance of Z.T., M.K., V.A.M., M.S. and C.H. Study supervision was conducted by M.S. and C.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The MIRA dataset was retrieved from the immuneACCESS© database [15]. Additional TCR CDR3 sequences from both CD4+ and CD8+ T cells were downloaded from McPAS [48] (version 4 January 2022), TCR3d [49] (version 13 January 2022) and VDJdb [50] (version 22 March 2022). All protein sequences were downloaded from UniProt [87]. The cross-reactivity exploration was achieved with custom Python scripts and Circos [88]. The three-dimensional structure of M1 and surface glycoprotein antigenic molecules was generated with Jmol, within Jalview software based on 1AA7 (A and B chain view) and 6 × 29 (A chain view) Protein Data Bank [89] entries for M1 and surface glycoprotein, respectively. The antigenic epitopes were aligned on reference sequences using the Clustal algorithm. The statistical, dimensionality reduction, ML and feature importance analyses were performed with in-house developed software based on Python’s scipy and scikit-learn as well as R’s rstudioapi, dplyr, plyr, GLDEX, TSDT, stats, stringr and ggplot2 libraries. The code and processed data used in this study can be found in gitlab.com/hyimmera/mira_analysis (accessed on 23 September 2022). Further inquiries can be directed to the corresponding authors.

Acknowledgments

We wish to acknowledge Lemonia Anagnostopoulos for proofreading the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. JHCRC. John Hopkins Coronavirus Resource Center; Johns Hopkins University School of Medicine: Baltimore, MD, USA. Available online: https://coronavirus.jhu.edu/covid-19-daily-video (accessed on 23 September 2022).
  2. Eastin, C.; Eastin, T. Clinical Characteristics of Coronavirus Disease 2019 in China. J. Emerg. Med. 2020, 58, 711–712. [Google Scholar] [CrossRef]
  3. Velavan, T.P.; Meyer, C.G. The COVID-19 Epidemic. Trop. Med. Int. Health 2020, 25, 278–280. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Lopez-Leon, S.; Wegman-Ostrosky, T.; Perelman, C.; Sepulveda, R.; Rebolledo, P.A.; Cuapio, A.; Villapol, S. More than 50 Long-Term Effects of COVID-19: A Systematic Review and Meta-Analysis. Sci. Rep. 2021, 11, 16144. [Google Scholar] [CrossRef] [PubMed]
  5. Whitley, R. Molnupiravir—A Step toward Orally Bioavailable Therapies for COVID-19. N. Engl. J. Med. 2021, 386, 592–593. [Google Scholar] [CrossRef] [PubMed]
  6. Gupta, A.; Gonzalez-Rojas, Y.; Juarez, E.; Crespo Casal, M.; Moya, J.; Falci, D.R.; Sarkis, E.; Solis, J.; Zheng, H.; Scott, N.; et al. Early Treatment for COVID-19 with SARS-CoV-2 Neutralizing Antibody Sotrovimab. N. Engl. J. Med. 2021, 385, 1941–1950. [Google Scholar] [CrossRef]
  7. Gottlieb, R.L.; Vaca, C.E.; Paredes, R.; Mera, J.; Webb, B.J.; Perez, G.; Oguchi, G.; Ryan, P.; Nielsen, B.U.; Brown, M.; et al. Early Remdesivir to Prevent Progression to Severe COVID-19 in Outpatients. N. Engl. J. Med. 2022, 386, 305–315. [Google Scholar] [CrossRef]
  8. Tregoning, J.S.; Flight, K.E.; Higham, S.L.; Wang, Z.; Pierce, B.F. Progress of the COVID-19 Vaccine Effort: Viruses, Vaccines and Variants versus Efficacy, Effectiveness and Escape. Nat. Rev. Immunol. 2021, 21, 626–636. [Google Scholar] [CrossRef] [PubMed]
  9. Minervina, A.A.; Komech, E.A.; Titov, A.; Bensouda Koraichi, M.; Rosati, E.; Mamedov, I.Z.; Franke, A.; Efimov, G.A.; Chudakov, D.M.; Mora, T.; et al. Longitudinal High-Throughput TCR Repertoire Profiling Reveals the Dynamics of T-Cell Memory Formation after Mild COVID-19 Infection. Elife 2021, 10, e63502. [Google Scholar] [CrossRef]
  10. Hanna, S.J.; Codd, A.S.; Gea-Mallorqui, E.; Scourfield, D.O.; Richter, F.C.; Ladell, K.; Borsa, M.; Compeer, E.B.; Moon, O.R.; Galloway, S.A.E.; et al. T Cell Phenotypes in COVID-19—A Living Review. Oxf. Open Immunol. 2020, 2, iqaa007. [Google Scholar] [CrossRef]
  11. Shomuradova, A.S.; Vagida, M.S.; Sheetikov, S.A.; Zornikova, K.V.; Kiryukhin, D.; Titov, A.; Peshkova, I.O.; Khmelevskaya, A.; Dianov, D.V.; Malasheva, M.; et al. SARS-CoV-2 Epitopes Are Recognized by a Public and Diverse Repertoire of Human T Cell Receptors. Immunity 2020, 53, 1245–1257.e5. [Google Scholar] [CrossRef]
  12. Chang, C.-M.; Feng, P.-H.; Wu, T.-H.; Alachkar, H.; Lee, K.-Y.; Chang, W.-C. Profiling of T Cell Repertoire in SARS-CoV-2-Infected COVID-19 Patients Between Mild Disease and Pneumonia. J. Clin. Immunol. 2021, 41, 1131–1145. [Google Scholar] [CrossRef] [PubMed]
  13. Li, L.; Chen, Q.; Han, X.; Shen, M.; Hu, C.; Chen, S.; Zhang, J.; Wang, Y.; Li, T.; Huang, J.; et al. T Cell Immunity Evaluation and Immunodominant Epitope T Cell Receptor Identification of Severe Acute Respiratory Syndrome Coronavirus 2 Spike Glycoprotein in COVID-19 Convalescent Patients. Front. Cell Dev. Biol. 2021, 9, 696662. [Google Scholar] [CrossRef] [PubMed]
  14. Wang, P.; Jin, X.; Zhou, W.; Luo, M.; Xu, Z.; Xu, C.; Li, Y.; Ma, K.; Cao, H.; Huang, Y.; et al. Comprehensive Analysis of TCR Repertoire in COVID-19 Using Single Cell Sequencing. Genomics 2021, 113, 456–462. [Google Scholar] [CrossRef] [PubMed]
  15. Nolan, S.; Vignali, M.; Klinger, M.; Dines, J.N.; Kaplan, I.M.; Svejnoha, E.; Craft, T.; Boland, K.; Pesesky, M.; Gittelman, R.M.; et al. A Large-Scale Database of T-Cell Receptor Beta (TCRβ) Sequences and Binding Associations from Natural and Synthetic Exposure to SARS-CoV-2. Res. Sq. 2020, 1–28. [Google Scholar] [CrossRef]
  16. Gittelman, R.M.; Lavezzo, E.; Snyder, T.M.; Zahid, H.J.; Elyanow, R.; Dalai, S.; Kirsch, I.; Baldo, L.; Manuto, L.; Franchin, E.; et al. Diagnosis and Tracking of Past SARS-CoV-2 Infection in a Large Study of Vo’, Italy through T-cell Receptor Sequeancing. medRxiv 2020, 9, 2020. [Google Scholar]
  17. Channappanavar, R.; Fett, C.; Zhao, J.; Meyerholz, D.K.; Perlman, S. Virus-Specific Memory CD8 T Cells Provide Substantial Protection from Lethal Severe Acute Respiratory Syndrome Coronavirus Infection. J. Virol. 2014, 88, 11034–11044. [Google Scholar] [CrossRef] [Green Version]
  18. Zhao, J.; Alshukairi, A.N.; Baharoon, S.A.; Ahmed, W.A.; Bokhari, A.A.; Nehdi, A.M.; Layqah, L.A.; Alghamdi, M.G.; Al Gethamy, M.M.; Dada, A.M.; et al. Recovery from the Middle East Respiratory Syndrome Is Associated with Antibody and T-Cell Responses. Sci. Immunol. 2017, 2, eaan5393. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Gallais, F.; Velay, A.; Wendling, M.J.; Nazon, C.; Partisani, M.; Sibilia, J.; Candon, S.; Fafi-Kremer, S. Intrafamilial Exposure to SARS-CoV-2 Associated with Cellular Immune Response without Seroconversion, France. Emerg. Infect. Dis. 2021, 27, 113–121. [Google Scholar] [CrossRef]
  20. Thieme, C.; Anft, M.; Paniskaki, K.; Blázquez Navarro, A.; Doevelaar, A.; Seibert, F.S.; Hölzer, B.; Konik, M.J.; Brenner, T.; Tempfer, C.; et al. The SARS-CoV-2 T-Cell Immunity Is Directed Against the Spike, Membrane, and Nucleocapsid Protein and Associated with COVID 19 Severity. medRxiv 2020. [Google Scholar] [CrossRef]
  21. Pinto, D.; Park, Y.J.; Beltramello, M.; Walls, A.C.; Tortorici, M.A.; Bianchi, S.; Jaconi, S.; Culap, K.; Zatta, F.; De Marco, A.; et al. Cross-Neutralization of SARS-CoV-2 by a Human Monoclonal SARS-CoV Antibody. Nature 2020, 583, 290–295. [Google Scholar] [CrossRef]
  22. Rydyznski Moderbacher, C.; Ramirez, S.I.; Dan, J.M.; Grifoni, A.; Hastie, K.M.; Weiskopf, D.; Belanger, S.; Abbott, R.K.; Kim, C.; Choi, J.; et al. Antigen-Specific Adaptive Immunity to SARS-CoV-2 in Acute COVID-19 and Associations with Age and Disease Severity. Cell 2020, 183, 996–1012. [Google Scholar] [CrossRef]
  23. Sekine, T.; Perez-Potti, A.; Rivera-Ballesteros, O.; Strålin, K.; Gorin, J.B.; Olsson, A.; Llewellyn-Lacey, S.; Kamal, H.; Bogdanovic, G.; Muschiol, S.; et al. Robust T Cell Immunity in Convalescent Individuals with Asymptomatic or Mild COVID-19. Cell 2020, 183, 158–168.e14. [Google Scholar] [CrossRef] [PubMed]
  24. Le Bert, N.; Tan, A.T.; Kunasegaran, K.; Tham, C.Y.L.; Hafezi, M.; Chia, A.; Chng, M.H.Y.; Lin, M.; Tan, N.; Linster, M.; et al. SARS-CoV-2-Specific T Cell Immunity in Cases of COVID-19 and SARS, and Uninfected Controls. Nature 2020, 584, 457–462. [Google Scholar] [CrossRef] [PubMed]
  25. Mateus, J.; Grifoni, A.; Tarke, A.; Sidney, J.; Ramirez, S.I.; Dan, J.M.; Burger, Z.C.; Rawlings, S.A.; Smith, D.M.; Phillips, E.; et al. Selective and Cross-Reactive SARS-CoV-2 T Cell Epitopes in Unexposed Humans. Science 2020, 370, 89–94. [Google Scholar] [CrossRef] [PubMed]
  26. Sette, A.; Crotty, S. Pre-Existing Immunity to SARS-CoV-2: The Knowns and Unknowns. Nat. Rev. Immunol. 2020, 20, 457–458. [Google Scholar] [CrossRef]
  27. Braun, J.; Loyal, L.; Frentsch, M.; Wendisch, D.; Georg, P.; Kurth, F.; Hippenstiel, S.; Dingeldey, M.; Kruse, B.; Fauchere, F.; et al. SARS-CoV-2-Reactive T Cells in Healthy Donors and Patients with COVID-19. Nature 2020, 587, 270–274. [Google Scholar] [CrossRef]
  28. Pennock, N.D.; White, J.T.; Cross, E.W.; Cheney, E.E.; Tamburini, B.A.; Kedl, R.M. T Cell Responses: Naive to Memory and Everything in Between. Adv. Physiol. Educ. 2013, 37, 273–283. [Google Scholar] [CrossRef] [Green Version]
  29. Petrova, G.; Ferrante, A.; Gorski, J. Cross-Reactivity of T Cells and Its Role in the Immune System. Crit. Rev. Immunol. 2012, 32, 349–372. [Google Scholar] [CrossRef] [Green Version]
  30. Welsh, R.M.; Che, J.W.; Brehm, M.A.; Selin, L.K. Heterologous Immunity between Viruses. Immunol. Rev. 2010, 235, 244–266. [Google Scholar] [CrossRef] [Green Version]
  31. Bangs, S.C.; Baban, D.; Cattan, H.J.; Li, C.K.-F.; McMichael, A.J.; Xu, X.-N. Human CD4+ Memory T Cells Are Preferential Targets for Bystander Activation and Apoptosis. J. Immunol. 2009, 182, 1962–1971. [Google Scholar] [CrossRef] [Green Version]
  32. Selin, L.K.; Varga, S.M.; Wong, I.C.; Welsh, R.M. Protective Heterologous Antiviral Immunity and Enhanced Immunopathogenesis Mediated by Memory T Cell Populations. J. Exp. Med. 1998, 188, 1705–1715. [Google Scholar] [CrossRef]
  33. Urbani, S.; Amadei, B.; Fisicaro, P.; Pilli, M.; Missale, G.; Bertoletti, A.; Ferrari, C. Heterologous T Cell Immunity in Severe Hepatitis C Virus Infection. J. Exp. Med. 2005, 201, 675–680. [Google Scholar] [CrossRef] [PubMed]
  34. Sharma, S.; Thomas, P.G. The Two Faces of Heterologous Immunity: Protection or Immunopathology. J. Leukoc. Biol. 2014, 95, 405–416. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Macdonald, W.A.; Chen, Z.; Gras, S.; Archbold, J.K.; Tynan, F.E.; Clements, C.S.; Bharadwaj, M.; Kjer-Nielsen, L.; Saunders, P.M.; Wilce, M.C.J.; et al. T Cell Allorecognition via Molecular Mimicry. Immunity 2009, 31, 897–908. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Wooldridge, L.; Ekeruche-Makinde, J.; van den Berg, H.A.; Skowera, A.; Miles, J.J.; Tan, M.P.; Dolton, G.; Clement, M.; Llewellyn-Lacey, S.; Price, D.A.; et al. A Single Autoimmune T Cell Receptor Recognizes More than a Million Different Peptides. J. Biol. Chem. 2012, 287, 1168–1177. [Google Scholar] [CrossRef] [Green Version]
  37. Christen, U.; Edelmann, K.H.; McGavern, D.B.; Wolfe, T.; Coon, B.; Teague, M.K.; Miller, S.D.; Oldstone, M.B.A.; von Herrath, M.G. A Viral Epitope That Mimics a Self Antigen Can Accelerate but Not Initiate Autoimmune Diabetes. J. Clin. Investig. 2004, 114, 1290–1298. [Google Scholar] [CrossRef] [Green Version]
  38. Reiser, J.-B.; Darnault, C.; Grégoire, C.; Mosser, T.; Mazza, G.; Kearney, A.; van der Merwe, P.A.; Fontecilla-Camps, J.C.; Housset, D.; Malissen, B. CDR3 Loop Flexibility Contributes to the Degeneracy of TCR Recognition. Nat. Immunol. 2003, 4, 241–247. [Google Scholar] [CrossRef]
  39. Ding, Y.H.; Baker, B.M.; Garboczi, D.N.; Biddison, W.E.; Wiley, D.C. Four A6-TCR/Peptide/HLA-A2 Structures That Generate Very Different T Cell Signals Are Nearly Identical. Immunity 1999, 11, 45–56. [Google Scholar] [CrossRef] [Green Version]
  40. Borbulevych, O.Y.; Piepenbrink, K.H.; Gloor, B.E.; Scott, D.R.; Sommese, R.F.; Cole, D.K.; Sewell, A.K.; Baker, B.M. T Cell Receptor Cross-Reactivity Directed by Antigen-Dependent Tuning of peptide-MHC Molecular Flexibility. Immunity 2009, 31, 885–896. [Google Scholar] [CrossRef] [Green Version]
  41. Cornberg, M.; Clute, S.C.; Watkin, L.B.; Saccoccio, F.M.; Kim, S.-K.; Naumov, Y.N.; Brehm, M.A.; Aslan, N.; Welsh, R.M.; Selin, L.K. CD8 T Cell Cross-Reactivity Networks Mediate Heterologous Immunity in Human EBV and Murine Vaccinia Virus Infections. J. Immunol. 2010, 184, 2825–2838. [Google Scholar] [CrossRef] [Green Version]
  42. Haanen, J.B.A.G.; Wolkers, M.C.; Kruisbeek, A.M.; Schumacher, T.N.M. Selective Expansion of Cross-Reactive Cd8+ Memory T Cells by Viral Variants. J. Exp. Med. 1999, 190, 1319–1328. [Google Scholar] [CrossRef] [PubMed]
  43. Clute, S.C.; Naumov, Y.N.; Watkin, L.B.; Aslan, N.; Sullivan, J.L.; Thorley-Lawson, D.A.; Luzuriaga, K.; Welsh, R.M.; Puzone, R.; Celada, F.; et al. Broad Cross-Reactive TCR Repertoires Recognizing Dissimilar Epstein-Barr and Influenza A Virus Epitopes. J. Immunol. 2010, 185, 6753–6764. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Sidhom, J.-W.; Baras, A.S. Deep Learning Identifies Antigenic Determinants of Severe SARS-CoV-2 Infection within T-Cell Repertoires. Sci. Rep. 2021, 11, 14275. [Google Scholar] [CrossRef] [PubMed]
  45. Shoukat, M.S.; Foers, A.D.; Woodmansey, S.; Evans, S.C.; Fowler, A.; Soilleux, E.J. Use of Machine Learning to Identify a T Cell Response to SARS-CoV-2. Cell Rep. Med. 2021, 2, 100192. [Google Scholar] [CrossRef] [PubMed]
  46. Klinger, M.; Pepin, F.; Wilkins, J.; Asbury, T.; Wittkop, T.; Zheng, J.; Moorhead, M.; Faham, M. Multiplex Identification of Antigen-Specific T Cell Receptors Using a Combination of Immune Assays and Immune Receptor Sequencing. PLoS ONE 2015, 10, e0141561. [Google Scholar] [CrossRef]
  47. Lefranc, M.-P.; Giudicelli, V.; Duroux, P.; Jabado-Michaloud, J.; Folch, G.; Aouinti, S.; Carillon, E.; Duvergey, H.; Houles, A.; Paysan-Lafosse, T.; et al. IMGT®, the international ImMunoGeneTics information system® 25 years on. Nucleic Acids Res. 2014, 43, D413–D422. [Google Scholar] [CrossRef] [Green Version]
  48. Tickotsky, N.; Sagiv, T.; Prilusky, J.; Shifrut, E.; Friedman, N. McPAS-TCR: A Manually Curated Catalogue of Pathology-Associated T Cell Receptor Sequences. Bioinformatics 2017, 33, 2924–2929. [Google Scholar] [CrossRef] [Green Version]
  49. Gowthaman, R.; Pierce, B.G. TCR3d: The T Cell Receptor Structural Repertoire Database. Bioinformatics 2019, 35, 5323–5325. [Google Scholar] [CrossRef] [Green Version]
  50. Bagaev, D.V.; Vroomans, R.M.A.; Samir, J.; Stervbo, U.; Rius, C.; Dolton, G.; Greenshields-Watson, A.; Attaf, M.; Egorov, E.S.; Zvyagin, I.V.; et al. VDJdb in 2019: Database Extension, New Analysis Infrastructure and a T-Cell Receptor Motif Compendium. Nucleic Acids Res. 2020, 48, D1057–D1062. [Google Scholar] [CrossRef]
  51. Sant, S.; Grzelak, L.; Wang, Z.; Pizzolla, A.; Koutsakos, M.; Crowe, J.; Loudovaris, T.; Mannering, S.I.; Westall, G.P.; Wakim, L.M.; et al. Single-Cell Approach to Influenza-Specific CD8+ T Cell Receptor Repertoires Across Different Age Groups, Tissues, and Following Influenza Virus Infection. Front. Immunol. 2018, 9, 1453. [Google Scholar] [CrossRef] [Green Version]
  52. Kaye, A.D.; Okeagu, C.N.; Pham, A.D.; Silva, R.A.; Hurley, J.J.; Arron, B.L.; Sarfraz, N.; Lee, H.N.; Ghali, G.E.; Gamble, J.W.; et al. Economic Impact of COVID-19 Pandemic on Healthcare Facilities and Systems: International Perspectives. Best Pract. Res. Clin. Anaesthesiol. 2021, 35, 293–306. [Google Scholar] [CrossRef] [PubMed]
  53. Nundy, S.; Ghosh, A.; Mesloub, A.; Albaqawy, G.A.; Alnaim, M.M. Impact of COVID-19 Pandemic on Socio-Economic, Energy-Environment and Transport Sector Globally and Sustainable Development Goal (SDG). J. Clean. Prod. 2021, 312, 127705. [Google Scholar] [CrossRef]
  54. Jasenosky, L.D.; Scriba, T.J.; Hanekom, W.A.; Goldfeld, A.E. T Cells and Adaptive Immunity to Mycobacterium Tuberculosis in Humans. Immunol. Rev. 2015, 264, 74–87. [Google Scholar] [CrossRef] [PubMed]
  55. Miller, A.; Reandelar, M.J.; Fasciglione, K.; Roumenova, V.; Li, Y.; Otazu, G.H. Correlation between Universal BCG Vaccination Policy and Reduced Mortality for COVID-19. medRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  56. Berg, M.K.; Yu, Q.; Salvador, C.E.; Melani, I.; Kitayama, S. Mandated Bacillus Calmette-Guérin (BCG) Vaccination Predicts Flattened Curves for the Spread of COVID-19. Sci. Adv. 2020, 6, eabc1463. [Google Scholar] [CrossRef]
  57. Charoenlap, S.; Piromsopa, K.; Charoenlap, C. Potential Role of Bacillus Calmette-Guérin (BCG) Vaccination in COVID-19 Pandemic Mortality: Epidemiological and Immunological Aspects. Asian Pac. J. Allergy Immunol. 2020, 38, 150–161. [Google Scholar]
  58. Escobar, L.E.; Molina-Cruz, A.; Barillas-Mury, C. BCG Vaccine Protection from Severe Coronavirus Disease 2019 (COVID-19). Proc. Natl. Acad. Sci. USA 2020, 117, 17720–17726. [Google Scholar] [CrossRef]
  59. Moorlag, S.J.C.F.M.; Arts, R.J.W.; van Crevel, R.; Netea, M.G. Non-Specific Effects of BCG Vaccine on Viral Infections. Clin. Microbiol. Infect. 2019, 25, 1473–1478. [Google Scholar] [CrossRef]
  60. Uthayakumar, D.; Paris, S.; Chapat, L.; Freyburger, L.; Poulet, H.; De Luca, K. Non-Specific Effects of Vaccines Illustrated Through the BCG Example: From Observations to Demonstrations. Front. Immunol. 2018, 9, 2869. [Google Scholar] [CrossRef] [Green Version]
  61. Saeed, S.; Quintin, J.; Kerstens, H.H.D.; Rao, N.A.; Aghajanirefah, A.; Matarese, F.; Cheng, S.-C.; Ratter, J.; Berentsen, K.; van der Ent, M.A.; et al. Epigenetic Programming of Monocyte-to-Macrophage Differentiation and Trained Innate Immunity. Science 2014, 345, 1251086. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Mehta, S.; Jeffrey, K.L. Beyond Receptors and Signaling: Epigenetic Factors in the Regulation of Innate Immunity. Immunol. Cell Biol. 2015, 93, 233–244. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Arts, R.J.W.; Moorlag, S.J.C.F.M.; Novakovic, B.; Li, Y.; Wang, S.-Y.; Oosting, M.; Kumar, V.; Xavier, R.J.; Wijmenga, C.; Joosten, L.A.B.; et al. BCG Vaccination Protects against Experimental Viral Infection in Humans through the Induction of Cytokines Associated with Trained Immunity. Cell Host Microbe 2018, 23, 89–100.e5. [Google Scholar] [CrossRef] [PubMed]
  64. Moulson, A.J.; Av-Gay, Y. BCG Immunomodulation: From the “hygiene Hypothesis” to COVID-19. Immunobiology 2021, 226, 152052. [Google Scholar] [CrossRef] [PubMed]
  65. Jehi, L.; Ji, X.; Milinovich, A.; Erzurum, S.; Rubin, B.P.; Gordon, S.; Young, J.B.; Kattan, M.W. Individualizing Risk Prediction for Positive Coronavirus Disease 2019 Testing: Results From 11,672 Patients. Chest 2020, 158, 1364–1375. [Google Scholar] [CrossRef] [PubMed]
  66. Noale, M.; Trevisan, C.; Maggi, S.; Antonelli Incalzi, R.; Pedone, C.; Di Bari, M.; Adorni, F.; Jesuthasan, N.; Sojic, A.; Galli, M.; et al. The Association between Influenza and Pneumococcal Vaccinations and SARS-CoV-2 Infection: Data from the EPICOVID19 Web-Based Survey. Vaccines 2020, 8, 471. [Google Scholar] [CrossRef]
  67. Pawlowski, C.; Puranik, A.; Bandi, H.; Venkatakrishnan, A.J.; Agarwal, V.; Kennedy, R.; O’Horo, J.C.; Gores, G.J.; Williams, A.W.; Halamka, J.; et al. Exploratory Analysis of Immunization Records Highlights Decreased SARS-CoV-2 Rates in Individuals with Recent non-COVID-19 Vaccinations. Sci. Rep. 2021, 11, 4741. [Google Scholar] [CrossRef]
  68. Fink, G.; Orlova-Fink, N.; Schindler, T.; Grisi, S.; Ferrer, A.P.; Daubenberger, C.; Brentani, A. Inactivated Trivalent Influenza Vaccination Is Associated with Lower Mortality among COVID-19 Patients in Brazil. BMJ Evid. Based Med. 2021, 26, 192–193. [Google Scholar] [CrossRef]
  69. Cohen, I.R. Antigenic Mimicry, Clonal Selection and Autoimmunity. J. Autoimmun. 2001, 16, 337–340. [Google Scholar] [CrossRef] [Green Version]
  70. Netea, M.G.; Domínguez-Andrés, J.; Barreiro, L.B.; Chavakis, T.; Divangahi, M.; Fuchs, E.; Joosten, L.A.B.; van der Meer, J.W.M.; Mhlanga, M.M.; Mulder, W.J.M.; et al. Defining Trained Immunity and Its Role in Health and Disease. Nat. Rev. Immunol. 2020, 20, 375–388. [Google Scholar] [CrossRef] [Green Version]
  71. Eldanasory, O.A.; Rabaan, A.A.; Al-Tawfiq, J.A. Can Influenza Vaccine Modify COVID-19 Clinical Course? Travel Med. Infect. Dis. 2020, 37, 101872. [Google Scholar] [CrossRef]
  72. Pallikkuth, S.; Williams, E.; Pahwa, R.; Hoffer, M.; Pahwa, S. Association of Flu Specific and SARS-CoV-2 Specific CD4 T Cell Responses in SARS-CoV-2 Infected Asymptomatic Heath Care Workers. Vaccine 2021, 39, 6019–6024. [Google Scholar] [CrossRef] [PubMed]
  73. Salem, M.L.; El-Hennawy, D. The Possible Beneficial Adjuvant Effect of Influenza Vaccine to Minimize the Severity of COVID-19. Med. Hypotheses 2020, 140, 109752. [Google Scholar] [CrossRef] [PubMed]
  74. Lehner, G.F.; Klein, S.J.; Zoller, H.; Peer, A.; Bellmann, R.; Joannidis, M. Correlation of Interleukin-6 with Epstein--Barr Virus Levels in COVID-19. Crit. Care 2020, 24, 657. [Google Scholar] [CrossRef] [PubMed]
  75. Nadeem, A.; Suresh, K.; Awais, H.; Waseem, S. Epstein-Barr Virus Coinfection in COVID-19. J. Investig. Med. High Impact Case Rep. 2021, 9, 23247096211040624. [Google Scholar] [CrossRef]
  76. Vigón, L.; García-Pérez, J.; Rodríguez-Mora, S.; Torres, M.; Mateos, E.; de la Osa, M.; Cervero, M.; Malo De Molina, R.; Navarro, C.; Murciano-Antón, M.A.; et al. Impaired Antibody-Dependent Cellular Cytotoxicity in a Spanish Cohort of Patients With COVID-19 Admitted to the ICU. Front. Immunol. 2021, 12, 742631. [Google Scholar] [CrossRef] [PubMed]
  77. Tan, S.S.; Chew, K.L.; Saw, S.; Jureen, R.; Sethi, S. Cross-Reactivity of SARS-CoV-2 with HIV Chemiluminescent Assay Leading to False-Positive Results. J. Clin. Pathol. 2021, 74, 614. [Google Scholar] [CrossRef]
  78. Salih, R.Q.; Salih, G.A.; Abdulla, B.A.; Ahmed, A.D.; Mohammed, H.R.; Kakamad, F.H.; Salih, A.M. False-Positive HIV in a Patient with SARS-CoV-2 Infection; a Case Report. Ann. Med. Surg. 2021, 71, 103027. [Google Scholar] [CrossRef]
  79. Kliger, Y.; Levanon, E.Y. Cloaked Similarity between HIV-1 and SARS-CoV Suggests an anti-SARS Strategy. BMC Microbiol. 2003, 3, 20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Mishra, N.; Kumar, S.; Singh, S.; Bansal, T.; Jain, N.; Saluja, S.; Kumar, R.; Bhattacharyya, S.; Palanichamy, J.K.; Mir, R.A.; et al. Cross-Neutralization of SARS-CoV-2 by HIV-1 Specific Broadly Neutralizing Antibodies and Polyclonal Plasma. PLoS Pathog. 2021, 17, e1009958. [Google Scholar] [CrossRef]
  81. Fendler, A.; Au, L.; Shepherd, S.T.C.; Byrne, F.; Cerrone, M.; Boos, L.A.; Rzeniewicz, K.; Gordon, W.; Shum, B.; Gerard, C.L.; et al. Functional Antibody and T Cell Immunity Following SARS-CoV-2 Infection, Including by Variants of Concern, in Patients with Cancer: The CAPTURE Study. Nat. Cancer 2021, 2, 1321–1337. [Google Scholar] [CrossRef]
  82. Mehandru, S.; Merad, M. Pathological Sequelae of Long-Haul COVID. Nat. Immunol. 2022, 23, 194–202. [Google Scholar] [CrossRef] [PubMed]
  83. Ehrenfeld, M.; Tincani, A.; Andreoli, L.; Cattalini, M.; Greenbaum, A.; Kanduc, D.; Alijotas-Reig, J.; Zinserling, V.; Semenova, N.; Amital, H.; et al. COVID-19 and Autoimmunity. Autoimmun. Rev. 2020, 19, 102597. [Google Scholar] [CrossRef] [PubMed]
  84. Galeotti, C.; Bayry, J. Autoimmune and Inflammatory Diseases Following COVID-19. Nat. Rev. Rheumatol. 2020, 16, 413–414. [Google Scholar] [CrossRef]
  85. Dotan, A.; Muller, S.; Kanduc, D.; David, P.; Halpert, G.; Shoenfeld, Y. The SARS-CoV-2 as an Instrumental Trigger of Autoimmunity. Autoimmun. Rev. 2021, 20, 102792. [Google Scholar] [CrossRef] [PubMed]
  86. Barzilai, O.; Ram, M.; Shoenfeld, Y. Viral Infection Can Induce the Production of Autoantibodies. Curr. Opin. Rheumatol. 2007, 19, 636–643. [Google Scholar] [CrossRef] [PubMed]
  87. UniProt Consortium. UniProt: The Universal Protein Knowledge Base in 2021. Nucleic Acids Res. 2021, 49, D480–D489. [Google Scholar] [CrossRef]
  88. Krzywinski, M.; Schein, J.; Birol, I.; Connors, J.; Gascoyne, R.; Horsman, D.; Jones, S.J.; Marra, M.A. Circos: An Information Aesthetic for Comparative Genomics. Genome Res. 2009, 19, 1639–1645. [Google Scholar] [CrossRef] [Green Version]
  89. Berman, H.M.; Westbrook, J.; Feng, Z.; Gilliland, G.; Bhat, T.N.; Weissig, H.; Shindyalov, I.N.; Bourne, P.E. The Protein Data Bank. Nucleic Acids Res. 2000, 28, 235–242. [Google Scholar] [CrossRef]
Figure 1. Overview of this study. (A) Outline of the Multiplexed Identification of T cell Receptor Antigen (MIRA) assay and the corresponding dataset available from the immunoACCESS© project web resource. (B) Analytic steps in this study, regarding the novel utilisation of the MIRA dataset for training Machine Learning algorithms that can highlight important SARS-CoV-2 antigens for distinguishing samples between healthy and COVID-19-convalescent cohorts. (C) Strategy for exploring T cell cross-reactivity between SARS-CoV-2 and other pathogens and diseases.
Figure 1. Overview of this study. (A) Outline of the Multiplexed Identification of T cell Receptor Antigen (MIRA) assay and the corresponding dataset available from the immunoACCESS© project web resource. (B) Analytic steps in this study, regarding the novel utilisation of the MIRA dataset for training Machine Learning algorithms that can highlight important SARS-CoV-2 antigens for distinguishing samples between healthy and COVID-19-convalescent cohorts. (C) Strategy for exploring T cell cross-reactivity between SARS-CoV-2 and other pathogens and diseases.
Biology 11 01531 g001
Figure 2. Exploratory analysis of the Multiplexed Identification of T cell Receptor Antigen (MIRA) dataset. (A) Number of samples in each MIRA cohort. (B) Per sample normalised number of unique T cell receptors (TCRs) in the healthy and convalescent cohorts. (C) Per sample normalised number of TCRs that recognise each SARS-CoV-2 antigen in the healthy and convalescent cohorts. (D) Projection of healthy and convalescent samples on the principal component analysis (PCA) space. Healthy and convalescent distributions in (B,C) were compared with the Mann-Whitney test.
Figure 2. Exploratory analysis of the Multiplexed Identification of T cell Receptor Antigen (MIRA) dataset. (A) Number of samples in each MIRA cohort. (B) Per sample normalised number of unique T cell receptors (TCRs) in the healthy and convalescent cohorts. (C) Per sample normalised number of TCRs that recognise each SARS-CoV-2 antigen in the healthy and convalescent cohorts. (D) Projection of healthy and convalescent samples on the principal component analysis (PCA) space. Healthy and convalescent distributions in (B,C) were compared with the Mann-Whitney test.
Biology 11 01531 g002
Figure 3. Evaluation of Machine Learning (ML) algorithms trained on the healthy and convalescent cohorts in the Multiplexed Identification of T cell Receptor Antigen (MIRA) dataset. (A) Balanced accuracy, precision, sensitivity, specificity and negative predictive value (NPV) of each algorithm after selecting a prediction score cut-off of 0.5. (B) Support Vector Machines (SVM) performance on multiple prediction score cut-offs. (C) Feature importance score after 50 permutations on all 20 randomly generated test sets. (D) ML algorithms’ performance after selecting only the important features for each algorithm and retraining.
Figure 3. Evaluation of Machine Learning (ML) algorithms trained on the healthy and convalescent cohorts in the Multiplexed Identification of T cell Receptor Antigen (MIRA) dataset. (A) Balanced accuracy, precision, sensitivity, specificity and negative predictive value (NPV) of each algorithm after selecting a prediction score cut-off of 0.5. (B) Support Vector Machines (SVM) performance on multiple prediction score cut-offs. (C) Feature importance score after 50 permutations on all 20 randomly generated test sets. (D) ML algorithms’ performance after selecting only the important features for each algorithm and retraining.
Biology 11 01531 g003
Figure 4. Exploration of the most common Multiplexed Identification of T cell Receptor Antigen (MIRA) T cell receptors (TCRs) in terms of clonal expansion and cross-reactivity. (A) Occurrence frequency and clonal expansion of all MIRA TCRs. The arrows point to the six most common TCRs (present in at least 11.5% of total number of subjects) that were further analysed in terms of clonal expansion in the two cohorts based on Mann-Whitney test (B). The statistical test could not be performed for some TCRs (denoted as p-val N/A). (C) Cohort distribution of the first and fourth most common MIRA TCRs that were found to be enriched in either cohort after applying Fisher’s exact test. (D) Secondary structure of surface glycoprotein (SARS-CoV-2) and Matrix protein 1 (M1) with cross-reactive sections, based on the most common MIRA TCR, highlighted with red colour. Locations highlighted with red colour consist of epitopes recognised by the cross-reactive TCRs and putatively reflect protein domains with similar structural or physicochemical properties. (Ε) Circular plot, as an alternative view of (D), depicting the cross-reactive property of the most common MIRA TCR that recognises epitopes from surface glycoprotein and M1. The inner and outer light-colored tracks represent the annotated domains.
Figure 4. Exploration of the most common Multiplexed Identification of T cell Receptor Antigen (MIRA) T cell receptors (TCRs) in terms of clonal expansion and cross-reactivity. (A) Occurrence frequency and clonal expansion of all MIRA TCRs. The arrows point to the six most common TCRs (present in at least 11.5% of total number of subjects) that were further analysed in terms of clonal expansion in the two cohorts based on Mann-Whitney test (B). The statistical test could not be performed for some TCRs (denoted as p-val N/A). (C) Cohort distribution of the first and fourth most common MIRA TCRs that were found to be enriched in either cohort after applying Fisher’s exact test. (D) Secondary structure of surface glycoprotein (SARS-CoV-2) and Matrix protein 1 (M1) with cross-reactive sections, based on the most common MIRA TCR, highlighted with red colour. Locations highlighted with red colour consist of epitopes recognised by the cross-reactive TCRs and putatively reflect protein domains with similar structural or physicochemical properties. (Ε) Circular plot, as an alternative view of (D), depicting the cross-reactive property of the most common MIRA TCR that recognises epitopes from surface glycoprotein and M1. The inner and outer light-colored tracks represent the annotated domains.
Biology 11 01531 g004
Figure 5. Cross-reactivity analysis of all Multiplexed Identification of T cell Receptor Antigen (MIRA) T cell receptors (TCRs). (A) Heatmap of unique MIRA complementarity-determining region 3 (CDR3) counts that exhibit cross-reactivity between SARS-CoV-2 (x-axis) and other pathogens and diseases (y-axis). The heatmap values correspond to the number of unique cross-reactive CDR3 sequences. (B) Circular plot that depicts the cross-reactivity of MIRA CDR3 regions between antigens that originate from SARS-CoV-2 and a selected subset of pathogens from (A). The inner and outer light-coloured tracks represent the annotated protein domains. Each connection represents the ability of a single CDR3 region to recognise a part of a SARS-CoV-2 antigen and a part of another pathogen’s protein. The connections are coloured based on their corresponding non-SARS-CoV-2 pathogens.
Figure 5. Cross-reactivity analysis of all Multiplexed Identification of T cell Receptor Antigen (MIRA) T cell receptors (TCRs). (A) Heatmap of unique MIRA complementarity-determining region 3 (CDR3) counts that exhibit cross-reactivity between SARS-CoV-2 (x-axis) and other pathogens and diseases (y-axis). The heatmap values correspond to the number of unique cross-reactive CDR3 sequences. (B) Circular plot that depicts the cross-reactivity of MIRA CDR3 regions between antigens that originate from SARS-CoV-2 and a selected subset of pathogens from (A). The inner and outer light-coloured tracks represent the annotated protein domains. Each connection represents the ability of a single CDR3 region to recognise a part of a SARS-CoV-2 antigen and a part of another pathogen’s protein. The connections are coloured based on their corresponding non-SARS-CoV-2 pathogens.
Biology 11 01531 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Georgakilas, G.K.; Galanopoulos, A.P.; Tsinaris, Z.; Kyritsi, M.; Mouchtouri, V.A.; Speletas, M.; Hadjichristodoulou, C. Machine-Learning-Assisted Analysis of TCR Profiling Data Unveils Cross-Reactivity between SARS-CoV-2 and a Wide Spectrum of Pathogens and Other Diseases. Biology 2022, 11, 1531. https://doi.org/10.3390/biology11101531

AMA Style

Georgakilas GK, Galanopoulos AP, Tsinaris Z, Kyritsi M, Mouchtouri VA, Speletas M, Hadjichristodoulou C. Machine-Learning-Assisted Analysis of TCR Profiling Data Unveils Cross-Reactivity between SARS-CoV-2 and a Wide Spectrum of Pathogens and Other Diseases. Biology. 2022; 11(10):1531. https://doi.org/10.3390/biology11101531

Chicago/Turabian Style

Georgakilas, Georgios K., Achilleas P. Galanopoulos, Zafeiris Tsinaris, Maria Kyritsi, Varvara A. Mouchtouri, Matthaios Speletas, and Christos Hadjichristodoulou. 2022. "Machine-Learning-Assisted Analysis of TCR Profiling Data Unveils Cross-Reactivity between SARS-CoV-2 and a Wide Spectrum of Pathogens and Other Diseases" Biology 11, no. 10: 1531. https://doi.org/10.3390/biology11101531

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop