Skip to main content

High throughput sequencing of T-cell receptor repertoire using dry blood spots

Abstract

Background

Immunology research, particularly next generation sequencing (NGS) of the immune T-cell receptor β (TCRβ) repertoire, has advanced progression in several fields, including treatment of various cancers and autoimmune diseases. This study aimed to identify the TCR repertoires from dry blood spots (DBS), a method that will help collecting real-world data for biomarker applications.

Methods

Finger-prick blood was collected onto a Whatman filter card. RNA was extracted from DBS of the filter card, and fully automated multiplex PCR was performed to generate a TCRβ chain library for next generation sequencing (NGS) analysis of unique CDR3s (uCDR3).

Results

We demonstrated that the dominant clonotypes from the DBS results recapitulated those found in whole blood. According to the statistical analysis and laboratory confirmation, 40 of 2-mm punch disks from the filter cards were enough to detect the shared top clones and have strong correlation in the uCDR3 discovery with whole blood. uCDR3 discovery was neither affected by storage temperatures (room temperature versus − 20 °C) nor storage durations (1, 14, and 28 days) when compared to whole blood. About 74–90% of top 50 uCDR3 clones of whole blood could also be detected from DBS. A low rate of clonotype sharing, 0.03–1.5%, was found among different individuals.

Conclusions

The DBS-based TCR repertoire profiling method is minimally invasive, provides convenient sampling, and incorporates fully automated library preparation. The system is sensitive to low RNA input, and the results are highly correlated with whole blood uCDR3 discovery allowing study scale-up to better understand the relationship and mutual influences between the immune and diseases.

Background

Immunology research, particularly next generation sequencing (NGS) of the immune T-cell receptor β (TCRβ) repertoire, has advanced progression in several fields including treatment of various cancers and autoimmune diseases [1]. The advantage of immunotherapy in cancers has led to increasing numbers of studies dedicated to exploring the impact and interaction between immunity and cancer cells. However, the complexity of the immune system in combination with the limitation of detection methods makes this subject difficult to research.

Currently, TCR repertoire sequencing is widely used to evaluate the immune system [2]. TCR repertoires of patients were explored in a variety of disorders—patients under cancer immunotherapy, autoimmune disease [rheumatoid arthritis (RA), ankylosing spondylitis], or subject to virus infection (hepatitis, human immunodeficiency virus (HIV)) [3,4,5,6,7]. Analyzing the TCR repertoire may help to gain a better understanding of the immune system features and of the etiology and progression of diseases, in particular those with unknown antigenic triggers. Rapid progress has been made in the deep profiling of TCR repertoires by using NGS to discover millions of sequences from the TCR repertoires [8]. There are different TCR library preparation for NGS, including: multiplex PCR, targeted enrichment methods, 5′Rapid amplification of cDNA ends (5′RACE) cDNA synthesis, template-switch, and nested PCR [9]. For TCR repertoire assessment, the current standard sample type is total RNA extracted from whole blood by venipuncture. However, the costs, participant burden, regulatory constraints, and logistics associated with venipuncture and RNA handling are major barriers to clinical application or community-based research on various diseases.

Dry blood spots (DBS) have been used broadly in disease screening, drug level monitoring, and infectious microorganism detection, such as the detection of HIV and plasmodium [10, 11]. There are two major advantages of this method. First, it is minimally invasive (requiring only 0.3 mL of capillary blood obtained by finger- or ear-prick) and a field-friendly alternative to venipuncture by professional medical staff. Second, it is low cost and convenient to handle in that the filter cards can be stored and shipped at room temperature without refrigeration for at least 1 month [12, 13].

The desire for TCR repertoire profiling in a variety of diseases and the advantages of a DBS-based method prompted us to establish a DBS-based TCR repertoire profiling method. Here, we report for the first time the use of a DBS sample in combination with a fully automated and closed system for TCRβ chain NGS library preparation. We have demonstrated that storage of DBS in filter cards up to 28 days at either room temperature or − 20 °C has no effect on unique CDR3 (uCDR3) discovery. Compared with whole blood, the DBS-based TCR repertoire profiling method is minimally invasive, provides convenience when sampling, is sensitive, and is highly correlated with whole blood in terms uCDR3 discovery, facilitating ease of incorporation and scale-up to studies which seek to explore the relationship between immunity in variety diseases.

Materials and methods

Subjects and sample collection

All enrolled subjects had written informed consent prior to the collection of whole blood or DBS via finger-prick. The study design and recruiting of patients were approved by the New England Independent Review Board® (NEIRB) (IRB Number: 14-378). All experiments were performed according to the relevant guidelines and regulations.

Whatman FTA 903 cards (Sigma-Aldrich Corp. MO, USA) were used for collection of capillary blood obtained by finger-prick. In addition, peripheral blood was collected via venipuncture for comparison. The Whatman cards were stored in zip-lock plastic bags with a desiccant (ULINE. Atlanta, USA), at room temperature or − 20 °C (for longer storage) for future studies.

RNA extraction from DBS and whole blood

RNA extraction from DBS was performed following the methods proposed by Karlsson et al. [14]. 2-mm disks were punched out from a DBS with a sterile Robbins true-cut disposable 2-mm biopsy punch (Robbins Instruments, Inc., NY, USA). The punch disks were then incubated in 700 μL RLT buffer plus β-mercaptoethanol (β-ME). The tubes were then incubated at 37 °C for 60 min in a thermomixer rotating at 1000 rpm. Following incubation, RNA was extracted with a RNeasy Micro Kit (Qiagen, Valencia, CA, USA) according to the manufacturer’s instruction. For whole blood, a RNeasy Mini Kit (Qiagen, Valencia, CA, USA) was used according to the manufacturer’s protocol. RNA concentrations were measured by spectrophotometry.

RT-PCR and sequencing of T-cell receptor β-chain CDR3 region

iRepertoire multiplex primer sets (iRepertoire, Inc. Huntsville, AL, USA) were used to amplify the CDR3 region of TCRβ chain by using RNA as template as described by Wang et al. [15]. The whole amplification process and library preparation process for NGS were fully automated in the iR-Procecessor and iR-Cassette (iRepertoire, Inc. Huntsville, AL, USA). Then, paired-end sequencing was performed on purified PCR products using an Illumina MiSeq v2 300-cycle Reagent Kit (Illumina Inc.), for an average read depth of 30,000 reads per sample.

Sequence data analysis

Raw cDNA sequences were first analyzed to identify V and J genes by using iR-map and visualized in iRweb (iRepertoire, Inc. AL, USA). Analyzed data from iRweb include peptide sequences, alignments to the international ImMunoGeneTics (IMGT) database, uCDR3, shared CDR3s, and V- and J-gene usage. Multiple alignments and hierarchical clustering of conserved amino acid sequences were analyzed as described by Wang et al. [15].

Computational error correction of bulk TCR sequences by replicates

Errors in sequencing resulting from PCR errors, PCR contamination and read error were mitigated according to the modified method of Glanville et al.’s study [16]. RNA samples were split into two reactions and processed as technical replicates. The coefficient of determination (R2) was calculated by linear regression to show the correlation between the replicates of TCRβ CDR3 frequencies prior to data analysis to exclude PCR and sequence errors.

TCR repertoire diversity analysis

The diversity of the TCR repertoire was calculated based on the diversity 50 (D50) value and the Diversity Index (DI) [17]. Briefly, D50 is the percent of dominant T cell clones that account for the cumulative 50% of the total CDR3s counted in the sample. The mathematical formulation of D50 is defined as follows:

$${\text{D}}50 = \left( {{\text{No}}.\;{\text{of uCDR}}3{\text{ that make up }}50\% {\text{ of the total reads}} \times 100} \right)/{\text{No}}.{\text{ of uCDR}}3{\text{s}}$$

The more diverse a library, the closer the value will be to 50.

The DI is defined mathematically as follows:

$${\text{Assume}}\;{\text{that the numbers }}\left( {\text{n}} \right){\text{ of uCDR3}}{:}\;r_{1} \, \ge \,r_{2} \, \ge \,r_{3} \, \ge \, \cdots r_{i} \, \ge \,r_{i + 1} \ge \cdots \, \ge \,r_{n}$$

where ri is the frequency of the i-th CDR3 and n is the total number of unique CDR3s.

$$x_{k} = \frac{k}{n},\;\;\;\;\;y_{k } = \frac{{\mathop {\sum\nolimits_{i = 1}^{k} {r_{i} } }\nolimits_{{}}^{{}} }}{{\sum\nolimits_{i = 1}^{n} {r_{i} } }}$$

In addition, tree maps and two-dimensional (2D) heat maps were used to reveal the diversity and characteristics of TCR repertoire. In a tree map, each rounded rectangle represents a unique entry: V–J combination of uCDR3, where the size of a spot denotes the relative frequency. 2D heat maps showed that the relative frequency of a consensus germline V-gene allele (as per alignment with the IMGT database) is plotted relative to the consensus germline J-gene allele. Therefore, it is immediately evident which V–J combination is used either frequently or infrequently by the color of the map.

Statistical analysis

Analysis was performed using the Statistical Package for the Social Sciences (SPSS) 17.0 software. Chi-squared test was performed on all categorical variables, with the exception of those with an expected frequency of < 5, which were analyzed by Fisher’s exact test. Two-sided Student’s t-tests was used for comparing the means of two independent variables. One-way analysis of variance (ANOVA) was used to analyze differences in the number of uCDR3s among the different storage durations.

Results

T-cell repertoire can be detected from the extracted RNA from DBS

The work flow of the DBS-based TCR repertoire discovery is shown in Fig. 1. RNA was extracted from DBS, followed by RT-PCR in a fully automated closed system to amplify the TCRβ library. After cDNA pooling of the TCRβ library and quantification, NGS sequencing was performed, and the data was analyzed as described previously.

Fig. 1
figure 1

The work flow of DBS based TCR repertoire discovery. Finger-prick blood is collected onto the filter card and extracted RNA from DBS is used as a template for amplification. The entire PCR amplification and library preparation are automated in a closed cassette. Post amplification, libraries are pooled and quantified then sequenced by NGS. Sequencing data is analyzed by iRweb’s bioinformatics platform

Since the percentage of dominant clones detected from DBS is highly dependent on TCR diversity of samples, we first calculated the number of punch disks required for dominant clone detection from whole TCR repertoire. Both confidence levels and intervals were calculated to compare TCR repertoire discovery between methods utilizing DBS versus 10 mL of whole blood. According to the normal range of differential count of complete blood cells, 20–35% of white blood cells are lymphocytes. T cells account for 70–85% of lymphocytes [18]. Therefore, there were about 700–2675 T cells per μL of whole blood for one man with 65 kg of body weight with normal white blood count range from 5000–9000 WBC/μL. According to the Whatman card specifications, one circle with a diameter of 0.5 in. can absorb 75–80 μL of blood; therefore, one 2-mm punch disk contains approximately 2 μL of blood. For individuals whose dominant clones comprise 60% of their total TCR repertoire, 80 μL of blood collected on 40 2-mm punch disks could provide 90% confidence level and 9% confidence interval. The volume of blood required for a given confidence level and interval is presented in Additional file 1: Table S1.

Having calculated the punch disk numbers required for TCR repertoire discovery, we compared between 3, 10, 40 and 80 punch disks. The average uCDR3 discovery for 3, 10, 40, and 80 disks were 169, 246, 1726, and 1581, respectively, with an average read depth of 30,000 reads per sample. The uCDR3 discovery from 40 disks was higher than those from 3 (p < 0.0001) or 10 punch disks (p < 0.0001). There was no significant statistical difference in average uCDR3 numbers between 40 and 80 punch disks (p = 0.520) (Fig. 2a). The tree maps and heat maps also revealed that the diversity and the identified combination of V- and J- gene segments were obviously higher in 40 or 80 punch disks than those in 3 or 10 punch disks (Fig. 2b, c). As a result, 40 disks were established as the optimum punch number for future studies.

Fig. 2
figure 2

T-cell repertoire is detected from the extracted DBS RNA. a uCDR3 discovery from different amounts of punch disks (**p < 0.01, by Student’s t-test). b TCR repertoire diversity is illustrated in tree maps where each rounded rectangle represents a unique entry: V–J–uCDR3 and the size of the spot denotes the relative frequency. c The relative frequency from low (black) to high (red) of V–J gene combinations, is illustrated in the 2D heat maps from varying amounts of punch disks (left to right: 3, 10, 40 and 80 disks). ns indicates no significant statistical difference

CDR3 frequency correlation between DBS and whole blood

The R2 of uCDR3 frequencies were calculated by linear regression between technical replicates to exclude PCR or sequencing error. The R2 of the regression for CDR3 discovery between technical replicates for DBS and whole blood samples were 0.990 and 0.995, respectively (Fig. 3a–d). uCDR3 frequency correlations were then determined to compare discovery between DBS and whole blood samples. uCDR3 discovery between DBS and whole blood revealed high regression correlation (R2 = 0.986 and R2 = 0.987, respectively) (Fig. 3e), with 94% and 92% of the top 50 and top 100 dominant clones from whole blood being detected from DBS, respectively (data not shown).

Fig. 3
figure 3

The comparison of TCRβ CDR3 between DBS and whole blood. a Tree maps of TCRβ CDR3 from duplicate DBS samples. b The coefficient of determination (R2) of linear regression in TCRβ CDR3 sequences from duplicate DBS samples. c Tree maps of TCRβ CDR3 from duplicate peripheral blood samples. d The R2 of linear regression in TCRβ CDR3 sequences from duplicate peripheral blood samples. e The R2 of linear regression in TCRβ CDR3 sequences between peripheral blood and DBS

Effect of storage temperature and duration on uCDR3 discovery

Despite the manufacturer’s suggested storage temperature of − 20 °C, prior studies show conflicting results for the effect of DBS storage temperature and duration on RNA preservation [19,20,21,22]. In this study, we compared the uCDR3 discovery between two storage temperatures (room temperature and − 20 °C). No significant difference was demonstrated in uCDR3 discovery between DBS samples stored at room temperature and − 20 °C for 14 days (p = 0.922) or 28 days (p = 0.700) (Fig. 4).

Fig. 4
figure 4

The comparison of uCDR3 discovery between different storage temperatures (room temperature versus − 20 °C). Student’s t test was used to determine that there was no significant difference in uCDR3 discovery between the two storage temperatures at both 14 days (p = 0.922) and 28 days (p = 0.700) storage duration. ns indicates no significant statistic difference

uCDR3 discovery from storage duration tests (1 day, 14 days and 28 days) was also compared in DBS samples with low, moderate, and high CDR3 diversities. Among different storage durations, no significant differences were found in uCDR3 discovery in DBS samples with low (p = 0.190), moderate (p = 0.077), and high (p = 0.857) uCDR3 diversities (Fig. 5), and this is true even under the two different storage temperatures tested, room temperature, and − 20 °C (Fig. 4). The uCDR3 frequency did not decrease over time. These results indicate that storage duration of 1, 14 and 28 days had no impact on uCDR3 discovery.

Fig. 5
figure 5

Comparison of uCDR3 discovery with storage durations of 1, 14 and 28 days from DBS samples with varying TCR diversities. There was no significant difference in TCRβ uCDR3 discovery between storage durations for samples with varying TCR diversities by one-way ANOVA analysis. ns indicates no significant statistic difference

We also demonstrated the relationship of TCRβ clonotype discovery between DBS and whole blood (Fig. 6). The R2 of the regression for CDR3 discovery between technical replicates for whole blood samples and DBS in different storage duration tests were 0.990 and 0.86–0.95, respectively. Long storage duration up to 28 days not only showed no impact on uCDR3 discovery, but also showed no impact on the correlation uCDR3 discovery between DBS and whole blood. In addition, the top 50 and 100 dominant CDR3 clones of whole blood could be detected by DBS for 62–88% and 45–69% of the dominant unique CDR3s at 28-day storage, respectively (Fig. 6).

Fig. 6
figure 6

Comparing uCDR3 discovery between different storage durations (1, 7, 14 and 28 days). The R2 of TCRβ CDR3 sequences was calculated by linear regression between the whole blood and DBS in different storage durations. The R2 of duplicate sequences showed high correlation in whole blood and DBS within different storage durations. The TCR diversities were shown by D50s, Diversity Indexes and tree maps. There was no significant difference in TCR diversity between storage durations. In addition, more than 74% of the top 50 dominant CDR3 clones of whole blood was detected by DBS. These results indicated that TCR uCDR3s of DBS had high correlation with those of whole blood. (R2: coefficient of determination calculated by linear regression; *The percentage of TCR dominant clones from whole blood could be detected from the DBS sample)

Relationship between TCR repertoire diversity and uCDR3 discovery

The diversity of each sample’s TCR repertoire can impact the correlation of uCDR3 discovery between DBS and whole blood samples. To explore this relationship, uCDR3 discovery was compared between DBS and whole blood samples known to have low, moderate, and high diversities of TCR repertoire. Table 1 shows the characteristics of TCRβ CDR3 between DBS and whole blood. The samples with low and moderate TCR diversity have a higher regression correlation between DBS and whole blood samples than those with higher TCR diversity (R2 = 0.989 for low diversity, 0.865 for moderate diversity and 0.806 for high diversify). 92% and 87% of the top 100 uCDR3 clones from whole blood could be detected in a DBS sample of low and moderate diversity, respectively, but only 74% of those in high diversity samples. For samples with high TCR diversity, the top 50 uCDR3 clones from whole blood could still be detected for 82% in a DBS sample. Although the diversity of samples had impact on the correlation of uCDR3 discovery, there was still a relatively high amount of shared uCDR3s between DBS sample and whole blood sample.

Table 1 The clinical characteristics of samples with low, moderate and high TCR repertoire diversities

uCDR3 sharing between samples

In order to determine whether the TCRβ CDR3s detected by this method are unique across individuals, the TCRβ CDR3 sharing was compared among different individuals. The sharing rate of uCDR3 between two individuals was low, ranging from 0.03% to 1.55% (Table 2 and Additional file 1: Table S2). The uCDR3 sharing rate between P5 and P6 is the highest; however, this is a parent–child relationship.

Table 2 The ratio of the shared uCDR3 reads divided by the total reads of the individual

Discussion

In this study, we reported for the first time the establishment of a DBS-based method as a minimally invasive alternative to venipuncture blood sampling for TCR repertoire studies. In addition to its cost effectiveness, we have demonstrated that the storage temperature and storage duration had no effect on uCDR3 discovery. The ease of sample collection and the convenience of storage for transportation purposes makes it a viable option for scale-up studies related to TCR repertoires.

The FTA cards entrap nucleic acids and are stabilized at room temperature, which allows long-term storage [20, 23, 24]. DNA extracted from the stored card yields adequate quality and amount for PCR amplification compared with that extracted from the fresh samples [25, 26]. The extracted DNA from the card can also be used for high-throughput molecular analysis [13, 27]. In addition, DBS also provides an alternative tool to formalin-fixed paraffin-embedded cell blocks for biobanking to detect epidermal growth factor receptor (EGFR) or Kirsten rat sarcoma 2 viral oncogene homolog (KRAS) mutation for personalized target therapy in lung cancer patients [28]. Therefore, this method can be utilized in multiple fields, not only cancer diagnosis but various disease etiologies. In addition to the stability of high-quality nucleic acids, the low cost of storage, and ease in transportation, FTA cards serve as an alternative to various traditional storage methods [20].

We established a method that uses RNA extracted from DBS for high-throughput sequencing. The entire process of RT-PCR amplification and library preparation of the immune repertoire was conducted automatically in cassettes. Automated amplification in a closed, contained environment throughout the entire process prevents contamination and improves the consistency of the results by avoiding error or variation, which may be introduced during manual processing. TCR repertoires derived from DBS are directly comparable with the dominant clones of TCR repertoire discovery from whole blood. As a result, DBS may be used as a screening tool for certain diseases or for genetic screening instead of whole blood.

Peluso et al. used FTA cards to store fine needle aspiration cytology samples of lymph node, and they amplified immunoglobulin heavy and light chains, TCR-β and γ chains by PCR from DNA [29]. The results showed that the amplified DNA from FTA cards is comparable to those of cryopreserved samples. However, a tremendous diversity of TCR repertoires results from extensive recombination, splicing and post-transcriptional processing to yield functional proteins. TCR repertoire discovery from RNA rather than genomic DNA prevents the inclusion of nonproductive or functionally irrelevant sequences [30]. There is currently no study report on the direct detection of the TCRβ repertoire from finger-pricked blood stored on FTA cards.

Although it is known that high-quality DNA can be extracted from the FTA cards [25, 31], there are few studies about RNA extracted from FTA cards compared with other preservation methods. Instability of RNA is a major concern with storage, temperature, and extraction methods. The World Health Organization (WHO) suggests that for HIV drug resistance, dried blood spot specimens should be transferred to − 20 °C or lower as soon as possible although they can be kept and/or transported at ambient temperature up to 14 days after collection [10]. However, Bertagolio et al. demonstrated that a 90% positive amplification rate was noted under the ambient storage temperature of DBS [19]. This is compatible with the current DBS study which shows comparable numbers of uCDR3s between disparate storage temperatures of room temperature and − 20 °C, even after 28 days of storage duration. Therefore, DBS can be stored and transported at room temperature without impacting TCRβ clonotype discovery.

Approximately 74–98% of the top 50 dominant clones from the whole blood are found in the DBS sample. This result indicates that the DBS-based immune repertoire profiling is very sensitive. According to the statistical analysis of sample sizes and the results of the current study, we are now confident that our DBS-based method can reflect whole blood TCRβ repertoire sampling with very good coverage and high sensitivity, i.e., enough shared top clonotypes and strong correlation in uCDR3 discovery.

The recent technological advances in the fields of immunology and immunotherapy hold promise for patients with cancers [32,33,34]. The TCR repertoire plays a pivotal factor in immunity [35]. For instance, TCR repertoire profiling has the potential to serve as a biomarker of treatment response in pancreatic ductal carcinoma patients who received immunotherapy [7]. Through recognition of major histocompatibility complex (MHC)-peptide complex of TCRs, T cells can be activated and specific T cell clones expanded to give a response to foreign pathogens or cancer cells [36].

Accumulating evidence suggests that immune cells play pivotal roles in a variety diseases and their overlapping regulatory mechanisms in addition to cancers, and TCR repertoire were also explored in different area. For example: Muraro et al. used TCR repertoire to follow the effect of autologous stem cell transplantation in multiple sclerosis patients [37]. For autoimmune diseases, the median TCR β-chain frequency in RA patients was increased tenfold, indicating marked contraction of the repertoire, and disease activity of RA was negatively correlated with the TCR repertoire diversity of CD4+ T cells [4, 38]. Even for infection disease, TCR repertoire was explored the function profile of T cells in HIV patients and the characteristics of immune landscape after injection of HIV vaccine [6, 39]. TCRβ-CDR3 were correlated with disease prognosis of hepatitis B and C [5, 40, 41]. DBS-based methodology provides a simple and quick approach to access changes of immunity and treatment effect in different disease areas.

Each individual has more than 109 T cells which express unique heterodimeric T-cell receptors as identified by the CDR3 peptide sequence of the TCRβ chain. This high diversity provides defense against foreign pathogens. Hou et al. reported that on average two individuals can share 4.85 ± 2.5% of their DNA sequences and 12.17 ± 0.81% of their CDR3 amino acid sequences [42]. According to the current study, less than 1% of uCDR3 clones are shared between two different people except when a parent–child relationship exists (Table 2). Putintseva et al. reported that the degree of overlap was always slightly higher for related individuals for all CDR3s, but this difference never approached a significant level compared to unrelated individuals [43]. In a parent–child relationship, 50% of their HLA alleles are shared. Exposure to the same pathogens may result in specific T cell expansion. Because V-gene usage is highly influenced by human leukocyte antigen (HLA)-type, certain V-genes may be preferentially selected by the immune system, ultimately allowing more shared CDR3s among closely related individuals. The high individual specificity detected by this method allows researchers to look for public CDR3s in a homogenous group or identify disease signatures, based on shared unique CDR3 peptides.

Figure 7 demonstrates the advantages of the DBS-based method for TCR repertoire analysis. DBS provides a minimally invasive sample collection alternative to venipuncture. It is convenient for storage and transportation, in contrast to the necessary of refrigeration of − 20 °C for some RNA-compatible vacutainer tubes, and cold chain shipping for whole blood. Then, high quality RNA can be extracted from DBS-based samples for future molecular amplification and analysis. The entire process of RT-PCR amplification and library preparation of the immune repertoire was conducted automatically in a fully closed cassette, which is easier than the labor-intensive manual process. It also helps improve result consistency and to prevent contamination of the laboratory environment with amplicon. Although other products (PAXgene Blood RNA Tube (IVD) (Qiagen), Tempus™ Blood RNA Tube, etc.) can be used for storage and transportation of RNA at room temperature, refrigeration is still necessary for long-term storage. In addition, they are more expensive and laborious than DBS-based processing due to the requirement for a phlebotomist and specialized extraction RNA kits.

Fig. 7
figure 7

Schematics for the comparison of the methods for TCR repertoire discovery by using a DBS sample and peripheral blood sample. DBS provides a minimally invasive sample collection and convenient alternative to venipuncture. High quality RNA is extracted from DBS-based samples for future molecular amplification and analysis. The entire process of multiplex PCR amplification and library preparation of the immune repertoire is automated in a fully closed cassette, which allows for a less labor-intensive manual process, improves result consistency and prevents contamination amplicon within the laboratory environment

Conclusions

The current DBS-based method shows high correlation to dominant TCRβ CDR3 clones detected from whole blood. The minimally invasive, cost effective, and convenient storage conditions will allow for the scale-up of surveillance studies or treatment response. More importantly, the development of such a tool allowed collecting of real-world data, i.e., when symptoms showed up, rather than doctors available for an appointment.

Abbreviations

NGS:

next generation sequencing

TCR:

T-cell receptor

5′RACE:

5′Rapid amplification of cDNA ends

DBS:

dry blood spots

RA:

rheumatoid arthritis

HIV:

human immunodeficiency virus

uCDR3:

unique CDR3s

β-ME:

β-mercaptoethanol

IMGT:

ImMunoGeneTics

DI:

Diversity Index

R 2 :

the coefficient of determination

2D:

two-dimensional

SPSS:

Statistical Package for the Social Sciences

ANOVA:

one-way analysis of variance

EGFR :

epidermal growth factor receptor

KRAS :

Kirsten rat sarcoma 2 viral oncogene homolog

WHO:

The World Health Organization

MHC:

major histocompatibility complex

References

  1. Davis MM, Tato CM, Furman D. Systems immunology: just getting started. Nat Immunol. 2017;18:725–32.

    Article  CAS  Google Scholar 

  2. Ichinohe T, Miyama T, Kawase T, Honjo Y, Kitaura K, Sato H, et al. Next-generation immune repertoire sequencing as a clue to elucidate the landscape of immune modulation by host–gut microbiome interactions. Front Immunol. 2018;9:668.

    Article  Google Scholar 

  3. Cui JH, Jin YB, Lin KR, Xiao P, Chen XP, Pan YM, et al. Characterization of peripheral blood TCR repertoire in patients with ankylosing spondylitis by high-throughput sequencing. Hum Immunol. 2018;79:485–90.

    Article  CAS  Google Scholar 

  4. Sakurai K, Ishigaki K, Shoda H, Nagafuchi Y, Tsuchida Y, Sumitomo S, et al. HLA-DRB1 shared epitope alleles and disease activity are correlated with reduced T cell receptor repertoire diversity in CD4+ T cells in rheumatoid arthritis. J Rheumatol. 2018;45:905–14.

    Article  CAS  Google Scholar 

  5. Ping Y, Song M, Wang M, Li Z, Zhang Y. CDR3 repertoire diversity of CD8+ T lymphocytes in patients with HCV. Cell Immunol. 2018;56:56. https://doi.org/10.1016/j.cellimm.2018.12.007 (Epub ahead of print).

    Article  CAS  Google Scholar 

  6. Wendel BS, Del Alcazar D, He C, Del Rio-Estrada PM, Aiamkitsumrit B, Ablanedo-Terrazas Y, et al. The receptor repertoire and functional profile of follicular T cells in HIV-infected lymph nodes. Sci Immunol. 2018;3:eaan8884. https://doi.org/10.1126/sciimmunol.aan8884.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Hopkins AC, Yarchoan M, Durham JN, Yusko EC, Rytlewski JA, Robins HS, et al. T cell receptor repertoire features associated with survival in immunotherapy-treated pancreatic ductal adenocarcinoma. JCI Insight. 2018;3. https://doi.org/10.1172/jci.insight.122092.

  8. Benichou J, Ben-Hamo R, Louzoun Y, Efroni S. Rep-Seq: uncovering the immunological repertoire through next-generation sequencing. Immunology. 2012;135:183–91.

    Article  CAS  Google Scholar 

  9. Rosati E, Dowds CM, Liaskou E, Henriksen EKK, Karlsen TH, Franke A. Overview of methodologies for T-cell receptor repertoire analysis. BMC Biotechnol. 2017;17:61.

    Article  Google Scholar 

  10. WHO. WHO manual for HIV drug resistance testing using dried blood spot specimens. Geneva: WHO; 2010.

    Google Scholar 

  11. Zainabadi K, Adams M, Han ZY, Lwin HW, Han KT, Ouattara A, et al. A novel method for extracting nucleic acids from dried blood spots for ultrasensitive detection of low-density Plasmodium falciparum and Plasmodium vivax infections. Malar J. 2017;16:377.

    Article  Google Scholar 

  12. Rahikainen AL, Palo JU, de Leeuw W, Budowle B, Sajantila A. DNA quality and quantity from up to 16 years old post-mortem blood stored on FTA cards. Forensic Sci Int. 2016;261:148–53.

    Article  CAS  Google Scholar 

  13. Saieg MA, Geddie WR, Boerner SL, Liu N, Tsao M, Zhang T, et al. The use of FTA cards for preserving unfixed cytological material for high-throughput molecular analysis. Cancer Cytopathol. 2012;120:206–14.

    Article  CAS  Google Scholar 

  14. Karlsson H, Guthenberg C, von Dobeln U, Kristenssson K. Extraction of RNA from dried blood on filter papers after long-term storage. Clin Chem. 2003;49:979–81.

    Article  CAS  Google Scholar 

  15. Wang C, Sanders CM, Yang Q, Schroeder HW Jr, Wang E, Babrzadeh F, et al. High throughput sequencing reveals a complex pattern of dynamic interrelationships among human T cell subsets. Proc Natl Acad Sci USA. 2010;107:1518–23.

    Article  CAS  Google Scholar 

  16. Glanville J, Huang H, Nau A, Hatton O, Wagar LE, Rubelt F, et al. Identifying specificity groups in the T cell receptor repertoire. Nature. 2017;547:94–8.

    Article  CAS  Google Scholar 

  17. Han J. Immunodiversity assessment method and its use. US20120183969 USA; 2012.

  18. Chernecky CC, Berger BJ, Laboratory tests and diagnostic procedures. 6th ed. In: Chernecky CC, Berger BJ, editors. St. Louis: Elsevier; 2013. p. 432–59.

  19. Bertagnolio S, Soto-Ramirez L, Pilon R, Rodriguez R, Viveros M, Fuentes L, et al. HIV-1 drug resistance surveillance using dried whole blood spots. Antivir Ther. 2007;12:107–13.

    CAS  PubMed  Google Scholar 

  20. da Cunha Santos G. FTA cards for preservation of nucleic acids for molecular assays: a review on the use of cytologic/tissue samples. Arch Pathol Lab Med. 2018;142:308–12.

    Article  Google Scholar 

  21. Garcia-Lerma JG, McNulty A, Jennings C, Huang D, Heneine W, Bremer JW. Rapid decline in the efficiency of HIV drug resistance genotyping from dried blood spots (DBS) and dried plasma spots (DPS) stored at 37 °C and high humidity. J Antimicrob Chemother. 2009;64:33–6.

    Article  CAS  Google Scholar 

  22. Nelson JA, Loftis AM, Kamwendo D, Fawzi WW, Taha TE, Goldenberg RL, et al. Nevirapine resistance in human immunodeficiency virus type 1-positive infants determined using dried blood spots stored for up to six years at room temperature. J Clin Microbiol. 2009;47:1209–11.

    Article  CAS  Google Scholar 

  23. Barth H, Morel A, Mougin C, Averous G, Legrain M, Fender M, et al. Long-term storage and safe retrieval of human papillomavirus DNA using FTA elute cards. J Virol Methods. 2016;229:60–5.

    Article  CAS  Google Scholar 

  24. Rajendram D, Ayenza R, Holder FM, Moran B, Long T, Shah HN. Long-term storage and safe retrieval of DNA from microorganisms for molecular analysis using FTA matrix cards. J Microbiol Methods. 2006;67:582–92.

    Article  CAS  Google Scholar 

  25. Peluso AL, Cascone AM, Lucchese L, Cozzolino I, Ieni A, Mignogna C, et al. Use of FTA cards for the storage of breast carcinoma nucleic acid on fine-needle aspiration samples. Cancer Cytopathol. 2015;123:582–92.

    Article  CAS  Google Scholar 

  26. Dobbs LJ, Madigan MN, Carter AB, Earls L. Use of FTA gene guard filter paper for the storage and transportation of tumor cells for molecular testing. Arch Pathol Lab Med. 2002;126:56–63.

    CAS  PubMed  Google Scholar 

  27. Kampmann ML, Buchard A, Borsting C, Morling N. High-throughput sequencing of forensic genetic samples using punches of FTA cards with buccal swabs. Biotechniques. 2016;61:149–51.

    Article  CAS  Google Scholar 

  28. da Cunha Santos G, Liu N, Tsao MS, Kamel-Reid S, Chin K, Geddie WR. Detection of EGFR and KRAS mutations in fine-needle aspirates stored on Whatman FTA cards: is this the tool for biobanking cytological samples in the molecular era? Cancer Cytopathol. 2010;118:450–6.

    Article  Google Scholar 

  29. Peluso AL, Cozzolino I, Bottiglieri A, Lucchese L, Di Crescenzo RM, Langella M, et al. Immunoglobulin heavy and light chains and T-cell receptor beta and gamma chains PCR assessment on cytological samples. A study comparing FTA cards and cryopreserved lymph node fine-needle cytology. Cytopathology. 2017;28:203–15.

    Article  CAS  Google Scholar 

  30. Calis JJ, Rosenberg BR. Characterizing immune repertoires by high throughput sequencing: strategies and applications. Trends Immunol. 2014;35:581–90.

    Article  CAS  Google Scholar 

  31. McClure MC, McKay SD, Schnabel RD, Taylor JF. Assessment of DNA extracted from FTA cards for use on the Illumina iSelect BeadChip. BMC Res Notes. 2009;2:107.

    Article  Google Scholar 

  32. Kaufman HL, Kirkwood JM, Hodi FS, Agarwala S, Amatruda T, Bines SD, et al. The Society for Immunotherapy of Cancer consensus statement on tumour immunotherapy for the treatment of cutaneous melanoma. Nat Rev Clin Oncol. 2013;10:588–98.

    Article  CAS  Google Scholar 

  33. Robert L, Harview C, Emerson R, Wang X, Mok S, Homet B, et al. Distinct immunological mechanisms of CTLA-4 and PD-1 blockade revealed by analyzing TCR usage in blood lymphocytes. Oncoimmunology. 2014;3:e29244.

    Article  Google Scholar 

  34. Robert L, Tsoi J, Wang X, Emerson R, Homet B, Chodon T, et al. CTLA4 blockade broadens the peripheral T-cell receptor repertoire. Clin Cancer Res. 2014;20:2424–32.

    Article  CAS  Google Scholar 

  35. Cooper MD, Alder MN. The evolution of adaptive immune systems. Cell. 2006;124:815–22.

    Article  CAS  Google Scholar 

  36. Castro CD, Luoma AM, Adams EJ. Coevolution of T-cell receptors with MHC and non-MHC ligands. Immunol Rev. 2015;267:30–55.

    Article  CAS  Google Scholar 

  37. Muraro PA, Robins H, Malhotra S, Howell M, Phippard D, Desmarais C, et al. T cell repertoire following autologous stem cell transplantation for multiple sclerosis. J Clin Invest. 2014;124:1168–72.

    Article  CAS  Google Scholar 

  38. Wagner UG, Koetz K, Weyand CM, Goronzy JJ. Perturbation of the T cell repertoire in rheumatoid arthritis. Proc Natl Acad Sci USA. 1998;95:14447–52.

    Article  CAS  Google Scholar 

  39. Wang L, Zhang W, Lin L, Li X, Saksena NK, Wu J, et al. A comprehensive analysis of the T and B lymphocytes repertoire shaped by HIV vaccines. Front Immunol. 2018;9:2194.

    Article  Google Scholar 

  40. Jiang Q, Zhao T, Zheng W, Zhou J, Wang H, Dong H, et al. Patient-shared TCRbeta-CDR3 clonotypes correlate with favorable prognosis in chronic hepatitis B. Eur J Immunol. 2018;48:1539–49.

    Article  CAS  Google Scholar 

  41. Yang G, Ou M, Chen H, Guo C, Chen J, Lin H, et al. Characteristic analysis of TCR beta-chain CDR3 repertoire for pre- and post-liver transplantation. Oncotarget. 2018;9:34506–19.

    PubMed  PubMed Central  Google Scholar 

  42. Hou X, Wang M, Lu C, Xie Q, Cui G, Chen J, et al. Analysis of the repertoire features of TCR beta chain CDR3 in human by high-throughput sequencing. Cell Physiol Biochem. 2016;39:651–67.

    Article  CAS  Google Scholar 

  43. Putintseva EV, Britanova OV, Staroverov DB, Merzlyak EM, Turchaninova MA, Shugay M, et al. Mother and child T cell receptor repertoires: deep profiling study. Front Immunol. 2013;4:463.

    Article  Google Scholar 

Download references

Authors’ contributions

SW, SL, JH and WP designed the study and patient specimen collection. SW, SL and WP completed the experimental process, data collection, literature search and generation of figures. SW, SL, XH, MBS; BB and MD wrote and edited the manuscript. HX and NH provided the experimental technical and laboratory facilities support. All authors reviewed the manuscript. All authors read and approved the final manuscript.

Acknowledgements

Not applicable.

Competing interests

WP, MBS, BB, MD, XH, SL, and JH are employees of iRepertoire. The other authors declare no competing interests.

Availability of data and materials

The data and materials supporting the conclusions of the current study are available from the corresponding author on reasonable request.

Consent for publication

Not applicable.

Ethics approval and consent to participate

The study design and recruiting of patients were approved by the New England Independent Review Board® (NEIRB) (IRB Number: 14-378). Written informed consent was obtained from all study participants. All experiments were performed according to the relevant guidelines and regulations.

Funding

The study did not receive any funding support.

Publisher’s Note

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

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Jian Han or Song Li.

Additional file

Additional file 1: Table S1.

The statistical results of sample (blood volume (µL)) calculation; Table S2. The numbers of shared CDR3 among the different people.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) 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

Wu, SG., Pan, W., Liu, H. et al. High throughput sequencing of T-cell receptor repertoire using dry blood spots. J Transl Med 17, 47 (2019). https://doi.org/10.1186/s12967-019-1796-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-019-1796-4

Keywords