Computational analyses show A-to-G mutations correlate with nascent mRNA hairpins at somatic hypermutation hotspots

Molecular mechanisms behind the generation of pro-oncogenic HIV-1 matrix protein p17 variants

HIV-1 matrix protein p17 variants (vp17s), characterized by amino acid insertions at the COOH-terminal region of the viral protein, have been recently identified and studied for their biological activity. Different from their wild-type counterpart (refp17), vp17s display a potent B cell growth and clonogenic activity. Recent data have highlighted the higher prevalence of vp17s in people living with HIV-1 (PLWH) with lymphoma compared with those without lymphoma, suggesting that vp17s may play a key role in lymphomagenesis. Molecular mechanisms involved in vp17 development are still unknown. Here we assessed the efficiency of HIV-1 Reverse Transcriptase (RT) in processing this genomic region and highlighted the existence of hot spots of mutation in Gag, at the end of the matrix protein and close to the matrix-capsid junction. This is possibly due to the presence of inverted repeats and palindromic sequences together with a high content of Adenine in the 322-342 nucleotide portion, which constrain HIV-1 RT to pause on the template. To define the recombinogenic properties of hot spots of mutation in the matrix gene, we developed plasmid vectors expressing Gag and a minimally modified Gag variant, and measured homologous recombination following cell co-nucleofection by next-generation sequencing. Data obtained allowed us to show that a wide range of recombination events occur in concomitance with the identified hot spots of mutation and that imperfect events may account for vp17s generation.

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.

Molecular Biology and Diversification of Human Retroviruses

Studies of retroviruses have led to many extraordinary discoveries that have advanced our understanding of not only human diseases, but also molecular biology as a whole. The most recognizable human retrovirus, human immunodeficiency virus type 1 (HIV-1), is the causative agent of the global AIDS epidemic and has been extensively studied. Other human retroviruses, such as human immunodeficiency virus type 2 (HIV-2) and human T-cell leukemia virus type 1 (HTLV-1), have received less attention, and many of the assumptions about the replication and biology of these viruses are based on knowledge of HIV-1. Existing comparative studies on human retroviruses, however, have revealed that key differences between these viruses exist that affect evolution, diversification, and potentially pathogenicity. In this review, we examine current insights on disparities in the replication of pathogenic human retroviruses, with a particular focus on the determinants of structural and genetic diversity amongst HIVs and HTLV.

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.

Blood and saliva-derived exomes from healthy Caucasian subjects do not display overt evidence of somatic mosaicism

Somatic mosaicism is a normal occurrence during development in the tissues and organs. As part of establishing a "healthy population "(HP) background or base-line, we investigated whether such mosaicism can be routinely detected in the circulating DNA secured from a rigorously designed healthy human liquid biopsy clinical trial (saliva, blood). We deployed next generation (NG) whole exome sequencing (WES) at median exome coverage rates of 97.2 % (-to-30x) and 70.0 % (-to-100x). We found that somatic mosaicism is not detectable by such standard bulk WES sequencing assays in saliva and blood DNA in 24 normal healthy Caucasians of both sexes from 18 to 60 years of age. We conclude that for circulating DNA using standard WES no novel somatic mutational variants can be detected in protein-coding regions of normal healthy subjects. This implies that the extent within normal tissues of somatic mosaicism must be at a lower level, below the detection threshold, for these circulating DNA WES read depths.

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.

Molecular model linking Th2 polarized M2 tumour-associated macrophages with deaminase-mediated cancer progression mutation signatures

A new and diverse range of somatic mutation signatures are observed in late-stage cancers, but the underlying reasons are not fully understood. We advance a "combinatorial association model" for deaminase binding domain (DBD) diversification to explain the generation of previously observed cancer-progression associated mutation signatures. We also propose that changes in the polarization of tumour-associated macrophages (TAMs) are accompanied by the expression of deaminases with a new and diverse range of DBDs, and thus accounting for the generation of new somatic mutation signatures. The mechanism proposed is molecularly reminiscent of combinatorial association of heavy (H) and light (L) protein chains following V(D)J recombination of immunoglobulin molecules (and similarly for protein chains in heterodimers α/β and γ/δ of V(D)Js of T Cell Receptors) required for pathogen antigen recognition by B cells and T cells, respectively. We also discuss whether extracellular vesicles (EVs) emanating from tumour enhancing M2-polarized macrophages represent a likely source of the de novo deaminase DBDs. We conclude that M2-polarized macrophages extruding EVs loaded with deaminase proteins or deaminase-specific transcription/translation regulatory factors and like information may directly trigger deaminase diversification within cancer cells, and thus account for the many new somatic mutation signatures that are indicative of cancer progression. This hypothesis now has a plausible evidentiary base, and it is worth direct testing in future investigations. A long-term objective would be to identify molecular biomarkers predicting cancer progression (or metastatic disease) and to support the development of new drug targets before metastatic pathways are activated.

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.

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.

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.

Biological significance of RNA editing in cells

RNA editing is one of the post-transcriptional RNA processes. RNA editing generates RNA and protein diversity in eukaryotes and results in specific amino acid substitutions, deletions, and changes in gene expression levels. Adenosine-to-inosine RNA editing represents the most important class of editing in human and affects function of many genes. The importance of balancing RNA modification levels across time and space is becoming increasingly evident. In this review, we overview the biological significance of RNA editing including RNA editing in tumorigenesis, RNA editing in neuronal tissues, RNA editing as a regulator of gene expression, and RNA editing in dsRNA-mediated gene silencing, which may increase our understanding of RNA biology.

Critical Analysis of Strand-Biased Somatic Mutation Signatures in TP53 versus Ig Genes, in Genome-Wide Data and the Etiology of Cancer

Previous analyses of rearranged immunoglobulin (Ig) variable genes (VDJs) concluded that the mechanism of Ig somatic hypermutation (SHM) involves the Ig pre-mRNA acting as a copying template resulting in characteristic strand biased somatic mutation patterns at A:T and G:C base pairs. We have since analysed cancer genome data and found the same mutation strand-biases, in toto or in part, in nonlymphoid cancers. Here we have analysed somatic mutations in a single well-characterised gene TP53. Our goal is to understand the genesis of the strand-biased mutation patterns in TP53—and in genome-wide data—that may arise by “endogenous” mechanisms as opposed to adduct-generated DNA-targeted strand-biased mutations caused by well-characterised “external” carcinogenic influences in cigarette smoke, UV-light, and certain dietary components. The underlying strand-biased mutation signatures in TP53, for many non-lymphoid cancers, bear a striking resemblance to the Ig SHM pattern. A similar pattern can be found in genome-wide somatic mutations in cancer genomes that have also mutated TP53. The analysis implies a role for base-modified RNA template intermediates coupled to reverse transcription in the genesis of many cancers. Thus Ig SHM may be inappropriately activated in many non-lymphoid tissues via hormonal and/or inflammation-related processes leading to cancer.

ADAR1 protein induces adenosine-targeted DNA mutations in senescent Bcl6 gene-deficient cells

Somatic mutations accumulate in senescent cells. Bcl6, which functions as a transcriptional repressor, has been identified as a potent inhibitor of cell senescence, but a role of Bcl6 in the accumulation of somatic mutations has remained unclear. Ig class-switch recombination simultaneously induces somatic mutations in an IgM class-switch (Ig-Sμ) region of IgG B cells. Surprisingly, mutations were detected in the Ig-Sμ region of Bcl6-deficient IgM B cells without class-switch recombination, and these mutations were mainly generated by conversion of adenosine to guanosine, suggesting a novel DNA mutator in the B cells. The ADAR1 (adenosine deaminase acting on RNA1) gene was overexpressed in Bcl6-deficient cells, and its promoter analysis revealed that ADAR1 is a molecular target of Bcl6. Exogenous ADAR1 induced adenosine-targeted DNA mutations in IgM B cells from ADAR1-transgenic mice and in wild-type mouse embryonic fibroblasts (MEFs). These mutations accumulated in senescent MEFs accompanied with endogenous ADAR1 expression, and the frequency in senescent Bcl6-deficient MEFs was higher than senescent wild-type MEFs. Thus, Bcl6 protects senescent cells from accumulation of adenosine-targeted DNA mutations induced by ADAR1.

Genesis of ancestral haplotypes: RNA modifications and reverse transcription-mediated polymorphisms

Understanding the genesis of the block haplotype structure of the genome is a major challenge. With the completion of the sequencing of the Human Genome and the initiation of the HapMap project the concept that the chromosomes of the mammalian genome are a mosaic, or patchwork, of conserved extended block haplotype sequences is now accepted by the mainstream genomics research community. Ancestral Haplotypes (AHs) can be viewed as a recombined string of smaller Polymorphic Frozen Blocks (PFBs). How have such variant extended DNA sequence tracts emerged in evolution? Here the relevant literature on the problem is reviewed from various fields of molecular and cell biology particularly molecular immunology and comparative and functional genomics. Based on our synthesis we then advance a testable molecular and cellular model. A critical part of the analysis concerns the origin of the strand biased mutation signatures in the transcribed regions of the human and higher primate genome, A-to-G versus T-to-C (ratio ∼ 1.5 fold) and C-to-T versus G-to-A (≥ 1.5 fold). A comparison and evaluation of the current state of the fields of immunoglobulin Somatic Hypermutation (SHM) and Transcription-Coupled DNA Repair focused on how mutations in newly synthesized RNA might be copied back to DNA thus accounting for some of the genome-wide strand biases (e.g., the A-to-G vs T-to-C component of the strand biased spectrum). We hypothesize that the genesis of PFBs and extended AHs occurs during mutagenic episodes in evolution (e.g., retroviral infections) and that many of the critical DNA sequence diversifying events occur first at the RNA level, e.g., recombination between RNA strings resulting in tandem and dispersed RNA duplications (retroduplications), RNA mutations via adenosine-to-inosine pre-mRNA editing events as well as error prone RNA synthesis. These are then copied back into DNA by a cellular reverse transcription process (also likely to be error-prone) that we have called "reverse transcription-mediated long DNA conversion." Finally we suggest that all these activities and others can be envisaged as being brought physically under the umbrella of special sites in the nucleus involved in transcription known as "transcription factories."

Nature, position, and frequency of mutations made in a single cycle of HIV-1 replication

There is considerable HIV-1 variation in patients. The extent of the variation is due to the high rate of viral replication, the high viral load, and the errors made during viral replication. Mutations can arise from errors made either by host DNA-dependent RNA polymerase II or by HIV-1 reverse transcriptase (RT), but the relative contributions of these two enzymes to the mutation rate are unknown. In addition, mutations in RT can affect its fidelity, but the effect of mutations in RT on the nature of the mutations that arise in vivo is poorly understood. We have developed an efficient system, based on existing technology, to analyze the mutations that arise in an HIV-1 vector in a single cycle of replication. A lacZalpha reporter gene is used to identify viral DNAs that contain mutations which are analyzed by DNA sequencing. The forward mutation rate in this system is 1.4 x 10(-5) mutations/bp/cycle, equivalent to the retroviral average. This rate is about 3-fold lower than previously reported for HIV-1 in vivo and is much lower than what has been reported for purified HIV-1 RT in vitro. Although the mutation rate was not affected by the orientation of lacZalpha, the sites favored for mutations (hot spots) in lacZalpha depended on which strand of lacZalpha was present in the viral RNA. The pattern of hot spots seen in lacZalpha in vivo did not match any of the published data obtained when purified RT was used to copy lacZalpha in vitro.

Alternative splicing and extensive RNA editing of human TPH2 transcripts

Brain serotonin (5-HT) neurotransmission plays a key role in the regulation of mood and has been implicated in a variety of neuropsychiatric conditions. Tryptophan hydroxylase (TPH) is the rate-limiting enzyme in the biosynthesis of 5-HT. Recently, we discovered a second TPH isoform (TPH2) in vertebrates, including man, which is predominantly expressed in brain, while the previously known TPH isoform (TPH1) is primarly a non-neuronal enzyme. Overwhelming evidence now points to TPH2 as a candidate gene for 5-HT-related psychiatric disorders. To assess the role of TPH2 gene variability in the etiology of psychiatric diseases we performed cDNA sequence analysis of TPH2 transcripts from human post mortem amygdala samples obtained from individuals with psychiatric disorders (drug abuse, schizophrenia, suicide) and controls. Here we show that TPH2 exists in two alternatively spliced variants in the coding region, denoted TPH2a and TPH2b. Moreover, we found evidence that the pre-mRNAs of both splice variants are dynamically RNA-edited in a mutually exclusive manner. Kinetic studies with cell lines expressing recombinant TPH2 variants revealed a higher activity of the novel TPH2B protein compared with the previously known TPH2A, whereas RNA editing was shown to inhibit the enzymatic activity of both TPH2 splice variants. Therefore, our results strongly suggest a complex fine-tuning of central nervous system 5-HT biosynthesis by TPH2 alternative splicing and RNA editing. Finally, we present molecular and large-scale linkage data evidencing that deregulated alternative splicing and RNA editing is involved in the etiology of psychiatric diseases, such as suicidal behaviour.

V-region mutation in vitro, in vivo, and in silico reveal the importance of the enzymatic properties of AID and the sequence environment

The somatic hypermutation of Ig variable regions requires the activity of activation-induced cytidine deaminase (AID) which has previously been shown to preferentially deaminate WRC (W = A/T, R = A/G) motif hot spots in in vivo and in vitro assays. We compared mutation profiles of in vitro assays for the 3' flanking intron of VhJ558-Jh4 region to previously reported in vivo profiles for the same region in the Msh2(-/-)Ung(-/-) mice that lack base excision and mismatch repair. We found that the in vitro and in vivo mutation profiles were highly correlated for the top (nontranscribed) strand, while for the bottom (transcribed) strand the correlation is far lower. We used an in silico model of AID activity to elucidate the relative importance of motif targeting in vivo. We found that the mutation process entails substantial complexity beyond motif targeting, a large part of which is captured in vitro. To elucidate the contribution of the sequence environment to the observed differences between the top and bottom strands, we analyzed intermutational distances. The bottom strand shows an approximately exponential distribution of distances in vivo and in vitro, as expected from a null model. However, the top strand deviates strongly from this distribution in that mutations approximately 50 nucleotides apart are greatly reduced, again both in vivo and in vitro, illustrating an important strand asymmetry. While we have confirmed that AID targeting of hot and cold spots is a key part of the mutation process, our results suggest that the sequence environment plays an equally important role.

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.

Length of the hypermutation motif DGYW/WRCH in the focus of statistical limits. Implications for a double-motif or extended motif recognition models

The motif DGYW/WRCH (Mh) and its frequently discussed simplified derivative GYW/WRC (Mhs) are involved in immunoglobulin (Ig) hypermutation. Both these motifs appear to be markedly shorter than the corresponding conventionally predicted minima of valid sequence lengths (MVSL). The same conclusion concerning both Mh and Mhs can also be obtained in the combined case including a less strict semi-empirically defined w-value and one nucleotide length tolerance related to MVSL. Such disagreement indicates considerably low information content in Mh and Mhs when evaluating these motifs as alphabetical structures (words). This fact raises a question of actually recognized structures (presumably longer than Mh and Mhs). Interestingly, both Mh and Mhs dimers or pairs of closely located Mh or Mhs achieve confirmation of length validity in the case of w=0.05, suggesting thus double-motif recognition as one of statistically consistent explanations. This possibility is also in agreement with the results of our model sequence study of mRNA derived from variable Ig gene sequences (rIgV) with respect to the most frequently occurring structures formed by motif overlaps in all model sequence sets. On the other hand, additional superior occurrence of motif pairs at a structurally important distance of a single DNA thread was found in the conserved domain (cd00099) related sequences of Elasmobranchii origin and less markedly in the corresponding human rIgV, but not in a randomly selected human subset of rIgV. The data are discussed with respect to statistical evaluation and structural properties of hypermutation motifs or the competent enzyme, i.e. activation-induced cytidine deaminase.

I. VH gene transcription creates stabilized secondary structures for coordinated mutagenesis during somatic hypermutation

During the adaptive immune response, antigen challenge triggers a million-fold increase in mutation rates in the variable-region antibody genes. The frequency of mutation is causally and directly linked to transcription, which provides ssDNA and drives supercoiling that stabilizes secondary structures containing unpaired, intrinsically mutable bases. Simulation analysis of transcription in VH5 reveals a dominant 65nt secondary structure in the non-transcribed strand containing six sites of mutable ssDNA that have also been identified independently in human B cell lines and in primary mouse B cells. This dominant structure inter-converts briefly with less stable structures and is formed repeatedly during transcription, due to periodic pauses and backtracking. In effect, this creates a stable yet dynamic "mutability platform" consisting of ever-changing patterns of unpaired bases that are simultaneously exposed and therefore able to coordinate mutagenesis. Such a complex of secondary structures may be the source of ssDNA for enzyme-based diversification, which ultimately results in high affinity antibodies.

Immunoglobulin somatic hypermutation

The immunoglobulin (Ig) repertoire achieves functional diversification through several somatic alterations of the Ig locus. One of these processes, somatic hypermutation (SHM), deposits point mutations into the variable region of the Ig gene to generate higher-affinity variants. Activation-induced cytidine deaminase (AID) converts cytidine to uridine to initiate the hypermutation process. Error-prone versions of DNA repair are believed to then process these lesions into a diverse spectrum of point mutations. We review the current understanding of the molecular mechanisms and regulation of SHM, and also discuss emerging ideas which merit further exploration.

RNA editing, DNA recoding and the evolution of human cognition

RNA editing appears to be the major mechanism by which environmental signals overwrite encoded genetic information to modify gene function and regulation, particularly in the brain. We suggest that the predominance of Alu elements in the human genome is the result of their evolutionary co-adaptation as a modular substrate for RNA editing, driven by selection for higher-order cognitive function. We show that RNA editing alters transcripts from loci encoding proteins involved in neural cell identity, maturation and function, as well as in DNA repair, implying a role for RNA editing not only in neural transmission and network plasticity but also in brain development, and suggesting that communication of productive changes back to the genome might constitute the molecular basis of long-term memory and higher-order cognition.