A comprehensive profiling of T- and B-lymphocyte receptor repertoires from a Chinese-origin rhesus macaque by high-throughput sequencing

scTCR-seq and HTS reveal a special novel TRBD2-TRBJ1 rearrangement in mammalian TRB CDR3 repertoire

Mammalian T cell receptor (TCR) beta-chain (TRB) V-D-J rearrangement mainly follows the "12/23 rule", and the "D-J rearrangement preceding the V-(D-J) rearrangement". Owing to the physical position of the D-J-C cluster in the TRB locus, the TRBD2 (D2) gene cannot directly perform inversional rearrangement or deletional/loop-out rearrangement with the TRBJ1 (J1) gene. Our previous studies revealed a single reverse TRBV30 (TRBV31 in mice) gene in the mammalian TRB locus, which can cause indirect rearrangement of the D2 gene and J1 gene; however, the mechanism and proportion involved in germline gene rearrangement are unknown. We obtained TRB CDR3 repertoires of thymus and peripheral tissues from humans and mice by HTS and scTCR-seq and found that 14% of the rearrangements in which the D2 gene is involved are D2-J1 rearrangements (D2-J2 rearrangements account for approximately 86%). The mechanism is that the reverse V30 gene preferentially performs inversional rearrangement with the D2 gene (V30-D2), leading to V30-D2-J1 rearrangement in humans, or the reverse V30 gene preferentially performs inversional rearrangement with the D1 gene (V30-D1), allowing the forward V genes (Vx) to perform Vx-D2-J1 rearrangement. We further found that D2-J1 rearrangements were present in more than 24% and more than 15% of the D2 gene rearrangements in rhesus monkeys and bats, respectively. Moreover, in bovine containing D1J1C1, D3J3C3, and D2J2C2 clusters, more than 11% D3-J1 and D2-J1 rearrangements and more than 22% D2-J3 rearrangements were found. This study provides a new perspective and feasible solution for further research on the significance of the special V-D-J recombination pattern in the mammalian TRB locus and the CDR3 repertoire formed by D2-J1 rearrangement.

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.

Immunization with germ line-targeting SOSIP trimers elicits broadly neutralizing antibody precursors in infant macaques

Adolescents are a growing population of people living with HIV. The period between weaning and sexual debut presents a low-risk window for HIV acquisition, making early childhood an ideal time for implementing an immunization regimen. Because the elicitation of broadly neutralizing antibodies (bnAbs) is critical for an effective HIV vaccine, our goal was to assess the ability of a bnAb B cell lineage-designed HIV envelope SOSIP (protein stabilized by a disulfide bond between gp120-gp41-named "SOS"-and an isoleucine-to-proline point mutation-named "IP"-at residue 559) to induce precursor CD4 binding site (CD4bs)-targeting bnAbs in early life. Infant rhesus macaques received either a BG505 SOSIP, based on the infant BG505 transmitted/founder virus, or the CD4bs germ line-targeting BG505 SOSIP GT1.1 (n = 5 per group). Although both strategies induced durable, high-magnitude plasma autologous virus neutralization responses, only GT1.1-immunized infants (n = 3 of 5) exhibited VRC01-like CD4bs bnAb precursor development. Thus, a multidose immunization regimen with bnAb lineage-designed SOSIPs shows promise for inducing early B cell responses with the potential to mature into protective HIV bnAbs before sexual debut.

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.

Germline-targeting SOSIP trimer immunization elicits precursor CD4 binding-site targeting broadly neutralizing antibodies in infant macaques

A vaccine that can achieve protective immunity prior to sexual debut is critical to prevent the estimated 410,000 new HIV infections that occur yearly in adolescents. As children living with HIV can make broadly neutralizing antibody (bnAb) responses in plasma at a faster rate than adults, early childhood is an opportune window for implementation of a multi-dose HIV immunization strategy to elicit protective immunity prior to adolescence. Therefore, the goal of our study was to assess the ability of a B cell lineage-designed HIV envelope SOSIP to induce bnAbs in early life. Infant rhesus macaques (RMs) received either BG505 SOSIP or the germline-targeting BG505 GT1.1 SOSIP (n=5/group) with the 3M-052-SE adjuvant at 0, 6, and 12 weeks of age. All infant RMs were then boosted with the BG505 SOSIP at weeks 26, 52 and 78, mimicking a pediatric immunization schedule of multiple vaccine boosts within the first two years of life. Both immunization strategies induced durable, high magnitude binding antibodies and plasma autologous virus neutralization that primarily targeted the CD4-binding site (CD4bs) or C3/465 epitope. Notably, three BG505 GT1.1-immunized infants exhibited a plasma HIV neutralization signature reflective of VRC01-like CD4bs bnAb precursor development and heterologous virus neutralization. Finally, infant RMs developed precursor bnAb responses at a similar frequency to that of adult RMs receiving a similar immunization strategy. Thus, a multi-dose immunization regimen with bnAb lineage designed SOSIPs is a promising strategy for inducing protective HIV bnAb responses in childhood prior to adolescence when sexual HIV exposure risk begins.

New insights into the germline genes and CDR3 repertoire of the TCRβ chain in Chiroptera

Introduction

Bats are recognized as natural reservoirs for many viruses, and their unique immune system enables them to coexist with these viruses without frequently exhibiting disease symptoms. However, the current understanding of the bat adaptive immune system is limited due to the lack of a database or tool capable of processing T-cell receptor (TCR) sequences for bats.

Methods

We performed germline gene annotation in three bat species using homologous genes and RSSs (Recombinational Signal Sequences) scanning method. Then we used the conserved C gene to construct the TCRβ chain receptor library of the Intermediate Horseshoe Bat. Bats' TCRβ data will be analyzed using MiXCR and constructed reference library.

Results

Regarding the annotation results, we found that the Pale Spear-nosed Bat has 37 members in the TRBV12 family, which is more than the total number of TRBV genes in the Greater Horseshoe Bat. The average number of unique TCRβ chain receptor sequences in each Intermediate Horseshoe Bat sample reached 24,904.

Discussion

The distinct variations in the distribution of TRBV genes among the three types of bats could have a direct impact on the diversity of the TCR repertoire, as evidenced by the presence of conserved amino acids that indicate the T-cell recognition of antigens in bats is MHC-restricted. The bats’ TCRβ repertoire is formed through the rearrangement of the V-D-J-C genes, with D-J/V-D deletions and insertions resulting in high diversity.

The methods and advances of adaptive immune receptors repertoire sequencing

The adaptive immune response is a powerful tool, capable of recognizing, binding to, and neutralizing a vast number of internal and external threats via T or B lymphatic receptors with widespread sets of antigen specificities. The emergence of high-throughput sequencing technology and bioinformatics provides opportunities for research in the fields of life sciences and medicine. The analysis and annotation for immune repertoire data can reveal biologically meaningful information, including immune prediction, target antigens, and effective evaluation. Continuous improvements of the immunological repertoire sequencing methods and analysis tools will help to minimize the experimental and calculation errors and realize the immunological information to meet the clinical requirements. That said, the clinical application of adaptive immune repertoire sequencing requires appropriate experimental methods and standard analytical tools. At the population cell level, we can acquire the overview of cell groups, but the information about a single cell is not obtained accurately. The information that is ignored may be crucial for understanding the heterogeneity of each cell, gene expression and drug response. The combination of high-throughput sequencing and single-cell technology allows us to obtain single-cell information with low-cost and high-throughput. In this review, we summarized the current methods and progress in this area.

Rhesus and cynomolgus macaque immunoglobulin heavy-chain genotyping yields comprehensive databases of germline VDJ alleles

Definition of the specific germline immunoglobulin (Ig) alleles present in an individual is a critical first step to delineate the ontogeny and evolution of antigen-specific antibody responses. Rhesus and cynomolgus macaques are important animal models for pre-clinical studies, with four main sub-groups being used: Indian- and Chinese-origin rhesus macaques and Mauritian and Indonesian cynomolgus macaques. We applied the (Ig) gene inference tool IgDiscover and performed extensive Sanger sequencing-based genomic validation to define germline VDJ alleles in these 4 sub-groups, comprising 45 macaques in total. There was allelic overlap between Chinese- and Indian-origin rhesus macaques and also between the two macaque species, which is consistent with substantial admixture. The island-restricted Mauritian cynomolgus population displayed the lowest number of alleles of the sub-groups, yet maintained high individual allelic diversity. These comprehensive databases of germline IGH alleles for rhesus and cynomolgus macaques provide a resource toward the study of B cell responses in these important pre-clinical models.

Systematic Profiling of Full-Length Ig and TCR Repertoire Diversity in Rhesus Macaque through Long Read Transcriptome Sequencing

The diversity of Ig and TCR repertoires is a focal point of immunological studies. Rhesus macaques (Macaca mulatta) are key for modeling human immune responses, placing critical importance on the accurate annotation and quantification of their Ig and TCR repertoires. However, because of incomplete reference resources, the coverage and accuracy of the traditional targeted amplification strategies for profiling rhesus Ig and TCR repertoires are largely unknown. In this study, using long read sequencing, we sequenced four Indian-origin rhesus macaque tissues and obtained high-quality, full-length sequences for over 6000 unique Ig and TCR transcripts, without the need for sequence assembly. We constructed, to our knowledge, the first complete reference set for the constant regions of all known isotypes and chain types of rhesus Ig and TCR repertoires. We show that sequence diversity exists across the entire variable regions of rhesus Ig and TCR transcripts. Consequently, existing strategies using targeted amplification of rearranged variable regions comprised of V(D)J gene segments miss a significant fraction (27-53% and 42-49%) of rhesus Ig/TCR diversity. To overcome these limitations, we designed new rhesus-specific assays that remove the need for primers conventionally targeting variable regions and allow single cell level Ig and TCR repertoire analysis. Our improved approach will enable future studies to fully capture rhesus Ig and TCR repertoire diversity and is applicable for improving annotations in any model organism.

Selection of potential cytokeratin-18 monoclonal antibodies following IGH repertoire evaluation in mice

Cytokeratin 18 (CK18), the main scaffold protein of keratinocyte, is distributed in epithelial cells. This structural protein maintains the integrity and continuity of epithelial tissue. Cytokeratin is also frequently used as an immunohistochemical marker of tumor growth. In recent years, immune repertoire (IR) evaluation using next-generation sequencing (NGS) have become increasingly efficient. Here we deep sequenced the mouse IR of the immunoglobulin heavy chain (IGH) after CK18 immunization. We comprehensively analyzed the IR based on complementarity determining region 3 (CDR3) abundance, germline gene usage polarization, clone diversity, and lineage. We found many convergence characteristics after CK18 immunization. Convergence represents a phenomenon that antigen stimulation or pathogen exposure induces the antigen specific clone expansion and enrichment. The convergence could be used for the immune evaluation and antibody screen. After immunization, the IGHV5 gene clusters became preponderant. The abundance and length of the most frequent CDR3 both increased, nevertheless the IR diversity level decreased. From the convergent IGH repertoires, we selected and expressed six antibodies with the most frequent CDR3s and IGH V-J combinations. The ELISA results suggested all screened six antibodies bound CK18 specifically. The most potential antibody had 9.424E-10M M affinity for the interaction with the CK18. Therefore, this is the NGS platform has been first used for anti-CK18 monoclonal antibodies (MAbs) discovery. These analyses methods could also be used for vaccine evaluation.

A Comprehensive Analysis of the T and B Lymphocytes Repertoire Shaped by HIV Vaccines

The exploitation of various human immunodeficiency virus type-1 (HIV-1) vaccines has posed great challenges for the researchers in precisely evaluating the vaccine-induced immune responses, however, the understanding of vaccination response suffers from the lack of unbiased characterization of the immune landscape. The rapid development of high throughput sequencing (HTS) makes it possible to scrutinize the extremely complicated immunological responses during vaccination. In the current study, three vaccines, namely N36, N51, and 5-Helix based on the HIV-1 gp41 pre-hairpin fusion intermediate were applied in rhesus macaques. We assessed the longitudinal vaccine responses using HTS, which delineated the evolutionary features of both T cell and B cell receptor repertoires with extreme diversities. Upon vaccination, we unexpectedly found significant discrepancies in the landscapes of T-cell and B-cell repertoires, together with the detection of significant class switching and the lineage expansion of the B cell receptor or immunoglobulin heavy chain (IGH) repertoire. The vaccine-induced expansions of lineages were further evaluated for mutation rate, lineage abundance, and lineage size features in their IGH repertoires. Collectively, these findings conclude that the N51 vaccine displayed superior performance in inducing the class-switch of B cell isotypes and promoting mutations of IgM B cells. In addition, the systematic HTS analysis of the immune repertoires demonstrates its wide applicability in enhancing the understanding of immunologic changes during pathogen challenge, and will guide the development, evaluation, and exploitation of new generation of diagnostic markers, immunotherapies, and vaccine strategies.

Profiling of the TCRβ repertoire in non-model species using high-throughput sequencing

In recent years, immune repertoire profiling with high-throughput sequencing (HTS) has advanced our understanding of adaptive immunity. However, fast progress in the field applied mostly to human and mouse research, with only few studies devoted to other model vertebrates. We present the first in-depth characterization of the T-cell receptor (TCR) repertoire in a non-model mammal (bank vole, Myodes glareolus), widely used in ecological and evolutionary research. We used RNA from spleens, 5′RACE and HTS to describe V and J segments of TCRβ, qualitatively characterize preferential V–J segment usage and CDR3 length distribution. Overall orthology to murine genes was preserved, with 11 J and 37 V genes found in voles (although 3 V genes lacked a close orthologue). Further, we implemented unique molecular identifiers for quantitative analysis of CDR3 repertoire with stringent error correction. A conservative, lower bound estimation of the TCRβ repertoire was similar to that found for mice (1.7–2.3 × 10⁵ clonotypes). We hope that by providing an easy-to-follow molecular protocol and on-line bioinformatics tools that do not require reference sequences (AmpliTCR and AmpliCDR3), we will encourage HTS immune repertoire profiling in other non-model vertebrates, thus opening new research avenues in e.g. comparative immunology, ecology and evolutionary biology.