Error-corrected next generation sequencing - Promises and challenges for genotoxicity and cancer risk assessment.

Error-corrected Next Generation Sequencing (ecNGS) is rapidly emerging as a valuable, highly sensitive and accurate method for detecting and characterizing mutations in any cell type, tissue or organism from which DNA can be isolated. Recent mutagenicity and carcinogenicity studies have used ecNGS to quantify drug-/chemical-induced mutations and mutational spectra associated with cancer risk. ecNGS has potential applications in genotoxicity assessment as a new readout for traditional models, for mutagenesis studies in 3D organotypic cultures, and for detecting off-target effects of gene editing tools. Additionally, early data suggest that ecNGS can measure clonal expansion of mutations as a mechanism-agnostic early marker of carcinogenic potential and can evaluate mutational load directly in human biomonitoring studies. In this review, we discuss promising applications, challenges, limitations, and key data initiatives needed to enable regulatory testing and adoption of ecNGS - including for advancing safety assessment, augmenting weight-of-evidence for mutagenicity and carcinogenicity mechanisms, identifying early biomarkers of cancer risk, and managing human health risk from chemical exposures.

A Brief Practical Guide to PCR.

Polymerase chain reaction (PCR) has been a powerful molecular biology tool since the mid-1980s. Millions of copies of specific sequence regions of DNA can be generated to allow the study of these regions. Fields that use this technology range from forensics to the experimental study of human biology. Standards for performing PCR and information tools to help design PCR protocols aid in successful implementation of PCR.

Assessment of Clonal Expansion Using CarcSeq Measurement of Lung Cancer Driver Mutations and Correlation With Mouse Strain- and Sex-Related Incidence of Spontaneous Lung Neoplasia.

Quantification of variation in levels of spontaneously occurring cancer driver mutations (CDMs) was developed to assess clonal expansion and predict future risk of neoplasm development. Specifically, an error-corrected next-generation sequencing method, CarcSeq, and a mouse CarcSeq panel (analogous to human and rat panels) were developed and used to quantify low-frequency mutations in a panel of amplicons enriched in hotspot CDMs. Mutations in a subset of panel amplicons, Braf, Egfr, Kras, Stk11, and Tp53, were related to incidence of lung neoplasms at 2 years. This was achieved by correlating median absolute deviation (MAD) from the overall median mutant fraction (MF) measured in the lung DNA of 16-week-old male and female, B6C3F1 and CD-1 mice (10 mice/sex/strain) with percentages of spontaneous alveolar/bronchioloalveolar adenomas and carcinomas reported in bioassay control groups. A total of 1586 mouse lung mutants with MFs >1 × 10-4 were recovered. The ratio of nonsynonymous to synonymous mutations was used to assess the proportion of recovered mutations conferring a positive selective advantage. The greatest ratio was observed in what is considered the most lung tumor-sensitive model examined, male B6C3F1 mice. Of the recurrent, nonsynonymous mouse mutations recovered, 55.5% have been reported in human tumors, with many located in or around the mouse equivalent of human cancer hotspot codons. MAD for the same subset of amplicons measured in normal human lung DNA samples showed a correlation of moderate strength and borderline significance with age (a cancer risk factor), as well as age-related cumulative lung cancer risk, suggesting MAD may inform species extrapolation.

CarcSeq Measurement of Rat Mammary Cancer Driver Mutations and Relation to Spontaneous Mammary Neoplasia.

The ability to deduce carcinogenic potential from subchronic, repeat dose rodent studies would constitute a major advance in chemical safety assessment and drug development. This study investigated an error-corrected NGS method (CarcSeq) for quantifying cancer driver mutations (CDMs) and deriving a metric of clonal expansion predictive of future neoplastic potential. CarcSeq was designed to interrogate subsets of amplicons encompassing hotspot CDMs applicable to a variety of cancers. Previously, normal human breast DNA was analyzed by CarcSeq and metrics based on mammary-specific CDMs were correlated with tissue donor age, a surrogate of breast cancer risk. Here we report development of parallel methodologies for rat. The utility of the rat CarcSeq method for predicting neoplastic potential was investigated by analyzing mammary tissue of 16-week-old untreated rats with known differences in spontaneous mammary neoplasia (Fischer 344, Wistar Han, and Sprague Dawley). Hundreds of mutants with mutant fractions ≥ 10-4 were quantified in each strain, most were recurrent mutations, and 42.5% of the nonsynonymous mutations have human homologs. Mutants in the mammary-specific target of the most tumor-sensitive strain (Sprague Dawley) showed the greatest nonsynonymous/synonymous mutation ratio, indicative of positive selection consistent with clonal expansion. For the mammary-specific target (Hras, Pik3ca, and Tp53 amplicons), median absolute deviation correlated with percentages of rats that develop spontaneous mammary neoplasia at 104 weeks (Pearson r = 1.0000, 1-tailed p = .0010). Therefore, this study produced evidence CarcSeq analysis of spontaneously occurring CDMs can be used to derive an early metric of clonal expansion relatable to long-term neoplastic outcome.

Quantification of cancer driver mutations in human breast and lung DNA using targeted, error-corrected CarcSeq.

There is a need for scientifically-sound, practical approaches to improve carcinogenicity testing. Advances in DNA sequencing technology and knowledge of events underlying cancer development have created an opportunity for progress in this area. The long-term goal of this work is to develop variation in cancer driver mutation (CDM) levels as a metric of clonal expansion of cells carrying CDMs because these important early events could inform carcinogenicity testing. The first step toward this goal was to develop and validate an error-corrected next-generation sequencing method to analyze panels of hotspot cancer driver mutations (hCDMs). The "CarcSeq" method that was developed uses unique molecular identifier sequences to construct single-strand consensus sequences for error correction. CarcSeq was used for mutational analysis of 13 amplicons encompassing >20 hotspot CDMs in normal breast, normal lung, ductal carcinomas, and lung adenocarcinomas. The approach was validated by detecting expected differences related to tissue type (normal vs. tumor and breast vs. lung) and mutation spectra. CarcSeq mutant fractions (MFs) correlated strongly with previously obtained ACB-PCR mutant fraction (MF) measurements from the same samples. A reconstruction experiment, in conjunction with other analyses, showed CarcSeq accurately quantifies MFs ≥10-4 . CarcSeq MF measurements were correlated with tissue donor age and breast cancer risk. CarcSeq MF measurements were correlated with deviation from median MFs analyzed to assess clonal expansion. Thus, CarcSeq is a promising approach to advance cancer risk assessment and carcinogenicity testing practices. Paradigms that should be investigated to advance this strategy for carcinogenicity testing are proposed.

Outgrowth of erlotinib-resistant subpopulations recapitulated in patient-derived lung tumor spheroids and organoids.

A model that recapitulates development of acquired therapeutic resistance is needed to improve oncology drug development and patient outcomes. To achieve this end, we established methods for the preparation and growth of spheroids from primary human lung adenocarcinomas, including methods to culture, passage, monitor growth, and evaluate changes in mutational profile over time. Primary lung tumor spheroids were cultured in Matrigel® with varying concentrations of erlotinib, a small molecule kinase inhibitor of epidermal growth factor receptor (EGFR) that is ineffective against KRAS mutant cells. Subtle changes in spheroid size and number were observed within the first two weeks of culture. Spheroids were cultured for up to 24 weeks, during which time interactions between different cell types, movement, and assembly into heterogeneous organoid structures were documented. Allele-specific competitive blocker PCR (ACB-PCR) was used to quantify low frequency BRAF V600E, KRAS G12D, KRAS G12V, and PIK3CA H1047R mutant subpopulations in tumor tissue residue (TR) samples and cultured spheroids. Mutant subpopulations, including multiple mutant subpopulations, were quite prevalent. Twelve examples of mutant enrichment were found in eight of the 14 tumors analyzed, based on the criteria that a statistically-significant increase in mutant fraction was observed relative to both the TR and the no-erlotinib control. Of the mutants quantified in erlotinib-treated cultures, PIK3CA H1047 mutant subpopulations increased most often (5/14 tumors), which is consistent with clinical observations. Thus, this ex vivo lung tumor spheroid model replicates the cellular and mutational tumor heterogeneity of human lung adenocarcinomas and can be used to assess the outgrowth of mutant subpopulations. Spheroid cultures with characterized mutant subpopulations could be used to investigate the efficacy of lung cancer combination therapies.

ACB-PCR Quantification of Low-Frequency Hotspot Cancer-Driver Mutations.

Allele-specific competitive blocker PCR (ACB-PCR) is a sensitive and quantitative approach for the selective amplification of a specific base substitution. Using the ACB-PCR technique, hotspot cancer-driver mutations (tumor-relevant mutations in oncogenes and tumor suppressor genes, which confer a selective growth advantage) are being developed as quantitative biomarkers of cancer risk. ACB-PCR employs a mutant-specific primer (with a 3'-penultimate mismatch relative to the mutant DNA sequence, but a double 3'-terminal mismatch relative to the wild-type DNA sequence) to selectively amplify rare mutant DNA molecules. A blocker primer having a non-extendable 3'-end and a 3'-penultimate mismatch relative to the wild-type DNA sequence, but a double 3'-terminal mismatch relative to the mutant DNA sequence is included in ACB-PCR to selectively repress amplification from abundant wild-type molecules. Consequently, ACB-PCR can quantify the level of a single base pair substitution mutation in a DNA population when present at a mutant:wild-type ratio of 1 × 10-5 or greater. Quantification of rare mutant alleles is achieved by parallel analysis of unknown samples and mutant fraction (MF) standards (defined mixtures of mutant and wild-type DNA sequences). The ability to quantify specific mutations with known association to cancer has several important applications in evaluating the carcinogenic potential of chemical exposures in rodent models. Further, the measurement of cancer-driver mutant subpopulations is important for precision cancer treatment (selecting the most appropriate targeted therapy and predicting the development of therapeutic resistance). This chapter provides a step-by-step description of the ACB-PCR methodology as it has been used to measure human PIK3CA codon 1047, CAT→CGT (H1047R) mutation.

Rationale and Roadmap for Developing Panels of Hotspot Cancer Driver Gene Mutations as Biomarkers of Cancer Risk.

Cancer driver mutations (CDMs) are necessary and causal for carcinogenesis and have advantages as reporters of carcinogenic risk. However, little progress has been made toward developing measurements of CDMs as biomarkers for use in cancer risk assessment. Impediments for using a CDM-based metric to inform cancer risk include the complexity and stochastic nature of carcinogenesis, technical difficulty in quantifying low-frequency CDMs, and lack of established relationships between cancer driver mutant fractions and tumor incidence. Through literature review and database analyses, this review identifies the most promising targets to investigate as biomarkers of cancer risk. Mutational hotspots were discerned within the 20 most mutated genes across the 10 deadliest cancers. Forty genes were identified that encompass 108 mutational hotspot codons overrepresented in the COSMIC database; 424 different mutations within these hotspot codons account for approximately 63,000 tumors and their prevalence across tumor types is described. The review summarizes literature on the prevalence of CDMs in normal tissues and suggests such mutations are direct and indirect substrates for chemical carcinogenesis, which occurs in a spatially stochastic manner. Evidence that hotspot CDMs (hCDMs) frequently occur as tumor subpopulations is presented, indicating COSMIC data may underestimate mutation prevalence. Analyses of online databases show that genes containing hCDMs are enriched in functions related to intercellular communication. In its totality, the review provides a roadmap for the development of tissue-specific, CDM-based biomarkers of carcinogenic potential, comprised of batteries of hCDMs and can be measured by error-correct next-generation sequencing. Environ. Mol. Mutagen. 61:152-175, 2020. Published 2019. This article is a U.S. Government work and is in the public domain in the USA. Environmental and Molecular Mutagenesis published by Wiley Periodicals, Inc. on behalf of Environmental Mutagen Society.

Low-Frequency Mutational Heterogeneity of Invasive Ductal Carcinoma Subtypes: Information to Direct Precision Oncology.

Information regarding the role of low-frequency hotspot cancer-driver mutations (CDMs) in breast carcinogenesis and therapeutic response is limited. Using the sensitive and quantitative Allele-specific Competitor Blocker PCR (ACB-PCR) approach, mutant fractions (MFs) of six CDMs (PIK3CA H1047R and E545K, KRAS G12D and G12V, HRAS G12D, and BRAF V600E) were quantified in invasive ductal carcinomas (IDCs; including ~20 samples per subtype). Measurable levels (i.e., ≥ 1 × 10-5, the lowest ACB-PCR standard employed) of the PIK3CA H1047R, PIK3CA E545K, KRAS G12D, KRAS G12V, HRAS G12D, and BRAF V600E mutations were observed in 34/81 (42%), 29/81 (36%), 51/81 (63%), 9/81 (11%), 70/81 (86%), and 48/81 (59%) of IDCs, respectively. Correlation analysis using available clinicopathological information revealed that PIK3CA H1047R and BRAF V600E MFs correlate positively with maximum tumor dimension. Analysis of IDC subtypes revealed minor mutant subpopulations of critical genes in the MAP kinase pathway (KRAS, HRAS, and BRAF) were prevalent across IDC subtypes. Few triple-negative breast cancers (TNBCs) had appreciable levels of PIK3CA mutation, suggesting that individuals with TNBC may be less responsive to inhibitors of the PI3K/AKT/mTOR pathway. These results suggest that low-frequency hotspot CDMs contribute significantly to the intertumoral and intratumoral genetic heterogeneity of IDCs, which has the potential to impact precision oncology approaches.

Fixation and Spread of Somatic Mutations in Adult Human Colonic Epithelium.

We investigated the means and timing by which mutations become fixed in the human colonic epithelium by visualizing somatic clones and mathematical inference. Fixation requires two sequential steps. First, one of approximately seven active stem cells residing within each colonic crypt has to be mutated. Second, the mutated stem cell has to replace neighbors to populate the entire crypt in a process that takes several years. Subsequent clonal expansion due to crypt fission is infrequent for neutral mutations (around 0.7% of all crypts undergo fission in a single year). Pro-oncogenic mutations subvert both stem cell replacement to accelerate fixation and clonal expansion by crypt fission to achieve high mutant allele frequencies with age. The benchmarking of these behaviors allows the advantage associated with different gene-specific mutations to be compared irrespective of the cellular mechanisms by which they are conferred.

Variation in organ-specific PIK3CA and KRAS mutant levels in normal human tissues correlates with mutation prevalence in corresponding carcinomas.

Large-scale sequencing efforts have described the mutational complexity of individual cancers and identified mutations prevalent in different cancers. As a complementary approach, allele-specific competitive blocker PCR (ACB-PCR) is being used to quantify levels of hotspot cancer driver mutations (CDMs) with high sensitivity, to elucidate the tissue-specific properties of CDMs, their occurrence as tumor cell subpopulations, and their occurrence in normal tissues. Here we report measurements of PIK3CA H1047R mutant fraction (MF) in normal colonic mucosa, normal lung, colonic adenomas, colonic adenocarcinomas, and lung adenocarcinomas. We report PIK3CA E545K MF measurements in those tissues, as well as in normal breast, normal thyroid, mammary ductal carcinomas, and papillary thyroid carcinomas. We report KRAS G12D and G12V MF measurements in normal colon. These MF measurements were integrated with previously published ACB-PCR data on KRAS G12D, KRAS G12V, and PIK3CA H1047R. Analysis of these data revealed a correlation between the degree of interindividual variability in these mutations (as log10 MF standard deviation) in normal tissues and the frequencies with which the mutations are detected in carcinomas of the corresponding organs in the COSMIC database. This novel observation has important implications. It suggests that interindividual variability in mutation levels of normal tissues may be used as a metric to identify mutations with critical early roles in tissue-specific carcinogenesis. Additionally, it raises the possibility that personalized cancer therapeutics, developed to target specifically activated oncogenic products, might be repurposed as prophylactic therapies to reduce the accumulation of cells carrying CDMs and, thereby, reduce future cancer risk. Environ. Mol. Mutagen. 58:466-476, 2017. © 2017 This article is a U.S. Government work and is in the public domain in the USA. Environmental and Molecular Mutagenesis published by Wiley Periodicals, Inc. on behalf of Environmental Mutagen Society.

Ovarian effects of prenatal exposure to benzo[a]pyrene: Roles of embryonic and maternal glutathione status.

Females deficient in the glutamate cysteine ligase modifier subunit (Gclm) of the rate-limiting enzyme in glutathione synthesis are more sensitive to ovarian follicle depletion and tumorigenesisby prenatal benzo[a]pyrene (BaP) exposure than Gclm+/+ mice. We investigated effects of prenatal exposure to BaP on reproductive development and ovarian mutations in Kras, a commonly mutated gene in epithelial ovarian tumors. Pregnantmice were dosed from gestational day 6.5 through 15.5 with 2mg/kg/day BaP or vehicle. Puberty onset occurred 5 days earlier in F1 daughters of all Gclm genotypes exposed to BaP compared to controls. Gclm+/- F1 daughters of Gclm+/- mothers and wildtype F1 daughters of wildtype mothers had similar depletion of ovarian follicles following prenatal exposure to BaP, suggesting that maternal Gclm genotype does not modify ovarian effects of prenatal BaP. We observed no BaP treatment or Gclm genotype related differences in ovarian Kras codon 12 mutations in F1 offspring.

Dose and temporal evaluation of ethylene oxide-induced mutagenicity in the lungs of male big blue mice following inhalation exposure to carcinogenic concentrations.

Ethylene oxide (EO) is a direct acting alkylating agent; in vitro and in vivo studies indicate that it is both a mutagen and a carcinogen. However, it remains unclear whether the mode of action (MOA) for cancer for EO is a mutagenic MOA, specifically via point mutation. To investigate the MOA for EO-induced mouse lung tumors, male Big Blue (BB) B6C3F1 mice (10/group) were exposed to EO by inhalation, 6 hr/day, 5 days/week for 4 (0, 10, 50, 100, or 200 ppm EO), 8, or 12 weeks (0, 100, or 200 ppm EO). Lung DNA samples were analyzed for cII mutant frequency (MF) at 4, 8 and 12 weeks of exposure; the mutation spectrum was analyzed for mutants from control and 200 ppm EO treatments. Although EO-induced cII MFs were 1.5- to 2.7-fold higher than the concurrent controls at 4 weeks, statistically significant increases in the cII MF were found only after 8 and 12 weeks of exposure and only at 200 ppm EO (P ≤ 0.05), which is twice the highest concentration used in the cancer bioassay. Consistent with the positive response, DNA sequencing of cII mutants showed a significant shift in the mutational spectra between control and 200 ppm EO following 8 and 12 week exposures (P ≤ 0.035), but not at 4 weeks. Thus, EO mutagenic activity in vivo was relatively weak and required higher than tumorigenic concentrations and longer than 4 weeks exposure durations. These data do not follow the classical patterns for a MOA mediated by point mutations. Environ. Mol. Mutagen. 58:122-134, 2017. © 2017 Wiley Periodicals, Inc.

Breast Cancer Heterogeneity Examined by High-Sensitivity Quantification of PIK3CA, KRAS, HRAS, and BRAF Mutations in Normal Breast and Ductal Carcinomas.

Mutant cancer subpopulations have the potential to derail durable patient responses to molecularly targeted cancer therapeutics, yet the prevalence and size of such subpopulations are largely unexplored. We employed the sensitive and quantitative Allele-specific Competitive Blocker PCR approach to characterize mutant cancer subpopulations in ductal carcinomas (DCs), examining five specific hotspot point mutations (PIK3CA H1047R, KRAS G12D, KRAS G12V, HRAS G12D, and BRAF V600E). As an approach to aid interpretation of the DC results, the mutations were also quantified in normal breast tissue. Overall, the mutations were prevalent in normal breast and DCs, with 9/9 DCs having measureable levels of at least three of the five mutations. HRAS G12D was significantly increased in DCs as compared to normal breast. The most frequent point mutation reported in DC by DNA sequencing, PIK3CA H1047R, was detected in all normal breast tissue and DC samples and was present at remarkably high levels (mutant fractions of 1.1 × 10(-3) to 4.6 × 10(-2)) in 4/10 normal breast samples. In normal breast tissue samples, PIK3CA mutation levels were positively correlated with age. However, the PIK3CA H1047R mutant fraction distributions for normal breast tissues and DCs were similar. The results suggest PIK3CA H1047R mutant cells have a selective advantage in breast, contribute to breast cancer susceptibility, and drive tumor progression during breast carcinogenesis, even when present as only a subpopulation of tumor cells.

Targeted therapies with companion diagnostics in the management of breast cancer: current perspectives.

Breast cancer is a multifaceted disease exhibiting both intertumoral and intratumoral heterogeneity as well as variable disease course. Over 2 decades of research has advanced the understanding of the molecular substructure of breast cancer, directing the development of new therapeutic strategies against these actionable targets. In vitro diagnostics, and specifically companion diagnostics, have been integral in the successful development and implementation of these targeted therapies, such as those directed against the human epidermal growth factor receptor 2. Lately, there has been a surge in the development, commercialization, and marketing of diagnostic assays to assist in breast cancer patient care. More recently, multigene signature assays, such as Oncotype DX, MammaPrint, and Prosigna, have been integrated in the clinical setting in order to tailor decisions on adjuvant endocrine and chemotherapy treatment. This review provides an overview of the current state of breast cancer management and the use of companion diagnostics to direct personalized approaches in the treatment of breast cancer.

Low-frequency KRAS mutations are prevalent in lung adenocarcinomas.

AIM: This study quantified low-frequency KRAS mutations in normal lung and lung adenocarcinomas, to understand their potential significance in the development of acquired resistance to EGFR-targeted therapies. MATERIALS & METHODS: Allele-specific Competitive Blocker-PCR was used to quantify KRAS codon 12 GAT (G12D) and GTT (G12V) mutation in 19 normal lung and 21 lung adenocarcinoma samples. RESULTS: Lung adenocarcinomas had KRAS codon 12 GAT and GTT geometric mean mutant fractions of 1.94 × 10-4 and 1.16 × 10-3, respectively. For 76.2% of lung adenocarcinomas, the level of KRAS mutation was greater than the upper 95% confidence interval of that in normal lung. CONCLUSION: KRAS mutant tumor subpopulations, not detectable by DNA sequencing, may drive resistance to EGFR blockade in lung adenocarcinoma patients.

ACB-PCR quantification of somatic oncomutation.

Allele-specific competitive blocker-polymerase chain reaction (ACB-PCR) is a sensitive approach for the selective amplification of an allele. Using the ACB-PCR technique, hotspot point mutations in oncogenes and tumor-suppressor genes (oncomutations) are being developed as quantitative biomarkers of cancer risk. ACB-PCR employs a mutant specific primer (with a 3'-penultimate mismatch relative to the mutant DNA sequence, but a double 3'-terminal mismatch relative to the wild-type DNA sequence) to selectively amplify rare mutant DNA molecules. A blocker primer (having a non-extendable 3'-end and with a 3'-penultimate mismatch relative to the wild-type DNA sequence, but a double 3'-terminal mismatch relative to the mutant DNA sequence) is included in ACB-PCR to selectively repress amplification from the abundant wild-type molecules. Consequently, ACB-PCR is capable of quantifying the level of a single basepair substitution mutation in a DNA population when present at a mutant:wild type ratio of 10(-5) or greater. Quantification of rare mutant alleles is achieved by parallel analysis of unknown samples and mutant fraction (MF) standards (defined mixtures of mutant and wild-type DNA sequences). The ability to quantify specific mutations with known association to cancer has several important applications, including evaluating the carcinogenic potential of chemical exposures in rodent models and in the diagnosis and treatment of cancer. This chapter provides a step-by-step description of the ACB-PCR methodology as it has been used to measure human KRAS codon 12 GGT to GAT mutation.

A subset of papillary thyroid carcinomas contain KRAS mutant subpopulations at levels above normal thyroid.

The molecular pathogenesis of papillary thyroid carcinoma (PTC) is largely attributed to chromosomal rearrangements and point mutations in genes within the MAPK pathway (i.e., BRAF and RAS). Despite KRAS being the 6th most frequently mutated gene for all cancers, the reported frequency in thyroid cancer is only 2%. This may be due, in part, to the use of insensitive mutation detection methods such as DNA sequencing. Therefore, using the sensitive and quantitative ACB-PCR approach, we quantified KRAS codon 12 GGT → GAT and GGT → GTT mutant fraction (MF) in 20 normal thyroid tissues, 17 primary PTC, 2 metastatic PTC, and 1 anaplastic thyroid carcinoma. We observed measurable levels of KRAS codon 12 GAT or GTT mutation in all normal thyroid tissues. For PTCs, 29.4% and 35.3% had KRAS codon 12 GAT and GTT MF above the 95% upper confidence interval for the corresponding MFs in normal thyroid. The highest observed KRAS codon 12 GTT MFs were associated with tumors with follicular characteristics and relatively high levels of tumor necrosis. The results indicate KRAS mutant subpopulations are present in a large number of thyroid tumors, a fact previously unrecognized. The presence of KRAS mutation may indicate a tumor with an aggressive phenotype, thus directing the course of clinical treatment.

Temporal changes in K-ras mutant fraction in lung tissue of big blue B6C3F1 mice exposed to ethylene oxide.

Ethylene oxide (EO) is a genotoxicant and a mouse lung carcinogen, but whether EO is carcinogenic through a mutagenic mode of action remains unclear. To investigate this question, 8-week-old male Big Blue B6C3F1 mice (10 mice/group) were exposed to EO by inhalation-6 h/day, 5 days/week for 4 weeks (0, 10, 50, 100, or 200 ppm EO) and 8 or 12 weeks (0, 100, or 200 ppm EO). Lung DNA samples were analyzed for levels of 3 K-ras codon 12 mutations (GGT→GAT, GGT→GTT, and GGT→TGT) using ACB-PCR. No measureable level of K-ras codon 12 TGT mutation was detected (ie, all lung mutant fractions [MFs] ≤ 10⁻5). Four weeks of inhalation of 100 ppm EO caused a significant increase in K-ras codon 12 GGT→GTT MF relative to controls, whereas 50, 100, and 200 ppm EO caused significant increases in K-ras codon 12 GGT→GAT MF. In addition, significant inverse correlations were observed between K-ras codon 12 GGT→GTT MF and cII mutant frequency in the lungs of the same mice exposed to 50, 100, or 200 ppm EO for 4 weeks. Surprisingly, 8 weeks of exposure to 100 and 200 ppm EO caused significant decreases in K-ras MFs relative to controls. Thus, the changes in K-ras MF as a function of cumulative EO dose were nonmonotonic and were consistent with EO causing early amplification of preexisting K-ras mutations, rather than induction of K-ras mutation through genotoxicity at codon 12. The possibility that these changes reflect K-ras mutant cell selection under varying degrees of oxidative stress is discussed.