Sequencing the peripheral blood B and T cell repertoire - Quantifying robustness and limitations

Tissue-resident memory T cells in human kidney transplants have alloreactive potential

The extent to which tissue-resident memory T (TRM) cells in transplanted organs possess alloreactivity is uncertain. This study investigates the alloreactive potential of TRM cells in kidney explants from 4 patients who experienced severe acute rejection leading to graft loss. Alloreactive T cell receptor (TCR) clones were identified in pretransplant blood samples through mixed lymphocyte reactions, followed by single-cell RNA and TCR sequencing of the proliferated recipient T cells. Subsequently, these TCR clones were traced in the TRM cells of kidney explants, which were also subjected to single-cell RNA and TCR sequencing. The proportion of recipient-derived TRM cells expressing an alloreactive TCR in the 4 kidney explants varied from 0% to 9%. Notably, these alloreactive TCRs were predominantly found among CD4+ and CD8+ TRM cells with an effector phenotype. Intriguingly, these clones were present not only in recipient-derived TRM cells but also in donor-derived TRM cells, constituting up to 4% of the donor population, suggesting the presence of self-reactive TRM cells. Overall, our study demonstrates that T cells with alloreactive potential present in the peripheral blood prior to transplantation can infiltrate the kidney transplant and adopt a TRM phenotype.

Fighting rare cancers: lessons from fibrolamellar hepatocellular carcinoma

The fight against rare cancers faces myriad challenges, including missed or wrong diagnoses, lack of information and diagnostic tools, too few samples and too little funding. Yet many advances in cancer biology, such as the realization that there are tumour suppressor genes, have come from studying well-defined, albeit rare, cancers. Fibrolamellar hepatocellular carcinoma (FLC), a typically lethal liver cancer, mainly affects adolescents and young adults. FLC is both rare, 1 in 5 million, and problematic to diagnose. From the paucity of data, it was not known whether FLC was one cancer or a collection with similar phenotypes, or whether it was genetically inherited or the result of a somatic mutation. A personal journey through a decade of work reveals answers to these questions and a road map of steps and missteps in our fight against a rare cancer.

The use of non-functional clonotypes as a natural calibrator for quantitative bias correction in adaptive immune receptor repertoire profiling

High-throughput sequencing of adaptive immune receptor repertoires is a valuable tool for receiving insights in adaptive immunity studies. Several powerful TCR/BCR repertoire reconstruction and analysis methods have been developed in the past decade. However, detecting and correcting the discrepancy between real and experimentally observed lymphocyte clone frequencies are still challenging. Here, we discovered a hallmark anomaly in the ratio between read count and clone count-based frequencies of non-functional clonotypes in multiplex PCR-based immune repertoires. Calculating this anomaly, we formulated a quantitative measure of V- and J-genes frequency bias driven by multiplex PCR during library preparation called Over Amplification Rate (OAR). Based on the OAR concept, we developed an original software for multiplex PCR-specific bias evaluation and correction named iROAR: immune Repertoire Over Amplification Removal (https://github.com/smiranast/iROAR). The iROAR algorithm was successfully tested on previously published TCR repertoires obtained using both 5' RACE (Rapid Amplification of cDNA Ends)-based and multiplex PCR-based approaches and compared with a biological spike-in-based method for PCR bias evaluation. The developed approach can increase the accuracy and consistency of repertoires reconstructed by different methods making them more applicable for comparative analysis.

The Clonal Diversity of Peripheral B Cell Receptor Immune Repertoire Impaired by Residual Malignant B Cells Predicts Treatment Efficacy in B Cell Lymphoma Patients

Germinal center (GC) is the vital locus for the evolution of naïve B cells into memory B and plasma cells, but also a hotbed for the proliferation of malignant B cells. We hypothesized that malignant B cells may locally or globally impact GCs to produce peripheral B cell receptor immune repertoire (BCR IR) with reduced clonal diversity. In this study, we first validated our hypothesis in a novel human in-vitro GC (hiGC) model. The addition of the diffuse large B cell lymphoma (DLBCL) cells to the hiGC culture attenuated the rate of diversity growth. For clinical validation, we collected samples from 17 DLBCL patients at various points during high-dose therapy and autologous stem cell rescue. The elimination and reestablishment of the patients’ lymphatic pool allowed us to unambiguously monitor the impact of tumor cells on the replenishment of the peripheral BCR IR. Compared to the nine patients who did not relapse after treatment, relapsed patients tended to have a slower rate of recovery regarding the clonal diversity of their peripheral BCR IR. Our results suggest a mechanistic and clinical connection between residual tumor cells and abnormal peripheral BCR IR, which may corelate with treatment efficacy in B cell lymphomas.

Optimization of peripheral blood volume for in silico reconstitution of the human B cell receptor repertoire

B cells recognize antigens via membrane‐expressed B‐cell receptors (BCR) and antibodies. Similar human BCR sequences are frequently found at a significantly higher frequency than that theoretically calculated. Patients infected with SARS‐CoV2 and HIV or with autoimmune diseases share very similar BCRs. Therefore, in silico reconstitution of BCR repertoires and identification of stereotypical BCR sequences related to human pathology have diagnostic potential. Furthermore, monitoring changes of clinically significant BCR sequences and isotype conversion has prognostic potential. For BCR repertoire analysis, peripheral blood (PB) is the most convenient source. However, the optimal human PB volume for in silico reconstitution of the BCR repertoire has not been studied in detail. Here, we sampled 5, 10, and 20 mL PB from the left arm and 40 mL PB from the right arm of two volunteers, reconstituted in silico PB BCR repertoires, and compared their composition. In both volunteers, PB sampling over 20 mL resulted in slight increases in functional unique sequences (FUSs) or almost no increase in repertoire diversity. All FUSs with a frequency above 0.08% or 0.03% in the 40 mL PB BCR repertoire were detected even in the 5 mL PB BCR repertoire from each volunteer. FUSs with a higher frequency were more likely to be found in BCR repertoires from reduced PB volume, and those coexisting in two repertoires showed a statistically significant correlation in frequency irrespective of sampled anatomical site. The correlation was more significant in higher‐frequency FUSs. These observations support the potential of BCR repertoire analysis for diagnosis.

ImmunoDataAnalyzer: a bioinformatics pipeline for processing barcoded and UMI tagged immunological NGS data

Background

Next-generation sequencing (NGS) is nowadays the most used high-throughput technology for DNA sequencing. Among others NGS enables the in-depth analysis of immune repertoires. Research in the field of T cell receptor (TCR) and immunoglobulin (IG) repertoires aids in understanding immunological diseases. A main objective is the analysis of the V(D)J recombination defining the structure and specificity of the immune repertoire. Accurate processing, evaluation and visualization of immune repertoire NGS data is important for better understanding immune responses and immunological behavior.

Results

ImmunoDataAnalyzer (IMDA) is a pipeline we have developed for automatizing the analysis of immunological NGS data. IMDA unites the functionality from carefully selected immune repertoire analysis software tools and covers the whole spectrum from initial quality control up to the comparison of multiple immune repertoires. It provides methods for automated pre-processing of barcoded and UMI tagged immune repertoire NGS data, facilitates the assembly of clonotypes and calculates key figures for describing the immune repertoire. These include commonly used clonality and diversity measures, as well as indicators for V(D)J gene segment usage and between sample similarity. IMDA reports all relevant information in a compact summary containing visualizations, calculations, and sample details, all of which serve for a more detailed overview. IMDA further generates an output file including key figures for all samples, designed to serve as input for machine learning frameworks to find models for differentiating between specific traits of samples.

Conclusions

IMDA constructs TCR and IG repertoire data from raw NGS reads and facilitates descriptive data analysis and comparison of immune repertoires. The IMDA workflow focus on quality control and ease of use for non-computer scientists. The provided output directly facilitates the interpretation of input data and includes information about clonality, diversity, clonotype overlap as well as similarity, and V(D)J gene segment usage. IMDA further supports the detection of sample swaps and cross-sample contamination that potentially occurred during sample preparation. In summary, IMDA reduces the effort usually required for immune repertoire data analysis by providing an automated workflow for processing raw NGS data into immune repertoires and subsequent analysis. The implementation is open-source and available on https://bioinformatics.fh-hagenberg.at/immunoanalyzer/.

Next Generation Imaging Techniques to Define Immune Topographies in Solid Tumors

In recent years, cancer immunotherapy experienced remarkable developments and it is nowadays considered a promising therapeutic frontier against many types of cancer, especially hematological malignancies. However, in most types of solid tumors, immunotherapy efficacy is modest, partly because of the limited accessibility of lymphocytes to the tumor core. This immune exclusion is mediated by a variety of physical, functional and dynamic barriers, which play a role in shaping the immune infiltrate in the tumor microenvironment. At present there is no unified and integrated understanding about the role played by different postulated models of immune exclusion in human solid tumors. Systematically mapping immune landscapes or "topographies" in cancers of different histology is of pivotal importance to characterize spatial and temporal distribution of lymphocytes in the tumor microenvironment, providing insights into mechanisms of immune exclusion. Spatially mapping immune cells also provides quantitative information, which could be informative in clinical settings, for example for the discovery of new biomarkers that could guide the design of patient-specific immunotherapies. In this review, we aim to summarize current standard and next generation approaches to define Cancer Immune Topographies based on published studies and propose future perspectives.

Topology and expressed repertoire of the Felis catus T cell receptor loci

Background

The domestic cat (Felis catus) is an important companion animal and is used as a large animal model for human disease. However, the comprehensive study of adaptive immunity in this species is hampered by the lack of data on lymphocyte antigen receptor genes and usage. The objectives of this study were to annotate the feline T cell receptor (TR) loci and to characterize the expressed repertoire in lymphoid organs of normal cats using high-throughput sequencing.

Results

The Felis catus TRG locus contains 30 genes: 12 TRGV, 12 TRGJ and 6 TRGC, the TRB locus contains 48 genes: 33 TRBV, 2 TRBD, 11 TRBJ, 2 TRBC, the TRD locus contains 19 genes: 11 TRDV, 2 TRDD, 5 TRDJ, 1 TRDC, and the TRA locus contains 127 genes: 62 TRAV, 64 TRAJ, 1 TRAC. Functional feline V genes form monophyletic clades with their orthologs, and clustering of multimember subgroups frequently occurs in V genes located at the 5′ end of TR loci. Recombination signal (RS) sequences of the heptamer and nonamer of functional V and J genes are highly conserved. Analysis of the TRG expressed repertoire showed preferential intra-cassette over inter-cassette rearrangements and dominant usage of the TRGV2–1 and TRGJ1–2 genes. The usage of TRBV genes showed minor bias but TRBJ genes of the second J-C-cluster were more commonly rearranged than TRBJ genes of the first cluster. The TRA/TRD V genes almost exclusively rearranged to J genes within their locus. The TRAV/TRAJ gene usage was relatively balanced while the TRD repertoire was dominated by TRDJ3.

Conclusions

This is the first description of all TR loci in the cat. The genomic organization of feline TR loci was similar to that of previously described jawed vertebrates (gnathostomata) and is compatible with the birth-and-death model of evolution. The large-scale characterization of feline TR genes provides comprehensive baseline data on immune repertoires in healthy cats and will facilitate the development of improved reagents for the diagnosis of lymphoproliferative diseases in cats. In addition, these data might benefit studies using cats as a large animal model for human disease.

Comprehensive RNA-Sequencing Analysis in Peripheral Blood Cells Reveals Differential Expression Signatures with Biomarker Potential for Idiopathic Membranous Nephropathy

To date, the clinical course of idiopathic membranous nephropathy (iMN) remains unclear and lacks direct and effective diagnostic methods. To better understand the host gene expression changes involved in the iMN process and identify the potential signatures for clinical diagnosis, we performed a whole genome-wide transcriptome profile of peripheral blood cells (PBC) from patients with iMN and healthy controls (HCs). A total of 188 differentially expressed genes (DEGs) were detected in patients with iMN versus HCs. Gene ontology (GO) functional enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that these DEGs were mainly correlated with protein targeting, ion homeostasis GO terms, and ribosome and phagosome pathways. The top 10 differentially expressed protein-coding genes with >2-fold changes and high expression levels were validated using quantitative real-time PCR, and showed high consistency with the high-throughput sequencing results. HLA-C, S100A8, and FTH1 genes were selected for further validation and showed the most significant difference between the iMN and HC group, indicating that they could be used as potential clinical diagnostic biomarkers. Our results provide novel potential diagnostic signatures for iMN and have important implications for better understanding the pathogenesis of iMN.