Sex-Biased Aging Effects on Ig Somatic Hypermutation Targeting

Ang Cui, Daniel G. Chawla and Steven H. Kleinstein
J Immunol January 1, 2021, 206 (1) 101-108; DOI: https://doi.org/10.4049/jimmunol.2000576

Key Points

  • SHM hot- and cold-spot biases are largely consistent between different age and sex groups.

  • A shift in SHM targeting is observed between young and older male subjects.

  • This shift in SHM targeting influences the Ig amino acid composition.

Abstract

Aged individuals, particularly males, display an impaired level of Ab response compared with their younger counterparts, yet the molecular mechanisms responsible for the discrepancy are not well understood. We hypothesize that some of this difference may be linked to B cell somatic hypermutation (SHM) targeting, including error-prone DNA repair activities that are crucial to Ab diversification. To examine the effects of aging on SHM targeting, we analyzed B cell Ig repertoire sequences from 27 healthy male and female human subjects aged 20–89. By studying mutation patterns based on 985,069 mutations obtained from 123,415 sequences, we found that the SHM mutability hierarchies on microsequence motifs (i.e., SHM hot/cold spots) are mostly consistent between different age and sex groups. However, we observed a lower frequency in mutations involving Phase II SHM DNA repair activities in older males, but not in females. We also observed, from a separate study, a decreased expression level of DNA mismatch repair genes involved in SHM in older individuals compared with younger individuals, with larger fold changes in males than in females. Finally, we showed that the balance between Phase I versus Phase II SHM activities impacts the resulting Ig phenotypes. Our results showed that the SHM process is altered in some older individuals, providing insights into observed clinical differences in immunologic responses between different age and sex groups.

Introduction

Age and sex are important sources of variation in immune responses. Older individuals tend to be more vulnerable to bacterial and viral infections, more prone to developing cancer and autoimmune diseases, and less likely to respond to vaccinations (1–3). The effects of aging on the immune system, however, are not uniform between sexes (4). Females tend to display a more robust immunologic response to infections and vaccinations, yet females are also at a higher risk of developing autoimmune diseases (5–8). Although the age and sex biases have long been observed, the molecular mechanisms causing these differences remain unclear. Gaining an understanding of the molecular mechanisms is crucial for designing more targeted therapies and vaccinations for different age and sex groups.

Among many immune cells that experience immunosenescence, B cells are of great interest because of their central role in humoral immunity. B cells produce Abs that defend us from a wide range of extracellular pathogens or bind to cells to elicit Ab-dependent cytotoxicity. Moreover, most vaccines for viral and bacterial infections to date target B cells and use their memory response to prevent future infections. Previous studies have shown that although the quantity of Abs does not alter significantly between young and aged individuals, older individuals display a decrease in the affinity, specificity, and diversity of the Abs they produce (9–12). Such observations prompted us to examine the mechanisms governing the Ab diversification process.

To generate high-affinity BCRs, which are the membrane-bound form of Abs, the Ig loci in B cells undergo a number of diversification steps, including V(D)J recombination and somatic hypermutation (SHM) (13). SHM, a process that occurs in germinal centers (GCs) in secondary lymphoid organs, diversifies BCRs by introducing point mutations into the Ig genes at a high rate (14). B cells with mutations leading to a high affinity to Ags survive and clonally expand in a selection process. SHM is initiated when the enzyme activation-induced cytidine deaminase (AID) converts cytosine (C) to uracil (U) in the Ig loci in B cells. Following this, three scenarios may occur. In Phase Ia, a simple replication leads to a C-to-thymine (T) or a guanine (G)-to-adenine (A) transition on the complement strand. In Phase Ib, the activation of UNG-induced short-patch base-excision repair (BER) pathway introduces mutation at the targeted C or G site (hereafter written as C/G site). Phase II results in the activation of long-patch BER and mismatch repair (MMR) pathways, leading to mutations on neighboring bases that can be on any nucleotide, including the A and T site (hereafter written as A/T site). Although stochastic, SHM preferentially targets certain DNA motifs (hot spots) while avoiding others (cold spots) (15). A previous study shows a decrease in AID expression level in older individuals, suggesting that a lower degree of SHM initiation may contribute to the compromised immunologic responses in these individuals (16). However, there has been a lack of investigation on sex biases or the impact of aging on SHM targeting, which dictates the types of mutations resulting from SHM.

With the advancement of high-throughput technologies, we can now profile immune repertoires using adaptive immune receptor repertoire sequencing (AIRR-seq) at an unprecedented depth (17, 18). Such technologies involve sequencing the Ig loci in DNA or RNA, enabling us to examine mutation patterns in BCRs in great detail. For example, it has been observed that aging is associated with longer CDR3 regions (the most diverse region in B cell Ig genes), an accumulation of more highly mutated IgM and IgG Ig genes, larger persistent clones in B cells in peripheral blood, and diminished intralineage diversification (19–21). Previously, we used such datasets to study SHM-targeting patterns across individuals and species at a high resolution (22, 23). We hypothesized that such approaches may also provide insights into the effects of age and sex in the SHM-targeting process.

In this study, we aimed to compare and contrast the SHM-targeting patterns between different age and sex groups. We analyzed high-throughput B cell Ig sequences from peripheral blood of 27 healthy individuals from two age groups: younger individuals (20–31 y old) and older individuals (61–89 y old) (19). As a validation, we examined mutations in Ig genes in whole-transcriptome RNA-sequencing (RNA-seq) data from 147 individuals aged 20–70 in the Genotype-Tissue Expression (GTEx) project (24). We examined mutations involved in Phase I versus Phase II SHM and found a significant decrease in Phase II–induced mutations in older males. Such observations lead to the hypothesis that older individuals, particularly males, experience alterations in SHM activities, which may contribute to their lower level of immunologic responses.

Materials and Methods

Repertoire-sequencing data processing and clonotype assignment

The high-throughput Ig repertoire AIRR-seq dataset was published in (19) by Wang et al. Raw high-throughput sequencing reads were quality controlled, assembled, and filtered using the REpertoire Sequencing TOolkit (pRESTO) (25). The IMGT/High-VQUEST tool was used to assign germline V(D)J segments and determine functionality (26). The Change-O command line tool was used to partition sequences into clonal groups (27). V(D)J sequences were assigned into clones based on identical IGHV gene, IGHJ gene, and junction length, with a weighted intraclonal distance threshold of 10 using the substitution probabilities previously described (28). Mutations were defined as nucleotides that were different from the germline sequence. A sliding-window approach was used to filter out sequences with highly clustered mutations, defined by six or more nucleotide differences in 10 consecutive nucleotides. Sequences with <5% mutation frequency were filtered out to reduce the effect of sequencing errors in the analysis.

Targeting and substitution models

The 1-mer substitution model was generated using the shazam R package version 0.1.1 in the Change-O suite (22, 27). In brief, the substitution matrix was defined by calculating the frequency of substitutions for each nucleotide. The 5-mer targeting model used for the construction of synonymous 5-mer models was previously described (22).

GTEx RNA-seq data processing

The raw whole-transcriptome RNA-seq data were collected by the GTEx project (v6) (24). The IgH locus was isolated and mapped to the germline Ig genes published by IMGT via IgBLAST (29). Reads with more than 65 nt (out of 75 nt) aligned were used for mutation analyses. To maximize the number of mutations we could obtain per individual, mutations from different tissues were aggregated for each individual. The A/T mutation frequency was normalized by the total number of A and T nucleotides in the germline sequences. Individuals with >0.8 or <0.2 A/T mutation frequency were considered outliers and excluded from the analysis.

Gene expression analysis

Gene expression data of healthy individuals of different age and sex groups were obtained from Gene Expression Omnibus accession number GSE65442 (30). Subject metadata were accessed from ImmPort SDY404. The baseline data of B cells prior to vaccination from the 2011–2012 season, which has at least four individuals in each age and sex group, were selected for the analysis. The gene expression profiles were quantile normalized using the R package limma. The genes involved in SHM repair activities, including AICDA, EXO1, MSH2, MSH6, PMS2, POLH, POLI, POLK, POLL, POLQ, REV1, UNG, and UNG2, were considered in the initial analysis (13). Probes with low detection signals (i.e., those with <90% of samples having a gene expression detection p value <0.01) were eliminated from downstream analyses.

Criteria for frailty

The operational definition of frailty was based on the following criteria: 1) low grip strength, 2) low energy, 3) slowed walking speed, 4) low physical activity, and 5) unintentional weight loss. Individuals with three to five symptoms were defined as frail.

Results

B cell Ig repertoire–sequencing data

To assess the effect of aging on SHM targeting, we analyzed a publicly available B cell Ig repertoire–sequencing (AIRR-seq) data from 27 healthy individuals (hereafter referred to as the “AIRR-seq” dataset) (19). The study includes 10 young adults (20–31 y old), five of whom are females, and 17 older adults (61–89 y old), 12 of whom are females. Within each age group, male and female subjects were of comparable ages (young: p = 0.65, older: p = 0.83; Welch t test).

High-throughput next-generation sequencing was used to profile the Ig sequences. For each individual, three replicates were obtained from different sequencing methods or platforms: genomic data sequenced at Stanford Genomic Center (DNA-GC), genomic data sequenced on the Roche 454 platform (DNA-Roche), and cDNA data from RNAs. We analyzed each replicate separately to prevent potential biases that arise from sequencing methods. A total of 123,415 unique sequences passed quality control (Materials and Methods) and were used for the downstream analyses (Table I). IMGT was used to determine the functionality of the Ig sequences (26). Nonfunctional sequences were defined as sequences containing out-of-frame junction region and were presumably not having participated in selection. To identify independent SHM events, we partitioned the sequences into groups that were likely to be clonally-related (Materials and Methods), and identified a set of distinct mutations from each clone. We collected a total of 985,069 SHMs, where mutations were defined as the nucleotides that were different from the inferred germline sequences in the V region. From the RNA replicate in which the isotype information is available, we found that the majority of the mutations were obtained from sequences of the IgG, IgA, and IgM isotypes. To study intrinsic SHM-targeting patterns in the absence of selection pressure, we analyzed two types of mutations from each replicate: 1) silent mutations (i.e., mutations that do not alter the amino acid sequence) from functional sequences and 2) both silent and replacement mutations from nonfunctional sequences (i.e., sequences that do not encode proteins).

View this table:
Table I. Sequencing data in the AIRR-seq dataset

Aging has little effect on SHM-targeting hierarchy on microsequence motifs

To assess whether aging influences the intrinsic biases of SHM targeting (i.e., hot/cold spots), we built targeting models based on microsequence motifs. To capture the classic SHM hot spot (WRC/GYW, where W = [A, T], R = [G, A], Y = [C, T]) on both the forward and reverse strands, the model was based on 5-mer motifs including two bases up- and two bases downstream of the mutated base (Materials and Methods) (22). The mutability of each 5-mer was defined as the probability of the central base being targeted by SHM relative to all other motifs. Because 5-mer mutability hierarchies are highly similar across individuals (22), we built a single targeting model for each age and sex group. We then used “hedgehog” plots to visualize the targeting biases (i.e., hot/cold-spot motifs) (Fig. 1A). In all age and sex groups, we observed high mutability in classic hot spots and low mutability in classic cold spots, consistent with our expectations. We compared the targeting models between age and sex groups and found that the mutability hierarchies were highly consistent between different ages and sexes (Spearman ρ > 0.92 in silent mutations in functional sequences, and ρ > 0.79 in both silent and replacement mutations in nonfunctional sequences in the RNA replicate) (Fig. 1B, 1C). The results suggest that age and sex have little effect on the hierarchy of SHM targeting on microsequence motifs.

FIGURE 1.
FIGURE 1.

SHM mutability hierarchy on microsequence motifs is consistent across age and sex groups. (A) SHM-targeting model for each age and sex group. Hedgehog plots depict the relative targeting of 5-mer motifs centered on the bases A, T, C, and G, with each center base in an individual circle. Classic hot spots are colored in red (WRC/GYW, where W = [A, T], R = [G, A], Y = [C, T]) and green (WA/TW), whereas cold spots are colored in blue (SYC/GRS, where S = [C, G]). (B) Scatterplot of 5-mer mutabilities between each age and sex group built on silent mutations in functional sequences. (C) Scatterplot of 5-mer mutabilities between each age and sex group built on all mutations in nonfunctional sequences. Each dot represents the mutability value of a 5-mer motif. R denotes Spearman correlation coefficient.

Aging is associated with a shift in SHM repair activities in males

In SHM, C/G nucleotides are targeted in both Phase I and Phase II, whereas A/T nucleotides are targeted in Phase II only. To detect whether there is a shift in SHM phases between age and sex groups, we examined the balance between mutations at C/G nucleotides versus mutations at A/T nucleotides. For each individual, we computed the proportion of mutations that occur on A/T bases, normalized by the number of each base in the germline sequences (hereafter written as A/T mutability). Across all technical replicates, we observed a significant decrease in A/T mutability in older males in nonfunctional sequences, with at least four of the five older males having lower mutations on A/T than all of the younger males (Fig. 2). We also observed a decrease in A/T mutability in older males in silent mutations from functional sequences, and the difference was statistically significant in the RNA replicate. In the male group, there may be a continuous decrease in A/T mutability with age (r < −0.70, p < 0.03, Pearson product-moment correlation in the RNA replicate) (Supplemental Fig. 1), although the level of continuity was unclear in the 31- to 60-y-old age group because of a lack of subjects. In females, the older individuals had comparable mutation patterns as the younger individuals. The older males also exhibited lower A/T mutabilities compared with older females across all replicates.

FIGURE 2.
FIGURE 2.

Older males display a lower level of A/T mutability in SHM. A/T mutability was calculated by normalizing the frequency of unselected mutations occurring on A or T nucleotides by the frequency of A or T bases in germline Ig sequences for each individual in the AIRR-seq dataset. Unselected mutations were obtained from silent mutations in functional sequences (Functional) or all mutations in nonfunctional sequences (Nonfunctional). Boxes represent the median and interquartile range of A/T mutabilities in each age and sex group. The analysis was performed in all three technical replicates: RNA (left), DNA-GC (middle), DNA-Roche (right). The p values were obtained from Welch t test.

We used an RNA-seq dataset published by the GTEx project to validate the results (hereafter referred to as the “GTEx RNA-seq” dataset) (24). This dataset contains RNA sequences from the entire transcriptome across many cell types, and therefore includes a limited number of mutated B cell Ig transcripts. We extracted the IgH locus and mapped the cDNA sequences to germline Ig genes using IgBLAST (Materials and Methods). The resulting data contained 11,312 mutations from 147 individuals. To maximize the number of mutations we can obtain per individual, sequences from different tissue sites were aggregated. We, again, calculated the A/T mutability by normalizing the proportion of mutations on A/T bases by the number of A/T nucleotides in the germline genes. Although the data contained a large amount of variation because of a small number of mutations that can be extracted from such types of data, we observed a similar trend as before: both sex groups displayed a decrease in A/T mutability with age, but the males had a much more prominent decrease (Pearson r = −0.24, p = 0.023 in males versus Pearson r = −0.13, p = 0.363 in females) (Supplemental Fig. 2). By applying linear models, we observed that A/T mutabilities in males decreased at a slightly faster rate (a mutability value of 0.5826 at age 0, decreasing by 0.0018 per year) compared with their female counterparts (a mutability value of 0.5291 at age 0, decreasing by 0.0011 per year). The decrease in A/T (versus C/G) mutability in older males suggests that these individuals have fewer mutations involved in Phase II SHM, which are induced by the activation of the long-patch BER and MSH2/MSH6 DNA MMR pathway that drives mutation spreading from the original AID-targeted C-to-U lesion to neighboring bases.

To determine whether there is also a shift between Phase Ia and Ib SHM repair pathways (i.e., simple replication versus UNG-induced short-patch BER at the AID-targeted C/G bases), we examined the frequency of C-to-T and G-to-A mutations (Phase Ia mutations) with respect to all mutations on C or G bases (Phase I mutations) in the AIRR-seq dataset. To observe this, we computed the frequency of each substitution (Supplemental Fig. 3A) and normalized them by the total number of substitutions from that base (Supplemental Fig. 3B). Although there appeared to be an insignificant elevation of Phase Ia mutations before normalization, this elevation disappeared after normalization. We tested the hypothesis on the GTEx RNA-seq dataset, and again did not see any correlation between age and Phase Ia mutation frequency. Thus, our analysis does not support a shift between Phase Ia and Phase Ib with age in these datasets.

Older males display a decrease in gene expression levels in molecules involved in SHM

The observation of a decreased mutation frequency occurring on A and T bases may suggest that the DNA repair mechanisms involved in SHM are altered in some older individuals. To assess whether DNA repair molecules involved in SHM have a different gene expression level, we analyzed the gene expression profiles of peripheral blood B cells collected from another cohort of healthy individuals (see Materials and Methods; hereafter referred to as the “gene expression” dataset). This cohort includes 16 younger (including seven male) and 14 older (including four male) subjects. B cells were sorted from blood samples and the gene expression levels of 47,260 probes were measured (30).

We examined the genes involved in SHM that were measured with significant detection level on the microarray (detection p < 0.01 for >90% of the samples), including MSH6, PMS2, REV1, and UNG. We compared the expression levels between the younger and older individuals in each sex group and found that these genes consistently displayed lower expression values in older males compared with younger males (Fig. 3). In particular, REV1 exhibited the largest gap between age groups (p = 0.01, Welch t test). In females, the expressions were comparable between age groups, although the older individuals had slightly lower expressions in most of the genes. Next, we examined whether the decrease in SHM repair activities was linked to health state. We divided the older individuals into “healthy” and “frail” categories based on the criteria described in Materials and Methods. There were a total of 11 healthy and three frail older individuals. Interestingly, the frail male showed lowest expression in multiple genes among the male group, although such a trend was not observed in the female group. These findings of gene expression changes in SHM molecules further demonstrate that older males may experience alterations in SHM-targeting activities.

FIGURE 3.
FIGURE 3.

Older males display a decrease in the expression of DNA repair genes involved in SHM. Log-transformed expressions of SHM repair genes with high detection levels in the gene expression dataset (described in Materials and Methods). Housekeeping genes were included as a control. Boxes represent the median and interquartile range of gene expressions in each age and sex group. The p values were obtained from Welch t test.

Balance between Phase I versus Phase II SHM repair activities impacts Ig phenotype

We hypothesized that an alteration in SHM targeting at the nucleotide level would affect the resulting Ig phenotypes and, consequently, contribute to a different immunological response. To detect the impact of a shift in A/T mutability on Ig amino acid composition, we examined the characteristics of amino acid substitutions by computationally translating the Ig nucleotide sequences into amino acids. We hypothesized that a decrease in A/T mutability would result in a decrease in the amino acid substitutions that occur because of mutations at A/T bases. Therefore, we examined a list of amino acid substitutions that would occur from a single nucleotide change unambiguously occurring at an A/T site (Table II). For example, the change from phenylalanine (Phe) to serine (Ser) was included in the list because the nucleotide T must be converted to C to induce this amino acid exchange in a single mutational event. In contrast, the change from serine to threonine was not included because it could be induced by mutations on T or G. There are 57 aa substitutions that match the criteria, with the most common ones being tyrosine (Tyr) to Phe or histidine (His), lysine (Lys) to arginine (Arg), Val to Ala, and Ser to Gly.

View this table:
Table II. Amino acid substitutions that are influenced by mutations on A/T nucleotides

For each individual, we calculated the proportion of such amino acid substitutions on functional sequences (i.e., sequences that encode BCRs) and compared it to the A/T mutability in nonfunctional sequences, which reflects the level of SHM-targeting activities not influenced by selection pressure. We found that a lower A/T mutability was significantly correlated with a lower frequency of amino acid substitutions that result from mutations at A/T sites for all replicates in the AIRR-seq dataset (RNA replicate: Spearman ρ = 0.82, p = 2 × 10−6) (Fig. 4, Supplemental Fig. 4). We observed that younger individuals clustered at the top right corner, which represents a higher A/T mutability and a higher frequency of amino acid substitutions due to such mutations at the nucleotide level. Older males and some older females occupied the bottom left corner, indicating a decrease in the proportion in amino acid substitutions involving A/T mutations. These results indicate that a shift in Phase I versus Phase II mutations affects the resulting Ig phenotype even if they are under selection pressure. As different amino acids have distinct properties (i.e., sizes, charges, and hydrophobicity), alterations in Ig amino acid composition would potentially impact immunologic responses, such as causing changes in Ig-binding affinity.

FIGURE 4.
FIGURE 4.

Changes in SHM targeting affects Ig phenotype. Comparison between A/T mutability found in nonfunctional sequences (reflecting intrinsic SHM-targeting biases) and the proportion of amino acid substitutions that are exclusively caused by A/T nucleotide changes in functional sequences (reflecting properties of amino acid sequences encoding functional BCRs) in the AIRR-seq dataset. Each point represents an individual; their corresponding age group is indicated by shade, and sex group is indicated by shape. The analysis was performed in all three technical replicates: RNA (shown in this figure), DNA-GC, and DNA-Roche (Supplemental Fig. 4).

Discussion

Elucidating age and sex discrepancies in immune responses at the molecular level can provide mechanistic insights to the differences in immune responses observed in clinical settings. Although it is well established that males and females of different age groups have different immune responses, the mechanisms responsible for the lack of Ab diversity and potency for certain groups are yet to be clarified. In this study, we analyzed high-throughput sequencing and gene expression datasets to study an Ab diversification process in different age and sex groups. We analyzed a large number of mutations gathered from high-throughput Ig-sequencing data, and used 1-mer and 5-mer microsequence-targeting frameworks to show a shift in A/T versus C/G targeting in older males. We used RNA-seq data from 147 individuals to validate the patterns observed. We also used gene expression profiles of an independent study to show that males have decreased levels of expression in several genes crucial to SHM-targeting activities. Finally, we showed that a shift in the balance between Phase I and Phase II SHM at the nucleotide level measured in the absence of selection pressure influences the amino acid compositions of functional Igs, potentially altering immunologic responses.

It has been observed that the Abs generated by older individuals have a lower level of affinity, specificity, and diversity. However, the molecular processes responsible for such declines have not been elucidated. Previous repertoire analyses have shown alterations in repertoire properties in older individuals, such as CDR3 lengths and clonal sizes. In this study, we investigated DNA repair pathways involved in SHM, processes upstream of these observations, and found that the mutations involved in different phases of SHM repair pathways are altered in some older individuals. Using one of the cohorts analyzed in this study (19), we previously found a decrease in WA hot spot targeting with age (31). However, unlike the sex-specific differences observed in this study in the targeting of A/T more generally, these hot spot differences were similar in males and females. Our finding that males exhibit much more prominent alterations in these pathways with age is consistent with the clinical finding that older men have a weaker response to vaccination and are more likely to develop infections or cancer. One hypothesis based on this study would be that the Phase II SHM process is less efficient in older males, resulting in a less optimal Ab repertoire in response to Ags.

Aside from aging, other factors may also affect SHM targeting, such as genetic factors and health state. Despite the fact that many factors can impact one’s immunity, we found significant correlations between A/T mutability and age in males. It is interesting that the frail male showed the largest decrease in the major SHM repair molecules, which might suggest that a lack of SHM repair activities is associated with weakened immune response. Future studies should examine the role of sex-specific molecules in the SHM process. For example, a study revealed the immunosuppressive function of testosterone in response to influenza vaccination, indicating that sex hormones could play a role (32).

It is interesting to note that many DNA repair enzymes and polymerases are involved in both Ab diversification and DNA damage repair in other parts of the genome. A previous study observed a correlation between microsatellite stability, particularly PMS2 expression, and types of mutations in B cell repertoires (33). In this study, we observed a slight reduction of PMS2 expression in older males, which may reflect their conclusion on the influence of MMR activity on the types of mutations generated through SHM. Future studies should examine the pleiotropic effects of DNA repair molecules. For example, individuals with deficiency in DNA repair genes may not only have an impaired Ab repertoire, but also a higher predisposition to diseases caused by DNA damages, such as cancer.

A limitation of the gene expression dataset is that the B cells were collected from peripheral blood instead of the GC where SHM takes place. Because of this, we found that some genes encoding important molecules in SHM, such as AID and pol-η, were not highly expressed, limiting our ability to study all genes involved in SHM. Additionally, other machineries, aside from gene expression of SHM-related molecules, may also influence SHM targeting. This prompts future studies to examine B cells in GCs to see what other DNA repair pathways contribute to the aging effects that we observed. A limitation of the GTEx RNA-seq dataset is the small number of B cell Ig genes one can extract from transcriptome-wide RNA-seq data on all cell types in a tissue. The AIRR-seq dataset has many more mutations per individual, although it has fewer subjects. However, the fact that both AIRR-seq and GTEx RNA-seq analyses yielded consistent findings improves our confidence in the sex-biased aging effects in SHM targeting. In the future, larger cohorts of individuals should be recruited to characterize the extent to which age and sex can influence SHM-targeting activities and which genes experience the most alteration in expression levels with age. Nevertheless, this study showed that age and sex could influence the balance between the two phases of SHM targeting, which is crucial to Ab diversification.

In conclusion, our analyses on BCR sequencing and B cell gene expression data showed alterations in SHM mutation patterns and gene expressions between age and sex groups. These results suggest an uneven rate of alteration in SHM DNA repair machineries between males and females across lifetime, providing insights into designing more targeted treatment and vaccination strategies for different age and sex groups.

Disclosures

S.H.K. receives consulting fees from Northrop Grumman.

Acknowledgments

We thank the Kleinstein laboratory members for feedback and technical assistance and Dr. Scott Boyd for sharing data. We thank the Yale University Biomedical High Performance Computing Center and the Broad Institute Computing Center for use of the computing resources.

Footnotes

  • This work was supported by grants from the National Institutes of Health (R03AI092379 and R01AI104739) to S.H.K. and a Natural Sciences and Engineering Research Council of Canada postgraduate fellowship (NSERC PGS-M) to A.C. The Yale University Biomedical High Performance Computing Center is funded by National Institutes of Health Grants RR19895 and RR029676-01. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

  • The online version of this article contains supplemental material.

  • Abbreviations used in this article:

    AID
    activation-induced cytidine deaminase
    AIRR-seq
    adaptive immune receptor repertoire sequencing
    BER
    base-excision repair
    GC
    germinal center
    GTEx
    Genotype-Tissue Expression
    MMR
    mismatch repair
    RNA-seq
    RNA-sequencing
    SHM
    somatic hypermutation.

  • Received May 22, 2020.
  • Accepted October 25, 2020.

References

  1. ↵
    1. Weinberger, B.,
    2. D. Herndler-Brandstetter,
    3. A. Schwanninger,
    4. D. Weiskopf,
    5. B. Grubeck-Loebenstein
    . 2008. Biology of immune responses to vaccines in elderly persons. Clin. Infect. Dis. 46: 1078–1084.
    OpenUrlCrossRefPubMed
    1. Castle, S. C.
    2000. Clinical relevance of age-related immune dysfunction. Clin. Infect. Dis. 31: 578–585.
    OpenUrlCrossRefPubMed
  3. ↵
    1. Castelo-Branco, C.,
    2. I. Soveral
    . 2014. The immune system and aging: a review. Gynecol. Endocrinol. 30: 16–22.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Márquez, E. J.,
    2. C.-H. Chung,
    3. R. Marches,
    4. R. J. Rossi,
    5. D. Nehar-Belaid,
    6. A. Eroglu,
    7. D. J. Mellert,
    8. G. A. Kuchel,
    9. J. Banchereau,
    10. D. Ucar
    . 2020. Sexual-dimorphism in human immune system aging. Nat. Commun. 11: 751.
    OpenUrlCrossRef
  5. ↵
    1. Beeson, P. B.
    1994. Age and sex associations of 40 autoimmune diseases. Am. J. Med. 96: 457–462.
    OpenUrlCrossRefPubMed
    1. Klein, S. L.,
    2. I. Marriott,
    3. E. N. Fish
    . 2015. Sex-based differences in immune function and responses to vaccination. Trans. R. Soc. Trop. Med. Hyg. 109: 9–15.
    OpenUrlCrossRefPubMed
    1. Fish, E. N.
    2008. The X-files in immunity: sex-based differences predispose immune responses. Nat. Rev. Immunol. 8: 737–744.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Gubbels Bupp, M. R.
    2015. Sex, the aging immune system, and chronic disease. Cell. Immunol. 294: 102–110.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Dunn-Walters, D. K.
    2016. The ageing human B cell repertoire: a failure of selection? Clin. Exp. Immunol. 183: 50–56.
    OpenUrl
    1. Dunn-Walters, D. K.,
    2. M. Banerjee,
    3. R. Mehr
    . 2003. Effects of age on antibody affinity maturation. Biochem. Soc. Trans. 31: 447–448.
    OpenUrlAbstract/FREE Full Text
    1. Weiskopf, D.,
    2. B. Weinberger,
    3. B. Grubeck-Loebenstein
    . 2009. The aging of the immune system. Transpl. Int. 22: 1041–1050.
    OpenUrlCrossRefPubMed
  12. ↵
    1. LeMaoult, J.,
    2. P. Szabo,
    3. M. E. Weksler
    . 1997. Effect of age on humoral immunity, selection of the B-cell repertoire and B-cell development. Immunol. Rev. 160: 115–126.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Peled, J. U.,
    2. F. L. Kuang,
    3. M. D. Iglesias-Ussel,
    4. S. Roa,
    5. S. L. Kalis,
    6. M. F. Goodman,
    7. M. D. Scharff
    . 2008. The biochemistry of somatic hypermutation. Annu. Rev. Immunol. 26: 481–511.
    OpenUrlCrossRefPubMed
  14. ↵
    1. McKean, D.,
    2. K. Huppi,
    3. M. Bell,
    4. L. Staudt,
    5. W. Gerhard,
    6. M. Weigert
    . 1984. Generation of antibody diversity in the immune response of BALB/c mice to influenza virus hemagglutinin. Proc. Natl. Acad. Sci. USA 81: 3180–3184.
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Rogozin, I. B.,
    2. N. A. Kolchanov
    . 1992. Somatic hypermutagenesis in immunoglobulin genes. II. Influence of neighbouring base sequences on mutagenesis. Biochim. Biophys. Acta 1171: 11–18.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Frasca, D.,
    2. A. M. Landin,
    3. S. C. Lechner,
    4. J. G. Ryan,
    5. R. Schwartz,
    6. R. L. Riley,
    7. B. B. Blomberg
    . 2008. Aging down-regulates the transcription factor E2A, activation-induced cytidine deaminase, and Ig class switch in human B cells. J. Immunol. 180: 5283–5290.
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Benichou, J.,
    2. R. Ben-Hamo,
    3. Y. Louzoun,
    4. S. Efroni
    . 2012. Rep-Seq: uncovering the immunological repertoire through next-generation sequencing. Immunology 135: 183–191.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Rubelt, F.,
    2. C. E. Busse,
    3. S. A. C. Bukhari,
    4. J. P. Bürckert,
    5. E. Mariotti-Ferrandiz,
    6. L. G. Cowell,
    7. C. T. Watson,
    8. N. Marthandan,
    9. W. J. Faison,
    10. U. Hershberg, et al, AIRR Community
    . 2017. Adaptive immune receptor repertoire community recommendations for sharing immune-repertoire sequencing data. Nat. Immunol. 18: 1274–1278.
    OpenUrlCrossRefPubMed
  19. ↵
    1. Wang, C.,
    2. Y. Liu,
    3. L. T. Xu,
    4. K. J. L. Jackson,
    5. K. M. Roskin,
    6. T. D. Pham,
    7. J. Laserson,
    8. E. L. Marshall,
    9. K. Seo,
    10. J.-Y. Lee, et al
    . 2014. Effects of aging, cytomegalovirus infection, and EBV infection on human B cell repertoires. J. Immunol. 192: 603–611.
    OpenUrlAbstract/FREE Full Text
    1. Jiang, N.,
    2. J. He,
    3. J. A. Weinstein,
    4. L. Penland,
    5. S. Sasaki,
    6. X. S. He,
    7. C. L. Dekker,
    8. N.-Y. Zheng,
    9. M. Huang,
    10. M. Sullivan, et al
    . 2013. Lineage structure of the human antibody repertoire in response to influenza vaccination. Sci. Transl. Med. 5: 171ra19.
  21. ↵
    1. de Bourcy, C. F. A.,
    2. C. J. L. Angel,
    3. C. Vollmers,
    4. C. L. Dekker,
    5. M. M. Davis,
    6. S. R. Quake
    . 2017. Phylogenetic analysis of the human antibody repertoire reveals quantitative signatures of immune senescence and aging. Proc. Natl. Acad. Sci. USA 114: 1105–1110.
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Yaari, G.,
    2. J. A. Vander Heiden,
    3. M. Uduman,
    4. D. Gadala-Maria,
    5. N. Gupta,
    6. J. N. Stern,
    7. K. C. O’Connor,
    8. D. A. Hafler,
    9. U. Laserson,
    10. F. Vigneault,
    11. S. H. Kleinstein
    . 2013. Models of somatic hypermutation targeting and substitution based on synonymous mutations from high-throughput immunoglobulin sequencing data. Front. Immunol. 4: 358.
    OpenUrlCrossRefPubMed
  23. ↵
    1. Cui, A.,
    2. R. Di Niro,
    3. J. A. Vander Heiden,
    4. A. W. Briggs,
    5. K. Adams,
    6. T. Gilbert,
    7. K. C. O’Connor,
    8. F. Vigneault,
    9. M. J. Shlomchik,
    10. S. H. Kleinstein
    . 2016. A model of somatic hypermutation targeting in mice based on high-throughput Ig sequencing data. J. Immunol. 197: 3566–3574.
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. GTEx Consortium
    . 2013. The genotype-tissue expression (GTEx) project. Nat. Genet. 45: 580–585.
    OpenUrlCrossRefPubMed
  25. ↵
    1. Vander Heiden, J. A.,
    2. G. Yaari,
    3. M. Uduman,
    4. J. N. H. Stern,
    5. K. C. O’Connor,
    6. D. A. Hafler,
    7. F. Vigneault,
    8. S. H. Kleinstein
    . 2014. pRESTO: a toolkit for processing high-throughput sequencing raw reads of lymphocyte receptor repertoires. Bioinformatics 30: 1930–1932.
    OpenUrlCrossRefPubMed
  26. ↵
    1. Alamyar, E.,
    2. P. Duroux,
    3. M.-P. Lefranc,
    4. V. Giudicelli
    . 2012. IMGT(®) tools for the nucleotide analysis of immunoglobulin (IG) and T cell receptor (TR) V-(D)-J repertoires, polymorphisms, and IG mutations: IMGT/V-QUEST and IMGT/HighV-QUEST for NGS. Methods Mol. Biol. 882: 569–604.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Gupta, N. T.,
    2. J. A. Vander Heiden,
    3. M. Uduman,
    4. D. Gadala-Maria,
    5. G. Yaari,
    6. S. H. Kleinstein
    . 2015. Change-O: a toolkit for analyzing large-scale B cell immunoglobulin repertoire sequencing data. Bioinformatics 31: 3356–3358.
    OpenUrlCrossRefPubMed
  28. ↵
    1. Smith, D. S.,
    2. G. Creadon,
    3. P. K. Jena,
    4. J. P. Portanova,
    5. B. L. Kotzin,
    6. L. J. Wysocki
    . 1996. Di- and trinucleotide target preferences of somatic mutagenesis in normal and autoreactive B cells. J. Immunol. 156: 2642–2652.
    OpenUrlAbstract
  29. ↵
    1. Ye, J.,
    2. N. Ma,
    3. T. L. Madden,
    4. J. M. Ostell
    . 2013. IgBLAST: an immunoglobulin variable domain sequence analysis tool. Nucleic Acids Res. 41(W1): W34–W40.
    OpenUrlCrossRefPubMed
  30. ↵
    1. Avey, S.,
    2. S. Mohanty,
    3. D. G. Chawla,
    4. H. Meng,
    5. T. Bandaranayake,
    6. I. Ueda,
    7. H. J. Zapata,
    8. K. Park,
    9. T. P. Blevins,
    10. S. Tsang, et al
    . 2020. Seasonal variability and shared molecular signatures of inactivated influenza vaccination in young and older adults. J. Immunol. 204: 1661–1673.
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. Hoehn, K. B.,
    2. J. A. Vander Heiden,
    3. J. Q. Zhou,
    4. G. Lunter,
    5. O. G. Pybus,
    6. S. H. Kleinstein
    . 2019. Repertoire-wide phylogenetic models of B cell molecular evolution reveal evolutionary signatures of aging and vaccination. Proc. Natl. Acad. Sci. USA 116: 22664–22672.
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Furman, D.,
    2. B. P. Hejblum,
    3. N. Simon,
    4. V. Jojic,
    5. C. L. Dekker,
    6. R. Thiébaut,
    7. R. J. Tibshirani,
    8. M. M. Davis
    . 2014. Systems analysis of sex differences reveals an immunosuppressive role for testosterone in the response to influenza vaccination. Proc. Natl. Acad. Sci. USA 111: 869–874.
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Rosner, K.,
    2. D. B. Winter,
    3. C. Kasmer,
    4. G. L. Skovgaard,
    5. R. E. Tarone,
    6. V. A. Bohr,
    7. P. J. Gearhart
    . 2001. Impact of age on hypermutation of immunoglobulin variable genes in humans. J. Clin. Immunol. 21: 102–115.
    OpenUrlCrossRefPubMed
Back to top

In this issue

Print
Download PDF
Article Alerts
Email Article
Citation Tools

Jump to section