Somatic hypermutation in immunity and cancer: Critical analysis of strand-biased and codon-context mutation signatures

Deaminase-Driven Reverse Transcription Mutagenesis in Oncogenesis: Critical Analysis of Transcriptional Strand Asymmetries of Single Base Substitution Signatures

This paper provides a critical analysis of the molecular mechanisms presently used to explain transcriptional strand asymmetries of single base substitution (SBS) signatures observed in cancer genomes curated at the Catalogue of Somatic Mutations in Cancer (COSMIC) database (Wellcome Trust Sanger Institute). The analysis is based on a deaminase-driven reverse transcriptase (DRT) mutagenesis model of cancer oncogenesis involving both the cytosine (AID/APOBEC) and adenosine (ADAR) mutagenic deaminases. In this analysis we apply what is known, or can reasonably be inferred, of the immunoglobulin somatic hypermutation (Ig SHM) mechanism to the analysis of the transcriptional stand asymmetries of the COSMIC SBS signatures that are observed in cancer genomes. The underlying assumption is that somatic mutations arising in cancer genomes are driven by dysregulated off-target Ig SHM-like mutagenic processes at non-Ig loci. It is reasoned that most SBS signatures whether of "unknown etiology" or assigned-molecular causation, can be readily understood in terms of the DRT-paradigm. These include the major age-related "clock-like" SBS5 signature observed in all cancer genomes sequenced and many other common subset signatures including SBS1, SBS3, SBS2/13, SBS6, SBS12, SBS16, SBS17a/17b, SBS19, SBS21, as well as signatures clearly arising from exogenous causation. We conclude that the DRT-model provides a plausible molecular framework that augments our current understanding of immunogenetic mechanisms driving oncogenesis. It accommodates both what is known about AID/APOBEC and ADAR somatic mutation strand asymmetries and provides a fully integrated understanding into the molecular origins of common COSMIC SBS signatures. The DRT-paradigm thus provides scientists and clinicians with additional molecular insights into the causal links between deaminase-associated genomic signatures and oncogenic processes.

Somatic mutation patterns at Ig and Non-Ig Loci

The reverse transcriptase (RT) model of immunoglobulin (Ig) somatic hypermutation (SHM) has received insufficient scientific attention. This is understandable given that DNA deamination mediated by activation-induced deaminase (AID), the initiating step of Ig SHM, has dominated experiments since 2002. We summarise some key history of the RT Ig SHM model dating to 1987. For example, it is now established that DNA polymerase η, the sole DNA repair polymerase involved in post-replication short-patch repair, is an efficient cellular RT. This implies that it is potentially able to initiate target site reverse transcription by RNA-directed DNA repair at AID-induced lesions. Recently, DNA polymerase θ has also been shown to be an efficient cellular RT. Since DNA polymerase θ plays no significant role in Ig SHM, it could serve a similar RNA-dependent DNA polymerase role as DNA polymerase η at non-Ig loci in the putative RNA-templated nucleotide excision repair of bulky adducts and other mutagenic lesions on the transcribed strand. A major yet still poorly recognised consequence of the proposed RT process in Ig SHM is the generation of significant and characteristic strand-biased mutation signatures at both deoxyadenosine/deoxythymidine and deoxyguanosine/deoxycytidine base pairs. In this historical perspective, we highlight how diagnostic strand-biased mutation signatures are detected in vivo during SHM at both Ig loci in germinal centre B lymphocytes and non-Ig loci in cancer genomes. These strand-biased signatures have been significantly obscured by technical issues created by improper use of the polymerase chain reaction technique. A heightened awareness of this fact should contribute to better data interpretation and somatic mutation pattern recognition both at Ig and non-Ig loci.

The transcription factor E2A can bind to and cleave single-stranded immunoglobulin heavy chain locus DNA

Class switch recombination (CSR) changes the constant region of the immunoglobulin heavy chain (IgH), and somatic hypermutation (SH) introduces point mutations in the variable regions of the antibody genes. Both these processes that optimize antibody responses of B lymphocytes are initiated by the enzyme Activation Induced cytidine Deaminase (AID). Here we have searched for CSR or SH coupled activities of the transcription factor E2A, since E2A is in a complex with AID and the transcription factors PAX5, ETS1 and IRF4 on key sequences of the Igh locus in B lymphocytes activated to CSR and SH. We report that E2A in contrast to other described transcription factors binds sequence specifically also to single-stranded DNA. The binding of E2A to single-stranded DNA has a strong sequence preference for one strand of a site in the intronic enhancer of the Igh locus. Furthermore, E2A was also found to cleave single-stranded DNA. The sequence profile of substrates cleaved by E2A is coupled to the sequences of substrates and products of AID, suggesting that E2A has a role not only in targeting of AID to switch regions and SH parts of antibody genes but also in cleavage of DNA at these sites.

Stochastic models of Mendelian and reverse transcriptional inheritance in state-structured cancer populations

Recent evidence suggests that a polyaneuploid cancer cell (PACC) state may play a key role in the adaptation of cancer cells to stressful environments and in promoting therapeutic resistance. The PACC state allows cancer cells to pause cell division and to avoid DNA damage and programmed cell death. Transition to the PACC state may also lead to an increase in the cancer cell’s ability to generate heritable variation (evolvability). One way this can occur is through evolutionary triage. Under this framework, cells gradually gain resistance by scaling hills on a fitness landscape through a process of mutation and selection. Another way this can happen is through self-genetic modification whereby cells in the PACC state find a viable solution to the stressor and then undergo depolyploidization, passing it on to their heritably resistant progeny. Here, we develop a stochastic model to simulate both of these evolutionary frameworks. We examine the impact of treatment dosage and extent of self-genetic modification on eco-evolutionary dynamics of cancer cells with aneuploid and PACC states. We find that under low doses of therapy, evolutionary triage performs better whereas under high doses of therapy, self-genetic modification is favored. This study generates predictions for teasing apart these biological hypotheses, examines the implications of each in the context of cancer, and provides a modeling framework to compare Mendelian and non-traditional forms of inheritance.

RNA-directed DNA repair and antibody somatic hypermutation

Somatic hypermutation at antibody loci affects both deoxyadenosine-deoxythymidine (A/T) and deoxycytidine-deoxyguanosine (C/G) pairs. Deamination of C to deoxyuridine (U) by activation-induced deaminase (AID) explains how mutation at C/G pairs is potentiated. Mutation at A/T pairs is triggered during the initial stages of repair of AID-generated U lesions and occurs through an as yet unknown mechanism in which polymerase η has a major role. Recent evidence confirms that human polymerase η can act as a reverse transcriptase. Here, we compare the popular suggestion of mutation at A/T pairs through nucleotide mispairing (owing to polymerase error) during short-patch repair synthesis with the alternative proposal of mutation at A/T pairs through RNA editing and RNA-directed DNA repair.

The efficient Lamarckian spread of life in the cosmos

In this Chapter we discuss the various mechanisms that are available for the possible transfer of cosmic microbial living systems from one cosmic habitat to another. With the 100 or so habitable planets that are now known to exist in our galaxy alone transfers of cometary dust carrying life including fragments of icy planetoids/asteroids would be expected to occur on a routine basis. It is thus easy to view the galaxy as a single connected "biosphere" of which our planet Earth is a minor component. The Hoyle-Wickramasinghe Panspermia paradigm provides a cogent biological rationale for the actual widespread existence of Lamarckian modes of inheritance in terrestrial systems (which we review here). Thus the Panspermia paradigm provides the raison d'etre for Lamarckian Inheritance. Under a terrestrially confined neoDarwinian viewpoint such an association may have been thought spurious in the past. Our aim here is to outline the main evidence for rapid terrestrial-based Lamarckian-based evolutionary hypermutation processes dependent on reverse transcription-coupled mechanisms among others. Such rapid adaptation mechanisms would be consistent with the effective cosmic spread of living systems. For example, a viable, or cryo-preserved, living system traveling through space in a protective matrix will of necessity need to adapt rapidly and proliferate on landing in a new cosmic niche. Lamarckian mechanisms thus come to the fore and supersede the slow (blind and random) genetic processes expected under neoDarwinian Earth centred theories.

A proposed reverse transcription mechanism for (CAG)n and similar expandable repeats that cause neurological and other diseases

The mechanism of (CAG)n repeat generation, and related expandable repeat diseases in non-dividing cells, is currently understood in terms of a DNA template-based DNA repair synthesis process involving hairpin stabilized slippage, local error-prone repair via MutSβ (MSH2-MSH3) hairpin protective stabilization, then nascent strand extension by DNA polymerases-β and -δ. We advance a very similar slipped hairpin-stabilized model involving MSH2-MSH3 with two key differences: the copying template may also be the nascent pre-mRNA with the repair pathway being mediated by the Y-family error-prone enzymes DNA polymerase-η and DNA polymerase-κ acting as reverse transcriptases. We argue that both DNA-based and RNA-based mechanisms could well be activated in affected non-dividing brain cells in vivo. Here, we compare the advantages of the RNA/RT-based model proposed by us as an adjunct to previously proposed models. In brief, our model depends upon dysregulated innate and adaptive immunity cascades involving AID/APOBEC and ADAR deaminases that are known to be involved in normal locus-specific immunoglobulin somatic hypermutation, cancer progression and somatic mutations at many off-target non-immunoglobulin sites across the genome: we explain how these processes could also play an active role in repeat expansion diseases at RNA polymerase II-transcribed genes.

Regulatory T cells and co-evolution of allele-specific MHC recognition by the TCR

What is the evolutionary mechanism for the TCR-MHC-conserved interaction? We extend Dembic's model (Dembic Z. In, Scand J Immunol e12806, 2019) of thymus positive selection for high-avidity anti-self-MHC Tregs among double (CD4 + CD8+)-positive (DP) developing thymocytes. This model is based on competition for self-MHC (+ Pep) complexes presented on cortical epithelium. Such T cells exit as CD4 + CD25+FoxP3 + thymic-derived Tregs (tTregs). The other positively selected DP T cells are then negatively selected on medulla epithelium removing high-avidity anti-self-MHC + Pep as T cells commit to CD4 + or CD8 + lineages. The process is likened to the competitive selection and affinity maturation in Germinal Centre for the somatic hypermutation (SHM) of rearranged immunoglobulin (Ig) variable region (V[D]Js) of centrocytes bearing antigen-specific B cell receptors (BCR). We now argue that the same direct SHM processes for TCRs occur in post-antigenic Germinal Centres, but now occurring in peripheral pTregs. This model provides a potential solution to a long-standing problem previously recognized by Cohn and others (Cohn M, Anderson CC, Dembic Z. In, Scand J Immunol e12790, 2019) of how co-evolution occurs of species-specific MHC alleles with the repertoire of their germline TCR V counterparts. We suggest this is not by 'blind', slow, and random Darwinian natural selection events, but a rapid structured somatic selection vertical transmission process. The pTregs bearing somatic TCR V mutant genes then, on arrival in reproductive tissues, can donate their TCR V sequences via soma-to-germline feedback as discussed in this journal earlier. (Steele EJ, Lindley RA. In, Scand J Immunol e12670, 2018) The high-avidity tTregs also participate in the same process to maintain a biased, high-avidity anti-self-MHC germline V repertoire.

Lamarck and Panspermia - On the Efficient Spread of Living Systems Throughout the Cosmos

We review the main lines of evidence (molecular, cellular and whole organism) published since the 1970s demonstrating Lamarckian Inheritance in animals, plants and microorganisms viz. the transgenerational inheritance of environmentally-induced acquired characteristics. The studies in animals demonstrate the genetic permeability of the soma-germline Weismann Barrier. The widespread nature of environmentally-directed inheritance phenomena reviewed here contradicts a key pillar of neo-Darwinism which affirms the rigidity of the Weismann Barrier. These developments suggest that neo-Darwinian evolutionary theory is in need of significant revision. We argue that Lamarckian inheritance strategies involving environmentally-induced rapid directional genetic adaptations make biological sense in the context of cosmic Panspermia allowing the efficient spread of living systems and genetic innovation throughout the Universe. The Hoyle-Wickramasinghe Panspermia paradigm also developed since the 1970s, unlike strictly geocentric neo-Darwinism provides a cogent biological rationale for the actual widespread existence of Lamarckian modes of inheritance - it provides its raison d'être. Under a terrestrially confined neo-Darwinian viewpoint such an association may have been thought spurious in the past. Our aim is to outline the conceptual links between rapid Lamarckian-based evolutionary hypermutation processes dependent on reverse transcription-coupled mechanisms among others and the effective cosmic spread of living systems. For example, a viable, or cryo-preserved, living system travelling through space in a protective matrix will need of necessity to rapidly adapt and proliferate on landing in a new cosmic niche. Lamarckian mechanisms thus come to the fore and supersede the slow (blind and random) genetic processes expected under a traditional neo-Darwinian evolutionary paradigm.

DNA damage pathways and B-cell lymphomagenesis

PURPOSE OF REVIEW

Recent lymphoma genome sequencing projects have shed light on the genomic landscape of indolent and aggressive lymphomas, as well as some of the molecular mechanisms underlying recurrent mutations and translocations in these entities. Here, we review these recent genomic discoveries, focusing on acquired DNA repair defects in lymphoma. In addition, we highlight recently identified actionable molecular vulnerabilities associated with recurrent mutations in chronic lymphocytic leukemia (CLL), which serves as a model entity.

RECENT FINDINGS

The results of several large lymphoma genome sequencing projects have recently been reported, including CLL, T-PLL and DLBCL. We align these discoveries with proposed mechanisms of mutation acquisition in B-cell lymphomas. Moreover, novel autochthonous mouse models of CLL have recently been generated and we discuss how these models serve as preclinical tools to drive the development of novel targeted therapeutic interventions. Lastly, we highlight the results of early clinical data on novel compounds targeting defects in the DNA damage response of CLL with a particular focus on deleterious ATM mutations.

SUMMARY

Defects in DNA repair pathways are selected events in cancer, including lymphomas. Specifically, ATM deficiency is associated with PARP1- and DNA-PKcs inhibitor sensitivity in vitro and in vivo.

Impact of DNA lesion repair, replication and formation on the mutational spectra of environmental carcinogens: Aflatoxin B1 as a case study

In a multicellular organism, somatic mutations represent a permanent record of the past chemical and biochemical perturbations experienced by a cell in its local microenvironment. Akin to a perpetual recording device, with every replication, genomic DNA accumulates mutations in patterns that reflect: i) the sequence context-dependent formation of DNA damage, due to environmental or endogenous reactive species, including spontaneous processes; ii) the activity of DNA repair pathways, which, depending on the type of lesion, can erase, ignore or exacerbate the mutagenic consequences of that DNA damage; and iii) the choice of replication machinery that synthesizes the nascent genomic copy. These three factors result in a richly contoured sequence context-dependent mutational spectrum that, from appearances, is distinct for most individual forms of DNA damage. Such a mutagenic legacy, if appropriately decoded, can reveal the local history of genome-altering events such as chemical or pathogen exposures, metabolic stress, and inflammation, which in turn can provide an indication of the underlying causes and mechanisms of genetic disease. Modern tools have positioned us to develop a deep mechanistic understanding of the cellular factors and pathways that modulate a mutational process and, in turn, provide opportunities for better diagnostic and prognostic biomarkers, better exposure risk assessment and even actionable therapeutic targets. The goal of this Perspective is to present a bottom-up, lesion-centric framework of mutagenesis that integrates the contributions of lesion replication, lesion repair and lesion formation to explain the complex mutational spectra that emerge in the genome following exposure to mutagens. The mutational spectra of the well-studied hepatocarcinogen aflatoxin B₁ are showcased here as specific examples, but the implications are meant to be generalizable.

Regulation of the DNA Repair Complex during Somatic Hypermutation and Class-Switch Recombination

B lymphocytes optimize Ab responses by somatic hypermutation (SH), which introduces point mutations in the variable regions of the Ab genes and by class-switch recombination (CSR), which changes the expressed C region exon of the IgH. These Ab diversification processes are initiated by the deaminating enzyme activation-induced cytidine deaminase followed by many DNA repair enzymes, ultimately leading to deletions and a high mutation rate in the Ab genes, whereas DNA lesions made by activation-induced cytidine deaminase are repaired with low error rate on most other genes. This indicates an advanced regulation of DNA repair. In this study, we show that initiation of Ab diversification in B lymphocytes of mouse spleen leads to formation of a complex between many proteins in DNA repair. We show also that BCR activation, which signals the end of successful SH, reduces interactions between some proteins in the complex and increases other interactions in the complex with varying kinetics. Furthermore, we show increased localization of SH- and CSR-coupled proteins on switch regions of the Igh locus upon initiation of SH/CSR and differential changes in the localization upon BCR signaling, which terminates SH. These findings provide early evidence for a DNA repair complex or complexes that may be of functional significance for carrying out essential roles in SH and/or CSR in B cells.

Cause of Cambrian Explosion - Terrestrial or Cosmic?

We review the salient evidence consistent with or predicted by the Hoyle-Wickramasinghe (H-W) thesis of Cometary (Cosmic) Biology. Much of this physical and biological evidence is multifactorial. One particular focus are the recent studies which date the emergence of the complex retroviruses of vertebrate lines at or just before the Cambrian Explosion of ∼500 Ma. Such viruses are known to be plausibly associated with major evolutionary genomic processes. We believe this coincidence is not fortuitous but is consistent with a key prediction of H-W theory whereby major extinction-diversification evolutionary boundaries coincide with virus-bearing cometary-bolide bombardment events. A second focus is the remarkable evolution of intelligent complexity (Cephalopods) culminating in the emergence of the Octopus. A third focus concerns the micro-organism fossil evidence contained within meteorites as well as the detection in the upper atmosphere of apparent incoming life-bearing particles from space. In our view the totality of the multifactorial data and critical analyses assembled by Fred Hoyle, Chandra Wickramasinghe and their many colleagues since the 1960s leads to a very plausible conclusion - life may have been seeded here on Earth by life-bearing comets as soon as conditions on Earth allowed it to flourish (about or just before 4.1 Billion years ago); and living organisms such as space-resistant and space-hardy bacteria, viruses, more complex eukaryotic cells, fertilised ova and seeds have been continuously delivered ever since to Earth so being one important driver of further terrestrial evolution which has resulted in considerable genetic diversity and which has led to the emergence of mankind.

High-throughput immune repertoire analysis with IGoR

High-throughput immune repertoire sequencing is promising to lead to new statistical diagnostic tools for medicine and biology. Successful implementations of these methods require a correct characterization, analysis, and interpretation of these data sets. We present IGoR (Inference and Generation Of Repertoires)-a comprehensive tool that takes B or T cell receptor sequence reads and quantitatively characterizes the statistics of receptor generation from both cDNA and gDNA. It probabilistically annotates sequences and its modular structure can be used to investigate models of increasing biological complexity for different organisms. For B cells, IGoR returns the hypermutation statistics, which we use to reveal co-localization of hypermutations along the sequence. We demonstrate that IGoR outperforms existing tools in accuracy and estimate the sample sizes needed for reliable repertoire characterization.

Reverse Transcriptase Mechanism of Somatic Hypermutation: 60 Years of Clonal Selection Theory

The evidence for the reverse transcriptase mechanism of somatic hypermutation is substantial and multifactorial. In this 60th anniversary year of the publication of Sir MacFarlane Burnet's Clonal Selection Theory, the evidence is briefly reviewed and updated.

ADAR deaminase A-to-I editing of DNA and RNA moieties of RNA:DNA hybrids has implications for the mechanism of Ig somatic hypermutation

The implications are discussed of recently published biochemical studies on ADAR-mediated A-to-I DNA and RNA deamination at RNA:DNA hybrids. The significance of these data are related to previous work on strand-biased and codon-context mutation signatures in B lymphocytes and cancer genomes. Those studies have established that there are two significant strand biases at A:T and G:C base pairs, A-site mutations exceed T-site mutations (A>>T) by 2.9 fold and G-site mutations exceed C-site mutations (G>>C) by 1.7 fold. Both these strand biases are inconsistent with alternative "DNA Deamination" mechanisms, yet are expected consequences of the RNA/RT-based "Reverse Transcriptase" mechanism of immunoglobulin (Ig) somatic hypermutation (SHM). The A-to-I DNA editing component at RNA:DNA hybrids that is likely to occur in Transcription Bubbles, while important, is of far lower A-to-I editing efficiency than in dsRNA substrates. The RNA moiety of RNA:DNA hybrids is also edited at similar lower frequencies relative to the editing rate at dsRNA substrates. Further, if the A-to-I DNA editing at RNA:DNA hybrids were the sole cause of A-to-I (read as A-to-G) mutation events for Ig SHM in vivo then the exact opposite strand biases at A:T base pairs (T>>A) of what is actually observed (A>>T) would be predicted. It is concluded that the strand-biased somatic mutation patterns at both A:T and G:C base pairs in vivo are best interpreted by the sequential steps of the RNA/RT-based mechanism. Further, the direct DNA A-to-I deamination at Transcription Bubbles is expected to contribute to the T-to-C component of the strand-biased Ig SHM spectrum.

BRILIA: Integrated Tool for High-Throughput Annotation and Lineage Tree Assembly of B-Cell Repertoires

The somatic diversity of antigen-recognizing B-cell receptors (BCRs) arises from Variable (V), Diversity (D), and Joining (J) (VDJ) recombination and somatic hypermutation (SHM) during B-cell development and affinity maturation. The VDJ junction of the BCR heavy chain forms the highly variable complementarity determining region 3 (CDR3), which plays a critical role in antigen specificity and binding affinity. Tracking the selection and mutation of the CDR3 can be useful in characterizing humoral responses to infection and vaccination. Although tens to hundreds of thousands of unique BCR genes within an expressed B-cell repertoire can now be resolved with high-throughput sequencing, tracking SHMs is still challenging because existing annotation methods are often limited by poor annotation coverage, inconsistent SHM identification across the VDJ junction, or lack of B-cell lineage data. Here, we present B-cell repertoire inductive lineage and immunosequence annotator (BRILIA), an algorithm that leverages repertoire-wide sequencing data to globally improve the VDJ annotation coverage, lineage tree assembly, and SHM identification. On benchmark tests against simulated human and mouse BCR repertoires, BRILIA correctly annotated germline and clonally expanded sequences with 94 and 70% accuracy, respectively, and it has a 90% SHM-positive prediction rate in the CDR3 of heavily mutated sequences; these are substantial improvements over existing methods. We used BRILIA to process BCR sequences obtained from splenic germinal center B cells extracted from C57BL/6 mice. BRILIA returned robust B-cell lineage trees and yielded SHM patterns that are consistent across the VDJ junction and agree with known biological mechanisms of SHM. By contrast, existing BCR annotation tools, which do not account for repertoire-wide clonal relationships, systematically underestimated both the size of clonally related B-cell clusters and yielded inconsistent SHM frequencies. We demonstrate BRILIA’s utility in B-cell repertoire studies related to VDJ gene usage, mechanisms for adenosine mutations, and SHM hot spot motifs. Furthermore, we show that the complete gene usage annotation and SHM identification across the entire CDR3 are essential for studying the B-cell affinity maturation process through immunosequencing methods.

Chemical Synthesis, Versatile Structures and Functions of Tailorable Adjuvants for Optimizing Oral Vaccination

Oral vaccines have become a recent focus because of their potential significance in disease prevention and therapy. In the development of oral vaccine-based therapeutics, synthetic materials with tailorable structures and versatile functions can act as antigen conveyers with adjuvant effects, reduce the time cost for vaccine optimization, and provide high security and enhanced immunity. This review presents an overview of the current status of tailoring synthetic adjuvants for oral vaccination, modification strategies for producing effectors with specific structures and functions, enhancement of immune-associated efficiencies, including the barrier-crossing capability to protect antigens in the gastrointestinal tract, coordination of the antigens penetrating mucosa and cell barriers, targeting of concentrated antigens to immune-associated cells, and direct stimulation of immune cells. Finally, we focus on prospective synthetic adjuvants that facilitate the use of oral vaccines via two approaches, namely, in vivo antigen expression and cancer immunotherapy.

Inosine-specific ribonuclease activity of natural variants of human endonuclease V

Adenine bases in DNA, RNA, and nucleotides are deaminated during normal metabolism via hydrolytic and nitrosative reactions. In RNA, the deaminated product inosine is resolved by human endonuclease V, and mice deficient in this enzyme are cancer-prone. We have now produced, purified, and characterized naturally occurring variants of human endonuclease V (V29I, R112Q, K114R, H141Y, and D201N). We found that H141Y, but not other variants, is catalytically impaired, suggesting that individuals homozygous for H141Y may be predisposed to disease.