Mutation rates of TGFBR2 and ACVR2 coding microsatellites in human cells with defective DNA mismatch repair

Microsatellite instability promotes colonic tumorigenesis through generating frameshift mutations at coding microsatellites of tumor suppressor genes, such as TGFBR2 and ACVR2. As a consequence, signaling through these TGFbeta family receptors is abrogated in DNA Mismatch repair (MMR)-deficient tumors. How these mutations occur in real time and mutational rates of these human coding sequences have not previously been studied. We utilized cell lines with different MMR deficiencies (hMLH1-/-, hMSH6-/-, hMSH3-/-, and MMR-proficient) to determine mutation rates. Plasmids were constructed in which exon 3 of TGFBR2 and exon 10 of ACVR2 were cloned +1 bp out of frame, immediately after the translation initiation codon of an enhanced GFP (EGFP) gene, allowing a -1 bp frameshift mutation to drive EGFP expression. Mutation-resistant plasmids were constructed by interrupting the coding microsatellite sequences, preventing frameshift mutation. Stable cell lines were established containing portions of TGFBR2 and ACVR2, and nonfluorescent cells were sorted, cultured for 7-35 days, and harvested for flow cytometric mutation detection and DNA sequencing at specific time points. DNA sequencing revealed a -1 bp frameshift mutation (A9 in TGFBR2 and A7 in ACVR2) in the fluorescent cells. Two distinct fluorescent populations, M1 (dim, representing heteroduplexes) and M2 (bright, representing full mutants) were identified, with the M2 fraction accumulating over time. hMLH1 deficiency revealed 11 (5.91 x 10(-4)) and 15 (2.18 x 10(-4)) times higher mutation rates for the TGFBR2 and ACVR2 microsatellites compared to hMSH6 deficiency, respectively. The mutation rate of the TGFBR2 microsatellite was approximately 3 times higher in both hMLH1 and hMSH6 deficiencies than the ACVR2 microsatellite. The -1 bp frameshift mutation rates of TGFBR2 and ACVR2 microsatellite sequences are dependent upon the human MMR background.     

NMR evidence of the stabilisation by the carcinogen N-2-acetylaminofluorene of a frameshift mutagenesis intermediate.

Two heteroduplexes d(C1A2C3T4C5G6C7A8C9A10C11)-d (G12T13G14T15G16G17A18G19T20G21) containing a bulged guanine either unmodified or modified with the carcinogen N-2-acetylaminofluorene (AAF) have been studied by nuclear magnetic resonance (NMR) as models of slipped mutagenic intermediates (SMI). Conformational equilibria are observed in both the unmodified and the AAF-modified heteroduplexes. The major conformation of the unmodified duplex is one where the extra guanine is stacked in the helix and the major conformation of the AAF-modified heteroduplex is one where the AAF is external to the helix. Unusual sugar proton chemical shifts of C5- and G6-AAF indicate that the AAF ring is pointing out in the 5' direction. A strong increase in the modified heteroduplex melting temperature (+15 degrees C) is observed. Moreover, in contrast to the unmodified heteroduplex, which shows extensive melting in the vicinity of the bulged guanine, the base pairs around the bulge in the AAF-modified heteroduplex remain paired at temperatures up to 30 degrees C. This exceptional stability of the site around the bulged modified guanine is suggested to be responsible for the high rate of -1 mutation induced by AAF at repetitive sequences.     

Effects of phosphate deficiency and sugars on expression of rab18 in Arabidopsis: hexokinase-dependent and okadaic acid-sensitive transduction of the sugar signal

Analbuminemia is a rare autosomal recessive disorder manifested by the absence or severe reduction of circulating serum albumin in homozygous subjects. In this report we describe a new molecular defect that caused the analbuminemic trait in a newborn of Iraqi origin. When the parents' DNA was analyzed, both subjects were found to be heterozygous for the same mutation found in the infant. All the 14 exon and flanking intron sequences of the albumin gene were amplified via PCR and screened for mutations by SSCP and heteroduplex analysis. A mutation in the DNA region encoding exon 1 and its flanking intron was revealed by the presence of a heteroduplex. The fragment, which was directly DNA sequenced, contains a previously unreported single nucleotide change, consisting in a G to A substitution at nucleotide 118 in the structural gene of the human protein. This mutation, involving the first base of intron 1, destroys the GT dinucleotide consensus sequence found at the 5' end of most intervening sequences and causes the defective pre-mRNA splicing responsible for the analbuminemic trait.     

Sequence-dependent modulation of frameshift mutagenesis at NarI-derived mutation hot spots.

The NarI sequence is known to be the strongest mutation hot spot for induced frameshift mutagenesis. Indeed, a single N-2-acetylaminofluorene (AAF) adduct induces -2 frameshift mutations (5'-GGCGAAFCC--> 5'-GGCC) more than 10(7)-fold over background mutagenesis in Escherichia coli. The mechanism of induction of the frameshift mutation involves a two nucleotide primer-template misalignment event during replication of the adduct-containing sequence. The slipped mutagenic intermediate (SMI) that is thus formed is strongly stabilised by the AAF residue. In order to understand the origin of the extreme susceptibility of this sequence to frameshift mutagenesis, we analysed AAF-induced mutagenesis at sequences 5'-NaGCGAAFCNb-3' containing the core dinucleotide GCGC repeat present in the NarI sequence flanked by variable nucleotides Na and Nb. The nature of nucleotide Nb was found to strongly modulate the frequency of induced -2 frameshift mutagenesis (up to 30 to 50-fold), while little if any effect could be attributed to nucleotide Na. The induction of -2 frameshifts, regardless of nucleotides Na and Nb, was found to be SOS-inducible but umuDC-independent as previously found for the authentic NarI sequence. The NarI sequence (GGCGCC) and sequence TGCGCA (Na=T, Nb=A) were found to be equally "hot" for -2 frameshift mutation induction compared to the sequence AGCGCT where induced mutagenesis was 30 to 50-fold lower.The analysis of replication events using constructions containing a strand marker across from the adduct site allowed us to demonstrate that the large difference in -2 frameshift mutagenesis is due to an intrinsic difference in the propensity of these sequences to slip during replication. How the nature of the nucleotide flanking the adduct on its 3'-side (Nb) differentially stabilises the SMI will be discussed in the light of recent structural data and theoretical models.     

Isolation and characterization of the human apolipoprotein A-II gene. Electron microscopic analysis of RNA:DNA hybrids, nucleotide sequence, identification of a polymorphic MspI site, and general structural organization of apolipoprotein genes

Three cloned apolipoprotein A-II genes were isolated from a human genomic cosmid library constructed in our laboratory. An approximately 3-kilobase HindIII insert containing the entire gene was analyzed by RNA:DNA hybridization and electron microscopy. The apo-A-II gene was found to consist of 4 exons and 3 intervening sequences (IVS), and the lengths of each exon and IVS were estimated by direct observation of the hybrids. The entire approximately 3-kilobase HindIII insert was sequenced. The 5' end of the gene was determined by primer extension. The DNA sequence confirms the presence of 4 exons and 3 IVS: exon 1, 34 nucleotides; exon 2, 76 nucleotides; exon 3, 133 nucleotides; exon 4, 230 nucleotides; IVS-I, 169 nucleotides; IVS-II, 299 nucleotides; and IVS-III, 396 nucleotides. A "TATA box" is located at position -29 from the CAP site. A "CAT box" is present at position -78. A "TG" element consisting of (TG)19 is identified at the 3' end of IVS-III. Furthermore, an enhancer core sequence, CTTTCCA, is identified at position -355 in the 5' flanking sequence. At positions -497 to -471 upstream from the CAP site is a stretch of 27 nucleotides that show high homology to stretches of 5' flanking sequences in the apo-C-II, apo-A-I, apo-E, and apo-C-III genes. An Alu dimer sequence is located approximately 300 nucleotides from the 3' end of the gene. Within this Alu sequence, we have identified a polymorphic MspI site. Restriction fragment length polymorphism involving this site has been previously shown to correlate with apo-A-II levels and high density lipoprotein structure. Analysis of conformation by Chou-Fasman analysis and by the helical hydrophobic moment of Eisenberg et al. (Eisenberg, D., Weiss, R. M., and Tergwillager, T. C. (1982). Nature (Lond.) 299, 371-374) indicates that in all of the 5 apolipoproteins characterized at the nucleotide level to date, i.e. apo-C-II, apo-A-II, apo-E, apo-A-I, and apo-C-III, the 2 IVS within the peptide coding regions of the gene tend to occur at regions corresponding to the surface of the polypeptide chain and divide the protein into distinct functional domains.     

Conserved signals in the 5' flanking region of eukaryotic nuclear tRNA genes.

The statistical analysis of 5' flanking regions of eukaryotic tRNA genes was done. The analysis of nucleotides in the sequence of fungi and invertebrates showed a high content of A and T in the flanking regions versus coding regions where G and C dominate. In contrast to these results in vertebrates sequences the preferences of any nucleotide in flanking regions was not observed. The analysis of tetrads showed five conserved signals: TTGT, (T/A)(T/A)ATA, A(C/T)(C/A)A in the tRNA genes of fungi, (A/T)TGA of invertebrates and (A/T)GAG of vertebrates. The analysis of 3' flanking regions did not show any conserved signals except well known poly-T tracks.     

A novel exon mutation in the human beta-hexosaminidase beta subunit gene affects 3' splice site selection.

The molecular basis of a dramatically decreased steady state level of beta-hexosaminidase beta subunit mRNA in a patient with juvenile Sandhoff disease was investigated. Nucleotide sequence analysis of the HEXB gene coding for the beta subunit revealed two single base substitutions, one in exon 2 (A to G, a known polymorphism) and the other in exon 11 (C to T). Analysis of the beta subunit mRNA species demonstrated activation of a cryptic splice site in exon 11 as well as skipping of the exon. A transfection assay using a chimeric gene containing intron 10 flanked by cDNA sequences carrying the mutation confirmed that the single base substitution located at position 8 of exon 11 inhibits the selection of the normal 3' splice site. The results demonstrate a new type of exon mutation affecting 3' splice site selection.     

Mutation screening of the EXT1 and EXT2 genes in patients with hereditary multiple exostoses.

Hereditary multiple exostoses (HME), the most frequent of all skeletal dysplasias, is an autosomal dominant disorder characterized by the presence of multiple exostoses localized mainly at the end of long bones. HME is genetically heterogeneous, with at least three loci, on 8q24.1 (EXT1), 11p11-p13 (EXT2), and 19p (EXT3). Both the EXT1 and EXT2 genes have been cloned recently and define a new family of potential tumor suppressor genes. This is the first study in which mutation screening has been performed for both the EXT1 and EXT2 genes prior to any linkage analysis. We have screened 17 probands with the HME phenotype, for alterations in all translated exons and flanking intronic sequences, in the EXT1 and EXT2 genes, by conformation-sensitive gel electrophoresis. We found the disease-causing mutation in 12 families (70%), 7 (41%) of which have EXT1 mutations and 5 (29%) EXT2 mutations. Together with the previously described 1-bp deletion in exon 6, which is present in 2 of our families, we report five new mutations in EXT1. Two are missense mutations in exon 2 (G339D and R340C), and the other three alterations (a nonsense mutation, a frameshift, and a splicing mutation) are likely to result in truncated nonfunctional proteins. Four new mutations are described in EXT2. A missense mutation (D227N) was found in 2 different families; the other three alterations (two nonsense mutations and one frameshift mutation) lead directly or indirectly to premature stop codons. The missense mutations in EXT1 and EXT2 may pinpoint crucial domains in both proteins and therefore give clues for the understanding of the pathophysiology of this skeletal disorder.     

Nucleotide sequence and genomic organization of a human T lymphocyte serine protease gene

We have determined the nucleotide sequence of 4508 base pairs of human genomic DNA which contain the human serine esterase gene from cytotoxic T lymphocytes (SECT) (equivalent to the 1-3E cDNA clone) and include 879 bp of 5' flanking DNA and 393 bp of 3' flanking DNA. The gene consists of five exons of 88, 148, 136, 261, and 257 nucleotides separated by four introns of 1043, 455, 205, and 643 nucleotides. The location of introns with respect to protein coding sequences in the SECT gene is identical to that of the human cathepsin G and murine granzyme B genes. Comparison of SECT gene exonic sequences to murine granzyme B-F cDNA sequences indicates similarities of 75 and 72% for granzymes B and C and 61, 59, and 61% for granzymes D, E, and F, respectively. The 5' flanking sequence of the SECT gene showed similarity only to the 5' flanking sequence of the murine granzyme B gene, indicating that these genes are homologous. Comparison of the SECT gene sequence to the human cathepsin G sequence indicated no similarity in the 5' flanking DNA although the exonic sequences show 64% sequence similarity overall and 45% sequence similarity in the respective 3' untranslated regions. These similarities suggest that the SECT and cathepsin G genes are members of the same family of serine protease genes. Evidence from high and low stringency Southern transfer analysis of human genomic DNA indicates the presence of another gene of at least 85% sequence similarity to the SECT gene.     

DNA base changes and RNA levels in N-acetoxy-2-acetylaminofluorene-induced dihydrofolate reductase mutants of Chinese hamster ovary cells

Formerly, we isolated a series of dihydrofolate reductase-deficient Chinese hamster ovary cell mutants that were induced by N-acetoxy-2-acetylaminofluorene. Deletions and complex gene rearrangements were detected in 28% of these mutants; 72% contained putative point mutations. In the present study, we have localized the putative point mutations in the 25,000 base dhfr gene by RNase heteroduplex mapping. Assignment of a position for each mutation was successful in 16 of 19 mutants studied. We cloned DNA fragments containing the mapped mutations from nine mutants into a bacteriophage lambda vector. In the case of 11 other mutants, DNA was amplified by the polymerase chain reaction procedure. Sequence analysis of cloned and amplified DNA confirmed the presence of point mutations. Most mutants (90%) carried base substitutions; the rest contained frameshift mutations. Of the point mutations, 75% were G.C to T.A transversions in either the dhfr coding sequence or at splice sites; transition G.C to A.T mutations were found in two mutants (10%). In one of these transition mutants, the base substitution occurred at the fifth base of the third intron. Of the frameshift mutations, one was a deletion of G.C pair and the other was an insertion of an A.T pair. Of the mapped mutants, 38% exhibited greatly reduced (approximately 10-fold) steady-state levels of dhfr mRNA. All eight sequenced mutants displaying this phenotype contained premature chain termination codons. Normal levels of dhfr mRNA were observed in five missense mutants and in five mutants carrying nonsense codons in the translated portion of exon VI. Taken together with the results of other mutagens at this locus, we conclude that the low dhfr mRNA phenotype is correlated with the presence of nonsense codons in exons II to V but not in the last exon of the dhfr gene.     

Trends in the 5' vs. 3' flanks of oligonucleotides in eukaryotic and prokaryotic genomes: the asymmetric roles played by cytosine and guanine.

Studies of sequence context preferences of oligonucleotides composed of (G/C)n and (A/T)m blocks (n + m = 3,4,5) unravel strong patterns. Comparisons of the 5' and 3' nearest neighbor doublets flanking these oligomers reveal the preference of (G/C)2 to be positioned immediately next to the (A/T)m block, enclosing it by (G/C) nucleotides rather than extending the (G/C)n block. That is, for a (G/C)n(A/T)m oligomer and a (G/C)2 doublet, (G/C)n(A/T)m(G/C)2 greater than (G/C)n + 2 (A/T)m. Similarly for an (A/T)m(G/C)n oligomer, (G/C)2(A/T)m(G/C)n greater than (A/T)m(G/C)n + 2. In an analogous manner, (A/T)2 flanking doublets prefer enclosing the (G/C)n blocks, although these patterns are weaker. Here we show a strong, direct relationship between the magnitude of the trends and the presence of Cs in the (G/C)n block in the (G/C)n(A/T)m oligomer, and the presence of Gs in the complementary (A/T)m(G/C)n oligomers. The trends are stronger in eukaryotic than in prokaryotic sequences. They are stronger for longer (G/C)n and shorter (A/T)m blocks. We suggest that the preference for (A/T)m to be enclosed by (G/C) rather than be flanked by them on only one side is related to DNA structure and DNA-protein interaction. Sequences of the (G/C)(A/T)(G/C) type may have more homogeneous minor groove geometry. In particular, the strong G vs. C asymmetry in the trends may be related to pyrimidine-purine junctions, possibly to CG sequences.     

Processing of mispaired and unpaired bases in heteroduplex DNA in E. coli

Bacteriophage lambda and phi X 174 DNAs, carrying sequenced mutations, have been used to construct in vitro defined species of heteroduplex DNA. Such heteroduplex DNAs were introduced by transfection, as single copies, into E. coli host cells. The progeny of individual heteroduplex molecules from each infective center was analyzed. The effect of the presence of GATC sequences (phi X 174 system) and of their methylation (lambda system) was tested. The following conclusions can be drawn: some mismatched base pairs trigger the process of mismatch repair, causing a localized strand-to-strand information transfer in heteroduplex DNA: transition mismatches G:T and A:C are efficiently repaired, whereas the six transversion mismatches are not always readily recognized and/or repaired. The recognition of transversion mismatches appears to depend on the neighbouring nucleotide sequence; single unpaired bases (frameshift mutation "mismatches") are recognized and repaired, some equally efficiently on both strands (longer and shorter), some more efficiently on the shorter (-1) strand; large non-homologies (about 800 bases) are not repaired by the Mut H, L, S, U system, but some other process repairs the non-homology with a relatively low efficiency; full methylation of GATC sequences inhibits mismatch repair on the methylated strand: this is the chemical basis of strand discrimination (old/new) in mismatch correction; unmethylated GATC sequences appear to improve mismatch repair of a G:T mismatch in phi X 174 DNA, but there may be some residual mismatch repair in GATC-free phi X 174, at least for some mismatches.(ABSTRACT TRUNCATED AT 250 WORDS)     

Conserved signals around the 5' splice sites in eukaryotic nuclear precursor mRNAs: G-runs are frequent in the introns and C in the exons near both 5' and 3' splice sites

It is known that the GT doublet is well conserved at the 5' exon/intron splice junction and is frequently embedded in the AGGT quartet. Although only the underlined G is invariable, splicing and ligation are accurately executed. In this work we search for additional conserved potential signals which may aid in 5' splice site recognition. Extensive searches which are not limited to a preconceived consensus sequence are carried out. We investigate the distributions of the 256 quartets in a 1000 nucleotide span around the 5' splice sites in approximately 1700 eukaryotic nuclear precursor mRNAs. Several potential signals are noted. Of particular interest are quartets containing runs of G, e.g., G4, G3T, G3C, G3A and AG3 in the intron immediately downstream and some C-containing quartets in the exon upstream of the 5' splice site. In an analogous calculation, (A)GGG(A) has also been found to be frequent in the intron, 60 nucleotides upstream and (A)CCC(A) in the exon downstream of the 3' splice site. These results are consistent with the recent indications that exon sequences may play a role in efficient splicing. Some models are proposed.     

Gaucher disease: A G+1----A+1 IVS2 splice donor site mutation causing exon 2 skipping in the acid beta-glucosidase mRNA.

Gaucher disease is the most frequent lysosomal storage disease and the most prevalent Jewish genetic disease. About 30 identified missense mutations are causal to the defective activity of acid beta-glucosidase in this disease. cDNAs were characterized from a moderately affected 9-year-old Ashkenazi Jewish Gaucher disease type 1 patient whose 80-year-old, enzyme-deficient, 1226G (Asn370----Ser [N370S]) homozygous grandfather was nearly asymptomatic. Sequence analyses revealed four populations of cDNAs with either the 1226G mutation, an exact exon 2 (delta EX2) deletion, a deletion of exon 2 and the first 115 bp of exon 3 (delta EX2-3), or a completely normal sequence. About 50% of the cDNAs were the delta EX2, the delta EX2-3, and the normal cDNAs, in a ratio of 6:3:1. Specific amplification and characterization of exon 2 and 5' and 3' intronic flanking sequences from the structural gene demonstrated clones with either the normal sequence or with a G+1----A+1 transition at the exon 2/intron 2 boundary. This mutation destroyed the splice donor consensus site (U1 binding site) for mRNA processing. This transition also was present at the corresponding exon/intron boundary of the highly homologous pseudogene. This new mutation, termed "IVS2 G+1----A+1," is the first splicing mutation described in Gaucher disease and accounted for about 3.4% of the Gaucher disease alleles in the Ashkenazi Jewish population. The occurrence of this "pseudogene"-type mutation in the structural gene indicates the role of acid beta-glucosidase pseudogene and structural gene rearrangements in the pathogenesis of this disease.