Recent advances in high-throughput sequencing present a new opportunity to deeply probe an organism's transcriptome. In this study, we used Illumina-based massively parallel sequencing to gain new insight into the transcriptome (RNA-Seq) of the human malaria parasite, Plasmodium falciparum. Using data collected at seven time points during the intraerythrocytic developmental cycle, we (i) detect novel gene transcripts; (ii) correct hundreds of gene models; (iii) propose alternative splicing events; and (iv) predict 5' and 3' untranslated regions. Approximately 70% of the unique sequencing reads map to previously annotated protein-coding genes. The RNA-Seq results greatly improve existing annotation of the P. falciparum genome with over 10% of gene models modified. Our data confirm 75% of predicted splice sites and identify 202 new splice sites, including 84 previously uncharacterized alternative splicing events. We also discovered 107 novel transcripts and expression of 38 pseudogenes, with many demonstrating differential expression across the developmental time series. Our RNA-Seq results correlate well with DNA microarray analysis performed in parallel on the same samples, and provide improved resolution over the microarray-based method. These data reveal new features of the P. falciparum transcriptional landscape and significantly advance our understanding of the parasite's red blood cell-stage transcriptome.

Tsetse flies, vectors of African trypanosomes, undergo viviparous reproduction (the deposition of live offspring). This reproductive strategy results in a large maternal investment and the deposition of a small number of progeny during a female's lifespan. The reproductive biology of tsetse has been studied on a physiological level; however the molecular analysis of tsetse reproduction requires deeper investigation. To build a foundation from which to base molecular studies of tsetse reproduction, a cDNA library was generated from female tsetse (Glossina morsitans morsitans) reproductive tissues and the intrauterine developmental stages. 3438 expressed sequence tags were sequenced and analyzed.

Improved molecular diagnostic methods for detection drug resistance in Mycobacterium tuberculosis (MTB) strains are required. Resistance to first- and second- line anti-tuberculous drugs has been associated with single nucleotide polymorphisms (SNPs) in particular genes. However, these SNPs can vary between MTB lineages therefore local data is required to describe different strain populations. We used whole genome sequencing (WGS) to characterize 37 extensively drug-resistant (XDR) MTB isolates from Pakistan and investigated 40 genes associated with drug resistance. Rifampicin resistance was attributable to SNPs in the rpoB hot-spot region. Isoniazid resistance was most commonly associated with the katG codon 315 (92%) mutation followed by inhA S94A (8%) however, one strain did not have SNPs in katG, inhA or oxyR-ahpC. All strains were pyrazimamide resistant but only 43% had pncA SNPs. Ethambutol resistant strains predominantly had embB codon 306 (62%) mutations, but additional SNPs at embB codons 406, 378 and 328 were also present. Fluoroquinolone resistance was associated with gyrA 91-94 codons in 81% of strains; four strains had only gyrB mutations, while others did not have SNPs in either gyrA or gyrB. Streptomycin resistant strains had mutations in ribosomal RNA genes; rpsL codon 43 (42%); rrs 500 region (16%), and gidB (34%) while six strains did not have mutations in any of these genes. Amikacin/kanamycin/capreomycin resistance was associated with SNPs in rrs at nt1401 (78%) and nt1484 (3%), except in seven (19%) strains. We estimate that if only the common hot-spot region targets of current commercial assays were used, the concordance between phenotypic and genotypic testing for these XDR strains would vary between rifampicin (100%), isoniazid (92%), flouroquinolones (81%), aminoglycoside (78%) and ethambutol (62%); while pncA sequencing would provide genotypic resistance in less than half the isolates. This work highlights the importance of expanded targets for drug resistance detection in MTB isolates.

Many mosquito species are naturally polymorphic for their abilities to transmit parasites, a feature which is of great interest for controlling vector-borne disease. Aedes aegypti, the primary vector of dengue and yellow fever and a laboratory model for studying lymphatic filariasis, is genetically variable for its capacity to harbor the filarial nematode Brugia malayi. The genome of Ae. aegypti is large and repetitive, making genome resequencing difficult and expensive. We designed exome captures to target protein-coding regions of the genome, and used association mapping in a wild Kenyan population to identify a single, dominant, sex-linked locus underlying resistance. This falls in a region of the genome where a resistance locus was previously mapped in a line established in 1936, suggesting that this polymorphism has been maintained in the wild for the at least 80 years. We then crossed resistant and susceptible mosquitoes to place both alleles of the gene into a common genetic background, and used RNA-seq to measure the effect of this locus on gene expression. We found evidence for Toll, IMD, and JAK-STAT pathway activity in response to early stages of B. malayi infection when the parasites are beginning to die in the resistant genotype. We also found that resistant mosquitoes express anti-microbial peptides at the time of parasite-killing, and that this expression is suppressed in susceptible mosquitoes. Together, we have found that a single resistance locus leads to a higher immune response in resistant mosquitoes, and we identify genes in this region that may be responsible for this trait.

Mycobacteria possess different type VII secretion (T7S) systems to secrete proteins across their unusual cell envelope. One of these systems, ESX-5, is only present in slow-growing mycobacteria and responsible for the secretion of multiple substrates. However, the role of ESX-5 substrates in growth and/or virulence is largely unknown. In this study, we show that esx-5 is essential for growth of both Mycobacterium marinum and Mycobacterium bovis. Remarkably, this essentiality can be rescued by increasing the permeability of the outer membrane, either by altering its lipid composition or by the introduction of the heterologous porin MspA. Mutagenesis of the first nucleotide-binding domain of the membrane ATPase EccC5 prevented both ESX-5-dependent secretion and bacterial growth, but did not affect ESX-5 complex assembly. This suggests that the rescuing effect is not due to pores formed by the ESX-5 membrane complex, but caused by ESX-5 activity. Subsequent proteomic analysis to identify crucial ESX-5 substrates confirmed that all detectable PE and PPE proteins in the cell surface and cell envelope fractions were routed through ESX-5. Additionally, saturated transposon-directed insertion-site sequencing (TraDIS) was applied to both wild-type M. marinum cells and cells expressing mspA to identify genes that are not essential anymore in the presence of MspA. This analysis confirmed the importance of esx-5, but we could not identify essential ESX-5 substrates, indicating that multiple of these substrates are together responsible for the essentiality. Finally, examination of phenotypes on defined carbon sources revealed that an esx-5 mutant is strongly impaired in the uptake and utilization of hydrophobic carbon sources. Based on these data, we propose a model in which the ESX-5 system is responsible for the transport of cell envelope proteins that are required for nutrient uptake. These proteins might in this way compensate for the lack of MspA-like porins in slow-growing mycobacteria.

Mycobacterium tuberculosis drug resistance (DR) challenges effective tuberculosis disease control. Current molecular tests examine limited numbers of mutations, and although whole genome sequencing approaches could fully characterise DR, data complexity has restricted their clinical application. A library (1,325 mutations) predictive of DR for 15 anti-tuberculosis drugs was compiled and validated for 11 of them using genomic-phenotypic data from 792 strains. A rapid online 'TB-Profiler' tool was developed to report DR and strain-type profiles directly from raw sequences. Using our DR mutation library, in silico diagnostic accuracy was superior to some commercial diagnostics and alternative databases. The library will facilitate sequence-based drug-susceptibility testing.

The apicomplexan parasite Theileria annulata transforms infected host cells, inducing uncontrolled proliferation and clonal expansion of the parasitized cell population. Shortly after sporozoite entry into the target cell, the surrounding host cell membrane is dissolved and an array of host cell microtubules (MTs) surrounds the parasite, which develops into the transforming schizont. The latter does not egress to invade and transform other cells. Instead, it remains tethered to host cell MTs and, during mitosis and cytokinesis, engages the cell's astral and central spindle MTs to secure its distribution between the two daughter cells. The molecular mechanism by which the schizont recruits and stabilizes host cell MTs is not known. MT minus ends are mostly anchored in the MT organizing center, while the plus ends explore the cellular space, switching constantly between phases of growth and shrinkage (called dynamic instability). Assuming the plus ends of growing MTs provide the first point of contact with the parasite, we focused on the complex protein machinery associated with these structures. We now report how the schizont recruits end-binding protein 1 (EB1), a central component of the MT plus end protein interaction network and key regulator of host cell MT dynamics. Using a range of in vitro experiments, we demonstrate that T. annulata p104, a polymorphic antigen expressed on the schizont surface, functions as a genuine EB1-binding protein and can recruit EB1 in the absence of any other parasite proteins. Binding strictly depends on a consensus SxIP motif located in a highly disordered C-terminal region of p104. We further show that parasite interaction with host cell EB1 is cell cycle regulated. This is the first description of a pathogen-encoded protein to interact with EB1 via a bona-fide SxIP motif. Our findings provide important new insight into the mode of interaction between Theileria and the host cell cytoskeleton.

Cholera infection continues to be a threat to global public health. The current cholera pandemic associated with Vibrio cholerae El Tor has now been ongoing for over half a century.

Tuberculosis (TB) caused by Mycobacterium tuberculosis (Mtb) is the second major cause of death from an infectious disease worldwide. Recent advances in DNA sequencing are leading to the ability to generate whole genome information in clinical isolates of M. tuberculosis complex (MTBC). The identification of informative genetic variants such as phylogenetic markers and those associated with drug resistance or virulence will help barcode Mtb in the context of epidemiological, diagnostic and clinical studies. Mtb genomic datasets are increasingly available as raw sequences, which are potentially difficult and computer intensive to process, and compare across studies. Here we have processed the raw sequence data (>1500 isolates, eight studies) to compile a catalogue of SNPs (n = 74,039, 63% non-synonymous, 51.1% in more than one isolate, i.e. non-private), small indels (n = 4810) and larger structural variants (n = 800). We have developed the PolyTB web-based tool (http://pathogenseq.lshtm.ac.uk/polytb) to visualise the resulting variation and important meta-data (e.g. in silico inferred strain-types, location) within geographical map and phylogenetic views. This resource will allow researchers to identify polymorphisms within candidate genes of interest, as well as examine the genomic diversity and distribution of strains. PolyTB source code is freely available to researchers wishing to develop similar tools for their pathogen of interest.

Malaria is a major public health problem that is actively being addressed in a global eradication campaign. Increased population mobility through international air travel has elevated the risk of re-introducing parasites to elimination areas and dispersing drug-resistant parasites to new regions. A simple genetic marker that quickly and accurately identifies the geographic origin of infections would be a valuable public health tool for locating the source of imported outbreaks. Here we analyse the mitochondrion and apicoplast genomes of 711 Plasmodium falciparum isolates from 14 countries, and find evidence that they are non-recombining and co-inherited. The high degree of linkage produces a panel of relatively few single-nucleotide polymorphisms (SNPs) that is geographically informative. We design a 23-SNP barcode that is highly predictive (~92%) and easily adapted to aid case management in the field and survey parasite migration worldwide.

Conjugative plasmids have been identified in a wide variety of different bacteria, ranging from proteobacteria to firmicutes, and conjugation is one of the most efficient routes for horizontal gene transfer. The most widespread mechanism of plasmid conjugation relies on different variants of the type IV secretion pathway. Here, we describe the identification of a novel type of conjugative plasmid that seems to be unique for mycobacteria. Interestingly, while this plasmid is efficiently exchanged between different species of slow-growing mycobacteria, including Mycobacterium tuberculosis, it could not be transferred to any of the fast-growing mycobacteria tested. Genetic analysis of the conjugative plasmid showed the presence of a locus containing homologues of three type IV secretion system components and a relaxase. In addition, a new type VII secretion locus was present. Using transposon insertion mutagenesis, we show that in fact both these secretion systems are essential for conjugation, indicating that this plasmid represents a new class of conjugative plasmids requiring two secretion machineries. This plasmid could form a useful new tool to exchange or introduce DNA in slow-growing mycobacteria.

Plasmodium falciparum causes most human malaria deaths, having prehistorically evolved from parasites of African Great Apes. Here we explore the genomic basis of P. falciparum adaptation to human hosts by fully sequencing the genome of the closely related chimpanzee parasite species P. reichenowi, and obtaining partial sequence data from a more distantly related chimpanzee parasite (P. gaboni). The close relationship between P. reichenowi and P. falciparum is emphasized by almost complete conservation of genomic synteny, but against this strikingly conserved background we observe major differences at loci involved in erythrocyte invasion. The organization of most virulence-associated multigene families, including the hypervariable var genes, is broadly conserved, but P. falciparum has a smaller subset of rif and stevor genes whose products are expressed on the infected erythrocyte surface. Genome-wide analysis identifies other loci under recent positive selection, but a limited number of changes at the host-parasite interface may have mediated host switching.

Schistosomiasis is the most important helminthic disease of humanity in terms of morbidity and mortality. Facile manipulation of schistosomes using lentiviruses would enable advances in functional genomics in these and related neglected tropical diseases pathogens including tapeworms, and including their non-dividing cells. Such approaches have hitherto been unavailable. Blood stream forms of the human blood fluke, Schistosoma mansoni, the causative agent of the hepatointestinal schistosomiasis, were infected with the human HIV-1 isolate NL4-3 pseudotyped with vesicular stomatitis virus glycoprotein. The appearance of strong stop and positive strand cDNAs indicated that virions fused to schistosome cells, the nucleocapsid internalized and the RNA genome reverse transcribed. Anchored PCR analysis, sequencing HIV-1-specific anchored Illumina libraries and Whole Genome Sequencing (WGS) of schistosomes confirmed chromosomal integration; >8,000 integrations were mapped, distributed throughout the eight pairs of chromosomes including the sex chromosomes. The rate of integrations in the genome exceeded five per 1,000 kb and HIV-1 integrated into protein-encoding loci and elsewhere with integration bias dissimilar to that of human T cells. We estimated ~ 2,100 integrations per schistosomulum based on WGS, i.e. about two or three events per cell, comparable to integration rates in human cells. Accomplishment in schistosomes of post-entry processes essential for HIV-1replication, including integrase-catalyzed integration, was remarkable given the phylogenetic distance between schistosomes and primates, the natural hosts of the genus Lentivirus. These enigmatic findings revealed that HIV-1 was active within cells of S. mansoni, and provided the first demonstration that HIV-1 can integrate into the genome of an invertebrate.

Drug discovery pipelines for the "neglected diseases" are now heavily populated with nitroheterocyclic compounds. Recently, the bicyclic nitro-compounds (R)-PA-824, DNDI-VL-2098 and delamanid have been identified as potential candidates for the treatment of visceral leishmaniasis. Using a combination of quantitative proteomics and whole genome sequencing of susceptible and drug-resistant parasites we identified a putative NAD(P)H oxidase as the activating nitroreductase (NTR2). Whole genome sequencing revealed that deletion of a single cytosine in the gene for NTR2 that is likely to result in the expression of a non-functional truncated protein. Susceptibility of leishmania was restored by reintroduction of the wild-type gene into the resistant line, which was accompanied by the ability to metabolise these compounds. Overexpression of NTR2 in wild-type parasites rendered cells hyper-sensitive to bicyclic nitro-compounds, but only marginally to the monocyclic nitro-drugs, nifurtimox and fexinidazole sulfone, known to be activated by a mitochondrial oxygen-insensitive nitroreductase (NTR1). Conversely, a double knockout NTR2 null cell line was completely resistant to bicyclic nitro-compounds and only marginally resistant to nifurtimox. Sensitivity was fully restored on expression of NTR2 in the null background. Thus, NTR2 is necessary and sufficient for activation of these bicyclic nitro-drugs. Recombinant NTR2 was capable of reducing bicyclic nitro-compounds in the same rank order as drug sensitivity in vitro. These findings may aid the future development of better, novel anti-leishmanial drugs. Moreover, the discovery of anti-leishmanial nitro-drugs with independent modes of activation and independent mechanisms of resistance alleviates many of the concerns over the continued development of these compound series.

Onchocerciasis (river blindness) is a neglected tropical disease that has been successfully targeted by mass drug treatment programs in the Americas and small parts of Africa. Achieving the long-term goal of elimination of onchocerciasis, however, requires additional tools, including drugs, vaccines, and biomarkers of infection. Here, we describe the transcriptome and proteome profiles of the major vector and the human host stages (L1, L2, L3, molting L3, L4, adult male, and adult female) of Onchocerca volvulus along with the proteome of each parasitic stage and of its Wolbachia endosymbiont (wOv). In so doing, we have identified stage-specific pathways important to the parasite's adaptation to its human host during its early development. Further, we generated a protein array that, when screened with well-characterized human samples, identified novel diagnostic biomarkers of O. volvulus infection and new potential vaccine candidates. This immunomic approach not only demonstrates the power of this postgenomic discovery platform but also provides additional tools for onchocerciasis control programs.

Human onchocerciasis is a serious neglected tropical disease caused by the filarial nematode Onchocerca volvulus that can lead to blindness and chronic disability. Control of the disease relies largely on mass administration of a single drug, and the development of new drugs and vaccines depends on a better knowledge of parasite biology. Here, we describe the chromosomes of O. volvulus and its Wolbachia endosymbiont. We provide the highest-quality sequence assembly for any parasitic nematode to date, giving a glimpse into the evolution of filarial parasite chromosomes and proteomes. This resource was used to investigate gene families with key functions that could be potentially exploited as targets for future drugs. Using metabolic reconstruction of the nematode and its endosymbiont, we identified enzymes that are likely to be essential for O. volvulus viability. In addition, we have generated a list of proteins that could be targeted by Federal-Drug-Agency-approved but repurposed drugs, providing starting points for anti-onchocerciasis drug development.

Burkholderia cenocepacia infection often leads to fatal cepacia syndrome in cystic fibrosis patients. However, antibiotic therapy rarely results in complete eradication of the pathogen due to its intrinsic resistance to many clinically available antibiotics. Recent attention has turned to the identification of essential genes as the proteins encoded by these genes may serve as potential targets for development of novel antimicrobials. In this study, we utilized TraDIS (Transposon Directed Insertion-site Sequencing) as a genome-wide screening tool to facilitate the identification of B. cenocepacia genes essential for its growth and viability. A transposon mutant pool consisting of approximately 500,000 mutants was successfully constructed, with more than 400,000 unique transposon insertion sites identified by computational analysis of TraDIS datasets. The saturated library allowed for the identification of 383 genes that were predicted to be essential in B. cenocepacia. We extended the application of TraDIS to identify conditionally essential genes required for in vitro growth and revealed an additional repertoire of 439 genes to be crucial for B. cenocepacia growth under nutrient-depleted conditions. The library of B. cenocepacia mutants can subsequently be subjected to various biologically related conditions to facilitate the discovery of genes involved in niche adaptation as well as pathogenicity and virulence.

Methods to reliably assess the accuracy of genome sequence data are lacking. Currently completeness is only described qualitatively and mis-assemblies are overlooked. Here we present REAPR, a tool that precisely identifies errors in genome assemblies without the need for a reference sequence. We have validated REAPR on complete genomes or de novo assemblies from bacteria, malaria and Caenorhabditis elegans, and demonstrate that 86% and 82% of the human and mouse reference genomes are error-free, respectively. When applied to an ongoing genome project, REAPR provides corrected assembly statistics allowing the quantitative comparison of multiple assemblies. REAPR is available at http://www.sanger.ac.uk/resources/software/reapr/.

Advances in both high-throughput sequencing and whole-genome amplification (WGA) protocols have allowed genomes to be sequenced from femtograms of DNA, for example from individual cells or from precious clinical and archived samples. Using the highly curated Caenorhabditis elegans genome as a reference, we have sequenced and identified errors and biases associated with Illumina library construction, library insert size, different WGA methods and genome features such as GC bias and simple repeat content. Detailed analysis of the reads from amplified libraries revealed characteristics suggesting that majority of amplified fragment ends are identical but inverted versions of each other. Read coverage in amplified libraries is correlated with both tandem and inverted repeat content, while GC content only influences sequencing in long-insert libraries. Nevertheless, single nucleotide polymorphism (SNP) calls and assembly metrics from reads in amplified libraries show comparable results with unamplified libraries. To utilize the full potential of WGA to reveal the real biological interest, this article highlights the importance of recognizing additional sources of errors from amplified sequence reads and discusses the potential implications in downstream analyses.

Cell-cycle progression and cell division in eukaryotes are governed in part by the cyclin family and their regulation of cyclin-dependent kinases (CDKs). Cyclins are very well characterised in model systems such as yeast and human cells, but surprisingly little is known about their number and role in Plasmodium, the unicellular protozoan parasite that causes malaria. Malaria parasite cell division and proliferation differs from that of many eukaryotes. During its life cycle it undergoes two types of mitosis: endomitosis in asexual stages and an extremely rapid mitotic process during male gametogenesis. Both schizogony (producing merozoites) in host liver and red blood cells, and sporogony (producing sporozoites) in the mosquito vector, are endomitotic with repeated nuclear replication, without chromosome condensation, before cell division. The role of specific cyclins during Plasmodium cell proliferation was unknown. We show here that the Plasmodium genome contains only three cyclin genes, representing an unusual repertoire of cyclin classes. Expression and reverse genetic analyses of the single Plant (P)-type cyclin, CYC3, in the rodent malaria parasite, Plasmodium berghei, revealed a cytoplasmic and nuclear location of the GFP-tagged protein throughout the lifecycle. Deletion of cyc3 resulted in defects in size, number and growth of oocysts, with abnormalities in budding and sporozoite formation. Furthermore, global transcript analysis of the cyc3-deleted and wild type parasites at gametocyte and ookinete stages identified differentially expressed genes required for signalling, invasion and oocyst development. Collectively these data suggest that cyc3 modulates oocyst endomitotic development in Plasmodium berghei.