Lineage-specific canonical and non-canonical activity of EZH2 in advanced prostate cancer subtypes

Enhancer of zeste homolog 2 (EZH2) is a histone methyltransferase and emerging therapeutic target that is overexpressed in most castration-resistant prostate cancers and implicated as a driver of disease progression and resistance to hormonal therapies. Here we define the lineage-specific action and differential activity of EZH2 in both prostate adenocarcinoma (PRAD) and neuroendocrine prostate cancer (NEPC) subtypes of advanced prostate cancer to better understand the role of EZH2 in modulating differentiation, lineage plasticity, and to identify mediators of response and resistance to EZH2 inhibitor therapy. Mechanistically, EZH2 modulates bivalent genes that results in upregulation of NEPC-associated transcriptional drivers (e.g., ASCL1) and neuronal gene programs, and leads to forward differentiation after targeting EZH2 in NEPC. Subtype-specific downstream effects of EZH2 inhibition on cell cycle genes support the potential rationale for co-targeting cyclin/CDK to overcome resistance to EZH2 inhibition.

Heterogeneity and transcriptional drivers of triple-negative breast cancer

Triple-negative breast cancer (TNBC) is a heterogeneous disease with limited treatment options. To characterize TNBC heterogeneity, we defined transcriptional, epigenetic, and metabolic subtypes and subtype-driving super-enhancers and transcription factors by combining functional and molecular profiling with computational analyses. Single-cell RNA sequencing revealed relative homogeneity of the major transcriptional subtypes (luminal, basal, and mesenchymal) within samples. We found that mesenchymal TNBCs share features with mesenchymal neuroblastoma and rhabdoid tumors and that the PRRX1 transcription factor is a key driver of these tumors. PRRX1 is sufficient for inducing mesenchymal features in basal but not in luminal TNBC cells via reprogramming super-enhancer landscapes, but it is not required for mesenchymal state maintenance or for cellular viability. Our comprehensive, large-scale, multiplatform, multiomics study of both experimental and clinical TNBC is an important resource for the scientific and clinical research communities and opens venues for future investigation.

Androgen receptor inhibition suppresses anti-tumor neutrophil response against bone metastatic prostate cancer via regulation of TβRI expression

Bone metastatic disease of prostate cancer (PCa) is incurable and progression in bone is largely dictated by tumor-stromal interactions in the bone microenvironment. We showed previously that bone neutrophils initially inhibit bone metastatic PCa growth yet metastatic PCa becomes resistant to neutrophil response. Further, neutrophils isolated from tumor-bone lost their ability to suppress tumor growth through unknown mechanisms. With this study, our goal was to define the impact of metastatic PCa on neutrophil function throughout tumor progression and to determine the potential of neutrophils as predictive biomarkers of metastatic disease. Using patient peripheral blood polymorphonuclear neutrophils (PMNs), we identified that PCa progression dictates PMN cell surface markers and gene expression, but not cytotoxicity against PCa. Importantly, we also identified a novel phenomenon in which second generation androgen deprivation therapy (ADT) suppresses PMN cytotoxicity via increased transforming growth factor beta receptor I (TβRI). High dose testosterone and genetic or pharmacologic TβRI inhibition rescued androgen receptor-mediated neutrophil suppression and restored neutrophil anti-tumor immune response. These studies highlight the ability to leverage standard-care ADT to generate neutrophil anti-tumor responses against bone metastatic PCa.

Telomerase Upregulation Induces Progression of Mouse BrafV600E-Driven Thyroid Cancers and Triggers Nontelomeric Effects

Mutations in the promoter of the telomerase reverse transcriptase (TERT) gene are the paradigm of a cross-cancer alteration in a noncoding region. TERT promoter mutations (TPM) are biomarkers of poor prognosis in cancer, including thyroid tumors. TPMs enhance TERT transcription, which is otherwise silenced in adult tissues, thus reactivating a bona fide oncoprotein. To study TERT deregulation and its downstream consequences, we generated a Tert mutant promoter mouse model via CRISPR/Cas9 engineering of the murine equivalent locus (Tert-123C>T) and crossed it with thyroid-specific BrafV600E-mutant mice. We also employed an alternative model of Tert overexpression (K5-Tert). Whereas all BrafV600E animals developed well-differentiated papillary thyroid tumors, 29% and 36% of BrafV600E+Tert-123C>T and BrafV600E+K5-Tert mice progressed to poorly differentiated cancers at week 20, respectively. Tert-upregulated tumors showed increased mitosis and necrosis in areas of solid growth, and older animals displayed anaplastic-like features, that is, spindle cells and macrophage infiltration. Murine TPM increased Tert transcription in vitro and in vivo, but temporal and intratumoral heterogeneity was observed. RNA-sequencing of thyroid tumor cells showed that processes other than the canonical Tert-mediated telomere maintenance role operate in these specimens. Pathway analysis showed that MAPK and PI3K/AKT signaling, as well as processes not previously associated with this tumor etiology, involving cytokine, and chemokine signaling, were overactivated. These models constitute useful preclinical tools to understand the cell-autonomous and microenvironment-related consequences of Tert-mediated progression in advanced thyroid cancers and other aggressive tumors carrying TPMs.

IMPLICATIONS

Telomerase-driven cancer progression activates pathways that can be dissected and perhaps therapeutically exploited.

Proteogenomic analysis of chemo-refractory high-grade serous ovarian cancer

To improve the understanding of chemo-refractory high-grade serous ovarian cancers (HGSOCs), we characterized the proteogenomic landscape of 242 (refractory and sensitive) HGSOCs, representing one discovery and two validation cohorts across two biospecimen types (formalin-fixed paraffin-embedded and frozen). We identified a 64-protein signature that predicts with high specificity a subset of HGSOCs refractory to initial platinum-based therapy and is validated in two independent patient cohorts. We detected significant association between lack of Ch17 loss of heterozygosity (LOH) and chemo-refractoriness. Based on pathway protein expression, we identified 5 clusters of HGSOC, which validated across two independent patient cohorts and patient-derived xenograft (PDX) models. These clusters may represent different mechanisms of refractoriness and implicate putative therapeutic vulnerabilities.

Comprehensive genetic and phenotype analysis of 95 individuals with mosaic tuberous sclerosis complex

Tuberous sclerosis complex (TSC) is a neurogenetic disorder due to loss-of-function TSC1 or TSC2 variants, characterized by tumors affecting multiple organs, including skin, brain, heart, lung, and kidney. Mosaicism for TSC1 or TSC2 variants occurs in 10%-15% of individuals diagnosed with TSC. Here, we report comprehensive characterization of TSC mosaicism by using massively parallel sequencing (MPS) of 330 TSC samples from a variety of tissues and fluids from a cohort of 95 individuals with mosaic TSC. TSC1 variants in individuals with mosaic TSC are much less common (9%) than in germline TSC overall (26%) (p < 0.0001). The mosaic variant allele frequency (VAF) is significantly higher in TSC1 than in TSC2, in both blood and saliva (median VAF: TSC1, 4.91%; TSC2, 1.93%; p = 0.036) and facial angiofibromas (median VAF: TSC1, 7.7%; TSC2 3.7%; p = 0.004), while the number of TSC clinical features in individuals with TSC1 and TSC2 mosaicism was similar. The distribution of mosaic variants across TSC1 and TSC2 is similar to that for pathogenic germline variants in general TSC. The systemic mosaic variant was not present in blood in 14 of 76 (18%) individuals with TSC, highlighting the value of analysis of multiple samples from each individual. A detailed comparison revealed that nearly all TSC clinical features are less common in individuals with mosaic versus germline TSC. A large number of previously unreported TSC1 and TSC2 variants, including intronic and large rearrangements (n = 11), were also identified.

Uncovering transcriptomic landscape alterations of CAN-2409 in in vitro and in vivo glioma models

Rationale

CAN-2409 is a locally delivered oncolytic therapy, which results in vaccination against the injected tumor. CAN-2409 consists of a non-replicating adenovirus armed with the Herpes virus thymidine kinase, which metabolizes ganciclovir into a phosphorylated nucleotide that is incorporated into the tumor cell's genome, thereby inflicting immunogenic cancer cell death. While CAN-2409's immunological impact has been well characterized, its effects on the tumor cells transcriptome remains unknown. We compared the transcriptomic landscape after treatment of glioblastoma models with CAN-2409 in vitro and in vivo to assess how the interplay with the tumor microenvironment influences CAN-2409-mediated transcriptome alterations.

Methods

We performed RNA-Seq with CAN-2409 treated patient-derived glioma stem-like cells and tumors of C57/BL6 mice and compared KEGG pathway usage and differential gene expression focusing on immune cell and cytokine profiles. T-cell -killing assays were performed to assess candidate effectors.

Results

PCA analysis showed distinct clustering of control and CAN-2409 samples under both conditions. KEGG pathway analysis revealed significant enrichment for p53 signaling and cell cycle pathway, with similar dynamics for key regulators of both pathways in vitro and in vivo, including MYC, CCNB1, PLK1 and CDC20. Selected alterations (PLK1 and CCNB1) were validated at the protein level. Cytokine expression analysis revealed upregulation of pro-inflammatory IL12a under both conditions; immune cell gene profiling showed reduction of myeloid associated genes. T-cell-killing assays showed increased killing in the presence of IL-12.

Conclusion

CAN-2409 significantly alters the transcriptome both in vitro and in vivo. Comparison of pathway enrichment revealed mutual and differential utilization of pathways under both conditions, suggesting a modulating influence on the cell cycle in tumor cells, and of the tumor microenvironment on the transcriptome in vivo. IL-12 synthesis likely depends on interactions with the tumor microenvironment, and it facilitates CAN-2409 cell killing. This dataset provides potential to understand resistance mechanisms and identify potential biomarkers for future studies.

Telomerase reactivation induces progression of mouse Braf V600E -driven thyroid cancers without telomere lengthening

Mutations in the promoter of the telomerase reverse transcriptase ( TERT ) gene are the paradigm of a cross-cancer alteration in a non-coding region. TERT promoter mutations (TPMs) are biomarkers of poor prognosis in several tumors, including thyroid cancers. TPMs enhance TERT transcription, which is otherwise silenced in adult tissues, thus reactivating a bona fide oncoprotein. To study TERT deregulation and its downstream consequences, we generated a Tert mutant promoter mouse model via CRISPR/Cas9 engineering of the murine equivalent locus (Tert -123C>T ) and crossed it with thyroid-specific Braf V600E -mutant mice. We also employed an alternative model of Tert overexpression (K5-Tert). Whereas all Braf V600E animals developed well-differentiated papillary thyroid tumors, 29% and 36% of Braf V600E +Tert -123C>T and Braf V600E +K5-Tert mice progressed to poorly differentiated thyroid cancers at week 20, respectively. Braf+Tert tumors showed increased mitosis and necrosis in areas of solid growth, and older animals from these cohorts displayed anaplastic-like features, i.e., spindle cells and macrophage infiltration. Murine Tert promoter mutation increased Tert transcription in vitro and in vivo , but temporal and intra-tumoral heterogeneity was observed. RNA-sequencing of thyroid tumor cells showed that processes other than the canonical Tert-mediated telomere maintenance role operate in these specimens. Pathway analysis showed that MAPK and PI3K/AKT signaling, as well as processes not previously associated with this tumor etiology, involving cytokine and chemokine signaling, were overactivated. Braf+Tert animals remained responsive to MAPK pathway inhibitors. These models constitute useful pre-clinical tools to understand the cell-autonomous and microenvironment-related consequences of Tert-mediated progression in advanced thyroid cancers and other aggressive tumors carrying TPMs.

Splicing factor SRSF1 controls autoimmune-related molecular pathways in regulatory T cells distinct from FoxP3

Regulatory T cells (Tregs) are vital for maintaining immune self-tolerance, and their impaired function leads to autoimmune disease. Mutations in FoxP3, the master transcriptional regulator of Tregs, leads to immune dysregulation, polyendocrinopathy, enteropathy, X-linked (IPEX) syndrome in humans and the early lethal "scurfy" phenotype with multi-organ autoimmune disease in mice. We recently identified serine/arginine-rich splicing factor 1 (SRSF1) as an indispensable regulator of Treg homeostasis and function. Intriguingly, Treg-conditional SRSF1-deficient mice exhibit early lethal systemic autoimmunity with multi-organ inflammation reminiscent of the scurfy mice. Importantly, SRSF1 is decreased in T cells from patients with the autoimmune disease systemic lupus erythematosus (SLE), and low SRSF1 levels inversely correlate with disease severity. Given that the Treg-specific deficiency of SRSF1 causes similarly profound autoimmune disease outcomes in mice as the deficiency/mutation in FoxP3, we aimed to evaluate the genes and molecular pathways controlled by these two indispensable regulatory proteins. We performed comparative bioinformatic analyses of transcriptomic profiles of Tregs from Srsf1-knockout mice and two Foxp3 mutant mice--the FoxP3-deficient ΔFoxp3 and the Foxp3 M370I mutant mice. We identified 132 differentially expressed genes (DEGs) unique to Srsf1-ko Tregs, 503 DEGs unique to Foxp3 M370I Tregs, and 1367 DEGs unique to ΔFoxp3 Tregs. Gene set enrichment and pathway analysis of DEGs unique to Srsf1-ko Tregs indicate that SRSF1 controls cytokine and immune response pathways. Conversely, FoxP3 controls pathways involved in DNA replication and cell cycle. Besides the distinct gene signatures, we identified only 30 shared genes between all three Treg mutants, mostly contributing to cytokine and immune defense pathways. Prominent genes included the chemokines CXCR6 and CCL1 and the checkpoint inhibitors FASLG and PDCD1. Thus, we demonstrate that SRSF1 and FoxP3 control common and distinct molecular pathways implicated in autoimmunity. Our analyses suggest that SRSF1 controls crucial immune functions in Tregs contributing to immune tolerance, and perturbations in its levels lead to systemic autoimmunity via mechanisms that are largely distinct from FoxP3.

Splicing factor SRSF1 is essential for CD8 T cell function and host antigen-specific viral immunity

Cytotoxic CD8 T cells are crucial for the host antigen-specific immune response to viral pathogens. Here we report the identification of an essential role for the serine/arginine-rich splicing factor (SRSF) 1 in CD8 T cell homeostasis and function. Specifically, SRSF1 is necessary for the maintenance of normal CD8 T lymphocyte numbers in the lymphoid compartment, and for the proliferative capacity and cytotoxic function of CD8 T cells. Furthermore, SRSF1 is required for antigen-specific IFN-γ cytokine responses in a viral infection challenge in mice. Transcriptomics analyses of Srsf1-deficient T cells reveal that SRSF1 controls proliferation, MAP kinase signaling and IFN signaling pathways. Mechanistically, SRSF1 controls the expression and activity of the Mnk2/p38-MAPK axis at the molecular level. Our findings reveal previously unrecognized roles for SRSF1 in the physiology and function of cytotoxic CD8 T lymphocytes and a potential molecular mechanism in viral immunopathogenesis.

Heme oxygenase-1 mitigates liver injury and fibrosis via modulation of LNX1/Notch1 pathway in myeloid cells

Activation of resident macrophages (Mϕ) and hepatic stellate cells is a key event in chronic liver injury. Mice with heme oxygenase-1 (HO-1; Hmox1)-deficient Mϕ (LysM-Cre:Hmox1 ^flfl ) exhibit increased inflammation, periportal ductular reaction, and liver fibrosis following bile duct ligation (BDL)-induced liver injury and increased pericellular fibrosis in NASH model. RiboTag-based RNA-sequencing profiling of hepatic HO-1-deficient Mϕ revealed dysregulation of multiple genes involved in lipid and amino acid metabolism, regulation of oxidative stress, and extracellular matrix turnover. Among these genes, ligand of numb-protein X1 (LNX1) expression is strongly suppressed in HO-1-deficient Mϕ. Importantly, HO-1 and LNX1 were expressed by hepatic Mϕ in human biliary and nonbiliary end-stage cirrhosis. We found that Notch1 expression, a downstream target of LNX1, was increased in LysM-Cre:Hmox1 ^flfl mice. In HO-1-deficient Mϕ treated with heme, transient overexpression of LNX1 drives M2-like Mϕ polarization. In summary, we identified LNX1/Notch1 pathway as a downstream target of HO-1 in liver fibrosis.

Ultrasensitive profiling of UV mutations identifies thousands of subclinical facial tumors in tuberous sclerosis complex

BACKGROUND

Tuberous Sclerosis Complex (TSC) is a neurogenetic syndrome due to loss-of-function mutations in TSC2 or TSC1, characterized by tumors at multiple body sites, including facial angiofibroma (FAF). Here, an ultrasensitive assessment of the extent and range of UV-induced mutations in TSC facial skin was performed.

METHODS

A Multiplex High-sensitivity PCR Assay (MHPA) was developed, enabling mutation detection at extremely low (<0.1%) variant allele frequencies (VAF).

RESULTS

MHPA assays were developed for both TSC2 and TP53, and applied to 81 samples, including 66 skin biopsies. UV-induced second hit mutation causing inactivation of TSC2 was pervasive in TSC facial skin with an average of 4.8 mutations per 2 mm biopsy at median VAF 0.08%, generating >150,000 incipient facial tumors (subclinical 'micro-FAFs') in the average TSC subject. The MHPA analysis also led to the identification of a refined UV-related indel signature and a recurrent complex mutation pattern, consisting of both a single or dinucleotide variant, and a 1-9 nt deletion, in cis.

CONCLUSION

TSC facial skin can be viewed as harboring a patchwork of clonal fibroblast proliferations (micro-FAF) with indolent growth, a small proportion of which develop into clinically observable FAF. Our observations also expand the spectrum of UV-related mutation signatures.

FUNDING

This work was supported by the TSC Alliance, Engles Family Fund for Research in TSC and LAM, and National Institutes of Health, National Heart, Lung, and Blood Institute [U01HL131022-04; Intramural Research Program].

A distinct core regulatory module enforces oncogene expression in KMT2A-rearranged leukemia

In this study, Harada et al. identified the transcription factors MEF2D and IRF8 as selective transcriptional dependencies of KMT2A-rearranged AML, where MEF2D displays partially redundant functions with its paralog, MEF2C. This study illustrates a mechanism of context-specific transcriptional addiction whereby a specific AML subclass depends on a highly specialized core regulatory module to directly enforce expression of common leukemia oncogenes.

Acute myeloid leukemia with KMT2A (MLL) rearrangements is characterized by specific patterns of gene expression and enhancer architecture, implying unique core transcriptional regulatory circuitry. Here, we identified the transcription factors MEF2D and IRF8 as selective transcriptional dependencies of KMT2A-rearranged AML, where MEF2D displays partially redundant functions with its paralog, MEF2C. Rapid transcription factor degradation followed by measurements of genome-wide transcription rates and superresolution microscopy revealed that MEF2D and IRF8 form a distinct core regulatory module with a narrow direct transcriptional program that includes activation of the key oncogenes MYC, HOXA9, and BCL2. Our study illustrates a mechanism of context-specific transcriptional addiction whereby a specific AML subclass depends on a highly specialized core regulatory module to directly enforce expression of common leukemia oncogenes.

Astrocyte-neuron crosstalk through Hedgehog signaling mediates cortical synapse development

Neuron-glia interactions play a critical role in the regulation of synapse formation and circuit assembly. Here we demonstrate that canonical Sonic hedgehog (Shh) pathway signaling in cortical astrocytes acts to coordinate layer-specific synaptic connectivity. We show that the Shh receptor Ptch1 is expressed by cortical astrocytes during development and that Shh signaling is necessary and sufficient to promote the expression of genes involved in regulating synaptic development and layer-enriched astrocyte molecular identity. Loss of Shh in layer V neurons reduces astrocyte complexity and coverage by astrocytic processes in tripartite synapses; conversely, cell-autonomous activation of Shh signaling in astrocytes promotes cortical excitatory synapse formation. Furthermore, Shh-dependent genes Lrig1 and Sparc distinctively contribute to astrocyte morphology and synapse formation. Together, these results suggest that Shh secreted from deep-layer cortical neurons acts to specialize the molecular and functional features of astrocytes during development to shape circuit assembly and function.

Xie et al. show that Sonic hedgehog expressed by cortical layer V neurons signals to Ptch1-expressing astrocytes where it is necessary to promote the expression of layer-enriched astrocyte gene networks that control synapse formation and astrocyte morphological complexity.

Splicing factor SRSF1 controls distinct molecular programs in regulatory and effector T cells implicated in systemic autoimmune disease

Systemic autoimmune diseases are characterized by hyperactive effector T cells (Teffs), aberrant cytokines and chemokines, and dysfunctional regulatory T cells (Tregs). We previously uncovered new roles for serine/arginine-rich splicing factor 1 (SRSF1) in the control of genes involved in T cell signaling and cytokine production in human T cells. SRSF1 levels are decreased in T cells from patients with systemic lupus erythematosus (SLE), and low levels correlate with severe disease. Moreover, T cell-conditional Srsf1-deficient mice recapitulate the autoimmune phenotype, exhibiting CD4 T cell hyperactivity, dysfunctional Tregs, systemic autoimmunity, and tissue inflammation. However, the role of SRSF1 in controlling molecular programs in Teffs and Tregs and how these pathways are implicated in autoimmunity is not known. Here, by comparative bioinformatics analysis, we demonstrate that SRSF1 controls largely distinct gene programs in Tregs and Teffs in vivo. SRSF1 regulates 189 differentially expressed genes (DEGs) unique to Tregs, 582 DEGs unique to Teffs, and 29 DEGs shared between both. Shared genes included IL-17A, IL-17F, CSF1, CXCL10, and CXCR4, and were highly enriched for inflammatory response and cytokine-cytokine receptor interaction pathways. SRSF1 controls distinct pathways in Tregs, which include chemokine signaling and immune cell differentiation, compared with pathways in Teffs, which include cytokine production, T cell homeostasis, and activation. We identified putative mRNA binding targets of SRSF1 which include CSF1, CXCL10, and IL-17F. Finally, comparisons with transcriptomics profiles from lupus-prone MRL/lpr mice reveal that SRSF1 controls genes and pathways implicated in autoimmune disease. The target genes of SRSF1 and putative binding targets we discovered, have known roles in systemic autoimmunity. Our findings suggest that SRSF1 controls distinct molecular pathways in Tregs and Teffs and aberrant SRSF1 levels may contribute to their dysfunction and immunopathogenesis of systemic autoimmune disease.

Multiomic analysis identifies CPT1A as a potential therapeutic target in platinum-refractory, high-grade serous ovarian cancer

Resistance to platinum compounds is a major determinant of patient survival in high-grade serous ovarian cancer (HGSOC). To understand mechanisms of platinum resistance and identify potential therapeutic targets in resistant HGSOC, we generated a data resource composed of dynamic (±carboplatin) protein, post-translational modification, and RNA sequencing (RNA-seq) profiles from intra-patient cell line pairs derived from 3 HGSOC patients before and after acquiring platinum resistance. These profiles reveal extensive responses to carboplatin that differ between sensitive and resistant cells. Higher fatty acid oxidation (FAO) pathway expression is associated with platinum resistance, and both pharmacologic inhibition and CRISPR knockout of carnitine palmitoyltransferase 1A (CPT1A), which represents a rate limiting step of FAO, sensitize HGSOC cells to platinum. The results are further validated in patient-derived xenograft models, indicating that CPT1A is a candidate therapeutic target to overcome platinum resistance. All multiomic data can be queried via an intuitive gene-query user interface (https://sites.google.com/view/ptrc-cell-line).

Performance assessment of DNA sequencing platforms in the ABRF Next-Generation Sequencing Study

Assessing the reproducibility, accuracy and utility of massively parallel DNA sequencing platforms remains an ongoing challenge. Here the Association of Biomolecular Resource Facilities (ABRF) Next-Generation Sequencing Study benchmarks the performance of a set of sequencing instruments (HiSeq/NovaSeq/paired-end 2 × 250-bp chemistry, Ion S5/Proton, PacBio circular consensus sequencing (CCS), Oxford Nanopore Technologies PromethION/MinION, BGISEQ-500/MGISEQ-2000 and GS111) on human and bacterial reference DNA samples. Among short-read instruments, HiSeq 4000 and X10 provided the most consistent, highest genome coverage, while BGI/MGISEQ provided the lowest sequencing error rates. The long-read instrument PacBio CCS had the highest reference-based mapping rate and lowest non-mapping rate. The two long-read platforms PacBio CCS and PromethION/MinION showed the best sequence mapping in repeat-rich areas and across homopolymers. NovaSeq 6000 using 2 × 250-bp read chemistry was the most robust instrument for capturing known insertion/deletion events. This study serves as a benchmark for current genomics technologies, as well as a resource to inform experimental design and next-generation sequencing variant calling.

Intra- and extra-cellular environments contribute to the fate of HIV-1 infection

HIV-1 entry into host cells leads to one of the following three alternative fates: (1) HIV-1 elimination by restriction factors, (2) establishment of HIV-1 latency, or (3) active viral replication in target cells. Here, we report the development of an improved system for monitoring HIV-1 fate at single-cell and population levels and show the diverse applications of this system to study specific aspects of HIV-1 fate in different cell types and under different environments. An analysis of the transcriptome of infected, primary CD4+ T cells that support alternative fates of HIV-1 identifies differential gene expression signatures in these cells. Small molecules are able to selectively target cells that support viral replication with no significant effect on viral latency. In addition, HIV-1 fate varies in different tissues following infection of humanized mice in vivo. Altogether, these studies indicate that intra- and extra-cellular environments contribute to the fate of HIV-1 infection.

Cross-Site Evaluation of Commercial Sanger Sequencing Chemistries

Sanger sequencing remains an essential tool utilized by researchers. Despite competition from commercial sequencing providers, many academic sequencing core facilities continue to offer these services based on a model of competitive pricing, knowledgeable technical support, and rapid turnaround time. In-house Sanger sequencing remains a viable core service and, until recently, Applied Biosystems BigDye Terminator chemistry was the only commercially available solution for Sanger DNA sequencing on Applied Biosystems (ABI) instruments; however, several new products employing novel dye chemistries and reaction configurations have entered the market. As a result, there is a need to benchmark the performance of these new chemistries on various DNA templates, including difficult-to-sequence templates, and their amenability to commonly employed cost-saving measures, such as dye dilution and reaction miniaturization. To evaluate these new reagents, a study was designed to compare the quality of Sanger sequencing data produced by ABI BigDye and commercially available kits from 2 other vendors using both control and difficult-to-sequence DNA templates under various reaction conditions. This study will serve as a valuable resource to core facilities conducting Sanger sequencing that wish to evaluate the use of an alternative chemistry in their sequencing core.

Multisite Evaluation of Next-Generation Methods for Small RNA Quantification

Small RNAs (smRNAs) are important regulators of many biologic processes and are now most frequently characterized using Illumina sequencing. However, although standard RNA sequencing library preparation has become routine in most sequencing facilities, smRNA sequencing library preparation has historically been challenging because of high input requirements, laborious protocols involving gel purifications, inability to automate, and a lack of benchmarking standards. Additionally, studies have suggested that many of these methods are nonlinear and do not accurately reflect the amounts of smRNAs in vivo. Recently, a number of new kits have become available that permit lower input amounts and less laborious, gel-free protocol options. Several of these new kits claim to reduce RNA ligase-dependent sequence bias through novel adapter modifications and to lessen adapter-dimer contamination in the resulting libraries. With the increasing number of smRNA kits available, understanding the relative strengths of each method is crucial for appropriate experimental design. In this study, we systematically compared 9 commercially available smRNA library preparation kits as well as NanoString probe hybridization across multiple study sites. Although several of the new methodologies do reduce the amount of artificially over- and underrepresented microRNAs (miRNAs), we observed that none of the methods was able to remove all of the bias in the library preparation. Identical samples prepared with different methods show highly varied levels of different miRNAs. Even so, many methods excelled in ease of use, lower input requirement, fraction of usable reads, and reproducibility across sites. These differences may help users select the most appropriate methods for their specific question of interest.

Therapeutic Targeting of the Secreted Lysophospholipase D Autotaxin Suppresses Tuberous Sclerosis Complex-Associated Tumorigenesis

Tuberous sclerosis complex (TSC) is an autosomal dominant disease characterized by multiorgan hamartomas, including renal angiomyolipomas and pulmonary lymphangioleiomyomatosis (LAM). TSC2 deficiency leads to hyperactivation of mTOR Complex 1 (mTORC1), a master regulator of cell growth and metabolism. Phospholipid metabolism is dysregulated upon TSC2 loss, causing enhanced production of lysophosphatidylcholine (LPC) species by TSC2-deficient tumor cells. LPC is the major substrate of the secreted lysophospholipase D autotaxin (ATX), which generates two bioactive lipids, lysophosphatidic acid (LPA) and sphingosine-1-phosphate (S1P). We report here that ATX expression is upregulated in human renal angiomyolipoma-derived TSC2-deficient cells compared with TSC2 add-back cells. Inhibition of ATX via the clinically developed compound GLPG1690 suppressed TSC2-loss associated oncogenicity in vitro and in vivo and induced apoptosis in TSC2-deficient cells. GLPG1690 suppressed AKT and ERK1/2 signaling and profoundly impacted the transcriptome of these cells while inducing minor gene expression changes in TSC2 add-back cells. RNA-sequencing studies revealed transcriptomic signatures of LPA and S1P, suggesting an LPA/S1P-mediated reprogramming of the TSC lipidome. In addition, supplementation of LPA or S1P rescued proliferation and viability, neutral lipid content, and AKT or ERK1/2 signaling in human TSC2-deficient cells treated with GLPG1690. Importantly, TSC-associated renal angiomyolipomas have higher expression of LPA receptor 1 and S1P receptor 3 compared with normal kidney. These studies increase our understanding of TSC2-deficient cell metabolism, leading to novel potential therapeutic opportunities for TSC and LAM. SIGNIFICANCE: This study identifies activation of the ATX-LPA/S1P pathway as a novel mode of metabolic dysregulation upon TSC2 loss, highlighting critical roles for ATX in TSC2-deficient cell fitness and in TSC tumorigenesis.

The KDM5A/RBP2 histone demethylase represses NOTCH signaling to sustain neuroendocrine differentiation and promote small cell lung cancer tumorigenesis

More than 90% of small cell lung cancers (SCLCs) harbor loss-of-function mutations in the tumor suppressor gene RB1 The canonical function of the RB1 gene product, pRB, is to repress the E2F transcription factor family, but pRB also functions to regulate cellular differentiation in part through its binding to the histone demethylase KDM5A (also known as RBP2 or JARID1A). We show that KDM5A promotes SCLC proliferation and SCLC's neuroendocrine differentiation phenotype in part by sustaining expression of the neuroendocrine transcription factor ASCL1. Mechanistically, we found that KDM5A sustains ASCL1 levels and neuroendocrine differentiation by repressing NOTCH2 and NOTCH target genes. To test the role of KDM5A in SCLC tumorigenesis in vivo, we developed a CRISPR/Cas9-based mouse model of SCLC by delivering an adenovirus (or an adeno-associated virus [AAV]) that expresses Cre recombinase and sgRNAs targeting Rb1, Tp53, and Rbl2 into the lungs of Lox-Stop-Lox Cas9 mice. Coinclusion of a KDM5A sgRNA decreased SCLC tumorigenesis and metastasis, and the SCLCs that formed despite the absence of KDM5A had higher NOTCH activity compared to KDM5A ^+/+ SCLCs. This work establishes a role for KDM5A in SCLC tumorigenesis and suggests that KDM5 inhibitors should be explored as treatments for SCLC.

Abstract 3647: Single cell RNA sequencing reveals mitogenic and progenitor gene programs inBRAF-rearranged pilocytic astrocytomas

Single-cell RNA sequencing (scRNAseq) has been performed across a range of intermediate to high-grade gliomas that each harbor multiple driver mutations. Pilocytic astrocytoma (PA), the most common childhood brain tumor, is a low-grade glioma. PAs frequently harbor oncogenic KIAA1549-BRAF fusions, and exhibit low rates of other driver mutations. While PAs exhibit a favorable prognosis compared to the higher-grade gliomas, treatment morbidity and tumor recurrence can represent major challenges for some PA patients. We performed scRNAseq of both tumor and non-tumor cells in newly-diagnosed PAs that contained KIAA1549-BRAF rearrangements. To confidently distinguish tumor cells from nontumor cells, we sorted cells by glial progenitor marker A2B5 status and profiled KIAA1549-BRAF fusion status of cells using a sensitive quantitative PCR-based assay. Results were validated using RNA in situ hybridization. When compared to higher-grade gliomas, a higher proportion of the PA tumor cells exhibited a differentiated, astrocyte-like phenotype. A smaller proportion of cells exhibited a progenitor-like phenotype with evidence of proliferation. These progenitor-like tumor cells expressed a mitogen-activated protein kinase (MAPK) program that was absent from higher-grade gliomas. Similar patterns of expression of genes associated with the astrocyte-like and MAPK gene programs were also seen in formalin fixed, paraffin embedded PA tissues using RNA in situ hybridization. Immune cells, especially microglia, comprised almost half of all cells in the PAs and accounted for differences in bulk expression profiles between tumor locations and subtypes. These single cell transcriptional data reveal a transcriptional developmental hierarchy in a pediatric low grade glioma that is skewed towards more mature brain cells compared to higher-grade gliomas. The results indicate that future analyses of bulk PA tissues should attempt to account for considerable infiltration by nontumor cells. Finally, the finding that a MAPK gene program is not uniformly expressed in PA tumor cells has implications for ongoing clinical investigations of therapies directed at the MAPK pathway for the treatment of PA.

Citation Format: Zachary J. Reitman, Brenton Paolella, Guillaume Bergthold, Kristine Pelton, Robert Jones, Sarah Becker, Claire E. Sinai, Haley Malkin, Ying Huang, Leslie Grimmett, Zachary T. Herbert, Yu Sun, Jessica Weatherbee, John Alberta, John Daley, Orit Rozenblatt-Rosen, Rosalind Segal, Daphne Haas-Kogan, Mariella G. Filbin, Mario L. Suva, Aviv Regev, Charles Stiles, Mark W. Kieran, Liliana Goumnerova, Keith L. Ligon, Alex K. Shalek, Pratiti Bandopadhayay, Rameen Beroukhim. Single cell RNA sequencing reveals mitogenic and progenitor gene programs inBRAF-rearranged pilocytic astrocytomas [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 3647.

Histone demethylase KDM6A directly senses oxygen to control chromatin and cell fate

Oxygen sensing is central to metazoan biology and has implications for human disease. Mammalian cells express multiple oxygen-dependent enzymes called 2-oxoglutarate (OG)-dependent dioxygenases (2-OGDDs), but they vary in their oxygen affinities and hence their ability to sense oxygen. The 2-OGDD histone demethylases control histone methylation. Hypoxia increases histone methylation, but whether this reflects direct effects on histone demethylases or indirect effects caused by the hypoxic induction of the HIF (hypoxia-inducible factor) transcription factor or the 2-OG antagonist 2-hydroxyglutarate (2-HG) is unclear. Here, we report that hypoxia promotes histone methylation in a HIF- and 2-HG-independent manner. We found that the H3K27 histone demethylase KDM6A/UTX, but not its paralog KDM6B, is oxygen sensitive. KDM6A loss, like hypoxia, prevented H3K27 demethylation and blocked cellular differentiation. Restoring H3K27 methylation homeostasis in hypoxic cells reversed these effects. Thus, oxygen directly affects chromatin regulators to control cell fate.

Cross-site comparison of ribosomal depletion kits for Illumina RNAseq library construction

Background

Ribosomal RNA (rRNA) comprises at least 90% of total RNA extracted from mammalian tissue or cell line samples. Informative transcriptional profiling using massively parallel sequencing technologies requires either enrichment of mature poly-adenylated transcripts or targeted depletion of the rRNA fraction. The latter method is of particular interest because it is compatible with degraded samples such as those extracted from FFPE and also captures transcripts that are not poly-adenylated such as some non-coding RNAs. Here we provide a cross-site study that evaluates the performance of ribosomal RNA removal kits from Illumina, Takara/Clontech, Kapa Biosystems, Lexogen, New England Biolabs and Qiagen on intact and degraded RNA samples.

Results

We find that all of the kits are capable of performing significant ribosomal depletion, though there are differences in their ease of use. All kits were able to remove ribosomal RNA to below 20% with intact RNA and identify ~ 14,000 protein coding genes from the Universal Human Reference RNA sample at >1FPKM. Analysis of differentially detected genes between kits suggests that transcript length may be a key factor in library production efficiency.

Conclusions

These results provide a roadmap for labs on the strengths of each of these methods and how best to utilize them.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4585-1) contains supplementary material, which is available to authorized users.

The association between germline BRCA2 variants and sensitivity to platinum-based chemotherapy among men with metastatic prostate cancer

BACKGROUND

Breast cancer 2 (BRCA2)-associated breast and ovarian cancers are sensitive to platinum-based chemotherapy. It is unknown whether BRCA2-associated prostate cancer responds favorably to such treatment.

METHODS

A retrospective analysis of a single-institution cohort of men with castration-resistant, metastatic prostate cancer was performed to determine the association between carrier status of pathogenic BRCA2 germline variants and prostate-specific antigen response to carboplatin-based chemotherapy. From 2001 through 2015, 8081 adult men with prostate cancer who had a consultation and/or underwent treatment at Dana-Farber Cancer Institute provided blood samples and consented to analyses of biologic material and clinical records. A subgroup of 141 men received at least 2 doses of carboplatin and docetaxel for castration-resistant disease (94% were also taxane refractory). These patients were categorized according to the absence or presence of pathogenic germline mutations in BRCA2 based on DNA sequencing from whole blood. The primary outcome was the response rate to carboplatin/docetaxel chemotherapy, defined according to a decline in prostate-specific antigen that exceeded 50% within 12 weeks of initiating this regimen. Associations between BRCA2 mutation status and response to carboplatin-based chemotherapy were tested using the Fisher exact test, with a 2-sided P value < .05 as the threshold for significance.

RESULTS

Pathogenic germline BRCA2 variants were observed in 8 of 141 men (5.7%; 95% confidence interval, 2.5%-10.9%). Six of 8 BRCA2 carriers (75%) experienced prostate-specific antigen declines >50% within 12 weeks, compared with 23 of 133 noncarriers (17%; absolute difference, 58%; 95% confidence interval, 27%-88%; P < .001). Prostate cancer cell lines functionally corroborated these clinical findings.

CONCLUSIONS

BRCA2-associated, castration-resistant prostate cancer is associated with a higher likelihood of response to carboplatin-based chemotherapy than non-BRCA2-associated prostate cancer. Cancer 2017;123:3532-9. © 2017 American Cancer Society.

The Public Repository of Xenografts Enables Discovery and Randomized Phase II-like Trials in Mice

More than 90% of drugs with preclinical activity fail in human trials, largely due to insufficient efficacy. We hypothesized that adequately powered trials of patient-derived xenografts (PDX) in mice could efficiently define therapeutic activity across heterogeneous tumors. To address this hypothesis, we established a large, publicly available repository of well-characterized leukemia and lymphoma PDXs that undergo orthotopic engraftment, called the Public Repository of Xenografts (PRoXe). PRoXe includes all de-identified information relevant to the primary specimens and the PDXs derived from them. Using this repository, we demonstrate that large studies of acute leukemia PDXs that mimic human randomized clinical trials can characterize drug efficacy and generate transcriptional, functional, and proteomic biomarkers in both treatment-naive and relapsed/refractory disease.