VDJ Gene Usage in IgM Repertoires of Rhesus and Cynomolgus Macaques

Macaques are frequently used to evaluate candidate vaccines and to study infection-induced antibody responses, requiring an improved understanding of their naïve immunoglobulin (IG) repertoires. Baseline gene usage frequencies contextualize studies of antigen-specific immune responses, providing information about how easily one may stimulate a response with a particular VDJ recombination. Studies of human IgM repertoires have shown that IG VDJ gene frequencies vary several orders of magnitude between the most and least utilized genes in a manner that is consistent across many individuals but to date similar analyses are lacking for macaque IgM repertoires. Here, we quantified VDJ gene usage levels in unmutated IgM repertoires of 45 macaques, belonging to two species and four commonly used subgroups: Indian and Chinese origin rhesus macaques and Indonesian and Mauritian origin cynomolgus macaques. We show that VDJ gene frequencies differed greatly between the most and least used genes, with similar overall patterns observed in macaque subgroups and individuals. However, there were also clear differences affecting the use of specific V, D and J genes. Furthermore, in contrast to humans, macaques of both species utilized IGHV4 family genes to a much higher extent and showed evidence of evolutionary expansion of genes of this family. Finally, we used the results to inform the analysis of a broadly neutralizing HIV-1 antibody elicited in SHIV-infected rhesus macaques, RHA1.V2.01, which binds the apex of the Env trimer in a manner that mimics the binding mode of PGT145. We discuss the likelihood that similar antibodies could be elicited in different macaque subgroups.

A Pan-Cancer Analysis of Tumor-Infiltrating B Cell Repertoires

Tumor-infiltrating B cells can play an important role in anti-tumor responses but their presence is not well understood. In this study, we extracted the B cell receptor repertoires from 9522 tumor and adjacent non-tumor samples across 28 tumor types in the Cancer Genome Atlas project and performed diversity and network analysis. We identified differences in diversity and network statistics across tumor types and subtypes and observed a trend towards increased clonality in primary tumors compared to adjacent non-tumor tissues. We also found significant associations between the repertoire features and mutation load, tumor stage, and age. Our V-gene usage analysis identified similar V-gene usage patterns in colorectal and endometrial cancers. Lastly, we evaluated the prognostic value of the repertoire features and identified significant associations with survival in seven tumor types. This study warrants further research into better understanding the role of tumor-infiltrating B cells across a wide range of tumor types.

Recent Advances in the Genetic of MALT Lymphomas

Mucosa-associated lymphoid tissue (MALT) lymphomas are a diverse group of lymphoid neoplasms with B-cell origin, occurring in adult patients and usually having an indolent clinical behavior. These lymphomas may arise in different anatomic locations, sharing many clinicopathological characteristics, but also having substantial variances in the aetiology and genetic alterations. Chromosomal translocations are recurrent in MALT lymphomas with different prevalence among different sites, being the 4 most common: t(11;18)(q21;q21), t(1;14)(p22;q32), t(14;18)(q32;q21), and t(3;14)(p14.1;q32). Several chromosomal numerical abnormalities have also been described, but probably represent secondary genetic events. The mutational landscape of MALT lymphomas is wide, and the most frequent mutations are: TNFAIP3, CREBBP, KMT2C, TET2, SPEN, KMT2D, LRP1B, PRDM1, EP300, TNFRSF14, NOTCH1/NOTCH2, and B2M, but many other genes may be involved. Similar to chromosomal translocations, certain mutations are enriched in specific lymphoma types. In the same line, variation in immunoglobulin gene usage is recognized among MALT lymphoma of different anatomic locations. In the last decade, several studies have analyzed the role of microRNA, transcriptomics and epigenetic alterations, further improving our knowledge about the pathogenic mechanisms in MALT lymphoma development. All these advances open the possibility of targeted directed treatment and push forward the concept of precision medicine in MALT lymphomas.

Cell type-specific analysis by single-cell profiling identifies a stable mammalian tRNA-mRNA interface and increased translation efficiency in neurons

The correlation between codon and anticodon pools influences the efficiency of translation, but whether differences exist in these pools across individual cells is unknown. We determined that codon usage and amino acid demand are highly stable across different cell types using available mouse and human single-cell RNA sequencing atlases. After demonstrating the robustness of ATAC-seq measurements for the analysis of tRNA gene usage, we quantified anticodon usage and amino acid supply in both mouse and human single-cell ATAC-seq atlases. We found that tRNA gene usage is overall coordinated across cell types, except in neurons, which clustered separately from other cell types. Integration of these datasets revealed a strong and statistically significant correlation between amino acid supply and demand across almost all cell types. Neurons have an enhanced translation efficiency over other cell types, driven by an increased supply of tRNAAla (AGC) anticodons. This results in faster decoding of the Ala-GCC codon, as determined by cell type-specific ribosome profiling, suggesting that the reduction of tRNAAla (AGC) anticodon pools may be implicated in neurological pathologies. This study, the first such examination of codon usage, anticodon usage, and translation efficiency resolved at the cell type level with single-cell information, identifies a conserved landscape of translation elongation across mammalian cellular diversity and evolution.

Biased TCR gene usage in citrullinated Tenascin C specific T-cells in rheumatoid arthritis

We aimed to search for common features in the autoreactive T cell receptor (TCR) repertoire in patients with rheumatoid arthritis (RA), focusing on the newly identified candidate antigen citrullinated Tenascin C (cit-TNC). Mononuclear cells from peripheral blood or synovial fluid of eight RA-patients positive for the RA-associated HLA-DRB1*04:01 allele were in-vitro cultured with recently identified citrullinated peptides from Tenascin C. Antigen-specific T cells were isolated using peptide-HLA tetramer staining and subsequently single-cell sequenced for paired alpha/beta TCR analyses by bioinformatic tools. TCRs were re-expressed for further studies of antigen-specificity and T cell responses. Autoreactive T cell lines could be grown out from both peripheral blood and synovial fluid. We demonstrate the feasibility of retrieving true autoreactive TCR sequences by validating antigen-specificity in T cell lines with re-expressed TCRs. One of the Tenascin C peptides, cit-TNC22, gave the most robust T cell responses including biased TCR gene usage patterns. The shared TCR-beta chain signature among the cit-TNC22-specific TCRs was evident in blood and synovial fluid of different patients. The identification of common elements in the autoreactive TCR repertoire gives promise to the possibility of both immune monitoring of the autoimmune components in RA and of future antigen- or TCR-targeted specific intervention in subsets of patients.

Single-cell RNA sequencing of Plasmodium vivax sporozoites reveals stage- and species-specific transcriptomic signatures

Plasmodium vivax sporozoites reside in the salivary glands of a mosquito before infecting a human host. Previous transcriptome-wide studies in populations of these forms were limited in their ability to elucidate cell-to-cell variation, thereby masking cellular states potentially important in understanding transmission outcomes. In this study, we performed transcription profiling on 9,947 P. vivax sporozoites to assess the extent to which they differ at single-cell resolution. We show that sporozoites residing in the mosquito's salivary glands exist in distinct developmental states, as defined by their transcriptomic signatures. Additionally, relative to P. falciparum, P. vivax displays overlapping and unique gene usage patterns, highlighting conserved and species-specific gene programs. Notably, distinguishing P. vivax from P. falciparum were a subset of P. vivax sporozoites expressing genes associated with translational regulation and repression. Finally, our comparison of single-cell transcriptomic data from P. vivax sporozoite and erythrocytic forms reveals gene usage patterns unique to sporozoites. In defining the transcriptomic signatures of individual P. vivax sporozoites, our work provides new insights into the factors driving their developmental trajectory and lays the groundwork for a more comprehensive P. vivax cell atlas.

Biased TCR gene usage in citrullinated Tenascin C specific T-cells in rheumatoid arthritis

Objectives To search for common features in the autoreactive T cell receptor (TCR)-repertoire in patients with rheumatoid arthritis (RA), focusing on the newly identified candidate antigen citrullinated Tenascin C (cit-TNC). Methods Mononuclear cells from peripheral blood or synovial fluid of eight RA-patients positive for the RA-associated HLA-DRB1*04:01 allele were in-vitro cultured with recently identified citrullinated peptides from Tenascin C. Antigen-specific T cells were isolated using peptide-HLA tetramer staining and subsequently single-cell sequenced for paired alpha/beta TCR analyses by bioinformatic tools. TCRs were re-expressed for further studies of antigen-specificity and T cell responses. Results Autoreactive T cell lines could be grown out from both peripheral blood and synovial fluid. We demonstrate the feasibility of retrieving true autoreactive TCR sequences by validating antigen-specificity in T cell lines with re-expressed TCRs. One of the Tenascin C peptides, cit-TNC22, gave the most robust T cell responses including biased TCR gene usage patterns. The shared TCR-beta chain signature among the citTNC22-specific TCRs was evident in blood and synovial fluid of different patients. Conclusion The identification of common elements in the autoreactive TCR repertoire gives promise to the possibility of both immune monitoring of the autoimmune components in RA and of future antigen- or TCR-targeted specific intervention in subsets of patients.

Fully Automated Workflows Quantify and Report Key T-Cell and B-Cell Receptor Biomarkers Relevant to Immuno-Oncology and Heme-Oncology Research

Background T-cell and B-cell repertoire analysis is used in oncology research, to understand the etiology of complex disease phenotypes, for the identification of biomarkers predictive of disease burden, outcome, and response to treatment, and for research in diagnosis and recurrence monitoring. Key predictors include secondary and tertiary repertoire features not reported by existing sequencing software solutions. For example, due to ongoing somatic hypermutation in mature B-cell receptors, the underlying sequence of a given clone can accumulate base differences and appear as several distinct clones with smaller frequencies, thereby hampering the ability of analysis software to detect its presence as a single dominant clone with the highest frequency. This has particularly detrimental implications for research in disorders such as follicular lymphoma and may require clonal lineage analysis for proper mitigation. Therefore, to aid the downstream analytics of biomarker identification and the study of complex disease, we developed fully automated analysis solutions that directly compute and report several key features (clonal lineage, amongst several others described below) pertinent to this area of research. Results We developed the Oncomine™ TCR Beta-SR, TCR Gamma-SR, BCR IGH-SR and BCR IGKL-SR workflows on Ion Reporter™ to characterize T-cell (β, γ chains) and B-cell (heavy and light (κ, δ) chains) repertoires. These workflows generate output tables and visualizations for primary repertoire features such as detected clones (viz., unique rearrangements in the receptor DNA sequence), their frequencies, as well as their somatic hypermutation levels in the case of B-cells (Figure 1a & 1b) for clonality assessment and rare clone detection. The software also quantifies and reports several secondary and tertiary repertoire features in a sample, such as clonal diversity, evenness of the clonal population, and B-cell lineage groupings useful in identifying related sub-clones. It includes spectratyping format plots to simultaneously assess the above features as a function of v-gene usage and CDR3 length combinations (Figure 1c & 1d), thereby providing users a complete snapshot of the repertoire, and also the capability to quickly determine CDR3 lengths and V-gene usage of highly expanded or mutated clones. A separate CDR3 lengths histogram is included, as well as a heatmap that depicts the distributions/intensity of Variable-Joining gene combinations (Figure 1e & 1f). Furthermore, the TCR workflows also report (i) convergence frequencies (fraction of clones with different nucleotide sequences, but identical amino acid sequences), and (ii) haplotype grouping for an analyzed sample, based on V-gene allele genotyping and clustering (Figure 1g). In addition, the long read Oncomine™ BCR IGH-LR workflow uniquely reports the isotype class for every detected clone, and includes a visualization of total reads, clones and lineages in the sample represented by isotype (Figure 1h). Conclusion The Oncomine™ immune repertoire workflows for T-cell and B-cell receptor sequencing were designed to be of high utility in distinct areas of malignancy research, and we expect them to greatly simplify complex downstream analyses. The unique capabilities of the workflows to automatically report secondary and tertiary repertoire features such as (i) clonal lineages for improved dominant clone detection in blood cancers, (ii) TCR clone convergence for prediction of response to immune checkpoint inhibitors [1,2], (iii) TCR haplotype grouping for evaluation of risk factors for autoimmunity and immune-related adverse events [3], and (iv) isotype classification in BCRs for studying pan-cancer immune evasion mechanisms, demonstrate the clear advantages of using these automated workflows over other existing solutions. For research use only. References 1) Looney TJ et al. (2020) TCR Convergence in Individuals Treated With Immune Checkpoint Inhibition for Cancer. Front. Immunol. 10:2985. 2) Naidus et al. (2021) Early changes in the circulating T cells are associated with clinical outcomes after PD-L1 blockade by durvalumab in advanced NSCLC patients. Cancer Immunology, Immunotherapy 70:2095-2102 3) Looney TJ et al. (2019) Haplotype Analysis of the T-Cell Receptor Beta (TCRB) Locus by Long-amplicon TCRB Repertoire Sequencing. Journal of Immunotherapy and Precision Oncology. 2 (4): 137-143. Figure 1 Figure 1. Disclosures Sarda: Thermo Fisher Scientific: Current Employment. Lowman: Thermo Fisher Scientific: Current Employment. Toro: Thermo Fisher Scientific: Current Employment. Pickle: Thermo Fisher Scientific: Current Employment. Looney: Thermo Fisher Scientific: Ended employment in the past 24 months; Singular Genomics: Current Employment. Hyland: Thermo Fisher Scientific: Current Employment.

Novel Allele Detection Tool Benchmark and Application With Antibody Repertoire Sequencing Dataset

Detailed knowledge of the diverse immunoglobulin germline genes is critical for the study of humoral immunity. Hundreds of alleles have been discovered by analyzing antibody repertoire sequencing (Rep-seq or Ig-seq) data via multiple novel allele detection tools (NADTs). However, the performance of these NADTs through antibody sequences with intrinsic somatic hypermutations (SHMs) is unclear. Here, we developed a tool to simulate repertoires by integrating the full spectrum features of an antibody repertoire such as germline gene usage, junctional modification, position-specific SHM and clonal expansion based on 2152 high-quality datasets. We then systematically evaluated these NADTs using both simulated and genuine Ig-seq datasets. Finally, we applied these NADTs to 687 Ig-seq datasets and identified 43 novel allele candidates (NACs) using defined criteria. Twenty-five alleles were validated through findings of other sources. In addition to the NACs detected, our simulation tool, the results of our comparison, and the streamline of this process may benefit further humoral immunity studies via Ig-seq.