Department of Biochemistry and Cell Biology, Stony Brook University, Stony Brook, NY, United States.
Department of Applied Mathematics and Statistics, Stony Brook University, Stony Brook, NY, United States.
Front Immunol. 2024 May 28;15:1407470. doi: 10.3389/fimmu.2024.1407470. eCollection 2024.
INTRODUCTION: Somatic hypermutation (SHM) of immunoglobulin variable (V) regions by activation induced deaminase (AID) is essential for robust, long-term humoral immunity against pathogen and vaccine antigens. AID mutates cytosines preferentially within WRCH motifs (where W=A or T, R=A or G and H=A, C or T). However, it has been consistently observed that the mutability of WRCH motifs varies substantially, with large variations in mutation frequency even between multiple occurrences of the same motif within a single V region. This has led to the notion that the immediate sequence context of WRCH motifs contributes to mutability. Recent studies have highlighted the potential role of local DNA sequence features in promoting mutagenesis of AGCT, a commonly mutated WRCH motif. Intriguingly, AGCT motifs closer to 5' ends of V regions, within the framework 1 (FW1) sub-region1, mutate less frequently, suggesting an SHM-suppressing sequence context. METHODS: Here, we systematically examined the basis of AGCT positional biases in human SHM datasets with DeepSHM, a machine-learning model designed to predict SHM patterns. This was combined with integrated gradients, an interpretability method, to interrogate the basis of DeepSHM predictions. RESULTS: DeepSHM predicted the observed positional differences in mutation frequencies at AGCT motifs with high accuracy. For the conserved, lowly mutating AGCT motifs in FW1, integrated gradients predicted a large negative contribution of 5'C and 3'G flanking residues, suggesting that a CAGCTG context in this location was suppressive for SHM. CAGCTG is the recognition motif for E-box transcription factors, including E2A, which has been implicated in SHM. Indeed, we found a strong, inverse relationship between E-box motif fidelity and mutation frequency. Moreover, E2A was found to associate with the V region locale in two human B cell lines. Finally, analysis of human SHM datasets revealed that naturally occurring mutations in the 3'G flanking residues, which effectively ablate the E-box motif, were associated with a significantly increased rate of AGCT mutation. DISCUSSION: Our results suggest an antagonistic relationship between mutation frequency and the binding of E-box factors like E2A at specific AGCT motif contexts and, therefore, highlight a new, suppressive mechanism regulating local SHM patterns in human V regions.
简介:通过激活诱导脱氨酶(AID)对免疫球蛋白可变(V)区进行体细胞超突变(SHM)是针对病原体和疫苗抗原产生强大、长期体液免疫的关键。AID 优先在 WRCH 基序(其中 W=A 或 T,R=A 或 G,H=A、C 或 T)内突变胞嘧啶。然而,人们一直观察到 WRCH 基序的突变率存在很大差异,即使在单个 V 区的同一基序的多个发生中,突变频率也有很大差异。这导致了这样一种观点,即 WRCH 基序的直接序列上下文有助于突变性。最近的研究强调了局部 DNA 序列特征在促进 AGCT 突变中的潜在作用,AGCT 是一个常见的突变 WRCH 基序。有趣的是,更靠近 V 区 5'端的 AGCT 基序,在框架 1(FW1)亚区 1 内,突变频率较低,表明存在 SHM 抑制序列上下文。
方法:在这里,我们使用深度学习 SHM(DeepSHM),一种旨在预测 SHM 模式的机器学习模型,系统地检查了人类 SHM 数据集中原位 AGCT 偏倚的基础。这与积分梯度(一种解释性方法)相结合,以探究 DeepSHM 预测的基础。
结果:DeepSHM 以高精度预测了 AGCT 基序中观察到的突变频率的位置差异。对于 FW1 中保守、突变率低的 AGCT 基序,积分梯度预测 5'C 和 3'G 侧翼残基的贡献很大,表明该位置的 CAGCTG 上下文对 SHM 具有抑制作用。CAGCTG 是 E 盒转录因子(包括 E2A)的识别基序,E2A 已被牵连到 SHM 中。事实上,我们发现 E 盒基序保真度与突变频率之间存在强烈的反比关系。此外,在两个人类 B 细胞系中发现 E2A 与 V 区位置相关。最后,对人类 SHM 数据集的分析表明,3'G 侧翼残基的自然突变有效地消除了 E 盒基序,与 AGCT 突变率的显著增加有关。
讨论:我们的结果表明,在特定的 AGCT 基序背景下,突变频率与 E 盒因子(如 E2A)的结合之间存在拮抗关系,因此突出了一种新的、抑制性机制,调节了人类 V 区的局部 SHM 模式。
Immunology. 2010-11-11
Annu Rev Immunol. 2022-4-26
Nucleic Acids Res. 2022-1-7
Proc Natl Acad Sci U S A. 2021-7-20
Gigascience. 2021-2-16
J Immunol. 2020-12-15