CCAR1 promotes DNA repair via alternative splicing.

DNA repair is directly performed by hundreds of core factors and indirectly regulated by thousands of others. We massively expanded a CRISPR inhibition and Cas9-editing screening system to discover factors indirectly modulating homology-directed repair (HDR) in the context of ∼18,000 individual gene knockdowns. We focused on CCAR1, a poorly understood gene that we found the depletion of reduced both HDR and interstrand crosslink repair, phenocopying the loss of the Fanconi anemia pathway. CCAR1 loss abrogated FANCA protein without substantial reduction in the level of its mRNA or that of other FA genes. We instead found that CCAR1 prevents inclusion of a poison exon in FANCA. Transcriptomic analysis revealed that the CCAR1 splicing modulatory activity is not limited to FANCA, and it instead regulates widespread changes in alternative splicing that would damage coding sequences in mouse and human cells. CCAR1 therefore has an unanticipated function as a splicing fidelity factor.

PnB Designer: a web application to design prime and base editor guide RNAs for animals and plants.

BACKGROUND: The rapid expansion of the CRISPR toolbox through tagging effector domains to either enzymatically inactive Cas9 (dCas9) or Cas9 nickase (nCas9) has led to several promising new gene editing strategies. Recent additions include CRISPR cytosine or adenine base editors (CBEs and ABEs) and the CRISPR prime editors (PEs), in which a deaminase or reverse transcriptase are fused to nCas9, respectively. These tools hold great promise to model and correct disease-causing mutations in animal and plant models. But so far, no widely-available tools exist to automate the design of both BE and PE reagents. RESULTS: We developed PnB Designer, a web-based application for the design of pegRNAs for PEs and guide RNAs for BEs. PnB Designer makes it easy to design targeting guide RNAs for single or multiple targets on a variant or reference genome from organisms spanning multiple kingdoms. With PnB Designer, we designed pegRNAs to model all known disease causing mutations available in ClinVar. Additionally, PnB Designer can be used to design guide RNAs to install or revert a SNV, scanning the genome with one CBE and seven different ABE PAM variants and returning the best BE to use. PnB Designer is publicly accessible at http://fgcz-shiny.uzh.ch/PnBDesigner/ CONCLUSION: With PnB Designer we created a user-friendly design tool for CRISPR PE and BE reagents, which should simplify choosing editing strategy and avoiding design complications.

Effect of progesterone on Candida albicans biofilm formation under acidic conditions: A transcriptomic analysis.

Vulvovaginal candidiasis (VVC) caused by Candida albicans is a common disease worldwide. A very important C. albicans virulence factor is its ability to form biofilms on epithelium and/or on intrauterine devices promoting VVC. It has been shown that VVC has a hormonal dependency and that progesterone affects virulence traits of C. albicans cells. To understand how the acidic environment (pH 4) and progesterone (either alone and in combination) modulate C. albicans response during formation of biofilm, a transcriptomic analysis was performed together with characterization of the biofilm properties. Compared to planktonic cells, acidic biofilm-cells exhibited major changes in their transcriptome, including modifications in the expression of 286 genes that were not previously associated with biofilm formation in C. albicans. The vast majority of the genes up-regulated in the acidic biofilm cells (including those uniquely identified in our study) are known targets of Sfl1, and consistently, Sfl1 deletion is herein shown to impair the formation of acidic biofilms (pH 4). Under the acidic conditions used, the presence of progesterone reduced C. albicans biofilm biomass and structural cohesion. Transcriptomic analysis of biofilms developed in the presence of progesterone led to the identification of 65 down-regulated genes including, among others, the regulator Tec1 and several of its target genes, suggesting that the function of this transcription factor is inhibited by the presence of the hormone. Additionally, progesterone reduced the susceptibility of biofilm cells to fluconazole, consistent with an up-regulation of efflux pumps. Overall, the results of this study show that progesterone modulates C. albicans biofilm formation and genomic expression under acidic conditions, which may have implications for C. albicans pathogenicity in the vaginal environment.

Proteogenomics produces comprehensive and highly accurate protein-coding gene annotation in a complete genome assembly of Malassezia sympodialis.

Complete and accurate genome assembly and annotation is a crucial foundation for comparative and functional genomics. Despite this, few complete eukaryotic genomes are available, and genome annotation remains a major challenge. Here, we present a complete genome assembly of the skin commensal yeast Malassezia sympodialis and demonstrate how proteogenomics can substantially improve gene annotation. Through long-read DNA sequencing, we obtained a gap-free genome assembly for M. sympodialis (ATCC 42132), comprising eight nuclear and one mitochondrial chromosome. We also sequenced and assembled four M. sympodialis clinical isolates, and showed their value for understanding Malassezia reproduction by confirming four alternative allele combinations at the two mating-type loci. Importantly, we demonstrated how proteomics data could be readily integrated with transcriptomics data in standard annotation tools. This increased the number of annotated protein-coding genes by 14% (from 3612 to 4113), compared to using transcriptomics evidence alone. Manual curation further increased the number of protein-coding genes by 9% (to 4493). All of these genes have RNA-seq evidence and 87% were confirmed by proteomics. The M. sympodialis genome assembly and annotation presented here is at a quality yet achieved only for a few eukaryotic organisms, and constitutes an important reference for future host-microbe interaction studies.

Multiple Origins of the Pathogenic Yeast Candida orthopsilosis by Separate Hybridizations between Two Parental Species.

Mating between different species produces hybrids that are usually asexual and stuck as diploids, but can also lead to the formation of new species. Here, we report the genome sequences of 27 isolates of the pathogenic yeast Candida orthopsilosis. We find that most isolates are diploid hybrids, products of mating between two unknown parental species (A and B) that are 5% divergent in sequence. Isolates vary greatly in the extent of homogenization between A and B, making their genomes a mosaic of highly heterozygous regions interspersed with homozygous regions. Separate phylogenetic analyses of SNPs in the A- and B-derived portions of the genome produces almost identical trees of the isolates with four major clades. However, the presence of two mutually exclusive genotype combinations at the mating type locus, and recombinant mitochondrial genomes diagnostic of inter-clade mating, shows that the species C. orthopsilosis does not have a single evolutionary origin but was created at least four times by separate interspecies hybridizations between parents A and B. Older hybrids have lost more heterozygosity. We also identify two isolates with homozygous genomes derived exclusively from parent A, which are pure non-hybrid strains. The parallel emergence of the same hybrid species from multiple independent hybridization events is common in plant evolution, but is much less documented in pathogenic fungi.

Genus-Wide Comparative Genomics of Malassezia Delineates Its Phylogeny, Physiology, and Niche Adaptation on Human Skin.

Malassezia is a unique lipophilic genus in class Malasseziomycetes in Ustilaginomycotina, (Basidiomycota, fungi) that otherwise consists almost exclusively of plant pathogens. Malassezia are typically isolated from warm-blooded animals, are dominant members of the human skin mycobiome and are associated with common skin disorders. To characterize the genetic basis of the unique phenotypes of Malassezia spp., we sequenced the genomes of all 14 accepted species and used comparative genomics against a broad panel of fungal genomes to comprehensively identify distinct features that define the Malassezia gene repertoire: gene gain and loss; selection signatures; and lineage-specific gene family expansions. Our analysis revealed key gene gain events (64) with a single gene conserved across all Malassezia but absent in all other sequenced Basidiomycota. These likely horizontally transferred genes provide intriguing gain-of-function events and prime candidates to explain the emergence of Malassezia. A larger set of genes (741) were lost, with enrichment for glycosyl hydrolases and carbohydrate metabolism, concordant with adaptation to skin's carbohydrate-deficient environment. Gene family analysis revealed extensive turnover and underlined the importance of secretory lipases, phospholipases, aspartyl proteases, and other peptidases. Combining genomic analysis with a re-evaluation of culture characteristics, we establish the likely lipid-dependence of all Malassezia. Our phylogenetic analysis sheds new light on the relationship between Malassezia and other members of Ustilaginomycotina, as well as phylogenetic lineages within the genus. Overall, our study provides a unique genomic resource for understanding Malassezia niche-specificity and potential virulence, as well as their abundance and distribution in the environment and on human skin.

Genefu: an R/Bioconductor package for computation of gene expression-based signatures in breast cancer.

UNLABELLED: Breast cancer is one of the most frequent cancers among women. Extensive studies into the molecular heterogeneity of breast cancer have produced a plethora of molecular subtype classification and prognosis prediction algorithms, as well as numerous gene expression signatures. However, reimplementation of these algorithms is a tedious but important task to enable comparison of existing signatures and classification models between each other and with new models. Here, we present the genefu R/Bioconductor package, a multi-tiered compendium of bioinformatics algorithms and gene signatures for molecular subtyping and prognostication in breast cancer. AVAILABILITY AND IMPLEMENTATION: The genefu package is available from Bioconductor. http://www.bioconductor.org/packages/devel/bioc/html/genefu.html Source code is also available on Github https://github.com/bhklab/genefu CONTACT: bhaibeka@uhnresearch.ca SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Comparative phenotypic analysis of the major fungal pathogens Candida parapsilosis and Candida albicans.

Candida parapsilosis and Candida albicans are human fungal pathogens that belong to the CTG clade in the Saccharomycotina. In contrast to C. albicans, relatively little is known about the virulence properties of C. parapsilosis, a pathogen particularly associated with infections of premature neonates. We describe here the construction of C. parapsilosis strains carrying double allele deletions of 100 transcription factors, protein kinases and species-specific genes. Two independent deletions were constructed for each target gene. Growth in >40 conditions was tested, including carbon source, temperature, and the presence of antifungal drugs. The phenotypes were compared to C. albicans strains with deletions of orthologous transcription factors. We found that many phenotypes are shared between the two species, such as the role of Upc2 as a regulator of azole resistance, and of CAP1 in the oxidative stress response. Others are unique to one species. For example, Cph2 plays a role in the hypoxic response in C. parapsilosis but not in C. albicans. We found extensive divergence between the biofilm regulators of the two species. We identified seven transcription factors and one protein kinase that are required for biofilm development in C. parapsilosis. Only three (Efg1, Bcr1 and Ace2) have similar effects on C. albicans biofilms, whereas Cph2, Czf1, Gzf3 and Ume6 have major roles in C. parapsilosis only. Two transcription factors (Brg1 and Tec1) with well-characterized roles in biofilm formation in C. albicans do not have the same function in C. parapsilosis. We also compared the transcription profile of C. parapsilosis and C. albicans biofilms. Our analysis suggests the processes shared between the two species are predominantly metabolic, and that Cph2 and Bcr1 are major biofilm regulators in C. parapsilosis.

Comparative genome analysis and gene finding in Candida species using CGOB.

The Candida Gene Order Browser (CGOB) was developed as a tool to visualize and analyze synteny relationships in multiple Candida species, and to provide an accurate, manually curated set of orthologous Candida genes for evolutionary analyses. Here, we describe major improvements to CGOB. The underlying structure of the database has been changed significantly. Genomic features are now based directly on genome annotations rather than on protein sequences, which allows non-protein features such as centromere locations in Candida albicans and tRNA genes in all species to be included. The data set has been expanded to 13 species, including genomes of pathogens (C. albicans, C. parapsilosis, C. tropicalis, and C. orthopsilosis), and those of xylose-degrading species with important biotechnological applications (C. tenuis, Scheffersomyces stipitis, and Spathaspora passalidarum). Updated annotations of C. parapsilosis, C. dubliniensis, and Debaryomyces hansenii have been incorporated. We discovered more than 1,500 previously unannotated genes among the 13 genomes, ranging in size from 29 to 3,850 amino acids. Poorly conserved and rapidly evolving genes were also identified. Re-analysis of the mating type loci of the xylose degraders suggests that C. tenuis is heterothallic, whereas both Spa. passalidarum and S. stipitis are homothallic. As well as hosting the browser, the CGOB website (http://cgob.ucd.ie) gives direct access to all the underlying genome annotations, sequences, and curated orthology data.

RamiGO: an R/Bioconductor package providing an AmiGO visualize interface.

SUMMARY: The R/Bioconductor package RamiGO is an R interface to AmiGO that enables visualization of Gene Ontology (GO) trees. Given a list of GO terms, RamiGO uses the AmiGO visualize API to import Graphviz-DOT format files into R, and export these either as images (SVG, PNG) or into Cytoscape for extended network analyses. RamiGO provides easy customization of annotation, highlighting of specific GO terms, colouring of terms by P-value or export of a simplified summary GO tree. We illustrate RamiGO functionalities in a genome-wide gene set analysis of prognostic genes in breast cancer. AVAILABILITY AND IMPLEMENTATION: RamiGO is provided in R/Bioconductor, is open source under the Artistic-2.0 License and is available with a user manual containing installation, operating instructions and tutorials. It requires R version 2.15.0 or higher. URL: http://bioconductor.org/packages/release/bioc/html/RamiGO.html

GeneSigDB: a manually curated database and resource for analysis of gene expression signatures.

GeneSigDB (http://www.genesigdb.org or http://compbio.dfci.harvard.edu/genesigdb/) is a database of gene signatures that have been extracted and manually curated from the published literature. It provides a standardized resource of published prognostic, diagnostic and other gene signatures of cancer and related disease to the community so they can compare the predictive power of gene signatures or use these in gene set enrichment analysis. Since GeneSigDB release 1.0, we have expanded from 575 to 3515 gene signatures, which were collected and transcribed from 1604 published articles largely focused on gene expression in cancer, stem cells, immune cells, development and lung disease. We have made substantial upgrades to the GeneSigDB website to improve accessibility and usability, including adding a tag cloud browse function, facetted navigation and a 'basket' feature to store genes or gene signatures of interest. Users can analyze GeneSigDB gene signatures, or upload their own gene list, to identify gene signatures with significant gene overlap and results can be viewed on a dynamic editable heatmap that can be downloaded as a publication quality image. All data in GeneSigDB can be downloaded in numerous formats including .gmt file format for gene set enrichment analysis or as a R/Bioconductor data file. GeneSigDB is available from http://www.genesigdb.org.

Inhibition of nonsense-mediated mRNA decay reduces the tumorigenicity of human fibrosarcoma cells.

Nonsense-mediated mRNA decay (NMD) is a eukaryotic RNA decay pathway with roles in cellular stress responses, differentiation, and viral defense. It functions in both quality control and post-transcriptional regulation of gene expression. NMD has also emerged as a modulator of cancer progression, although available evidence supports both a tumor suppressor and a pro-tumorigenic role, depending on the model. To further investigate the role of NMD in cancer, we knocked out the NMD factor SMG7 in the HT1080 human fibrosarcoma cell line, resulting in suppression of NMD function. We then compared the oncogenic properties of the parental cell line, the SMG7-knockout, and a rescue cell line in which we re-introduced both isoforms of SMG7. We also tested the effect of a drug inhibiting the NMD factor SMG1 to distinguish NMD-dependent effects from putative NMD-independent functions of SMG7. Using cell-based assays and a mouse xenograft tumor model, we showed that suppression of NMD function severely compromises the oncogenic phenotype. Molecular pathway analysis revealed that NMD suppression strongly reduces matrix metalloprotease 9 (MMP9) expression and that MMP9 re-expression partially rescues the oncogenic phenotype. Since MMP9 promotes cancer cell migration and invasion, metastasis and angiogenesis, its downregulation may contribute to the reduced tumorigenicity of NMD-suppressed cells. Collectively, our results highlight the potential value of NMD inhibition as a therapeutic approach.

Epitope-engineered human hematopoietic stem cells are shielded from CD123-targeted immunotherapy.

Targeted eradication of transformed or otherwise dysregulated cells using monoclonal antibodies (mAb), antibody-drug conjugates (ADC), T cell engagers (TCE), or chimeric antigen receptor (CAR) cells is very effective for hematologic diseases. Unlike the breakthrough progress achieved for B cell malignancies, there is a pressing need to find suitable antigens for myeloid malignancies. CD123, the interleukin-3 (IL-3) receptor alpha-chain, is highly expressed in various hematological malignancies, including acute myeloid leukemia (AML). However, shared CD123 expression on healthy hematopoietic stem and progenitor cells (HSPCs) bears the risk for myelotoxicity. We demonstrate that epitope-engineered HSPCs were shielded from CD123-targeted immunotherapy but remained functional, while CD123-deficient HSPCs displayed a competitive disadvantage. Transplantation of genome-edited HSPCs could enable tumor-selective targeted immunotherapy while rebuilding a fully functional hematopoietic system. We envision that this approach is broadly applicable to other targets and cells, could render hitherto undruggable targets accessible to immunotherapy, and will allow continued posttransplant therapy, for instance, to treat minimal residual disease (MRD).

survcomp: an R/Bioconductor package for performance assessment and comparison of survival models.

SUMMARY: The survcomp package provides functions to assess and statistically compare the performance of survival/risk prediction models. It implements state-of-the-art statistics to (i) measure the performance of risk prediction models; (ii) combine these statistical estimates from multiple datasets using a meta-analytical framework; and (iii) statistically compare the performance of competitive models.

Recursive Editing improves homology-directed repair through retargeting of undesired outcomes.

CRISPR-Cas induced homology-directed repair (HDR) enables the installation of a broad range of precise genomic modifications from an exogenous donor template. However, applications of HDR in human cells are often hampered by poor efficiency, stemming from a preference for error-prone end joining pathways that yield short insertions and deletions. Here, we describe Recursive Editing, an HDR improvement strategy that selectively retargets undesired indel outcomes to create additional opportunities to produce the desired HDR allele. We introduce a software tool, named REtarget, that enables the rational design of Recursive Editing experiments. Using REtarget-designed guide RNAs in single editing reactions, Recursive Editing can simultaneously boost HDR efficiencies and reduce undesired indels. We also harness REtarget to generate databases for particularly effective Recursive Editing sites across the genome, to endogenously tag proteins, and to target pathogenic mutations. Recursive Editing constitutes an easy-to-use approach without potentially deleterious cell manipulations and little added experimental burden.

Controlled Cycling and Quiescence Enables Efficient HDR in Engraftment-Enriched Adult Hematopoietic Stem and Progenitor Cells.

Genome editing often takes the form of either error-prone sequence disruption by non-homologous end joining (NHEJ) or sequence replacement by homology-directed repair (HDR). Although NHEJ is generally effective, HDR is often difficult in primary cells. Here, we use a combination of immunophenotyping, next-generation sequencing, and single-cell RNA sequencing to investigate and reprogram genome editing outcomes in subpopulations of adult hematopoietic stem and progenitor cells. We find that although quiescent stem-enriched cells mostly use NHEJ, non-quiescent cells with the same immunophenotype use both NHEJ and HDR. Inducing quiescence before editing results in a loss of HDR in all cell subtypes. We develop a strategy of controlled cycling and quiescence that yields a 6-fold increase in the HDR/NHEJ ratio in quiescent stem cells ex vivo and in vivo. Our results highlight the tension between editing and cellular physiology and suggest strategies to manipulate quiescent cells for research and therapeutic genome editing.

ATF4 Regulates MYB to Increase γ-Globin in Response to Loss of ß-Globin.

ß-Hemoglobinopathies can trigger rapid production of red blood cells in a process known as stress erythropoiesis. Cellular stress prompts differentiating erythroid precursors to express high levels of fetal γ-globin. However, the mechanisms underlying γ-globin production during cellular stress are still poorly defined. Here, we use CRISPR-Cas genome editing to model the stress caused by reduced levels of adult ß-globin. We find that decreased ß-globin is sufficient to induce robust re-expression of γ-globin, and RNA sequencing (RNA-seq) of differentiating isogenic erythroid precursors implicates ATF4 as a causal regulator of this response. ATF4 binds within the HBS1L-MYB intergenic enhancer and regulates expression of MYB, a known γ-globin regulator. Overall, the reduction of ATF4 upon ß-globin knockout decreases the levels of MYB and BCL11A. Identification of ATF4 as a key regulator of globin compensation adds mechanistic insight to the poorly understood phenomenon of stress-induced globin compensation and could inform strategies to treat hemoglobinopathies.

Integrated genomic analysis reveals mutated ELF3 as a potential gallbladder cancer vaccine candidate.

Gallbladder cancer (GBC) is an aggressive gastrointestinal malignancy with no approved targeted therapy. Here, we analyze exomes (n = 160), transcriptomes (n = 115), and low pass whole genomes (n = 146) from 167 gallbladder cancers (GBCs) from patients in Korea, India and Chile. In addition, we also sequence samples from 39 GBC high-risk patients and detect evidence of early cancer-related genomic lesions. Among the several significantly mutated genes not previously linked to GBC are ETS domain genes ELF3 and EHF, CTNNB1, APC, NSD1, KAT8, STK11 and NFE2L2. A majority of ELF3 alterations are frame-shift mutations that result in several cancer-specific neoantigens that activate T-cells indicating that they are cancer vaccine candidates. In addition, we identify recurrent alterations in KEAP1/NFE2L2 and WNT pathway in GBC. Taken together, these define multiple targetable therapeutic interventions opportunities for GBC treatment and management.

The Indian cobra reference genome and transcriptome enables comprehensive identification of venom toxins.

Snakebite envenoming is a serious and neglected tropical disease that kills ~100,000 people annually. High-quality, genome-enabled comprehensive characterization of toxin genes will facilitate development of effective humanized recombinant antivenom. We report a de novo near-chromosomal genome assembly of Naja naja, the Indian cobra, a highly venomous, medically important snake. Our assembly has a scaffold N50 of 223.35 Mb, with 19 scaffolds containing 95% of the genome. Of the 23,248 predicted protein-coding genes, 12,346 venom-gland-expressed genes constitute the 'venom-ome' and this included 139 genes from 33 toxin families. Among the 139 toxin genes were 19 'venom-ome-specific toxins' (VSTs) that showed venom-gland-specific expression, and these probably encode the minimal core venom effector proteins. Synthetic venom reconstituted through recombinant VST expression will aid in the rapid development of safe and effective synthetic antivenom. Additionally, our genome could serve as a reference for snake genomes, support evolutionary studies and enable venom-driven drug discovery.