Kathryn Crouch

EuPathDB: the eukaryotic pathogen genomics database resource

The Eukaryotic Pathogen Genomics Database Resource (EuPathDB, http://eupathdb.org) is a collection of databases covering 170+ eukaryotic pathogens (protists & fungi), along with relevant free-living and non-pathogenic species, and select pathogen hosts. To facilitate the discovery of meaningful biological relationships, the databases couple preconfigured searches with visualization and analysis tools for comprehensive data mining via intuitive graphical interfaces and APIs. All data are analyzed with the same workflows, including creation of gene orthology profiles, so data are easily compared across data sets, data types and organisms. EuPathDB is updated with numerous new analysis tools, features, data sets and data types. New tools include GO, metabolic pathway and word enrichment analyses plus an online workspace for analysis of personal, non-public, large-scale data. Expanded data content is mostly genomic and functional genomic data while new data types include protein microarray, metabolic pathways, compounds, quantitative proteomics, copy number variation, and polysomal transcriptomics. New features include consistent categorization of searches, data sets and genome browser tracks; redesigned gene pages; effective integration of alternative transcripts; and a EuPathDB Galaxy instance for private analyses of a users data. Forthcoming upgrades include user workspaces for private integration of data with existing EuPathDB data and improved integration and presentation of host–pathogen interactions.

Characterization of the immunoglobulin repertoire of the spiny dogfish (Squalus acanthias).

The cartilaginous fish (chimeras, sharks, skates and rays) are the oldest group relative to mammals in which an adaptive immune system founded upon immunoglobulins has been found. In this manuscript we characterize the immunoglobulins of the spiny dogfish (Squalus acanthias) at both the molecular and expressed protein levels. Despite the presence of hundreds of IgM clusters in this species the serum levels of this isotype are comparatively low. However, analysis of cDNA sequences and serum protein suggests microheterogeneity in the IgM heavy chains and supports the proposal that different clusters are preferentially used in the two forms (monomer or pentamer) of this isotype. We also found that the IgNAR isotype in this species exists in a previously unknown multimeric format in serum. Finally, we identified a new form of the IgW isotype (the shark IgD orthologue), in which the leader is spliced directly to the first constant domain, resulting in a molecule lacking an antigen-binding domain.

TrypanoCyc : a community-led biochemical pathways database for Trypanosoma brucei

The metabolic network of a cell represents the catabolic and anabolic reactions that interconvert small molecules (metabolites) through the activity of enzymes, transporters and non-catalyzed chemical reactions. Our understanding of individual metabolic networks is increasing as we learn more about the enzymes that are active in particular cells under particular conditions and as technologies advance to allow detailed measurements of the cellular metabolome. Metabolic network databases are of increasing importance in allowing us to contextualise data sets emerging from transcriptomic, proteomic and metabolomic experiments. Here we present a dynamic database, TrypanoCyc (http://www.metexplore.fr/trypanocyc/), which describes the generic and condition-specific metabolic network of Trypanosoma brucei, a parasitic protozoan responsible for human and animal African trypanosomiasis. In addition to enabling navigation through the BioCyc-based TrypanoCyc interface, we have also implemented a network-based representation of the information through MetExplore, yielding a novel environment in which to visualise the metabolism of this important parasite.

FungiDB: An Integrated Bioinformatic Resource for Fungi and Oomycetes.

FungiDB (fungidb.org) is a free online resource for data mining and functional genomics analysis for fungal and oomycete species. FungiDB is part of the Eukaryotic Pathogen Genomics Database Resource (EuPathDB, eupathdb.org) platform that integrates genomic, transcriptomic, proteomic, and phenotypic datasets, and other types of data for pathogenic and nonpathogenic, free-living and parasitic organisms. FungiDB is one of the largest EuPathDB databases containing nearly 100 genomes obtained from GenBank, Aspergillus Genome Database (AspGD), The Broad Institute, Joint Genome Institute (JGI), Ensembl, and other sources. FungiDB offers a user-friendly web interface with embedded bioinformatics tools that support custom in silico experiments that leverage FungiDB-integrated data. In addition, a Galaxy-based workspace enables users to generate custom pipelines for large-scale data analysis (e.g., RNA-Seq, variant calling, etc.). This review provides an introduction to the FungiDB resources and focuses on available features, tools, and queries and how they can be used to mine data across a diverse range of integrated FungiDB datasets and records.

The ATR kinase of Trypanosoma brucei links DNA damage signalling and monoallelic control of surface antigen gene expression during antigenic variation

To evade mammalian immunity, Trypanosoma brucei switches the variant surface glycoprotein (VSG) expressed on its surface. Key to this reaction are controls exerted to ensure only one of many subtelomeric multigene VSG expression sites are transcribed at a time. DNA repair activities have to date been implicated only in catalysis of VSG switching by recombination, not transcriptional control. However, how VSG switching is signalled to guide the appropriate reaction, or to integrate switching into parasite growth, is unknown. Here we show that loss of ATR, a DNA damage signalling protein kinase, is lethal and causes increased nuclear genome lesions. ATR depletion also causes expression of mixed VSGs on the cell surface, increased transcription of genes from silent expression sites, and altered localisation of RNA Polymerase I and VEX1, factors involved in VSG transcription. The work therefore reveals that VSG expression control is mediated by a nuclear DNA damage signalling factor.

Validation of the protein kinase PfCLK3 as a multi-stage cross species malarial drug target

The requirement for next generation anti-malarials to be both curative and transmission blockers necessitate the identification of molecular pathways essential for viability of both asexual and sexual parasite life stages. Here we identify a selective inhibitor to the Plasmodium falciparum protein kinase PfCLK3 which we use in combination with chemogenetics, whole genome sequencing and transcriptomics to validate PfCLK3 as a druggable target acting at multiple parasite life stages. Consistent with the proposed role of PfCLK3 as a regulator of RNA splicing, inhibition results in the down-regulation of >400 genes essential for parasite survival. Through this mechanism, blocking PfCLK3 activity not only results in rapid killing of asexual blood stage parasites but is also effective on sporozoites and gametocytes as well as showing parasiticidal activity in all Plasmodium species tested. Hence, our data establishes PfCLK3 as a target with the potential to deliver both symptomatic treatment and transmission blocking in malaria.

RibonucleaseH1-targeted R-loops in surface antigen gene expression sites can direct trypanosome immune evasion

Switching of the Variant Surface Glycoprotein (VSG) in Trypanosoma brucei provides a crucial host immune evasion strategy that is catalysed both by transcription and recombination reactions, each operating within specialised telomeric VSG expression sites (ES). VSG switching is likely triggered by events focused on the single actively transcribed ES, from a repertoire of around 15, but the nature of such events is unclear. Here we show that RNA-DNA hybrids, called R-loops, form preferentially within sequences termed the 70 bp repeats in the actively transcribed ES, but spread throughout the active and inactive ES in the absence of RNase H1, which degrades R-loops. Loss of RNase H1 also leads to increased levels of VSG coat switching and replication-associated genome damage, some of which accumulates within the active ES. This work indicates VSG ES architecture elicits R-loop formation, and that these RNA-DNA hybrids connect T. brucei immune evasion by transcription and recombination. Author summary All pathogens must survive eradication by the host immune response in order to continue infections and be passed on to a new host. Changes in the proteins expressed on the surface of the pathogen, or on the surface of the cells the pathogen infects, is a widely used strategy to escape immune elimination. Understanding how this survival strategy, termed antigenic variation, operates in any pathogen is critical, both to understand interaction between the pathogen and host and disease progression. A key event in antigenic variation is the initiation of the change in expression of the surface protein gene, though how this occurs has been detailed in very few pathogens. Here we examine how changes in expression of the surface coat of the African trypanosome, which causes sleeping sickness disease, are initiated. We reveal that specialised nucleic acid structures, termed R-loops, form around the expressed trypanosome surface protein gene and increase in abundance after mutation of an enzyme that removes them, leading to increased changes in the surface coat in trypanosome cells that are dividing. We therefore shed light on the earliest acting events in trypanosome antigenic variation.

Genome-wide mapping reveals conserved and diverged R-loop activities in the unusual genetic landscape of the African trypanosome genome

Abstract R-loops are stable RNA–DNA hybrids that have been implicated in transcription initiation and termination, as well as in telomere maintenance, chromatin formation, and genome replication and instability. RNA Polymerase (Pol) II transcription in the protozoan parasite Trypanosoma brucei is highly unusual: virtually all genes are co-transcribed from multigene transcription units, with mRNAs generated by linked trans-splicing and polyadenylation, and transcription initiation sites display no conserved promoter motifs. Here, we describe the genome-wide distribution of R-loops in wild type mammal-infective T. brucei and in mutants lacking RNase H1, revealing both conserved and diverged functions. Conserved localization was found at centromeres, rRNA genes and retrotransposon-associated genes. RNA Pol II transcription initiation sites also displayed R-loops, suggesting a broadly conserved role despite the lack of promoter conservation or transcription initiation regulation. However, the most abundant sites of R-loop enrichment were within the regions between coding sequences of the multigene transcription units, where the hybrids coincide with sites of polyadenylation and nucleosome-depletion. Thus, instead of functioning in transcription termination the most widespread localization of R-loops in T. brucei suggests a novel correlation with pre-mRNA processing. Finally, we find little evidence for correlation between R-loop localization and mapped sites of DNA replication initiation.

A Leishmania infantum genetic marker associated with miltefosine treatment failure for visceral leishmaniasis

Background Miltefosine has been used successfully to treat visceral leishmaniasis (VL) in India, but it was unsuccessful for VL in a clinical trial in Brazil. Methods To identify molecular markers that predict VL treatment failure whole genome sequencing of 26 L. infantum isolates, from cured and relapsed patients allowed a GWAS analysis of SNPs, gene and chromosome copy number variations. Findings A strong association was identified (p = 0·0005) between the presence of a genetically stable L. infantumMiltefosine Sensitivity Locus (MSL), and a positive response to miltefosine treatment. The risk of treatment failure increased 9·4-fold (95% CI 2·11–53·54) when an isolate did not have the MSL. The complete absence of the MSL predicted miltefosine failure with 0·92 (95% CI 0·65–0·996) sensitivity and 0·78 (95% CI 0·52–0·92) specificity. A genotyping survey of L. infantum (n = 157) showed that the frequency of MSL varies in a cline from 95% in North East Brazil to <5% in the South East. The MSL was found in the genomes of all L. infantum and L. donovani sequenced isolates from the Old World (n = 671), where miltefosine can have a cure rate higher than 93%. Interpretation Knowledge on the presence or absence of the MSL in L. infantum will allow stratification of patients prior to treatment, helping to establish better therapeutic strategies for VL treatment. Fund CNPq, FAPES, GCRF MRC and Wellcome Trust.