Skip to main content

V-J combinations of T-cell receptor predict responses to erythropoietin in end-stage renal disease patients



Anemia is common among end-stage renal disease (ESRD) patients who undergone hemodialysis. The total reduction of red blood cell (RBC) count is associated with poor prognosis in these patients. Although erythropoietin (EPO) has been used as an effective treatment for ESRD patients with anemia, a large number of patients still present poor responses to EPO treatment.


We measured T-cell receptor sequencing profiles, including length of complementarity-deteremining region 3 (CDR3), intra- and inter-group (EPO resistant vs. responsive) clonotype diversity, V(D)J usage profiles and V-J combinations from ESRD patients and to investigate the correlation between these features and EPO treatment efficacy.


Our results revealed statistical significance in the top 3 ~ 15 most abundant joint distributions of Vβ/Jβ among the two groups, suggesting the importance of V or J gene utilization in the EPO response of ESRD patients.


In summary, we provided evidence addressing the potential correlation between the immune repertoire and EPO response in ESRD patients.

Trial registration

TMU-JIRB 201309026. Registered 16 October 2013.


End-stage renal disease (ESRD), a final-stage kidney disease, is an important health issue worldwide. Increasing ESRD populations have become a substantial burden on global healthcare expenditures [1]. In addition, the rapid decline in renal function, defined as the glomerular filtration rate (GFR), is associated with an increased risk of cardiovascular and all-cause mortality [2]. According to KDIGO 2012 guidelines, the progression from chronic kidney disease (CKD) to ESRD is characterized by grade 5 CKD requiring renal replacement therapy such as kidney transplantation or hemodialysis. The United States Renal Data System (USRDS) Annual Data Report of 2014 indicated that Mexico (466.5 per million population), Taiwan (449.7 per million population), and the United States (358.7 per million population) are with the highest incidence rate of ESRD in 2012. Besides, Taiwan (2902.1 per million population), Japan (2365.2 per million population), and the United States (1975.5 per million population) are with the highest ESRD prevalence in 2012. In particular in Taiwan, a CKD epidemic area, around 60,000 ~ 70,000 ESRD patients receive dialysis each year, and consume 5.89% of the National Health Insurance budget [3].

The impairment of the renal endocrine role in ESRD patients can lead to anemia [4], and the incidence of secondary anemia increases as the GFR declines. In a healthy condition, renal secretion of erythropoietin (EPO) plays a crucial role in red blood cell (RBC) proliferation and differentiation, and insufficient EPO production in ESRD patients leads to both decreasing RBC production and a 30% ~ 70% shorter RBC lifespan [5]. In 2013, Lin et al. reported a significant decrease in the hemoglobin (Hb) level in CKD patients compared to non-CKD samples in both elderly (≥60 years old) and non-elderly (>60 years old) subjects using a large Taiwanese cohort (n = 3352). As anemia may lead to a significant increase of cardiovascular disease risk and faster renal function decline, managing anemia is therefore crucial for the quality of life of ESRD patients. As a good management of anemia will improve the clinical outcomes of ESRD patients, erythropoiesis-stimulating agents (ESAs) are therefore widely used in ESRD-associated anemia. However, some ESRD patients show poor responsibility to EPO treatment. The resistance to EPO was reported to be associated with adverse prognoses and poorer clinical outcomes [6, 7]. Patients with ESRD have increased levels of several cytokines that negatively regulate erythropoiesis, e.g. interleukin (IL)-1, IL-6, tumor necrosis factor (TNF)-α, and interferon [IFN]-ɣ [8, 9]. The potential mechanisms include the direct inhibition to the growth of erythroid precursor, induction of apoptosis, and disrupting iron metabolism. Besides, it is well known that ESRD patients possess disruption of the T-cell-mediated immune responses [10]. Expansion of the memory T-cell population, impairment of T-cell functions, and diminished regulatory T-cells were observed in patients with renal failure [11, 12]. Zal et al. indicated that clonal expansion of the CD4+/CD28 T-cell subpopulation may contribute to poorer cardiovascular outcomes in ESRD patients who receive hemodialysis [13]. Although the strong associations between the T-cell-mediated immune reactions and ESRD have been characterized, the potential mechanisms that influence T cell activation, cytokines production and EPO responses in ESRD patients remain unclear. Immune repertoire sequencing has recently emerged as a powerful method for profiling of the TCR repertoire [14]. This TCR sequencing method quantifying individual clonotypes by sequencing gene segments on variable (V), diversity (D) and joining (J) regions on the TCR, and further depicting the diversity of TCR repertoires. Using this approach, we hypothesized that TCR diversity and clonality among ESRD patients is the key to link the aberrant cytokine production and eventually contribute to the EPO responses.


Subject recruitment

Peripheral blood mononuclear cells (PBMCs) of samples were collected after obtaining informed consent from subjects at Taipei Medical University Hospital. The samples were collected from patients with written informed consent. The experiment was approved by the institutional review board at the Taipei Medical University and the methods were carried out in accordance with relevant guidelines and regulations. According to the KDOQI Clinical Practice Guidelines, ESRD patients were diagnosed according to eGFR decline of <15 ml/min/1.73m2 over 3 or more months. ESRD patients often exhibit anemia, and recombinant human EPO therapy is required to recover the hemoglobin (Hb) levels to at least 11 g/dL. However, some patients have to use a higher dosage of EPO than others to achieve Hb level targets. In this study, we recruited seven ESRD patients, whose eGFR had declined by <15 ml/min/1.73m2 and who were under regular hemodialysis for 3 months or more, and had also used EPO treatment. Three patients with an erythropoietin resistance index (ERI) of >5 U/kg/week (range, 11.6 ~ 12.3 U/kg/week) were classified into an EPO-nonresponsive group, and 4 EPO responders were defined by an ERI of <5.0 U/kg/week (range, 0.84 ~ 4.81 U/kg/week). This criterion resulted in an ERI grouping the same as in previous studies [15,16,17].

Complementary (c)DNA synthesis and a two-step polymerase chain reaction (PCR)

Whole blood was diluted 1:1 with phosphate-buffered saline (PBS), added gently to the top of Histopaque (Sigma), and centrifuged at 400 ×g for 30 min. The peripheral blood mononuclear cells (PBMCs) were harvested after centrifuging and washing with PBS. Sequentially, the RNeasy Mini Kit (Qiagen) and QIAshredder spin column (Qiagen) were applied to extract the RNA from PBMCs. A reverse-transcription (RT)-PCR was performed with an One-Step RT-PCR kit (Qiagen) plus iRepertoire Primers (iRepertoire) according to the manufacturer’s instructions. cDNA templates were amplified by performing a Multiplex PCR with a Multiplex PCR kit (Qiagen) and Illumina Communal Primers of the HTBI-M type (iRepertoire). After library preparation, reads were obtained using the Illumina MiSeq System (Illumina).

Raw sequencing data preprocessing

Raw paired-end sequencing data were de-multiplexed by sample barcodes to the FASTA format. Then, sample indices were trimmed while removing the adapter (Additional file 1). The paired-end reads were further overlapped. Non-overlapping reads were discarded in a downstream analysis. When overlapping, the base with the higher phred score was taken, thereby improving the overall overlapping read quality. The merged sequencing data were filtered by the following criteria: reads with a length of <50 base pairs (bp); reads with any nucleotide with a phred quality score of <20; and reads with <50% nucleotide bases that met the criterion of a phred quality score of >20. Ultimately, this resulted in overlapping reads with ~350 bp, which was sufficient to cover the entire complementarity-determining region (CDR) 3. In addition, the read length distribution was visualized using the seeFastqPlot function in the systemPipeR package (Additional file 2).

VDJ mapping and clonotype identification

The resulting FASTA files were subjected to CDR3β extraction, V, D, and J regions mapping, clonotype identification of PCR, and sequencing errors with using MiTCR software [18] with the following parameters: ‘species = hs’ and ‘gene = TRB’. A phred quality score of >25 was set as the threshold value for each base within the CDR3 region. In addition, we used the ‘eliminate these errors’ algorithm for the PCR and sequencing error correction, which corrects a maximal number of accumulated errors, while sacrificing 1% ~ 3% of the TCR β diversity, according to our experiments (data not shown). The mappability of each sample was calculated by the number of mapped reads divided by the total number of reads. After clonotype clusterization, clonotype enumeration was conducted by summing the count of reads that coded for the same clonotype sequence. The percentage of clonotypes was calculated by the number of clonotypes divided by the number of nucleotide clones; and the percentage of in-frame clonotypes was calculated by the number of in-frame clonotypes divided by the total number.

TCR β repertoire data analysis

In-frame reads were used for all downstream analyses. The inverse Simpson’s index (1/D) was applied to identify a sample’s clonotype diversity with the vegan package. In addition, to avoid the influence of rare clonotypes, we also assessed the diversity of individual sample repertoires by calculating the proportion of the top 5 most-abundant clonotypes. This parameter is robust for sequencing depth, and therefore allows comparisons between individuals. We also evaluated TRBV and TRBJ gene usage to enable profiling of repertoire manifestations of EPO responsiveness in ESRD patients.

Statistical analysis

R software and Bioconductor were applied for all downstream statistical tests and visualization tasks using in-house R scripts.


Quality assessment of the TCR β repertoire sequencing data in ESRD patients

ESRD samples from seven patients were analyzed. Four of seven samples were defined as EPO-responsive (R) based on the erythropoietin resistance index (ERI) while the remaining three were defined as EPO-nonresponsive (NR) (Table 1). The ERI was calculated by the weekly weighted-adjusted EPO dose (U/kg/week) divided by the hemoglobin level (g/dL). The ERI of patients showed a significant difference between the EPO-responsive and EPO-nonresponsive groups (two-sample Fisher-Pitman permutation p = 0.017), with EPO responders and EPO non-responders having mean ERI values of 2.75 and 11.9, respectively (Fig. 1a). An association analysis of ESRD patients’ clinical features also revealed a significant correlation between the hemoglobin level and EPO responsiveness (two-sample Fisher-Pitman permutation p = 0.048, Table 2 and Fig. 1b).

Table 1 Patient characteristics of end-stage renal disease (ESRD) samples
Fig. 1
figure 1

a The erythropoietin resistance index (ERI) levels in both erythropoietin (EPO) resistant and responsive groups (p = 0.017). b Correlation between end-stage renal disease patients’ clinical features and erythropoietin (EPO) responsiveness. Significant correlations between clinical features and EPO responsive status are denoted by “*”. Y-axis, mean square contingency coefficient

Table 2 Correlation between end-stage renal disease patients’ clinical features and erythropoietin (EPO) responsiveness

To evaluate the analysis pipeline for the repertoire analysis, we exploited the 2 following parameters: (1) overlapping proportion and (2) mappability of reads. In the data-preprocessing step, individual TCR β cDNA sequences were subjected to overlapping. The robustness of overlapping step (Additional file 3) was found by a mean of 83.16% of overlap, and the small variation (standard deviation (SD) = 4.77%). In addition, after VDJ region determination, a mappability of 89.35% ± 4.91% further revealed the reliability of this approach (Additional file 4). The depths of CDR3 amino acid sequences (clonotypes) of individuals were also calculated. Only in-frame clone types, with nucleotide lengths that were divisible by 3, were subjected to all downstream analyses (Additional file 5).

Repertoire diversity and EPO responsiveness in ESRD patients

To assess clonal expansions of TCR β in ESRD patients, we first depicted the β-CDR3 length distribution in seven ESRD patients (Fig. 2a, c and Additional file 6). As shown, substantial oligoclonal enrichment was observed in samples R1, R2, NR1, and NR3. The length of the most frequent β-CDR3 sequences ranging from 13 to 15 amino acids in these patients. In addition, the sequence pattern of the first 4 residues was very similar to PBMC of healthy individuals (Fig. 2b, d) [19].

Fig. 2
figure 2

a Spectratype plot showing the distribution of T cell receptor (TCR) β complementarity determining region (CDR)3 length and their frequencies in sample R1. b 13-amino-acid CDR3 sequences of sample R1 were selected for WebLogo visualization. c Spectratype plot showing the distribution of TCR β CDR3 length and their frequencies in sample NR1. d 13-amino-acid CDR3 sequences of sample NR1 were selected for WebLogo visualization

We then determined the most frequent clonotype by percentage in both EPO responders and non-responders (Additional file 7a). We found that the most frequent clonotypes comprised 17.8%, 14.26%, 2.01%, 7.80%, 9.50%, 2.81% and 8.64% of reads in samples R1, R2, R3, R4, NR1, NR2, and NR3, respectively. In addition, we found 2, 2, 0, 1, 3, 0, and 2 clonotypes with frequencies of >5% in samples R1, R2, R3, R4, NR1, NR2, and NR3, respectively (Table 3).

Table 3 Most frequently emerging T-cell receptor (TCR) β sequences that constituted a proportion greater than 5% of all sequences

We further quantified the overall diversity of clonotypes using the inverse Simpson’s index (1/D), a mathematical index for assessing diversity for which the higher value means the greater the diversity. No significant difference in diversity indices between EPO responders and non-responders was observed (two-sample Fisher-Pitman permutation p = 0.642, Additional file 7b).

To further explore the difference between the TCR β repertoires in EPO responders and non-responders in terms of the clone type pattern, the sum of proportions of the top 5 clones was employed to compare the diversity in each sample (Fig. 3a). The 5 most abundant TCR β clone types constructed 22.1% ± 10.9% of sequences in EPO responders, compared to 18.8% ± 12.1% in EPO nonresponders (two-sample Fisher-Pitman permutation p = 0.682). We further calculated the pairwise clonotype repertoire diversity index (RDI) [20] among inter-group and intra-group samples, however, no significant difference was found across different groups (Kruskal-Wallis rank sum p = 0.792, Fig. 3b). Taken together, we concluded that no significant correlation between TCR β diversity and EPO responsiveness was observed.

Fig. 3
figure 3

Clonal diversity of T cells in end-stage renal disease patients. a T-cell receptor (TCR) repertoires are depicted to show the individual diversity in erythropoietin (EPO) responders compared to EPO non-responders. Colors represent the abundance of each clone in percentages. b Comparison of pairwise repertoire diversity indexes (RDIs) of clonotype distribution across sample groups

Characterization of gene usage of TRBV and TRBJ in ESRD patients

We next assess the gene usage of TRBV and TRBJ in EPO responders and non-responders. We observed that except several overexpression of specific TRBV genes, most ESRD patients shared a similar pattern of TRBV usage (Fig. 4a), regardless their responsiveness to EPO treatment. In addition, these patients also shared very similar TRBJ gene usage profiles (Fig. 4b). To validate this observation, we further conducted PCA analysis using TRBV and TRBJ gene usage frequencies. As result shown, EPO responders and non-responders were indistinguishable on the first and second principal component, suggesting the similar frequency of TRBV (Fig. 4c) and TRBJ (Fig. 4d) gene usage across two groups. Therefore, our results suggested that EPO responders and non-responders shared very similar TRBV and TRBJ gene usage profiles, despite some patient-specific gene enrichment pattern (especially on TRBV gene) was observed.

Fig. 4
figure 4

a Heatmap showing TRBV gene expression in all ESRD samples. b Heatmap showing TRBJ gene expression in all ESRD samples. c Figure depicting principal component analysis (PCA) result of TRBV gene usage across EPO responders and non-responders. d Figure depicting PCA result of TRBJ gene usage across EPO responders and non-responders

Specific joint TRBV/TRBJ2–7 gene usage was enriched in EPO responders

To clarify TRBV and TRBJ gene usage profiles in ESRD patients, we assessed the joint distribution between genes usages of Vβ and Jβ (Fig. 5a and Additional file 8). The top frequency of TRBV/TRBJ comprised 4.30%, 3.92%, 4.48%, 4.52%, 1.94%, 4.42% and 1.87% in samples R1, R2, R3, R4, NR1, NR2, and NR3, respectively. The high abundance (>3%) of joint Vβ/Jβ usage was observed in R1 (TRBV20–1/TRBJ2–1, TRBV20–1/TRBJ2–7, TRBV29–1/TRBJ2–7, TRBV29–1/TRBJ1–1, TRBV29–1/TRBJ2–1, TRBV29–1/TRBJ2–7), R2 (TRBV5–1/TRBJ1–1, TRBV5–1/TRBJ2–1, TRBV5–1/TRBJ2–7, TRBV6–1/TRBJ2–1, and TRBV6–1/TRBJ2–7), R3 (TRBV28/TRBJ2–1, TRBV28/TRBJ2–7), R4 (TRBV6–5/TRBJ2–1, TRBV6–5/TRBJ2–7), and NR2 (TRBV13/TRBJ2–7). However, we found no high abundance joint Vβ/Jβ usage in NR1 and NR3. In short, we observed high abundance rearrangements in all EPO responders (4/4, 100.0%) in contrast to 1 EPO non-responders (1/3, 33.3%; Chi-square p = 7.6 × 10−9).

Fig. 5
figure 5

a 3D bar plots quantifying joint distributions of TRBV and TRBJ in R1 patients. The x-axis represent TRBV genes, and the z-axis represent TRBJ genes. Y-axis represents frequency of each joint Vβ/Jβ rearrangement. b Correlations between the sum of top 1 ~ 30 most abundant joint distributions of Vβ/Jβ to erythropoietin (EPO) responsiveness in ESRD patients

To further assess the potential clinical utility of joint Vβ/Jβ profiles, we then asked whether joint Vβ/Jβ distributions may predict the EPO responsiveness of ESRD patients. We quantified the top abundant joint distributions of Vβ/Jβ, and this revealed a significant difference in the sum of top 3 to top 15 joint Vβ/Jβ distribution proportions between EPO-responsive and EPO-nonresponsive ESRD patients (two-sample Fisher-Pitman permutation p < 0.05, Fig. 5b). Together, we demonstrated different Vβ/Jβ usage patterns among EPO responders and non-responders, providing evidence for the correlation between TCR-mediated immune aberrations and EPO responsiveness (Fig. 6).

Fig. 6
figure 6

Schematic showing the TCR repertoire analysis focused on EPO responsiveness in ESRD patients. The x axis indicated the summation of top V-J combination(s) and the y axis indicated the significant value (−log10[P value]) regarding the sum of top V-J combination(s) across EPO responders and EPO nonresponders


In this study, we identified high abundance rearrangements in all EPO responders compared to EPO non-responders, and showed that EPO responders had a higher sum of the top 3 to 15 sequence proportions than EPO non-responders.

The immune profiles of ESRD patients can be very complicated, as both suppression and activation of the immune system may coexist in patients [21]. Chronic immune system activation, which implies expansion of T-lymphocytes and elevated levels of pro-inflammatory cytokines (including interleukin [IL]-1, tumor necrosis factor [TNF]-α, and interferon [IFN]-ɣ) and anti-inflammatory cytokines (including IL-10 and IL-13), which were secreted by regulatory T cells or other immune suppressive cell populations, may contribute to the poor responsiveness to EPO [22]. These cytokines can affect the process of erythropoiesis and are associated with low EPO sensitivity. In the aspect of T-cell phenotypes, Cooper et al. conducted flow cytometry experiments and revealed significant CD28 loss of both CD4+ and CD8+ T cells in poor EPO responders, which implied that the T cells may lose the costimulatory factors for full activation [23]. But whether the first activation signal induced by T cell receptor was turned on is still unknown. Therefore, it is important to investigate the TCR constituents (especially clonality) that may contribute to beneficial EPO treatment responses via adaptive immune effects. Nest generation sequencing technique provided us a very good platform to study T-cell receptor profiles in ESRD patients, which were poorly reported before.

We separately demonstrated the results of TCR repertoire analyses in diversity and specificity aspects to depict the hallmarks of TCR-mediated adaptive immunity in ESRD patients who received EPO treatment. As shown, we found no significant difference in TCR β repertoire diversity in either EPO responders or non-responders. Diversity index that measures biodiversity in ecological systems has been adopted to quantify the TCR diversity [24, 25]. Since the elevated IL-10 was observed in EPO non-responders compared to the normal controls, we hypothesize that those enriched clones in our non-responder group such as NR1 and NR3 may indicate the clonal expansion of Treg population, although the further studies are needed to verify this. Meanwhile, we noted the individual difference of intra-group TCR β repertoire diversity in both patient groups, for which inverse Simpson’s diversity indexes ranged 25 ~ 420 in EPO responders and 47 ~ 641 in EPO non-responders, indicating distinct immune environments among individual ESRD patients.

In the current study, we clearly showed that partial expansion of T cells was observed in all EPO responders which might imply the initiation of adaptive immune response induced by TCR. However, whether those expanded T-cell clones also secreted cytotoxic cytokines and factors need to be examined in future studies. Due to the lack of the human leukocyte antigen (HLA) types information of patients, identified different common CDR3 sequences in ESRD patients may indicate the different antigens loaded in different HLA molecules among patients. Our pilot study will strongly inspire the following researches to correlate HLA, immune-phenotypes and cytokines together to understand the complicated immune environment in ESRD patients.

As shown in Fig. 2 we identified highly enriched Vβ/Jβ rearrangements in EPO responder group, which may represent that EPO responders have their own preference on V and J usages to generate the specific CDR3 regions for immune responses. In addition, although the enriched Vβ/Jβ rearrangement was also detected in NR2 who is in non-responder group, this patient has very diverse CDR3 clonal types. In other words, the enriched Vβ/Jβ rearrangement was not correlated with enriched CDR3 sequences in this non-responder. This interesting result implies that the mechanism of the complicated TCR rearrangement for immune responses may be quite different between EPO responder group and non-responder group.

We mentioned several limitations of our study in the following: Firstly, the sample size is relatively small to make a strong conclusion, and may limit the power of the study. However, we still observed the enrichment of particular V-J combinations in EPO-resistant in compared to EPO-responsive ESRD patients. Therefore, we may expect a more concrete conclusion regarding the observed correlations in larger sample cohort. In considering that there are so many factors that may confound the observed correlation of TCR profiles across two groups, we assessed the possible confounders and found no significant difference between these covariates across EPO responders and non-responders. However, the correlations between TCR clonotypes and immune-related factors (pro-inflammatory cytokines and anti-inflammatory cytokines) in ESRD patients remains to be elucidated. In addition, the human leukocyte antigen (HLA) gene profiles have not been assessed in our study, which may also limit the generalizability of the results to the whole population.


In summary, we demonstrated a DNA sequencing-based method to characterize T-cell immunity in ESRD patients and measured the T-cell repertoire complexity to reveal similarities and differences between T-cell repertoire profiles underlying immune-pathological mechanisms of EPO responsiveness in ESRD patients (Fig. 6). As TCR repertoire behavior may affect the immune responses [26], our study provided a better understanding of the clonality of T-lymphocytes in ESRD patients that may explain the role of TCR-mediated adaptive immunity in EPO responsiveness.





Complementarity-determining region 3


Chronic kidney disease




Erythropoietin resistance index


Erythropoiesis-stimulating agents


End-stage renal disease


Glomerular filtration rate




Human leukocyte antigen


High-throughput sequencing






Peripheral blood mononuclear cells


Phosphate-buffered saline


Principal component analysis


Polymerase chain reaction


Red blood cell


Repertoire diversity index




Standard deviation


T-cell receptor


Tumor necrosis factor


United States Renal Data System


  1. Krishnan M, Franco E, McMurray S, Petra E, Nissenson AR. ESRD special needs plans: a proof of concept for integrated care. Nephrol News Issues. 2014;28(12):30. 32, 34-36

    PubMed  Google Scholar 

  2. The global issue of kidney disease. Lancet 2013, 382(9887):101.

  3. Chan TC, Fan IC, Liu MS, Su MD, Chiang PH. Addressing health disparities in chronic kidney disease. Int J Environ Res Public Health. 2014;11(12):12848–65.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Lankhorst CE, Wish JB. Anemia in renal disease: diagnosis and management. Blood Rev. 2010;24(1):39–47.

    Article  CAS  PubMed  Google Scholar 

  5. Ly J, Marticorena R, Donnelly S. Red blood cell survival in chronic renal failure. Am J Kidney Dis. 2004;44(4):715–9.

    Article  PubMed  Google Scholar 

  6. Szczech LA, Barnhart HX, Inrig JK, Reddan DN, Sapp S, Califf RM, Patel UD, Singh AK. Secondary analysis of the CHOIR trial epoetin-alpha dose and achieved hemoglobin outcomes. Kidney Int. 2008;74(6):791–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Messana JM, Chuang CC, Turenne M, Wheeler J, Turner J, Sleeman K, Tedeschi P, Hirth R. Association of quarterly average achieved hematocrit with mortality in dialysis patients: a time-dependent comorbidity-adjusted model. Am J Kidney Dis. 2009;53(3):503–12.

    Article  PubMed  Google Scholar 

  8. Icardi A, Paoletti E, De Nicola L, Mazzaferro S, Russo R, Cozzolino M. Renal anaemia and EPO hyporesponsiveness associated with vitamin D deficiency: the potential role of inflammation. Nephrol Dial Transplant. 2013;28(7):1672–9.

    Article  CAS  PubMed  Google Scholar 

  9. Honda H, Kimata N, Wakai K, Akizawa T. Associations among epoetin therapy, inflammation, nutritional status, and mortality in patients on hemodialysis. J Ren Nutr. 2014;24(5):322–9.

    Article  CAS  PubMed  Google Scholar 

  10. Meijers RW, Betjes MG, Baan CC, Litjens NH. T-cell ageing in end-stage renal disease patients: assessment and clinical relevance. World J Nephrol. 2014;3(4):268–76.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Betjes MG, Litjens NH. Chronic kidney disease and premature ageing of the adaptive immune response. Curr Urol Rep. 2015;16(1):471.

    Article  PubMed  Google Scholar 

  12. Litjens NH, Huisman M, van den Dorpel M, Betjes MG. Impaired immune responses and antigen-specific memory CD4+ T cells in hemodialysis patients. J Am Soc Nephrol. 2008;19(8):1483–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Zal B, Chitalia N, Ng YS, Trieu V, Javed S, Warrington R, Kaski JC, Banerjee D, Baboonian C. Killer cell immunoglobulin receptor profile on CD4(+) CD28(−) T cells and their pathogenic role in non-dialysis-dependent and dialysis-dependent chronic kidney disease patients. Immunology. 2015;145(1):105–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Klinger M, Kong K, Moorhead M, Weng L, Zheng J, Faham M. Combining next-generation sequencing and immune assays: a novel method for identification of antigen-specific T cells. PLoS One. 2013;8(9):e74231.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Chung S, Song HC, Shin SJ, Ihm SH, Park CS, Kim HY, Yang CW, Kim YS, Choi EJ, Kim YK. Relationship between erythropoietin resistance index and left ventricular mass and function and cardiovascular events in patients on chronic hemodialysis. Hemodial Int Symp Home Hemodial. 2012;16(2):181–7.

    Article  Google Scholar 

  16. Lopez-Gomez JM, Portoles JM, Aljama P. Factors that condition the response to erythropoietin in patients on hemodialysis and their relation to mortality. Kidney Int Suppl. 2008;111:S75–81.

    Article  CAS  Google Scholar 

  17. Locatelli F, Andrulli S, Memoli B, Maffei C, Del Vecchio L, Aterini S, De Simone W, Mandalari A, Brunori G, Amato M, et al. Nutritional-inflammation status and resistance to erythropoietin therapy in haemodialysis patients. Nephrol Dial Transplant. 2006;21(4):991–8.

    Article  CAS  PubMed  Google Scholar 

  18. Bolotin DA, Shugay M, Mamedov IZ, Putintseva EV, Turchaninova MA, Zvyagin IV, Britanova OV, Chudakov DM. MiTCR: software for T-cell receptor sequencing data analysis. Nat Methods. 2013;10(9):813–4.

    Article  CAS  PubMed  Google Scholar 

  19. Freeman JD, Warren RL, Webb JR, Nelson BH, Holt RA. Profiling the T-cell receptor beta-chain repertoire by massively parallel sequencing. Genome Res. 2009;19(10):1817–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Rubelt F, Bolen CR, McGuire HM, Vander Heiden JA, Gadala-Maria D, Levin M, Euskirchen GM, Mamedov MR, Swan GE, Dekker CL, et al. Individual heritable differences result in unique cell lymphocyte receptor repertoires of naive and antigen-experienced cells. Nat Commun. 2016;7:11112.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Macdougall IC. Could anti-inflammatory cytokine therapy improve poor treatment outcomes in dialysis patients? Nephrol Dial Transplant. 2004;19(Suppl 5):V73–8.

    Article  CAS  PubMed  Google Scholar 

  22. Macdougall IC, Cooper AC. Erythropoietin resistance: the role of inflammation and pro-inflammatory cytokines. Nephrolo Dial Transplant. 2002;17(Suppl 11):39–43.

    Article  CAS  Google Scholar 

  23. Cooper AC, Breen CP, Vyas B, Ochola J, Kemeny DM, Macdougall IC. Poor response to recombinant erythropoietin is associated with loss of T-lymphocyte CD28 expression and altered interleukin-10 production. Nephrol Dial Transplant. 2003;18(1):133–40.

    Article  CAS  PubMed  Google Scholar 

  24. Liu X, Venkataraman G, Lin J, Kiyotani K, Smith S, Montoya M, Nakamura Y, Kline J. Highly clonal regulatory T-cell population in follicular lymphoma - inverse correlation with the diversity of CD8 T cells. Oncoimmunology. 2015;4(5):e1002728.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Yew PY, Alachkar H, Yamaguchi R, Kiyotani K, Fang H, Yap KL, Liu HT, Wickrema A, Artz A, van Besien K, et al. Quantitative characterization of T-cell repertoire in allogeneic hematopoietic stem cell transplant recipients. Bone Marrow Transplant. 2015;50(9):1227–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Nikolich-Zugich J, Slifka MK, Messaoudi I. The many important facets of T-cell repertoire diversity. Nat Rev Immunol. 2004;4(2):123–32.

    Article  CAS  PubMed  Google Scholar 

Download references


I am grateful to Prof. Yusuke Nakamura (University of Chicago, USA) for reading the manuscript and providing valuable comments.


This work was supported by Grants from Ministry of Science and Technology, Taiwan (MY3NSC101–2320-B-038-029-MY3; MOST105-2628-B-038-001-MY4) and from Taipei Medical University, Taiwan (105TMU-CIT-02-3).

Availability of data and materials

Not applicable.

Author information

Authors and Affiliations



MSW, KCC, WCC and WCC conceived basic concept. HSCW, and YWH performed computational analysis. CMC and YWH conducted library preparation. HSCW, CMC, YWH, XL, MSW and WCC wrote the manuscript. MSW, KCC and WCC supervised and coordinated this research. All authors discussed and interpreted results of this manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Mai-Szu Wu or Wei-Chiao Chang.

Ethics declarations

Ethics approval and consent to participate

Samples were collected after obtaining informed consent from subjects at Taipei Medical University Hospital. The samples were collected from patients with written informed consent. The experiment was approved by the institutional review board at the Taipei Medical University and the methods were carried out in accordance with relevant guidelines and regulations.

Consent for publication

We accept the conditions of submission and the BioMed Central Copyright and License Agreement.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Schematic showing steps of T-cell receptor (TCR) β repertoire sequencing data processing. (DOCX 25 kb)

Additional file 2:

Overlapping read length distribution plot of end-stage renal disease patients. QC, quality control. (DOCX 260 kb)

Additional file 3:

Summary of repertoire sequencing output and overlapping proportions in end-stage renal disease samples. (DOCX 14 kb)

Additional file 4:

Mappability to variable/diversity/joining (VDJ) regions in end-stage renal disease patients. (DOCX 13 kb)

Additional file 5:

Summary of T-cell receptor (TCR) β clonotype statistics. (DOCX 14 kb)

Additional file 6:

T-cell receptor (TCR) β complementarity determining region (CDR)3 size spectratype plots of amino acid clonotypes from end-stage renal disease patients. The CDR3 length is defined as the number of amino acids between the conserved positions CDR3, i.e., cysteine (Cys)104 and phenylalanine (Phe)118. The color of the bar is annotated by the number of corresponding nucleotide read counts in the samples. (DOCX 155 kb)

Additional file 7:

Quantification of T-cell repertoire diversity in end-stage renal disease (ESRD) patients. a The clonotype distribution plots of ESRD patients. The x-axis and y-axis are log10 scaled. The value of 1/D represents the T-cell receptor (TCR) β repertoire diversity. b Comparison of TCR repertoire diversity (log-scaled) between erythropoietin (EPO) responders and EPO resistants (p = 0.642) (DOCX 232 kb)

Additional file 8:

3D bar plots quantifying joint distributions of TRBV and TRBJ in end-stage renal disease (ESRD) patients. (DOCX 517 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wong, H.SC., Chang, CM., Kao, CC. et al. V-J combinations of T-cell receptor predict responses to erythropoietin in end-stage renal disease patients. J Biomed Sci 24, 43 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: