Literature DB >> 32434680

AID in Antibody Diversification: There and Back Again.

Yuqing Feng1, Noé Seija2, Javier M Di Noia3, Alberto Martin4.   

Abstract

Activation-Induced cytidine Deaminase (AID) initiates affinity maturation and isotype switching by deaminating deoxycytidines within immunoglobulin genes, leading to somatic hypermutation (SHM) and class switch recombination (CSR). AID thus potentiates the humoral response to clear pathogens. Marking the 20th anniversary of the discovery of AID, we review the current understanding of AID function. We discuss AID biochemistry and how error-free forms of DNA repair are co-opted to prioritize mutagenesis over accuracy during antibody diversification. We discuss the regulation of DNA double-strand break (DSB) repair pathways during CSR. We describe genomic targeting of AID as a multilayered process involving chromatin architecture, cis- and trans-acting factors, and determining mutagenesis - distinct from AID occupancy at loci that are spared from mutation.
Copyright © 2020 Elsevier Ltd. All rights reserved.

Entities:  

Year:  2020        PMID: 32434680      PMCID: PMC7183997          DOI: 10.1016/j.it.2020.04.009

Source DB:  PubMed          Journal:  Trends Immunol        ISSN: 1471-4906            Impact factor:   16.687


AID as a Central Factor to Generating Antibody Repertoires

Antibodies are essential to clearing viral, bacterial, and parasitic infections, as well as neutralizing toxins. The ability of Activation-Induced cytidine Deaminase (AID) to enhance antibody affinity through somatic hypermutation (SHM) (see Glossary) and alter the effector function of antibodies through class switch recombination (CSR) is vital to these tasks (see later). In addition, the generation of a long-lasting antigen-specific antibody memory response is crucial to establishing immunity against pathogens, and AID plays a key role in this endeavor by engraving memories of past pathogen experience through genetic modifications of immunoglobulin (Ig) genes. Here, we review the most essential findings on AID since its discovery two decades ago by providing a timely update on the molecular mechanisms of AID in B cell antibody production. As this review was written during the period that the world was engulfed by the COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-COV-2), we hope it will provide background to inspire vaccine design and the generation of therapeutic antibodies against SARS-CoV-2, as well as inevitable future pathogens.

General Features of Secondary Antibody Diversification Mechanisms

Vertebrates have evolved elegant genetic mechanisms to produce a repertoire of antibodies capable of recognizing and neutralizing a potentially unlimited number of foreign pathogens. Prior to antigen exposure, site-specific recombination of germline Ig variable (V), diversity (D), and joining (J) segments during B cell development provides the primary antibody repertoire [1]. The processes of SHM, gene conversion (GCV), and CSR further diversify this antibody repertoire upon B cell-cognate antigen interactions. In antigen-activated B cells, SHM introduces point mutations in the Ig V region of both Ig heavy (Igh) and light chains to enable the process of affinity maturation in the germinal center of secondary lymphoid organs. GCV is an alternative mechanism of altering antibody specificity and affinity in select vertebrate species, including rabbits and birds. GCV involves the production of single-stranded (ss) DNA breaks or DSBs within the rearranged V(D)J region, triggering a homologous recombination-like process that repairs the lesion using Ig V pseudogenes. This produces templated sequence alterations within the V region. CSR is a region-specific deletional recombination event whereby the exons encoding for the IgM constant region within the Igh are replaced with a downstream constant region without changing the antigen-binding specificity, thereby altering the effector functions of the antibody. While all three processes are distinct genetic transactions, they share some common features, such as transcription of their target sequence [2], and preferential mutation of RGYW/WRCY motifs (where W=A/T, R=A/G, Y=C/T), wherein the underlined G/C is mutated [3]. While not demonstrated until the discovery of AID, these and other specific similarities suggested a common inception to each of these processes.

Discovery of AID and Its Molecular Mechanism in Antibody Diversification

AID was identified by Tasuku Honjo’s group in 1999 [4], and shown to be essential for SHM and CSR in both mice and humans in 2000 [5,6], and later for GCV in chicken DT40 B lymphoma cell line [7]. Variable lymphocyte receptor diversification in ancient jawless vertebrates also depends on an AID-like cytidine deaminases named CDA1 and CDA2 [8], indicating that diversification of antigen receptors by deamination is an ancestral mechanism in vertebrates. AID was originally proposed to be an RNA editing enzyme due to sequence homology to its paralog, the RNA-cytosine deaminase apolipoprotein B mRNA editing catalytic polypeptide 1 (APOBEC 1) [4]. However, the DNA deamination model which proposes that AID promotes antibody diversity by deaminating deoxycytidine (dC) to deoxyuridine (dU) within Ig genes is currently supported by significant experimental evidence [9]. Purified AID preferentially mutates dC within the WRCY motif [10,11]. As this motif is preferentially mutated during SHM in vivo [12], this result provides strong evidence supporting the DNA deamination model for AID function. Biochemical data suggest that AID is catalytically inefficient [13,14]. The low catalytic rate has been proposed to act as an evolutionary safeguard to limit the mutagenic activity of AID [11]. Moreover, the accessibility of AID’s catalytic pocket appears to be a key determinant of AID’s catalytic rate [15]. AID has been suggested to act on ssDNA by inducing multiple deaminations per binding event [10]. A single-molecule resolution Förster resonance energy transfer (FRET) study visualized the ssDNA scanning motion of AID, and suggested that the enzyme was able to move bidirectionally while remaining bound to the same ssDNA [16]. This behavior of AID on ssDNA may play a role in enhancing mutagenic repair by overwhelming faithful DNA repair (see later). Modeling AID on pre-existing structures of APOBEC enzymes [15] and crystallizing truncated AID variants [17,18] have provided valuable structural details. X-ray structure analysis with human monomeric AID reveals that AID prefers to bind and deaminate structured substrates, such as G-quadruplex (G4) structures that can form at the switch (S) region, over linear ssDNA substrates [18., 19., 20.] (Box 1 ). However, as such structures are rare, especially over diverse Ig V exons, it is likely that other structures are also acted upon by AID, such as short patches of melted DNA [21]. These X-ray structure studies suggested that AID contains a bifurcated nucleic acid binding surface [15,18]: one required for catalysis, and another positively charged surface in AID helix 6, thought to play a role in recognizing structured substrates. A similar bifurcated nucleic acid binding surface has been identified in APOBEC3H, mediating dimerization [22]. Besides facilitating substrate binding, specific arginine residues in AID helix 6 may underpin a licensing mechanism that links AID to transcription elongation and warrants productive targeting [23]. Nevertheless, biochemical, biophysical and structural approaches have provided valuable knowledge on AID targeting mechanisms during SHM and CSR (see later). A detailed summary of the structure–function description of AID can be found in [9]. The transcriptional units encoding the different constant domains of Igh have revealed roles for RNA structures and RNA processing in CSR and AID targeting. These elements are transcribed from a cytokine-driven promoter, which produces a sterile transcript to allow AID access. The primary transcript contains a large intron encompassing the S region (1–12 kb), made up of repeats rich in WRCY motifs. The G-rich composition of the S regions enables the formation of secondary structures such as DNA–RNA hybrids producing R loops and G4 quadruplexes during transcription [20]. One model envisages that cotranscriptional formation of DNA G4 on the nontemplate strand ssDNA stabilizes R loops on the template strand, while processing of these R loops by the RNA exosome exposes ssDNA and permits deamination of both strands [20]. Others have proposed a role for R loops in CSR resolution, rather than on attracting AID activity [20]. AID shows high activity on G4 or branched DNA substrates in biochemical assays, allegedly due to a bifurcated substrate-binding surface that can recognize ssDNA overhangs on these structures [18]. G4 DNA additionally promotes AID oligomerization in vitro, which would facilitate DSB formation by increasing the deamination density at the S regions [18]; however, this remains to be tested in vivo. Another model proposes that AID binds to G4-RNA formed in the debranched S-region intron, which would then guide AID back to the DNA [104], presumably by base complementarity. This model could account for the link of CSR to RNA splicing [110] and the interaction of AID with splicing factors such as heterogeneous nuclear ribonucleoproteins or polypyrimidine tract-binding proteins [9]. However, the mechanism by which G4-RNA guides and hands off AID to the DNA is not clear. The unwinding of G4-RNA in S-region transcripts by the DEAD-box RNA helicase 1 (DDX1) to allow the formation of DNA:RNA R-loop hybrids at the Igh by base complementarity was recently proposed as a mechanism [105]. However, ablation of DDX1 reduced CSR only by 25% when analyzed per cell division, indicating that this mechanism is not essential for CSR. Additional aspects of this model remain to be demonstrated, including evidence of AID transfer from RNA to DNA in vivo, and whether this transfer participates in targeting AID for SHM at the Ig V or at off-target loci. Alt-text: Box 1

Co-opted DNA Repair Pathways Involved in SHM

Uracils introduced in DNA are normally repaired faithfully, restoring the original DNA sequence (Figure 1 ). This is not the case during SHM, GCV, or CSR, whereby the initiating AID-induced lesion promotes further mutagenesis that benefits affinity maturation by providing a wide range of amino acid changes, or causes DNA breaks that trigger recombination events. Why faithful repair of DNA uracils is subverted during SHM has been the subject of intense investigation. The majority of the uracils in the genome are faithfully repaired by the base excision repair (BER) pathway: excision of uracils by uracil DNA glycosylase (UNG); incision at the resulting abasic sites by apurinic/apyrimidinic endonuclease 1; processing of the abasic sites by DNA polymerase β, followed by sealing of the final nick by DNA ligase 3 [24]. Uracils escaping detection by UNG or single-strand-selective monofunctional uracil-DNA glycosylase 1 (SMUG1) can be copied over by the replication machinery, producing transition mutations at G:C base pairs (Figure 1) [25,26]. However, during SHM, the later steps of BER are circumvented, for reasons that are not fully clear (see Outstanding Questions). Replication across the noninstructive abasic site by the translesional DNA polymerases REV1 and Pol η generates both transition and transversion mutations at G:C pairs [27,28] (Figure 1). Indeed, mice deficient in UNG (Ung -/-) exhibit reduced transversion mutations at G:C pairs in Ig V regions relative to wild-type (WT) mice, with a minor impact on mutagenesis at A:T base pairs [29].
Figure 1

Multifaceted Uracil Lesion Processing in Somatic Hypermutation (SHM).

AID-Induced dU:dG mismatches can be resolved by various pathways, leading to a diverse mutation spectrum at G:C and A:T base pairs in immunoglobulin variable (Ig V) regions. AID-induced uracil lesions can be repaired faithfully by canonical BER and MMR, leading to no SHM. Replicating over uracils by replicative polymerases generates transition mutations at G:C base pairs. UNG, and to a lesser extent, SMUG1, from the BER pathway, detect and excise uracils from DNA, creating abasic sites. Replicating over abasic sites by translesion synthesis polymerase from short-patch ncBER pathway generates transition and transversion mutations at G:C base pairs. Noncanonical MMR pathway or long-patch BER can excise single-stranded DNA surrounding either the dU:dG mismatch, or the abasic site, respectively, to initiate patch excision that leads to the engagement of translesion synthesis polymerases, generating both transition and transversion mutations at A:T base pairs. Abbreviations: AID, Activation-induced Cytidine Deaminase; ncBER, noncanonical base excision repair; ncMMR, noncanonical mismatch repair; SMUG1, single-strand selective monofunctional uracil DNA glycosylase; UNG, uracil DNA glycosylase. This figure was created using BioRender (https://biorender.com/).

Multifaceted Uracil Lesion Processing in Somatic Hypermutation (SHM). AID-Induced dU:dG mismatches can be resolved by various pathways, leading to a diverse mutation spectrum at G:C and A:T base pairs in immunoglobulin variable (Ig V) regions. AID-induced uracil lesions can be repaired faithfully by canonical BER and MMR, leading to no SHM. Replicating over uracils by replicative polymerases generates transition mutations at G:C base pairs. UNG, and to a lesser extent, SMUG1, from the BER pathway, detect and excise uracils from DNA, creating abasic sites. Replicating over abasic sites by translesion synthesis polymerase from short-patch ncBER pathway generates transition and transversion mutations at G:C base pairs. Noncanonical MMR pathway or long-patch BER can excise single-stranded DNA surrounding either the dU:dG mismatch, or the abasic site, respectively, to initiate patch excision that leads to the engagement of translesion synthesis polymerases, generating both transition and transversion mutations at A:T base pairs. Abbreviations: AID, Activation-induced Cytidine Deaminase; ncBER, noncanonical base excision repair; ncMMR, noncanonical mismatch repair; SMUG1, single-strand selective monofunctional uracil DNA glycosylase; UNG, uracil DNA glycosylase. This figure was created using BioRender (https://biorender.com/). Mismatch repair (MMR) also mediates repair of uracil lesions, detected as U:G mispairs. Replication-associated canonical MMR (cMMR) is traditionally a high-fidelity process that preserves genomic stability. cMMR begins with the recognition of a base mispair by the MutSα heterodimer complex (MSH2/MSH6), followed by the recruitment of MutLα complex (MLH1/PMS2) and the activation of the endonuclease activity of PMS2, and removal of a track of the DNA strand by exonuclease EXO1 [30,31]. The resynthesis of the ssDNA gap by high-fidelity replicative polymerases restores the original DNA sequences. Instead, during SHM, a subverted MMR pathway termed noncanonical MMR (ncMMR) using the translesion DNA polymerase η is largely responsible for mutations at A:T base pairs [31., 32., 33., 34., 35., 36., 37.] (Figure 1). Mutagenic ncMMR is largely linked to AID lesions, as non-AID targets are faithfully repaired by MMR [38]. Recent work by a number of groups has shed light on ncMMR, whereby the UNG/SMUG1 and MMR pathways converge and work collaboratively to facilitate A:T and G:C mutagenesis [39., 40., 41., 42., 43.]. This view is supported by findings that both UNG and ncMMR are able to function in the G1 phase of the cell cycle [41,42], where AID activity peaks [44]. Hence, the introduction of an abasic site within the EXO1 excision track induces translesional repair and mutations at A:T base pairs. This can arise in a number of ways: EXO1 excises the wrong strand of an AID-induced G:U mismatch, thereby exposing dU on the ssDNA track to uracil excision by UNG or SMUG1. Alternatively, AID disrupts cMMR by deaminating a dC within the EXO1 excision track, triggering UNG/SMUG1 to create an abasic site that leads to mutagenic translesional repair.

AID-Initiated DSB Formation and Repair during CSR

DSB formation in the Ig S region is a prerequisite for CSR. Ig S regions are highly enriched with AGCT motifs, a palindromic variant of the WRCY motif (Figure 2A) [45]. Similar to SHM, the generation of DSBs during CSR also co-opts the activities of UNG and MMR, as Ung -/- Msh2 -/- mice show abolished CSR relative to WT mice [46]. MMR- and UNG-induced nicks and gaps on opposite DNA strands are converted into DSBs in donor and acceptor S regions, which are synapsed and ligated by DSB repair mechanisms, completing the CSR process. DSB repair during CSR primarily relies on canonical nonhomologous end joining (NHEJ), which promotes DSB end joining with no or minimal sequence microhomology. Alternative end joining (Alt-EJ) plays a minor role in CSR and utilizes short regions of microhomology between both DNA ends to join DSBs [47] (Figure 2B). DSB repair by NHEJ involves limited DNA end processing, whereas increased resection and microhomology usage are two features of repair by Alt-EJ. As such, mechanisms that regulate DNA end resection promote usage of one DSB repair pathway over the other (see later). The structure of the AID-induced DSB has also been suggested to provide a decision of which DSB repair pathway to use [48,49]. The stochastic activity of AID on both strands of DNA leads to a mixed spectrum of DSB structures in the S region. Accordingly, staggered DSBs are the preferred substrates for Alt-EJ, whereas blunt DSBs are preferentially repaired by the NHEJ pathway [49]. Indeed, CSR serves as an attractive physiological model for studying various DNA repair pathways at the same time, as all of BER, MMR, NHEJ, and Alt-EJ play key roles in DSB formation and repair.
Figure 2

AID-Mediated DSB Formation within S Regions Drives CSR.

(A) Sequence feature of S regions in mouse Ig heavy chain locus: S regions preceding different constant region exons are highly enriched with AID deamination hotspot 5′-AGCT-3′, a palindromic variant representing the canonical RGYW/WRCY sequences. Only Sμ and Sα are shown here for simplicity. (B) DSB ends are joined through NHEJ and Alt-EJ to complete CSR. NHEJ mediates DSB EJ with no or minimal microhomology, whereas Alt-EJ is facilitated by complementary base-paring interactions between microhomologous sequences present between the two synapsed S regions DSBs. Black bars represent homologous sequences near DSB ends. Abbreviations: 53BP1, p53-binding protein 1; AID, Activation-induced Cytidine Deaminase; Alt-EJ, alternative end joining; CSR, class switch recombination; CTIP, CtBP-interacting protein; DNA-PKcs, DNA-dependent protein kinase, catalytic subunit; DSB, double-strand break; ES cell specific; HMCES, 5-Hydroxymethylcytosine binding, ES-cell-specific; KU, Ku 70/80 complex; LIG3, ligase 3; LIG4, ligase 4; MRN, MRE11–RAD50–NBN complex; NHEJ, nonhomologous end joining; RIF1, replication timing regulatory factor 1; S, switch; SHIELDIN, protein complex composed of SHLD1, SHLD2, SHLD3, and REV 7; XLF, XRCC4-like factor; XRCC 1/4, X-ray repair cross-complementing protein 1/4. This figure was created using BioRender (https://biorender.com/).

AID-Mediated DSB Formation within S Regions Drives CSR. (A) Sequence feature of S regions in mouse Ig heavy chain locus: S regions preceding different constant region exons are highly enriched with AID deamination hotspot 5′-AGCT-3′, a palindromic variant representing the canonical RGYW/WRCY sequences. Only Sμ and Sα are shown here for simplicity. (B) DSB ends are joined through NHEJ and Alt-EJ to complete CSR. NHEJ mediates DSB EJ with no or minimal microhomology, whereas Alt-EJ is facilitated by complementary base-paring interactions between microhomologous sequences present between the two synapsed S regions DSBs. Black bars represent homologous sequences near DSB ends. Abbreviations: 53BP1, p53-binding protein 1; AID, Activation-induced Cytidine Deaminase; Alt-EJ, alternative end joining; CSR, class switch recombination; CTIP, CtBP-interacting protein; DNA-PKcs, DNA-dependent protein kinase, catalytic subunit; DSB, double-strand break; ES cell specific; HMCES, 5-Hydroxymethylcytosine binding, ES-cell-specific; KU, Ku 70/80 complex; LIG3, ligase 3; LIG4, ligase 4; MRN, MRE11–RAD50–NBN complex; NHEJ, nonhomologous end joining; RIF1, replication timing regulatory factor 1; S, switch; SHIELDIN, protein complex composed of SHLD1, SHLD2, SHLD3, and REV 7; XLF, XRCC4-like factor; XRCC 1/4, X-ray repair cross-complementing protein 1/4. This figure was created using BioRender (https://biorender.com/).

NHEJ Cascade

The Ku70/80 complex, XRCC4, and Ligase 4 (Lig4) represent the core NHEJ factors based on their exclusive role in end joining. Upon AID-induced S-region DSB formation, the Mre11–Rad50–Nbs1 (MRN) sensor–-kinase complex binds DSBs and recruits the ataxia telangiectasia-mutated (ATM) kinase to the break site. ATM phosphorylates numerous substrates, leading to cell cycle arrest, as well as phosphorylating histone H2AX, which stimulates a cascade of events leading to the recruitment of 53BP1 to chromatin flanking the DSB [50,51]. 53BP1 promotes NHEJ and CSR by preventing DSB end resection [52] through the recruitment of Rif1 [53., 54., 55., 56.] and Shieldin, a protein complex composed of Shld1, Shld2, Shld3, and Rev7 [57., 58., 59., 60., 61., 62.] (Figure 2B). How the Shieldin complex prevents DNA end resection is still not clear, but it requires the binding of Shld2 to ssDNA exposed at DSBs [63] (see Outstanding Questions). After end processing, Lig4 and XRCC4 ligate the two DNA ends [64].

Alt-EJ

DSBs generated during CSR are also repaired by the Alt-EJ pathway [47]. Factors that are implicated in Alt-EJ include CtIP [65,66], which seems to have a role in end processing, XRCC1, ligase 1, and ligase 3, which are involved in end joining [67,68]. HMCES (5-Hydroxymethylcytosine binding, ES-cell-specific) acts as an abasic site sensor in ssDNA [69], and also functions in the Alt-EJ pathway during CSR [70]. HMCES-defective mouse primary B cells and CH12F3 B cells have shown reduced microhomology usage in switch junctions relative to WT cells [70]. This finding has been postulated to be mechanistically due to HMCES’s ability to bind ssDNA overhangs, likely caused by strand resection of DSB at the S regions, through its SOS-response-associated-peptidase domain [71].

AID Targeting

The mechanism by which AID preferentially mutates Ig genes is still incompletely understood (see Outstanding questions). This issue has implications for immunity as well as for the oncogenic mutations introduced by AID outside the Ig loci [72]. AID targeting entails successive steps of increasing specificity: entering the nucleus, associating quite extensively to chromatin, and mutating Ig loci at a high frequency, in addition to a small fraction of other genes at low (but significant) frequency.

Entering the Nucleus

Several mechanisms regulate the amount of AID protein per B cell present at steady state, as previously reviewed [72,73]. AID is a nuclear–cytoplasmic shuttling protein largely excluded from the nucleus at steady state. Restraining nuclear residence limits AID activity, which may balance Ig diversification and off-target AID activity. The small proportion of AID present in the nucleus during the cell cycle interphase reflects a dynamic equilibrium between opposing mechanisms. AID is retained in the cytoplasm; a heat shock protein 90–DnaJ homolog subfamily A member 1 (HSP90–DNAJA1) complex first stabilizes AID [9]. This molecular chaperone stabilizes proteins until they acquire functional proficiency. For AID, this step is defined by its transition to another cytoplasmic complex that contains elongation factor 1α (eEF1A) and holds the bulk of functional AID, preventing its diffusion into the nucleus [74]. The maturation step required for this transition, and whether this process is regulated, remain unknown [74]. In addition, nuclear AID is either continuously exported by chromosomal maintenance 1 (CRM1) or targeted for proteasomal degradation [9]. Three mechanisms can allow AID into the nucleus. First, AID continuously enters the nucleus through a nuclear localization signal recognized by karyopherins [9]. Second, AID is trapped in the nucleus when the nuclear membrane reassembles at the end of mitosis [44]. This transient accumulation correlates with the time when more AID-induced mutations are observed, and is consistent with SHM and CSR being initiated in the G1 phase [72]. Third, AID can transiently and stochastically accumulate in the nucleus in short pulses [75]. This provocative observation was made in non-B cells (i.e., human fibrosarcoma line HT1080 and SV40-transformed fibroblast line GM639) overexpressing AID–GFP, so the pulses could not be correlated with SHM or CSR. Unlike AID pulses, AID nuclear import and entrapment during mitosis have both been confirmed in B cells by several groups [23,44,76., 77., 78.]. Yet, it is unclear why a nuclear import mechanism for AID would evolve if passive trapping was sufficient to access the nucleus for SHM and CSR. While different AID fractions could conceivably access the nucleus by different pathways, this issue remains unknown, but certainly warrants further research.

Associating with Chromatin and Transcription

AID most preferentially mutates Ig genes but it also deaminates other genes (termed off-target genes) in normal B cells. In addition, the available data are consistent with AID occupying up to thousands of loci, of which just ~450 genes have been found mutated in normal B cells [23,79., 80., 81., 82.] (Boxes 2 and 3 ). These findings suggest the existence of AID susceptible and nonsusceptible loci [72]. If the occupancy of a locus by AID does not necessarily result in mutagenesis, what determines AID activity after it reaches chromatin? It is important to identify the loci occupied by AID to understand its targeting mechanisms. The only available analysis of endogenous AID chromatin occupancy was performed by ChIP-seq in primary activated mouse B cells, which identified >5000 loci occupied by AID [23,82]. Others reanalyzed the data and reported no AID-specific peaks [111], which was rebutted, indicating that AID ChIP-seq requires specific analyses [112]. ChiP-seq peaks are defined by signal enrichment over background. The signal of AID ChIP is typically low, as only ~5% of the B cell AID is chromatin-associated according to biochemical fractionation [23,83]. The background is primarily determined by the specificity of the antibody. The AID ChIP-seq study reported two experiments using different antibodies with similar results. They also used primary B cells from AID knockout mice as negative controls [23,82]. Chromatin fragmentation by mechanical shearing can be a source of nonspecific signals. Highly transcribed DNA regions are hypersensitive to sonication, resulting in an over-representation of these regions in the background. AID ChIP-seq revealed 5514 genes occupied in WT cells, of which 396 were shared with AID-deficient samples, and were also regions that were hypersensitive to sonication, (e.g., the Igh locus). Given the preference of AID for accessible chromatin, this study distinguished specific AID peaks from the background in AID-deficient B cells by considering this hypersensitivity effect for normalizing the signal [23,82,112]. While independent confirmation of these data is desirable, a large number of occupied loci is consistent with the broadly distributed chromatin-bound endogenous AID observed in the CH12F3 mouse B cell line by immunofluorescence, and by the biochemical abundance of chromatin-bound AID [23,83]. A stepwise process of AID targeting might harmonize these observations with the much smaller number of genes that AID has been found to mutate [79., 80., 81.,102,113] (Figure I). Nuclear AID could broadly associate to the chromatin through ribonucleoproteins [72], as suggested by RNA dependency [23]. Interaction with paused RNAPII–SPT5 could allow accumulation around the TSS, consistent with the CHIP-seq profile of endogenous AID in B cells [82]; this would then enable its access to selected gene bodies based on the mechanism discussed in Figure 3 in the main text. This model is analogous to how transcription or DNA repair factors find their target sequences by scanning chromatin, and provides a rationale for the broad association of AID with chromatin and its presence around the TSS of unmutated genes.
Figure I

Stepwise Model for Productive Deamination of Selected AID Off-Targets.

While AID may be able to associate to chromatin in different ways, it could accumulate at paused promoter regions of many genes by virtue of its interaction with SPT5, but access only a subset of those genes bodies, for instance by a licensing mechanism and/or the abundance of ssDNA substrate. Abbreviations: AID, Activation-induced Cytidine Deaminase; RNAPII, RNA polymerase II.

Figure 3

Multilayered Productive AID Targeting in B cells.

(A) AID activity is restricted to type A chromatin compartments, which are accessible and transcriptionally active, versus transcriptionally inactive B compartments. (B) AID activity is also limited to some TADs of DNA, often affecting highly transcribed genes regulated by superenhancers, with clusters by long-range interactions between regulatory sequences. The Ig and some off-target genes additionally harbor cis-acting elements named DIVAC (diversification activators) that attract AID activity. (C, D) AID can be recruited to promoter or promoter-proximal paused RNAPII via SPT5. (C) AID accesses the gene body of susceptible loci through associations with the elongating RNAPII complex, via a putative licensing mechanism. At these loci, several mechanisms can provide ssDNA AID substrates (yellow boxes): DNA supercoiling, DNA secondary structures, and RNA processing of sense and antisense RNA. Transcriptional stalling, caused by nucleic acid secondary structures or by early termination due to antisense transcription, can also expose and favor AID jumping to ssDNA, possibly helped by the regulated association of AID to the ssDNA binding protein RPA. Stalling and SPT5 association may also allow time for AID to act on the small transcription bubble within the RNAPII. (D) At loci that are not susceptible to mutation, AID may fail its coupling with transcription elongation, and/or lack access to ssDNA substrates. Abbreviations: AID, Activation-induced Cytidine Deaminase; ncRNA, noncoding RNA; NELF, negative elongation factor (causes RNAPII pausing); PAF, RNAPII-associated factor complex; RNAPII, RNA polymerase II; ssDNA, single-stranded DNA; TAD, topologically associated domain; TSS, transcription start site.

Stepwise Model for Productive Deamination of Selected AID Off-Targets. While AID may be able to associate to chromatin in different ways, it could accumulate at paused promoter regions of many genes by virtue of its interaction with SPT5, but access only a subset of those genes bodies, for instance by a licensing mechanism and/or the abundance of ssDNA substrate. Abbreviations: AID, Activation-induced Cytidine Deaminase; RNAPII, RNA polymerase II. Alt-text: Box 2 Three strategies have been used to identify AID off-target genes. First, sequencing of select genes revealed a mutation pattern resembling SHM in human B cell lymphoma samples expressing AID, notably diffuse large B-cell lymphoma (DLBCL) [114]. The identity of these relatively few genes can be biased by selection of mutations advantageous for the tumor cells and, accordingly, include many proto-oncogenes or tumor suppressors. Second, mid- and high-throughput sequencing analyses have been undertaken in normal, activated, or germinal center B cells, in mouse models. These efforts still sampled a small proportion of the genome, selecting genes based on defined criteria [79,80]. These approaches have also revealed that some AID-deaminated genes are subsequently repaired by UNG-BER and/or MMR with high efficiency [79,80]. Third, a series of studies have relied on methods that detect and/or trap AID-induced DSB [102,113,115]. These require off-target mutations to be converted to DSB, likely biasing the identification towards genes with the highest deamination load and/or displaying features that facilitate DSB. The latter could include high transcription and/or AID occupancy, enrichment in the WRCY motifs preferred by AID, the formation of secondary structures, or a combination thereof. For instance, G4 DNA may be preferred by AID and may amplify its activity by inducing AID oligomerization (Box 1) [18]. The overlap between the lists of AID off-target genes determined by each method is only partial, demonstrating that each approach has a different bias, which must be considered in AID targeting models. The implementation of unsupervised analysis of mutational signatures in large cancer sequence datasets has expanded the potential off-targets of AID in human samples, but this represents correlative data [116,117]. It would be possible to use a similar unbiased approach to obtain the full repertoire of AID off-targets in normal B cells by comparing germinal center B cells obtained from AID-proficient and -deficient mice, perhaps in combination with UNG and/or MMR deficiencies. Alt-text: Box 3 Biochemical fractionation indicates that ~50% of nuclear AID is associated with chromatin in mouse primary B cells and the CH12F3 B cell line [23,83]. The phosphorylated form of AID at Ser38 (pS38), a protein kinase A substrate, is enriched in the chromatin-associated AID fraction [83], which allows its interaction with the ssDNA-binding replication protein A (RPA), boosting CSR and SHM [84,85]. Notably, pS38 seems dispensable for AID activity on Myc and increasing pS38 increases Igh–Myc translocations without affecting CSR [86]. This suggests that pS38, and potentially other post-translational modifications, can provide context-dependent regulation of AID activity at different loci. The Ig and off-target genes share some features that suggest similarities between the underlying mutation mechanism. One basic common requirement for AID targeting is transcription [72]. AID associates with the RNA polymerase II (RNAPII) complex and with multiple transcription, splicing, and RNA processing factors, as well as with RNA–DNA hybrid structures arising from transcription, which have all been implicated in AID targeting [9,72]. As these are all ubiquitous factors, none by themselves can explain productive AID targeting; thus, they must be integrated into a comprehensive picture that in conjunction with chromatin architecture, could account for AID specificity.

Chromatin Architecture

Chromatin architecture has emerged as an important determinant of genome-wide mutagenic AID targeting. AID-susceptible loci are associated with highly accessible chromatin and significantly correlate with histone marks typical of active enhancers (i.e., H3K27Ac and H3K4me1) and transcription elongation (i.e., H3K36me3 and H3K79me2) [87]. AID activity is enriched at promoter-proximal exons and introns, but can also be detected at enhancers that form regulatory clusters with other promoters and enhancers through long-range interactions [81]. Most AID-susceptible clusters regulated by these superenhancers are highly interconnected and show greater accessibility than other regions (Figure 3A, B). However, the proportion of AID off-targets regulated by superenhancers ranges between 40 and 70% [80,81,88]. The discrepancy may originate from the different techniques used to detect AID off-target activity (Box 3). Multilayered Productive AID Targeting in B cells. (A) AID activity is restricted to type A chromatin compartments, which are accessible and transcriptionally active, versus transcriptionally inactive B compartments. (B) AID activity is also limited to some TADs of DNA, often affecting highly transcribed genes regulated by superenhancers, with clusters by long-range interactions between regulatory sequences. The Ig and some off-target genes additionally harbor cis-acting elements named DIVAC (diversification activators) that attract AID activity. (C, D) AID can be recruited to promoter or promoter-proximal paused RNAPII via SPT5. (C) AID accesses the gene body of susceptible loci through associations with the elongating RNAPII complex, via a putative licensing mechanism. At these loci, several mechanisms can provide ssDNA AID substrates (yellow boxes): DNA supercoiling, DNA secondary structures, and RNA processing of sense and antisense RNA. Transcriptional stalling, caused by nucleic acid secondary structures or by early termination due to antisense transcription, can also expose and favor AID jumping to ssDNA, possibly helped by the regulated association of AID to the ssDNA binding protein RPA. Stalling and SPT5 association may also allow time for AID to act on the small transcription bubble within the RNAPII. (D) At loci that are not susceptible to mutation, AID may fail its coupling with transcription elongation, and/or lack access to ssDNA substrates. Abbreviations: AID, Activation-induced Cytidine Deaminase; ncRNA, noncoding RNA; NELF, negative elongation factor (causes RNAPII pausing); PAF, RNAPII-associated factor complex; RNAPII, RNA polymerase II; ssDNA, single-stranded DNA; TAD, topologically associated domain; TSS, transcription start site. An overlapping of superenhancers with gene bodies of susceptible loci can define the site of mutagenesis, at least in part, by producing antisense transcription that clashes with transcription from the gene promoter. This convergent transcription can stall the RNAPII, causing early termination, and/or recruit RNA processing factors to favor AID activity by exposing ssDNA on both strands [72,89] (Figure 3C). In transcription-linked in vitro deamination assays, the transcript would normally protect the template strand from deamination, as the nontemplate strand is usually mutated. However, purified AID can mutate both strands of an Ig V region when transcribed in the presence of SPT4 and SPT5, a heterodimeric complex that regulates RNAPII processivity [90]. This suggests that AID has an intrinsic ability to mutate both strands under those conditions, but it is likely that in the more complex context of chromatin, some target genes require convergent transcription to expose both DNA strands to AID [88,89]. Chromatin organization into topologically associated domains (TADs) can contribute to AID activity by bringing together distant DNA regions, as exemplified by chromatin extrusion to ensure CSR directionality [91]. Deletional CSR requires long-range interactions between the donor (Eμ) and acceptor (Eacc) S-region enhancers, and the Igh 3′regulatory region (3′RR) [92]. The loop extrusion mechanism envisages cohesin-loading cis elements in the Igh, drawing the 3′RR and Eμ together to form a CSR center. Transcription of an acceptor S region induced by cytokine signaling would expose Eacc cohesin-binding sites, driving chromatin extrusion for juxtaposing the Eacc with the CSR center, and thus arranging distal AID-induced DSBs in a direction favorable for deletional CSR. Additionally, long-range inter-TAD interactions promoted by noncoding RNA can contribute to AID function. Specifically, disrupting transcription by CRISPR/Cas9 on a TAD containing a long noncoding RNA that is 2.6 Mb away from the Igh is sufficient to halve the efficiency of CSR to IgA in mouse Peyer’s patch B cells and in the CH12F3 B cell line [93].

Gene Deamination Susceptibility

Since neither transcription, AID occupancy, nor the presence of superenhancers is sufficient to predict mutation susceptibility of a gene, additional features must contribute to define susceptibility to AID. The presence of cis elements, quality of transcription, formation of DNA/RNA structures, and a licensing mechanism coupling AID to transcription elongation have been proposed [94] (reviewed later).

Cis Elements

Ig loci are privileged for SHM [95], implying that AID is preferentially targeted to and/or more active at those loci. The cause of this preference is still unclear but is independent of the Ig promoter. Ig loci have strong enhancers and are highly transcribed, and at least the Igh locus harbors a superenhancer [81]. Ig enhancers contain defined evolutionarily conserved elements, known as DIVAC elements (diversification activators), that are sufficient to attract SHM to a juxtaposed transcribed sequence [96]. DIVACs can affect mutability at the TAD level [97]. Thus, the insertion of a DIVAC into an otherwise nontargeted TAD can attract SHM, as shown in the Ramos B cell line [97]. In this study, adjacent TADs were spared, additionally suggesting a topologically delineated confinement of AID mutagenesis [97]. The same study found that some AID off-target genes have associated enhancers with DIVAC activity [97], which may explain why they are colaterally mutated by AID. The mechanism of DIVAC activity is unclear. Evidence in the DT40 and Ramos B cell lines suggests that recruitment of different transcription factors to DIVACs enhances SHM without an effect on transcription [96,97]. Nevertheless, their relative contribution remains to be dissected. DIVACs are transcribed, producing enhancer RNA (eRNA) from both strands [97]. Bidirectional transcription, either convergent (as discussed above) or divergent transcription, and eRNA processing by the RNA exosome [89] can contribute to the AID targeting activity of the DIVAC (Figure 3C,D). Although eRNA has not been shown to be a universal feature of AID targets, it is possible that some intragenic AID off-target sites are due to bidirectional transcription driven by overlapping intragenic enhancers [88,89].

Transcription Quality and Nucleic Acid Structures

RNAPII and its processivity factor SPT5 are required for AID activity and colocalize with AID genome-wide [23,98]. SPT5 mediates the AID–RNAPII interaction and is sufficient to recruit AID to chromatin [23,98]. The function of SPT5 is to initially cause promoter-proximal pausing of RNAPII [99], and is thus enriched immediately downstream from the transcription start site (TSS). AID shows a similar genome-wide pattern by ChIP-seq in primary B cells, consistent with its accumulation at SPT5-rich sites (Box 2) [82,98]. During transcription elongation, SPT5 remains associated with and acts as a positive factor for elongating RNAPII; thus, both accumulate inside gene bodies upon transcriptional stalling, a phenomenon considered key for AID activity [90,98,100,101]. Accordingly, stalled RNAPII and SPT5 are found within gene bodies at the Ig V and S regions, as well as at some AID off-targets [72,100,102], which could recruit AID to stalled transcription [98,100]. However, SPT5 is still abundant at the Ig constant regions and at many genes that are occupied by AID, but spared from mutation [80,81,88,98]. Since AID is an inefficient enzyme in vitro, RNAPII stalling via SPT5 or RNA–DNA structures might not only expose ssDNA but also provide time for AID to work [15,90] (Figure 3C). DNA–RNA R loops or G4 DNA formed during transcription can stall RNAPII [103], and in fact, are preferred AID substrates in vitro [20]; this might in turn contribute to AID susceptibility at certain loci. Another targeting mechanism has been proposed in which G4-RNA formed in the intron of the S-region transcript binds AID and brings it to the cDNA [104,105] (Box 1). However, while R loops and G4 can form in the S region and some AID off-targets such as Myc [19,20,89], their presence at the Ig V and at other AID off-targets, is unlikely and uncertain, respectively; therefore, other mechanisms must be at play for SHM [20] (Box 1). Transcriptional stalling may cause premature termination, another potential source of ssDNA demonstrated at the Ig V [101]. In addition, transcription-induced DNA supercoiling also provides ssDNA substrate to AID [106].

Coupling to Transcription Elongation

The identification of catalytically active AID mutants that can occupy the Ig promoter, but not the gene body, and fail to enable SHM or CSR, has suggested a model in which AID is recruited by SPT5 to promoter-proximal-paused RNAPII; in this model, AID gains access to the gene body by traveling with transcription elongation [23]. Supporting the link between AID activity and transcription elongation, factors that interact with RNAPII after the pause release, namely SPT5, RNAPII-associated factor (PAF) and the histone chaperones SPT6 and FAcilitates Chromatin Transcription (FACT), are important for SHM and/or CSR, and biochemically copurify with AID [72,107,108]. Whether histone chaperones facilitate AID activity at the Ig by their role in transcription or by defining histone modifications, is unclear. Of note, the RNA exosome has been linked to elongating RNAPII via SPT6 in Drosophila cells [109] which may facilitate AID activity. A more precise understanding of how the transition from transcriptional pausing/stalling to elongation has an impact on AID activity is required. A speculative model integrating several observations would have AID accessing the gene body with elongating RNAPII, but only being able to mutate genes in which the transcriptional and genomic contexts cause simultaneous RNAPII stalling and ssDNA exposure (Figure 3C and Outstanding Questions). In addition, the dissociation of AID from Ig loci to prevent mutation of constant regions must be somehow regulated – a step that has received little attention.

Concluding Remarks

Research carried out over the past two decades has provided a fundamental understanding of B cell-mediated adaptive immunity by uncovering the biological functions of AID, including its biochemical and structural properties, as well as the DNA repair pathways downstream of AID, and the regulatory mechanisms that facilitate productive AID targeting at Ig loci. Although there is no evidence of a specific AID targeting factor, the implication of ubiquitous transcription-associated factors, and the verification that several factors contribute to, but are not essential for SHM or CSR, suggests a multilayered AID targeting mechanism in which chromatin architecture, transcription quality, and AID regulation all contribute to define whether a locus is mutated or not, and to what extent. This is a fundamental mechanism to further understand, as it has direct implications in affinity maturation, and might eventually be manipulated to improve the efficiency of immunizations. In addition, future endeavors should be directed towards developing a better understanding of the biochemical events involved in antibody diversification processes, including enzymology and mechanisms of action during SHM and CSR. What is the precise mechanism leading to mutagenic repair of AID-induced uracil residues within Ig loci? Does the Shieldin complex in CSR have a similar mechanism of action as Shelterin-mediated telomere end processing? Is the enzymology for the SHM and CSR processes completely understood? How much each of the proposed AID targeting mechanisms contributes to SHM, CSR, and off-target mutagenesis? Which of these mechanisms are in fact linked and which are independent of each other? What exact combination of genomic features and factors promotes AID activity after occupancy? What is the mechanistic role of G4 and R-loop structures in AID targeting to Ig and other genes? Is there a coordination between the transcriptional landscape of susceptible genes and AID post-translational modifications or licensing? Alt-text: Outstanding Questions
  116 in total

1.  DNA polymerase eta is an A-T mutator in somatic hypermutation of immunoglobulin variable genes.

Authors:  X Zeng; D B Winter; C Kasmer; K H Kraemer; A R Lehmann; P J Gearhart
Journal:  Nat Immunol       Date:  2001-06       Impact factor: 25.606

2.  AID-RNA polymerase II transcription-dependent deamination of IgV DNA.

Authors:  Phuong Pham; Sohail Malik; Chiho Mak; Peter C Calabrese; Robert G Roeder; Myron F Goodman
Journal:  Nucleic Acids Res       Date:  2019-11-18       Impact factor: 16.971

3.  AID is required to initiate Nbs1/gamma-H2AX focus formation and mutations at sites of class switching.

Authors:  Michel C Nussenzweig; André Nussenzweig; Simone Petersen; Rafael Casellas; Bernardo Reina-San-Martin; Hua Tang Chen; Michael J Difilippantonio; Patrick C Wilson; Leif Hanitsch; Arkady Celeste; Masamichi Muramatsuk; Duane R Pilch; Christophe Redon; Thomas Ried; William M Bonner; Tasuku Honjo
Journal:  Nature       Date:  2001-12-06       Impact factor: 49.962

4.  Noncoding RNA transcription alters chromosomal topology to promote isotype-specific class switch recombination.

Authors:  Gerson Rothschild; Wanwei Zhang; Junghyun Lim; Pankaj Kumar Giri; Brice Laffleur; Yiyun Chen; Mingyan Fang; Yuling Chen; Lekha Nair; Zhi-Ping Liu; Haiteng Deng; Lennart Hammarström; Jiguang Wang; Uttiya Basu
Journal:  Sci Immunol       Date:  2020-02-07

5.  AID binds to transcription-induced structures in c-MYC that map to regions associated with translocation and hypermutation.

Authors:  Michelle L Duquette; Phuong Pham; Myron F Goodman; Nancy Maizels
Journal:  Oncogene       Date:  2005-09-01       Impact factor: 9.867

Review 6.  The biochemistry of activation-induced deaminase and its physiological functions.

Authors:  Mani Larijani; Alberto Martin
Journal:  Semin Immunol       Date:  2012-06-12       Impact factor: 11.130

7.  The concerted action of Msh2 and UNG stimulates somatic hypermutation at A . T base pairs.

Authors:  Darina Frieder; Mani Larijani; Cathy Collins; Marc Shulman; Alberto Martin
Journal:  Mol Cell Biol       Date:  2009-07-13       Impact factor: 4.272

8.  A source of the single-stranded DNA substrate for activation-induced deaminase during somatic hypermutation.

Authors:  Xiaohua Wang; Manxia Fan; Susan Kalis; Lirong Wei; Matthew D Scharff
Journal:  Nat Commun       Date:  2014-06-13       Impact factor: 14.919

9.  Sequence-Intrinsic Mechanisms that Target AID Mutational Outcomes on Antibody Genes.

Authors:  Leng-Siew Yeap; Joyce K Hwang; Zhou Du; Robin M Meyers; Fei-Long Meng; Agnė Jakubauskaitė; Mengyuan Liu; Vinidhra Mani; Donna Neuberg; Thomas B Kepler; Jing H Wang; Frederick W Alt
Journal:  Cell       Date:  2015-11-12       Impact factor: 41.582

10.  Spt5 accumulation at variable genes distinguishes somatic hypermutation in germinal center B cells from ex vivo-activated cells.

Authors:  Robert W Maul; Zheng Cao; Lakshmi Venkataraman; Carol A Giorgetti; Joan L Press; Yves Denizot; Hansen Du; Ranjan Sen; Patricia J Gearhart
Journal:  J Exp Med       Date:  2014-10-06       Impact factor: 14.307

View more
  33 in total

1.  Ig Enhancers Increase RNA Polymerase II Stalling at Somatic Hypermutation Target Sequences.

Authors:  Alina Tarsalainen; Yaakov Maman; Fei-Long Meng; Minna K Kyläniemi; Anni Soikkeli; Paulina Budzyńska; Jessica J McDonald; Filip Šenigl; Frederic W Alt; David G Schatz; Jukka Alinikula
Journal:  J Immunol       Date:  2021-12-03       Impact factor: 5.422

2.  Controllable genome editing with split-engineered base editors.

Authors:  Kiara N Berríos; Niklaus H Evitt; Rachel A DeWeerd; Diqiu Ren; Meiqi Luo; Aleksia Barka; Tong Wang; Caroline R Bartman; Yemin Lan; Abby M Green; Junwei Shi; Rahul M Kohli
Journal:  Nat Chem Biol       Date:  2021-10-18       Impact factor: 15.040

3.  SHLD1 is dispensable for 53BP1-dependent V(D)J recombination but critical for productive class switch recombination.

Authors:  Estelle Vincendeau; Wenming Wei; Xuefei Zhang; Cyril Planchais; Wei Yu; Hélène Lenden-Hasse; Thomas Cokelaer; Juliana Pipoli da Fonseca; Hugo Mouquet; David J Adams; Frederick W Alt; Stephen P Jackson; Gabriel Balmus; Chloé Lescale; Ludovic Deriano
Journal:  Nat Commun       Date:  2022-06-28       Impact factor: 17.694

4.  APOBEC3A regulates transcription from interferon-stimulated response elements.

Authors:  Manabu Taura; John A Frank; Takehiro Takahashi; Yong Kong; Eriko Kudo; Eric Song; Maria Tokuyama; Akiko Iwasaki
Journal:  Proc Natl Acad Sci U S A       Date:  2022-05-12       Impact factor: 12.779

5.  The uncharacterized SANT and BTB domain-containing protein SANBR inhibits class switch recombination.

Authors:  Simin Zheng; Allysia J Matthews; Numa Rahman; Kayleigh Herrick-Reynolds; Emily Sible; Jee Eun Choi; Alec Wishnie; Yan Kee Ng; Daniela Rhodes; Stephen J Elledge; Bao Q Vuong
Journal:  J Biol Chem       Date:  2021-04-06       Impact factor: 5.157

Review 6.  Genome edited B cells: a new frontier in immune cell therapies.

Authors:  Geoffrey L Rogers; Paula M Cannon
Journal:  Mol Ther       Date:  2021-09-24       Impact factor: 11.454

7.  Role of Dot1L and H3K79 methylation in regulating somatic hypermutation of immunoglobulin genes.

Authors:  Zhi Duan; Linda B Baughn; Xiaohua Wang; Yongwei Zhang; Varun Gupta; Thomas MacCarthy; Matthew D Scharff; Guojun Yu
Journal:  Proc Natl Acad Sci U S A       Date:  2021-07-20       Impact factor: 11.205

Review 8.  AID in Chronic Lymphocytic Leukemia: Induction and Action During Disease Progression.

Authors:  Pablo Oppezzo; Marcelo Navarrete; Nicholas Chiorazzi
Journal:  Front Oncol       Date:  2021-05-10       Impact factor: 6.244

Review 9.  The Ultimate (Mis)match: When DNA Meets RNA.

Authors:  Benoit Palancade; Rodney Rothstein
Journal:  Cells       Date:  2021-06-08       Impact factor: 7.666

10.  AID and APOBECs as Multifaceted Intrinsic Virus-Restricting Factors: Emerging Concepts in the Light of COVID-19.

Authors:  Anastasia Meshcheryakova; Peter Pietschmann; Philip Zimmermann; Igor B Rogozin; Diana Mechtcheriakova
Journal:  Front Immunol       Date:  2021-07-01       Impact factor: 7.561

View more

北京卡尤迪生物科技股份有限公司 © 2022-2023.