A statistical approach for identifying differential distributions in single-cell RNA-seq experiments

The ability to quantify cellular heterogeneity is a major advantage of single-cell technologies. However, statistical methods often treat cellular heterogeneity as a nuisance. We present a novel method to characterize differences in expression in the presence of distinct expression states within and among biological conditions. We demonstrate that this framework can detect differential expression patterns under a wide range of settings. Compared to existing approaches, this method has higher power to detect subtle differences in gene expression distributions that are more complex than a mean shift, and can characterize those differences. The freely available R package scDD implements the approach. Electronic supplementary material The online version of this article (doi:10.1186/s13059-016-1077-y) contains supplementary material, which is available to authorized users

High‐throughput identification of RNA nuclear enrichment sequences

Abstract In the post‐genomic era, thousands of putative noncoding regulatory regions have been identified, such as enhancers, promoters, long noncoding RNAs (lncRNAs), and a cadre of small peptides. These ever‐growing catalogs require high‐throughput assays to test their functionality at scale. Massively parallel reporter assays have greatly enhanced the understanding of noncoding DNA elements en masse. Here, we present a massively parallel RNA assay (MPRNA) that can assay 10,000 or more RNA segments for RNA‐based functionality. We applied MPRNA to identify RNA‐based nuclear localization domains harbored in lncRNAs. We examined a pool of 11,969 oligos densely tiling 38 human lncRNAs that were fused to a cytosolic transcript. After cell fractionation and barcode sequencing, we identified 109 unique RNA regions that significantly enriched this cytosolic transcript in the nucleus including a cytosine‐rich motif. These nuclear enrichment sequences are highly conserved and over‐represented in global nuclear fractionation sequencing. Importantly, many of these regions were independently validated by single‐molecule RNA fluorescence in situ hybridization. Overall, we demonstrate the utility of MPRNA for future investigation of RNA‐based functionalities

Differential substrate use in EGF‐ and oncogenic KRAS‐stimulated human mammary epithelial cells

Many metabolic phenotypes in cancer cells are also characteristic of proliferating nontransformed mammalian cells, and attempts to distinguish between phenotypes resulting from oncogenic perturbation from those associated with increased proliferation are limited. Here, we examined the extent to which metabolic changes corresponding to oncogenic KRAS expression differed from those corresponding to epidermal growth factor (EGF)-driven proliferation in human mammary epithelial cells (HMECs). Removal of EGF from culture medium reduced growth rates and glucose/glutamine consumption in control HMECs despite limited changes in respiration and fatty acid synthesis, while the relative contribution of branched-chain amino acids to the TCA cycle and lipogenesis increased in the near-quiescent conditions. Most metabolic phenotypes measured in HMECs expressing mutant KRAS were similar to those observed in EGF-stimulated control HMECs that were growing at comparable rates. However, glucose and glutamine consumption as well as lactate and glutamate production were lower in KRAS-expressing cells cultured in media without added EGF, and these changes correlated with reduced sensitivity to GLUT1 inhibitor and phenformin treatment. Our results demonstrate the strong dependence of metabolic behavior on growth rate and provide a model to distinguish the metabolic influences of oncogenic mutations and nononcogenic growth

Transparency and reproducibility in artificial intelligence

Breakthroughs in artificial intelligence (AI) hold enormous potential as it can automate complex tasks and go even beyond human performance. In their study, McKinney et al. showed the high potential of AI for breast cancer screening. However, the lack of methods’ details and algorithm code undermines its scientific value. Here, we identify obstacles hindering transparent and reproducible AI research as faced by McKinney et al., and provide solutions to these obstacles with implications for the broader field

Transparency and reproducibility in artificial intelligence

Breakthroughs in artificial intelligence (AI) hold enormous potential as it can automate complex tasks and go even beyond human performance. In their study, McKinney et al. showed the high potential of AI for breast cancer screening. However, the lack of methods’ details and algorithm code undermines its scientific value. Here, we identify obstacles hindering transparent and reproducible AI research as faced by McKinney et al., and provide solutions to these obstacles with implications for the broader field

CDK4/6 inhibition reprograms the breast cancer enhancer landscape by stimulating AP-1 transcriptional activity

Pharmacologic inhibitors of cyclin-dependent kinases 4 and 6 (CDK4/6) were designed to induce cancer cell cycle arrest. Recent studies have suggested that these agents also exert other effects, influencing cancer cell immunogenicity, apoptotic responses, and differentiation. Using cell-based and mouse models of breast cancer together with clinical specimens, we show that CDK4/6 inhibitors induce remodeling of cancer cell chromatin characterized by widespread enhancer activation, and that this explains many of these effects. The newly activated enhancers include classical super-enhancers that drive luminal differentiation and apoptotic evasion, as well as a set of enhancers overlying endogenous retroviral elements that is enriched for proximity to interferon-driven genes. Mechanistically, CDK4/6 inhibition increases the level of several Activator Protein-1 (AP-1) transcription factor proteins, which are in turn implicated in the activity of many of the new enhancers. Our findings offer insights into CDK4/6 pathway biology and should inform the future development of CDK4/6 inhibitors

Reversal of viral and epigenetic HLA class I repression in Merkel cell carcinoma

Cancers avoid immune surveillance through an array of mechanisms, including perturbation of HLA class I antigen presentation. Merkel cell carcinoma (MCC) is an aggressive, HLA-I–low, neuroendocrine carcinoma of the skin often caused by the Merkel cell polyomavirus (MCPyV). Through the characterization of 11 newly generated MCC patient-derived cell lines, we identified transcriptional suppression of several class I antigen presentation genes. To systematically identify regulators of HLA-I loss in MCC, we performed parallel, genome-scale, gain- and loss-of-function screens in a patient-derived MCPyV-positive cell line and identified MYCL and the non-canonical Polycomb repressive complex 1.1 (PRC1.1) as HLA-I repressors. We observed physical interaction of MYCL with the MCPyV small T viral antigen, supporting a mechanism of virally mediated HLA-I suppression. We further identify the PRC1.1 component USP7 as a pharmacologic target to restore HLA-I expression in MCC

Prostate cancer reactivates developmental epigenomic programs during metastatic progression

Epigenetic processes govern prostate cancer (PCa) biology, as evidenced by the dependency of PCa cells on the androgen receptor (AR), a prostate master transcription factor. We generated 268 epigenomic datasets spanning two state transitions—from normal prostate epithelium to localized PCa to metastases—in specimens derived from human tissue. We discovered that reprogrammed AR sites in metastatic PCa are not created de novo; rather, they are prepopulated by the transcription factors FOXA1 and HOXB13 in normal prostate epithelium. Reprogrammed regulatory elements commissioned in metastatic disease hijack latent developmental programs, accessing sites that are implicated in prostate organogenesis. Analysis of reactivated regulatory elements enabled the identification and functional validation of previously unknown metastasis-specific enhancers at HOXB13, FOXA1 and NKX3-1. Finally, we observed that prostate lineage-specific regulatory elements were strongly associated with PCa risk heritability and somatic mutation density. Examining prostate biology through an epigenomic lens is fundamental for understanding the mechanisms underlying tumor progression