Unique Molecular Identifiers (UMIs) are random oligonucleotide barcodes that are increasingly used in high-throughput sequencing experiments. Through a UMI, identical copies arising from distinct molecules can be distinguished from those arising through PCR amplification of the same molecule. However, bioinformatic methods to leverage the information from UMIs have yet to be formalized. In particular, sequencing errors in the UMI sequence are often ignored or else resolved in an ad hoc manner. We show that errors in the UMI sequence are common and introduce network-based methods to account for these errors when identifying PCR duplicates. Using these methods, we demonstrate improved quantification accuracy both under simulated conditions and real iCLIP and single-cell RNA-seq data sets. Reproducibility between iCLIP replicates and single-cell RNA-seq clustering are both improved using our proposed network-based method, demonstrating the value of properly accounting for errors in UMIs. These methods are implemented in the open source UMI-tools software package.
CpG islands (CGIs) are associated with most mammalian gene promoters. A subset of CGIs act as polycomb response elements (PREs) and are recognized by the polycomb silencing systems to regulate expression of genes involved in early development. How CGIs function mechanistically as nucleation sites for polycomb repressive complexes remains unknown. Here we discover that KDM2B (FBXL10) specifically recognizes non-methylated DNA in CGIs and recruits the polycomb repressive complex 1 (PRC1). This contributes to histone H2A lysine 119 ubiquitylation (H2AK119ub1) and gene repression. Unexpectedly, we also find that CGIs are occupied by low levels of PRC1 throughout the genome, suggesting that the KDM2B-PRC1 complex may sample CGI-associated genes for susceptibility to polycomb-mediated silencing. These observations demonstrate an unexpected and direct link between recognition of CGIs by KDM2B and targeting of the polycomb repressive system. This provides the basis for a new model describing the functionality of CGIs as mammalian PREs.
17 November 2003 in Proceedings of the National Academy of Sciences
New antifungal agents are urgently required to combat life-threatening infections caused by opportunistic fungal pathogens likeCandida albicans. The manipulation of endogenous fungal programmed cell death responses could provide a basis for future therapies. Here we assess the physiology of death inC. albicansin response to environmental stresses (acetic acid and hydrogen peroxide) and an antifungal agent (amphotericin B). Exposure ofC. albicansto 40-60 mM acetic acid, 5-10 mM hydrogen peroxide, or 4-8 μg·ml-1amphotericin B produced cellular changes reminiscent of mammalian apoptosis. Nonviable cells that excluded propidium iodide displayed the apoptotic marker phosphatidylserine (as shown by annexin-V-FITC labeling), were terminal deoxynucleotidyltransferase-mediated dUTP nick end labeling (TUNEL)-positive (indicating nuclease-mediated double-strand DNA breakage), and produced reactive oxygen species. Ultrastructural changes in apoptotic cells included chromatin condensation and margination, separation of the nuclear envelope, and nuclear fragmentation.C. albicanscells treated at higher doses of these compounds showed cellular changes characteristic of necrosis. Necrotic cells displayed reduced TUNEL staining, a lack of surface phosphatidylserine, limited reactive oxygen species production, and an inability to exclude propidium iodide. Necrotic cells lacked defined nuclei and showed extensive intracellular vacuolization. Apoptosis inC. albicanswas associated with an accumulation of cells in the G2/M phase of the cell cycle, and under some apoptosis-inducing conditions, significant proportions of yeast cells switched to hyphal growth before dying. This is a demonstration of apoptosis in a medically important fungal pathogen.
We introduce alevin, a fast end-to-end pipeline to process droplet-based single-cell RNA sequencing data, performing cell barcode detection, read mapping, unique molecular identifier (UMI) deduplication, gene count estimation, and cell barcode whitelisting. Alevin’s approach to UMI deduplication considers transcript-level constraints on the molecules from which UMIs may have arisen and accounts for both gene-unique reads and reads that multimap between genes. This addresses the inherent bias in existing tools which discard gene-ambiguous reads and improves the accuracy of gene abundance estimates. Alevin is considerably faster, typically eight times, than existing gene quantification approaches, while also using less memory.
Androgen-deprivation therapy (ADT) is standard treatment for locally advanced or metastatic prostate cancer (PCa). Many patients develop castration resistance (castration-resistant PCa [CRPC]) after approximately 2-3 yr, with a poor prognosis. The molecular mechanisms underlying CRPC progression are unclear. To undertake quantitative tumour transcriptome profiling prior to and following ADT to identify functionally important androgen-regulated pathways or genes that may be reactivated in CRPC. RNA sequencing (RNA-seq) was performed on tumour-rich, targeted prostatic biopsies from seven patients with locally advanced or metastatic PCa before and approximately 22 wk after ADT initiation. Differentially regulated genes were identified in treatment pairs and further investigated by quantitative reverse transcription-polymerase chain reaction (qRT-PCR) on cell lines and immunohistochemistry on a separate CRPC patient cohort. Functional assays were used to determine the effect of pathway modulation on cell phenotypes. We searched for gene expression changes affecting key cell signalling pathways that may be targeted as proof of principle in a CRPC in vitro cell line model. We identified ADT-regulated signalling pathways, including the Wnt/β-catenin signalling pathway, and observed overexpression of β-catenin in a subset of CRPC by immunohistochemistry. We validated 6 of 12 (50%) pathway members by qRT-PCR on LNCaP/LNCaP-AI cell RNAs, of which 4 (67%) demonstrated expression changes consistent with RNA-seq data. We show that the tankyrase inhibitor XAV939 (which promotes β-catenin degradation) reduced androgen-independent LNCaP-AI cell line growth compared with androgen-responsive LNCaP cells via an accumulation of cell proportions in the G0/G1 phase and reduction in the S and G2/M phases. Our biopsy protocol did not account for tumour heterogeneity, and pathway inhibition was limited to pharmacologic approaches. RNA-seq of paired PCa samples revealed ADT-regulated signalling pathways. Proof-of-principle inhibition of the Wnt/β-catenin signalling pathway specifically delays androgen-independent PCa cell cycle progression and proliferation and warrants further investigation as a potential target for therapy for CRPC.