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.

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.

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.

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.

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.

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.

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.