Single Circulating-Tumor-Cell-Targeted Sequencing to Identify Somatic Variants in Liquid Biopsies in Non-Small-Cell Lung Cancer Patients.

Non-small-cell lung cancer (NSCLC) accounts for most cancer-related deaths worldwide. Liquid biopsy by a blood draw to detect circulating tumor cells (CTCs) is a tool for molecular profiling of cancer using single-cell and next-generation sequencing (NGS) technologies. The aim of the study was to identify somatic variants in single CTCs isolated from NSCLC patients by targeted NGS. Thirty-one subjects (20 NSCLC patients, 11 smokers without cancer) were enrolled for blood draws (7.5 mL). CTCs were identified by immunofluorescence, individually retrieved, and DNA-extracted. Targeted NGS was performed to detect somatic variants (single-nucleotide variants (SNVs) and insertions/deletions (Indels)) across 65 oncogenes and tumor suppressor genes. Cancer-associated variants were classified using OncoKB database. NSCLC patients had significantly higher CTC counts than control smokers (p = 0.0132; Mann-Whitney test). Analyzing 23 CTCs and 13 white blood cells across seven patients revealed a total of 644 somatic variants that occurred in all CTCs within the same subject, ranging from 1 to 137 per patient. The highest number of variants detected in ≥1 CTC within a patient was 441. A total of 18/65 (27.7%) genes were highly mutated. Mutations with oncogenic impact were identified in functional domains of seven oncogenes/tumor suppressor genes (NF1, PTCH1, TP53, SMARCB1, SMAD4, KRAS, and ERBB2). Single CTC-targeted NGS detects heterogeneous and shared mutational signatures within and between NSCLC patients. CTC single-cell genomics have potential for integration in NSCLC precision oncology.

Beyond Circulating Tumor Cell Enumeration: Cell-Based Liquid Biopsy to Assess Protein Biomarkers and Cancer Genomics Using the RareCyte® Platform.

The practice of medicine has steadily employed less invasive methods to obtain information derived from the tumor to guide clinical management of patients. Liquid biopsy-the sampling of blood-is a non-invasive method for generating information previously only available from tissue biopsies of the tumor mass. Analysis of fragmented circulating tumor DNA in the plasma is clinically used to identify actionable mutations and detect residual or recurrent disease. Plasma analysis cannot, however, assess cancer phenotypes, including the expression of drug targets and protein biomarkers. Circulating tumor cells (CTCs) are intact cancer cells that have entered the blood that have the potential for distant metastasis. While enumeration of CTCs is prognostic of outcome, recently developed technology allows for the interrogation of protein biomarkers on CTCs that could be predictive of response. Furthermore, since CTCs contain intact whole cancer genomes, isolating viable CTCs detected during therapy could provide a rational approach to assessing mutational profiles of resistance. Identification, characterization and molecular analysis of CTCs together will advance the capacity of liquid biopsy to meet the requirements of twenty-first century medicine.

Identification and characterization of effusion tumor cells (ETCs) from remnant pleural effusion specimens.

BACKGROUND: Cancer is a leading cause of death worldwide, and patients may have advanced disease when diagnosed. Targeted therapies guided by molecular subtyping of cancer can benefit patients significantly. Pleural effusions are frequently observed in patients with metastatic cancer and are routinely removed for therapeutic purposes; however, effusion specimens have not been recognized as typical substrates for clinical molecular testing because of frequent low tumor cellularity. METHODS: Excess remnant pleural effusion samples (N = 25) from 21 patients with and without suspected malignancy were collected at Stanford Health Care between December 2019 and November 2020. Samples were processed into ThinPrep slides and underwent novel effusion tumor cell (ETC) analysis. The ETC results were compared with the original clinical diagnoses for accuracy. A subset of confirmed ETCs was further isolated and processed for molecular profiling to identify cancer driver mutations. All samples were obtained with Institutional Review Board approval. RESULTS: The authors established novel quantitative standards to identify ETCs and detected epithelial malignancy with 89.5% sensitivity and 100% specificity in the pleural effusion samples. Molecular profiling of confirmed ETCs (pools of 5 cells evaluated) revealed key pathogenic mutations consistent with clinical molecular findings. CONCLUSIONS: In this study, the authors developed a novel ETC-testing assay that detected epithelial malignancies in pleural effusions with high sensitivity and specificity. Molecular profiling of 5 ETCs showed promising concordance with the clinical molecular findings. To promote cancer subtyping and guide treatment, this ETC-testing assay will need to be validated in larger patient cohorts to facilitate integration into cytologic workflow.

Lower mitochondrial DNA content but not increased mutagenesis associates with decreased base excision repair activity in brains of AD subjects.

Accumulation of oxidative mitochondrial DNA (mtDNA) damage and impaired base excision repair (BER) in brains have been associated with Alzheimer's disease (AD). However, it is still not clear how these affect mtDNA stability, as reported levels of mtDNA mutations in AD are conflicting. Thus, we investigated whether alterations in BER correlate with mtDNA instability in AD using postmortem brain samples from cognitively normal AD subjects and individuals who show neuropathological features of AD, but remained cognitively normal (high-pathology control). To date, no data on DNA repair and mtDNA stability are available for these individuals. BER activities, mtDNA mutations, and mtDNA copy number were measured in the nuclear and mitochondrial extracts. Significantly lower uracil DNA glycosylase activity was detected in nuclear and mitochondrial extracts from AD subjects, while apurinic/apyrimidinic endonuclease activity was similar in all groups. Although mtDNA mutation frequency was similar in all groups, mtDNA copy number was significantly decreased in the temporal cortex of AD brains but not of high-pathology control subjects. Our results show that lower mitochondrial uracil DNA glycosylase activity does not result in increased mutagenesis, but rather in depletion of mtDNA in early-affected brain regions during AD development.

Multiple Molecular Mechanisms Rescue mtDNA Disease in C. elegans.

Genetic instability of the mitochondrial genome (mtDNA) plays an important role in human aging and disease. Thus far, it has proven difficult to develop successful treatment strategies for diseases that are caused by mtDNA instability. To address this issue, we developed a model of mtDNA disease in the nematode C. elegans, an animal model that can rapidly be screened for genes and biological pathways that reduce mitochondrial pathology. These worms recapitulate all the major hallmarks of mtDNA disease in humans, including increased mtDNA instability, loss of respiration, reduced neuromuscular function, and a shortened lifespan. We found that these phenotypes could be rescued by intervening in numerous biological pathways, including IGF-1/insulin signaling, mitophagy, and the mitochondrial unfolded protein response, suggesting that it may be possible to ameliorate mtDNA disease through multiple molecular mechanisms.

Cyclophosphamide leads to persistent deficits in physical performance and in vivo mitochondria function in a mouse model of chemotherapy late effects.

Fatigue is the symptom most commonly reported by long-term cancer survivors and is increasingly recognized as related to skeletal muscle dysfunction. Traditional chemotherapeutic agents can cause acute toxicities including cardiac and skeletal myopathies. To investigate the mechanism by which chemotherapy may lead to persistent skeletal muscle dysfunction, mature adult mice were injected with a single cyclophosphamide dose and evaluated for 6 weeks. We found that exposed mice developed a persistent decrease in treadmill running time compared to baseline (25.7±10.6 vs. 49.0±16.8 min, P = 0.0012). Further, 6 weeks after drug exposure, in vivo parameters of mitochondrial function remained below baseline including maximum ATP production (482.1 ± 48.6 vs. 696.2 ± 76.6, P = 0.029) and phosphocreatine to ATP ratio (3.243 ± 0.1 vs. 3.878 ± 0.1, P = 0.004). Immunoblotting of homogenized muscles from treated animals demonstrated a transient increase in HNE adducts 1 week after exposure that resolved by 6 weeks. However, there was no evidence of an oxidative stress response as measured by quantitation of SOD1, SOD2, and catalase protein levels. Examination of mtDNA demonstrated that the mutation frequency remained comparable between control and treated groups. Interestingly, there was evidence of a transient increase in NF-ĸB p65 protein 1 day after drug exposure as compared to saline controls (0.091±0.017 vs. 0.053±0.022, P = 0.033). These data suggest that continued impairment in muscle and mitochondria function in cyclophosphamide-treated animals is not linked to persistent oxidative stress and that alternative mechanisms need to be considered.

Age-related accumulation of phosphorylated mitofusin 2 protein in retinal ganglion cells correlates with glaucoma progression.

Dysregulation of axonal bioenergetics is likely a key mechanism in the initiation and progression of age-related neurodegenerative diseases. Glaucoma is a quintessential neurodegenerative disorder characterized by progressive deterioration of the optic nerve (ON) and eventual death of retinal ganglion cells (RGCs). Age and elevation of intraocular pressure are key risk factors in glaucoma, but the common early hallmarks of decreased axonal transport and increased bioenergetic vulnerability likely underlie disease initiation. We examined the correlation between bioenergetics and axonal transport with mitochondrial mutation frequency and post-translational modifications of mitofusin 2 (Mfn2) in RGCs during glaucoma progression. No increase in the frequency of mtDNA mutations was detected, but we observed significant shifts in mitochondrial protein species. Mfn2 is a fusion protein that functions in mitochondrial biogenesis, maintenance, and mitochondrial transport. We demonstrate that Mfn2 accumulates selectively in RGCs during glaucomatous degeneration, that two novel states of Mfn2 exist in retina and ON, and identify a phosphorylated form that selectively accumulates in RGCs, but is absent in ON. Phosphorylation of Mfn2 is correlated with higher ubiquitination, and failure of the protein to reach the ON. Together, these data suggest that post-translational modification of Mfn2 is associated with its dysregulation during a window of metabolic vulnerability that precedes glaucomatous degeneration. Future work to either manipulate expression of Mfn2 or to prevent its degradation could have therapeutic value in the treatment of neurodegenerative diseases where long-tract axons are vulnerable.

Massively parallel digital transcriptional profiling of single cells.

Characterizing the transcriptome of individual cells is fundamental to understanding complex biological systems. We describe a droplet-based system that enables 3' mRNA counting of tens of thousands of single cells per sample. Cell encapsulation, of up to 8 samples at a time, takes place in ∼6 min, with ∼50% cell capture efficiency. To demonstrate the system's technical performance, we collected transcriptome data from ∼250k single cells across 29 samples. We validated the sensitivity of the system and its ability to detect rare populations using cell lines and synthetic RNAs. We profiled 68k peripheral blood mononuclear cells to demonstrate the system's ability to characterize large immune populations. Finally, we used sequence variation in the transcriptome data to determine host and donor chimerism at single-cell resolution from bone marrow mononuclear cells isolated from transplant patients.

Mitochondrial DNA exhibits resistance to induced point and deletion mutations.

The accumulation of somatic mitochondrial DNA (mtDNA) mutations contributes to the pathogenesis of human disease. Currently, mitochondrial mutations are largely considered results of inaccurate processing of its heavily damaged genome. However, mainly from a lack of methods to monitor mtDNA mutations with sufficient sensitivity and accuracy, a link between mtDNA damage and mutation has not been established. To test the hypothesis that mtDNA-damaging agents induce mtDNA mutations, we exposed MutaTMMouse mice to benzo[a]pyrene (B[a]P) or N-ethyl-N-nitrosourea (ENU), daily for 28 consecutive days, and quantified mtDNA point and deletion mutations in bone marrow and liver using our newly developed Digital Random Mutation Capture (dRMC) and Digital Deletion Detection (3D) assays. Surprisingly, our results demonstrate mutagen treatment did not increase mitochondrial point or deletion mutation frequencies, despite evidence both compounds increase nuclear DNA mutations and demonstrated B[a]P adduct formation in mtDNA. These findings contradict models of mtDNA mutagenesis that assert the elevated rate of mtDNA mutation stems from damage sensitivity and abridged repair capacity. Rather, our results demonstrate induced mtDNA damage does not readily convert into mutation. These findings suggest robust mitochondrial damage responses repress induced mutations after mutagen exposure.

Cells Comprising the Prostate Cancer Microenvironment Lack Recurrent Clonal Somatic Genomic Aberrations.

UNLABELLED: Prostate cancer-associated stroma (CAS) plays an active role in malignant transformation, tumor progression, and metastasis. Molecular analyses of CAS have demonstrated significant changes in gene expression; however, conflicting evidence exists on whether genomic alterations in benign cells comprising the tumor microenvironment (TME) underlie gene expression changes and oncogenic phenotypes. This study evaluates the nuclear and mitochondrial DNA integrity of prostate carcinoma cells, CAS, matched benign epithelium and benign epithelium-associated stroma by whole-genome copy-number analyses, targeted sequencing of TP53, and FISH. Array comparative genomic hybridization (aCGH) of CAS revealed a copy-neutral diploid genome with only rare and small somatic copy-number aberrations (SCNA). In contrast, several expected recurrent SCNAs were evident in the adjacent prostate carcinoma cells, including gains at 3q, 7p, and 8q, and losses at 8p and 10q. No somatic TP53 mutations were observed in CAS. Mitochondrial DNA (mtDNA) extracted from carcinoma cells and stroma identified 23 somatic mtDNA mutations in neoplastic epithelial cells, but only one mutation in stroma. Finally, genomic analyses identified no SCNAs, LOH, or copy-neutral LOH in cultured cancer-associated fibroblasts, which are known to promote prostate cancer progression in vivo IMPLICATIONS: The gene expression changes observed in prostate cancer-adjacent stroma and the attendant contribution of the stroma to the development and progression of prostate cancer are not due to frequent or recurrent genomic alterations in the TME.

Targeted single molecule mutation detection with massively parallel sequencing.

Next-generation sequencing (NGS) technologies have transformed genomic research and have the potential to revolutionize clinical medicine. However, the background error rates of sequencing instruments and limitations in targeted read coverage have precluded the detection of rare DNA sequence variants by NGS. Here we describe a method, termed CypherSeq, which combines double-stranded barcoding error correction and rolling circle amplification (RCA)-based target enrichment to vastly improve NGS-based rare variant detection. The CypherSeq methodology involves the ligation of sample DNA into circular vectors, which contain double-stranded barcodes for computational error correction and adapters for library preparation and sequencing. CypherSeq is capable of detecting rare mutations genome-wide as well as those within specific target genes via RCA-based enrichment. We demonstrate that CypherSeq is capable of correcting errors incurred during library preparation and sequencing to reproducibly detect mutations down to a frequency of 2.4 × 10(-7) per base pair, and report the frequency and spectra of spontaneous and ethyl methanesulfonate-induced mutations across the Saccharomyces cerevisiae genome.

The metabolome regulates the epigenetic landscape during naive-to-primed human embryonic stem cell transition.

For nearly a century developmental biologists have recognized that cells from embryos can differ in their potential to differentiate into distinct cell types. Recently, it has been recognized that embryonic stem cells derived from both mice and humans exhibit two stable yet epigenetically distinct states of pluripotency: naive and primed. We now show that nicotinamide N-methyltransferase (NNMT) and the metabolic state regulate pluripotency in human embryonic stem cells (hESCs).  Specifically, in naive hESCs, NNMT and its enzymatic product 1-methylnicotinamide are highly upregulated, and NNMT is required for low S-adenosyl methionine (SAM) levels and the H3K27me3 repressive state. NNMT consumes SAM in naive cells, making it unavailable for histone methylation that represses Wnt and activates the HIF pathway in primed hESCs. These data support the hypothesis that the metabolome regulates the epigenetic landscape of the earliest steps in human development.

Subacute calorie restriction and rapamycin discordantly alter mouse liver proteome homeostasis and reverse aging effects.

Calorie restriction (CR) and rapamycin (RP) extend lifespan and improve health across model organisms. Both treatments inhibit mammalian target of rapamycin (mTOR) signaling, a conserved longevity pathway and a key regulator of protein homeostasis, yet their effects on proteome homeostasis are relatively unknown. To comprehensively study the effects of aging, CR, and RP on protein homeostasis, we performed the first simultaneous measurement of mRNA translation, protein turnover, and abundance in livers of young (3 month) and old (25 month) mice subjected to 10-week RP or 40% CR. Protein abundance and turnover were measured in vivo using (2) H3 -leucine heavy isotope labeling followed by LC-MS/MS, and translation was assessed by polysome profiling. We observed 35-60% increased protein half-lives after CR and 15% increased half-lives after RP compared to age-matched controls. Surprisingly, the effects of RP and CR on protein turnover and abundance differed greatly between canonical pathways, with opposite effects in mitochondrial (mt) dysfunction and eIF2 signaling pathways. CR most closely recapitulated the young phenotype in the top pathways. Polysome profiles indicated that CR reduced polysome loading while RP increased polysome loading in young and old mice, suggesting distinct mechanisms of reduced protein synthesis. CR and RP both attenuated protein oxidative damage. Our findings collectively suggest that CR and RP extend lifespan in part through the reduction of protein synthetic burden and damage and a concomitant increase in protein quality. However, these results challenge the notion that RP is a faithful CR mimetic and highlight mechanistic differences between the two interventions.

Skeletal muscle mitochondrial DNA deletions are not increased in CuZn-superoxide dismutase deficient mice.

Mitochondrial DNA (mtDNA) deletion mutations are proposed contributors to aging-related muscle fiber loss and atrophy, but evidence of a causal role for these mutations in muscle aging is lacking. Elucidating the etiology of in vivo mtDNA deletion mutations will help to better understand and test the possible roles of these mutations in aging. The implication of mtDNA mutations in aging is based on the susceptibility of mtDNA to oxidative damage by reactive oxygen species (ROS) due to residing in mitochondria, the primary source of endogenous ROS. Cells possess many pathways for neutralizing ROSs, including a variety of superoxide dismutases (SOD). Mice lacking CuZnSOD (Sod1(-/-) mice) have high levels of oxidative damage in many tissues including skeletal muscle and are a model for testing the role of oxidative damage in the formation of mtDNA deletion mutations. The increased DNA oxidative damage in Sod1(-/-) mice is associated with increased mtDNA deletion mutations in a variety of tissues, but skeletal muscle mtDNA mutations have not been reported. We hypothesized that a life-long absence of mouse muscle CuZnSOD would increase mtDNA deletion mutation frequency and focal accumulation of these mutations in aging mouse skeletal muscle. Focal accumulations of mtDNA deletion mutations were detected by histochemical staining for cytochrome c oxidase (cytOX) activity and detection of cytOX-negative fibers, a marker of focal mtDNA mutation accumulation, within approximately 20,000 muscle fibers through a distance of 1000µm. Total DNA was extracted from intervening unstained sections and mtDNA deletion mutation frequency was measured by a droplet digital PCR. Droplet digital PCR quantification of mtDNA deletion mutations showed no difference in mtDNA deletion mutation frequency in Sod1(-/-) mouse muscle compared to wild-type mice and we observed no significant increase in the number of cytOX-negative muscle fibers, in Sod1(-/-) mice compared to wild-type mice. These data demonstrate that not all changes in cellular oxidative stress are linked to mtDNA deletion mutations and shift the focus to other etiologies for these mutations that need to be clarified to better test their possible role in aging.

Targeted enrichment and high-resolution digital profiling of mitochondrial DNA deletions in human brain.

Due largely to the inability to accurately quantify and characterize de novo deletion events, the mechanisms underpinning the pathogenic expansion of mtDNA deletions in aging and neuromuscular disorders remain poorly understood. Here, we outline and validate a new tool termed 'Digital Deletion Detection' (3D) that allows for high-resolution analysis of rare deletions occurring at frequencies as low as 1 × 10(-8) . 3D is a three-step process that includes targeted enrichment for deletion-bearing molecules, single-molecule partitioning of genomes into thousands of droplets for direct quantification via droplet digital PCR, and breakpoint characterization using massively parallel sequencing. Using 3D, we interrogated over 8 billion mitochondrial genomes to analyze the age-related dynamics of mtDNA deletions in human brain tissue. We demonstrate that the total deletion load increases with age, while the total number and diversity of unique deletions remain constant. Our data provide support for the hypothesis that expansion of pre-existing mutations is the primary factor contributing to age-related accumulation of mtDNA deletions.

Digital genomic quantification of tumor-infiltrating lymphocytes.

Infiltrating T lymphocytes are frequently found in malignant tumors and are suggestive of a host cancer immune response. Multiple independent studies have documented that the presence and quantity of tumor-infiltrating lymphocytes (TILs) are strongly correlated with increased survival. However, because of methodological factors, the exact effect of TILs on prognosis has remained enigmatic, and inclusion of TILs in standard prognostic panels has been limited. For example, some reports enumerate all CD3(+) cells, some count only cytotoxic CD8(+) T cells, and the criteria used to score tumors as TIL-positive or TIL-negative are inconsistent among studies. To address this limitation, we introduce a robust digital DNA-based assay, termed QuanTILfy, to count TILs and assess T cell clonality in tissue samples, including tumors. We demonstrate the clonal specificity of this approach by the diagnosis of T cell acute lymphoblastic leukemia and the accurate, sensitive, and highly reproducible measurement of TILs in primary and metastatic ovarian cancer. Our experiments demonstrate an association between higher TIL counts and improved survival among women with ovarian cancer, and are consistent with previous observations that the immune response against ovarian cancer is a meaningful and independent prognostic factor. Surprisingly, the TIL repertoire is diverse for all tumors in the study with no notable oligoclonal expansions. Furthermore, because variability in the measurement and characterization of TILs has limited their clinical utility as biomarkers, these results highlight the significant translational potential of a robust, standardizable DNA-based assay to assess TILs in a variety of cancer types.

Decreased mitochondrial DNA mutagenesis in human colorectal cancer.

Genome instability is regarded as a hallmark of cancer. Human tumors frequently carry clonally expanded mutations in their mitochondrial DNA (mtDNA), some of which may drive cancer progression and metastasis. The high prevalence of clonal mutations in tumor mtDNA has commonly led to the assumption that the mitochondrial genome in cancer is genetically unstable, yet this hypothesis has not been experimentally tested. In this study, we directly measured the frequency of non-clonal (random) de novo single base substitutions in the mtDNA of human colorectal cancers. Remarkably, tumor tissue exhibited a decreased prevalence of these mutations relative to adjacent non-tumor tissue. The difference in mutation burden was attributable to a reduction in C:G to T:A transitions, which are associated with oxidative damage. We demonstrate that the lower random mutation frequency in tumor tissue was also coupled with a shift in glucose metabolism from oxidative phosphorylation to anaerobic glycolysis, as compared to non-neoplastic colon. Together these findings raise the intriguing possibility that fidelity of mitochondrial genome is, in fact, increased in cancer as a result of a decrease in reactive oxygen species-mediated mtDNA damage.

Roles of DNA polymerase I in leading and lagging-strand replication defined by a high-resolution mutation footprint of ColE1 plasmid replication.

DNA polymerase I (pol I) processes RNA primers during lagging-strand synthesis and fills small gaps during DNA repair reactions. However, it is unclear how pol I and pol III work together during replication and repair or how extensive pol I processing of Okazaki fragments is in vivo. Here, we address these questions by analyzing pol I mutations generated through error-prone replication of ColE1 plasmids. The data were obtained by direct sequencing, allowing an accurate determination of the mutation spectrum and distribution. Pol I's mutational footprint suggests: (i) during leading-strand replication pol I is gradually replaced by pol III over at least 1.3 kb; (ii) pol I processing of Okazaki fragments is limited to ∼20 nt and (iii) the size of Okazaki fragments is short (∼250 nt). While based on ColE1 plasmid replication, our findings are likely relevant to other pol I replicative processes such as chromosomal replication and DNA repair, which differ from ColE1 replication mostly at the recruitment steps. This mutation footprinting approach should help establish the role of other prokaryotic or eukaryotic polymerases in vivo, and provides a tool to investigate how sequence topology, DNA damage, or interactions with protein partners may affect the function of individual DNA polymerases.

Generation, function, and prognostic utility of somatic mitochondrial DNA mutations in cancer.

Exciting new studies are increasingly strengthening the link between mitochondrial mutagenesis and tumor progression. Here we provide a comprehensive review and meta-analysis of studies reporting on mitochondrial DNA mutations in common human cancers. We discuss possible mechanisms by which mitochondrial DNA mutations may influence carcinogenesis, outline important caveats for interpreting the detected mutations--particularly differentiating causality from association--and suggest how new mutational assays may help resolve fundamental controversies in the field and delineate the origin and expansion of neoplastic cell lineages. Finally, we discuss the potential clinical utility of mtDNA mutations for improving the sensitivity of early cancer diagnosis, rapidly detecting cancer recurrence, and predicting the disease outcome.

Molecularly evolved thymidylate synthase inhibits 5-fluorodeoxyuridine toxicity in human hematopoietic cells.

Thymidylate synthase (TS) inhibitors, such as 5-fluorouracil (5-FU) and 5-fluorodeoxyuridine (5-FUdR), are amongst the most frequently used chemotherapeutic drugs available, although their efficacy is often limited by myelotoxicity. An emerging strategy for overcoming bone marrow toxicity involves ex vivo genetic transfer of drug resistance to autologous hematopoietic progenitor cells, followed by reimplantation of the transfected cells before chemotherapy. Here we establish that expression of mutant TS genes, selected from millions of engineered variants, renders human hematopoietic cells resistant to 5-FUdR, and identify the most efficacious variant for gene therapeutic rescue of drug-induced myelosuppression.