MADDD-seq, a novel massively parallel sequencing tool for simultaneous detection of DNA damage and mutations.

Our genome is exposed to a wide variety of DNA-damaging agents. If left unrepaired, this damage can be converted into mutations that promote carcinogenesis or the development of genetically inherited diseases. As a result, researchers and clinicians require tools that can detect DNA damage and mutations with exceptional sensitivity. In this study, we describe a massively parallel sequencing tool termed Mutation And DNA Damage Detection-seq (MADDD-seq) that is capable of detecting O6-methyl guanine lesions and mutations simultaneously, with a single assay. To illustrate the dual capabilities of MADDD-seq, we treated WT and DNA repair deficient yeast cells with the DNA-damaging agent MNNG and tracked DNA lesions and mutations over a 24-h time period. This approach allowed us to identify thousands of DNA adducts and mutations in a single sequencing run and gain deep insight into the kinetics of DNA repair and mutagenesis.

Defining cellular population dynamics at single-cell resolution during prostate cancer progression.

Advanced prostate malignancies are a leading cause of cancer-related deaths in men, in large part due to our incomplete understanding of cellular drivers of disease progression. We investigate prostate cancer cell dynamics at single-cell resolution from disease onset to the development of androgen independence in an in vivo murine model. We observe an expansion of a castration-resistant intermediate luminal cell type that correlates with treatment resistance and poor prognosis in human patients. Moreover, transformed epithelial cells and associated fibroblasts create a microenvironment conducive to pro-tumorigenic immune infiltration, which is partially androgen responsive. Androgen-independent prostate cancer leads to significant diversification of intermediate luminal cell populations characterized by a range of androgen signaling activity, which is inversely correlated with proliferation and mRNA translation. Accordingly, distinct epithelial populations are exquisitely sensitive to translation inhibition, which leads to epithelial cell death, loss of pro-tumorigenic signaling, and decreased tumor heterogeneity. Our findings reveal a complex tumor environment largely dominated by castration-resistant luminal cells and immunosuppressive infiltrates.

An In Vivo Model of Human Macrophages in Metastatic Melanoma.

Despite recent therapeutic progress, advanced melanoma remains lethal for many patients. The composition of the immune tumor microenvironment (TME) has decisive impacts on therapy response and disease outcome, and high-dimensional analyses of patient samples reveal the heterogeneity of the immune TME. Macrophages infiltrate TMEs and generally associate with tumor progression, but the underlying mechanisms are incompletely understood. Because experimental systems are needed to elucidate the functional properties of these cells, we developed a humanized mouse model reconstituted with human immune cells and human melanoma. We used two strains of recipient mice, supporting or not supporting the development of human myeloid cells. We found that human myeloid cells favored metastatic spread of the primary tumor, thereby recapitulating the cancer-supportive role of macrophages. We next analyzed the transcriptome of human immune cells infiltrating tumors versus other tissues. This analysis identified a cluster of myeloid cells present in the TME, but not in other tissues, which do not correspond to canonical M2 cells. The transcriptome of these cells is characterized by high expression of glycolytic enzymes and multiple chemokines and by low expression of gene sets associated with inflammation and adaptive immunity. Compared with humanized mouse results, we found transcriptionally similar myeloid cells in patient-derived samples of melanoma and other cancer types. The humanized mouse model described here thus complements patient sample analyses, enabling further elucidation of fundamental principles in melanoma biology beyond M1/M2 macrophage polarization. The model can also support the development and evaluation of candidate antitumor therapies.

Astrocytic laminin-211 drives disseminated breast tumor cell dormancy in brain.

Although dormancy is thought to play a key role in the metastasis of breast tumor cells to the brain, our knowledge of the molecular mechanisms regulating disseminated tumor cell (DTC) dormancy in this organ is limited. Here using serial intravital imaging of dormant and metastatic triple-negative breast cancer lines, we identify escape from the single-cell or micrometastatic state as the rate-limiting step towards brain metastasis. We show that every DTC occupies a vascular niche, with quiescent DTCs residing on astrocyte endfeet. At these sites, astrocyte-deposited laminin-211 drives DTC quiescence by inducing the dystroglycan receptor to associate with yes-associated protein, thereby sequestering it from the nucleus and preventing its prometastatic functions. These findings identify a brain-specific mechanism of DTC dormancy and highlight the need for a more thorough understanding of tumor dormancy to develop therapeutic approaches that prevent brain metastasis.

Spatial transcriptomics at subspot resolution with BayesSpace.

Recent spatial gene expression technologies enable comprehensive measurement of transcriptomic profiles while retaining spatial context. However, existing analysis methods do not address the limited resolution of the technology or use the spatial information efficiently. Here, we introduce BayesSpace, a fully Bayesian statistical method that uses the information from spatial neighborhoods for resolution enhancement of spatial transcriptomic data and for clustering analysis. We benchmark BayesSpace against current methods for spatial and non-spatial clustering and show that it improves identification of distinct intra-tissue transcriptional profiles from samples of the brain, melanoma, invasive ductal carcinoma and ovarian adenocarcinoma. Using immunohistochemistry and an in silico dataset constructed from scRNA-seq data, we show that BayesSpace resolves tissue structure that is not detectable at the original resolution and identifies transcriptional heterogeneity inaccessible to histological analysis. Our results illustrate BayesSpace's utility in facilitating the discovery of biological insights from spatial transcriptomic datasets.

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.

Towards precision prevention: Technologies for identifying healthy individuals with high risk of disease.

The rise of advanced technologies for characterizing human populations at the molecular level, from sequence to function, is shifting disease prevention paradigms toward personalized strategies. Because minimization of adverse outcomes is a key driver for treatment decisions for diseased populations, developing personalized therapy strategies represent an important dimension of both precision medicine and personalized prevention. In this commentary, we highlight recently developed enabling technologies in the field of DNA damage, DNA repair, and mutagenesis. We propose that omics approaches and functional assays can be integrated into population studies that fuse basic, translational and clinical research with commercial expertise in order to accelerate personalized prevention and treatment of cancer and other diseases linked to aberrant responses to DNA damage. This collaborative approach is generally applicable to efforts to develop data-driven, individualized prevention and treatment strategies for other diseases. We also recommend strategies for maximizing the use of biological samples for epidemiological studies, and for applying emerging technologies to clinical applications.

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.

Epithelial mesenchymal-like transition occurs in a subset of cells in castration resistant prostate cancer bone metastases.

TGFß is a known driver of epithelial-mesenchymal transition (EMT) which is associated with tumor aggressiveness and metastasis. However, EMT has not been fully explored in clinical specimens of castration-resistant prostate cancer (CRPC) metastases. To assess EMT in CRPC, gene expression analysis was performed on 149 visceral and bone metastases from 62 CRPC patients and immunohistochemical analysis was performed on 185 CRPC bone and visceral metastases from 42 CRPC patients. In addition, to assess the potential of metastases to seed further metastases the mitochondrial genome was sequenced at different metastatic sites in one patient. TGFß was increased in bone versus visceral metastases. While primarily cytoplasmic; nuclear and cytoplasmic Twist were significantly higher in bone than in visceral metastases. Slug and Zeb1 were unchanged, with the exception of nuclear Zeb1 being significantly higher in visceral metastases. Importantly, nuclear Twist, Slug, and Zeb1 were only present in a subset of epithelial cells that had an EMT-like phenotype. Underscoring the relevance of EMT-like cells, mitochondrial sequencing revealed that metastases could seed additional metastases in the same patient. In conclusion, while TGFß expression and EMT-associated protein expression is present in a considerable number of CRPC visceral and bone metastases, nuclear Twist, Slug, and Zeb1 localization and an EMT-like phenotype (elongated nuclei and cytoplasmic compartment) was only present in a small subset of CRPC bone metastases. Mitochondrial sequencing from different metastases in a CRPC patient provided evidence for the seeding of metastases from previously established metastases, highlighting the biological relevance of EMT-like behavior in CRPC metastases.

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.

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.

High-throughput sequencing of T-cell receptors reveals a homogeneous repertoire of tumour-infiltrating lymphocytes in ovarian cancer.

The cellular adaptive immune system mounts a response to many solid tumours mediated by tumour-infiltrating T lymphocytes (TILs). Basic measurements of these TILs, including total count, show promise as prognostic markers for a variety of cancers, including ovarian and colorectal. In addition, recent therapeutic advances are thought to exploit this immune response to effectively fight melanoma, with promising studies showing efficacy in additional cancers. However, many of the basic properties of TILs are poorly understood, including specificity, clonality, and spatial heterogeneity of the T-cell response. We utilize deep sequencing of rearranged T-cell receptor beta (TCRB) genes to characterize the basic properties of TILs in ovarian carcinoma. Due to somatic rearrangement during T-cell development, the TCR beta chain sequence serves as a molecular tag for each T-cell clone. Using these sequence tags, we assess similarities and differences between infiltrating T cells in discretely sampled sections of large tumours and compare to T cells from peripheral blood. Within the limits of sensitivity of our assay, the TIL repertoires show strong similarity throughout each tumour and are distinct from the circulating T-cell repertoire. We conclude that the cellular adaptive immune response within ovarian carcinomas is spatially homogeneous and distinct from the T-cell compartment of peripheral blood.

Simultaneous digital quantification and fluorescence-based size characterization of massively parallel sequencing libraries.

Due to the high cost of failed runs and suboptimal data yields, quantification and determination of fragment size range are crucial steps in the library preparation process for massively parallel sequencing (or next-generation sequencing). Current library quality control methods commonly involve quantification using real-time quantitative PCR and size determination using gel or capillary electrophoresis. These methods are laborious and subject to a number of significant limitations that can make library calibration unreliable. Herein, we propose and test an alternative method for quality control of sequencing libraries using droplet digital PCR (ddPCR). By exploiting a correlation we have discovered between droplet fluorescence and amplicon size, we achieve the joint quantification and size determination of target DNA with a single ddPCR assay. We demonstrate the accuracy and precision of applying this method to the preparation of sequencing libraries.

Heterogeneity of tumor-induced gene expression changes in the human metabolic network.

Reprogramming of cellular metabolism is an emerging hallmark of neoplastic transformation. However, it is not known how the expression of metabolic genes in tumors differs from that in normal tissues, or whether different tumor types exhibit similar metabolic changes. Here we compare expression patterns of metabolic genes across 22 diverse types of human tumors. Overall, the metabolic gene expression program in tumors is similar to that in the corresponding normal tissues. Although expression changes of some metabolic pathways (e.g., upregulation of nucleotide biosynthesis and glycolysis) are frequently observed across tumors, expression changes of other pathways (e.g., oxidative phosphorylation) are very heterogeneous. Our analysis also suggests that the expression changes of some metabolic genes (e.g., isocitrate dehydrogenase and fumarate hydratase) may enhance or mimic the effects of recurrent mutations in tumors. On the level of individual biochemical reactions, many hundreds of metabolic isoenzymes show significant and tumor-specific expression changes. These isoenzymes are potential targets for anticancer therapy.

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.

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.

Quantification of random mutations in the mitochondrial genome.

Mitochondrial DNA (mtDNA) mutations contribute to the pathology of a number of age-related disorders, including Parkinson disease [A. Bender et al., Nat. Genet. 38 (2006) 515,Y. Kraytsberg et al., Nat. Genet. 38 (2006) 518], muscle-wasting [J. Wanagat, Z. Cao, P. Pathare, J.M. Aiken, FASEB J. 15 (2001) 322], and the metastatic potential of cancers [K. Ishikawa et al., Science 320 (2008) 661]. The impact of mitochondrial DNA mutations on a wide variety of human diseases has made it increasingly important to understand the mechanisms that drive mitochondrial mutagenesis. In order to provide new insight into the etiology and natural history of mtDNA mutations, we have developed an assay that can detect mitochondrial mutations in a variety of tissues and experimental settings [M. Vermulst et al., Nat. Genet. 40 (2008) 4, M. Vermulst et al., Nat. Genet. 39 (2007) 540]. This methodology, termed the Random Mutation Capture assay, relies on single-molecule amplification to detect rare mutations among millions of wild-type bases [J.H. Bielas, L.A. Loeb, Nat. Methods 2 (2005) 285], and can be used to analyze mitochondrial mutagenesis to a single base pair level in mammals.