Cancerous hyper-mutagenesis in p53 genes is possibly associated with transcriptional bypass of DNA lesions

Molecular mechanisms of RNA polymerase II transcription elongation elucidated by kinetic network models

Transcription elongation cycle (TEC) of RNA polymerase II (Pol II) is a process of adding a nucleoside triphosphate to the growing messenger RNA chain. Due to the long timescale events in Pol II TEC, an advanced computational technique, such as Markov State Model (MSM), is needed to provide atomistic mechanism and reaction rates. The combination of MSM and experimental results can be used to build a kinetic network model (KNM) of the whole TEC. This review provides a brief protocol to build MSM and KNM of the whole TEC, along with the latest findings of MSM and other computational studies of Pol II TEC. Lastly, we offer a perspective on potentially using a sequence dependent KNM to predict genome-wide transcription error.

Accurate RNA consensus sequencing for high-fidelity detection of transcriptional mutagenesis-induced epimutations

Transcriptional mutagenesis (TM) due to misincorporation during RNA transcription can result in mutant RNAs, or epimutations, that generate proteins with altered properties. TM has long been hypothesized to play a role in aging, cancer, and viral and bacterial evolution. However, inadequate methodologies have limited progress in elucidating a causal association. We present a high-throughput, highly accurate RNA sequencing method to measure epimutations with single-molecule sensitivity. Accurate RNA consensus sequencing (ARC-seq) uniquely combines RNA barcoding and generation of multiple cDNA copies per RNA molecule to eliminate errors introduced during cDNA synthesis, PCR, and sequencing. The stringency of ARC-seq can be scaled to accommodate the quality of input RNAs. We apply ARC-seq to directly assess transcriptome-wide epimutations resulting from RNA polymerase mutants and oxidative stress.

Stationary-Phase Mutagenesis in Stressed Bacillus subtilis Cells Operates by Mfd-Dependent Mutagenic Pathways

In replication-limited cells of Bacillus subtilis, Mfd is mutagenic at highly transcribed regions, even in the absence of bulky DNA lesions. However, the mechanism leading to increased mutagenesis through Mfd remains currently unknown. Here, we report that Mfd may promote mutagenesis in nutritionally stressed B. subtilis cells by coordinating error-prone repair events mediated by UvrA, MutY and PolI. Using a point-mutated gene conferring leucine auxotrophy as a genetic marker, it was found that the absence of UvrA reduced the Leu⁺ revertants and that a second mutation in mfd reduced mutagenesis further. Moreover, the mfd and polA mutants presented low but similar reversion frequencies compared to the parental strain. These results suggest that Mfd promotes mutagenic events that required the participation of NER pathway and PolI. Remarkably, this Mfd-dependent mutagenic pathway was found to be epistatic onto MutY; however, whereas the MutY-dependent Leu⁺ reversions required Mfd, a direct interaction between these proteins was not apparent. In summary, our results support the concept that Mfd promotes mutagenesis in starved B. subtilis cells by coordinating both known and previously unknown Mfd-associated repair pathways. These mutagenic processes bias the production of genetic diversity towards highly transcribed regions in the genome.

Visualizing translocation dynamics and nascent transcript errors in paused RNA polymerases in vivo

BACKGROUND

Transcription elongation is frequently interrupted by pausing signals in DNA, with downstream effects on gene expression. Transcription errors also induce prolonged pausing, which can lead to a destabilized genome by interfering with DNA replication. Mechanisms of pausing associated with translocation blocks and misincorporation have been characterized in vitro, but not in vivo.

RESULTS

We investigate the pausing pattern of RNA polymerase (RNAP) in Escherichia coli by a novel approach, combining native elongating transcript sequencing (NET-seq) with RNase footprinting of the transcripts (RNET-seq). We reveal that the G-dC base pair at the 5' end of the RNA-DNA hybrid interferes with RNAP translocation. The distance between the 5' G-dC base pair and the 3' end of RNA fluctuates over a three-nucleotide width. Thus, the G-dC base pair can induce pausing in post-translocated, pre-translocated, and backtracked states of RNAP. Additionally, a CpG sequence of the template DNA strand spanning the active site of RNAP inhibits elongation and induces G-to-A errors, which leads to backtracking of RNAP. Gre factors efficiently proofread the errors and rescue the backtracked complexes. We also find that pausing events are enriched in the 5' untranslated region and antisense transcription of mRNA genes and are reduced in rRNA genes.

CONCLUSIONS

In E. coli, robust transcriptional pausing involves RNAP interaction with G-dC at the upstream end of the RNA-DNA hybrid, which interferes with translocation. CpG DNA sequences induce transcriptional pausing and G-to-A errors.

UV signature mutations

Sequencing complete tumor genomes and exomes has sparked the cancer field's interest in mutation signatures for identifying the tumor's carcinogen. This review and meta-analysis discusses signatures and their proper use. We first distinguish between a mutagen's canonical mutations—deviations from a random distribution of base changes to create a pattern typical of that mutagen—and the subset of signature mutations, which are unique to that mutagen and permit inference backward from mutations to mutagen. To verify UV signature mutations, we assembled literature datasets on cells exposed to UVC, UVB, UVA, or solar simulator light (SSL) and tested canonical UV mutation features as criteria for clustering datasets. A confirmed UV signature was: ≥60% of mutations are C→T at a dipyrimidine site, with ≥5% CC→TT. Other canonical features such as a bias for mutations on the nontranscribed strand or at the 3' pyrimidine had limited application. The most robust classifier combined these features with criteria for the rarity of non-UV canonical mutations. In addition, several signatures proposed for specific UV wavelengths were limited to specific genes or species; UV's nonsignature mutations may cause melanoma BRAF mutations; and the mutagen for sunlight-related skin neoplasms may vary between continents.

Direct assessment of transcription fidelity by high-resolution RNA sequencing

Cancerous and aging cells have long been thought to be impacted by transcription errors that cause genetic and epigenetic changes. Until now, a lack of methodology for directly assessing such errors hindered evaluation of their impact to the cells. We report a high-resolution Illumina RNA-seq method that can assess noncoded base substitutions in mRNA at 10⁻⁴–10⁻⁵ per base frequencies in vitro and in vivo. Statistically reliable detection of changes in transcription fidelity through ∼10³ nt DNA sites assures that the RNA-seq can analyze the fidelity in a large number of the sites where errors occur. A combination of the RNA-seq and biochemical analyses of the positions for the errors revealed two sequence-specific mechanisms that increase transcription fidelity by Escherichia coli RNA polymerase: (i) enhanced suppression of nucleotide misincorporation that improves selectivity for the cognate substrate, and (ii) increased backtracking of the RNA polymerase that decreases a chance of error propagation to the full-length transcript after misincorporation and provides an opportunity to proofread the error. This method is adoptable to a genome-wide assessment of transcription fidelity.

Transcriptional de-repression and Mfd are mutagenic in stressed Bacillus subtilis cells

In recent years, it has been proposed that conflicts between transcription and active chromosomal replication engender genome instability events. Furthermore, transcription elongation factors have been reported to prevent conflicts between transcription and replication and avoid genome instability. Here, we examined transcriptional de-repression as a genetic diversity-producing agent and showed, through the use of physiological and genetic means, that transcriptional de-represssion of a leuC defective allele leads to accumulation of Leu(+) mutations. We also showed, by using riboswitches that activate transcription in conditions of tyrosine or methionine starvation, that the effect of transcriptional de-repression of the leuC construct on the accumulation of Leu(+) mutations was independent of selection. We examined the role of Mfd, a transcription elongation factor involved in DNA repair, in this process and showed that proficiency of this factor promotes mutagenic events. These results are in stark contrast to previous reports in Escherichia coli, which showed that Mfd prevents replication fork collapses. Because our assays place cells under non-growing conditions, by starving them for two amino acids, we surmised that the Mfd mutagenic process associated with transcriptional de-repression does not result from conflicts with chromosomal replication. These results raise the interesting concept that transcription elongation factors may serve two functions in cells. In growing conditions, these factors prevent the generation of mutations, while in stress or non-growing conditions they mediate the production of genetic diversity.

Transcriptional mutagenesis: causes and involvement in tumour development

The majority of human cells do not multiply continuously but are quiescent or slow-replicating and devote a large part of their energy to transcription. When DNA damage in the transcribed strand of an active gene is bypassed by a RNA polymerase, they can miscode at the damaged site and produce mutant transcripts. This process is known as transcriptional mutagenesis and, as discussed in this Perspective, could lead to the production of mutant proteins and might therefore be important in tumour development.

Transcription-associated mutation in Bacillus subtilis cells under stress

Adaptive (stationary phase) mutagenesis is a phenomenon by which nondividing cells acquire beneficial mutations as a response to stress. Although the generation of adaptive mutations is essentially stochastic, genetic factors are involved in this phenomenon. We examined how defects in a transcriptional factor, previously reported to alter the acquisition of adaptive mutations, affected mutation levels in a gene under selection. The acquisition of mutations was directly correlated to the level of transcription of a defective leuC allele placed under selection. To further examine the correlation between transcription and adaptive mutation, we placed a point-mutated allele, leuC427, under the control of an inducible promoter and assayed the level of reversion to leucine prototrophy under conditions of leucine starvation. Our results demonstrate that the level of Leu(+) reversions increased significantly in parallel with the induced increase in transcription levels. This mutagenic response was not observed under conditions of exponential growth. Since transcription is a ubiquitous biological process, transcription-associated mutagenesis may influence evolutionary processes in all organisms.

[Mapping DNA damage to understand somatic mutagenesis]

Somatic mutation theory explains how DNA damage can lead to the malignant transformation of cells. It therefore elucidates the connection between genotoxic agents and cancers. Mutational spectra, which tend to be characteristic of a cancer type, are available for certain genes like p53 which is frequently mutated in tumors. A mutational spectrum could therefore be the signature of the genotoxic agent(s) at the origin of the malignant transformation. Ligation-mediated PCR (LMPCR) is a genomic sequencing method that can be used for the mapping of DNA damage at nucleotide resolution. Such a mapping can then be compared to a mutational spectrum to test the hypothesis that implies one agent can cause mutations into one cancer type. LMPCR has been used this way to map DNA damage generated by different UV wavelengths. The frequently damaged sites following UVB irradiation correlate with the mutational spectrum of p53 in skin cancer. Similarly, BPDE, the activated form of the benzo[a]pyrene present in tobacco smoke, generates frequent adducts at sites corresponding to mutation hotspots of p53 in lung cancers. Still, the correlation between BPDE damage sites and p53 mutations is not perfect and this suggests a role of other genotoxic substances that are also present in tobacco smoke, such as the nitrosamine NNK. Finally, and beyond this objective of better understanding somatic mutagenesis, LMPCR is commonly used whenever DNA damage frequency and/or repair is to be investigated.

Stationary phase mutagenesis in B. subtilis: a paradigm to study genetic diversity programs in cells under stress

One of the experimental platforms to study programs increasing genetic diversity in cells under stressful or nondividing conditions is adaptive mutagenesis, also called stationary phase mutagenesis or stress-induced mutagenesis. In some model systems, there is evidence that mutagenesis occurs in genes that are actively transcribed. Some of those genes may be actively transcribed as a result of environmental stress giving the appearance of directed mutation. That is, cells under conditions of starvation or other stresses accumulate mutations in transcribed genes, including those transcribed because of the selective pressure. An important question concerns how, within the context of stochastic processes, a cell biases mutation to genes under selection pressure? Because the mechanisms underlying DNA transactions in prokaryotic cells are well conserved among the three domains of life, these studies are likely to apply to the examination of genetic programs in eukaryotes. In eukaryotes, increasing genetic diversity in differentiated cells has been implicated in neoplasia and cell aging. Historically, Escherichia coli has been the paradigm used to discern the cellular processes driving the generation of adaptive mutations; however, examining adaptive mutation in Bacillus subtilis has contributed new insights. One noteworthy contribution is that the B. subtilis' ability to accumulate chromosomal mutations under conditions of starvation is influenced by cell differentiation and transcriptional derepression, as well as by proteins homologous to transcription and repair factors. Here we revise and discuss concepts pertaining to genetic programs that increase diversity in B. subtilis cells under nutritional stress.

G:C > A:T mutations and potential epigenetic regulation of p53 in breast cancer

Analysis of germline p53 mutations in breast cancer reveals that the Li-Fraumeni and Li-Fraumeni-like syndromes are mostly related to the loss of initiation codon 133 of regulatory TP53 isoforms (Delta133p53). In eight codons of exons 5-8 (including 133), mutations are frequent in Li-Fraumeni-related, but scarce in sporadic breast cancer, while in six more codons they are frequent both in familial and sporadic breast cancers. At the proximity of these codons, we observed in somatic mutation databases, 16 codons (minihotspots mostly in exons 7, 8) which undergo frequent G:C > A:T transitions (non-CpG) in all sporadic cancers. In addition, in sporadic breast cancer we observed 35 adjacent codons in which the following types of mutation are observed: frequent G:C > A:T transitions at CCs/GGs, frequent silent mutations in exons 5,6 and suppressed nonsense mutations (5 codons, few records). Non-CpG G:C > A:T transitions in the 35 codons are rare in familial cancers (p53, BRCA1, or BRCA2-related), but frequent in sporadic cancers in organs where Li-Fraumeni-related carcinogenesis is common e.g. adrenal cortex, soft tissues. These data are in support of the following tissue-specific processes: in sporadic breast cancer (sarcomas etc.), loss of methylation sites (in 35 codons mostly next to codon 133), might lead to loss of silencing of TP53 isoforms which are suppressed in these tissues. On the contrary, "spreading" of cytosine methylation (asymmetric) in a G:C-rich region next to common hotspots (codons 238-252 in minihotspots) and mutagenesis probably destabilizes all tissues. Frequent C > T activation at non-CpG is also observed in prostate sporadic cancer, which similarly to breast, undergoes age-related crisis. The above data reveal that tissue-specific epigenetic regulatory mechanisms might be involved in p53 instability.

Novel role of mfd: effects on stationary-phase mutagenesis in Bacillus subtilis

Previously, using a chromosomal reversion assay system, we established that an adaptive mutagenic process occurs in nongrowing Bacillus subtilis cells under stress, and we demonstrated that multiple mechanisms are involved in generating these mutations (41, 43). In an attempt to delineate how these mutations are generated, we began an investigation into whether or not transcription and transcription-associated proteins influence adaptive mutagenesis. In B. subtilis, the Mfd protein (transcription repair coupling factor) facilitates removal of RNA polymerase stalled at transcriptional blockages and recruitment of repair proteins to DNA lesions on the transcribed strand. Here we demonstrate that the loss of Mfd has a depressive effect on stationary-phase mutagenesis. An association between Mfd mutagenesis and aspects of transcription is discussed.

Predicting the oncogenicity of missense mutations reported in the International Agency for Cancer Research (IARC) mutation database on p53

Many mutation databases, comprising thousands of reported mutations, are available. Often the clinical significance of the reported mutations is unknown. In this study we developed an algorithm that allows prediction of the clinical significance of missense mutations reported in a mutation database. Nonsense mutations are used as a referent group for this assessment. We used the International Association for Research on Cancer (IARC) mutation database on TP53 to implement the algorithm. First, on the basis of published data [Nachman MW, Crowell SL. 2000. Genetics 156:297-304], we ascribed mutation rates to every single nucleotide substitution (SNS) in the core domain of the TP53 gene. Second, for every possible SNS we computed the expected number of missense mutations, under the assumption that missense mutations are as oncogenic as nonsense ones. The natural logarithm of the ratio of the observed to the expected number of missense mutations (LR) was used as a quantitative measure of oncogenicity (i.e., the ability of a mutation to produce cancer). We estimated the relative oncogenicity of all missense mutations reported in the IARC p53 mutation database, and constructed a profile of oncogenicity of the missense mutations along the DNA-binding region of p53.

Methylation and repeats in silent and nonsense mutations of p53

All exonic CG sequences in p53 are methylated; this epigenetic modification is correlated with frequent G:C-->A:T transitions in p53. Recent reports reveal the presence in p53 of non-CG methylation in CC and CCC sequences, complementary to sites of selective guanosine adduct formation (GG and GGG), and the association of genetic instability with methylation at repetitive sequences. We presently investigated the distribution of methylation sites and repetitive elements in silent and nonsense p53 mutations (2051) among the IARC's TP53 somatic mutation database for exons 5-8. Silent mutations are nonrandom, but mostly involve G:C-->A:T transitions (62%); in particular C-->T mutations (39% of all silent mutations) are mostly correlated with CC and CCC sequences, while G-->A mutations with GG sequences. Sequence analysis of all non-G:C-->A:T silent mutations reveals the frequent formation of new methylation sites (CG), new CCC and GGG sequences in the resulting sequence, refinement of symmetry elements at interrupted microsatellite-like sequences and formation of small repeats (55.3%). The G:C-->A:T silent mutations characterize cancers associated with cigarette smoking (e.g. bladder or lung and bronchus cancer versus colorectal cancer); on the contrary, non-G:C-->A:T silent mutations have similar frequencies in most cancers. Nonsense mutations in exons 5-8, all resulting in mutants lacking amino acids 307-393, which are crucial for p53 activity, were also analyzed. The frequency of nonsense mutations is higher at methylated sites or repeats 1-2 nucleotides removed from methylation sites. Frameshift mutations are also more frequent at repeated sequences. The frequent G:C-->A:T silent mutations could indicate that CC and CCC sequences of exons 5-8 are occasionally targets of non-CpG methylation of cytosine. This process of de novo methylation in the presence of microsatellite-like sequences and small repeats might influence the genetic stability of a variety of genes.

Repositioning-dependent fate of duplicate genes

Gene duplication is the main source of evolutionary novelties. However, the problem with duplicates is that the purifying selection overlooks deleterious mutations in the redundant sequence, which therefore, instead of gaining a new function, often degrades into a functionless pseudogene. This risk of functional loss instead of gain is much higher for small populations of higher organisms with a slow and complex development. We propose that it is the epigenetic tissue/stage-complementary silencing of duplicates that makes them exposable to the purifying selection, thus saving them from pseudogenization and opening the way towards new function(s). Our genome-wide analyses of gene duplicates in several eukaryotic species combined with the phylogenetic comparison of vertebrate alpha- and beta-globin gene clusters strongly support this epigenetic complementation (EC) model. The distinctive condition for a new duplicate to survive by the EC mechanism seems to be its repositioning to an ectopic site, which is accompanied by changes in the rate and direction of mutagenesis. The most distinguished in this respect is the human genome. In this review, we extend and discuss the data on the EC- and repositioning-dependent fate of gene duplicates with the special emphasis on the problem of detecting brief postduplication period of adaptive evolution driven by positive selection. Accordingly, we propose a new CpG-focused measure of selection that is insensitive to translocation-caused biases in mutagenesis.

Non-CpG cytosine methylation of p53 exon 5 in non-small cell lung carcinoma

Non-CpG methylation of cytosine residues, a mechanism associated with regulation of gene expression, has not been investigated in human cancer until now. Analysis of the p53 exon 5 mutation spectrum in mutation databases for lung cancer reveals frequent GC>AT transitions, several of which occur at non-CpG sequences. To investigate the involvement of cytosine methylation in this mutagenesis process, we analyzed the methylation profile of p53 exon 5, in lung carcinoma. In this report, we present evidence that extensive clustered non-CpG methylation is observed in three regions of this exon, namely the sequences spanning codons 156-159, 175-179 and the 3' splice site, as well as in scattered CpA sequences. This methylation pattern was verified using direct methylation sequencing, and a two-stage methylation-specific PCR assay (MSP), designed for the detection of methylation in a GC rich region (oligo C sequence, of codons 175-179) of exon 5. The results from this MSP assay reveal that DNA from cancerous specimens was more heavily methylated in non-CpG cytosines, compared to that from non-cancerous lung tissue of cancer patients (14/19 cancerous and 6/19 non-cancerous, respectively). DNA isolated from human leucocytes and some non-cancerous specimens (2/19) was free of non-CpG methylation. Careful analysis of the mutations reported in p53 mutation databases also provides corroborating evidence that the high incidence of GC>AT mutations in the p53 gene, observed in lung cancer, might also be related to non-CpG methylation, as well as to the overall increase of methylation sites in this locus.

Clusters of mutations from transient hypermutability

Collections of mutants usually contain more mutants bearing multiple mutations than expected from the mutant frequency and a random distribution of mutations. This excess is seen in a variety of organisms and also after DNA synthesis in vitro. The excess is unlikely to originate in mutator mutants but rather from transient hypermutability resulting from a perturbation of one of the many transactions that maintain genetic fidelity. The multiple mutations are sometimes clustered and sometimes randomly distributed. We model some spectra as populations comprising a majority with a low mutation frequency and a minority with a high mutation frequency. In the case of mutants produced in vitro by a bacteriophage RB69 mutator DNA polymerase, mutants with two mutations are in approximately 10-fold excess and mutants with three mutations are in even greater excess. However, phenotypically undetectable mutations seen only as hitchhikers with detectable mutations are approximately 5-fold more frequent than mutants bearing detectable mutations, indicating that they arose in a subpopulation with a higher mutation frequency. Excess multiple mutations may contribute critically to carcinogenesis and to adaptive mutation, including the adaptations of pathogens as they move from host to host. In the case of the rapidly mutating riboviruses, the viral population appears to be composed of a majority with a mutation frequency substantially lower than the average and a minority with a huge mutational load.

Position-associated GC asymmetry of gene duplicates

It is well known that repositioning of a gene often exerts a strong impact on its own expression and whole development. Here we report the results of genome-wide analyses suggesting that repositioning may also radically change the evolutionary fate of gene duplicates. As an indicator of these changes, we used the GC content of gene pairs which originated by duplication. This indicator turned out to be duplicate-asymmetric, which means that genes in a pair differ significantly in GC content despite their apparent origin from a common ancestor. Such an asymmetry necessarily implies that after duplication two originally identical genes mutated in opposite directions-toward GC-rich and GC-poor content, respectively. In mammalian genomes, this trend is definitely associated with presumably methylated hypermutable CpG sites, and in a typical GC-asymmetric gene pair, its two member genes are embedded in GC-contrasting isochores. However, we unexpectedly found similar significant GC asymmetry in fish, fly, worm, and yeast. This means that neither methylation alone nor methylation in combination with isochores can be counted as a primary cause of the GC asymmetry; rather they represent specific realizations of some universal principle of genome evolution. Remarkably, genes from pairs with the greatest GC asymmetry tend to be on different chromosomes, suggesting that the mutational difference between gene duplicates is associated with translocation of a new gene to a different place in the genome, whereas GC symmetric pairs demonstrate the opposite tendency. A recently emerged extra gene copy is usually on the same chromosome as is its parent but quickly, by 0.05 substitution per synonymous site, either has perished or occupies a different chromosome. During this earliest posttranslocation period, the ratio of nonsynonymous/synonymous base substitutions is unusually high, suggesting a rapid adaptive evolution of novel functions. In a general context of evolution by gene duplication, our interpretation of this position-dependent GC asymmetry between duplicated genes is that evolution of redundant genes toward a new function has often been associated with their very early, postduplication repositioning in the genome, with a concomitant abrupt change in epigenetic control of tissue/stage-specific expression and an increase in the mutation rate. Of eight eukaryotic genomes studied, the most distinguished in this respect is the human genome.

Genotoxicity of environmental tobacco smoke: a review

Environmental tobacco smoke (ETS), or second-hand smoke, is a widespread contaminant of indoor air in environments where smoking is not prohibited. It is a significant source of exposure to a large number of substances known to be hazardous to human health. Numerous expert panels have concluded that there is sufficient evidence to classify involuntary smoking (or passive smoking) as carcinogenic to humans. According to the recent evaluation by the International Agency for Research on Cancer, involuntary smoking causes lung cancer in never-smokers with an excess risk in the order of 20% for women and 30% for men. The present paper reviews studies on genotoxicity and related endpoints carried out on ETS since the mid-1980s. The evidence from in vitro studies demonstrates induction of DNA strand breaks, formation of DNA adducts, mutagenicity in bacterial assays and cytogenetic effects. In vivo experiments in rodents have shown that exposure to tobacco smoke, whole-body exposure to mainstream smoke (MS), sidestream smoke (SS), or their mixture, causes DNA single strand breaks, aromatic adducts and oxidative damage to DNA, chromosome aberrations and micronuclei. Genotoxicity of transplacental exposure to ETS has also been reported. Review of human biomarker studies conducted among non-smokers with involuntary exposure to tobacco smoke indicates presence of DNA adducts, urinary metabolites of carcinogens, urinary mutagenicity, SCEs and hypoxanthine-guanine phosphoribosyltransferase (HPRT) gene mutations (in newborns exposed through involuntary smoking of the mother). Studies on human lung cancer from smokers and never-smokers involuntarily exposed to tobacco smoke suggest occurrence of similar kinds of genetic alterations in both groups. In conclusion, these overwhelming data are compatible with the current knowledge on the mechanisms of carcinogenesis of tobacco-related cancers, occurring not only in smokers but with a high biological plausibility also in involuntary smokers.

Missense mutations in cancer suppressor gene TP53 are colocalized with exonic splicing enhancers (ESEs)

Mutation databases can be viewed as footprints of functional organization of a gene and thus can be used to infer its functional organization. We studied the association of exonic splicing enhancers (ESEs) with missense mutations in the tumor suppressor gene TP53 using the International Agency for Research on Cancer (IARC) mutation database. The goals of the study were: (i) to verify the hypothesis that deleterious missense mutations are colocalized with ESEs; (ii) to identify potentially functional ESE sites in the open reading frame (ORF) of the TP53. If some sequence functions as a splicing enhancer, then nucleotide substitutions in the site will disturb splicing, abrogate p53 function, and cause an increased susceptibility to cancer. Therefore, among cancers showing p53 mutations, more missense mutations are expected within functional ESE sites as compared to non-functional ESE motifs. Using several statistical tests, we found that missense mutations in TP53 are strongly colocalized with ESEs, and that only a small fraction of ESE sites contributes to the association. There are usually one or two ESEs per exon showing a statistically significant association with missense mutations--so-called significant ESE sites. In many respects significant ESE sites are different from those that do not show association with missense mutations. We found that positions of significant ESE sites are codon-dependent--significant ESEs preferentially start from the first position of a codon, whereas non-significant ESEs show no position dependence. Significant ESEs showed a more limited set of sequences compared to non-significant ESEs. These findings suggest that there is a limited number of missense mutations that influence ESE sites and our analysis provides further insight into the types of sites that harbor exonic enhancer elements.