GENTLE: a novel bioinformatics tool for generating features and building classifiers from T cell repertoire cancer data
BMC Bioinformatics volume 24, Article number: 32 (2023)
In the global effort to discover biomarkers for cancer prognosis, prediction tools have become essential resources. TCR (T cell receptor) repertoires contain important features that differentiate healthy controls from cancer patients or differentiate outcomes for patients being treated with different drugs. Considering, tools that can easily and quickly generate and identify important features out of TCR repertoire data and build accurate classifiers to predict future outcomes are essential.
This paper introduces GENTLE (GENerator of T cell receptor repertoire features for machine LEarning): an open-source, user-friendly web-application tool that allows TCR repertoire researchers to discover important features; to create classifier models and evaluate them with metrics; and to quickly generate visualizations for data interpretations. We performed a case study with repertoires of TRegs (regulatory T cells) and TConvs (conventional T cells) from healthy controls versus patients with breast cancer. We showed that diversity features were able to distinguish between the groups. Moreover, the classifiers built with these features could correctly classify samples (‘Healthy’ or ‘Breast Cancer’)from the TRegs repertoire when trained with the TConvs repertoire, and from the TConvs repertoire when trained with the TRegs repertoire.
The paper walks through installing and using GENTLE and presents a case study and results to demonstrate the application’s utility. GENTLE is geared towards any researcher working with TCR repertoire data and aims to discover predictive features from these data and build accurate classifiers. GENTLE is available on https://github.com/dhiego22/gentle and https://share.streamlit.io/dhiego22/gentle/main/gentle.py.
Identifying high-quality biomarkers in cancer data is a formidable challenge, but TCR repertoires have shown to be a useful source in surmounting this obstacle [1, 2]. TCRs are generated by a VDJ (variable, diversity, joining) recombination process that can generate a potential diversity of 1019 unique TCRs . This process yields two protein chains consisting primarily of alpha and beta chains; in approximately 10 percent of cases, they consist of gamma and delta chains . The beta chain uniquely contains the diversity (D) gene segment, which includes greater diversity–specifically in the CDR3 (Complementarity Determining Region 3) region where the D gene segment is located. Despite its complexity, a sample of these TCR repertoires generates substantial insights into immune system behavior, thus providing important information regarding the choice of therapy .
Many features of the repertoire, such as diversity, clonality, motifs, and network presentation, are often used as metrics to distinguish between cancer patients who may respond to a certain drug or simply between sick and healthy individuals [6,7,8]. Table 1 provides additional information about the terms introduced above. Originally used in ecology, diversity metrics can be adapted to characterize the TCR repertoire. The Shannon and Simpson indices are the most commonly used metrics , although other metrics, including Hill numbers, Pielou, and Gini, also serve as strategies [6, 9]. We previously identified  significant differences of 1 minus Pielou (referred to as clonality in the article; see their supplemental information file) and Simpson indices between control and mice with mammary tumors after immunotherapy. Network-based approaches may consider each TCR as a node and utilize distance metrics to build edges between the nodes. In another work , we demonstrated that a network representation could stratify control and transgenic mice blood samples. A longitudinal analysis of TCR repertoire networks showed variations in the density index of the network. The Levenshtein metric was employed to represent the editing distances used to build the networks. Using physicochemical motifs, Ostmeyer and collaborators  distinguished tumor tissue from patient-matched healthy tissue in the same organ by extracting 4-mers from CDR3 sequences and using their frequencies to build high-accuracy classifiers for breast cancer and colorectal cancer. Additionally, Wang and colleagues  identified 11 structural motifs to distinguish long-term survivors from short-term survivors with nasopharyngeal carcinoma.
Recent technological breakthroughs have advanced the entire field . In tandem with progress in immunotherapy, great improvements in NGS (Next Generation Sequencing) options now allow the sequencing of longer reads in large quantities and at an affordable price . These improvements have contributed to the gradually increasing availability of TCR-based data, while data analysis tools are still limited. TCR repertoire-oriented tools can be used in two types of analysis: low-level or high-level. Low-level tasks involve raw data processing, while high-level tasks extract information from processed data . IMGT/HighV-QUEST , IgBLAST , MiXCR , and MiTCR  are some of the most cited tools for low-level tasks, while Table 2 summarizes the most cited tools for high-level tasks and depicts their main features. The experimental design should guide the adequate tool selection according to a specific purpose. GENTLE generates features and allows users to create classifiers and evaluate their predictive power on experimental samples, thus meeting the need for fast and easy-to-use data analysis tools developed using–and based on–easily accessible sources .
In summary, we have developed GENTLE, the first tool to offer a machine learning pipeline for TCR repertoire cancer data: it enables the generation of features, normalization methods, feature selection algorithms, classifiers construction, and evaluation metrics for internal and external validation.
The code for GENTLE was written in Python 3 and can be run on version 3.9 or higher. We used Streamlit to construct the GUI (Graphical User Interface), Pandas for data manipulation, Plotly for data visualization, and scikit-learn for many machine learning algorithms . Streamlit was chosen for its speed and simplicity in implementing an application front-end and for easily and freely deploying and sharing the application with the community. Streamlit also makes available some functions that do not run unless a specific parameter has changed. In light of this, we incorporated these functions to avoid redundant calculations and to keep the program running fast. The source code of GENTLE is available on https://github.com/dhiego22/gentle; the README documentation is succinct and clear for users, and the program can be easily installed with virtualenv or docker.
The input into GENTLE is a.csv file format. For files that surpass the maximum size supported (200 MB), the.csv file can be zipped and uploaded in a zip format. The file must be a dataframe in which the rows represent the samples of the experiment, and the columns represent the TCRs with one additional column representing the label of each sample (e.g., case/control). The values should represent the counts of the TCRs in the samples. We provide examples of the input data in the Github repository.
After uploading the data, four options will appear in the sidebar; these options are the feature dimensions to be analyzed. It is important to emphasize that when uploading a different dataframe, one should always erase the cache from the options menu in the top-right corner of the screen. The Diversity metric calculates popular diversity measurements widely used in ecology such as richness, Shannon, Simpson, inverse Simpson, Pielou, one minus Pielou, hill numbers and Gini indices. The Network metric will use the TCR sequences as nodes and calculate a Levenshtein distance of two, according to , to create edges between the nodes. After creating the networks, features like the number of nodes and edges, density, clustering coefficient, transitivity and connected components are calculated. The Motif metric calculates the frequency of contiguous letters specifically, 2-mers, 3-mers and 4-mers. Finally, GENTLE gives the option to use six different dimensionality reduction methods: PCA (Principal Component Analysis) , t-SNE (t-distributed Stochastic Neighbor Embedding) , UMAP (Uniform Manifold Approximation and Projection) , ICA (Independent Component Analysis) , SVD (Singular Value Decomposition) , and ISOMAP (Isometric Mapping) . In addition, when exporting the features as a dataframe, GENTLE offers three normalizing options in the sidebar. The first option is standard normalization, which converts the data to an average value of 0.0 and a standard deviation of 1.0. The second option is min–max normalization, which linearly converts the data such that the minimum value is − 1.0 and the maximum value is 1.0. The third option is robust scaler, which subtracts the median value and linearly scales the data based on the interquartile range. Normalization is optional, but it can significantly impact the algorithmic performance .
GENTLE provides four feature selection methods: Pearson’s correlation, Ridge, XGBoost, and mRMR. Upon selecting the methods, a dataframe will be created with a rank of the features with the greatest predictive power, according to each method, where zero means the feature was not selected by the method; one means the feature was selected as the most predictive; two, the second-most predictive; and so on. It is important to emphasize that Pearson's method considers only one feature at a time when defining its predictive power; this way, the two most predictive features can be so correlated that their combination may not improve a classifier's predictive power when trained together. In contrast, the other methods will consider the combination of the features, which means the two most predictive features are the combination of two features that will produce the most predictive classifier if trained together. In addition, for visualization purposes, one can choose two featuresto display a 2D scatter plot, or three featuresto display a rotating 3D scatter plot.
For the next step, one can perform the classification and validate the predictive power of the selected features. Four classifiers can be chosen: GNB (Gaussian Naive Bayes), LDA (Linear Discriminant Analysis), LR (Logistic Regression), and DT (Decision tree). A radar plot will be generated, representing the five main scoring methods for classifiers: accuracy, precision, recall, F1, and AUC (Area Under Curve) ROC (Receiver Operating Characteristics) curve.
Finally, one can upload another dataset for external validation purposes. This ultimate step produces a confusion matrix and a radar plot with accuracy, precision, recall, F1-score and AUC ROC scoring methods as explained above.
Each dataframe generated can be downloaded (in.csv file format), along with the networks, the charts (in.png image format), and the classifier model (in pickle format). The networks created can also be visualized. GENTLE can be used for educational purposes due to its user-friendly interface and simplicity. This tool is particularly useful in providing fast feedback when analyzing a TCR repertoire and its features. Figure 1 summarizes the main steps for using GENTLE and understanding its capabilities. There is also a concise walkthrough with screenshots in 11 steps available in the Additional file 1—A Walkthrough of GENTLE and the methods from GENTLE’s flow are summarized in the Additional file 2—Diversity metrics, network metrics, dimensionality reduction, classifiers and scoring metrics.
Algorithms for feature selection and classification
Pearson’s correlation is a widely applied strategy to find the features most related to a target and to eliminate any redundancies between features . It provides a rank of the most correlated features according to a given target; it does not consider a set of features for prediction but rather considers features independently. For example, if feature A and feature B have the highest correlation values but are alike, we can achieve the same results with only one. mRMR (minimum redundancy maximum relevance) circumvents this issue by selecting features with high predictive power which are simultaneously different from one another; i.e., it selects the smallest relevant subset of features . Ridge regression (also known as L2 regularization) is the ideal method to tackle the overfitting issue, given its regularization use. This method determines variables with zero effect on the data without wasting predictive information. XGBoost (Extreme Gradient Boosting) uses ensembles of decision tree methods, like gradient boosting, to estimate the importance of the features when training a predictive model. It is an appropriate method when working with large datasets; moreover, it can reflect complex interactions between the features . Here, we adopt algorithms belonging to the three main classes of feature selection methods: filter-based represented by Pearson’s correlation, embedded represented by Ridge and XGBoost, and wrapper represented by mRMR.
To build the classifiers, we chose four methods based on their transparency, speed, and capacity to compete with sophisticated methods (e.g., neural networks, ensemble methods) and degree of relevance in the literature. In terms of transparency, we are referring to internal processes used and the various weighted factors that remain unknown, commonly known as a ‘black box’ algorithm. Naive Bayes could classify between malignant and benign breast cancer using the Breast Cancer Wisconsin Data Set, providing fast and accurate results . Linear Discriminant Analysis is widely used in the biomedical field to separate groups according to disease and response to treatments. It is an efficient approach to dealing with the small sample size problem . Logistic Regression was applied to calculate the log-likelihood ratio of radiological response to anti-PD1 (programmed cell death protein 1) therapy TCR repertoire data of patients with metastatic melanoma . Decision Tree was the base method used to investigate features in 15 patients with NSCLC (non-small cell lung cancer) using a combination of exome, transcriptome, and TCR repertoire data .
In summary, the feature selection methods and the classifier models can tackle diverse issues, perform rapidly, and deliver competitive and accurate results compared to other state-of-the-art approaches.
Results and discussion
To perform a case study workflow, we used a public dataset from the TCRdb that can be accessed from the link http://bioinfo.life.hust.edu.cn/TCRdb/#/download (access the project PRJNA297261). The website contains seven TCR repertoire projects, but only two of them contain at least two conditions: ‘Healthy’ or ‘Breast Cancer’. From these two projects, the immunoSEQ20 project is unbalanced, as it contains 60 ‘Breast Cancer’ samples and only three ‘Healthy’ samples. Considering that, we decided to move forward with the PRJNA297261 project, which was balanced. Notice that the dataset here is only used for demonstration purposes and not for its medical/biological merit. We did not produce this dataset. This database provides preprocessed dataframes that require less processing to fit GENTLE’s input. The original project labels each sample with the condition of ‘Breast Cancer’ or ‘Healthy’ tissues; the project also labels each sample with the cell type of ‘TRegs’ or ‘TConvs’. We made available on Github the script in which we processed the data extracted from the website, turning it into a dataframe that serves as input for GENTLE. This script splits the original project (PRJNA297261) into two dataframes; one with only the TRegs samples, and the other one with only the TConvs samples. Each dataframe labels its samples as ‘Healthy’ or ‘Breast Cancer’. The processed data is summarized in Table 3 and is available on the Github page. The original data comprises TCR beta chain repertoire of regulatory and conventional T cells in peripheral blood from breast cancer patients and healthy individuals. TRegs are important for the regulation of immune response, including TConvs, which can differentiate into effector cells and respond to non-self antigens. Although TRegs and TConvs have different functions they can descend from common clones .
In this case study, we identified the features that can separate the healthy from the sick samples in both datasets, and used them to build classifiers. To avoid overfitting, we limited the number of features according to the minimum number of samples found in each category minus one . When building the classifiers, we considered two features: one less from the three healthy patients from the TRegs and in both categories from the TConvs. We trained one model with the TRegs dataset using the TConvs dataset as the test/unseen data; we also made a switch by training one model with the TConvs dataset and using the TRegs dataset as the holdout/unseen data, thus allowing us to analyze the predictive power of each feature for the TCR repertoire.
By analyzing all four dimensions (diversity, network, motif, and dimensionality reduction) using the TConvs dataset, we could construct scatter plots that portray a separation between the healthy and the sick samples. The features in the x and y-axis of each scatter plot were chosen based on the feature selection method’s choice for the most predictive features. Both Simpson and Shannon indices from the diversity features were able to separate the healthy and sick samples (Fig. 2A), wherein the sick samples had higher values for the Shannon index and lower values for the Simpson index; the opposite was true for the healthy samples. The density and the number of arrows could also accurately separate the samples (Fig. 2B), in which the healthy samples had lower density values and more arrows from the built networks than the samples with breast cancer. The frequencies of the motifs ‘VS’ (valine followed by a serine) and ‘SV’ (serine followed by a valine) could likewise separate the samples (Fig. 2C); both frequencies were more common in healthy patients than in sick patients. Many dimensional reduction methods were also able to separate the samples, but we depicted an interesting scenario in which the IC1 feature did not, but the IC2 feature completely separated the samples (Fig. 2D). The features generated from these repertoires clearly distinguished the sick and healthy samples as seen in the scatter plots in Fig. 2A–D.
The TRegs dataset could not distinguish the healthy from the sick samples as the TConvs did. The only feature that perfectly separated the samples was the Shannon index (Fig. 2E). Again, density and number of arrows proved to be the most predictive features from the built networks (Fig. 2F). Although they did not completely separate all the samples, they showed a slight tendency towards separation. The combination of the motifs ‘GG’ (two guanines together) and ‘SQ’ (serine followed by a glutamine) shown in (Fig. 2G), can be considered to separate the samples. The dimension reduction methods had difficulty accomplishing this task; PC1 and PC2 are depicted in (Fig. 2H).
It is important to emphasize that the feature selection algorithms chose the features shown in each scatter plot as the ones with the highest predictive power. Many features from the diversity dimension (omitted here) could distinguish the ‘Healthy’ and ‘Breast Cancer’ samples from the TConvs dataset, while even the most predictive features from the network, motif, and dimension reduction had difficulty distinguishing the samples from the TRegs dataset.
Based on the exploratory analysis we performed (shown in Fig. 2), we decided to build and evaluate the classifiers only with the selected features by the feature selection methods. The Shannon and Simpson indexes were the only features able to build classifiers with high scores of the internal validation and which could classify all the samples correctly from the holdout dataset (see Fig. 3). For the TConvs dataset, all the built models had scores close to 1 with the internal validation, and all of them could predict the samples from the TRegs dataset perfectly. The only model trained with the TRegs dataset that could classify the TConvs dataset perfectly was the Decision Tree.
In sum, this analysis demonstrates the similarity between TRegs and TConvs when analyzing their diversity features, and it portrays how these features can be predictive when analyzing TCR repertoires of breast cancer. Many studies–to cite a few, [2, 51, 52], and –have opted to explore TCR repertoires in their diversity dimension (specifically, the Shannon index). In corroboration with our results, diversity is an essential feature for TCR repertoire analyses and should not be omitted from any relevant discussion. Although this analysis was performed with only a few samples, we believe that it provides an example of the strength of the tool. With the addition of new datasets, more analyses can be done to strengthen the underlying hypothesis regarding information in the repertoire. Both cell types (TRegs and TConvs) are influenced by the same signaling pathways that dictate their development, differentiation, and function . Therefore, understanding TRegs’ and TConvs’ similarities can be a promising pathway to new therapeutic approaches.
This paper introduces GENTLE, a platform designed to help researchers easily and swiftly analyze their TCR repertoire data. GENTLE provides visualization capabilities and a user-friendly interface.It serves as the first graphical web tool to to incorporate feature selection methods to identify important features built within the TCR repertoire. It also makes available a set of diverse machine learning methods to generate models for classification purposes. The platform makes it possible to compare the performance of the classifiers through the main evaluation metrics for binary classification and also offers metrics for external validation. All data generated by GENTLE can be downloaded for further analysis. As an open-source web application, GENTLE provides researchers tools to analyze data efficiently and to be able to extract biomarkers and build classifiers that could positively improve treatment prospects across healthcare. Our case study showed that diversity features, such as the Shannon and Simpson indices, can be important biomarkers for healthy and sick patients with breast cancer when analyzing their repertoires of TRegs and Tconvs. For future works, we will add features generated from time series data, as some insights can only be gleaned upon analyzing changes in the repertoires over time. Moreover, we will add more options for classifiers and metrics of feature selection aligned more closely with time series data.
Availability and requirements
Project name: GENTLE. Project home page: https://github.com/dhiego22/gentle & https://share.streamlit.io/dhiego22/gentle/main/gentle.py. Operating system(s): Linux, Windows, Mac. Programming language: Python 3.9 + . Other requirements: Streamlit, Plotly, Pandas, Sckit-learn. License: MIT License. Any restrictions to use by non-academics: None.
Availability of data and materials
The data used in this work was extracted from the TCRdb, a database that contains a plethora of TCR repertoires of many cancer and cell types. The data can be accessed through the link http://bioinfo.life.hust.edu.cn/TCRdb/#/download (project PRJNA297261).
GENerator of T cell receptor repertoire features for machine LEarning
T cell receptor
Variable, diversity, joining genes
Graphical user interface
Regulatory T cells
Conventional T cells
Next generation sequencing
Complementarity determining region 3
Principal component analysis
T-distributed stochastic neighbor embedding
Uniform Manifold approximation and projection
Independent component analysis
Singular value decomposition
Gaussian Naive Bayes
Linear discriminant analysis
Area under curve
Receiver operating characteristics
Non-small cell lung cancer
Extreme gradient boosting
Minimum redundancy maximum relevance
Programmed cell death protein 1
Kumagai S, Togashi Y, Kamada T, Sugiyama E, Nishinakamura H, Takeuchi Y, et al. The PD-1 expression balance between effector and regulatory T cells predicts the clinical efficacy of PD-1 blockade therapies. Nat Immunol. 2020;21:1346–58. https://doi.org/10.1038/s41590-020-0769-3.
Zhao J, Chen AX, Gartrell RD, Silverman AM, Aparicio L, Chu T, et al. Immune and genomic correlates of response to anti-PD-1 immunotherapy in glioblastoma. Nat Med. 2019;25:462–9. https://doi.org/10.1038/s41591-019-0349-y.
Pai JA, Satpathy AT. High-throughput and single-cell T cell receptor sequencing technologies. Nat Methods. 2021;18:881–92. https://doi.org/10.1038/s41592-021-01201-8.
Girardi M. Immunosurveillance and immunoregulation by gammadelta T cells. J Invest Dermatol. 2006;126:25–31. https://doi.org/10.1038/sj.jid.5700003.
Arnaout RA, Prak ETL, Schwab N, Rubelt F. Adaptive immune receptor repertoire community. The future of blood testing is the immunome. Front Immunol. 2021;12:626793. https://doi.org/10.3389/fimmu.2021.626793.
Chiffelle J, Genolet R, Perez MA, Coukos G, Zoete V, Harari A. T-cell repertoire analysis and metrics of diversity and clonality. Curr Opin Biotechnol. 2020;65:284–95. https://doi.org/10.1016/j.copbio.2020.07.010.
Valkiers S, de Vrij N, Gielis S, Verbandt S, Ogunjimi B, Laukens K, et al. Recent advances in T-cell receptor repertoire analysis: bridging the gap with multimodal single-cell RNA sequencing. ImmunoInformatics. 2022;5:100009. https://doi.org/10.1016/j.immuno.2022.100009.
Katayama Y, Kobayashi TJ. Comparative study of repertoire classification methods reveals data efficiency of -mer feature extraction. Front Immunol. 2022;13:797640. https://doi.org/10.3389/fimmu.2022.797640.
Kidman J, Principe N, Watson M, Lassmann T, Holt RA, Nowak AK, et al. Characteristics of TCR repertoire associated with successful immune checkpoint therapy responses. Front Immunol. 2020;11:587014. https://doi.org/10.3389/fimmu.2020.587014.
Philip H, Snir T, Gordin M, Shugay M, Zilberberg A, Efroni S. A T cell repertoire timestamp is at the core of responsiveness to CTLA-4 blockade. IScience. 2021;24:102100. https://doi.org/10.1016/j.isci.2021.102100.
Priel A, Gordin M, Philip H, Zilberberg A, Efroni S. Network representation of T-cell repertoire—a novel tool to analyze immune response to cancer formation. Front Immunol. 2018;9:2913. https://doi.org/10.3389/fimmu.2018.02913.
Ostmeyer J, Christley S, Toby IT, Cowell LG. Biophysicochemical motifs in T-cell receptor sequences distinguish repertoires from tumor-infiltrating lymphocyte and adjacent healthy tissue. Cancer Res. 2019;79:1671–80. https://doi.org/10.1158/0008-5472.CAN-18-2292.
Wang G, Mudgal P, Wang L, Shuen TWH, Wu H, Alexander PB, et al. TCR repertoire characteristics predict clinical response to adoptive CTL therapy against nasopharyngeal carcinoma. Oncoimmunology. 2021;10:1955545. https://doi.org/10.1080/2162402X.2021.1955545.
Naylor K, Li G, Vallejo AN, Lee W-W, Koetz K, Bryl E, et al. The Influence of age on T cell generation and TCR diversity. J Immunol. 2005;174:7446–52. https://doi.org/10.4049/jimmunol.174.11.7446.
Mahe E, Pugh T, Kamel-Reid S. T cell clonality assessment: past, present and future. J Clin Pathol. 2018;71:195–200. https://doi.org/10.1136/jclinpath-2017-204761.
Bagaev DV, Vroomans RMA, Samir J, Stervbo U, Rius C, Dolton G, et al. VDJdb in 2019: database extension, new analysis infrastructure and a T-cell receptor motif compendium. Nucleic Acids Res. 2020;48:D1057–62. https://doi.org/10.1093/nar/gkz874.
Joshi K, Milighetti M, Chain BM. Application of T cell receptor (TCR) repertoire analysis for the advancement of cancer immunotherapy. Curr Opin Immunol. 2022;74:1–8. https://doi.org/10.1016/j.coi.2021.07.006.
Zhang Y, Yang X, Zhang Y, Zhang Y, Wang M, Ou JX, et al. Tools for fundamental analysis functions of TCR repertoires: a systematic comparison. Brief Bioinform. 2020;21:1706–16. https://doi.org/10.1093/bib/bbz092.
Alamyar E, Duroux P, Lefranc M-P, Giudicelli V. IMGT® Tools for the Nucleotide Analysis of Immunoglobulin (IG) and T Cell Receptor (TR) V-(D)-J Repertoires, Polymorphisms, and IG Mutations: IMGT/V-QUEST and IMGT/HighV-QUEST for NGS. Methods in Molecular BiologyTM. 2012;569–604. https://doi.org/10.1007/978-1-61779-842-9_32.
Ye J, Ma N, Madden TL, Ostell JM. IgBLAST: an immunoglobulin variable domain sequence analysis tool. Nucleic Acids Res. 2013;41:W34-40. https://doi.org/10.1093/nar/gkt382.
Bolotin DA, Poslavsky S, Mitrophanov I, Shugay M, Mamedov IZ, Putintseva EV, et al. MiXCR: software for comprehensive adaptive immunity profiling. Nat Methods. 2015;12:380–1. https://doi.org/10.1038/nmeth.3364.
Bolotin DA, Shugay M, Mamedov IZ, Putintseva EV, Turchaninova MA, Zvyagin IV, et al. MiTCR: software for T-cell receptor sequencing data analysis. Nat Methods. 2013;10:813–4. https://doi.org/10.1038/nmeth.2555.
Henry VJ, Bandrowski AE, Pepin A-S, Gonzalez BJ, Desfeux A. OMICtools: an informative directory for multi-omic data analysis. Database. 2014. https://doi.org/10.1093/database/bau069.
Pavlović M, Scheffer L, Motwani K, Kanduri C, Kompova R, Vazov N, et al. The immuneML ecosystem for machine learning analysis of adaptive immune receptor repertoires. Nat Mach Intell. 2021;3:936–44. https://doi.org/10.1038/s42256-021-00413-z.
Sturm G, Szabo T, Fotakis G, Haider M, Rieder D, Trajanoski Z, et al. Scirpy: a Scanpy extension for analyzing single-cell T-cell receptor-sequencing data. Bioinformatics. 2020;36:4817–8. https://doi.org/10.1093/bioinformatics/btaa611.
Popov A, ivan-immunomind, MVolobueva, Nazarov VI, Bot I, Rumynskiy E, et al. Immunarch 0.6.8: an R package for painless bioinformatics analysis of T-cell and B-cell immune repertoires. Zenodo; 2022. 10.5281/ZENODO.3367200.
Morin A, Kwan T, Ge B, Letourneau L, Ban M, Tandre K, et al. Immunoseq: the identification of functionally relevant variants through targeted capture and sequencing of active regulatory regions in human immune cells. BMC Med Genomics. 2016;9:59. https://doi.org/10.1186/s12920-016-0220-7.
Vander Heiden JA, Yaari G, Uduman M, Stern JNH, O’Connor KC, Hafler DA, et al. pRESTO: a toolkit for processing high-throughput sequencing raw reads of lymphocyte receptor repertoires. Bioinformatics. 2014;30:1930–2. https://doi.org/10.1093/bioinformatics/btu138.
Gupta NT, Vander Heiden JA, Uduman M, Gadala-Maria D, Yaari G, Kleinstein SH. Change-O: a toolkit for analyzing large-scale B cell immunoglobulin repertoire sequencing data. Bioinformatics. 2015;31:3356–8. https://doi.org/10.1093/bioinformatics/btv359.
Shugay M, Bagaev DV, Turchaninova MA, Bolotin DA, Britanova OV, Putintseva EV, et al. VDJtools: Unifying post-analysis of T cell receptor repertoires. PLoS Comput Biol. 2015;11:e1004503. https://doi.org/10.1371/journal.pcbi.1004503.
Schattgen SA, Guion K, Crawford JC, Souquette A, Barrio AM, Stubbington MJT, et al. Integrating T cell receptor sequences and transcriptional profiles by clonotype neighbor graph analysis (CoNGA). Nat Biotechnol. 2022;40:54–63. https://doi.org/10.1038/s41587-021-00989-2.
Borcherding N, Bormann NL, Kraus G. scRepertoire: An R-based toolkit for single-cell immune receptor analysis. F1000Res. 2020;9:47. https://doi.org/10.12688/f1000research.22139.2.
Weber CR, Rubio T, Wang L, Zhang W, Robert PA, Akbar R, et al. Reference-based comparison of adaptive immune receptor repertoires. Cell Rep Methods. 2022;2:100269. https://doi.org/10.1016/j.crmeth.2022.100269.
Garreta R, Moncecchi G. Learning Scikit-Learn: Machine Learning in Python. Packt Pub Limited; 2013. https://jmlr.csail.mit.edu/papers/v12/pedregosa11a.html. Accessed 21 Nov 2022.
Wold S, Esbensen K, Geladi P. Principal component analysis. Chemometrics Intellig Lab Syst. 1987;2:37–52. https://doi.org/10.1016/0169-7439(87)80084-9.
Soni J, Prabakar N, Upadhyay H. Visualizing high-dimensional data using t-distributed stochastic neighbor embedding algorithm. Princ Data Sci. 2020. https://doi.org/10.1007/978-3-030-43981-1_9.
Sainburg T, McInnes L, Gentner TQ. Parametric UMAP embeddings for representation and semisupervised learning. Neural Comput. 2021;33:2881–907. https://doi.org/10.1162/neco_a_01434.
Hyvärinen A, Oja E. Independent component analysis: algorithms and applications. Neural Netw. 2000;13:411–30. https://doi.org/10.1016/s0893-6080(00)00026-5.
Van Loan CF. Generalizing the Singular Value Decomposition. SIAM J Numer Anal. 2006. https://doi.org/10.1137/0713009.10.1137/0713009.
Tenenbaum JB, de Silva V, Langford JC. A global geometric framework for nonlinear dimensionality reduction. Science. 2000;290:2319–23. https://doi.org/10.1126/science.290.5500.2319.
Singh D, Singh B. Investigating the impact of data normalization on classification performance. Appl Soft Comput. 2020;97:105524. https://doi.org/10.1016/j.asoc.2019.105524.
Gnanadesikan R, Kettenring JR, Tsao SL. Weighting and selection of variables for cluster analysis. J Classif. 1995;12:113–36. https://doi.org/10.1007/bf01202271.
Ding C, Peng H. Minimum redundancy feature selection from microarray gene expression data. J Bioinform Comput Biol. 2005;3:185–205. https://doi.org/10.1142/s0219720005001004.
Alsahaf A, Petkov N, Shenoy V, Azzopardi G. A framework for feature selection through boosting. Expert Syst Appl. 2022;187:115895. https://doi.org/10.1016/j.eswa.2021.115895.
Wood A, Shpilrain V, Najarian K, Kahrobaei D. Private naive bayes classification of personal biomedical data: application in cancer data analysis. Comput Biol Med. 2019;105:144–50. https://doi.org/10.1016/j.compbiomed.2018.11.018.
Sharma A, Paliwal KK. Linear discriminant analysis for the small sample size problem: an overview. Int J Mach Learn Cybern. 2015;6:443–54. https://doi.org/10.1007/s13042-013-0226-9.
Valpione S, Mundra PA, Galvani E, Campana LG, Lorigan P, De Rosa F, et al. The T cell receptor repertoire of tumor infiltrating T cells is predictive and prognostic for cancer survival. Nat Commun. 2021. https://doi.org/10.1038/s41467-021-24343-x.
Jia Q, Wu W, Wang Y, Alexander PB, Sun C, Gong Z, et al. Local mutational diversity drives intratumoral immune heterogeneity in non-small cell lung cancer. Nat Commun. 2018;9:5361. https://doi.org/10.1038/s41467-018-07767-w.
Wolf KJ, Emerson RO, Pingel J, Buller RM, DiPaolo RJ. Conventional and regulatory CD4+ T cells that share identical TCRs are derived from common clones. PLoS ONE. 2016;11:e0153705. https://doi.org/10.1371/journal.pone.0153705.
Hua J, Xiong Z, Lowey J, Suh E, Dougherty ER. Optimal number of features as a function of sample size for various classification rules. Bioinformatics. 2005;21:1509–15. https://doi.org/10.1093/bioinformatics/bti171.
Ji S, Li J, Chang L, Zhao C, Jia R, Tan Z, et al. Peripheral blood T-cell receptor repertoire as a predictor of clinical outcomes in gastrointestinal cancer patients treated with PD-1 inhibitor. Clin Transl Oncol. 2021;23:1646–56. https://doi.org/10.1007/s12094-021-02562-4.
Cui J-H, Lin K-R, Yuan S-H, Jin Y-B, Chen X-P, Su X-K, et al. TCR repertoire as a novel indicator for immune monitoring and prognosis assessment of patients with cervical Cancer. Front Immunol. 2018;9:2729. https://doi.org/10.3389/fimmu.2018.02729.
Chaara W, Gonzalez-Tort A, Florez L-M, Klatzmann D, Mariotti-Ferrandiz E, Six A. RepSeq data representativeness and robustness assessment by Shannon entropy. Front Immunol. 2018;9:1038. https://doi.org/10.3389/fimmu.2018.01038.
Camirand G, Riella LV. Treg-centric view of immunosuppressive drugs in transplantation: a balancing act. Am J Transplant. 2017;17:601–10. https://doi.org/10.1111/ajt.14029.
The authors thank Isa Goldberg, Or Malca, Beatriz Stransky, Thiago Felipe, Emmanuel Barbosa, Leonardo Capitani, students and professors of the Federal University of Rio Grande do Norte and Bar-Ilan University for the encouragement, support, discussions and insights generated along the development of this work.
This work is supported in part by funds from the Brazilian Funding agency CAPES—National Coordination of High Education Personnel Formation Programs (Grants Numbers 88887.161820/2017-0, 88887.469283/2019-00 and 88887.600071/2021-0). The APC was funded by the Federal University of Rio Grande do Norte. This research was supported by NPAD/UFRN.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
César Rennó-Costa, Sol Efroni share senior authorship
About this article
Cite this article
Andrade, D.S., Terrematte, P., Rennó-Costa, C. et al. GENTLE: a novel bioinformatics tool for generating features and building classifiers from T cell repertoire cancer data. BMC Bioinformatics 24, 32 (2023). https://doi.org/10.1186/s12859-023-05155-w