The signature of somatic hypermutation appears to be written into the germline IgV segment repertoire

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.

Germline V repertoires: Origin, maintenance, diversification

In our view, Melvin Cohn (Scand J Immunol. 2018;87:e12640) has set out the logical guidelines towards a resolution of the very real enigma of the selectability of vertebrate germline Ig V repertoires under the current evolutionary paradigm…" A somatically derived repertoire scrambles this (germline VL + VH) substrate so that its specificities are lost, making it un-selectable in the germline. Consequently, evolution faced an incompatibility." It is argued here in Reply that a reverse transcriptase-based soma-to-germline process (S->G) targeting germline V segment arrays goes some considerable way to resolving fundamental contradictions on the origin, maintenance and then real-time adaptive diversification of these limited sets of V segments encoded within various V repertoire arrays.

APOBEC and ADAR deaminases may cause many single nucleotide polymorphisms curated in the OMIM database

Cytosine and adenosine deamination events (DNA, RNA substrates) account for most codon-context Targeted Somatic Mutation (TSM) patterns observed in immunoglobulin (Ig) somatic hypermutation (SHM), and in cancer exomes following Ig-SHM-like responses. TSM refers to the process of somatic mutagenesis involving deamination events that results on a dominant type of mutation (e.g., C-to-T), and co-incident at a particular motif (e.g., WRC), and preferentially targeting the first, second or third nucleotide position within the mutated codon (e.g. MC1, MC2 or MC3, read 5-prime to 3-prime). It is now widely accepted that if left uncorrected, the accumulation of uncorrected TSMs involving the deaminases, may lead to a diagnosis of cancer or other degenerative disease. Our hypothesis is that many missense, nonsense and synonymous single nucleotide polymorphisms (SNPs) associated with clinically significant diseases may have arisen in the population by similar highly targeted deamination events. The OMIM database was searched for disease-associated SNPs on the X chromosome, and for all chromosomes. The nucleotide substitution patterns for disease-associated SNPs were analyzed by the TSM method to identify the likely deaminase source for C-to-U (C-to-T/G-to-A) and A-to-I (A-to-G/T-to-C) derived gene mutations preferentially targeting known sequence motifs associated with the deaminases: AID, APOBEC3G, APOBEC3B and ADAR 1/2. Of the 789 OMIM SNPs analysed. In both data sets, the mutation targeting preferences within the mutated codon reveal a statistically significant bias (p < 0.001). The results imply that a deamination of C-site and A-site targets are written into the human germline for the chromosome wide exomic SNPs analysed. This is consistent with previously observed mutation patterns arising in cancer genomes and hypermutated Ig genes during SHM. The results imply that similar types of deaminase-mediated molecular processes that occur in somatic hypermutation and cancer, may be contributing causative drivers of human exomic SNPs.

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.

Soma-to-germline feedback is implied by the extreme polymorphism at IGHV relative to MHC: The manifest polymorphism of the MHC appears greatly exceeded at Immunoglobulin loci, suggesting antigen-selected somatic V mutants penetrate Weismann's Barrier

Soma-to-germline feedback is forbidden under the neo-Darwinian paradigm. Nevertheless, there is a growing realization it occurs frequently in immunoglobulin (Ig) variable (V) region genes. This is a surprising development. It arises from a most unlikely source in light of the exposure of co-author EJS to the haplotype data of RL Dawkins and others on the polymorphism of the Major Histocompatibility Complex, which is generally assumed to be the most polymorphic region in the genome (spanning ∼4 Mb). The comparison between the magnitude of MHC polymorphism with estimates for the human heavy chain immunoglobulin V locus (spanning ∼1 Mb), suggests IGHV could be many orders of magnitude more polymorphic than the MHC. This conclusion needs airing in the literature as it implies generational churn and soma-to-germline gene feedback. Pedigree-based experimental strategies to resolve the IGHV issue are outlined.

A unified approach to the evolutionary consequences of genetic and nongenetic inheritance

Inheritance-the influence of ancestors on the phenotypes of their descendants-translates natural selection into evolutionary change. For the past century, inheritance has been conceptualized almost exclusively as the transmission of DNA sequence variation from parents to offspring in accordance with Mendelian rules, but advances in cell and developmental biology have now revealed a rich array of inheritance mechanisms. This empirical evidence calls for a unified conception of inheritance that combines genetic and nongenetic mechanisms and encompasses the known range of transgenerational effects, including the transmission of genetic and epigenetic variation, the transmission of plastic phenotypes (acquired traits), and the effects of parental environment and genotype on offspring phenotype. We propose a unified theoretical framework based on the Price equation that can be used to model evolution under an expanded inheritance concept that combines the effects of genetic and nongenetic inheritance. To illustrate the utility and generality of this framework, we show how it can be applied to a variety of scenarios, including nontransmissible environmental noise, maternal effects, indirect genetic effects, transgenerational epigenetic inheritance, RNA-mediated inheritance, and cultural inheritance.

Common chromosomal fragile sites (CFS) may be involved in normal and traumatic cognitive stress memory consolidation and altered nervous system immunity

Previous reports of specific patterns of increased fragility at common chromosomal fragile sites (CFS) found in association with certain neurobehavioural disorders did not attract attention at the time due to a shift towards molecular approaches to delineate neuropsychiatric disorder candidate genes. Links with miRNA, altered methylation and the origin of copy number variation indicate that CFS region characteristics may be part of chromatinomic mechanisms that are increasingly linked with neuroplasticity and memory. Current reports of large-scale double-stranded DNA breaks in differentiating neurons and evidence of ongoing DNA demethylation of specific gene promoters in adult hippocampus may shed new light on the dynamic epigenetic changes that are increasingly appreciated as contributing to long-term memory consolidation. The expression of immune recombination activating genes in key stress-induced memory regions suggests the adoption by the brain of this ancient pattern recognition and memory system to establish a structural basis for long-term memory through controlled chromosomal breakage at highly specific genomic regions. It is furthermore considered that these mechanisms for management of epigenetic information related to stress memory could be linked, in some instances, with the transfer of the somatically acquired information to the germline. Here, rearranged sequences can be subjected to further selection and possible eventual retrotranscription to become part of the more stable coding machinery if proven to be crucial for survival and reproduction. While linkage of cognitive memory with stress and fear circuitry and memory establishment through structural DNA modification is proposed as a normal process, inappropriate activation of immune-like genomic rearrangement processes through traumatic stress memory may have the potential to lead to undesirable activation of neuro-inflammatory processes. These theories could have a significant impact on the interpretation of risks posed by heredity and the environment and the search for neuropsychiatric candidate genes.

Mechanism of somatic hypermutation: critical analysis of strand biased mutation signatures at A:T and G:C base pairs

The DNA sequence data of the somatic hypermutation (SHM) field published since 1984 has been critically reviewed. The analysis has revealed three strand biased mutation signatures. The first concerns the mutations generated at G:C base pairs in mice genetically deficient in uracil-DNA glycosylase and MSH2-MSH6-mediated mismatch repair. Such mice display the AID deaminase footprint and here C mutations exceed G mutations at least 1.5-fold. This supports earlier and more recent studies claiming that dC-to-dU deaminations occur preferentially in the single stranded DNA regions of the displaced nontranscribed strand (NTS) during transcription. The second concerns the signature generated in immunised mice where G mutations exceed C mutations by at least 1.7-fold. This is a newly identified strand bias which has previously gone undetected. It is consistent with the polynucleotide polymerisation signature of RNA polymerase II copying the template DNA strand carrying AID-mediated lesions generated at C bases, viz. uracils and abasic sites. A reverse transcription step would then need to intervene to fix the mutation pattern in DNA. The third concerns the long recognised strand biased signature generated in normal aged or actively immunised mice whereby A mutations exceed T mutations by two- to three-fold. It is argued that this pattern is best understood as a combination of adenosine-to-inosine (A-to-I) RNA editing followed by a reverse transcription step fixing the A-to-G, as well as A-to-T and A-to-C, as strand biased mutation signatures in DNA. The reasons why the AID-linked RNA polymerase II mutation signature had previously gone undetected are discussed with regard to limitations of standard PCR-based SHM assay techniques. It is concluded that the most economical SHM mechanism involves both DNA and RNA deaminations coupled to a reverse transcription process, most likely involving DNA polymerase eta acting in its reverse transcriptase mode. Experimental approaches to differentiate this RNA-based model from the standard DNA deamination model are discussed.

Substitution patterns in alleles of immunoglobulin V genes in humans and mice

Immunoglobulins (Igs) constitute a subfamily of rapidly evolving proteins. It is postulated that this characteristic is due mainly to the participation of these proteins in highly diverse functions of recognition and defense. Although this vision of rapid evolution in Igs is widely accepted, various studies have demonstrated that diverse and contradictory forces not yet completely understood converge in the evolution of these receptors. In a recent study of the substitution patterns in the alleles that form the human IGHV locus, we found that the variation in genetic and structural information does not occur homogeneously among the different genes, nor among the regions and positions conforming said locus. In view of these results and of the importance of a better understanding of the basic evolutionary process in specific receptors (such as Igs) for both immunology and molecular evolution, it is important to explore the nature of the diversification process in these proteins in detail. In this work, therefore, we analyzed the substitution patterns in all the alleles reported for loci IGKV and IGLV in humans and mice, and we compared the results with those previously observed in the human IGHV locus. We found that the process of evolutionary variation of the Igs reflect the diversity of selective pressures operating on the different loci, genes, sub-regions and positions; for example, diversification through substitution is generally centered on CDRs, but only few positions inside the CDRs were frequently substituted. In spite of this general tendency, it is possible to observe differences in the degree of diversification among loci, families and genes. These tendencies to modify only certain attributes of IGV genes seem to be in agreement with differential strategies associated with the restrictions of the molecular immune recognition mechanism. The complexity of the evolutionary patterns observed in this study leads us to think that the predispositions observed herein may also be due in part to processes of DNA dynamics.

Structural analysis of substitution patterns in alleles of human immunoglobulin VH genes

The diversity in repertoires of antibodies (Abs) needed in response to the antigen challenge is produced by evolutionary and somatic processes. The mechanisms operating at a somatic level have been studied in great detail. In contrast, neither the mechanisms nor the strategies of diversification at an evolutionary level have yet been understood in similar detail. Particularly, the substitution patterns in alleles of immunoglobulin genes (Igs) have not been systematically studied. Furthermore, there is a scarcity of studies which link the analysis at a genetic level of the diversification of repertoires with the structural consequences at the protein level of the changes in DNA information. For the purpose of systematically characterizing the strategies of evolutionary diversification through sequence variation at alleles, in this work, we built a database for all the alleles of the IGHV locus in humans reported until now. Based on these data, we performed diverse analyses of substitution patterns and linked these results with studies at the protein level. We found that the sequence diversification in different alleles does not operate with equal intensity for all V genes. Our studies, both of the number of substitutions and of the type of amino acid change per sub-segment of the V-REGION evidenced differences in the selective pressure to which these regions are exposed. The implications of these results for understanding the evolutionary diversification strategies, as well as for the somatic generation of antibody repertoires are discussed.

A 30 kb region of the Epstein-Barr virus genome is colinear with the rearranged human immunoglobulin gene loci: implications for a "ping-pong evolution" model for persisting viruses and their hosts. A review

The left part of the Epstein-Barr virus (EBV) genome exhibits a strong colinearity of structural and functional elements with the immunoglobulin (Ig) gene loci which is only partially reflected in nucleotide sequence homologies. We propose that this colinearity may be the result of an inter-dependent co-evolution of the immunoglobulin loci together with EBV. Our observation could help elucidating the mechanisms of somatic hypermutation, explaining the ability of EBV to accidentally cause tumors, and shedding more light on the general mechanisms of viral and organismal evolution. We suggest that persisting viruses served as a complement for the organismal germline like in a ping-pong game and outline The Ping-Pong Evolution Hypothesis.

Molecular evolution of V(H)9 germline genes isolated from DBA, BALB, 129 and C57BL mouse strains and sublines

We have used the polymerase chain reaction (PCR) in an attempt to clone and sequence the exons and hitherto unavailable contiguous flanks of all members of the small V(H) 9 germline gene family from inbred mouse strains and sublines that have had a common ancestry within the last century, and to analyze the molecular evolution of these sequences. Fifteen genuine germline genes were isolated (designated V(H) 9.1 through V(H) 9.15) from strains and sublines of DBA, BALB, 129 and C57BL inbred mice. Of the 15 genuine isolates, nine are novel: seven sequences from DBA strains and sublines ( V(H) 9.3 to V(H) 9.9) and two sequences from C57BL strains ( V(H) 9.13 and V(H) 9.14). We have identified sequencing errors and PCR recombinant artefacts in previously published sequences. We detected no sequence divergence of individual genes shared by the strains and sublines studied. However, we isolated two genes from DBA strains and sublines, V(H) 9.1 and V(H) 9.3, that differ only by five nucleotides encoding three amino acid changes that are concentrated within a 33 nucleotide (11 codon) region. Of these 11 codons, eight encode a putative antigen binding site. There were no differences in the remaining 733 nucleotides sequenced (including both 5' and 3' flanking regions). Potential explanations for the generation of V(H) 9.1 and V(H) 9.3 are discussed.

A model of the early evolution of soma-to-germline feedback

The V-genes of the immunoglobulin locus in vertebrates code for a part of the heavy and light chain variable regions of antibodies and are extremely variable. Steele (1979) has developed a theory that explains the evolution of adaptive immune response by a soma-to-germline flow of cDNAs derived from somatically mutated V-genes. Here we model the early evolution of soma-to-germline feedback in a population living in a changing viral environment in terms of the dynamics of an initially rare genetic modifier that controls transfer of V-genes to germ cells' DNA. It is shown that a modifier invades the population and creates a great variety of V-genes if the environment follows stepwise temporal changes, i.e. a soma-to-germline feedback machinery evolves in a population if newly derived V-alleles still play a role in protecting the population against foreign antigens in some following generations. The distribution of the age of V-genes evolves to a bell-shaped curve the width and the maximum of which depend mainly on selection strength. Two phases of modifier evolution are distinguished. In the first phase, the dynamics are slow while the number of different V-genes is small. In the second phase, when a sufficiently large number of different V-genes is created, the modifier increases faster in frequency. Linkage of V-genes and the modifier enhances the rate of evolution.

Analysis of antibodies of known structure suggests a lack of correspondence between the residues in contact with the antigen and those modified by somatic hypermutation

Forty unique murine antibody-antigen complexes determined at 2.5 A or less resolution are analyzed to determine whether the residues in direct contact with the antigen are modified by somatic hypermutation. This was done by taking advantage of the recent characterization of the pool of Vkappa germline genes of the mouse. The average number of residues in contact with the antigen in the V(L) gene, which contains the CDRL-1, CDRL-2, and all but one residue of CDRL-3, was six. The average number of somatic mutations was similar (around five). However, as many as 53% of the antibodies did not show somatic replacements of residues in contact with the antigen. Another 28% had only one. Overall, the frequency of antibodies with increasing number of somatic replacements in residues in contact with the antigen decreased exponentially. A possible explanation of this finding is that mutations in the contacting residues have an adverse effect on the antigen-antibody interaction. This implies that most of the observed mutations are those remaining after negative (purifying) selection. Therefore, efficient strategies of site-directed mutagenesis to improve the affinity of antibodies should be focused on residues other than those directly interacting with the antigen.

Testing the reverse transcriptase model of somatic mutation

Somatic hypermutation of the variable (V) regions of rearranged immunoglobulin genes leads to antibody affinity maturation. Although this process has been extensively studied, the mechanisms responsible for these multiple point mutations are still elusive. One mechanism that was proposed over 10 years ago by Steele and Pollard was that an intrinsic reverse transcriptase (RT) copies the nascent mRNA creating the large number of observed point mutations due to its high error rate. A cDNA copy of the mutated V region would then replace the endogenous DNA through a gene conversion-like event, thus integrating these point mutations into the genome. This model of hypermutation would account for the very high mutation rate, the presence of hotspots, strand bias, the requirement for transcription and localization of mutation within the immunoglobulin V region. Using AZT and ddC to inhibit endogenous RTs, we have assayed for somatic mutation using a murine in vivo model. Somatic mutation occurred at similar frequencies and with the same characteristics with or without treatment of RT inhibitors, suggesting that standard reverse transcription is not required for antibody V region hypermutation in the mouse.

Molecular genetics and pathogenesis of autosomal dominant polycystic kidney disease

Autosomal dominant polycystic kidney disease (ADPKD) is a common and systemic disease characterized by formation of focal cysts. Of the three potential causes of cysts, downstream obstruction, compositional changes in extracellular matrix, and proliferation of partially dedifferentiated cells, evidence strongly supports the latter as the primary abnormality. In the vast majority of cases, the disease is caused by mutations in PKD1 or PKD2, and appears to be recessive at the cellular level. Somatic second hits in the normal allele of cells containing the germ line mutation initiate or accelerate formation of cysts. The intrinsically high frequency of somatic second hits in epithelia appears to be sufficient to explain the frequent occurrence of somatic second hits in the disease-causing genes. PKD1 and PKD2 encode a putative adhesive/ion channel regulatory protein and an ion channel, respectively. The two proteins interact directly in vitro. Their cellular and subcellular localization suggest that they may also function independently in a common signaling pathway that may involve the membrane skeleton and that links cell-cell and cell-matrix adhesion to the development of cell polarity.

The reverse transcriptase model of somatic hypermutation

The evidence supporting the reverse transcriptase model of somatic hypermutation is critically reviewed. The model provides a coherent explanation for many apparently unrelated findings. We also show that the somatic hypermutation pattern in the human BCL-6 gene can be interpreted in terms of the reverse transcriptase model and the notion of feedback of somatically mutated sequences to the germline over evolutionary time.

The genetic variability of the VH genes in follicular lymphoma: the impact of the hypermutation mechanism

Follicular lymphoma (FL) cells have inherited an activated hypermutation mechanism from their origin of germinal centre B cells. Based on today's knowledge of the intrinsic properties related to this mechanism and the VH base composition, reconsideration of previous reports should be made on a broader range of samples. The present study examined the mutation pattern of the VH genes expressed by 55 cases of FL. FL VH genes showed evidence of antigenic selection in 30% of cases with 88% carrying a functional sIg and 78.2% showing intraclonal variation. VH family and gene segment utilization was found to be roughly similar to that of normal B lymphocytes. FL VH genes revealed extensive variations. 17% of the VH exons harboured a total of five deletions, three duplications and two insertions as compared to the most homologous germline counterpart. The VH genes of one tumour displayed three populations with varying CDR3 length at diagnosis. At relapse, emergence of a differently mutated gene, additional mutations reminiscent of ongoing mutations or no variation was prominent. From this study the heterogeneity of FLs is well established and ongoing mutations are seen in the scope of the activated status of the hypermutation mechanism rather than antigen-stimulated tumour growth.

Selection in a T-dependent primary humoral response: new insights from polypeptide models

Regulatory mechanisms involved in the induction and progression of T-dependent humoral responses have been extensively delineated using a variety of haptens as model antigens. However, several unanswered questions remain with respect to those elicited by structurally more complex molecules. Our own laboratory has been pursuing this latter aspect using designed synthetic peptides as model systems. The cumulative results indeed support that humoral responses to such antigens involve several additional layers of regulation, beyond that identified with haptens. At the first level, the multiplicity of antigenic determinants recognized by the preimmune B-cell pool is soon subject to competitive pressures that restrict, both at the level of repertoire and epitope, fine specificities of early activated clonotypes. Selection at this stage is on the basis of affinity for epitope, which, in turn, is under thermodynamic control. This selected B-cell subset proceeds to populate germinal centers, where further optimization--by way of somatic hypermutation followed by clonal selection--is in favor of increased on-rates of antigen binding. Thus, contrary to findings with hapten antigens, maturation of antibody responses to polypeptides occurs in two discrete, but sequential, stages. The first is for B cells with optimum affinity for the corresponding epitope. This is then followed by further improvement on the basis of increased on-rates of antigen/epitope binding. It is a combination of these two processes which results in the high fidelity of antibodies produced in the secondary response.

Alternative Splicing and Hypermutation of a Nonproductively Rearranged TCR α-Chain in a T Cell Hybridoma

Like Ig genes, TCR genes are formed by somatic rearrangements of noncontiguous genomic V, J, and C regions. Unlike Ig genes, somatic hypermutation of TCR V regions is an infrequent event. We describe the occurrence of spontaneous hypermutation in a nonproductively rearranged TCR α-chain gene in a clonal T cell hybridoma that had lost its productively rearranged α-chain. The mutating hybridoma was eventually supplanted in culture by a nonmutating variant that had restored an open reading frame in the nonproductively rearranged TCR α-chain through the use of cryptic splice sites in the Vα region. Evidence is presented for the presence of cDNA reverse transcripts of the TCR α-chain within the hybridoma, suggesting a role for reverse transcriptase in the generation of mutations.