SeqDoC: rapid SNP and mutation detection by direct comparison of DNA sequence chromatograms

Mastering DNA chromatogram analysis in Sanger sequencing for reliable clinical analysis

BACKGROUND

Sanger dideoxy sequencing is vital in clinical analysis due to its accuracy, ability to analyze genetic markers like SNPs and STRs, capability to generate reliable DNA profiles, and its role in resolving complex clinical cases. The precision and robustness of Sanger sequencing contribute significantly to the scientific basis of clinical investigations. Though the reading of chromatograms seems to be a routine step, many errors conducted in PCR may lead to consequent limitations in the readings of AGCT peaks. These errors are possibly associated with improper DNA amplification and its subsequent interpretation of DNA sequencing files, such as noisy peaks, artifacts, and confusion between double-peak technical errors, heterozygosity, and double infection potentials. Thus, it is not feasible to read nucleic acid sequences without giving serious attention to these technical problems. To ensure the accuracy of DNA sequencing outcomes, it is also imperative to detect and rectify technical challenges that may lead to misinterpretation of the DNA sequence, resulting in errors and incongruities in subsequent analyses.

SHORT CONCLUSION

This overview sheds light on prominent technical concerns that can emerge prior to and during the interpretation of DNA chromatograms in Sanger sequencing, along with offering strategies to address them effectively. The significance of identifying and tackling these technical limitations during the chromatogram analysis is underscored in this review. Recognizing these concerns can aid in enhancing the quality of downstream analyses for Sanger sequencing results, which holds notable improvement in accuracy, reliability, and ability to provide crucial genetic information in clinical analysis.

Emergence of Ancylostoma caninum parasites with the benzimidazole resistance F167Y polymorphism in the US dog population

BACKGROUND

Anthelmintic resistance to benzimidazole has been detected in the canine hookworm, Ancylostoma caninum. Benzimidazole resistance is believed to have developed originally in greyhounds, but has also been detected in non-greyhound pet dogs. The aim of this study was to validate a probe-based allele-specific real-time PCR tests for the F167Y polymorphism on the β-tubulin isotype-1 gene and to determine the geographic distribution.

METHODS

Allele-specific real-time PCR tests were established and validated to detect the codon 167 polymorphism in the Ancylostoma caninum β-tubulin isotype-1gene. Additionally, real-time PCR tests were validated for Ancylostoma spp. and Uncinaria stenocephala. Two nucleic acid extraction protocols were validated including mechanical disruption of parasite structures in stool. The frequency of the F167Y single nucleotide polymorphism (SNP) was determined in hookworm confirmed stool samples. Samples with the resistant 167Y genotype were confirmed by β-tubulin gene sequencing and allele frequencies were determined.

RESULTS

The Ancylostoma spp. and A. caninum F167Y allele-specific real-time PCR tests were highly sensitive and specific when tested against synthetic DNA, spiked samples, and characterized parasites. Using an optimized total nucleic acid extraction protocol, 54 of 511 (10.6%) were found to contain the benzimidazole resistance allele. All 55 samples containing hookworms with the resistance mutation were confirmed by β-tubulin gene sequencing. The majority of resistant hookworms (44 resistant, 183 tested; 24.4%) originated from Florida, five from California (103 tested, 4.9%), three from Idaho (40 tested, 7.5%), two from Nevada (22 tested, 9.1%), and one sample from Hawaii (13 tested, 7.7%). Resistant genotypes were found in 14 different dog breeds including eight in Greyhounds. Allele-frequency determination revealed resistance allele frequencies between 1 and 100% with 58% above 50%.

CONCLUSIONS

This data strongly supports recent findings of benzimidazole resistant canine hookworms present throughout the general US pet dog population.

A novel DNA methylation motif identified in Bacillus pumilus BA06 and possible roles in the regulation of gene expression

Single-molecule real-time (SMRT) sequencing can be used to identify a wide variety of chemical modifications of the genome, such as methylation. Here, we applied this approach to identify N6-methyl-adenine (m6A) and N4-methyl-cytosine (m4C) modification in the genome of Bacillus pumilus BA06. A typical methylation recognition motif of the type I restriction-modification system (R-M), 5'-TC^m6AN8TTGG-3'/3'-AGTN8^m6AACC-5', was identified. We confirmed that this motif was a new type I methylation site using REBASE analysis and that it was recognized by a type I R-M system, Bpu6ORFCP, according to methylation sensitivity assays in vivo and vitro. Furthermore, we found that deletion of the R-M system Bpu6ORFCP induced transcriptional changes in many genes and led to increased gene expression in pathways related to ABC transporters, sulfur metabolism, ribosomes, cysteine and methionine metabolism and starch and sucrose metabolism, suggesting that the R-M system in B. pumilus BA06 has other significant biological functions beyond protecting the B. pumilus BA06 genome from foreign DNA.

Single-Nucleotide-Resolution Computing and Memory in Living Cells

The ability to process and store information in living cells is essential for developing next-generation therapeutics and studying biology in situ. However, existing strategies have limited recording capacity and are challenging to scale. To overcome these limitations, we developed DOMINO, a robust and scalable platform for encoding logic and memory in bacterial and eukaryotic cells. Using an efficient single-nucleotide-resolution Read-Write head for DNA manipulation, DOMINO converts the living cells' DNA into an addressable, readable, and writable medium for computation and storage. DOMINO operators enable analog and digital molecular recording for long-term monitoring of signaling dynamics and cellular events. Furthermore, multiple operators can be layered and interconnected to encode order-independent, sequential, and temporal logic, allowing recording and control over the combination, order, and timing of molecular events in cells. We envision that DOMINO will lay the foundation for building robust and sophisticated computation-and-memory gene circuits for numerous biotechnological and biomedical applications.

Recurrent triploidy due to a failure to complete maternal meiosis II: whole-exome sequencing reveals candidate variants

Triploidy is a relatively common cause of miscarriage; however, recurrent triploidy has rarely been reported. A healthy 34-year-old woman was ascertained because of 18 consecutive miscarriages with triploidy found in all 5 karyotyped losses. Molecular results in a sixth loss were also consistent with triploidy. Genotyping of markers near the centromere on multiple chromosomes suggested that all six triploid conceptuses occurred as a result of failure to complete meiosis II (MII). The proband's mother had also experienced recurrent miscarriage, with a total of 18 miscarriages. Based on the hypothesis that an inherited autosomal-dominant maternal predisposition would explain the phenotype, whole-exome sequencing of the proband and her parents was undertaken to identify potential candidate variants. After filtering for quality and rarity, potentially damaging variants shared between the proband and her mother were identified in 47 genes. Variants in genes coding for proteins implicated in oocyte maturation, oocyte activation or polar body extrusion were then prioritized. Eight of the most promising candidate variants were confirmed by Sanger sequencing. These included a novel change in the PLCD4 gene, and a rare variant in the OSBPL5 gene, which have been implicated in oocyte activation upon fertilization and completion of MII. Several variants in genes coding proteins playing a role in oocyte maturation and early embryonic development were also identified. The genes identified may be candidates for the study in other women experiencing recurrent triploidy or recurrent IVF failure.

Genetic variants and increased expression of Parascaris equorum P-glycoprotein-11 in populations with decreased ivermectin susceptibility

Macrocyclic lactones (MLs) represent the major drug class for control of parasitic infections in humans and animals. However, recently reports of treatment failures became more frequent. In addition to human and ruminant parasitic nematodes this also is the case for the horse-nematode Parascaris equorum. Nevertheless, to date the molecular basis of ML resistance is still not understood. Unspecific resistance mechanisms involving transporters such as P-glycoproteins (Pgps) are expected to contribute to ML resistance in nematodes. Here, complete sequences of two P. equorum Pgps were cloned and identified as orthologs of Caenorhabditis elegans Ppg-11 and an unnamed Caenorhabditis briggsae Pgp designated as Pgp-16 using phylogenetic analysis. Quantitative real-time PCR was used to compare expression between tissues. Significantly higher PeqPgp-11 expression was found in the gut for both genders, whereas for PeqPgp-16 the body wall was identified as predominant expression site. Furthermore, Pgps were analyzed regarding their participation in resistance development. Using SeqDoC analyses, Pgp-sequences of P. equorum populations with different ML susceptibility were compared. This approach revealed three single nucleotide polymorphisms (SNPs) causing missense mutations in the PeqPgp-11 sequence which correlated with decreased ML susceptibility. However, no resistance associated differences in mRNA expression levels were detected between embryonated eggs of these populations. In contrast, comparison of two pre-adult groups with different ivermectin (IVM) susceptibility revealed the presence of the three SNPs and in addition statistically significant PeqPgp-11 overexpression in the group of worms with reduced susceptibility. These results indicate that Pgp-11 might be involved in IVM resistance in P. equorum as it shows increased expression in an IVM exposed life-cycle stage of an IVM resistant population as well as occurrence of putatively resistance associated SNPs in populations with reduced IVM susceptibility. These SNPs are promising diagnostic candidates for detection of ML resistance with potential also for other parasitic nematode species.

Potential contribution of P-glycoproteins to macrocyclic lactone resistance in the cattle parasitic nematode Cooperia oncophora

Resistance against macrocyclic lactones such as ivermectin is widespread among parasitic gastrointestinal nematodes of small ruminants and is rapidly increasing in cattle parasites. ABC transporters of the subfamily B, the so-called P-glycoproteins (Pgps) have been frequently implicated in ivermectin resistance and are a major cause of multi-drug resistance in protozoa and helminths. The Pgp inhibitor verapamil (VPL) dramatically enhanced susceptibility of the cattle parasitic nematode Cooperia oncophora to ivermectin in vitro as measured in a larval developmental assay and a larval migration inhibition assay using third stage larvae. Moreover, VPL completely restored susceptibility to ivermectin in a resistant isolate resulting in virtually identical dose-response curves of susceptible and resistant isolates in the presence of VPL. Further characterisation of the molecular mechanisms resulting in Pgp-mediated ivermectin resistance is still hampered by the lack of molecular and biochemical information for Pgps of parasitic nematodes. Using PCR with degenerate primers, fragments of four different C. oncophora Pgps could be amplified and the Conpgp-2, previously implicated in macrocyclic lactone resistance in Haemonchus contortus, and Conpgp-3 full-length cDNAs were obtained by RACE PCR. The pgp sequences presented here contribute important data required to systematically screen resistant C. oncophora isolates for up- or down-regulation of Pgps and for the detection of single nucleotide polymorphisms in Pgps to detect selection of specific Pgp alleles by anthelmintics as early as possible.

Rare allelic forms of PRDM9 associated with childhood leukemogenesis

One of the most rapidly evolving genes in humans, PRDM9, is a key determinant of the distribution of meiotic recombination events. Mutations in this meiotic-specific gene have previously been associated with male infertility in humans and recent studies suggest that PRDM9 may be involved in pathological genomic rearrangements. In studying genomes from families with children affected by B-cell precursor acute lymphoblastic leukemia (B-ALL), we characterized meiotic recombination patterns within a family with two siblings having hyperdiploid childhood B-ALL and observed unusual localization of maternal recombination events. The mother of the family carries a rare PRDM9 allele, potentially explaining the unusual patterns found. From exomes sequenced in 44 additional parents of children affected with B-ALL, we discovered a substantial and significant excess of rare allelic forms of PRDM9. The rare PRDM9 alleles are transmitted to the affected children in half the cases; nonetheless there remains a significant excess of rare alleles among patients relative to controls. We successfully replicated this latter observation in an independent cohort of 50 children with B-ALL, where we found an excess of rare PRDM9 alleles in aneuploid and infant B-ALL patients. PRDM9 variability in humans is thought to influence genomic instability, and these data support a potential role for PRDM9 variation in risk of acquiring aneuploidies or genomic rearrangements associated with childhood leukemogenesis.

A substitution in the ligand binding domain of the porcine glucocorticoid receptor affects activity of the adrenal gland

Glucocorticoids produced in the adrenal cortex under the control of the hypothalamic-pituitary axis play a vital role in the maintenance of basal and stress-related homeostasis and influence health and well-being. To identify loci affecting regulation of the hypothalamic-pituitary-adrenal (HPA) axis in the pig we performed a genome-wide association study for two parameters of acute and long-term adrenal activity: plasma cortisol level and adrenal weight. We detected a major quantitative trait locus at the position of the glucocorticoid receptor gene (NR3C1) – a key regulator of HPA axis activity. To determine the causal variant(s), we resequenced the coding region of NR3C1 and found three missense single nucleotide polymorphisms (SNPs). SNP c.1829C>T, leading to a p.Ala610Val substitution in the ligand binding domain, showed large (about 0.6× and 1.2× phenotypic standard deviations for cortisol level and adrenal weight, respectively), and highly significant (2.1E-39≤log10(1/p)≤1.7E+0) negative effects on both traits. We were able to replicate the association in three commercial pig populations with different breed origins. We analyzed effects of the p.Ala610Val substitution on glucocorticoid-induced transcriptional activity of porcine glucocorticoid receptor (GR) in vitro and determined that the substitution introduced by SNP c.1829C>T increased sensitivity of GR by about two-fold. Finally, we found that non-coding polymorphisms in linkage disequilibrium with SNP c.1829C>T have only a minor effect on the expression of NR3C1 in tissues related to the HPA axis. Our findings provide compelling evidence that SNP c.1829C>T in porcine NR3C1 is a gain-of-function mutation with a major effect on the activity of the adrenal gland. Pigs carrying this SNP could provide a new animal model to study neurobiological and physiological consequences of genetically based GR hypersensitivity and adrenal hypofunction.

Methylation by a unique α-class N4-cytosine methyltransferase is required for DNA transformation of Caldicellulosiruptor bescii DSM6725

Thermophilic microorganisms capable of using complex substrates offer special advantages for the conversion of lignocellulosic biomass to biofuels and bioproducts. Members of the gram-positive bacterial genus Caldicellulosiruptor are anaerobic thermophiles with optimum growth temperatures between 65°C and 78°C and are the most thermophilic cellulolytic organisms known. In fact, they efficiently use biomass non-pretreated as their sole carbon source and in successive rounds of application digest 70% of total switchgrass substrate. The ability to genetically manipulate these organisms is a prerequisite to engineering them for use in conversion of these complex substrates to products of interest as well as identifying gene products critical for their ability to utilize non-pretreated biomass. Here, we report the first example of DNA transformation of a member of this genus, C. bescii. We show that restriction of DNA is a major barrier to transformation (in this case apparently absolute) and that methylation with an endogenous unique α-class N4-Cytosine methyltransferase is required for transformation of DNA isolated from E. coli. The use of modified DNA leads to the development of an efficient and reproducible method for DNA transformation and the combined frequencies of transformation and recombination allow marker replacement between non-replicating plasmids and chromosomal genes providing the basis for rapid and efficient methods of genetic manipulation.

SNPmplexViewer--toward a cost-effective traceability system

Background

Beef traceability has become mandatory in many regions of the world and is typically achieved through the use of unique numerical codes on ear tags and animal passports. DNA-based traceability uses the animal's own DNA code to identify it and the products derived from it. Using SNaPshot, a primer-extension-based method, a multiplex of 25 SNPs in a single reaction has been practiced for reducing the expense of genotyping a panel of SNPs useful for identity control.

Findings

To further decrease SNaPshot's cost, we introduced the Perl script SNPmplexViewer, which facilitates the analysis of trace files for reactions performed without the use of fluorescent size standards. SNPmplexViewer automatically aligns reference and target trace electropherograms, run with and without fluorescent size standards, respectively. SNPmplexViewer produces a modified target trace file containing a normalised trace in which the reference size standards are embedded. SNPmplexViewer also outputs aligned images of the two electropherograms together with a difference profile.

Conclusions

Modified trace files generated by SNPmplexViewer enable genotyping of SnaPshot reactions performed without fluorescent size standards, using common fragment-sizing software packages. SNPmplexViewer's normalised output may also improve the genotyping software's performance. Thus, SNPmplexViewer is a general free tool enabling the reduction of SNaPshot's cost as well as the fast viewing and comparing of trace electropherograms for fragment analysis. SNPmplexViewer is available at http://cowry.agri.huji.ac.il/cgi-bin/SNPmplexViewer.cgi.

Whole-exome sequencing of pediatric acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL), the most common malignant disorder in childhood, is typically associated with numerical chromosomal aberrations, fusion genes or small focal deletions, thought to represent important pathogenetic events in the development of the leukemia. Mutations, such as single nucleotide changes, have also been reported in childhood ALL, but these have only been studied by sequencing a small number of candidate genes. Herein, we report the first unbiased sequencing of the whole exome of two cases of pediatric ALL carrying the ETV6/RUNX1 (TEL/AML1) fusion gene (the most common genetic subtype) and corresponding normal samples. A total of 14 somatic mutations were identified, including four and seven protein-altering nucleotide substitutions in each ALL. Twelve mutations (86%) occurred in genes previously described to be mutated in other types of cancer, but none was found to be recurrent in an extended series of 29 ETV6/RUNX1-positive ALLs. The number of single nucleotide mutations was similar to the number of copy number alterations as detected by single nucleotide polymorphism arrays. Although the true pathogenetic significance of the mutations must await future functional evaluations, this study provides a first estimate of the mutational burden at the genetic level of t(12;21)-positive childhood ALL.

Absence of SYCP3 mutations in women with recurrent miscarriage with at least one trisomic miscarriage

Mutations within the coding regions of the synaptonemal complex gene SYCP3 have previously been reported in women with recurrent miscarriage. The present study found no mutations in any of the coding exons or the intron/exon boundaries among 50 recurrent miscarriage patients with at least one documented trisomic miscarriage, suggesting that mutations in SYCP3 do not contribute significantly to risk for recurrent miscarriage through maternal meiotic nondisjunction.

Mutant generation in vertebrate model organisms by TILLING

TILLING (Targeting Induced Local Lesions IN Genomes) is a popular reverse genetic approach that has been successfully applied in several genetic model organisms such as zebrafish, rat, Drosophila, Arabidopsis, or medaka. In contrast to classical targeted knockout technologies that work in mice by directly targeting a gene of interest, TILLING follows an indirect strategy. The first step of the TILLING pipeline is the generation of a TILLING library that consists of large numbers of mutagenized individuals. In a second step, these individuals are screened for mutations in any gene of interest. Screening is performed by PCR amplification of specific exons from each individual of a library followed by mutation detection. This could be done, for example, by direct re-sequencing of PCR fragments or alternatively, by CEL1 endonuclease-mediated mutation discovery. Individuals carrying potentially deleterious point mutations are isolated from the library and mutant lines are established. TILLING allows the identification of a whole range of point mutations, covering nonsense, splice site, and missense mutations in only one screening round, because the generation of mutations by mutagenesis as well as the screening tools is not biased. Potential knockout mutations are initially the mutations of choice, but TILLING screens can also be used to isolate allelic series of point mutations ranging from complete null phenotypes to hypomorphic or even dominant-negative or conditional alleles. These allelic series can be helpful for a comprehensive functional analysis of a gene of interest. TILLING is applicable to any kind of genetically tractable model organism, as long as this model organism is amenable to chemical mutagenesis, and genomic sequence information for a gene of interest is available. This chapter describes the design and pipeline of a TILLING facility as we are currently operating it for zebrafish in Dresden. Protocols for mutation detection by direct re-sequencing are described in detail. However, alternatives to this pipeline do exist and will be mentioned briefly.

Comparative view of in silico DNA sequencing analysis tools

DNA sequencing is an important tool for discovery of genetic variants. The task of detecting single-nucleotide variants is complicated by noise and sequencing artifacts in sequencing data. Several in silico tools have been developed to assist this process. These tools interpret the raw chromatogram data and perform a specialized base-calling and quality-control assessment procedure to identify variants. The approach used to identify variants differs between the tools, with some specific to SNPs and other for Indels. The choice of a tool is guided by the design of the sequencing project and the nature of the variant to be discovered. In this chapter, these tools are compared to facilitate the choice of a tool used for variant discovery.

Anthelmintic resistance: markers for resistance, or susceptibility?

The Consortium for Anthelmintic Resistance and Susceptibility (CARS) brings together researchers worldwide, with a focus of advancing knowledge of resistance and providing information on detection methods and treatment strategies. Advances in this field suggest mechanisms and features of resistance that are shared among different classes of anthelmintic. Benzimidazole resistance is characterized by specific amino acid substitutions in beta-tubulin. If present, these substitutions increase in frequency upon drug treatment and lead to treatment failure. In the laboratory, sequence substitutions in ion-channels can contribute to macrocyclic lactone resistance, but there is little evidence that they are significant in the field. Changes in gene expression are associated with resistance to several different classes of anthelmintic. Increased P-glycoprotein expression may prevent drug access to its site of action. Decreased expression of ion-channel subunits and the loss of specific receptors may remove the drug target. Tools for the identification and genetic analysis of parasitic nematodes and a new online database will help to coordinate research efforts in this area. Resistance may result from a loss of sensitivity as well as the appearance of resistance. A focus on the presence of anthelmintic susceptibility may be as important as the detection of resistance.

Linkage analysis of the C4A/C4B copy number variation and polymorphisms of the adjacent steroid 21-hydroxylase gene in a healthy population

Genes encoding the steroid 21-hydroxylase (CYP21A2) and the complement component C4 proteins (C4A and C4B) are located in the MHC region in a strongly linked structure named RCCX module. Previous studies found that carriers of C4B gene deficiency (C4B*Q0) have higher risk for cardiovascular diseases. A potential explanation is that lacking the C4B gene may result in altered function of the neighboring CYP21A2 gene. Therefore we sequenced the CYP21A2 gene in 96 healthy individuals to identify polymorphisms and to characterize their linkage pattern. Fifty-three variations were detected including a new one which alters the TATA-box of the gene. Only three known mutations (V281L, Q318X and R479L) associated with congenital adrenal hyperplasia, were found in 7, 2 and 1 subjects, respectively. Linkage analysis revealed that some variations exhibit strong correlation with the C4 copy number polymorphism and constituents of the MHC III region. Rare alleles of three polymorphisms were identified as components of the 8.1 ancestral haplotype. Haplotyping and family study confirmed that the variant alleles of two intronic SNPs were constituents of haplotype blocks lacking the C4B gene. These results suggest that variations of CYP21A2 gene can be involved in disease associations of the 8.1 haplotype and the C4B*Q0 genotype.

Single nucleotide polymorphism (SNP) markers for benzimidazole resistance in veterinary nematodes

Resistance to the benzimidazole class of anthelmintics in nematodes of veterinary importance has a long history. Research into the mechanisms responsible for this resistance is subsequently at a more advanced stage than for other classes of anthelmintics. The principal mechanism of resistance to benzimidazoles is likely to involve changes in the primary structure of β-tubulins, the building blocks of microtubules. Specifically, point mutations in the β-tubulin isotype 1 gene leading to amino acid substitutions in codons 167, 198, and 200 of the protein have been associated with resistance in nematodes. These single nucleotide polymorphisms offer a means of detecting the presence of resistance within populations. In this mini-review, we focus on the prevalence and importance of these polymorphisms in three groups of nematodes: trichostrongylids, cyathostomins, and hookworms. A brief overview of existing strategies for genotyping single nucleotide polymorphisms is also presented. The CARS initiative hopes to exploit these known polymorphisms to further our understanding of the phenomenon of BZ resistance.

Rapid identification of single nucleotide substitutions using SeqDoC

Identification and characterization of nucleotide substitutions in DNA sequences for single nucleotide polymorphism or point mutation detection can be a time consuming and sometimes inaccurate process, particularly in relatively low-throughput situations where fully automated solutions may not be appropriate. SeqDoC provides a simple web-based application to simplify this identification process, by using direct subtractive comparison of the raw sequence traces to highlight differences characteristic of nucleotide substitutions. Sequencing artefacts, such as variable peak separation and signal strength, are compensated for with moving window normalisation functions, whereas the signal to noise ratio of the comparison trace is greatly enhanced by applying an algorithm to emphasise features associated with nucleotide substitutions. Analysis of the output is simple and intuitive, permitting rapid identification of points of difference between the reference and test sequence traces.

Spontaneous and UV radiation-induced multiple metastatic melanomas in Cdk4R24C/R24C/TPras mice

Human melanoma susceptibility is often characterized by germ-line inactivating CDKN2A (INK4A/ARF) mutations, or mutations that activate CDK4 by preventing its binding to and inhibition by INK4A. We have previously shown that a single neonatal UV radiation (UVR) dose delivered to mice that carry melanocyte-specific activation of Hras (TPras) increases melanoma penetrance from 0% to 57%. Here, we report that activated Cdk4 cooperates with activated Hras to enhance susceptibility to melanoma in mice. Whereas UVR treatment failed to induce melanomas in Cdk4(R24C/R24C) mice, it greatly increased the penetrance and decreased the age of onset of melanoma development in Cdk4(R24C/R24C)/TPras animals compared with TPras alone. This increased penetrance was dependent on the threshold of Cdk4 activation as Cdk4(R24C/+)/TPras animals did not show an increase in UVR-induced melanoma penetrance compared with TPras alone. In addition, Cdk4(R24C/R24C)/TPras mice invariably developed multiple lesions, which occurred rarely in TPras mice. These results indicate that germ-line defects abrogating the pRb pathway may enhance UVR-induced melanoma. TPras and Cdk4(R24C/R24C)/TPras tumors were comparable histopathologically but the latter were larger and more aggressive and cultured cells derived from such melanomas were also larger and had higher levels of nuclear atypia. Moreover, the melanomas in Cdk4(R24C/R24C)/TPras mice, but not in TPras mice, readily metastasized to regional lymph nodes. Thus, it seems that in the mouse, Hras activation initiates UVR-induced melanoma development whereas the cell cycle defect introduced by mutant Cdk4 contributes to tumor progression, producing more aggressive, metastatic tumors.