Analysis of primary Wilms tumors for WT1 RNA processing alterations

Wilms tumor (WT) is an embryonal renal tumor with a heterogeneous genetic etiology that serves as a valuable model for studying tumorigenesis. Biallelic inactivation of the tumor suppressor gene WT1, a zinc-finger transcriptional regulator located at 11p13, is critical for the development of some Wilms tumors. Interestingly, WT1 genomic analysis has demonstrated mutations in less than 20% of WT cases. This suggests either other genes play a more major role in Wilms tumorigenesis or WT1 is functionally altered by mechanisms other than DNA mutation. Previous observations in rat and in WT xenograft cell lines have suggested that abnormal WT1 RNA processing (exon 6 RNA editing and aberrant exon 2 splicing, respectively) is a potential mechanism of altering WT1 function in the absence of a WT1 DNA mutation. However, the role of this abnormal RNA processing has not previously been assessed in primary Wilms tumors. ^ To test the hypothesis that abnormal WT1 RNA processing is a mechanism of WT1alteration during tumor development, WT1 RNA from 85 primary tumors was analyzed using reverse transcription and polymerase chain reaction amplification (RT-PCR). Although no evidence for WT1 RNA editing was observed, variable levels (5% to 50%) of aberrant WT1 exon 2 splicing were detected for 11 tumors in the absence of a detectable WT1 DNA mutation. Also, alteration of normal WT1 alternative splicing, observed as RNA isoform loss, was detected in five tumors with no apparent WT1 genomic alteration, although no consistent pattern of RNA isoform loss was detected. This abnormal WT1 splicing, detected by either loss of exon 2 from some of the transcripts or loss of RNA isoforms, is statistically correlated with relapse (p = 0.005). These studies demonstrate that abnormal WT1 RNA processing is not a common mechanism of abrogating normal WT1 function in primary tumors. However, in those cases in which abnormal WTI splicing is present, these data indicate that it may serve as a useful prognostic marker for relapse in WT patients. ^

DOUBLE-STRAND BREAK REPAIR PATHWAYS IN DNA STRUCTURE-INDUCED GENETIC INSTABILITY

Genetic instability in mammalian cells can occur by many different mechanisms. In the absence of exogenous sources of DNA damage, the DNA structure itself has been implicated in genetic instability. When the canonical B-DNA helix is naturally altered to form a non-canonical DNA structure such as a Z-DNA or H-DNA, this can lead to genetic instability in the form of DNA double-strand breaks (DSBs) (1, 2). Our laboratory found that the stability of these non-B DNA structures was different in mammals versus Escherichia coli (E.coli) bacteria (1, 2). One explanation for the difference between these species may be a result of how DSBs are repaired within each species. Non-homologous end-joining (NHEJ) is primed to repair DSBs in mammalian cells, while bacteria that lack NHEJ (such as E.coli), utilize homologous recombination (HR) to repair DSBs. To investigate the role of the error-prone NHEJ repair pathway in DNA structure-induced genetic instability, E.coli cells were modified to express genes to allow for a functional NHEJ system under different HR backgrounds. The Mycobacterium tuberculosis NHEJ sufficient system is composed of Ku and Ligase D (LigD) (3). These inducible NHEJ components were expressed individually and together in E.coli cells, with or without functional HR (RecA/RecB), and the Z-DNA and H-DNA-induced mutations were characterized. The Z-DNA structure gave rise to higher mutation frequencies compared to the controls, regardless of the DSB repair pathway(s) available; however, the type of mutants produced after repair was greatly dictated on the available DSB repair system, indicated by the shift from 2% large-scale deletions in the total mutant population to 24% large-scale deletions when NHEJ was present (4). This suggests that NHEJ has a role in the large deletions induced by Z-DNA-forming sequences. H-DNA structure, however, did not exhibit an increase in mutagenesis in the newly engineered E.coli environment, suggesting the involvement of other factors in regulating H-DNA formation/stability in bacterial cells. Accurate repair by established DNA DSB repair pathways is essential to maintain the stability of eukaryotic and prokaryotic genomes and our results suggest that an error-prone NHEJ pathway was involved in non-B DNA structure-induced mutagenesis in both prokaryotes and eukaryotes.

Evolutionary conservation of residues in vertebrate DNA polymerase N conferring low fidelity and bypass activity.

POLN is a nuclear A-family DNA polymerase encoded in vertebrate genomes. POLN has unusual fidelity and DNA lesion bypass properties, including strong strand displacement activity, low fidelity favoring incorporation of T for template G and accurate translesion synthesis past a 5S-thymine glycol (5S-Tg). We searched for conserved features of the polymerase domain that distinguish it from prokaryotic pol I-type DNA polymerases. A Lys residue (679 in human POLN) of particular interest was identified in the conserved 'O-helix' of motif 4 in the fingers sub-domain. The corresponding residue is one of the most important for controlling fidelity of prokaryotic pol I and is a nonpolar Ala or Thr in those enzymes. Kinetic measurements show that K679A or K679T POLN mutant DNA polymerases have full activity on nondamaged templates, but poorly incorporate T opposite template G and do not bypass 5S-Tg efficiently. We also found that a conserved Tyr residue in the same motif not only affects sensitivity to dideoxynucleotides, but also greatly influences enzyme activity, fidelity and bypass. Protein sequence alignment reveals that POLN has three specific insertions in the DNA polymerase domain. The results demonstrate that residues have been strictly retained during evolution that confer unique bypass and fidelity properties on POLN.

Mutation analysis of the PVRL1 gene in caucasians with nonsyndromic cleft lip/palate.

Nonsyndromic cleft lip with or without cleft palate (nsCL/P, MIM 119530) is perhaps the most common major birth defect. Homozygous PVRL1 loss-of-function mutations result in an autosomal recessive CL/P syndrome, CLPED1, and a PVRL1 nonsense mutation is associated with sporadic nsCL/P in Northern Venezuela. To address the more general role of PVRL1 variation in risk of nsCL/P, we carried out mutation analysis of PVRL1 in North American and Australian nsCL/P cases and population-matched controls. We identified a total of 15 variants, 5 of which were seen in both populations and 1 of which, an in-frame insertion at Glu442, was more frequent in patients than in controls in both populations, though the difference was not statistically significant. Another variant, which is specific to the PVRL1 beta (HIgR) isoform, S447L, was marginally associated with nsCL/P in North American Caucasian patients, but not in Australian patients, and overall variants that affect the beta-isoform were significantly more frequent among North American patients. One Australian patient had a splice junction mutation of PVRL1. Our results suggest that PVRL1 may play a minor role in susceptibility to the occurrence of nsCL/P in some Caucasian populations, and that variation involving the beta (HIgR) isoform might have particular importance for risk of orofacial clefts. Nevertheless, these results underscore the need for studies that involve very large numbers when assessing the possible role of rare variants in risk of complex traits such as nsCL/P.

Severe aortic and arterial aneurysms associated with a TGFBR2 mutation.

BACKGROUND: A 24-year-old man presented with previously diagnosed Marfan's syndrome. Since the age of 9 years, he had undergone eight cardiovascular procedures to treat rapidly progressive aneurysms, dissection and tortuous vascular disease involving the aortic root and arch, the thoracoabdominal aorta, and brachiocephalic, vertebral, internal thoracic and superior mesenteric arteries. Throughout this extensive series of cardiovascular surgical repairs, he recovered without stroke, paraplegia or renal impairment. INVESTIGATIONS: CT scans, arteriogram, genetic mutation screening of transforming growth factor beta receptors 1 and 2. DIAGNOSIS: Diffuse and rapidly progressing vascular disease in a patient who met the diagnostic criteria for Marfan's syndrome, but was later rediagnosed with Loeys-Dietz syndrome. Genetic testing also revealed a de novo mutation in transforming growth factor beta receptor 2. MANAGEMENT: Regular cardiovascular surveillance for aneurysms and dissections, and aggressive surgical treatment of vascular disease.

Detection and sequence analysis of a nickel-induced mutation in a retroviral model system

We have developed a novel way to assess the mutagenicity of environmentally important metal carcinogens, such as nickel, by creating a positive selection system based upon the conditional expression of a retroviral transforming gene. The target gene is the v-mos gene in MuSVts110, a murine retrovirus possessing a growth temperature dependent defect in expression of the transforming gene due to viral RNA splicing. In normal rat kidney cells infected with MuSVts110 (6m2 cells), splicing of the MuSVts110 RNA to form the mRNA from which the transforming protein, p85$\sp{\rm gag-mos}$, is translated is growth-temperature dependent, occurring at 33 C and below but not at 39 C and above. This splicing "defect" is mediated by cis-acting viral sequences. Nickel chloride treatment of 6m2 cells followed by growth at 39 C, allowed the selection of "revertant" cells which constitutively express p85$\sp{\rm gag-mos}$ due to stable changes in the viral RNA splicing phenotype, suggesting that nickel, a carcinogen whose mutagenicity has not been well established, could induce mutations in mammalian genes. We also show by direct sequencing of PCR-amplified integrated MuSVts110 DNA from a 6m2 nickel-revertant cell line that the nickel-induced mutation affecting the splicing phenotype is a cis-acting 70-base duplication of a region of the viral DNA surrounding the 3$\sp\prime$ splice site. These findings provide the first example of the molecular basis for a nickel-induced DNA lesion and establish the mutagenicity of this potent carcinogen. ^

Estimation of male-to-female ratios of mutation rate in primates and rodents

Haldane (1935) developed a method for estimating the male-to-female ratio of mutation rate ($\alpha$) by using sex-linked recessive genetic disease, but in six different studies using hemophilia A data the estimates of $\alpha$ varied from 1.2 to 29.3. Direct genomic sequencing is a better approach, but it is laborious and not readily applicable to non-human organisms. To study the sex ratios of mutation rate in various mammals, I used an indirect method proposed by Miyata et al. (1987). This method takes advantage of the fact that different chromosomes segregate differently between males and females, and uses the ratios of mutation rate in sequences on different chromosomes to estimate the male-to-female ratio of mutation rate. I sequenced the last intron of ZFX and ZFY genes in 6 species of primates and 2 species of rodents; I also sequenced the partial genomic sequence of the Ube1x and Ube1y genes of mice and rats. The purposes of my study in addition to estimation of $\alpha$'s in different mammalian species, are to test the hypothesis that most mutations are replication dependent and to examine the generation-time effect on $\alpha$. The $\alpha$ value estimated from the ZFX and ZFY introns of the six primate specise is ${\sim}$6. This estimate is the same as an earlier estimate using only 4 species of primates, but the 95% confidence interval has been reduced from (2, 84) to (2, 33). The estimate of $\alpha$ in the rodents obtained from Zfx and Zfy introns is ${\sim}$1.9, and that deriving from Ube1x and Ube1y introns is ${\sim}$2. Both estimates have a 95% confidence interval from 1 to 3. These two estimates are very close to each other, but are only one-third of that of the primates, suggesting a generation-time effect on $\alpha$. An $\alpha$ of 6 in primates and 2 in rodents are close to the estimates of the male-to-female ratio of the number of germ-cell divisions per generation in humans and mice, which are 6 and 2, respectively, assuming the generation time in humans is 20 years and that in mice is 5 months. These findings suggest that errors during germ-cell DNA replication are the primary source of mutation and that $\alpha$ decreases with decreasing length of generation time. ^

Functional studies on DNA double-strand break repair proteins (MmRad51, Ku80) in mice

RecA in Escherichia coli and it's homologue, ScRad51 in Saccharomyces cerevisiae, play important roles in recombinational repair. ScRad51 homologues have been discovered in a wide range of organisms including Schizosaccharomyces pombe, lily, chicken, mouse and human. To date there is no direct evidence to describe that mouse Rad51(MmRad51) is involved in DNA double-strand break repair. In order to elucidate the role of MmRad51 in vivo, it was mutated by the embryonic stem (ES) cell/gene targeting technology in mice. The mutant embryos arrested in development shortly after implantation. There was a decrease in cell proliferation followed by programmed cell death, and trophectoderm-derived cells were sensitive to $\gamma$-radiation. Severe chromosome loss was observed in most mitotically dividing cells. The mutant embryos lived longer and developed further in a p53 mutant background; however, double-mutant embryonic fibroblasts failed to proliferate in tissue culture, reflecting the embryos limited life span. Based on these data, MmRad51 repairs DNA damage induced by $\gamma$-radiation, is needed to maintain euplody, and plays an important role in proliferating cells.^ Ku is a heterodimer of 70 and 80 kDs subunit, which binds to DNA ends and other altered DNA structures such as hairpins, nicks, and gaps. In addition, Ku is required for DNA-PK activity through a direct association. Although the biochemical properties of Ku and DNA-PKcs have been characterized in cells, their physiological functions are not clear. In order to understand the function of Ku in vivo, we generated mice homozygous for a mutation of the Ku80 gene. Ku80-deficient mice, like scid mice, showed severe immunodeficiency due to a impairment of V(D)J recombination. Mutant mice were semiviable and runted, cells derived from mutant embryos displayed hypersensitivity to $\gamma$-radiation, a decreased growth rate, a slow entry into S phase, altered colony size distributions, and a short life span. Based on these results, mutant cells and mice appeared to prematurely age. ^

Role of {\it ERCC1\/} in DNA repair and genetic recombination

The ERCC1 (Excision Repair Cross-Complementing-1) gene is the presumptive mammalian homolog of the Saccharomyces cerevisiae RAD10 gene. In mammalian NER, the Ercc1/XpF complex functions as an endonuclease that specifically recognizes 5$\sp\prime$ double-strand-3$\sp\prime$ single-strand structures. In yeast, the analogous function is performed by the Rad1/Rad10 complex. These observations and the conservation of amino acid homology between the Rad1 and XpF and the Rad10 and Ercc1 proteins has led to a general assumption of functional homology between these genes.^ In addition to NER, the Rad1/Rad10 endonuclease complex is also required in certain specialized mitotic recombination pathways in yeast. However, a similiar requirement for the endonuclease function of the Ercc1/XpF complex during genetic recombination in mammalian cells has not been directly demonstrated. The experiments performed in these studies were designed to determine if ERCC1 deficiency would produce recombination-deficient phenotypes in CHO cells similar to those observed in RAD10 deletion mutants, including: (1) decreased single-reciprocal exchange recombination, and (2) inability to process 5$\sp\prime$ sequence heterology in recombination intermediates.^ Specifically, these studies describe: (1) The isolation and characterization of the ERCC1 locus of Chinese hamster ovary cells; (2) The production of an ERCC1 null mutant cell line by targeted knock-out of the endogenous ERCC1 gene in a Chinese hamster ovary cell line, CHO-ATS49tg, which contains an endogenous locus, APRT, suitable as a chromosomal target for homologous recombination; (3) The characterization of mutant ERCC1 alleles from a panel of Chinese hamster ovary cell ERCC1 mutants derived by conventional mutagenesis; (4) An investigation of the effects of ERCC1 mutation on mitotic recombination through targeting of the APRT locus in an ERCC1 null background.^ The results of these studies strongly suggest that the role of ERCC1 in homologous recombination in mammalian cells is analogous to that of the yeast RAD10 gene. ^

MATHEMATICAL STUDIES ON THE EVOLUTIONARY CHANGE OF DNA SEQUENCES

With the aim of understanding the mechanism of molecular evolution, mathematical problems on the evolutionary change of DNA sequences are studied. The problems studied and the results obtained are as follows: (1) Estimation of evolutionary distance between nucleotide sequences. Studying the pattern of nucleotide substitution for the case of unequal substitution rates, a new mathematical formula for estimating the average number of nucleotide substitutions per site between two homologous DNA sequences is developed. It is shown that this formula has a wider applicability than currently available formulae. A statistical method for estimating the number of nucleotide changes due to deletion and insertion is also developed. (2) Biases of the estimates of nucleotide substitutions obtained by the restriction enzyme method. The deviation of the estimate of nucleotide substitutions obtained by the restriction enzyme method from the true value is investigated theoretically. It is shown that the amount of the deviation depends on the nucleotides in the recognition sequence of the restriction enzyme used, unequal rates of substitution among different nucleotides, and nucleotide frequences, but the primary factor is the unequal rates of nucleotide substitution. When many different kinds of enzymes are used, however, the amount of average deviation is generally small. (3) Distribution of restriction fragment lengths. To see the effect of undetectable restriction fragments and fragment differences on the estimate of nucleotide differences, the theoretical distribution of fragment lengths is studied. This distribution depends on the type of restriction enzymes used as well as on the relative frequencies of four nucleotides. It is shown that undetectability of small fragments or fragment differences gives a serious underestimate of nucleotide substitutions when the length-difference method of estimation is used, but the extent of underestimation is small when the site-difference method is used. (4) Evolutionary relationships of DNA sequences in finite populations. A mathematical theory on the expected evolutionary relationships among DNA sequences (nucleons) randomly chosen from the same or different populations is developed under the assumption that the evolutionary change of nucleons is determined solely by mutation and random genetic drift. . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author). UMI ^

STUDIES ON THE PATTERNS OF NUCLEOTIDE AND AMINO ACID SUBSTITUTION (ELECTROPHORESIS, MUTATION RATE, SELECTION, BASE, COMPOSITION)

Theoretical and empirical studies were conducted on the pattern of nucleotide and amino acid substitution in evolution, taking into account the effects of mutation at the nucleotide level and purifying selection at the amino acid level. A theoretical model for predicting the evolutionary change in electrophoretic mobility of a protein was also developed by using information on the pattern of amino acid substitution. The specific problems studied and the main results obtained are as follows: (1) Estimation of the pattern of nucleotide substitution in DNA nuclear genomes. The pattern of point mutations and nucleotide substitutions among the four different nucleotides are inferred from the evolutionary changes of pseudogenes and functional genes, respectively. Both patterns are non-random, the rate of change varying considerably with nucleotide pair, and that in both cases transitions occur somewhat more frequently than transversions. In protein evolution, substitution occurs more often between amino acids with similar physico-chemical properties than between dissimilar amino acids. (2) Estimation of the pattern of nucleotide substitution in RNA genomes. The majority of mutations in retroviruses accumulate at the reverse transcription stage. Selection at the amino acid level is very weak, and almost non-existent between synonymous codons. The pattern of mutation is very different from that in DNA genomes. Nevertheless, the pattern of purifying selection at the amino acid level is similar to that in DNA genomes, although selection intensity is much weaker. (3) Evaluation of the determinants of molecular evolutionary rates in protein-coding genes. Based on rates of nucleotide substitution for mammalian genes, the rate of amino acid substitution of a protein is determined by its amino acid composition. The content of glycine is shown to correlate strongly and negatively with the rate of substitution. Empirical formulae, called indices of mutability, are developed in order to predict the rate of molecular evolution of a protein from data on its amino acid sequence. (4) Studies on the evolutionary patterns of electrophoretic mobility of proteins. A theoretical model was constructed that predicts the electric charge of a protein at any given pH and its isoelectric point from data on its primary and quaternary structures. Using this model, the evolutionary change in electrophoretic mobilities of different proteins and the expected amount of electrophoretically hidden genetic variation were studied. In the absence of selection for the pI value, proteins will on the average evolve toward a mildly basic pI. (Abstract shortened with permission of author.) ^

DNA nucleotide excision repair gene single nucleotide polymorphisms and hereditary nonpolyposis colorectal cancer

Hereditary nonpolyposis colorectal cancer (HNPCC) is an autosomal dominant disease caused by germline mutations in DNA mismatch repair(MMR) genes. The nucleotide excision repair(NER) pathway plays a very important role in cancer development. We systematically studied interactions between NER and MMR genes to identify NER gene single nucleotide polymorphism (SNP) risk factors that modify the effect of MMR mutations on risk for cancer in HNPCC. We analyzed data from polymorphisms in 10 NER genes that had been genotyped in HNPCC patients that carry MSH2 and MLH1 gene mutations. The influence of the NER gene SNPs on time to onset of colorectal cancer (CRC) was assessed using survival analysis and a semiparametric proportional hazard model. We found the median age of onset for CRC among MMR mutation carriers with the ERCC1 mutation was 3.9 years earlier than patients with wildtype ERCC1(median 47.7 vs 51.6, log-rank test p=0.035). The influence of Rad23B A249V SNP on age of onset of HNPCC is age dependent (likelihood ratio test p=0.0056). Interestingly, using the likelihood ratio test, we also found evidence of genetic interactions between the MMR gene mutations and SNPs in ERCC1 gene(C8092A) and XPG/ERCC5 gene(D1104H) with p-values of 0.004 and 0.042, respectively. An assessment using tree structured survival analysis (TSSA) showed distinct gene interactions in MLH1 mutation carriers and MSH2 mutation carriers. ERCC1 SNP genotypes greatly modified the age onset of HNPCC in MSH2 mutation carriers, while no effect was detected in MLH1 mutation carriers. Given the NER genes in this study play different roles in NER pathway, they may have distinct influences on the development of HNPCC. The findings of this study are very important for elucidation of the molecular mechanism of colon cancer development and for understanding why some mutation carriers of the MSH2 and MLH1 gene develop CRC early and others never develop CRC. Overall, the findings also have important implications for the development of early detection strategies and prevention as well as understanding the mechanism of colorectal carcinogenesis in HNPCC. ^

Role of telomere dysfunction, DNA damage response and p53 mutations in tumorigenesis and aging

The ends of eukaryotic chromosomes are protected by specialized ribonucleoprotein structures termed telomeres. Telomeres protect chromosomes from end-to-end fusions, inappropriate repair and degradation. Disruption of this complex activates an ATM/ATR DNA damage response (DDR) pathway. One component of the complex is the Protection Of Telomeres 1 (POT1) protein, an evolutionarily conserved protein which binds single-stranded 3' overhang and is required for both chromosomal end protection and telomere length regulation. The mouse contains two POT1 orthologs, Pot1a and Pot1b. Here we show that both proteins colocalize with telomeres through interaction with the adapter protein TPP1. In addition, compared to Pot1a, the OB-folds of Pot1b possess less sequence specificity for telomeres. Disruption of POT1 proteins result in telomere dysfunction and activation of an ATR-dependent DDR at telomeres, suggesting that this response is normally suppressed by POT1 binding to the single-stranded G-overhang. ^ Telomeres are maintained by telomerase, and its absence in somatic cells results in telomere progressive loss that triggers the activation of p53. Telomere dysfunction initiates genomic instability and induces both p53-dependent replicative senescence and apoptosis to suppress tumorigenesis. In the absence of functional p53, this genomic instability promotes cancer. It was previously not known which aspect of the p53 dependent DNA damage response is important to suppress tumorigenesis initiated by dysfunctional telomeres. The p53R172P knock-in mouse, which is unable to induce apoptosis but retains intact cell cycle arrest/cellular senescence pathways, allowed us to examine whether p53-dependent apoptosis is a major tumor suppression pathway initiated in the setting of telomere dysfunction. Spontaneous tumorigenesis remains potently suppressed in late generation telomerase null mice possessing the p53P/P mutation. These results suggest that suppression of spontaneous tumorigenesis initiated by dysfunctional telomeres requires activation of a p53-dependent senescence pathway. In addition, we used another knock-in mouse model with a p53R172H (p53H) point mutation to test the hypothesis that telomere dysfunction promotes chromosomal instability and accelerates the onset of tumorigenesis in vivo in the setting of this most common gain-of-function mutation in the human Li Fraumeni cancer syndrome. We unexpectedly observed that telomerase null mice possessing dysfunctional telomeres in the setting of the p53H/+ mutation develop significantly fewer tumors, die prematurely and exhibit higher level of cellular senescence, apoptosis and elevated genomic instability compared to telomerase intact p53H/+ and telomerase null p53+/+ mice. These contrasting results thus link cancer and aging to the functional status of telomeres and the integrity of the p53 pathway. ^

NRAS and BRAF mutations in primary cutaneous melanoma: A comparison of mutation rates between radial and vertical growth phases in individual tumors

Primary cutaneous melanoma is a cancer arising from melanocytes in the skin. In recent decades the incidence of this malignancy has increased significantly. Mortality rates are high for patients with tumors measuring over a few millimeters in thickness. Response rates to conventional radiation and chemotherapy are very low in patients with metastatic melanoma. New therapies targeting melanoma’s aberrant cell signaling pathways such as the MAP Kinase pathway are being developed. Mutations of NRAS and BRAF genes are quite common in cutaneous melanoma and lead to constitutive activation of the MAP Kinase pathway. This study tests the hypothesis that NRAS and BRAF mutations increase as a tumor progresses from the noninvasive radial growth phase (RGP) to the invasive vertical growth phase (VGP). Laser capture microdissection was used to obtain separate, pure tumor DNA samples from the RGP and VGP of thirty primary cutaneous melanomas. PCR was used to amplify NRAS exon 2 and BRAF exon 15 tumor DNA. The amplified DNA was sequenced and analyzed for mutations. An overall mutation rate of 74% was obtained for the twenty-three melanomas in which there were complete sequence results. With the exception of one melanoma NRAS and BRAF mutations were mutually exclusive. All seven NRAS exon 2 mutations involved codon 61. Three of these melanomas had mutations in both the RGP and VGP. The remaining four tumors were wild type for NRAS exon 2 in the RGP but mutated in the VGP. Of the fifteen BRAF exon 15 mutated melanomas all but one involved codon 600. Twelve of the fifteen BRAF exon 15 mutations were the T1799A type. Nine of the fifteen BRAF mutated tumors had the same mutation in both the RGP and VGP. Five of fifteen melanomas had wild type RGP DNA and BRAF exon 15 mutated VGP DNA. A single melanoma had BRAF exon 15 mutated DNA in the RGP and wild type DNA in the VGP. Overall, these results suggest a trend toward the acquisition of NRAS and BRAF mutations as cutaneous melanomas change from a noninvasive to an invasive, potentially deadly cancer.^