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

Multi-compartmental diversification of neutralizing antibody lineages dissected in SARS-CoV-2 spike-immunized macaques

The continued evolution of SARS-CoV-2 underscores the need to understand qualitative aspects of the humoral immune response elicited by spike immunization. Here, we combine monoclonal antibody (mAb) isolation with deep B cell receptor (BCR) repertoire sequencing of rhesus macaques immunized with prefusion-stabilized spike glycoprotein. Longitudinal tracing of spike-sorted B cell lineages in multiple immune compartments demonstrates increasing somatic hypermutation and broad dissemination of vaccine-elicited B cells in draining and non-draining lymphoid compartments, including the bone marrow, spleen and, most notably, periaortic lymph nodes. Phylogenetic analysis of spike-specific monoclonal antibody lineages identified through deep repertoire sequencing delineates extensive intra-clonal diversification that shaped neutralizing activity. Structural analysis of the spike in complex with a broadly neutralizing mAb provides a molecular basis for the observed differences in neutralization breadth between clonally related antibodies. Our findings highlight that immunization leads to extensive intra-clonal B cell evolution where members of the same lineage can both retain the original epitope specificity and evolve to recognize additional spike variants not previously encountered.

Evolution of T cells in the cancer-resistant naked mole-rat

Naked mole-rats (NMRs) are best known for their extreme longevity and cancer resistance, suggesting that their immune system might have evolved to facilitate these phenotypes. Natural killer (NK) and T cells have evolved to detect and destroy cells infected with pathogens and to provide an early response to malignancies. While it is known that NMRs lack NK cells, likely lost during evolution, little is known about their T-cell subsets in terms of the evolution of the genes that regulate their function, their clonotypic diversity, and the thymus where they mature. Here we find, using single-cell transcriptomics, that NMRs have a large circulating population of γδT cells, which in mice and humans mostly reside in peripheral tissues and induce anti-cancer cytotoxicity. Using single-cell-T-cell-receptor sequencing, we find that a cytotoxic γδT-cell subset of NMRs harbors a dominant clonotype, and that their conventional CD8 αβT cells exhibit modest clonotypic diversity. Consistently, perinatal NMR thymuses are considerably smaller than those of mice yet follow similar involution progression. Our findings suggest that NMRs have evolved under a relaxed intracellular pathogenic selective pressure that may have allowed cancer resistance and longevity to become stronger targets of selection to which the immune system has responded by utilizing γδT cells.

Three immunizations with Novavax's protein vaccines increase antibody breadth and provide durable protection from SARS-CoV-2

The immune responses to Novavax's licensed NVX-CoV2373 nanoparticle Spike protein vaccine against SARS-CoV-2 remain incompletely understood. Here, we show in rhesus macaques that immunization with Matrix-MTM adjuvanted vaccines predominantly elicits immune events in local tissues with little spillover to the periphery. A third dose of an updated vaccine based on the Gamma (P.1) variant 7 months after two immunizations with licensed NVX-CoV2373 resulted in significant enhancement of anti-spike antibody titers and antibody breadth including neutralization of forward drift Omicron variants. The third immunization expanded the Spike-specific memory B cell pool, induced significant somatic hypermutation, and increased serum antibody avidity, indicating considerable affinity maturation. Seven months after immunization, vaccinated animals controlled infection by either WA-1 or P.1 strain, mediated by rapid anamnestic antibody and T cell responses in the lungs. In conclusion, a third immunization with an adjuvanted, low-dose recombinant protein vaccine significantly improved the quality of B cell responses, enhanced antibody breadth, and provided durable protection against SARS-CoV-2 challenge.

Host immunity associated with spontaneous suppression of viremia in therapy-naïve young rhesus macaques following neonatal SHIV infection

IMPORTANCE

Despite the advent of highly active anti-retroviral therapy, people are still dying from HIV-related causes, many of whom are children, and a protective vaccine or cure is needed to end the HIV pandemic. Understanding the nature and activation states of immune cell subsets during infection will provide insights into the immunologic milieu associated with viremia suppression that can be harnessed via therapeutic strategies to achieve a functional cure, but these are understudied in pediatric subjects. We evaluated humoral and adaptive host immunity associated with suppression of viremia in rhesus macaques infected soon after birth with a pathogenic SHIV. The results from our study provide insights into the immune cell subsets and functions associated with viremia control in young macaques that may translate to pediatric subjects for the design of future anti-viral strategies in HIV-1-infected infants and children and contribute to an understudied area of HIV-1 pathogenesis in pediatric subjects.

FLAIRR-Seq: A Method for Single-Molecule Resolution of Near Full-Length Antibody H Chain Repertoires

Current Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) using short-read sequencing strategies resolve expressed Ab transcripts with limited resolution of the C region. In this article, we present the near-full-length AIRR-seq (FLAIRR-seq) method that uses targeted amplification by 5' RACE, combined with single-molecule, real-time sequencing to generate highly accurate (99.99%) human Ab H chain transcripts. FLAIRR-seq was benchmarked by comparing H chain V (IGHV), D (IGHD), and J (IGHJ) gene usage, complementarity-determining region 3 length, and somatic hypermutation to matched datasets generated with standard 5' RACE AIRR-seq using short-read sequencing and full-length isoform sequencing. Together, these data demonstrate robust FLAIRR-seq performance using RNA samples derived from PBMCs, purified B cells, and whole blood, which recapitulated results generated by commonly used methods, while additionally resolving H chain gene features not documented in IMGT at the time of submission. FLAIRR-seq data provide, for the first time, to our knowledge, simultaneous single-molecule characterization of IGHV, IGHD, IGHJ, and IGHC region genes and alleles, allele-resolved subisotype definition, and high-resolution identification of class switch recombination within a clonal lineage. In conjunction with genomic sequencing and genotyping of IGHC genes, FLAIRR-seq of the IgM and IgG repertoires from 10 individuals resulted in the identification of 32 unique IGHC alleles, 28 (87%) of which were previously uncharacterized. Together, these data demonstrate the capabilities of FLAIRR-seq to characterize IGHV, IGHD, IGHJ, and IGHC gene diversity for the most comprehensive view of bulk-expressed Ab repertoires to date.

Large-Scale Complete Sequencing and Haplotyping of 1-10 kb DNA Molecules Using Short Massively Parallel Reads

In this chapter, we describe a simple, low-cost method for making many copies of a single DNA molecule (1-10 kb in length) as a concatemer on a long DNA strand. This can enable applications requiring high-quality contiguous sequence and haplotype data from long single DNA molecules at large scale.

Alternative splicing and genetic variation of mhc-e: implications for rhesus cytomegalovirus-based vaccines

Rhesus cytomegalovirus (RhCMV)-based vaccination against Simian Immunodeficiency virus (SIV) elicits MHC-E-restricted CD8+ T cells that stringently control SIV infection in ~55% of vaccinated rhesus macaques (RM). However, it is unclear how accurately the RM model reflects HLA-E immunobiology in humans. Using long-read sequencing, we identified 16 Mamu-E isoforms and all Mamu-E splicing junctions were detected among HLA-E isoforms in humans. We also obtained the complete Mamu-E genomic sequences covering the full coding regions of 59 RM from a RhCMV/SIV vaccine study. The Mamu-E gene was duplicated in 32 (54%) of 59 RM. Among four groups of Mamu-E alleles: three ~5% divergent full-length allele groups (G1, G2, G2_LTR) and a fourth monomorphic group (G3) with a deletion encompassing the canonical Mamu-E exon 6, the presence of G2_LTR alleles was significantly (p = 0.02) associated with the lack of RhCMV/SIV vaccine protection. These genomic resources will facilitate additional MHC-E targeted translational research.

Identification and Tracking of Alloreactive T Cell Clones in Rhesus Macaques Through the RM-scTCR-Seq Platform

T cell receptor (TCR) clonotype tracking is a powerful tool for interrogating T cell mediated immune processes. New methods to pair a single cell's transcriptional program with its TCR identity allow monitoring of T cell clonotype-specific transcriptional dynamics. While these technologies have been available for human and mouse T cells studies, they have not been developed for Rhesus Macaques (RM), a critical translational organism for autoimmune diseases, vaccine development and transplantation. We describe a new pipeline, 'RM-scTCR-Seq', which, for the first time, enables RM specific single cell TCR amplification, reconstruction and pairing of RM TCR's with their transcriptional profiles. We apply this method to a RM model of GVHD, and identify and track in vitro detected alloreactive clonotypes in GVHD target organs and explore their GVHD driven cytotoxic T cell signature. This novel, state-of-the-art platform fundamentally advances the utility of RM to study protective and pathogenic T cell responses.

Single-Cell-Based High-Throughput Ig and TCR Repertoire Sequencing Analysis in Rhesus Macaques

Recent advancements in microfluidics and high-throughput sequencing technologies have enabled recovery of paired H and L chains of Igs and VDJ and VJ chains of TCRs from thousands of single cells simultaneously in humans and mice. Despite rhesus macaques being one of the most well-studied model organisms for the human adaptive immune response, high-throughput single-cell immune repertoire sequencing assays are not yet available due to the complexity of these polyclonal receptors. We used custom primers that capture all known rhesus macaque Ig and TCR isotypes and chains that are fully compatible with a commercial solution for single-cell immune repertoire profiling. Using these rhesus-specific assays, we sequenced Ig and TCR repertoires in >60,000 cells from cryopreserved rhesus PBMCs, splenocytes, and FACS-sorted B and T cells. We were able to recover every Ig isotype and TCR chain, measure clonal expansion in proliferating T cells, and pair Ig and TCR repertoires with gene expression profiles of the same single cells. Our results establish the ability to perform high-throughput immune repertoire analysis in rhesus macaques at the single-cell level.

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.

Different adjuvanted pediatric HIV envelope vaccines induced distinct plasma antibody responses despite similar B cell receptor repertoires in infant rhesus macaques

Different HIV vaccine regimens elicit distinct plasma antibody responses in both human and nonhuman primate models. Previous studies in human and non-human primate infants showed that adjuvants influenced the quality of plasma antibody responses induced by pediatric HIV envelope vaccine regimens. We recently reported that use of the 3M052-SE adjuvant and longer intervals between vaccinations are associated with higher magnitude of antibody responses in infant rhesus macaques. However, the impact of different adjuvants in HIV vaccine regimens on the developing infant B cell receptor (BCR) repertoire has not been studied. This study evaluated whether pediatric HIV envelope vaccine regimens with different adjuvants induced distinct antigen-specific memory B cell repertoires and whether specific immunoglobulin (Ig) immunogenetic characteristics are associated with higher magnitude of plasma antibody responses in vaccinated infant rhesus macaques. We utilized archived preclinical pediatric HIV vaccine studies PBMCs and tissue samples from 19 infant rhesus macaques immunized either with (i) HIV Env protein with a squalene adjuvant, (ii) MVA-HIV and Env protein co-administered using a 3-week interval, (iii) MVA-HIV prime/ protein boost with an extended 6-week interval between immunizations, or (iv) with HIV Env administered with 3M-052-SE adjuvant. Frequencies of vaccine-elicited HIV Env-specific memory B cells from PBMCs and tissues were similar across vaccination groups (frequency range of 0.06–1.72%). There was no association between vaccine-elicited antigen-specific memory B cell frequencies and plasma antibody titer or avidity. Moreover, the epitope specificity and Ig immunogenetic features of vaccine-elicited monoclonal antibodies did not differ between the different vaccine regimens. These data suggest that pediatric HIV envelope vaccine candidates with different adjuvants that previously induced higher magnitude and quality of plasma antibody responses in infant rhesus macaques were not driven by distinct antigen-specific memory BCR repertoires.

Comprehensive transcriptome characterization of Grus japonensis using PacBio SMRT and Illumina sequencing

The red-crowned crane (Grus japonensis) is an endangered species distributed across southeast Russia, northeast China, Korea, and Japan. Here, we sequenced for the first time the full-length unreferenced transcriptome of red-crowned crane mixed samples using a PacBio Sequel platform. A total of 359,136 circular consensus sequences (CCS) were obtained via clustering to remove redundancy. A total of 303,544 full-length non-chimeric sequences were identified by judging whether CCS contained 5' and 3' adapters, and the poly(A) tail. Eight samples were sequenced using Illumina, and PacBio sequencing data were corrected according to the collected Illumina data to obtain more accurate full-length transcripts. A total of 4,100 long non-coding RNAs, 13,115 simple sequences repeat loci and 29 transcription factor families were identified. The expression of lncRNAs and TFs in pancreas was lowest comparing with other tissues. Many enriched immune-related transmission pathways (MHC and IL receptors) were identified in the spleen. This study will contribute to a better understanding of the gene structure and post-transcriptional regulatory network, and provide references for future studies on red-crowned cranes.

Ultra-accurate microbial amplicon sequencing with synthetic long reads

BACKGROUND

Out of the many pathogenic bacterial species that are known, only a fraction are readily identifiable directly from a complex microbial community using standard next generation DNA sequencing. Long-read sequencing offers the potential to identify a wider range of species and to differentiate between strains within a species, but attaining sufficient accuracy in complex metagenomes remains a challenge.

METHODS

Here, we describe and analytically validate LoopSeq, a commercially available synthetic long-read (SLR) sequencing technology that generates highly accurate long reads from standard short reads.

RESULTS

LoopSeq reads are sufficiently long and accurate to identify microbial genes and species directly from complex samples. LoopSeq perfectly recovered the full diversity of 16S rRNA genes from known strains in a synthetic microbial community. Full-length LoopSeq reads had a per-base error rate of 0.005%, which exceeds the accuracy reported for other long-read sequencing technologies. 18S-ITS and genomic sequencing of fungal and bacterial isolates confirmed that LoopSeq sequencing maintains that accuracy for reads up to 6 kb in length. LoopSeq full-length 16S rRNA reads could accurately classify organisms down to the species level in rinsate from retail meat samples, and could differentiate strains within species identified by the CDC as potential foodborne pathogens.

CONCLUSIONS

The order-of-magnitude improvement in length and accuracy over standard Illumina amplicon sequencing achieved with LoopSeq enables accurate species-level and strain identification from complex- to low-biomass microbiome samples. The ability to generate accurate and long microbiome sequencing reads using standard short read sequencers will accelerate the building of quality microbial sequence databases and removes a significant hurdle on the path to precision microbial genomics. Video abstract.