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.

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.

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.

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.

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.

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.

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.

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.

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.

Long term rapamycin treatment improves mitochondrial DNA quality in aging mice.

Age-induced mitochondrial DNA deletion mutations may underlie cell loss and tissue aging. Rapamycin extends mouse lifespan and modulates mitochondrial quality control. We hypothesized that reduced deletion mutation abundance may contribute to rapamycin's life extension effects. To test this hypothesis, genetically heterogeneous male and female mice were treated with rapamycin, compounded in chow at 14 or 42 ppm, from 9 months to 22 months of age. Mice under a 40% dietary restriction were included as a control known to protect mtDNA quality. To determine if chronic rapamycin treatment affects mitochondrial DNA quality, we assayed mtDNA deletion frequency and electron transport chain deficient fiber abundances in mouse quadriceps muscle. At 42 ppm rapamycin, we observed a 57% decrease in deletion frequency, a 2.8-fold decrease in ETC deficient fibers, and a 3.4-fold increase in the number of mice without electron transport chain deficient fibers. We observed a similar trend with the 14 ppm dose. DR significantly decreased ETC deficient fiber abundances with a trend toward lower mtDNA deletion frequency. The effects of rapamycin treatment on mitochondrial DNA quality were greatest in females at the highest dose. Rapamycin treatment at 14 ppm did not affect muscle mass or function. Dietary restriction also reduced deletion frequency and ETC deficient fibers. These data support the concept that the lifespan extending effects of rapamycin treatment result from enhanced mitochondrial DNA quality.

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.

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.

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.

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.