Ian Sudbery

Sequencing depth and coverage: key considerations in genomic analyses

Sequencing technologies have placed a wide range of genomic analyses within the capabilities of many laboratories. However, sequencing costs often set limits to the amount of sequences that can be generated and, consequently, the biological outcomes that can be achieved from an experimental design. In this Review, we discuss the issue of sequencing depth in the design of next-generation sequencing experiments. We review current guidelines and precedents on the issue of coverage, as well as their underlying considerations, for four major study designs, which include de novo genome sequencing, genome resequencing, transcriptome sequencing and genomic location analyses (for example, chromatin immunoprecipitation followed by sequencing (ChIP–seq) and chromosome conformation capture (3C)).

Apoptosis induced by environmental stresses and amphotericin B in Candida albicans

New antifungal agents are urgently required to combat life-threatening infections caused by opportunistic fungal pathogens like Candida albicans. The manipulation of endogenous fungal programmed cell death responses could provide a basis for future therapies. Here we assess the physiology of death in C. albicans in response to environmental stresses (acetic acid and hydrogen peroxide) and an antifungal agent (amphotericin B). Exposure of C. albicans to 40-60 mM acetic acid, 5-10 mM hydrogen peroxide, or 4-8 μg·ml-1 amphotericin 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. albicans cells 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 in C. albicans was 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.

KDM2B links the Polycomb Repressive Complex 1 (PRC1) to recognition of CpG islands

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. DOI: http://dx.doi.org/10.7554/eLife.00205.001

Next-generation sequencing of advanced prostate cancer treated with androgen-deprivation therapy.

BACKGROUND 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. OBJECTIVE 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. DESIGN, SETTING, AND PARTICIPANTS 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. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS 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. RESULTS AND LIMITATIONS 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. CONCLUSIONS 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.

UMI-tools: modeling sequencing errors in Unique Molecular Identifiers to improve quantification accuracy

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.

Next-generation sequencing of vertebrate experimental organisms

Next-generation sequencing technologies are revolutionizing biology by allowing for genome-wide transcription factor binding-site profiling, transcriptome sequencing, and more recently, whole-genome resequencing. While it is currently not possible to generate complete de novo assemblies of higher-vertebrate genomes using next-generation sequencing, improvements in sequence read lengths and throughput, coupled with new assembly algorithms for large data sets, will soon make this a reality. These developments will in turn spawn a revolution in how genomic data are used to understand genetics and how model organisms are used for disease gene discovery. This review provides an overview of the current next-generation sequencing platforms and the newest computational tools for the analysis of next-generation sequencing data. We also describe how next-generation sequencing may be applied in the context of vertebrate model organism genetics.

Deep short-read sequencing of chromosome 17 from the mouse strains A/J and CAST/Ei identifies significant germline variation and candidate genes that regulate liver triglyceride levels

Genome sequences are essential tools for comparative and mutational analyses. Here we present the short read sequence of mouse chromosome 17 from the Mus musculus domesticus derived strain A/J, and the Mus musculus castaneus derived strain CAST/Ei. We describe approaches for the accurate identification of nucleotide and structural variation in the genomes of vertebrate experimental organisms, and show how these techniques can be applied to help prioritize candidate genes within quantitative trait loci.

CGAT: computational genomics analysis toolkit

Summary: Computational genomics seeks to draw biological inferences from genomic datasets, often by integrating and contextualizing next-generation sequencing data. CGAT provides an extensive suite of tools designed to assist in the analysis of genome scale data from a range of standard file formats. The toolkit enables filtering, comparison, conversion, summarization and annotation of genomic intervals, gene sets and sequences. The tools can both be run from the Unix command line and installed into visual workflow builders, such as Galaxy. Availability: The toolkit is freely available from http://github.com/CGATOxford/cgat Contact: andreas.heger@dpag.ox.ac.uk

Cell cycle-independent phospho-regulation of Fkh2 during hyphal growth regulates Candida albicans pathogenesis.

The opportunistic human fungal pathogen, Candida albicans, undergoes morphological and transcriptional adaptation in the switch from commensalism to pathogenicity. Although previous gene-knockout studies have identified many factors involved in this transformation, it remains unclear how these factors are regulated to coordinate the switch. Investigating morphogenetic control by post-translational phosphorylation has generated important regulatory insights into this process, especially focusing on coordinated control by the cyclin-dependent kinase Cdc28. Here we have identified the Fkh2 transcription factor as a regulatory target of both Cdc28 and the cell wall biosynthesis kinase Cbk1, in a role distinct from its conserved function in cell cycle progression. In stationary phase yeast cells 2D gel electrophoresis shows that there is a diverse pool of Fkh2 phospho-isoforms. For a short window on hyphal induction, far before START in the cell cycle, the phosphorylation profile is transformed before reverting to the yeast profile. This transformation does not occur when stationary phase cells are reinoculated into fresh medium supporting yeast growth. Mass spectrometry and mutational analyses identified residues phosphorylated by Cdc28 and Cbk1. Substitution of these residues with non-phosphorylatable alanine altered the yeast phosphorylation profile and abrogated the characteristic transformation to the hyphal profile. Transcript profiling of the phosphorylation site mutant revealed that the hyphal phosphorylation profile is required for the expression of genes involved in pathogenesis, host interaction and biofilm formation. We confirmed that these changes in gene expression resulted in corresponding defects in pathogenic processes. Furthermore, we identified that Fkh2 interacts with the chromatin modifier Pob3 in a phosphorylation-dependent manner, thereby providing a possible mechanism by which the phosphorylation of Fkh2 regulates its specificity. Thus, we have discovered a novel cell cycle-independent phospho-regulatory event that subverts a key component of the cell cycle machinery to a role in the switch from commensalism to pathogenicity.

Phosphoregulation of Nap1 Plays a Role in Septin Ring Dynamics and Morphogenesis in Candida albicans

ABSTRACT Nap1 has long been identified as a potential septin regulator in yeasts. However, its function and regulation remain poorly defined. Here, we report functional characterization of Nap1 in the human-pathogenic fungus Candida albicans. We find that deletion of NAP1 causes constitutive filamentous growth and changes of septin dynamics. We present evidence that Nap1’s cellular localization and function are regulated by phosphorylation. Phos-tag gel electrophoresis revealed that Nap1 phosphorylation is cell cycle dependent, exhibiting the lowest level around the time of bud emergence. Mass spectrometry identified 10 phosphoserine and phosphothreonine residues in a cluster near the N terminus, and mutation of these residues affected Nap1’s localization to the septin ring and cellular function. Nap1 phosphorylation involves two septin ring-associated kinases, Cla4 and Gin4, and its dephosphorylation occurs at the septin ring in a manner dependent on the phosphatases PP2A and Cdc14. Furthermore, the nap1Δ/Δ mutant and alleles carrying mutations of the phosphorylation sites exhibited greatly reduced virulence in a mouse model of systemic candidiasis. Together, our findings not only provide new mechanistic insights into Nap1’s function and regulation but also suggest the potential to target Nap1 in future therapeutic design. IMPORTANCE Septins are conserved filament-forming GTPases involved in a wide range of cellular events, such as cytokinesis, exocytosis, and morphogenesis. In Candida albicans, the most prevalent human fungal pathogen, septin functions are indispensable for its virulence. However, the molecular mechanisms by which septin structures are regulated are poorly understood. In this study, we deleted NAP1, a gene encoding a putative septin regulator, in C. albicans and found that cells lacking NAP1 showed abnormalities in morphology, invasive growth, and septin ring dynamics. We identified a conserved N-terminal phosphorylation cluster on Nap1 and demonstrated that phosphorylation at these sites regulates Nap1 localization and function. Importantly, deletion of NAP1 or mutation in the N-terminal phosphorylation cluster strongly reduced the virulence of C. albicans in a mouse model of systemic infection. Thus, this study not only provides mechanistic insights into septin regulation but also suggests Nap1 as a potential antifungal target. Septins are conserved filament-forming GTPases involved in a wide range of cellular events, such as cytokinesis, exocytosis, and morphogenesis. In Candida albicans, the most prevalent human fungal pathogen, septin functions are indispensable for its virulence. However, the molecular mechanisms by which septin structures are regulated are poorly understood. In this study, we deleted NAP1, a gene encoding a putative septin regulator, in C. albicans and found that cells lacking NAP1 showed abnormalities in morphology, invasive growth, and septin ring dynamics. We identified a conserved N-terminal phosphorylation cluster on Nap1 and demonstrated that phosphorylation at these sites regulates Nap1 localization and function. Importantly, deletion of NAP1 or mutation in the N-terminal phosphorylation cluster strongly reduced the virulence of C. albicans in a mouse model of systemic infection. Thus, this study not only provides mechanistic insights into septin regulation but also suggests Nap1 as a potential antifungal target.