VDJML: a file format with tools for capturing the results of inferring immune receptor rearrangements

Analysis of CDR3 Sequences from T-Cell Receptor β in Acute Respiratory Distress Syndrome

Acute Respiratory Distress Syndrome (ARDS) is an illness that typically develops in people who are significantly ill or have serious injuries. ARDS is characterized by fluid build-up that occurs in the alveoli. T-cells are implicated as playing a role in the modulation of the aberrant response leading to excessive tissue damage and, eventually, ARDS. Complementarity Determining Region 3 (CDR3) sequences derived from T-cells are key players in the adaptive immune response. This response is governed by an elaborate specificity for distinct molecules and the ability to recognize and vigorously respond to repeated exposures to the same molecules. Most of the diversity in T-cell receptors (TCRs) is contained in the CDR3 regions of the heterodimeric cell-surface receptors. For this study, we employed the novel technology of immune sequencing to assess lung edema fluid. Our goal was to explore the landscape of CDR3 clonal sequences found within these samples. We obtained more than 3615 CDR3 sequences across samples in the study. Our data demonstrate that: (1) CDR3 sequences from lung edema fluid exhibit distinct clonal populations, and (2) CDR3 sequences can be further characterized based on biochemical features. Analysis of these CDR3 sequences offers insight into the CDR3-driven T-cell repertoire of ARDS. These findings represent the first step towards applications of this technology with these types of biological samples in the context of ARDS.

Immune Repertoire Analysis on High-Performance Computing Using VDJServer V1: A Method by the AIRR Community

AIRR-seq data sets are usually large and require specialized analysis methods and software tools. A typical Illumina MiSeq sequencing run generates 20-30 million 2 × 300 bp paired-end sequence reads, which roughly corresponds to 15 GB of sequence data to be processed. Other platforms like NextSeq, which is useful in projects where the full V gene is not needed, create about 400 million 2 × 150 bp paired-end reads. Because of the size of the data sets, the analysis can be computationally expensive, particularly the early analysis steps like preprocessing and gene annotation that process the majority of the sequence data. A standard desktop PC may take 3-5 days of constant processing for a single MiSeq run, so dedicated high-performance computational resources may be required.VDJServer provides free access to high-performance computing (HPC) at the Texas Advanced Computing Center (TACC) through a graphical user interface (Christley et al. Front Immunol 9:976, 2018). VDJServer is a cloud-based analysis portal for immune repertoire sequence data that provides access to a suite of tools for a complete analysis workflow, including modules for preprocessing and quality control of sequence reads, V(D)J gene assignment, repertoire characterization, and repertoire comparison. Furthermore, VDJServer has parallelized execution for tools such as IgBLAST, so more compute resources are utilized as the size of the input data grows. Analysis that takes days on a desktop PC might take only a few hours on VDJServer. VDJServer is a free, publicly available, and open-source licensed resource. Here, we describe the workflow for performing immune repertoire analysis on VDJServer's high-performance computing.

Single-Cell Sequencing Reveals the Transcriptome and TCR Characteristics of pTregs and in vitro Expanded iTregs

Regulatory T cells (Tregs) play a critical role in the maintenance of immune tolerance and tumor evasion. However, the relative low proportion of these cells in peripheral blood and tissues has hindered many studies. We sought to establish a rapamycin-based in vitro Treg expansion procedure in patients diagnosed with colorectal cancer and perform single-cell sequencing to explore the characteristics of Treg cells. CD25+ cells enriched from peripheral blood mononuclear cells (PBMC) of colorectal tumor patients were cultured in X-VIVO15 medium, supplemented with 5% human AB serum, L-glutamine, rapamycin, interleukin-2 (IL-2), and Dynabeads human Treg expander for 21 days to expand Tregs. Treg cells with satisfactory phenotype and function were successfully expanded from CD4+CD25+ cells in patients with colorectal cancer. The median expansion fold was 75 (range, 20-105-fold), and >90.0% of the harvest cells were CD4+CD25+CD127^dim/- cells. The ratio of CD4+CD25+Foxp3+ cells exceeded 60%. Functional assays showed that iTregs significantly inhibited CD8+T cell proliferation in vitro. Single-cell sequencing showed that the transcriptome of pTreg (CD4+CD25+CD127^dim/- cells isolated from PBMC of colorectal cancer patients) and iTreg (CD4+CD25+CD127^dim/- cells expanded in vitro according to the above regimen) cells were interlaced. pTregs exhibited enhanced suppressive function, whereas iTregs exhibited increased proliferative capacity. TCR repertoire analysis indicated minimal overlap between pTregs and iTregs. Pseudo-time trajectory analysis of Tregs revealed that pTregs were a continuum composed of three main branches: activated/effector, resting and proliferative Tregs. In contrast, in vitro expanded iTregs were a mixture of proliferating and activated/effector cells. The expression of trafficking receptors was also different in pTregs and iTregs. Various chemokine receptors were upregulated in pTregs. Activated effector pTregs overexpressed the chemokine receptor CCR10, which was not expressed in iTregs. The chemokine CCL28 was overexpressed in colorectal cancer and associated with poor prognosis. CCR10 interacted with CCL28 to mediate the recruitment of Treg into tumors and accelerated tumor progression. Depletion of CCR10+Treg cells from tumor microenvironment (TME) could be used as an effective treatment strategy for colorectal cancer patients. Our data distinguished the transcriptomic characteristics of different subsets of Treg cells and revealed the context-dependent functions of different populations of Treg cells, which was crucial to the development of alternative therapeutic strategies for Treg cells in autoimmune disease and cancer.

AIRR Community Standardized Representations for Annotated Immune Repertoires

Increased interest in the immune system's involvement in pathophysiological phenomena coupled with decreased DNA sequencing costs have led to an explosion of antibody and T cell receptor sequencing data collectively termed "adaptive immune receptor repertoire sequencing" (AIRR-seq or Rep-Seq). The AIRR Community has been actively working to standardize protocols, metadata, formats, APIs, and other guidelines to promote open and reproducible studies of the immune repertoire. In this paper, we describe the work of the AIRR Community's Data Representation Working Group to develop standardized data representations for storing and sharing annotated antibody and T cell receptor data. Our file format emphasizes ease-of-use, accessibility, scalability to large data sets, and a commitment to open and transparent science. It is composed of a tab-delimited format with a specific schema. Several popular repertoire analysis tools and data repositories already utilize this AIRR-seq data format. We hope that others will follow suit in the interest of promoting interoperable standards.

VDJServer: A Cloud-Based Analysis Portal and Data Commons for Immune Repertoire Sequences and Rearrangements

Background

Recent technological advances in immune repertoire sequencing have created tremendous potential for advancing our understanding of adaptive immune response dynamics in various states of health and disease. Immune repertoire sequencing produces large, highly complex data sets, however, which require specialized methods and software tools for their effective analysis and interpretation.

Results

VDJServer is a cloud-based analysis portal for immune repertoire sequence data that provide access to a suite of tools for a complete analysis workflow, including modules for preprocessing and quality control of sequence reads, V(D)J gene segment assignment, repertoire characterization, and repertoire comparison. VDJServer also provides sophisticated visualizations for exploratory analysis. It is accessible through a standard web browser via a graphical user interface designed for use by immunologists, clinicians, and bioinformatics researchers. VDJServer provides a data commons for public sharing of repertoire sequencing data, as well as private sharing of data between users. We describe the main functionality and architecture of VDJServer and demonstrate its capabilities with use cases from cancer immunology and autoimmunity.

Conclusion

VDJServer provides a complete analysis suite for human and mouse T-cell and B-cell receptor repertoire sequencing data. The combination of its user-friendly interface and high-performance computing allows large immune repertoire sequencing projects to be analyzed with no programming or software installation required. VDJServer is a web-accessible cloud platform that provides access through a graphical user interface to a data management infrastructure, a collection of analysis tools covering all steps in an analysis, and an infrastructure for sharing data along with workflows, results, and computational provenance. VDJServer is a free, publicly available, and open-source licensed resource.

Reproducibility and Reuse of Adaptive Immune Receptor Repertoire Data

High-throughput sequencing (HTS) of immunoglobulin (B-cell receptor, antibody) and T-cell receptor repertoires has increased dramatically since the technique was introduced in 2009 (1–3). This experimental approach explores the maturation of the adaptive immune system and its response to antigens, pathogens, and disease conditions in exquisite detail. It holds significant promise for diagnostic and therapy-guiding applications. New technology often spreads rapidly, sometimes more rapidly than the understanding of how to make the products of that technology reliable, reproducible, or usable by others. As complex technologies have developed, scientific communities have come together to adopt common standards, protocols, and policies for generating and sharing data sets, such as the MIAME protocols developed for microarray experiments. The Adaptive Immune Receptor Repertoire (AIRR) Community formed in 2015 to address similar issues for HTS data of immune repertoires. The purpose of this perspective is to provide an overview of the AIRR Community’s founding principles and present the progress that the AIRR Community has made in developing standards of practice and data sharing protocols. Finally, and most important, we invite all interested parties to join this effort to facilitate sharing and use of these powerful data sets (join@airr-community.org).

T Cell Receptor Profiling in Type 1 Diabetes

PURPOSE OF REVIEW

The genetic susceptibility and dominant protection for type 1 diabetes (T1D) associated with human leukocyte antigen (HLA) haplotypes, along with minor risk variants, have long been thought to shape the T cell receptor (TCR) repertoire and eventual phenotype of autoreactive T cells that mediate β-cell destruction. While autoantibodies provide robust markers of disease progression, early studies tracking autoreactive T cells largely failed to achieve clinical utility.

RECENT FINDINGS

Advances in acquisition of pancreata and islets from T1D organ donors have facilitated studies of T cells isolated from the target tissues. Immunosequencing of TCR α/β-chain complementarity determining regions, along with transcriptional profiling, offers the potential to transform biomarker discovery. Herein, we review recent studies characterizing the autoreactive TCR signature in T1D, emerging technologies, and the challenges and opportunities associated with tracking TCR molecular profiles during the natural history of T1D.

The Repertoire Dissimilarity Index as a method to compare lymphocyte receptor repertoires

BACKGROUND

The B and T cells of the human adaptive immune system leverage a highly diverse repertoire of antigen-specific receptors to protect the human body from pathogens. The sequencing and analysis of immune repertoires is emerging as an important tool to understand immune responses, whether beneficial or harmful (in the case of autoimmunity). However, methods for studying these repertoires, and for directly comparing different immune repertoires, are lacking.

RESULTS

In this paper, we present a non-parametric method for directly comparing sequencing repertoires, with the goal of rigorously quantifying differences in V, D, and J gene segment utilization. This method, referred to as the Repertoire Dissimilarity Index (RDI), uses a bootstrapped subsampling approach to account for variance in sequencing depth, and, coupled with a data simulation approach, allows for direct quantification of the average variation between repertoires. We use the RDI method to recapitulate known differences in the formation of the CD4+ and CD8+ T cell repertoires, and further show that antigen-driven activation of naïve CD8+ T cells is more selective than in the CD4+ repertoire, resulting in a more specialized CD8+ memory repertoire.

CONCLUSIONS

We prove that the RDI method is an accurate and versatile method for comparisons of immune repertoires. The RDI method has been implemented as an R package, and is available for download through Bitbucket.

Reproducibility and Reuse of Adaptive Immune Receptor Repertoire Data

High-throughput sequencing (HTS) of immunoglobulin (B-cell receptor, antibody) and T-cell receptor repertoires has increased dramatically since the technique was introduced in 2009 (1-3). This experimental approach explores the maturation of the adaptive immune system and its response to antigens, pathogens, and disease conditions in exquisite detail. It holds significant promise for diagnostic and therapy-guiding applications. New technology often spreads rapidly, sometimes more rapidly than the understanding of how to make the products of that technology reliable, reproducible, or usable by others. As complex technologies have developed, scientific communities have come together to adopt common standards, protocols, and policies for generating and sharing data sets, such as the MIAME protocols developed for microarray experiments. The Adaptive Immune Receptor Repertoire (AIRR) Community formed in 2015 to address similar issues for HTS data of immune repertoires. The purpose of this perspective is to provide an overview of the AIRR Community's founding principles and present the progress that the AIRR Community has made in developing standards of practice and data sharing protocols. Finally, and most important, we invite all interested parties to join this effort to facilitate sharing and use of these powerful data sets (join@airr-community.org).