The complete 685-kilobase DNA sequence of the human beta T cell receptor locus

A multiomic atlas identifies a treatment-resistant, bone marrow progenitor-like cell population in T cell acute lymphoblastic leukemia

Refractoriness to initial chemotherapy and relapse after remission are the main obstacles to curing T cell acute lymphoblastic leukemia (T-ALL). While tumor heterogeneity has been implicated in treatment failure, the cellular and genetic factors contributing to resistance and relapse remain unknown. Here we linked tumor subpopulations with clinical outcome, created an atlas of healthy pediatric hematopoiesis and applied single-cell multiomic analysis to a diverse cohort of 40 T-ALL cases. We identified a bone marrow progenitor (BMP)-like leukemia subpopulation associated with treatment failure and poor overall survival. The single-cell-derived molecular signature of BMP-like blasts predicted poor outcome across multiple subtypes of T-ALL and revealed that NOTCH1 mutations additively drive T-ALL blasts away from the BMP-like state. Through in silico and in vitro drug screenings, we identified a therapeutic vulnerability of BMP-like blasts to apoptosis-inducing agents including venetoclax. Collectively, our study establishes multiomic signatures for rapid risk stratification and targeted treatment of high-risk T-ALL.

Bulk T cell repertoire sequencing (TCR-Seq) is a powerful technology for understanding inflammation-mediated diseases

Multiple alterations in the T cell repertoire were identified in many chronic inflammatory diseases such as inflammatory bowel disease, multiple sclerosis, and rheumatoid arthritis, suggesting that T cells might, directly or indirectly, be implicated in these pathologies. This has sparked a deep interest in characterizing disease-associated T cell clonotypes as well as in identifying and quantifying their contribution to the pathophysiology of different autoimmune and inflammation-mediated diseases. Bulk T cell repertoire sequencing (TCR-Seq) has emerged as a powerful method to profile the T cell repertoire of a sample in a high throughput fashion. Given the increasing utilization of TCR-Seq, we aimed here to provide a comprehensive, up-to-date review of the method, its extensions, and its ability to investigate chronic and autoimmune diseases. Specifically, we started by introducing the immunological basis of TCR repertoire generation and features, followed by discussing different experimental approach to perform TCR-Seq, then we describe different methods and frameworks for analyzing the generated datasets. Subsequently, different experimental techniques for investigating the antigenicity of T cell clonotypes are described. Lastly, we discuss recent studies that utilized TCR-Seq to understand different inflammation-mediated diseases, discuss fallbacks of the technology and potential future directions to overcome current limitations.

HTS and scRNA-seq revealed that the location and RSS quality of the mammalian TRBV and TRBJ genes impact biased rearrangement

The quality of Recombination signal sequences (RSSs), location, and genetics of mammalian V, D, and J genes synergistically affect the recombination frequency of genes; however, the specific regulatory mechanism and efficiency have not been elucidated. By taking advantage of single-cell RNA-sequencing (scRNA-seq) and high-throughput sequencing (HTS) to investigate V(D)J rearrangement characteristics in the CDR3 repertoire, we found that the distal and proximal V genes (or J genes) "to D" gene were involved in rearrangement significantly more frequently than the middle V genes (or J genes) in the TRB locus among various species, including Primates (human and rhesus monkey), Rodentia (BALB/c, C57BL/6, and Kunming mice), Artiodactyla (buffalo), and Chiroptera (Rhinolophus affinis). The RSS quality of the V and J genes affected their frequency in rearrangement to varying degrees, especially when the V-RSSs with recombination signal information content (RIC) score < -45 significantly reduced the recombination frequency of the V gene. The V and J genes that were "away from D" had the dual advantages of recombinant structural accessibility and relatively high-quality RSSs, which promoted their preferential utilization in rearrangement. The quality of J-RSSs formed during mammalian evolution was apparently greater than that of V-RSSs, and the D-J distance was obviously shorter than that of V-D, which may be one of the reasons for guaranteeing that the "D-to-J preceding V-to-DJ rule" occurred when rearranged. This study provides a novel perspective on the mechanism and efficiency of V-D-J rearrangement in the mammalian TRB locus, as well as the biased utilization characteristics and application of V and J genes in the initial CDR3 repertoire.

Phenotypic and functional characterization of posoleucel, a multivirus-specific T cell therapy for the treatment and prevention of viral infections in immunocompromised patients

BACKGROUND

Deficits in T cell immunity translate into increased risk of severe viral infection in recipients of solid organ and hematopoietic cell transplants. Thus, therapeutic strategies that employ the adoptive transfer of virus-specific T cells are being clinically investigated to treat and prevent viral diseases in these highly immunocompromised patients. Posoleucel is an off-the-shelf multivirus-specific T cell investigational product for the treatment and prevention of infections due to adenovirus, BK virus, cytomegalovirus, Epstein-Barr virus, human herpesvirus 6 or JC virus.

METHODS

Herein we perform extensive characterization of the phenotype and functional profile of posoleucel to illustrate the cellular properties that may contribute to its in vivo activity.

RESULTS AND CONCLUSIONS

Our results demonstrate that posoleucel is enriched for central and effector memory CD4+ and CD8+ T cells with specificity for posoleucel target viruses and expressing a broad repertoire of T cell receptors. Antigen-driven upregulation of cell-surface molecules and production of cytokine and effector molecules indicative of proliferation, co-stimulation, and cytolytic potential demonstrate the specificity of posoleucel and its potential to mount a broad, polyfunctional, and effective Th1-polarized antiviral response upon viral exposure. We also show the low risk for off-target and nonspecific effects as evidenced by the enrichment of posoleucel in memory T cells, low frequency of naive T cells, and lack of demonstrated alloreactivity in vitro. The efficacy of posoleucel is being explored in four placebo-controlled clinical trials in transplant recipients to treat and prevent viral infections (NCT05179057, NCT05305040, NCT04390113, NCT04605484).

The reverse TRBV30 gene of mammals: a defect or superiority in evolution?

At the 3' end of the C2 gene in the mammalian TRB locus, a distinct reverse TRBV30 gene (named TRBV31 in mice) has been conserved throughout evolution. In the fully annotated TRB locus of 14 mammals (including six orders), we observed noteworthy variations in the localization and quality of the reverse V30 genes and Recombination Signal Sequences (RSSs) in the gene trees of 13 mammals. Conversely, the forward V29 genes and RSSs were generally consistent with the species tree of their corresponding species. This finding suggested that the evolution of the reverse V30 gene was not synchronous and likely played a crucial role in regulating adaptive immune responses. To further investigate this possibility, we utilized single-cell TCR sequencing (scTCR-seq) and high-throughput sequencing (HTS) to analyze TCRβ CDR3 repertoires from both central and peripheral tissues of Primates (Homo sapiens and Macaca mulatta), Rodentia (Mus musculus: BALB/c, C57BL/6, and Kunming mice), Artiodactyla (Bos taurus and Bubalus bubalis), and Chiroptera (Rhinolophus affinis and Hipposideros armige). Our investigation revealed several novel observations: (1) The reverse V30 gene exhibits classical rearrangement patterns adhering to the '12/23 rule' and the 'D-J rearrangement preceding the V-(D-J) rearrangement'. This results in the formation of rearranged V30-D2J2, V30-D1J1, and V30-D1J2. However, we also identified 'special rearrangement patterns' wherein V30-D rearrangement preceding D-J rearrangement, giving rise to rearranged V30-D2-J1 and forward Vx-D2-J. (2) Compared to the 'deletional rearrangement' (looping out) of forward V1-V29 genes, the reverse V30 gene exhibits preferential utilization with 'inversional rearrangement'. This may be attributed to the shorter distance between the V30 gene and D gene and the 'inversional rearrangement' modes. In summary, in the mammalian TRB locus, the reverse V30 gene has been uniquely preserved throughout evolution and preferentially utilized in V(D)J recombination, potentially serving a significant role in adaptive immunity. These results will pave the way for novel and specialized research into the mechanisms, efficiency, and function of V(D)J recombination in mammals.

An ancient enzyme finds a new home: Prevalence and neofunctionalization of trypsin in marine phytoplankton

Trypsin is an ancient protease best known as a digestive enzyme in animals, and traditionally believed to be absent in plants and protists. However, our recent studies have revealed its wide presence and important roles in marine phytoplankton. Here, to gain a better understanding on the importance of trypsin in phytoplankton, we further surveyed the distribution, diversity, evolution and potential ecological roles of trypsin in global ocean phytoplankton. Our analysis indicated that trypsin is widely distributed both taxonomically and geographically in marine phytoplankton. Furthermore, by systematic comparative analyses we found that algal trypsin could be classified into two subfamilies (trypsin I and trypsin II) and exhibited highly duplicated and diversified during evolution. We also observed markedly different domain sequences and organizations between and within the subfamilies, suggesting potential neofunctionalization. Diatoms contain both subfamilies of trypsin, with higher numbers of genes and more environment-responsive expression of trypsin than other lineages. The duplication and subsequent neofunctionalization of the trypsin family may be important in diatoms for adapting to dynamical environmental conditions, contributing to diatoms' dominance in the coastal oceans. This work advances our knowledge on the distribution and neofunctionalization of this ancient enzyme and creates a new window of research on phytoplankton biology.

The PRSS3P2 and TRY7 deletion copy number variant modifies risk for chronic pancreatitis

BACKGROUND

PRSS1 and PRSS2 constitute the only functional copies of a tandemly-arranged five-trypsinogen-gene cluster (i.e., PRSS1, PRSS3P1, PRSS3P2, TRY7 and PRSS2) on chromosome 7q35. Variants in PRSS1 and PRSS2, including missense and copy number variants (CNVs), have been reported to predispose to or protect against chronic pancreatitis (CP). We wondered whether a common trypsinogen pseudogene deletion CNV (that removes two of the three trypsinogen pseudogenes, PRSS3P2 and TRY7) might be associated with CP causation/predisposition.

METHODS

We analyzed the common PRSS3P2 and TRY7 deletion CNV in a total of 1536 CP patients and 3506 controls from France, Germany, India and Japan by means of quantitative fluorescent multiplex polymerase chain reaction.

RESULTS

We demonstrated that the deletion CNV variant was associated with a protective effect against CP in the French, German and Japanese cohorts whilst a trend toward the same association was noted in the Indian cohort. Meta-analysis under a dominant model yielded a pooled odds ratio (OR) of 0.68 (95% confidence interval (CI) 0.52-0.89; p = 0.005) whereas an allele-based meta-analysis yielded a pooled OR of 0.84 (95% CI 0.77-0.92; p = 0.0001). This protective effect is explicable by reference to the recent finding that the still functional PRSS3P2/TRY7 pseudogene enhancers upregulate pancreatic PRSS2 expression.

CONCLUSIONS

The common PRSS3P2 and TRY7 deletion CNV was associated with a reduced risk for CP. This finding provides additional support for the emerging view that dysregulated PRSS2 expression represents a discrete mechanism underlying CP predisposition or protection.

High-quality read-based phasing of cystic fibrosis cohort informs genetic understanding of disease modification

Phasing of heterozygous alleles is critical for interpretation of cis-effects of disease-relevant variation. We sequenced 477 individuals with cystic fibrosis (CF) using linked-read sequencing, which display an average phase block N50 of 4.39 Mb. We use these samples to construct a graph representation of CFTR haplotypes, demonstrating its utility for understanding complex CF alleles. These are visualized in a Web app, CFTbaRcodes, that enables interactive exploration of CFTR haplotypes present in this cohort. We perform fine-mapping and phasing of the chr7q35 trypsinogen locus associated with CF meconium ileus, an intestinal obstruction at birth associated with more severe CF outcomes and pancreatic disease. A 20-kb deletion polymorphism and a PRSS2 missense variant p.Thr8Ile (rs62473563) are shown to independently contribute to meconium ileus risk (p = 0.0028, p = 0.011, respectively) and are PRSS2 pancreas eQTLs (p = 9.5 × 10⁻⁷ and p = 1.4 × 10⁻⁴, respectively), suggesting the mechanism by which these polymorphisms contribute to CF. The phase information from linked reads provides a putative causal explanation for variation at a CF-relevant locus, which also has implications for the genetic basis of non-CF pancreatitis, to which this locus has been reported to contribute.

Attentive Variational Information Bottleneck for TCR-peptide interaction prediction

MOTIVATION

We present a multi-sequence generalization of Variational Information Bottleneck and call the resulting model Attentive Variational Information Bottleneck (AVIB). Our AVIB model leverages multi-head self-attention to implicitly approximate a posterior distribution over latent encodings conditioned on multiple input sequences. We apply AVIB to a fundamental immuno-oncology problem: predicting the interactions between T-cell receptors (TCRs) and peptides.

RESULTS

Experimental results on various datasets show that AVIB significantly outperforms state-of-the-art methods for TCR-peptide interaction prediction. Additionally, we show that the latent posterior distribution learned by AVIB is particularly effective for the unsupervised detection of out-of-distribution amino acid sequences.

AVAILABILITY AND IMPLEMENTATION

The code and the data used for this study are publicly available at: https://github.com/nec-research/vibtcr.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Targeted long-read sequencing facilitates phased diploid assembly and genotyping of the human T cell receptor alpha, delta, and beta loci

T cell receptors (TCRs) recognize peptide fragments presented by the major histocompatibility complex (MHC) and are critical to T cell-mediated immunity. Recent data have indicated that genetic diversity within TCR-encoding gene regions is underexplored, limiting understanding of the impact of TCR loci polymorphisms on TCR function in disease, even though TCR repertoire signatures (1) are heritable and (2) associate with disease phenotypes. To address this, we developed a targeted long-read sequencing approach to generate highly accurate haplotype resolved assemblies of the TCR beta (TRB) and alpha/delta (TRA/D) loci, facilitating the genotyping of all variant types, including structural variants. We validate our approach using two mother-father-child trios and 5 unrelated donors representing multiple populations. This resulted in improved genotyping accuracy and the discovery of 84 undocumented V, D, J, and C alleles, demonstrating the utility of this framework for improving our understanding of TCR diversity and function in disease.

Attention-aware contrastive learning for predicting T cell receptor-antigen binding specificity

MOTIVATION

It has been proven that only a small fraction of the neoantigens presented by major histocompatibility complex (MHC) class I molecules on the cell surface can elicit T cells. This restriction can be attributed to the binding specificity of T cell receptor (TCR) and peptide-MHC complex (pMHC). Computational prediction of T cells binding to neoantigens is a challenging and unresolved task.

RESULTS

In this paper, we proposed an attention-aware contrastive learning model, ATMTCR, to infer the TCR-pMHC binding specificity. For each TCR sequence, we used a transformer encoder to transform it to latent representation, and then masked a percentage of amino acids guided by attention weights to generate its contrastive view. Compared to fully-supervised baseline model, we verified that contrastive learning-based pretraining on large-scale TCR sequences significantly improved the prediction performance of downstream tasks. Interestingly, masking a percentage of amino acids with low attention weights yielded best performance compared to other masking strategies. Comparison experiments on two independent datasets demonstrated our method achieved better performance than other existing algorithms. Moreover, we identified important amino acids and their positional preference through attention weights, which indicated the potential interpretability of our proposed model.

On TCR binding predictors failing to generalize to unseen peptides

Several recent studies investigate TCR-peptide/-pMHC binding prediction using machine learning or deep learning approaches. Many of these methods achieve impressive results on test sets, which include peptide sequences that are also included in the training set. In this work, we investigate how state-of-the-art deep learning models for TCR-peptide/-pMHC binding prediction generalize to unseen peptides. We create a dataset including positive samples from IEDB, VDJdb, McPAS-TCR, and the MIRA set, as well as negative samples from both randomization and 10X Genomics assays. We name this collection of samples TChard. We propose the hard split, a simple heuristic for training/test split, which ensures that test samples exclusively present peptides that do not belong to the training set. We investigate the effect of different training/test splitting techniques on the models’ test performance, as well as the effect of training and testing the models using mismatched negative samples generated randomly, in addition to the negative samples derived from assays. Our results show that modern deep learning methods fail to generalize to unseen peptides. We provide an explanation why this happens and verify our hypothesis on the TChard dataset. We then conclude that robust prediction of TCR recognition is still far for being solved.

Profiling genes encoding the adaptive immune receptor repertoire with gAIRR Suite

Adaptive immune receptor repertoire (AIRR) is encoded by T cell receptor (TR) and immunoglobulin (IG) genes. Profiling these germline genes encoding AIRR (abbreviated as gAIRR) is important in understanding adaptive immune responses but is challenging due to the high genetic complexity. Our gAIRR Suite comprises three modules. gAIRR-seq, a probe capture-based targeted sequencing pipeline, profiles gAIRR from individual DNA samples. gAIRR-call and gAIRR-annotate call alleles from gAIRR-seq reads and annotate whole-genome assemblies, respectively. We gAIRR-seqed TRV and TRJ of seven Genome in a Bottle (GIAB) DNA samples with 100% accuracy and discovered novel alleles. We also gAIRR-seqed and gAIRR-called the TR and IG genes of a subject from both the peripheral blood mononuclear cells (PBMC) and oral mucosal cells. The calling results from these two cell types have a high concordance (99% for all known gAIRR alleles). We gAIRR-annotated 36 genomes to unearth 325 novel TRV alleles and 29 novel TRJ alleles. We could further profile the flanking sequences, including the recombination signal sequence (RSS). We validated two structural variants for HG002 and uncovered substantial differences of gAIRR genes in references GRCh37 and GRCh38. gAIRR Suite serves as a resource to sequence, analyze, and validate germline TR and IG genes to study various immune-related phenotypes.

Tetraploid Ancestry Provided Atlantic Salmon With Two Paralogue Functional T Cell Receptor Beta Regions Whereof One Is Completely Novel

Protective cellular immune responses have been difficult to study in fish, due to lack of basic understanding of their T cell populations, and tools to study them. Cellular immunity is thus mostly ignored in vaccination and infection studies compared to humoral responses. High throughput sequencing, as well as access to well assembled genomes, now advances studies of cellular responses. Here we have used such resources to describe organization of T cell receptor beta genes in Atlantic salmon. Salmonids experienced a unique whole genome duplication approximately 94 million years ago, which provided these species with many functional duplicate genes, where some duplicates have evolved new functions or sub-functions of the original gene copy. This is also the case for T cell receptor beta, where Atlantic salmon has retained two paralogue T cell receptor beta regions on chromosomes 01 and 09. Compared to catfish and zebrafish, the genomic organization in both regions is unique, each chromosomal region organized with dual variable- diversity- joining- constant genes in a head to head orientation. Sequence identity of the chromosomal constant sequences between TRB01 and TRB09 is suggestive of rapid diversification, with only 67 percent as opposed to the average 82-90 percent for other duplicated genes. Using virus challenged samples we find both regions expressing bona fide functional T cell receptor beta molecules. Adding the 292 variable T cell receptor alpha genes to the 100 variable TRB genes from 14 subgroups, Atlantic salmon has one of the most diverse T cell receptor alpha beta repertoire of any vertebrate studied so far. Perhaps salmonid cellular immunity is more advanced than we have imagined.

A Decentralized Kidney Transplant Biopsy Classifier for Transplant Rejection Developed Using Genes of the Banff-Human Organ Transplant Panel

Introduction

A decentralized and multi-platform-compatible molecular diagnostic tool for kidney transplant biopsies could improve the dissemination and exploitation of this technology, increasing its clinical impact. As a first step towards this molecular diagnostic tool, we developed and validated a classifier using the genes of the Banff-Human Organ Transplant (B-HOT) panel extracted from a historical Molecular Microscope^® Diagnostic system microarray dataset. Furthermore, we evaluated the discriminative power of the B-HOT panel in a clinical scenario.

Materials and Methods

Gene expression data from 1,181 kidney transplant biopsies were used as training data for three random forest models to predict kidney transplant biopsy Banff categories, including non-rejection (NR), antibody-mediated rejection (ABMR), and T-cell-mediated rejection (TCMR). Performance was evaluated using nested cross-validation. The three models used different sets of input features: the first model (B-HOT Model) was trained on only the genes included in the B-HOT panel, the second model (Feature Selection Model) was based on sequential forward feature selection from all available genes, and the third model (B-HOT+ Model) was based on the combination of the two models, i.e. B-HOT panel genes plus highly predictive genes from the sequential forward feature selection. After performance assessment on cross-validation, the best-performing model was validated on an external independent dataset based on a different microarray version.

Results

The best performances were achieved by the B-HOT+ Model, a multilabel random forest model trained on B-HOT panel genes with the addition of the 6 most predictive genes of the Feature Selection Model (ST7, KLRC4-KLRK1, TRBC1, TRBV6-5, TRBV19, and ZFX), with a mean accuracy of 92.1% during cross-validation. On the validation set, the same model achieved Area Under the ROC Curve (AUC) of 0.965 and 0.982 for NR and ABMR respectively.

Discussion

This kidney transplant biopsy classifier is one step closer to the development of a decentralized kidney transplant biopsy classifier that is effective on data derived from different gene expression platforms. The B-HOT panel proved to be a reliable highly-predictive panel for kidney transplant rejection classification. Furthermore, we propose to include the aforementioned 6 genes in the B-HOT panel for further optimization of this commercially available panel.

Naive and memory T cells TCR-HLA-binding prediction

T cells recognize antigens through the interaction of their T cell receptor (TCR) with a peptide-major histocompatibility complex (pMHC) molecule. Following thymic-positive selection, TCRs in peripheral naive T cells are expected to bind MHC alleles of the host. Peripheral clonal selection is expected to further increase the frequency of antigen-specific TCRs that bind to the host MHC alleles. To check for a systematic preference for MHC-binding T cells in TCR repertoires, we developed Natural Language Processing-based methods to predict TCR-MHC binding independently of the peptide presented for Class I MHC alleles. We trained a classifier on published TCR-pMHC binding pairs and obtained a high area under curve (AUC) of over 0.90 on the test set. However, when applied to TCR repertoires, the accuracy of the classifier dropped. We thus developed a two-stage prediction model, based on large-scale naive and memory TCR repertoires, denoted TCR HLA-binding predictor (CLAIRE). Since each host carries multiple human leukocyte antigen (HLA) alleles, we first computed whether a TCR on a CD8 T cell binds an MHC from any of the host Class-I HLA alleles. We then performed an iteration, where we predict the binding with the most probable allele from the first round. We show that this classifier is more precise for memory than for naïve cells. Moreover, it can be transferred between datasets. Finally, we developed a CD4-CD8 T cell classifier to apply CLAIRE to unsorted bulk sequencing datasets and showed a high AUC of 0.96 and 0.90 on large datasets. CLAIRE is available through a GitHub at: https://github.com/louzounlab/CLAIRE, and as a server at: https://claire.math.biu.ac.il/Home.

IMGT®Homo sapiens IG and TR Loci, Gene Order, CNV and Haplotypes: New Concepts as a Paradigm for Jawed Vertebrates Genome Assemblies

IMGT®, the international ImMunoGeneTics information system®, created in 1989, by Marie-Paule Lefranc (Université de Montpellier and CNRS), marked the advent of immunoinformatics, a new science which emerged at the interface between immunogenetics and bioinformatics for the study of the adaptive immune responses. IMGT® is based on a standardized nomenclature of the immunoglobulin (IG) and T cell receptor (TR) genes and alleles from fish to humans and on the IMGT unique numbering for the variable (V) and constant (C) domains of the immunoglobulin superfamily (IgSF) of vertebrates and invertebrates, and for the groove (G) domain of the major histocompatibility (MH) and MH superfamily (MhSF) proteins. IMGT® comprises 7 databases, 17 tools and more than 25,000 pages of web resources for sequences, genes and structures, based on the IMGT Scientific chart rules generated from the IMGT-ONTOLOGY axioms and concepts. IMGT® reference directories are used for the analysis of the NGS high-throughput expressed IG and TR repertoires (natural, synthetic and/or bioengineered) and for bridging sequences, two-dimensional (2D) and three-dimensional (3D) structures. This manuscript focuses on the IMGT®Homo sapiens IG and TR loci, gene order, copy number variation (CNV) and haplotypes new concepts, as a paradigm for jawed vertebrates genome assemblies.

T cell receptor beta germline variability is revealed by inference from repertoire data

BACKGROUND

T and B cell receptor (TCR, BCR) repertoires constitute the foundation of adaptive immunity. Adaptive immune receptor repertoire sequencing (AIRR-seq) is a common approach to study immune system dynamics. Understanding the genetic factors influencing the composition and dynamics of these repertoires is of major scientific and clinical importance. The chromosomal loci encoding for the variable regions of TCRs and BCRs are challenging to decipher due to repetitive elements and undocumented structural variants.

METHODS

To confront this challenge, AIRR-seq-based methods have recently been developed for B cells, enabling genotype and haplotype inference and discovery of undocumented alleles. However, this approach relies on complete coverage of the receptors' variable regions, whereas most T cell studies sequence a small fraction of that region. Here, we adapted a B cell pipeline for undocumented alleles, genotype, and haplotype inference for full and partial AIRR-seq TCR data sets. The pipeline also deals with gene assignment ambiguities, which is especially important in the analysis of data sets of partial sequences.

RESULTS

From the full and partial AIRR-seq TCR data sets, we identified 39 undocumented polymorphisms in T cell receptor Beta V (TRBV) and 31 undocumented 5 ' UTR sequences. A subset of these inferences was also observed using independent genomic approaches. We found that a single nucleotide polymorphism differentiating between the two documented T cell receptor Beta D2 (TRBD2) alleles is strongly associated with dramatic changes in the expressed repertoire.

CONCLUSIONS

We reveal a rich picture of germline variability and demonstrate how a single nucleotide polymorphism dramatically affects the composition of the whole repertoire. Our findings provide a basis for annotation of TCR repertoires for future basic and clinical studies.

Identifying COVID-19-Specific Transcriptomic Biomarkers with Machine Learning Methods

COVID-19, a severe respiratory disease caused by a new type of coronavirus SARS-CoV-2, has been spreading all over the world. Patients infected with SARS-CoV-2 may have no pathogenic symptoms, i.e., presymptomatic patients and asymptomatic patients. Both patients could further spread the virus to other susceptible people, thereby making the control of COVID-19 difficult. The two major challenges for COVID-19 diagnosis at present are as follows: (1) patients could share similar symptoms with other respiratory infections, and (2) patients may not have any symptoms but could still spread the virus. Therefore, new biomarkers at different omics levels are required for the large-scale screening and diagnosis of COVID-19. Although some initial analyses could identify a group of candidate gene biomarkers for COVID-19, the previous work still could not identify biomarkers capable for clinical use in COVID-19, which requires disease-specific diagnosis compared with other multiple infectious diseases. As an extension of the previous study, optimized machine learning models were applied in the present study to identify some specific qualitative host biomarkers associated with COVID-19 infection on the basis of a publicly released transcriptomic dataset, which included healthy controls and patients with bacterial infection, influenza, COVID-19, and other kinds of coronavirus. This dataset was first analysed by Boruta, Max-Relevance and Min-Redundancy feature selection methods one by one, resulting in a feature list. This list was fed into the incremental feature selection method, incorporating one of the classification algorithms to extract essential biomarkers and build efficient classifiers and classification rules. The capacity of these findings to distinguish COVID-19 with other similar respiratory infectious diseases at the transcriptomic level was also validated, which may improve the efficacy and accuracy of COVID-19 diagnosis.

Contribution of T Cell Receptor Alpha and Beta CDR3, MHC Typing, V and J Genes to Peptide Binding Prediction

Introduction

Predicting the binding specificity of T Cell Receptors (TCR) to MHC-peptide complexes (pMHCs) is essential for the development of repertoire-based biomarkers. This affinity may be affected by different components of the TCR, the peptide, and the MHC allele. Historically, the main element used in TCR-peptide binding prediction was the Complementarity Determining Region 3 (CDR3) of the beta chain. However, recently the contribution of other components, such as the alpha chain and the other V gene CDRs has been suggested. We use a highly accurate novel deep learning-based TCR-peptide binding predictor to assess the contribution of each component to the binding.

Methods

We have previously developed ERGO-I (pEptide tcR matchinG predictiOn), a sequence-based T-cell receptor (TCR)-peptide binding predictor that employs natural language processing (NLP) -based methods. We improved it to create ERGO-II by adding the CDR3 alpha segment, the MHC typing, V and J genes, and T cell type (CD4+ or CD8+) as to the predictor. We then estimate the contribution of each component to the prediction.

Results and Discussion

ERGO-II provides for the first time high accuracy prediction of TCR-peptide for previously unseen peptides. For most tested peptides and all measures of binding prediction accuracy, the main contribution was from the beta chain CDR3 sequence, followed by the beta chain V and J and the alpha chain, in that order. The MHC allele was the least contributing component. ERGO-II is accessible as a webserver at http://tcr2.cs.biu.ac.il/ and as a standalone code at https://github.com/IdoSpringer/ERGO-II.

The T Cell Receptor (TRB) Locus in Tursiops truncatus: From Sequence to Structure of the Alpha/Beta Heterodimer in the Human/Dolphin Comparison

The bottlenose dolphin (Tursiops truncatus) belongs to the Cetartiodactyla and, similarly to other cetaceans, represents the most successful mammalian colonization of the aquatic environment. Here we report a genomic, evolutionary, and expression study of T. truncatus T cell receptor beta (TRB) genes. Although the organization of the dolphin TRB locus is similar to that of the other artiodactyl species, with three in tandem D-J-C clusters located at its 3' end, its uniqueness is given by the reduction of the total length due essentially to the absence of duplications and to the deletions that have drastically reduced the number of the germline TRBV genes. We have analyzed the relevant mature transcripts from two subjects. The simultaneous availability of rearranged T cell receptor α (TRA) and TRB cDNA from the peripheral blood of one of the two specimens, and the human/dolphin amino acids multi-sequence alignments, allowed us to calculate the most likely interactions at the protein interface between the alpha/beta heterodimer in complex with major histocompatibility class I (MH1) protein. Interacting amino acids located in the complementarity-determining region according to IMGT numbering (CDR-IMGT) of the dolphin variable V-alpha and beta domains were identified. According to comparative modelization, the atom pair contact sites analysis between the human MH1 grove (G) domains and the T cell receptor (TR) V domains confirms conservation of the structure of the dolphin TR/pMH.

Prediction of Specific TCR-Peptide Binding From Large Dictionaries of TCR-Peptide Pairs

Current sequencing methods allow for detailed samples of T cell receptors (TCR) repertoires. To determine from a repertoire whether its host had been exposed to a target, computational tools that predict TCR-epitope binding are required. Currents tools are based on conserved motifs and are applied to peptides with many known binding TCRs. We employ new Natural Language Processing (NLP) based methods to predict whether any TCR and peptide bind. We combined large-scale TCR-peptide dictionaries with deep learning methods to produce ERGO (pEptide tcR matchinG predictiOn), a highly specific and generic TCR-peptide binding predictor. A set of standard tests are defined for the performance of peptide-TCR binding, including the detection of TCRs binding to a given peptide/antigen, choosing among a set of candidate peptides for a given TCR and determining whether any pair of TCR-peptide bind. ERGO reaches similar results to state of the art methods in these tests even when not trained specifically for each test. The software implementation and data sets are available at https://github.com/louzounlab/ERGO. ERGO is also available through a webserver at: http://tcr.cs.biu.ac.il/.

Evolution of the T-Cell Receptor (TR) Loci in the Adaptive Immune Response: The Tale of the TRG Locus in Mammals

T lymphocytes are the principal actors of vertebrates' cell-mediated immunity. Like B cells, they can recognize an unlimited number of foreign molecules through their antigen-specific heterodimer receptors (TRs), which consist of αβ or γδ chains. The diversity of the TRs is mainly due to the unique organization of the genes encoding the α, β, γ, and δ chains. For each chain, multi-gene families are arranged in a TR locus, and their expression is guaranteed by the somatic recombination process. A great plasticity of the gene organization within the TR loci exists among species. Marked structural differences affect the TR γ (TRG) locus. The recent sequencing of multiple whole genome provides an opportunity to examine the TR gene repertoire in a systematic and consistent fashion. In this review, we report the most recent findings on the genomic organization of TRG loci in mammalian species in order to show differences and similarities. The comparison revealed remarkable diversification of both the genomic organization and gene repertoire across species, but also unexpected evolutionary conservation, which highlights the important role of the T cells in the immune response.

The effects of CTX damage or inhibition of bone marrow hematopoiesis and GM-CSF stimulation of bone marrow hematopoiesis on the peripheral blood TCRβ CDR3 repertoire of BALB/c mice

Objective: This paper aims to investigate the dynamic changes of the T-cell receptor (TCR) β complementarity-determining region 3 (CDR3) repertoire during cyclophosphamide or Cytoxan (CTX) damage or inhibition of bone marrow hematopoiesis caused by a reduction of peripheral blood white blood cells (WBCs) in BALB/c mice.Methods: We analyze TCR CDR3 repertoire of BALB/c mice including (1) NS control group (2) CTX damage group (3) CTX damage + GM-CSF recovery group (4) CTX damage + auto-recovery group.Results: The number of WBCs in the CTX group is significantly lower than that in the NS group and after GM-CSF injection, the GM-CSF group is higher than that in the NS group. The diversity of the CTX damage group is the highest and there is a significant difference in high-frequency clonal proliferation between the CTX damage group and CTX damage + GM-CSF recovery group compared with the NS control group. In addition, the numbers of unique productive CDR3 overlapping numbers in the four experimental groups are similar.Conclusions: These data reveal that CTX significantly reduced the number of WBCs and ratio of high-frequency TCR CDR3 sequences, and indirectly increased the diversity of the TCR CDR3 repertoire. GM-CSF quickly restored the number of WBCs, and partially restored changes in the TCR CDR3 repertoire induced by CTX. Results from monitoring the dynamic changes of the TCR CDR3 repertoire can be used to assess the effects of CTX and GM-CSF on the function of peripheral blood T cells and to explore the possible underlying mechanisms.

Genomic organization of the chicken TCRβ locus originated by duplication of a Vβ segment combined with a trypsinogen gene

Based on the latest assembly of the red jungle fowl (Gallus gallus) genome sequence, we characterized the detailed genomic organization of the T cell receptor beta (TCRβ) locus of chicken. The chicken TCRβ locus spans approximately 210 kb, and is organized in a typical translocon organization as previously reported. Within this locus, a total of 16 germline Vβ gene segments were classified into three subgroups, containing 11, four, and one members, respectively. Phylogenetic analysis revealed that the chicken Vβ3.1 segment was homologous with the duck Vβ1 subgroup, and further clustered with Vβ  segments from reptiles but not amphibians. We also identified nine protease serine 1 (PRSS1) and three protease serine 2 (PRSS2) genes, which were interspersed within the chicken TCRβ locus. Dot-plot analysis of the chicken TCRβ locus against itself revealed that the 5' part of the locus had arisen through a series of tandem duplication events. The homology units were composed of one Vβ1 segment followed by a PRSS1 gene, or one Vβ2 segment followed by a PRSS2 gene. This duplication pattern, in which the Vβ segments and trypsinogen genes form a duplication unit, was unique to TCRβ loci of chicken and duck, but not observed in TCRβ loci of other tetrapods studied thus far. By analyzing the cloned TCRβ cDNA sequences, we found that the usage pattern of Vβ segments was consistent with the results of previous studies. These studies showed that members of the Vβ1 subgroup are preferentially utilized in V-D-J recombination. Furthermore, we found that the Vβ3.1 segment participated into V-D-J recombination, but at a very low frequency. The length distribution of the chicken complementarity-determining region 3β (CDR3β) showed a tendency similar to that observed for the duck.

Worldwide genetic variation of the IGHV and TRBV immune receptor gene families in humans

This article presents a comprehensive study of the IGHV and TRBV gene families in a globally diverse sample of humans and shows that the two gene families exhibit starkly different patterns of variation.

The immunoglobulin heavy variable (IGHV) and T cell beta variable (TRBV) loci are among the most complex and variable regions in the human genome. Generated through a process of gene duplication/deletion and diversification, these loci can vary extensively between individuals in copy number and contain genes that are highly similar, making their analysis technically challenging. Here, we present a comprehensive study of the functional gene segments in the IGHV and TRBV loci, quantifying their copy number and single-nucleotide variation in a globally diverse sample of 109 (IGHV) and 286 (TRBV) humans from over a 100 populations. We find that the IGHV and TRBV gene families exhibit starkly different patterns of variation. In addition to providing insight into the different evolutionary paths of the IGHV and TRBV loci, our results are also important to the adaptive immune repertoire sequencing community, where the lack of frequencies of common alleles and copy number variants is hampering existing analytical pipelines.

The Florida manatee (Trichechus manatus latirostris) T cell receptor loci exhibit V subgroup synteny and chain-specific evolution

The Florida manatee (Trichechus manatus latirostris) has limited diversity in the immunoglobulin heavy chain. We therefore investigated the antigen receptor loci of the other arm of the adaptive immune system: the T cell receptor. Manatees are the first species from Afrotheria, a basal eutherian superorder, to have an in-depth characterization of all T cell receptor loci. By annotating the genome and expressed transcripts, we found that each chain has distinct features that correlates to their individual functions. The genomic organization also plays a role in modulating sequence conservation between species. There were extensive V subgroup synteny blocks in the TRA and TRB loci between T. m. latirostris and human. Increased genomic locus complexity correlated to increased locus synteny. We also identified evidence for a VHD pseudogene for the first time in a eutherian mammal. These findings emphasize the value of including species within this basal eutherian radiation in comparative studies.

Analysis of TCRβ and TCRγ genes in Chinese alligator provides insights into the evolution of TCR genes in jawed vertebrates

All jawed vertebrates have four T cell receptor (TCR) chains that are expressed by thymus-derived lymphocytes and play a major role in animal immune defence. However, few studies have investigated the TCR chains of crocodilians compared with those of birds and mammals, despite their key evolutionary position linking amphibians, reptiles, birds and mammals. Here, employing an Alligator sinensis genomic bacterial artificial chromosome (BAC) library and available genome data, we characterized the genomic organization, evolution and expression of TRB and TRG loci in Alligator sinensis. According to the sequencing data, the Alligator sinensis TRB locus spans approximately 500 Kb of genomic DNA containing two D-J-C clusters and 43 V gene segments and is organized as Vβ₍₃₉₎-pJβ1-pCβ1-pDβ1-Dβ2- Jβ2₍₁₂₎-Cβ2-Vβ₍₄₎, whereas the TRG locus spans 115 Kb of DNA genomic sequence consisting of 18 V gene segments, nine J gene segments and one C gene segment and is organized in a classical translocon pattern as Vγ₍₁₈₎-Jγ₍₉₎-Cγ. Moreover, syntenic analysis of TRB and TRG chain loci suggested a high degree of conserved synteny in the genomic regions across mammals, birds and Alligator sinensis. By analysing the cloned TRB/TRG cDNA, we identified the usage pattern of V families in the expressed TRB and TRG. An analysis of the junctions of the recombined VJ revealed the presence of N and P nucleotides in both expressed TRB and TRG sequences. Phylogenetic analysis revealed that TRB and TRG loci possess distinct evolutionary patterns. Most Alligator sinensis V subgroups have closely related orthologues in chicken and duck, and a small number of Alligator sinensis V subgroups have orthologues in mammals, which supports the hypothesis that crocodiles are the closest relatives of birds and mammals. Collectively, these data provide insights into TCR gene evolution in vertebrates and improve our understanding of the Alligator sinensis immune system.

The Bayesian optimist's guide to adaptive immune receptor repertoire analysis

Probabilistic modeling is fundamental to the statistical analysis of complex data. In addition to forming a coherent description of the data-generating process, probabilistic models enable parameter inference about given datasets. This procedure is well developed in the Bayesian perspective, in which one infers probability distributions describing to what extent various possible parameters agree with the data. In this paper, we motivate and review probabilistic modeling for adaptive immune receptor repertoire data then describe progress and prospects for future work, from germline haplotyping to adaptive immune system deployment across tissues. The relevant quantities in immune sequence analysis include not only continuous parameters such as gene use frequency but also discrete objects such as B-cell clusters and lineages. Throughout this review, we unravel the many opportunities for probabilistic modeling in adaptive immune receptor analysis, including settings for which the Bayesian approach holds substantial promise (especially if one is optimistic about new computational methods). From our perspective, the greatest prospects for progress in probabilistic modeling for repertoires concern ancestral sequence estimation for B-cell receptor lineages, including uncertainty from germline genotype, rearrangement, and lineage development.

Ci8 short, a novel LPS-induced peptide from the ascidian Ciona intestinalis, modulates responses of the human immune system

The selective modulation of immunity is an emerging concept driven by the vast advances in our understanding of this crucial host defense system. Invertebrates have raised researchers' interest as potential sources of new bioactive molecules owing to their antibacterial, anticancer and immunomodulatory activities. A LipoPolySaccharide (LPS) challenge in the ascidian Ciona intestinalis generates the transcript, Ci8 short, with cis-regulatory elements in the 3' UTR region that are essential for shaping innate immune responses. The derived amino acidic sequence in silico analysis showed specific binding to human Major Histocompatibility Complex (MHC) Class I and Class II alleles. The role of Ci8 short peptide was investigated in a more evolved immune system using human Peripheral Blood Mononuclear Cells (PBMCs) as in vitro model. The biological activities of this molecule include the activation of 70kDa TCR ζ chain Associated Protein kinase (ZAP-70) and T Cell Receptor (TCR) Vβ oligo clonal selection on CD4⁺ T lymphocytes as well as increased proliferation and IFN-γ secretion. Furthermore Ci8 short affects CD4⁺/CD25^high induced regulatory T cells (iTreg) subset selection which co-expressed the functional markers TGF-β1/Latency Associated Protein (LAP) and CD39/CD73. This paper describes a new molecule that modulates important responses of the human adaptive immune system.

Analysis of the CDR3 length repertoire and the diversity of T cell receptor α and β chains in swine CD4+ and CD8+ T lymphocytes

The T cell receptor (TCR) is a complex heterodimer that recognizes fragments of antigens as peptides and binds to major histocompatibility complex molecules. The TCR α and β chains possess three hypervariable regions termed complementarity determining regions (CDR1, 2 and 3). CDR3 is responsible for recognizing processed antigen peptides. Immunoscope spectratyping is a simple technique for analyzing CDR3 polymorphisms and sequence length diversity, in order to investigate T cell function and the pattern of TCR utilization. The present study employed this technique to analyze CDR3 polymorphisms and the sequence length diversity of TCR α and β chains in porcine CD4⁺ and CD8⁺ T cells. Polymerase chain reaction products of 19 TCR α variable regions (AV) and 20 TCR β variable regions (BV) gene families obtained from the CD4⁺ and CD8⁺ T cells revealed a clear band following separation by 1.5% agarose gel electrophoresis, and each family exhibited >8 bands following separation by 6% sequencing gel electrophoresis. CDR3 spectratyping of all identified TCR AV and BV gene families in the sorted CD4⁺ and CD8⁺ T cells by GeneScan, demonstrated a standard Gaussian distribution with >8 peaks. CDR3 in CD4⁺ and CD8⁺ T cells demonstrated different expression patterns. The majority of CDR3 recombined in frame and the results revealed that there were 10 and 14 amino acid discrepancies between the longest and shortest CDR3 lengths in specific TCR AV and TCR BV gene families, respectively. The results demonstrated that CDR3 polymorphism and length diversity demonstrated different expression and utilization patterns in CD4⁺ and CD8⁺ T cells. These results may facilitate future research investigating the porcine TCR CDR3 gene repertoire as well as the functional complexity and specificity of the TCR molecule.

Alcohol Consumption and Cigarette Smoking are Important Modifiers of the Association Between Acute Pancreatitis and the PRSS1-PRSS2 Locus in Men

OBJECTIVES

The present study was designed to investigate whether the susceptibility to acute pancreatitis (AP) attributable to polymorphism rs10273639 at the PRSS1-PRSS2 locus is dependent on alcohol consumption and cigarette smoking.

METHODS

A total of 603 unrelated Russian individuals including 304 patients with physician-diagnosed AP and 299 sex- and age-matched healthy controls have been recruited for the study. A polymorphism rs10273639 (-408C>T) of PRSS1-PRSS2 was genotyped by TaqMan-based assay.

RESULTS

A variant allele -408T (P = 0.003) and genotypes -408CT plus TT (P = 0.002) were associated with decreased AP risk only in men. The odds ratios for AP in the CC homozygotes versus the variant genotypes were 1.95 [95% confidence interval (CI), 0.65-5.85; P = 0.23], 1.72 (95% CI, 0.93-3.20; P = 0.08), and 2.37 (95% CI, 1.09-5.13; P = 0.03) for men who consumed up to 28, 29 to 59, and more than 60 alcohol drinks a week, respectively. Cigarette smokers with the -408CC genotype had an increased risk of AP (odds ratio, 2.07; 95% CI, 1.25-3.42; P = 0.004), whereas nonsmoker carriers did not have a disease risk (odds ratio, 1.48; 95% CI, 0.58-3.82; P = 0.42).

CONCLUSIONS

We confirmed a robust association of polymorphism rs10273639 at PRSS1-PRSS2 with AP in the Russian population. The present study is the first to show that relationship between the locus and disease is significantly modified by alcohol consumption and cigarette smoking.

A comprehensive analysis of the germline and expressed TCR repertoire in White Peking duck

Recently, many immune-related genes have been extensively studied in ducks, but relatively little is known about their TCR genes. Here, we determined the germline and expressed repertoire of TCR genes in White Peking duck. The genomic organization of the duck TCRα/δ, TCRγ and unconventional TCRδ2 loci are highly conserved with their counterparts in mammals or chickens. By contrast, the duck TCRβ locus is organized in an unusual pattern, (Vβ)_n-Dβ-(Jβ)₂-Cβ1-(Jβ)₄-Cβ2, which differs from the tandem-aligned clusters in mammals or the translocon organization in some teleosts. Excluding the first exon encoding the immunoglobulin domain, the subsequent exons of the two Cβ show significant diversity in nucleotide sequence and exon structure. Based on the nucleotide sequence identity, 49 Vα, 30 Vδ, 13 Vβ and 15 Vγ unique gene segments are classified into 3 Vα, 5 Vδ, 4 Vβ and 6 Vγ subgroups, respectively. Phylogenetic analyses revealed that most duck V subgroups, excluding Vβ1, Vγ5 and Vγ6, have closely related orthologues in chicken. The coding joints of all cDNA clones demonstrate conserved mechanisms that are used to increase junctional diversity. Collectively, these data provide insight into the evolution of TCRs in vertebrates and improve our understanding of the avian immune system.

Copy number variations exploration of multiple genes in Graves' disease

BACKGROUND

Few previous published papers reported copy number variations of genes could affect the predisposition of Graves' disease (GD). Herein, the aim of this study was to explore the association between copy number variations (CNV) profile and GD.

METHODS

The preliminary copy number microarray used to screen copy number variant genes was performed in 6 GD patients. Five CNV candidate genes (CFH, CFHR1, KIAA0125, UGT2B15, and UGT2B17) were then validated in an independent set of samples (50 GD patients and 50 matched healthy ones) by the Accucopy assay method. The CNV of the other 2 genes TRY6 and CCL3L1 was investigated in 144 GD patients and 144 healthy volunteers by the definitive genotyping technique using the Taqman quantitative polymerase-chain-reaction (Taqman qPCR). TRY6 gene-associated single nucleotide polymorphism (SNP), rs13230029, was genotyped by the PCR-ligase detection reaction (LDR) in 675 GD patients and 898 healthy controls.

RESULTS

There were no correlation of the gene copy number (GCN) of CFH, CFHR1, KIAA0125, UGT2B15, and UGT2B17 with GD. In comparison with that of controls, the GCN distribution of TRY6 and CCL3L1 in GD patients did not show significantly differ (P > 0.05). Furthermore, TRY6-related polymorphism (rs13230029) showed no difference between GD patients and controls. No correlation was found between CNV or SNP genotype and clinical phenotypes. Generally, there were no link of the copy numbers of several genes, including CFH, CFHR1, KIAA0125, UGT2B15, UGT2B17, TRY6, and CCL3L1 to GD.

CONCLUSION

Our results clearly indicated that the copy number variations of multiple genes, namely CFH, CFHR1, KIAA0125, UGT2B15, UGT2B17, TRY6, and CCL3L1, were not associated with the development of GD.

TRIg: a robust alignment pipeline for non-regular T-cell receptor and immunoglobulin sequences

BACKGROUND

T cells and B cells are essential in the adaptive immunity via expressing T cell receptors and immunoglogulins respectively for recognizing antigens. To recognize a wide variety of antigens, a highly diverse repertoire of receptors is generated via complex recombination of the receptor genes. Reasonably, frequencies of the recombination events have been shown to predict immune diseases and provide insights into the development of immunity. The field is further boosted by high-throughput sequencing and several computational tools have been released to analyze the recombined sequences. However, all current tools assume regular recombination of the receptor genes, which is not always valid in data prepared using a RACE approach. Compared to the traditional multiplex PCR approach, RACE is free of primer bias, therefore can provide accurate estimation of recombination frequencies. To handle the non-regular recombination events, a new computational program is needed.

RESULTS

We propose TRIg to handle non-regular T cell receptor and immunoglobulin sequences. Unlike all current programs, TRIg does alignments to the whole receptor gene instead of only to the coding regions. This brings new computational challenges, e.g., ambiguous alignments due to multiple hits to repetitive regions. To reduce ambiguity, TRIg applies a heuristic strategy and incorporates gene annotation to identify authentic alignments. On our own and public RACE datasets, TRIg correctly identified non-regularly recombined sequences, which could not be achieved by current programs. TRIg also works well for regularly recombined sequences.

CONCLUSIONS

TRIg takes into account non-regular recombination of T cell receptor and immunoglobulin genes, therefore is suitable for analyzing RACE data. Such analysis will provide accurate estimation of recombination events, which will benefit various immune studies directly. In addition, TRIg is suitable for studying aberrant recombination in immune diseases. TRIg is freely available at https://github.com/TLlab/trig .

Quantification of a Selective Expansion of T Cell Receptor Vβ by Superantigen Using Real-Time PCR

Selective expansion of T cells bearing specific T cell receptor Vβ segments is a hallmark of superantigens. Analyzing Vβ specificity of superantigens is important for characterizing newly discovered superantigens and understanding differential T cell responses to each toxin. Here, we describe a real-time PCR method using SYBR green I and primers specific to Cβ and Vβ genes for an absolute quantification. The established method was applied to quantify a selective expansion of T cell receptor Vβ expansion by superantigens and generated accurate, reproducible, and comparable results.

Internal Duplications of DH, JH, and C Region Genes Create an Unusual IgH Gene Locus in Cattle

It has been suspected for many years that cattle possess two functional IgH gene loci, located on Bos taurus autosome (BTA) 21 and BTA11, respectively. In this study, based on fluorescence in situ hybridization and additional experiments, we showed that all functional bovine IgH genes were located on BTA21, and only a truncated μCH2 exon was present on BTA11. By sequencing of seven bacterial artificial chromosome clones screened from a Hostein cow bacterial artificial chromosome library, we generated a 678-kb continuous genomic sequence covering the bovine IGHV, IGHD, IGHJ, and IGHC genes, which are organized as IGHVn-IGHDn-IGHJn-IGHM1-(IGHDP-IGHV3-IGHDn)3-IGHJn-IGHM2-IGHD-IGHG3-IGHG1-IGHG2-IGHE-IGHA. Although both of two functional IGHM genes, IGHM1 and IGHM2, can be expressed via independent VDJ recombinations, the IGHM2 can also be expressed through class switch recombination. Likely because more IGHD segments can be involved in the expression of IGHM2, the IGHM2 gene was shown to be dominantly expressed in most tissues throughout different developmental stages. Based on the length and identity of the coding sequence, the 23 IGHD segments identified in the locus could be divided into nine subgroups (termed IGHD1 to IGHD9). Except two members of IGHD9 (14 nt in size), all other functional IGHD segments are longer than 30 nt, with the IGHD8 gene (149 bp) to be the longest. These remarkably long germline IGHD segments play a pivotal role in generating the exceptionally great H chain CDR 3 length variability in cattle.

High-content molecular profiling of T-cell therapy in oncology

Recent clinical data have revealed the remarkable potential for T-cell-modulating agents to induce potent and durable responses in a subset of cancer patients. In this review, we discuss molecular approaches, platforms, and strategies that enable a broader interrogation of the activity of agents that modulate the activity of tumor-specific T cells, to more comprehensively understand how and why the agents succeed and fail, as well as examples of data sets generated in clinical trials that have provided important insights into the biological activity of T-cell therapies and that support further rational development of this exciting treatment modality.

Cross-Reactivity of TCR Repertoire: Current Concepts, Challenges, and Implication for Allotransplantation

Being able to track donor reactive T cells during the course of organ transplantation is a key to improve the graft survival, to prevent graft dysfunction, and to adapt the immunosuppressive regimen. The attempts of transplant immunologists have been for long hampered by the large size of the alloreactive T cell repertoire. Understanding how self-TCR can interact with allogeneic MHC is a key to critically appraise the different assays available to analyze the TCR Vβ repertoire usage. In this report, we will review conceptually and experimentally the process of cross-reactivity. We will then highlight what can be learned from allotransplantation, a situation of artificial cross-reactivity. Finally, the low- and high-resolution techniques to characterize the TCR Vβ repertoire usage in transplantation will be critically discussed.

Gene conversion between cationic trypsinogen (PRSS1) and the pseudogene trypsinogen 6 (PRSS3P2) in patients with chronic pancreatitis

Mutations of the human cationic trypsinogen gene (PRSS1) are frequently found in association with hereditary pancreatitis. The most frequent variants p.N29I and p.R122H are recognized as disease-causing mutations. Three pseudogene paralogs in the human trypsinogen family, including trypsinogen 6 (PRSS3P2), carry sequence variations in exon 3 that mimic the p.R122H mutation. In routine genetic testing of patients with chronic pancreatitis, we identified in two unrelated individuals similar gene conversion events of 24-71 nucleotides length between exon 3 of the PRSS1 (acceptor) and PRSS3P2 (donor) genes. The converted allele resulted in three nonsynonymous alterations c.343T>A (p.S115T), c.347G>C (p.R116P), and c.365_366delinsAT (p.R122H). Functional analysis of the conversion triple mutant revealed markedly increased autoactivation resulting in high and sustained trypsin activity in the presence of chymotrypsin C. This activation phenotype was identical to that of the p.R122H mutant. In addition, cellular secretion of the triple mutant from transfected HEK 293T cells was increased about twofold and this effect was attributable to mutation p.R116P. Our observations confirm and extend the notion that recombination events between members of the trypsinogen family can generate high-risk PRSS1 alleles. The pathogenic phenotype of the novel conversion is explained by a unique combination of increased trypsinogen activation and secretion.

Genomic analysis offers insights into the evolution of the bovine TRA/TRD locus

Background

The TRA/TRD locus contains the genes for V(D)J somatic rearrangement of TRA and TRD chains expressed by αβ and γδ T cells respectively. Previous studies have demonstrated that the bovine TRA/TRD locus contains an exceptionally large number of TRAV/TRDV genes. In this study we combine genomic and transcript analysis to provide insights into the evolutionary development of the bovine TRA/TRD locus and the remarkable TRAV/TRDV gene repertoire.

Results

Annotation of the UMD3.1 assembly identified 371 TRAV/TRDV genes (distributed in 42 subgroups), 3 TRDJ, 6 TRDD, 62 TRAJ and single TRAC and TRDC genes, most of which were located within a 3.5 Mb region of chromosome 10. Most of the TRAV/TRDV subgroups have multiple members and several have undergone dramatic expansion, most notably TRDV1 (60 genes). Wide variation in the proportion of pseudogenes within individual subgroups, suggest that differential ‘birth’ and ‘death’ rates have been used to form a functional bovine TRAV/TRDV repertoire which is phylogenetically distinct from that of humans and mice. The expansion of the bovine TRAV/TRDV gene repertoire has predominantly been achieved through a complex series of homology unit (regions of DNA containing multiple gene) replications. Frequent co-localisation within homology units of genes from subgroups with low and high pseudogene proportions suggest that replication of homology units driven by evolutionary selection for the former may have led to a ‘collateral’ expansion of the latter. Transcript analysis was used to define the TRAV/TRDV subgroups available for recombination of TRA and TRD chains and demonstrated preferential usage of different subgroups by the expressed TRA and TRD repertoires, indicating that TRA and TRD selection have had distinct impacts on the evolution of the TRAV/TRDV repertoire.

Conclusion

Both TRA and TRD selection have contributed to the evolution of the bovine TRAV/TRDV repertoire. However, our data suggest that due to homology unit duplication TRD selection for TRDV1 subgroup expansion may have substantially contributed to the genomic expansion of several TRAV subgroups. Such data demonstrate how integration of genomic and transcript data can provide a more nuanced appreciation of the evolutionary dynamics that have led to the dramatically expanded bovine TRAV/TRDV repertoire.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-994) contains supplementary material, which is available to authorized users.

Genetics of acute and chronic pancreatitis: An update

Progress made in identifying the genetic susceptibility underlying acute and chronic pancreatitis has benefitted the clinicians in understanding the pathogenesis of the disease in a better way. The identification of mutations in cationic trypsinogen gene (PRSS1 gene; functional gain mutations) and serine protease inhibitor kazal type 1 (SPINK1 gene; functional loss mutations) and other potential susceptibility factors in genes that play an important role in the pancreatic secretory functions or response to inflammation during pancreatic injury has changed the current concepts and understanding of a complex multifactorial disease like pancreatitis. An individual’s susceptibility to the disease is governed by genetic factors in combination with environmental factors. Candidate gene and genetic linkage studies have identified polymorphisms in cationic trypsinogen (PRSS1), SPINK1, cystic fibrosis trans-membrane conductance regulator (CFTR), Chymotrypsinogen C (CTRC), Cathepsin B (CTSB) and calcium sensing receptor (CASR). Individuals with polymorphisms in the mentioned genes and other as yet identified genes are at an enhanced risk for the disease. Recently, polymorphisms in genes other than those involved in “intra-pancreatic trypsin regulatory mechanism” namely Claudin-2 (CLDN2) and Carboxypeptidase A1 (CPA1) gene have also been identified for their association with pancreatitis. With ever growing number of studies trying to identify the genetic susceptibility in the form of single nucleotide polymorphisms, this review is an attempt to compile the available information on the topic.

Human cationic trypsinogen (PRSS1) variants and chronic pancreatitis

Variations in the serine protease 1 (PRSS1) gene encoding human cationic trypsinogen have been conclusively associated with autosomal dominant hereditary pancreatitis and sporadic nonalcoholic chronic pancreatitis. Most high-penetrance PRSS1 variants increase intrapancreatic trypsin activity by stimulating trypsinogen autoactivation and/or by inhibiting chymotrypsin C-dependent trypsinogen degradation. Alternatively, some PRSS1 variants can cause trypsinogen misfolding, which results in intracellular retention and degradation with consequent endoplasmic reticulum stress. However, not all PRSS1 variants are pathogenic, and clinical relevance of rare variants is often difficult to ascertain. Here we review the PRSS1 variants published since 1996 and discuss their functional properties and role in chronic pancreatitis.