Combining genotypes and T cell receptor distributions to infer genetic loci determining V(D)J recombination probabilities

Evolutionary trade-offs constraining the MHC gene expansion: beyond simple TCR depletion model

The immune system is as much shaped by the pressure of pathogens as it is by evolutionary trade-offs that constrain its structure and function. A perfect example comes from the major histocompatibility complex (MHC), molecules that initiate adaptive immune response by presentation of foreign antigens to T cells. The remarkable, population-level polymorphism of MHC genes is assumed to result mainly from a co-evolutionary arms race between hosts and pathogens, while the limited, within-individual number of functional MHC loci is thought to be the consequence of an evolutionary trade-off between enhanced pathogen recognition and excessive T cell depletion during negative selection in the thymus. Certain mathematical models and infection studies suggest that an intermediate individual MHC diversity would thus be optimal. A recent, more direct test of this hypothesis has shown that the effects of MHC diversity on T-cell receptor (TCR) repertoires may differ between MHC classes, supporting the depletion model only for MHC class I. Here, we used the bank vole (Myodes=Cletronomys glareolus), a rodent species with variable numbers of expressed MHC genes, to test how an individual MHC diversity influences the proportions and TCR repertoires of responding T cell subsets. We found a non-linear relationship between MHC diversity and T cell proportions (with intermediate MHC numbers coinciding with the largest T cell proportions), perhaps reflecting an optimality effect of balanced positive and negative thymic selection. The association was strongest for the relationship between MHC class I and splenic CD8+ T cells. The CD8+ TCR richness alone was unaffected by MHC class I diversity, suggesting that MHC class I expansion may be limited by decreasing T cell counts, rather than by direct depletion of TCR richness. In contrast, CD4+ TCR richness was positively correlated with MHC class II diversity, arguing against a universal TCR depletion. It also suggests that different evolutionary forces or trade-offs may limit the within-individual expansion of the MHC class II loci.

Immunosequencing and Profiling of T Cells at the Maternal-Fetal Interface of Women with Preterm Labor and Chronic Chorioamnionitis

T cells are implicated in the pathophysiology of preterm labor and birth, the leading cause of neonatal morbidity and mortality worldwide. Specifically, maternal decidual T cells infiltrate the chorioamniotic membranes in chronic chorioamnionitis (CCA), a placental lesion considered to reflect maternal anti-fetal rejection, leading to preterm labor and birth. However, the phenotype and TCR repertoire of decidual T cells in women with preterm labor and CCA have not been investigated. In this study, we used phenotyping, TCR sequencing, and functional assays to elucidate the molecular characteristics and Ag specificity of T cells infiltrating the chorioamniotic membranes in women with CCA who underwent term or preterm labor. Phenotyping indicated distinct enrichment of human decidual effector memory T cell subsets in cases of preterm labor with CCA without altered regulatory T cell proportions. TCR sequencing revealed that the T cell repertoire of CCA is characterized by increased TCR richness and decreased clonal expansion in women with preterm labor. We identified 15 clones associated with CCA and compared these against established TCR databases, reporting that infiltrating T cells may possess specificity for maternal and fetal Ags, but not common viral Ags. Functional assays demonstrated that choriodecidual T cells can respond to maternal and fetal Ags. Collectively, our findings provide, to our knowledge, novel insight into the complex processes underlying chronic placental inflammation and further support a role for effector T cells in the mechanisms of disease for preterm labor and birth. Moreover, this work further strengthens the contribution of adaptive immunity to the syndromic nature of preterm labor and birth.

IGHV allele similarity clustering improves genotype inference from adaptive immune receptor repertoire sequencing data

In adaptive immune receptor repertoire analysis, determining the germline variable (V) allele associated with each T- and B-cell receptor sequence is a crucial step. This process is highly impacted by allele annotations. Aligning sequences, assigning them to specific germline alleles, and inferring individual genotypes are challenging when the repertoire is highly mutated, or sequence reads do not cover the whole V region. Here, we propose an alternative naming scheme for the V alleles, as well as a novel method to infer individual genotypes. We demonstrate the strengths of the two by comparing their outcomes to other genotype inference methods. We validate the genotype approach with independent genomic long-read data. The naming scheme is compatible with current annotation tools and pipelines. Analysis results can be converted from the proposed naming scheme to the nomenclature determined by the International Union of Immunological Societies (IUIS). Both the naming scheme and the genotype procedure are implemented in a freely available R package (PIgLET https://bitbucket.org/yaarilab/piglet). To allow researchers to further explore the approach on real data and to adapt it for their uses, we also created an interactive website (https://yaarilab.github.io/IGHV_reference_book).

Robust and prototypical immune responses toward COVID-19 vaccine in First Nations peoples are impacted by comorbidities

High-risk groups, including Indigenous people, are at risk of severe COVID-19. Here we found that Australian First Nations peoples elicit effective immune responses to COVID-19 BNT162b2 vaccination, including neutralizing antibodies, receptor-binding domain (RBD) antibodies, SARS-CoV-2 spike-specific B cells, and CD4+ and CD8+ T cells. In First Nations participants, RBD IgG antibody titers were correlated with body mass index and negatively correlated with age. Reduced RBD antibodies, spike-specific B cells and follicular helper T cells were found in vaccinated participants with chronic conditions (diabetes, renal disease) and were strongly associated with altered glycosylation of IgG and increased interleukin-18 levels in the plasma. These immune perturbations were also found in non-Indigenous people with comorbidities, indicating that they were related to comorbidities rather than ethnicity. However, our study is of a great importance to First Nations peoples who have disproportionate rates of chronic comorbidities and provides evidence of robust immune responses after COVID-19 vaccination in Indigenous people.

Statistical inference reveals the role of length, GC content, and local sequence in V(D)J nucleotide trimming

To appropriately defend against a wide array of pathogens, humans somatically generate highly diverse repertoires of B cell and T cell receptors (BCRs and TCRs) through a random process called V(D)J recombination. Receptor diversity is achieved during this process through both the combinatorial assembly of V(D)J-genes and the junctional deletion and insertion of nucleotides. While the Artemis protein is often regarded as the main nuclease involved in V(D)J recombination, the exact mechanism of nucleotide trimming is not understood. Using a previously published TCRβ repertoire sequencing data set, we have designed a flexible probabilistic model of nucleotide trimming that allows us to explore various mechanistically interpretable sequence-level features. We show that local sequence context, length, and GC nucleotide content in both directions of the wider sequence, together, can most accurately predict the trimming probabilities of a given V-gene sequence. Because GC nucleotide content is predictive of sequence-breathing, this model provides quantitative statistical evidence regarding the extent to which double-stranded DNA may need to be able to breathe for trimming to occur. We also see evidence of a sequence motif that appears to get preferentially trimmed, independent of GC-content-related effects. Further, we find that the inferred coefficients from this model provide accurate prediction for V- and J-gene sequences from other adaptive immune receptor loci. These results refine our understanding of how the Artemis nuclease may function to trim nucleotides during V(D)J recombination and provide another step toward understanding how V(D)J recombination generates diverse receptors and supports a powerful, unique immune response in healthy humans.

A genome-wide association study for allergen component sensitizations identifies allergen component-specific and allergen protein group-specific associations

Background

Allergic diseases are some of the most common diseases worldwide. Genome-wide association studies (GWASs) have been conducted to elucidate the genetic factors of allergic diseases. However, no GWASs for allergen component sensitization have been performed.

Objective

We sought to detect genetic variants associated with differences in immune responsiveness against allergen components.

Methods

The participants of the present study were recruited from the Tokyo Children's Health, Illness, and Development study, and allergen component-specific IgE level at age 9 years was measured by means of allergen microarray immunoassays. We performed GWASs for allergen component sensitization against each allergen (single allergen component sensitization, number of allergen components analyzed, n = 31), as well as against allergen protein families (allergen protein group sensitization, number of protein groups analyzed, n = 16).

Results

We performed GWAS on 564 participants of the Tokyo Children's Health, Illness, and Development study and found associations between Amb a 1 sensitization and the immunoglobulin heavy-chain variable gene on chromosome 14 and between Phl p 1 sensitization and the HLA class II region on chromosome 6 (P < 5.0 × 10-8). A GWAS-significant association was also observed between the HLA class II region and profilin sensitization (P < 5.0 × 10-8).

Conclusions

Our data provide the first demonstration of genetic risk for allergen component sensitization and show that this genetic risk is related to immune response genes including immunoglobulin heavy-chain variable gene and HLA.

Targeted long-read sequencing facilitates phased diploid assembly and genotyping of the human T cell receptor alpha, delta, and beta loci

T cell receptors (TCRs) recognize peptide fragments presented by the major histocompatibility complex (MHC) and are critical to T cell-mediated immunity. Recent data have indicated that genetic diversity within TCR-encoding gene regions is underexplored, limiting understanding of the impact of TCR loci polymorphisms on TCR function in disease, even though TCR repertoire signatures (1) are heritable and (2) associate with disease phenotypes. To address this, we developed a targeted long-read sequencing approach to generate highly accurate haplotype resolved assemblies of the TCR beta (TRB) and alpha/delta (TRA/D) loci, facilitating the genotyping of all variant types, including structural variants. We validate our approach using two mother-father-child trios and 5 unrelated donors representing multiple populations. This resulted in improved genotyping accuracy and the discovery of 84 undocumented V, D, J, and C alleles, demonstrating the utility of this framework for improving our understanding of TCR diversity and function in disease.