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.

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.

Rapidly inducible Cas9 and DSB-ddPCR to probe editing kinetics.

We developed a chemically inducible Cas9 (ciCas9) and a droplet digital PCR assay for double-strand breaks (DSB-ddPCR) to investigate the kinetics of Cas9-mediated generation and repair of DSBs in cells. ciCas9 is a rapidly activated, single-component Cas9 variant engineered by replacing the protein's REC2 domain with the BCL-xL protein and fusing an interacting BH3 peptide to the C terminus. ciCas9 can be tunably activated by a compound that disrupts the BCL-xL-BH3 interaction within minutes. DSB-ddPCR demonstrates time-resolved, highly quantitative, and targeted measurement of DSBs. Combining these tools facilitated an unprecedented exploration of the kinetics of Cas9-mediated DNA cleavage and repair. We find that sgRNAs targeting different sites generally induce cleavage within minutes and repair within 1 or 2 h. However, we observe distinct kinetic profiles, even for proximal sites, and this suggests that target sequence and chromatin state modulate cleavage and repair kinetics.

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.

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.

Mitochondrial point mutations do not limit the natural lifespan of mice.

Whether mitochondrial mutations cause mammalian aging, or are merely correlated with it, is an area of intense debate. Here, we use a new, highly sensitive assay to redefine the relationship between mitochondrial mutations and age. We measured the in vivo rate of change of the mitochondrial genome at a single-base pair level in mice, and we demonstrate that the mutation frequency in mouse mitochondria is more than ten times lower than previously reported. Although we observed an 11-fold increase in mitochondrial point mutations with age, we report that a mitochondrial mutator mouse was able to sustain a 500-fold higher mutation burden than normal mice, without any obvious features of rapidly accelerated aging. Thus, our results strongly indicate that mitochondrial mutations do not limit the lifespan of wild-type mice.

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.

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.

Tumor Regression Following Engineered Polyomavirus-Specific T Cell Therapy in Immune Checkpoint Inhibitor-Refractory Merkel Cell Carcinoma.

Although immune check-point inhibitors (CPIs) revolutionized treatment of Merkel cell carcinoma (MCC), patients with CPI-refractory MCC lack effective therapy. More than 80% of MCC express T-antigens encoded by Merkel cell polyomavirus, which is an ideal target for T-cell receptor (TCR)-based immunotherapy. However, MCC often repress HLA expression, requiring additional strategies to reverse the downregulation for allowing T cells to recognize their targets. We identified TCRMCC1 that recognizes a T-antigen epitope restricted to human leukocyte antigen (HLA)-A*02:01. Seven CPI-refractory metastatic MCC patients received CD4 and CD8 T cells transduced with TCRMCC1 (TTCR-MCC1) preceded either by lymphodepleting chemotherapy or an HLA-upregulating regimen (single-fraction radiation therapy (SFRT) or systemic interferon gamma (IFNγ)) with concurrent avelumab. Two patients who received preceding SFRT and IFNγ respectively experienced tumor regression. One experienced regression of 13/14 subcutaneous lesions with 1 'escape' lesion and the other had delayed tumor regression in all lesions after initial progression. Although TTCR-MCC1 cells with an activated phenotype infiltrated tumors including the 'escape' lesion, all progressing lesions transcriptionally lacked HLA expression. While SFRT/IFNγ did not immediately upregulate tumor HLA expression, a secondary endogenous antigen-specific T cell infiltrate was detected in one of the regressing tumors and associated with HLA upregulation, indicating in situ immune responses have the potential to reverse HLA downregulation. Indeed, supplying a strong co-stimulatory signal via a CD200R-CD28 switch receptor allows TTCR-MCC1 cells to control HLA-downregulated MCC cells in a xenograft mouse model, upregulating HLA expression. Our results demonstrate the potential of TCR gene therapy for metastatic MCC and propose a next strategy for overcoming epigenetic downregulation of HLA in MCC.

A random mutation capture assay to detect genomic point mutations in mouse tissue.

Herein, a detailed protocol for a random mutation capture (RMC) assay to measure nuclear point mutation frequency in mouse tissue is described. This protocol is a simplified version of the original method developed for human tissue that is easier to perform, yet retains a high sensitivity of detection. In contrast to assays relying on phenotypic selection of reporter genes in transgenic mice, the RMC assay allows direct detection of mutations in endogenous genes in any mouse strain. Measuring mutation frequency within an intron of a transcribed gene, we show this assay to be highly reproducible. We analyzed mutation frequencies from the liver tissue of animals with a mutation within the intrinsic exonuclease domains of the two major DNA polymerases, δ and ε. These mice exhibited significantly higher mutation frequencies than did wild-type animals. A comparison with a previous analysis of these genotypes in Big Blue mice revealed the RMC assay to be more sensitive than the Big Blue assay for this application. As RMC does not require analysis of a particular gene, simultaneous analysis of mutation frequency at multiple genetic loci is feasible. This assay provides a versatile alternative to transgenic mouse models for the study of mutagenesis in vivo.

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.

Aging neural progenitor cells have decreased mitochondrial content and lower oxidative metabolism.

Although neurogenesis occurs in discrete areas of the adult mammalian brain, neural progenitor cells (NPCs) produce fewer new neurons with age. To characterize the molecular changes that occur during aging, we performed a proteomic comparison between primary-cultured NPCs from the young adult and aged mouse forebrain. This analysis yielded changes in proteins necessary for cellular metabolism. Mitochondrial quantity and oxygen consumption rates decrease with aging, although mitochondrial DNA in aged NPCs does not have increased mutation rates. In addition, aged cells are resistant to the mitochondrial inhibitor rotenone and proliferate in response to lowered oxygen conditions. These results demonstrate that aging NPCs display an altered metabolic phenotype, characterized by a coordinated shift in protein expression, subcellular structure, and metabolic physiology.

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.

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.

Targeting an alternate Wilms' tumor antigen 1 peptide bypasses immunoproteasome dependency.

Designing effective antileukemic immunotherapy will require understanding mechanisms underlying tumor control or resistance. Here, we report a mechanism of escape from immunologic targeting in an acute myeloid leukemia (AML) patient, who relapsed 1 year after immunotherapy with engineered T cells expressing a human leukocyte antigen A*02 (HLA-A2)-restricted T cell receptor (TCR) specific for a Wilms' tumor antigen 1 epitope, WT1126-134 (TTCR-C4). Resistance occurred despite persistence of functional therapeutic T cells and continuous expression of WT1 and HLA-A2 by the patient's AML cells. Analysis of the recurrent AML revealed expression of the standard proteasome, but limited expression of the immunoproteasome, specifically the beta subunit 1i (ß1i), which is required for presentation of WT1126-134. An analysis of a second patient treated with TTCR-C4 demonstrated specific loss of AML cells coexpressing ß1i and WT1. To determine whether the WT1 protein continued to be processed and presented in the absence of immunoproteasome processing, we identified and tested a TCR targeting an alternative, HLA-A2-restricted WT137-45 epitope that was generated by immunoproteasome-deficient cells, including WT1-expressing solid tumor lines. T cells expressing this TCR (TTCR37-45) killed the first patients' relapsed AML resistant to WT1126-134 targeting, as well as other primary AML, in vitro. TTCR37-45 controlled solid tumor lines lacking immunoproteasome subunits both in vitro and in an NSG mouse model. As proteasome composition can vary in AML, defining and preferentially targeting these proteasome-independent epitopes may maximize therapeutic efficacy and potentially circumvent AML immune evasion by proteasome-related immunoediting.

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.

Distinct populations of antigen-specific tissue-resident CD8+ T cells in human cervix mucosa.

The ectocervix is part of the lower female reproductive tract (FRT), which is susceptible to sexually transmitted infections (STIs). Comprehensive knowledge of the phenotypes and T cell receptor (TCR) repertoire of tissue-resident memory T cells (TRMs) in the human FRT is lacking. We took single-cell RNA-Seq approaches to simultaneously define gene expression and TCR clonotypes of the human ectocervix. There were significantly more CD8+ than CD4+ T cells. Unsupervised clustering and trajectory analysis identified distinct populations of CD8+ T cells with IFNGhiGZMBloCD69hiCD103lo or IFNGloGZMBhiCD69medCD103hi phenotypes. Little overlap was seen between their TCR repertoires. Immunofluorescence staining showed that CD103+CD8+ TRMs were preferentially localized in the epithelium, whereas CD69+CD8+ TRMs were distributed evenly in the epithelium and stroma. Ex vivo assays indicated that up to 14% of cervical CD8+ TRM clonotypes were HSV-2 reactive in HSV-2-seropositive persons, reflecting physiologically relevant localization. Our studies identified subgroups of CD8+ TRMs in the human ectocervix that exhibited distinct expression of antiviral defense and tissue residency markers, anatomic locations, and TCR repertoires that target anatomically relevant viral antigens. Optimization of the location, number, and function of FRT TRMs is an important approach for improving host defenses to STIs.

LOH-proficient embryonic stem cells: a model of cancer progenitor cells?

Cancers are thought to originate in stem cells through the accumulation of multiple mutations. Some of these mutations result in a loss of heterozygosity (LOH). A recent report demonstrates that exposure of mouse embryonic stem cells to nontoxic amounts of mutagens triggers a marked increase in the frequency of LOH. Thus, mutagen induction of LOH in embryonic stem cells suggests a new pathway to account for the multiple homozygous mutations in human tumors. This induction could mimic early mutagenic events that generate cancers in human tissue stem cells.

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.

High information content assays for genetic toxicology testing: A report of the International Workshops on Genotoxicity Testing (IWGT).

We live in an era of 'big data', where the volume, velocity, and variety of the data being generated is increasingly influencing the way toxicological sciences are practiced. With this in mind, a workgroup was formed for the 2017 International Workshops on Genotoxicity Testing (IWGT) to consider the use of high information content data in genetic toxicology assessments. Presentations were given on adductomics, global transcriptional profiling, error-reduced single-molecule sequencing, and cellular phenotype-based assays, which were identified as methodologies that are relevant to present-day genetic toxicology assessments. Presenters and workgroup members discussed the state of the science for these methodologies, their potential use in genetic toxicology, current limitations, and the future work necessary to advance their utility and application. The session culminated with audience-assisted SWOT (strength, weakness, opportunities, and threats) analyses. The summary report described herein is structured similarly. A major conclusion of the workgroup is that while conventional regulatory genetic toxicology testing has served the public well over the last several decades, it does not provide the throughput that has become necessary in modern times, and it does not generate the mechanistic information that risk assessments ideally take into consideration. The high information content assay platforms that were discussed in this session, as well as others under development, have the potential to address aspect(s) of these issues and to meet new expectations in the field of genetic toxicology.