Molecular analysis of mutagenesis in mammalian cells

Mammalian cells are constantly facing various types of mutagens. However, due to the high complexity of the cell genome, the molecular analysis of mutagenesis has not yet been possible. Therefore, we have used simian virus 40 (SV40) as a biological and molecular probe to characterize mutagenesis at the nucleotide level. By using a reversion assay from a temperature-sensitive phenotype towards a wild-type phenotype, we have analysed mutagenesis induced by u.v.-light and by apurinic sites (Ap sites). We report here experiments allowing us to quantify and to compare the mutagenic efficiency of various DNA lesions measured on the SV40 genome. The Ap sites are very mutagenic in this type of assay. The molecular analysis of u.v.-induced mutagenesis reveals that mutations correspond to single base-pair substitutions always located opposite Py-Py lesions. The mutations are almost equally distributed between transition and transversion types, and between the 5' and the 3' side of the Py-Py targets. These results demonstrate for the first time in animal cells the existence of targeted mutations induced by u.v.-light. We propose therefore, the use of SV40 as an efficient biological and molecular probe for assaying mutagenic pathways in mammalian cells.     

Co-evolutionary analysis of domains in interacting proteins reveals insights into domain-domain interactions mediating protein-protein interactions

Recent advances in functional genomics have helped generate large-scale high-throughput protein interaction data. Such networks, though extremely valuable towards molecular level understanding of cells, do not provide any direct information about the regions (domains) in the proteins that mediate the interaction. Here, we performed co-evolutionary analysis of domains in interacting proteins in order to understand the degree of co-evolution of interacting and non-interacting domains. Using a combination of sequence and structural analysis, we analyzed protein-protein interactions in F1-ATPase, Sec23p/Sec24p, DNA-directed RNA polymerase and nuclear pore complexes, and found that interacting domain pair(s) for a given interaction exhibits higher level of co-evolution than the non-interacting domain pairs. Motivated by this finding, we developed a computational method to test the generality of the observed trend, and to predict large-scale domain-domain interactions. Given a protein-protein interaction, the proposed method predicts the domain pair(s) that is most likely to mediate the protein interaction. We applied this method on the yeast interactome to predict domain-domain interactions, and used known domain-domain interactions found in PDB crystal structures to validate our predictions. Our results show that the prediction accuracy of the proposed method is statistically significant. Comparison of our prediction results with those from two other methods reveals that only a fraction of predictions are shared by all the three methods, indicating that the proposed method can detect known interactions missed by other methods. We believe that the proposed method can be used with other methods to help identify previously unrecognized domain-domain interactions on a genome scale, and could potentially help reduce the search space for identifying interaction sites.     

Transition-transversion bias is not universal: a counter example from grasshopper pseudogenes

Comparisons of the DNA sequences of metazoa show an excess of transitional over transversional substitutions. Part of this bias is due to the relatively high rate of mutation of methylated cytosines to thymine. Postmutation processes also introduce a bias, particularly selection for codon-usage bias in coding regions. It is generally assumed, however, that there is a universal bias in favour of transitions over transversions, possibly as a result of the underlying chemistry of mutation. Surprisingly, this underlying trend has been evaluated only in two types of metazoan, namely Drosophila and the Mammalia. Here, we investigate a third group, and find no such bias. We characterize the point substitution spectrum in Podisma pedestris, a grasshopper species with a very large genome. The accumulation of mutations was surveyed in two pseudogene families, nuclear mitochondrial and ribosomal DNA sequences. The cytosine-guanine (CpG) dinucleotides exhibit the high transition frequencies expected of methylated sites. The transition rate at other cytosine residues is significantly lower. After accounting for this methylation effect, there is no significant difference between transition and transversion rates. These results contrast with reports from other taxa and lead us to reject the hypothesis of a universal transition/transversion bias. Instead we suggest fundamental interspecific differences in point substitution processes.     

Specificity of Escherichia coli mutD and mutL mutator strains

The products of the mutD and mutL genes of Escherichia coli are involved in proofreading by DNA polymerase III and DNA adenine MTase (Dam)-dependent mismatch repair, respectively. We have used the plasmid-borne bacteriophage P22 mnt gene as a target to determine the types of mutations produced in mutL25 and mutD5 strains. Of 60 mutations identified from mutL25 cells, 52 were transition mutations and of these the AT----GC subset predominated (40 out of 52). The majority of AT----GC mutations were found at the same three sites (hotspots). In contrast, transversion mutations (47 out of 76) were found about twice as frequently as transitions (28 out of 76) from mutD5 bacteria. Two hotspots were identified but at different sites than those in the mutL25 cells. These results suggest that the proofreading function of DNA polymerase III primarily repairs potential transversion mutations while Dam-dependent mismatch repair rectifies potential transition mutations.     

蕨类植物rpoC1内含子缺失及其分子进化速率

rpoC1基因编码RNA聚合酶β°亚基蛋白, 在转录过程中与DNA模板结合, 与β亚基形成的β-β°亚基复合体构成RNA合成的催化中心。以rpoC1基因为研究对象, 在贝叶斯因子大于20的条件下, 用HyPhy软件位点模型检测到3个正选择位点和541个负选择位点; 用PAML软件位点模型检测到10个正选择位点, 其中3个位点的后验概率超过99%。此外, 基于最大似然法构建64种蕨类植物的系统发育树, 结合HyPhy软件分析rpoC1基因的转换率、颠换率、转换率/颠换率、同义替换率、非同义替换率以及同义替换率/非同义替换率, 探讨rpoC1基因内含子丢失与分子进化速率的关系。结果表明, rpoC1基因内含子缺失对转换率、颠换率以及非同义替换率有一定影响。     

Mutational spectrum in the recent human genome inferred by single nucleotide polymorphisms

So far, there is no genome-wide estimation of the mutational spectrum in humans. In this study, we systematically examined the directionality of the point mutations and maintenance of GC content in the human genome using approximately 1.8 million high-quality human single nucleotide polymorphisms and their ancestral sequences in chimpanzees. The frequency of C-->T (G-->A) changes was the highest among all mutation types and the frequency of each type of transition was approximately fourfold that of each type of transversion. In intergenic regions, when the GC content increased, the frequency of changes from G or C increased. In exons, the frequency of G:C-->A:T was the highest among the genomic categories and contributed mainly by the frequent mutations at the CpG sites. In contrast, mutations at the CpG sites, or CpG-->TpG/CpA mutations, occurred less frequently in the CpG islands relative to intergenic regions with similar GC content. Our results suggest that the GC content is overall not in equilibrium in the human genome, with a trend toward shifting the human genome to be AT rich and shifting the GC content of a region to approach the genome average. Our results, which differ from previous estimates based on limited loci or on the rodent lineage, provide the first representative and reliable mutational spectrum in the recent human genome and categorized genomic regions.     

The complementary neighborhood patterns and methylation-to-mutation likelihood structures of 15,110 single-nucleotide polymorphisms in the bovine genome

Bayesian analysis was performed to examine the single-nucleotide polymorphism (SNPs) neighborhood patterns in cattle using 15,110 SNPs, each with a flanking sequence of 500 bp. Our analysis confirmed three well-known features reported in plants and/or other animals: (1) the transition is the most abundant type of SNPs, accounting for 69.8% in cattle; (2) the transversion occurs most frequently (38.56%) in cattle when the A + T content equals two at their immediate adjacent sites; and (3) C <--> T and A <--> G transitions have reverse complementary neighborhood patterns and so do A <--> C and G <--> T transversions. Our study also revealed several novel SNP neighborhood patterns that have not been reported previously. First, cattle and humans share an overall SNP pattern, indicating a common mutation system in mammals. Second, unlike C <--> T/A <--> G and A <--> C/G <--> T, the true neighborhood patterns for A <--> T and C <--> G might remain mysterious because the sense and antisense sequences flanking these mutations are not actually recognizable. Third, among the reclassified four types of SNPs, the neighborhood ratio between A + T and G + C was quite different. The ratio was lowest for C <--> G, but increased for C <--> T/A <--> G, further for A <--> C/G <--> T, and the most for A <--> T. Fourth, when two immediate adjacent sites provide structures for CpG, it significantly increased transitions compared to the structures without the CpG. Finally, unequal occurrence between A <--> G and C <--> T in five paired neighboring structures indicates that the methylation-induced deamination reactions were responsible for approximately 20% of total transitions. In addition, conversion can occur at both CpG sites and non-CpG sites. Our study provides new insights into understanding molecular mechanisms of mutations and genome evolution.     

The oncoprotein Bcl-3 can facilitate NF-kappa B-mediated transactivation by removing inhibiting p50 homodimers from select kappa B sites.

Previously we have proposed a role for Bcl-3 in facilitating transactivation through kappa B sites by counteracting the inhibitory effects of bound, non-transactivating homodimers of the p50 subunit of NF-kappa B. Such homodimers are abundant for example in nuclei of unstimulated primary T cells. Here we extend the model and provide new evidence which fulfills a number of predictions. (i) Bcl-3 preferentially targets p50 homodimers over NF-kappa B heterodimers since the homodimers are completely dissociated from kappa B sites at concentrations of Bcl-3 which do not affect NF-kappa B. (ii) Select kappa B sites associate very strongly and stably with p50 homodimers, completely preventing binding by NF-kappa B. Such kappa B sites are likely candidates for regulation by p50 homodimers and Bcl-3. (iii) Bcl-3 and p50 can be co-localized in the nucleus, a requirement for active removal of homodimers from their binding sites in vivo. (iv) The ankyrin repeat domain of Bcl-3 is sufficient for the reversal of p50 homodimer-mediated inhibition, correlating with the ability of this domain alone to inhibit p50 binding to kappa B sites in vitro. Our data support the model that induction of nuclear Bcl-3 may be required during cellular stimulation to actively remove stably bound p50 homodimers from certain kappa B sites in order to allow transactivating NF-kappa B complexes to engage. This exact mechanism is demonstrated with in vitro experiments.     

Repair of a Mismatch Is Influenced by the Base Composition of the Surrounding Nucleotide Sequence

Heteroduplexes with single base pair mismatches of known sequence were prepared by annealing separated strands of bacteriophage lambda DNA and used to transfect Escherichia coli. A series of transition (G:T and A:C) and transversion (G:A and C:T) mismatches located throughout most of the bacteriophage lambda cI gene has been examined. The results suggest that the transition mismatches are generally better repaired than the transversion mismatches and that, at least for the transversion mismatches studied, repair efficiency increases with increasing G:C content in the neighboring nucleotide sequence. This specificity of the E. coli mismatch repair system can account, in part, for the similar frequencies of base substitution mutations throughout the E. coli genome.     

Influence of uvrB and pKM101 on the spectrum of spontaneous, UV- and gamma-ray-induced base substitutions that revert hisG46 in Salmonella typhimurium

Oligonucleotide probes were used to identify base substitutions in 1089 revertants of hisG46 in Salmonella typhimurium that arose spontaneously or following irradiation with UV- or gamma-rays. The hisG46 allele, carrying a mutant CCC codon (Pro) in place of the wild-type codon CTC (Leu69) reverted via 6 distinguishable mutational events--C to T transitions at codon sites 1 or 2, C to A or C to G transversions at codon site 1, C to A at codon site 2, and an extragenic suppressor mutation. The distribution of hisG46 revertants differed among treatments and was influenced by the DNA-repair capacity of the bacteria. Plasmid pKM101 enhanced the frequencies of both spontaneous and induced mutations; transversion events were enhanced more efficiently by pKM101 than were transition events. Compared to Uvr+ bacteria, Uvr- bacteria had higher frequencies of spontaneous and induced mutations; transition mutations were enhanced more efficiently than were transversion mutations. The influence of DNA-repair activities on the mutational spectra provides some insights on the origins of spontaneous and UV-induced mutations.     

Thermodynamic and phylogenetic prediction of RNA secondary structures in the coding region of hepatitis C virus

The existence and functional importance of RNA secondary structure in the replication of positive-stranded RNA viruses is increasingly recognized. We applied several computational methods to detect RNA secondary structure in the coding region of hepatitis C virus (HCV), including thermodynamic prediction, calculation of free energy on folding, and a newly developed method to scan sequences for covariant sites and associated secondary structures using a parsimony-based algorithm. Each of the prediction methods provided evidence for complex RNA folding in the core- and NS5B-encoding regions of the genome. The positioning of covariant sites and associated predicted stem-loop structures coincided with thermodynamic predictions of RNA base pairing, and localized precisely in parts of the genome with marked suppression of variability at synonymous sites. Combined, there was evidence for a total of six evolutionarily conserved stem-loop structures in the NS5B-encoding region and two in the core gene. The virus most closely related to HCV, GB virus-B (GBV-B) also showed evidence for similar internal base pairing in its coding region, although predictions of secondary structures were limited by the absence of comparative sequence data for this virus. While the role(s) of stem-loops in the coding region of HCV and GBV-B are currently unknown, the structure predictions in this study could provide the starting point for functional investigations using recently developed self-replicating clones of HCV.     

Use of supF, an Escherichia coli tyrosine suppressor tRNA gene, as a mutagenic target in shuttle-vector plasmids

The Escherichia coli tyrosine amber suppressor tRNA gene, supF, has been utilized as a mutagenic target in several shuttle-vector plasmids. Data on mutagenic inactivation of suppressor activity was obtained from induced mutagenesis experiments with plasmids pZ189 and p3AC, and from studies on alterations of the supF gene transduced into E. coli. 162 single or tandem base-substitution mutations that reduce or eliminate suppressor activity were identified at 86 sites within 158 base pairs. The 2 transition and 4 transversion mutations possible in double-stranded DNA were all detectable. At 56 sites two different inactivating mutations were found; and at 20 sites all 3 possible base substitution mutations inactivated suppressor function. Most of the mutations were clustered within the mature tRNA region: 144 of the base-substitution mutations were found at 74 sites within the 85-bp mature tRNA region. Insertions of 1 or 2 bases at 4 sites and deletions of 1 to 3 bases at 15 sites were found to inactivate supF function. A few silent mutations which do not inactivate suppressor function were found: single base-substitutions at 4 sites, 14 pairs of silent double mutations, and a large deletion including the promoter region. The supF gene is thus an extremely sensitive target for mutagenic inactivation in shuttle-vector plasmids.     

Computational analysis of hydrogen bonds in protein-RNA complexes for interaction patterns

Structural analysis of protein-RNA complexes is labor-intensive, yet provides insight into the interaction patterns between a protein and RNA. As the number of protein-RNA complex structures reported has increased substantially in the last few years, a systematic method is required for automatically identifying interaction patterns. This paper presents a computational analysis of the hydrogen bonds in the most representative set of protein-RNA complexes. The analysis revealed several interesting interaction patterns. (1) While residues in the beta-sheets favored unpaired nucleotides, residues in the helices showed no preference and residues in turns favored paired nucleotides. (2) The backbone hydrogen bonds were more dominant than the base hydrogen bonds in the paired nucleotides, but the reverse was observed in the unpaired nucleotides. (3) The protein-RNA complexes contained more paired nucleotides than unpaired nucleotides, but the unpaired nucleotides were observed more frequently interacting with the proteins. And (4) Arg-U, Thr-A, Lys-A, and Asn-U were the most frequently observed pairs. The interaction patterns discovered from the analysis will provide us with useful information in predicting the structure of the RNA binding protein and the structure of the protein binding RNA.     

Variation in the Pattern of Nucleotide Substitution Across Sites

A model of nucleotide substitution that allows the transition/transversion rate bias to vary across sites was constructed. We examined the fit of this model using likelihood-ratio tests by analyzing 13 protein coding genes and 1 pseudogene. Likelihood-ratio testing indicated that a model that allows variation in the transition/transversion rate bias across sites provided a significant improvement in fit for most protein coding genes but not for the pseudogene. When the analysis was repeated with parameters estimated separately for first, second, and third codon positions, strong heterogeneity was uncovered for the first and second codon positions; the variation in the transition/transversion rate was generally weaker at the third codon position. The transition rate bias and branch lengths are underestimated when variation in the transition/transversion rate was not accommodated, suggesting that it may be important to accommodate variation in the pattern of nucleotide substitution for accurate estimation of evolutionary parameters. Received: 4 November 1997 / Accepted: 19 May 1998     

Heterogeneity of the transition/transversion ratio in Drosophila and Hominidae genomes

Mutation rate varies between sites in the genome. Part of this variation can be explained by well-recognized short nucleotide contexts, but a large component of this variation remains cryptic. We used data on interspecies divergence and intraspecies polymorphism in Drosophila and Hominidae to analyze variation of the average rate of the 12 possible kinds of single-nucleotide mutations and in the transition/transversion ratio κ at single-nucleotide resolution. Both the average mutation rate and κ vary by a factor of ~3 between nucleotide sites. The characteristic scale of variation in κ is up to at least ~30 nucleotides in Drosophila and ~5 nucleotides in Hominidae. Genome segments with locally elevated mutation rates possess lower values of κ; however, a substantial fraction of variation in κ cannot be directly explained by the local mutation rates.     

p53 tumor suppressor gene mutations in fibroblast-like synoviocytes from erosion synovium and non-erosion synovium in rheumatoid arthritis

Abnormalities in the p53 tumor suppressor gene have been detected in rheumatoid arthritis (RA) and could contribute to the pathogenesis of chronic disease. To determine whether synoviocytes from invasive synovium in RA have an increased number of mutations compared with non-erosion synoviocytes, p53 cDNA subclones from fibroblast-like synoviocytes (FLS) derived from erosion and non-erosion sites of the same synovium were examined in patients requiring total joint replacement. Ten erosion FLS lines and nine non-erosion FLS lines were established from nine patients with RA. Exons 5-10 from 209 p53 subclones were sequenced (114 from erosion FLS, 95 from non-erosion FLS). Sixty percent of RA FLS cell lines and 8.6% of the p53 subclones isolated from FLS contained p53 mutations. No significant differences were observed between the erosion and non-erosion FLS with regard to the frequency or type of p53 mutation. The majority of the mutations were missense transition mutations, which are characteristic of oxidative damage. In addition, paired intact RA synovium and cultured FLS from the same joints were evaluated for p53 mutations. Matched synovium and cultured synoviocytes contained p53 mutations, although there was no overlap in the specific mutations identified in the paired samples. Clusters of p53 mutations in subclones were detected in some FLS, including one in codon 249, which is a well-recognized 'hot spot' associated with cancer. Our data are consistent with the hypothesis that p53 mutations are randomly induced by genotoxic exposure in small numbers of RA synoviocytes localized to erosion and non-erosion regions of RA synovium. The determining factor for invasiveness might be proximity to bone or cartilage rather than the presence of a p53 mutation.