key: cord-0859523-tiuq9q9o
authors: Assmann, Jorn L.J.C.; Kolijn, P. Martijn; Schrijver, Benjamin; van Gammeren, Adriaan J.; Loth, Daan W.; Ermens, Ton A.A.M.; Dik, Willem A.; van der Velden, Vincent H.J.; Langerak, Anton W.
title: TRB sequences targeting ORF1a/b are associated with disease severity in hospitalized COVID‐19 patients
date: 2021-04-13
journal: J Leukoc Biol
DOI: 10.1002/jlb.6covcra1120-762r
sha: c3a2a09f2e0383e516ce1906f40ef32918f30390
doc_id: 859523
cord_uid: tiuq9q9o

The potential protective or pathogenic role of the adaptive immune response to SARS‐CoV‐2 infection has been vigorously debated. While COVID‐19 patients consistently generate a T lymphocyte response to SARS‐CoV‐2 antigens, evidence of significant immune dysregulation in these patients continues to accumulate. In this study, next generation sequencing of the T cell receptor beta chain (TRB) repertoire was conducted in hospitalized COVID‐19 patients to determine if immunogenetic differences of the TRB repertoire contribute to disease course severity. Clustering of highly similar TRB CDR3 amino acid sequences across COVID‐19 patients yielded 781 shared TRB sequences. The TRB sequences were then filtered for known associations with common diseases such as EBV and CMV. The remaining sequences were cross‐referenced to a publicly accessible dataset that mapped COVID‐19 specific TCRs to the SARS‐CoV‐2 genome. We identified 158 SARS‐CoV‐2 specific TRB sequences belonging to 134 clusters in our COVID‐19 patients. Next, we investigated 113 SARS‐CoV‐2 specific clusters binding only one peptide target in relation to disease course. Distinct skewing of SARS‐CoV‐2 specific TRB sequences toward the nonstructural proteins (NSPs) encoded within ORF1a/b of the SARS‐CoV‐2 genome was observed in clusters associated with critical disease course when compared to COVID‐19 clusters associated with a severe disease course. These data imply that T‐lymphocyte reactivity towards peptides from NSPs of SARS‐CoV‐2 may not constitute an effective adaptive immune response and thus may negatively affect disease severity.

response in COVID-19, alongside the established role of the humoral response. 9, 10, 11, 12 The T-cell receptor beta (TRB) repertoire is tremendously diverse.

V-D-J recombination can theoretically result in the formation of 10 12 TCRs, allowing recognition of essentially all peptide antigens presented by human leukocyte antigen (HLA) molecules. 13 The TRB chain is most important in determining TCR diversity, the CDR3 region in particular is highly variable and heavily involved in binding to peptide HLA complexes. Accordingly, the CDR3 region of TRB sequences can be considered a unique barcode for T-lymphocytes targeting a particular antigen. 14 Upon initiating a cellular immune response, T-lymphocytes undergo clonal expansion, 15 which might be represented through a decreased diversity of the TRB gene repertoire and/or appearance of epitopespecific (e.g., COVID-19 specific) clonotypes. 16 Clonotypes are mostly defined as sequences acquired by usage of the same V-, D-, and Jgene and identical CDR3 region, and can be readily identified by next-generation sequencing (NGS). However, it has been shown that epitope-specific TRBs often share conserved CDR3 features, without being fully identical. To this extent, epitope-specific TRBs often form clusters of clonotypes that share similar sequence features and therefore CDR3 regions. 2 Until now, studies regarding the TRB repertoire in COVID-19 patients were largely descriptive 17, 18 and a detailed analysis of the dynamics of the TRB repertoire in relation to disease course (i.e., severe vs. critical disease) has, to our knowledge, not been performed. We hypothesize that immunogenetic differences in the TRB repertoire may contribute to a severe (hospitalized) vs critical (ICU and/or death) disease course and that SARS-CoV-2 epitope-specific TRBs directed towards peptides from nonstructural proteins might be involved in the defective immune response that might lead to critical disease outcome.

To test this hypothesis, we employed NGS of the TRB repertoire and characterized TRB sequences of interest which have highly similar TRB CDR3 amino acid sequences derived from TRB sequences enhanced in COVID-19 patients. These sequences were clustered using the TEIRE-SIAS algorithm. To determine which of the enhanced sequences in COVID-19 patients were SARS-CoV-2 epitope-specific, we utilized publicly available datasets of TRB-antigen associations. 19, 20 2

We determined the temporal dynamics of the TRB repertoire in a total of 87 COVID-19 patients of whom 46 had a severe and 41 a critical disease course. No significant differences were observed in the cumulative frequency of the top 10 clonotypes (i.e., sum of frequencies of 10 most abundant clonotypes) in severe and critical COVID-19 patients at baseline and longitudinally ( Figure 1A and B). Furthermore, the TRB Shannon diversity score did not differ between the severe and critical outcome groups ( Figure 1C and D). Given the lack of overt differences in TRB repertoire between severe and critical COVID- 19, we postulate that the T-lymphocyte response to COVID-19 may be mediated by several low frequency clonotypes, rather than by a single immuno-dominant clonotype. This idea is supported by other studies, where multiple smaller COVID-19 specific clonotypes were found in individual patients. 20 Formally, we cannot entirely exclude that clonal T-lymphocytes have homed to the lungs and therefore may be underrepresented in the peripheral blood, as has been described in studies exploring paired TCR sequencing of bronchial alveolar lavage fluid and peripheral blood. 21 The key to new insights based on the current dataset would be to establish which of the TRB clonotypes identified are specific for SARS-CoV-2. Public clonotypes, identical clonotypes shared across individuals, are a well-described phenomenon in both healthy and diseased subjects. In fact, mechanisms such as convergent recombination result in a far more limited TCR repertoire as theoretically possible.

T-lymphocyte responses towards the same pathogen are often driven through the same antigens across individuals, resulting in epitopespecific clonotypes (enhanced sequences) which have been described

for several viral diseases, including influenza, EBV, and CMV. 2, 22 Therefore, we (and others) postulated that such enhanced sequences could be found in COVID-19 patients as well. 23 In our cohort, we identified 781 highly similar TRB CDR3 amino acid sequences representing 483 clusters in the 87 patients ( Figure 2A ).

Most of the clusters (443) consisted of TRB sequences shared in 2 or 3 patients, with a minority (41) shared by ≥4 patients (Supplementary Table 1 ). We defined a cluster as "severe" or "critical" if all TRB sequences within that cluster were derived from patients with the same disease course. In total, we observed 96 severe clusters and 117 critical clusters, while the remaining 268 clusters were not associated with a particular outcome (Supplementary Table 1 ). To filter out shared public TRB sequences targeting common diseases such as EBV and CMV, we queried the VDJ database (VDJdb), an online repository for known TCR-antigen interactions. 19 We discovered 54 of the 781 shared TRB sequences had previously been associated with an antigen from a disease other than COVID-19 in the VDJdb, suggesting that the prevalence of these TRB sequences was likely not related to COVID-19, albeit the potential for cross-reactivity cannot be completely excluded ( Figure 2A ). 19 As the remaining 727 shared TRB sequences were not associated with any antigens from previously studied diseases in the VDJdb, the prevalence of these TRB sequences in COVID-19 patients may be the result of stimulation by SARS-CoV-2 peptides. To further explore this, we utilized a dataset from a recently published study available through medRxiv that described the TRB repertoire of over 1500 COVID-19 patients and 3500 controls. 20 The authors stimulated Next, we selected the 113 SARS-CoV-2 specific clusters that exclusively targeted peptides from either structural or nonstructural proteins (NSPs) and screened for potential associations of clusters with disease course. We identified 19 clusters that were exclusively found in patients with a severe disease course, while 23 clusters were exclusively found in patients with a critical disease course (Supplemen-tary Table 2b ). Analysis of these clusters revealed a distinct skewing towards NSPs and accessory proteins of the SARS-CoV-2 genome for clusters associated with a critical disease course, while SARS-CoV-2 specific TRB sequences in clusters associated with severe disease course target peptides from structural proteins such as nucleocapsid phosphoprotein (N) and surface glycoprotein (S) (P = 0.026; Figure 2B ).

Upon further review, the skewing of the clusters associated with a critical disease course towards NSPs and accessory proteins appeared primarily mediated by the NSPs of ORF1a/b of the SARS-CoV-2 genome (P < 0.001) ( Figure 2B; Supplementary Table 2b ).

The observed skewing is interesting as multiple studies regarding SARS-CoV-1 have previously indicated that peptides from N and S proteins can induce strong SARS-CoV-1 specific T-lymphocyte responses. [28] [29] [30] Additionally, the memory response to SARS-CoV-1 is almost exclusively directed towards peptides from N proteins. 31 Furthermore, all candidate SARS-CoV-2 vaccines of major pharmaceutical companies are also primarily based on S proteins. 32 -2 0 Severe outcome F I G U R E 2 Mapping of shared highly similar TRB CDR3 amino acid sequences to the SARS-CoV-2 genome. A Venn diagram of shared sequences between highly similar sequences identified in our cohort (TEIRESIAS), TRB sequences stored in the VDJ database of known antigen-TCR associations (VDJdb) and TRB sequences previously associated with SARS-CoV-2 antigens (MIRA) (A). Clusters associated with a critical disease course ((B), P < 0.001) target peptides from nonstructural proteins (NSPs) significantly (P = 0.026), or NSPs encoded in ORF1a/b significantly (P < 0.0001) more frequently compared to clusters associated with a severe disease course, for a detailed overview of the clusters see SupplementaryTable 2b. Temporal dynamics of SARS-CoV-2-specific TRB sequences directed towards peptides from structural and nonstructural proteins in patients with severe (C) or critical (D) COVID-19 disease course. Severe group n = 46, Critical group n = 41. P-values are calculated using Fisher's Exact Test T-lymphocyte response in convalescent COVID-19 patients, these responses were an order of magnitude lower than responses towards peptides from N and S proteins. 5 Furthermore, a recent publication reported that early induction of SARS-CoV-2-specific T-lymphocytes is present in patients with mild disease and accelerated viral clearance.

The T-lymphocytes targeted primarily the N and S proteins, as well as the accessory proteins of ORF7/8. 35 Altogether, our data indicate that an immune response to the nonstructural proteins of ORF1a/b of SARS-CoV-2 may be associated with an inability to sufficiently combat SARS-CoV-2 infection and consequently lead to a critical disease course.

We investigated the relationship between disease duration and the antigenic targets of the TRB repertoire (structural vs. NSPs of ORF1a/b). In the first 16 days after symptom onset, a consistent low frequency response of < 1% of the total TRB repertoire could be observed in 37 of 46 patients in the severe outcome group, primarily targeted at peptides derived from structural proteins ( Figure 2C ).

For the patients with a critical outcome, a longer follow up was available, both due to closer monitoring in the ICU and due to prolonged hospital stay. We observed that clonotypes targeting peptides from nonstructural proteins of ORF1a/b are already present early during COVID-19 infection in patients with a critical outcome, in sharp contrast to the severe outcome group (Figure 2C and D) . While clonotypes targeting structural antigens could be detected up to 24 days after symptom onset, all clonotypes detected after that point targeted peptides of NSPs of ORF1a/b ( Figure 2D ). Based on these observations, we hypothesize there may be some potential for clonotypes specific for peptides from nonstructural proteins of ORF1a/b in COVID-19 as a marker for critical disease.

Collectively our data illustrate a clear role for the T-lymphocyte response in COVID-19. Epitope-specific highly similar clonotypes can be found in the TRB repertoire of COVID-19 patients with both severe and critical disease course. We observed that SARS-CoV-2 specific TRB sequences in clusters associated with a critical disease course, targeted the nonstructural proteins of SARS-CoV-2 significantly more frequently than clusters associated with a severe disease course, which tended to target structural proteins. Therefore, T-lymphocyte reactivity towards peptides from NSPs of ORF1a/b may be related to a critical disease course. These results may be the first step towards new insights into the mechanisms behind COVID-19 disease severity, although our data need to be validated in larger cohorts that also include non-hospitalized COVID-19 patients with a mild disease course.

Eighty-seven patients with a positive SARS-CoV-2 PCR test and a clin- 

Total cellular DNA was isolated from whole blood and 1000 ng was amplified through multiplex PCR of TRBV-TRBD-TRBJ gene rearrangements following the BIOMED-2 protocol. 36 

Raw FASTQ files were uploaded to the ARResT/Interrogate immunoprofiler for annotation and initial exploratory analysis. 37 The full productive TRB repertoire for each sample was downloaded in .csv format and processed in R (R Core Team, 2020) during further analysis. 38 Shannon diversity index was calculated using the vegan package. 39 A summary of relevant QC metrics is depicted in Supplementary Table 5 .

For more details please refer to our entry in the GEO database under accession number GSE161810.

Clonotypes below 50 reads were excluded from clustering to reduce the influence of technical sequencing errors. Clustering was performed using the TEIRESIAS algorithm, 24 a commonly used algorithm for pattern recognition in biological studies. The TEIRESIAS algorithm was initially developed by the Bioinformatics and Pattern Discovery group at the IBM Computational Biology Center 23 and later adapted for CDR3 amino-acid sequences in chronic lymphocytic leukemia. 24 We in turn adapted the identity and similarity thresholds to 80% and 100% respectively, to accommodate for the use of TCR data and to increase the likelihood that the identified highly similar TCRs indeed target the same epitope. Clustering was performed based on CDR3 amino acid sequence alone, TRBV and TRBJ gene identity was not taken into account for clustering.

The MIRA dataset was obtained from the ImmuneCODE data release. 20 However, in this raw dataset of the MIRA results, 3178 out of 120,703 unique CDR3 amino acid TRB sequences targeted more than one peptide from the COVID-19 genome, with some TRB targeting as many as 15. As we deemed it unlikely that such promiscuous binding was entirely accurate, we applied an additional filtering step to these 3178 TRB sequences. In this filtering step, we discarded all peptide-TRB sequence associations that were reported in only one patient, while keeping the peptide-TRB associations reported in multiple patients, reducing this subset to 2647 TRB sequences with robust antigen associations. While this conservative approach certainly risks filtering out some true-positive hits, we considered it more important to prevent any doubts introduced by the potential false-positives in the MIRA-dataset. The 2647 TRB sequences and their peptide associations were merged back into the original MIRA dataset.

Both the MIRA dataset and our own clusters were screened for known TCR-antigen associations through the VDJdatabase (VDJdb) and matches were tagged.

The CDR3 amino acid sequences for each entry in the MIRA dataset were compared with the CDR3 amino acid sequences identified during clustering. All matches can be found in Supplementary Table   2a . For downstream analysis clusters matching with both the MIRA dataset and an unrelated target (autoantigen, other pathogen) in the VDJdb were filtered. Clusters with multiple targets within the same class (structural proteins vs nonstructural proteins) were grouped and labeled as targeting that class. Clusters not exclusively associated with an outcome or targeting peptides from both nonstructural and structural proteins were filtered.

Frequency of a clonotype indicates the % reads of the total TRB repertoire in the patient sample. Cumulative frequency was calculated by summing the frequency of the 10 largest productive clonotypes in the repertoire. Cumulative frequency of the top 10 clonotypes and Shannon diversity score between outcome groups was compared through

Wilcoxon signed-rank test. For the comparison of peptide targets between outcome groups, Fisher's Exact tests were applied on 2 × 2 contingency tables. All statistical analysis was performed in R (R Core Team, 2020).

Quantifiable predictive features define epitope-specific T cell receptor repertoires

Estimation of the asymptomatic ratio of novel coronavirus infections (COVID-19)

Care for critically ill patients with COVID-19

Targets of T cell responses to SARS-CoV-2 coronavirus in humans with COVID-19 disease and unexposed individuals

Dysregulation of immune response in patients with COVID-19 in Wuhan, China. Clinical Infectious Diseases

Functional exhaustion of antiviral lymphocytes in COVID-19 patients

Extensive longitudinal immune profiling reveals sustained innate immune activaton in COVID-19 patients with unfavorable outcome

Memory T cell responses targeting the SARS coronavirus persist up to 11 years post-infection

T cell-dependent B cell activation. Annual review of immunology

Detection of SARS-CoV-2-specific humoral and cellular immunity in COVID-19 convalescent individuals. Immunity

Immune cell profiling of COVID-19 patients in the recovery stage by single-cell sequencing

The many important facets of T-cell repertoire diversity

How TCRs bind MHCs, peptides, and coreceptors

A natural immunological adjuvant enhances T cell clonal expansion through a CD28-dependent, interleukin (IL)-2-independent mechanism

Age-related appearance of a CMV-specific high-avidity CD8+ T cell clonotype which does not occur in young adults

Next-generation sequencing of T and B cell receptor repertoires from COVID-19 patients showed signatures associated with severity of disease

Longitudinal high-throughput TCR repertoire profiling reveals the dynamics of T-cell memory formation after mild COVID-19 infection

VDJdb: a curated database of T-cell receptor sequences with known antigen specificity

Magnitude and dynamics of the T-cell response to SARS-CoV-2 infection at both individual and population levels. medRxiv

The differential immune responses to COVID-19 in peripheral and lung revealed by single-cell RNA sequencing. medRxiv

Identifying specificity groups in the T cell receptor repertoire

Combinatorial pattern discovery in biological sequences: the TEIRESIAS algorithm

A different ontogenesis for chronic lymphocytic leukemia cases carrying stereotyped antigen receptors: molecular and computational evidence

Multiplex identification of antigenspecific T cell receptors using a combination of immune assays and immune receptor sequencing

Analyzing the Mycobacterium tuberculosis immune response by T-cell receptor clustering with GLIPH2 and genome-wide antigen screening

A large-scale database of T-cell receptor beta (TCRβ) sequences and binding associations from natural and synthetic exposure to SARS-CoV

A DNA vaccine induces SARS coronavirus neutralization and protective immunity in mice

Effects of a SARS-associated coronavirus vaccine in monkeys

Generation and characterization of DNA vaccines targeting the nucleocapsid protein of severe acute respiratory syndrome coronavirus

SARS-CoV-2-specific T cell immunity in cases of COVID-19 and SARS, and uninfected controls

Safety and immunogenicity of the ChAdOx1 nCoV-19 vaccine against SARS-CoV-2: a preliminary report of a phase 1/2, single-blind, randomised controlled trial

Single-shot Ad26 vaccine protects against SARS-CoV-2 in rhesus macaques

COVID-19 vaccine BNT162b1 elicits human antibody and TH 1 T-cell responses

Early induction of functional SARS-CoV-2 specific T cells associates with rapid viral clearance and mild disease in COVID-19 patients

Design and standardization of PCR primers and protocols for detection of clonal immunoglobulin and T-cell receptor gene recombinations in suspect lymphoproliferations: report of the BIOMED-2 Concerted Action BMH4-CT98-3936

ARResT/Interrogate: an interactive immunoprofiler for IG/TR NGS data

R: A language and environment for statistical computing

Vegan: community ecology package. R package version

We gratefully acknowledge Mr. P.J. Hengeveld for constructive discussion. 

The authors declare no conflict of interest.

Shared last author;

The data that support the findings of this study are available from the corresponding author upon reasonable request. TRB NGS data are uploaded in the GEO database and are available under accession number GSE161810.

The study was performed in accordance with the guidelines for sharing of patient data of observational scientific research in case of exceptional health situations, as issued by the Commission on Codes of Conduct of the Foundation Federation of Dutch Medical Scientific Societies (https://www.federa.org/federa-english).

Additional supporting information may be found online in the Supporting Information section at the end of the article.