Single-cell analysis of immune and stroma cell remodeling in clear cell renal cell carcinoma primary tumors and bone metastatic lesions.

BACKGROUND: Despite therapeutic advances, once a cancer has metastasized to the bone, it represents a highly morbid and lethal disease. One third of patients with advanced clear cell renal cell carcinoma (ccRCC) present with bone metastasis at the time of diagnosis. However, the bone metastatic niche in humans, including the immune and stromal microenvironments, has not been well-defined, hindering progress towards identification of therapeutic targets. METHODS: We collected fresh patient samples and performed single-cell transcriptomic profiling of solid metastatic tissue (Bone Met), liquid bone marrow at the vertebral level of spinal cord compression (Involved), and liquid bone marrow from a different vertebral body distant from the tumor site but within the surgical field (Distal), as well as bone marrow from patients undergoing hip replacement surgery (Benign). In addition, we incorporated single-cell data from primary ccRCC tumors (ccRCC Primary) for comparative analysis. RESULTS: The bone marrow of metastatic patients is immune-suppressive, featuring increased, exhausted CD8 + cytotoxic T cells, T regulatory cells, and tumor-associated macrophages (TAM) with distinct transcriptional states in metastatic lesions. Bone marrow stroma from tumor samples demonstrated a tumor-associated mesenchymal stromal cell population (TA-MSC) that appears to be supportive of epithelial-to mesenchymal transition (EMT), bone remodeling, and a cancer-associated fibroblast (CAFs) phenotype. This stromal subset is associated with poor progression-free and overall survival and also markedly upregulates bone remodeling through the dysregulation of RANK/RANKL/OPG signaling activity in bone cells, ultimately leading to bone resorption. CONCLUSIONS: These results provide a comprehensive analysis of the bone marrow niche in the setting of human metastatic cancer and highlight potential therapeutic targets for both cell populations and communication channels.

Dissecting the immune suppressive human prostate tumor microenvironment via integrated single-cell and spatial transcriptomic analyses.

The treatment of low-risk primary prostate cancer entails active surveillance only, while high-risk disease requires multimodal treatment including surgery, radiation therapy, and hormonal therapy. Recurrence and development of metastatic disease remains a clinical problem, without a clear understanding of what drives immune escape and tumor progression. Here, we comprehensively describe the tumor microenvironment of localized prostate cancer in comparison with adjacent normal samples and healthy controls. Single-cell RNA sequencing and high-resolution spatial transcriptomic analyses reveal tumor context dependent changes in gene expression. Our data indicate that an immune suppressive tumor microenvironment associates with suppressive myeloid populations and exhausted T-cells, in addition to high stromal angiogenic activity. We infer cell-to-cell relationships from high throughput ligand-receptor interaction measurements within undissociated tissue sections. Our work thus provides a highly detailed and comprehensive resource of the prostate tumor microenvironment as well as tumor-stromal cell interactions.

A transcriptional metastatic signature predicts survival in clear cell renal cell carcinoma.

Clear cell renal cell carcinoma (ccRCC) is the most common type of kidney cancer in adults. When ccRCC is localized to the kidney, surgical resection or ablation of the tumor is often curative. However, in the metastatic setting, ccRCC remains a highly lethal disease. Here we use fresh patient samples that include treatment-naive primary tumor tissue, matched adjacent normal kidney tissue, as well as tumor samples collected from patients with bone metastases. Single-cell transcriptomic analysis of tumor cells from the primary tumors reveals a distinct transcriptional signature that is predictive of metastatic potential and patient survival. Analysis of supporting stromal cells within the tumor environment demonstrates vascular remodeling within the endothelial cells. An in silico cell-to-cell interaction analysis highlights the CXCL9/CXCL10-CXCR3 axis and the CD70-CD27 axis as potential therapeutic targets. Our findings provide biological insights into the interplay between tumor cells and the ccRCC microenvironment.

Long-term Oncologic Impact of Positive Anterior and Posterior Surgical Margins After Radical Prostatectomy.

OBJECTIVE: The objective of this study was to evaluate the impact of the anterior/posterior status of positive surgical margin (PSM) on long-term outcomes after radical prostatectomy for prostate cancer. PATIENTS AND METHODS: We included 391 consecutive PSM patients after radical prostatectomy between 1993 and 2007 excluding cases with multiple location PSM or lack of anterior/posterior status data. The oncologic impact of anterior-PSM and posterior-PSM were examined by Kaplan-Meier analysis and the Cox proportional hazards model. RESULTS: There were 115 cases (29.4%) with apex-PSM, 257 cases (65.7%) with peripheral PSM, and 19 cases (4.9%) with bladder neck PSM. Among the 257 peripheral PSM cases, 58 cases (22.6%) were with anterior-PSM, 174 cases (67.7%) were with posterior-PSM, and 25 cases (9.7%) were with both anterior and posterior PSM. Over a median follow-up of 12.6 years, patients with anterior-PSM, especially those with low to intermediate Gleason score (≤7), showed a biochemical recurrence (BCR) prognosis similar to those with apex-PSM. In contrast, patients with posterior-PSM showed significantly higher BCR risk on both univariate and multivariate analyses when compared with those with apex-PSM. No impact on metastasis-free survival or overall survival was observed. CONCLUSIONS: In our study, we found that prostate cancer patients with anterior-PSM showed a more favorable BCR prognosis similar to those with apex-PSM when comparing to patients with posterior-PSM. Our study results may help physicians to choose different treatment options for patients diagnosed with different PSM status including considering further adjuvant treatment for patients with posterior-PSM.

MicroRNAs Cause Accelerated Decay of Short-Tailed Target mRNAs.

MicroRNAs (miRNAs) specify the recruitment of deadenylases to mRNA targets. Despite this recruitment, we find that miRNAs have almost no effect on steady-state poly(A)-tail lengths of their targets in mouse fibroblasts, which motivates the acquisition of pre-steady-state measurements of the effects of miRNAs on tail lengths, mRNA levels, and translational efficiencies. Effects on translational efficiency are minimal compared to effects on mRNA levels, even for newly transcribed target mRNAs. Effects on target mRNA levels accumulate as the mRNA population approaches steady state, whereas effects on tail lengths peak for recently transcribed target mRNAs and then subside. Computational modeling of this phenomenon reveals that miRNAs cause not only accelerated deadenylation of their targets but also accelerated decay of short-tailed target molecules. This unanticipated effect of miRNAs largely prevents short-tailed target mRNAs from accumulating despite accelerated target deadenylation. The net result is a nearly imperceptible change to the steady-state tail-length distribution of targeted mRNAs.

The Dynamics of Cytoplasmic mRNA Metabolism.

For all but a few mRNAs, the dynamics of metabolism are unknown. Here, we developed an experimental and analytical framework for examining these dynamics for mRNAs from thousands of genes. mRNAs of mouse fibroblasts exit the nucleus with diverse intragenic and intergenic poly(A)-tail lengths. Once in the cytoplasm, they have a broad (1000-fold) range of deadenylation rate constants, which correspond to cytoplasmic lifetimes. Indeed, with few exceptions, degradation appears to occur primarily through deadenylation-linked mechanisms, with little contribution from either endonucleolytic cleavage or deadenylation-independent decapping. Most mRNA molecules degrade only after their tail lengths fall below 25 nt. Decay rate constants of short-tailed mRNAs vary broadly (1000-fold) and are larger for short-tailed mRNAs that have previously undergone more rapid deadenylation. This coupling helps clear rapidly deadenylated mRNAs, enabling the large range in deadenylation rate constants to impart a similarly large range in stabilities.

Assessment of 5-year overall survival in bladder cancer patients with incidental prostate cancer identified at radical cystoprostatectomy.

OBJECTIVE: To investigate the oncological impact of incidental prostate cancer (iPCa) found during radical cystoprostatectomy (RCP) on overall survival (OS) prognosis of urothelial carcinoma of the bladder (BCa). PATIENTS AND METHODS: A total of 122 RCP cases resected between 2002 and 2012 at our center were included for study. Survival of BCa patient was compared using the Kaplan-Meier method and the log-rank test. Cox proportional hazards regression models were used to analyze the impact of iPCa on the 5-year overall mortality of BCa patients after RCP. RESULTS: Among the 122 BCa cases that underwent RCP, 38 cases (31.1%) had iPCa, in which, 17 cases (44.7%) were identified as clinically significant iPCa (csPCa). BCa patients with iPCa were older (71 vs 64 years, p = 0.004) and had higher preoperative PSA level (3.1 ng/mL vs 1.4 ng/mL, p = 0.017) when compared to those without iPCa. Cases with iPCa showed a more favorable 5-year OS than cases without iPCa, although this difference did not reach statistical significance (p = 0.219). When excluding the higher risk cases with Gleason score (GS) ≥ 4 + 3 and/or preoperative PSA > 10 ng/mL, BCa patients with iPCa showed a significantly longer OS than cases without iPCa on univariate analysis (p = 0.044), but not on multivariate analysis (p = 0.125). CONCLUSION: Our results demonstrated that the frequent findings of low-risk iPCa in BCa patients could indicate the potential possibility of shared pathogenesis pathways between iPCa and BCa. Future study with a larger cohort is warranted to validate this result.

Impact of Multifocality and Multilocation of Positive Surgical Margin After Radical Prostatectomy on Predicting Oncological Outcome.

OBJECTIVE: To assess the impact of focality and location of positive surgical margins (PSM) on long-term outcomes after radical prostatectomy (RP) for prostate cancer (PCa), including biochemical recurrence (BCR), metastasis and overall mortality. PATIENTS AND METHODS: From a total of 2796 cases of RP between 1993 and 2007 in our single hospital, 476 cases with PSMs were identified and included in this study. PSM location was categorized into apex, peripheral, and bladder neck. Survival was estimated using the Kaplan-Meier method. Cox proportional hazard regression models were used to analyze the impact of PSM focality and location status on oncologic survival. RESULTS: Of these 476 cases with PSMs, 335 (70.4%) cases were with single focal (sF) PSMs and 141 (29.6%) cases were with multifocal (mF) PSMs. Furthermore, 406 (85.3%) cases were found to have single location (sL) PSMs, and 70 (14.7%) cases were with multilocation (mL) PSMs. The median follow-up was 12.9 years. mF-PSMs and mL-PSMs showed significant impact on increased BCR risk on univariate analysis, and mL-PSMs remained significant on multivariate analysis (P = .048). Furthermore, the combination of multifocality and multilocation showed added prognostic value on predicting BCR-free survival, but not on metastasis-free survival or overall survival. CONCLUSION: The presence of mF-PSMs and mL-PSMs, and especially the combination of both, demonstrated significant impact on BCR prognosis. Patients with apex sLsF-PSMs were less likely to have BCR when compared with all those with non-apex sLsF-PSMs. These results should be considered when evaluating patients for adjuvant therapy.

Predicting microRNA targeting efficacy in Drosophila.

BACKGROUND: MicroRNAs (miRNAs) are short regulatory RNAs that derive from hairpin precursors. Important for understanding the functional roles of miRNAs is the ability to predict the messenger RNA (mRNA) targets most responsive to each miRNA. Progress towards developing quantitative models of miRNA targeting in Drosophila and other invertebrate species has lagged behind that of mammals due to the paucity of datasets measuring the effects of miRNAs on mRNA levels. RESULTS: We acquired datasets suitable for the quantitative study of miRNA targeting in Drosophila. Analyses of these data expanded the types of regulatory sites known to be effective in flies, expanded the mRNA regions with detectable targeting to include 5' untranslated regions, and identified features of site context that correlate with targeting efficacy in fly cells. Updated evolutionary analyses evaluated the probability of conserved targeting for each predicted site and indicated that more than a third of the Drosophila genes are preferentially conserved targets of miRNAs. Based on these results, a quantitative model was developed to predict targeting efficacy in insects. This model performed better than existing models, and it drives the most recent version, v7, of TargetScanFly. CONCLUSIONS: Our evolutionary and functional analyses expand the known scope of miRNA targeting in flies and other insects. The existence of a quantitative model that has been developed and trained using Drosophila data will provide a valuable resource for placing miRNAs into gene regulatory networks of this important experimental organism.

The influence of microRNAs and poly(A) tail length on endogenous mRNA-protein complexes.

BACKGROUND: All mRNAs are bound in vivo by proteins to form mRNA-protein complexes (mRNPs), but changes in the composition of mRNPs during posttranscriptional regulation remain largely unexplored. Here, we have analyzed, on a transcriptome-wide scale, how microRNA-mediated repression modulates the associations of the core mRNP components eIF4E, eIF4G, and PABP and of the decay factor DDX6 in human cells. RESULTS: Despite the transient nature of repressed intermediates, we detect significant changes in mRNP composition, marked by dissociation of eIF4G and PABP, and by recruitment of DDX6. Furthermore, although poly(A)-tail length has been considered critical in post-transcriptional regulation, differences in steady-state tail length explain little of the variation in either PABP association or mRNP organization more generally. Instead, relative occupancy of core components correlates best with gene expression. CONCLUSIONS: These results indicate that posttranscriptional regulatory factors, such as microRNAs, influence the associations of PABP and other core factors, and do so without substantially affecting steady-state tail length.

mRNA poly(A)-tail changes specified by deadenylation broadly reshape translation in Drosophila oocytes and early embryos.

Because maturing oocytes and early embryos lack appreciable transcription, posttranscriptional regulatory processes control their development. To better understand this control, we profiled translational efficiencies and poly(A)-tail lengths throughout Drosophila oocyte maturation and early embryonic development. The correspondence between translational-efficiency changes and tail-length changes indicated that tail-length changes broadly regulate translation until gastrulation, when this coupling disappears. During egg activation, relative changes in poly(A)-tail length, and thus translational efficiency, were largely retained in the absence of cytoplasmic polyadenylation, which indicated that selective poly(A)-tail shortening primarily specifies these changes. Many translational changes depended on PAN GU and Smaug, and these changes were largely attributable to tail-length changes. Our results also revealed the presence of tail-length-independent mechanisms that maintained translation despite tail-length shortening during oocyte maturation, and prevented essentially all translation of bicoid and several other mRNAs before egg activation. In addition to these fundamental insights, our results provide valuable resources for future studies.

Poly(A)-tail profiling reveals an embryonic switch in translational control.

Poly(A) tails enhance the stability and translation of most eukaryotic messenger RNAs, but difficulties in globally measuring poly(A)-tail lengths have impeded greater understanding of poly(A)-tail function. Here we describe poly(A)-tail length profiling by sequencing (PAL-seq) and apply it to measure tail lengths of millions of individual RNAs isolated from yeasts, cell lines, Arabidopsis thaliana leaves, mouse liver, and zebrafish and frog embryos. Poly(A)-tail lengths were conserved between orthologous mRNAs, with mRNAs encoding ribosomal proteins and other 'housekeeping' proteins tending to have shorter tails. As expected, tail lengths were coupled to translational efficiencies in early zebrafish and frog embryos. However, this strong coupling diminished at gastrulation and was absent in non-embryonic samples, indicating a rapid developmental switch in the nature of translational control. This switch complements an earlier switch to zygotic transcriptional control and explains why the predominant effect of microRNA-mediated deadenylation concurrently shifts from translational repression to mRNA destabilization.

Extensive alternative polyadenylation during zebrafish development.

The post-transcriptional fate of messenger RNAs (mRNAs) is largely dictated by their 3' untranslated regions (3' UTRs), which are defined by cleavage and polyadenylation (CPA) of pre-mRNAs. We used poly(A)-position profiling by sequencing (3P-seq) to map poly(A) sites at eight developmental stages and tissues in the zebrafish. Analysis of over 60 million 3P-seq reads substantially increased and improved existing 3' UTR annotations, resulting in confidently identified 3' UTRs for >79% of the annotated protein-coding genes in zebrafish. mRNAs from most zebrafish genes undergo alternative CPA, with those from more than a thousand genes using different dominant 3' UTRs at different stages. These included one of the poly(A) polymerase genes, for which alternative CPA reinforces its repression in the ovary. 3' UTRs tend to be shortest in the ovaries and longest in the brain. Isoforms with some of the shortest 3' UTRs are highly expressed in the ovary, yet absent in the maternally contributed RNAs of the embryo, perhaps because their 3' UTRs are too short to accommodate a uridine-rich motif required for stability of the maternal mRNA. At 2 h post-fertilization, thousands of unique poly(A) sites appear at locations lacking a typical polyadenylation signal, which suggests a wave of widespread cytoplasmic polyadenylation of mRNA degradation intermediates. Our insights into the identities, formation, and evolution of zebrafish 3' UTRs provide a resource for studying gene regulation during vertebrate development.

Identification of small molecule inhibitors of pyruvate kinase M2.

A common feature of tumors arising from diverse tissue types is a reliance on aerobic glycolysis for glucose metabolism. This metabolic difference between cancer cells and normal cells could be exploited for therapeutic benefit in patients. Cancer cells universally express the M2 isoform of the glycolytic enzyme pyruvate kinase (PKM2), and previous work has demonstrated that PKM2 expression is necessary for aerobic glycolysis and cell proliferation in vivo. Because most normal tissues express an isoform of pyruvate kinase other than PKM2, selective targeting of PKM2 provides an opportunity to target cell metabolism for cancer therapy. PKM2 has an identical catalytic site as the related M1 splice variant (PKM1). However, isoform selective inhibition is possible as PKM2 contains a unique region for allosteric regulation. We have screened a library of greater than 1,00,000 small molecules to identify such inhibitors. The inhibitors identified for PKM2 fell primarily into three distinct structural classes. The most potent PKM2 inhibitor resulted in decreased glycolysis and increased cell death following loss of growth factor signaling. At least part of this effect was due to on-target PKM2 inhibition as less cell death was observed in cells engineered to express PKM1. These data suggest that isoform selective inhibition of PKM2 with small molecules is feasible and support the hypothesis that inhibition of glucose metabolism in cancer cells is a viable strategy to treat human malignancy.

Ribosomal synthesis of N-methyl peptides.

N-methyl amino acids (N-Me AAs) are a common component of nonribosomal peptides (NRPs), a class of natural products from which many clinically important therapeutics are obtained. N-Me AAs confer peptides with increased conformational rigidity, membrane permeability, and protease resistance. Hence, these analogues are highly desirable building blocks in the ribosomal synthesis of unnatural peptide libraries, from which functional, NRP-like molecules may be identified. By supplementing a reconstituted Escherichia coli translation system with specifically aminoacylated total tRNA that has been chemically methylated, we have identified three N-Me AAs (N-Me Leu, N-Me Thr, and N-Me Val) that are efficiently incorporated into peptides by the ribosome. Moreover, we have demonstrated the synthesis of peptides containing up to three N-Me AAs, a number comparable to that found in many NRP drugs. With improved incorporation efficiency and translational fidelity, it may be possible to synthesize combinatorial libraries of peptides that contain multiple N-Me AAs. Such libraries could be subjected to in vitro selection methods to identify drug-like, high-affinity ligands for protein targets of interest.