Profiling of T Cell Receptor β-Chain Complimentary Determining Regions 3 Repertoire in Subarachnoid Hemorrhage Patients Using High-Throughput Sequencing

Article information

J Korean Neurosurg Soc. 2021;64(4):505-513

Publication date (electronic) : 2021 June 29

doi : https://doi.org/10.3340/jkns.2020.0214

Bong Jun Kim^1,*, Jun Hyong Ahn^2,*, Dong Hyuk Youn¹, Jin Pyeong Jeon^,1,2,3

¹Institute of New Frontier Stroke Research, Hallym University College of Medicine, Chuncheon, Korea

²Department of Neurosurgery, Hallym University College of Medicine, Chuncheon, Korea

³Genetic and Research Inc., Chuncheon, Korea

Address for reprints : Jin Pyeong Jeon Department of Neurosurgery, Hallym University College of Medicine, 77 Sakju-ro, Chuncheon 24253, Korea Tel : +82-33-240-5171, Fax : +82-33-240-9970, E-mail : jjs6553@daum.net

*Bong Jun Kim and Jun Hyong Ahn contributed equally to this work and should be considered co-first authors.

Received 2020 July 27; Revised 2020 September 18; Accepted 2020 October 7.

Abstract

INTRODUCTION

The adaptive immune response is mediated by diverse receptors on the surface of T cells that allow for recognition of pathogens [8]. Accordingly, identification of the T cell receptor (TCR) repertoire can be used to monitor the host immune response. Characterization of TCR repertoires have proven useful in cancer and autoimmune disease patients who received immunotherapy [34]. Sims et al. [30] sequenced TCRs of glioma αand β-chains. They showed that more divergent VJ-independent repertoires are associated with tumors than non-neoplastic tissues. Yu et al. [34] evaluated the associations between TCR repertoire and symptom remission in systemic lupus erythematosus (SLE) patients treated with glucocorticoids. They suggested that monitoring of the TCR complimentary determining region (CDR) 3 is related to prognosis in SLE, suggesting that it mirrors therapeutic efficacy. T cells are also actively involved in atherosclerosis [14,22]. Jia et al. [14] reported epigenetic suppression of forkhead box that was related with a decrease in regulatory T cells, resulting in an increased risk of atherosclerosis. In addition, the T cell response is involved in the initiation of plaque disruption [10].

SAH occurs spontaneously via ruptured aneurysm, arteriovenous malformation (AVM), or traumatic brain injury. The mortality rate of SAH due to ruptured cerebral aneurysm exceeded 20% [21,33]. Although an aneurysm itself may be successfully occluded, neurologic complications can occur due to direct brain injury and vascular dysregulation [26]; such complications are associated with poor neurologic outcomes if not promptly treated. Consequently, SAH patients require multimodal monitoring in the neurointensive care unit for early detection of neurologic complications and subsequent optimal intervention [26]. In attempts to understand the pathophysiology of SAH, many researchers have focused on genomic and proteomic analyses; by comparison, the adaptive immune response has not been well studied [1]. Recently, we reported the TCR repertoire for the first time in SAH patients with and without delayed cerebral ischemia (DCI) [18]. In addition, we further evaluated the T cell immunology based on SAH severity to delineate SAH pathogenesis.

MATERIALS AND METHODS

RESULTS

Baseline characteristics of the study

TCRB CDR3 repertoires from six SAH patients were analyzed and compared at two different time points, admission and 7-day follow-up. Three patients presented with goodgrade SAH, while the three others had poor-grade SAH. The median age was 53.5 years (range, 41–77). Four of the patients were female (66.7%). Most aneurysms (n=5, 83.3%) were located in the anterior circulation. Two patients experienced severe DCI and recovered well after chemical angioplasty. Three patients with good-grade SAH recovered well without significant disability. However, one out of three patients (33.3%) exhibited poor neurologic outcome defined by mRS score of 3 (Table 1).

Table 1.

Clinical characteristics of the enrolled patients

Sequence, clonality, and CDR3 lengths of TCRB

Variables such as total and unique sequences, clonality, and frequency of the TCRB CDR3 repertoire were measured (Table 2). Clonality using the Shannon index at admission and follow-up were 0.059 (0.037–0.038) and 0.027 (0.014–0.082) (p=0.699), respectively. Clonality of poor-grade (0.025 [0.011–0.038]) SAH was significantly lower than that of good-grade SAH (0.095 [0.079–0.101]; p=0.026; Fig. 1A).

Table 2.

TCRB CDR3 repertoire – total sequences, unique sequences, and clonality – in patients with subarachnoid hemorrhage at two times, admission and follow-up

Fig. 1.

Comparison of TCRB clonality (A) and CDR3 length repertoire (B) according to SAH severity, good- vs. poor-grade. Good-grade SAH patients had higher clonality and shorter CDR3 length than did poor-grade SAH patients. CDR : complimentary determining regions, TCRB : T cell receptor β-chain, SAH : subarachnoid hemorrhage.

Read frequency of ≤0.005% to <0.05% dominated the clonal expansion, but the changes between admission and follow-up did not show a definite association (Supplementary Table 1). We ranked TCRB clonotype frequencies and compared repertoire diversity using the Simpson diversity index (D) [5], where a score close to one indicates high diversity (Supplementary Fig. 1). Overall, the diversity indices ranged from 0.993 to 0.999 at admission and from 0.987 to 0.999 at follow-up. Poor-grade SAH patients had diversity scores of 0.998 and 0.999, which were higher than the scores of good-grade SAH patients.

CDR3 lengths ranged from 15 to 78 nucleotides, with a peak at 42 in good-grade and at 45 in poor-grade SAH patients. The five most common CDR3 lengths were 36, 39, 42, 45, and 48 nucleotides. The results are calculated as the average of five individuals. Good-grade SAH cases tended to have shorter CDR lengths compared to poor-grade SAH (Fig. 1B). The ten most expanded sequences were presented and compared in terms of amino acid levels (Fig. 2) [32]. The distribution of the top ten most common sequences was different in each patient. Overall, poor-grade SAH patients had a lower proportion of the top ten most common sequences than did good-grade SAH patients.

Fig. 2.

Quantification of TCRB CDR3 sequences in each patient between the day of admission and the 7-day follow-up. The 10 most common sequences are presented in different colors. Patients who presented with poor-grade SAH tended to show a low proportion of the top 10 most common sequences compared to those with good-grade SAH at admission. Adm : admission, F/U : follow-up, TCRB : T cell receptor β-chain, CDR : complimentary determining regions, SAH : subarachnoid hemorrhage.

Distribution of TCRBV and TCRBJ

We analyzed the distribution of the unique clonotypes in both the TCRBV and TCRBJ gene segments. Overall, TCRBV gene segments from the follow-up day had different clonotypes than did those from the admission day, but the distribution of TCRBJ gene segments was similar between days. Compared to good-grade SAH, poor-grade SAH tended to exhibit a great number of different clonotypes of the TCRBV gene segments (Fig. 3). We analyzed the relative frequencies of V-J combinations using different colors representing disease severity. We obtained 754 annotated V-J pairs. TCRBV21-01/TCRBJ01-05, TCRBV21-01/TCRBJ01-04, and TCRBV21-01/TCRBJ02-04 were the most commonly increased V-J pairs in poor-grade SAH compared to good-grade SAH patients (Supplementary Fig. 2). We further compared the V-J combinations at two different time points (admission and 7-day follow-up) in good-grade and poor-grade SAH patients, respectively. The results showed that TCRBV19-01/TCRBJ02-04 was the most increased V-J pair and TCRBV28-01/TCRBJ02-04 was the most decreased V-J pair in good-grade SAH. Regarding poorgrade SAH patients, the most increased and decreased V-J pairs were TCRBV28-01/TCRBJ02-06 and TCRBV30-01/TCRBJ02-04, respectively (Fig. 4).

Fig. 3.

Distribution and change in the unique clonotypes of the TCRBV and TCRBJ genes in patients with SAH. Compared to TCRBJ, TCRBV showed different distributions between the day of admission and the 7-day follow-up. Different colors are used to represent gene frequencies in each patient. TCRBJ : T cell receptor β chain J segment, TCRBV : T cell receptor β chain V segment, Adm : admission, F/U : follow-up, SAH : subarachnoid hemorrhage.

Fig. 4.

Comparison of average TCRB V-J gene utilization at admission and during the 7-day follow-up according to SAH severity of good (A) and poor (B)-grade. V gene and J gene segments are arranged on the x- and y-axis. Different colors indicate the differences in gene frequencies of poor-grade SAH patients compared with good-grade SAH patients. TCRB : T cell receptor β-chain, SAH : subarachnoid hemorrhage.

DISCUSSION

T cells recognize antigens with the major histocompatibility complex-TCR complex. The TCR repertoire mainly relies on diversity in the α- and β-chains [8]. Among the hypervariable regions of CDRs chains, CDR3 has been targeted for high-throughput immune sequencing due to its diversity [8,32]. Chang et al. [5] compared TCRB repertoires in rheumatoid arthritis patients who received different medications. In their study, the diversity of the TCRB repertoire varied according to different biological medications and was inversely correlated with disease activity [5]. Sun et al. [31] reported that TCRB diversity was lower overall in neonatal sepsis than in controls, with different clonotypes and V-J gene combinations in the two. In addition, shortening of CDR3 was suggested to increase autoimmune disease risk by increased self-recognition during thymic T cell development in type 1 diabetus mellitus [9]. Therefore, the TCR repertoire can function as a mirror of the health of the host, a marker of disease severity, or a marker of biological response for immunotherapy [2].

It has been postulated that SAH can spark innate and adaptive immune responses [25]. The study of immune response in SAH has generally focused on pro-inflammatory cytokines, which are responsible for brain edema after hemorrhagic insult, and treatments to alleviate the condition [6]. Hemorrhage leak into the subarachnoid space can trigger recruitment of leukocytes and inflammatory responses [15]. SAH patients exhibited increased levels of interleukin (IL)-6 and soluble IL-2 receptor in the cerebrospinal fluid [24]. Chen et al. [6] reported a beneficial effect of melatonin on preservation of tight junction proteins in the blood-brain barrier, resulting in decreased brain edema in an experimental SAH model in rats. In their study, the levels of pro-inflammatory cytokines such as IL-6, IL-1β, and tumor necrosis factor-α were suppressed after melatonin injection. Moraes et al. [25] also reported higher levels of innate immune cells such as CD14⁺⁺, CD16⁺, and CD69⁺ following SAH. However, the adaptive immune response, and TCRB repertoire in particular, has not been well studied in SAH patients. In this study, we explored the composition and variation of the TCRB CDR3 repertoire to understand T cell-mediated adaptive immunity after SAH. TCRB CDR3 repertoires are dynamic in nature according to disease severity. Poor-grade SAH patients showed significantly decreased clonality and increased diversity as compared with good-grade SAH patients, as well as longer CDR3 amino acid lengths. Interestingly, disease severity was associated with differences in VJ gene usage of TCRB and CDR3 lengths, which can affect TCR rearrangement. Therefore, further studies on the role of individual clonotypes in the development of clinical complications such as DCI or SAH-induced immunocompromised states are necessary.

In this study, we only included SAH patients who received endovascular coil embolization though the transfemoral approach, not those with craniotomy and clipping. The host immune response to tissue injury can affect systemic inflammatory and anti-inflammatory responses [23]. Marik and Flemmer [23] reported increased T-helper 2 cells which can result in impaired cell-mediated immunity after physical injury. Compared to coil embolization, clipping is associated with concerns about infection of the surgical site and tissue injury. Accordingly, we only sequenced and analyzed the TCRB repertoire in SAH patients with coiling to reduce inherent bias from different treatment methods. Therefore, patients who underwent craniotomy and clipping may have different TCRB CDR3 repertoires following SAH, which requires additional investigation.

The study limitations are as follows. First, the small sample size may be a concern, although this is the largest series to date involving SAH patients. Most previous studies enrolled less than 10 patients due to economic concerns [5,12]. Our previous study analyzing the TCRB repertoire in DCI also included five SAH patients [18]. Nevertheless, economic concerns cannot justify the results of the study. Accordingly, further studies are needed to investigate the adaptive immune response in a large number of the SAH patients to validate our findings. Second, we did not evaluate the TCR α-chain. Although most research has been done on either the α- or β-chain, one at a time [8], quantification of both the α- and β-chains can more accurately determine the TCR repertoire. Third, we have focused on the clinical grading system of H-H grade, but not Fisher grade. T cell adaptive immune response via extravasation of blood can trigger additional inflammatory response. We have analyzed and compared the average TCRB V-J gene utilization according to Fisher grades, II and III (A) and IV (B) (Supplementary Fig. 3). TCRBV19-01/TCRBJ02-04 was the most increased V-J pair in Fisher grades II and III in the 7-day follow-up and TCRBV28-01/TCRBJ02-04 was the most decreased V-J pair. In case of Fisher grade IV, TCRBV28-01/TCRBJ02-06 and TCRBV30-01/TCRBJ02-04 were the most increased and decreased V-J pairs, respectively. Nevertheless, the small number of enrolled patients may be a concern for the interpretation, suggesting the need for further TCRB studies focusing on Fisher grade. Fourth, different drug regimens, which possibly affect immunity, cannot accurately reflect the adaptive immune response. Accordingly, additional multicenter studies based on drug regimen are needed to evaluate TCRB repertoire.

CONCLUSION

The TCRB repertoire has a dynamic nature following SAH. Poor-grade SAH patients displayed a reduction in clonality and an increase in diversity and CDR3 length. Monitoring TCRB CDR3 repertoires can provide a better understanding of the adaptive immune response and could be useful for identification of therapeutic targets in SAH patients.

Supplementary Materials

Article information Continued