Characterising genome architectures using genome decomposition analysis.

Genome architecture describes how genes and other features are arranged in genomes. These arrangements reflect the evolutionary pressures on genomes and underlie biological processes such as chromosomal segregation and the regulation of gene expression. We present a new tool called Genome Decomposition Analysis (GDA) that characterises genome architectures and acts as an accessible approach for discovering hidden features of a genome assembly. With the imminent deluge of high-quality genome assemblies from projects such as the Darwin Tree of Life and the Earth BioGenome Project, GDA has been designed to facilitate their exploration and the discovery of novel genome biology. We highlight the effectiveness of our approach in characterising the genome architectures of single-celled eukaryotic parasites from the phylum Apicomplexa and show that it scales well to large genomes.

Antibody-independent mechanisms regulate the establishment of chronic Plasmodium infection.

Malaria is caused by parasites of the genus Plasmodium. All human-infecting Plasmodium species can establish long-lasting chronic infections1-5, creating an infectious reservoir to sustain transmission1,6. It is widely accepted that the maintenance of chronic infection involves evasion of adaptive immunity by antigenic variation7. However, genes involved in this process have been identified in only two of five human-infecting species: Plasmodium falciparum and Plasmodium knowlesi. Furthermore, little is understood about the early events in the establishment of chronic infection in these species. Using a rodent model we demonstrate that from the infecting population, only a minority of parasites, expressing one of several clusters of virulence-associated pir genes, establishes a chronic infection. This process occurs in different species of parasites and in different hosts. Establishment of chronicity is independent of adaptive immunity and therefore different from the mechanism proposed for maintenance of chronic P. falciparum infections7-9. Furthermore, we show that the proportions of parasites expressing different types of pir genes regulate the time taken to establish a chronic infection. Because pir genes are common to most, if not all, species of Plasmodium10, this process may be a common way of regulating the establishment of chronic infections.

Population, genetic, and antigenic diversity of the apicomplexan Eimeria tenella and their relevance to vaccine development.

The phylum Apicomplexa includes serious pathogens of humans and animals. Understanding the distribution and population structure of these protozoan parasites is of fundamental importance to explain disease epidemiology and develop sustainable controls. Predicting the likely efficacy and longevity of subunit vaccines in field populations relies on knowledge of relevant preexisting antigenic diversity, population structure, the likelihood of coinfection by genetically distinct strains, and the efficiency of cross-fertilization. All four of these factors have been investigated for Plasmodium species parasites, revealing both clonal and panmictic population structures with exceptional polymorphism associated with immunoprotective antigens such as apical membrane antigen 1 (AMA1). For the coccidian Toxoplasma gondii only genomic diversity and population structure have been defined in depth so far; for the closely related Eimeria species, all four variables are currently unknown. Using Eimeria tenella, a major cause of the enteric disease coccidiosis, which exerts a profound effect on chicken productivity and welfare, we determined population structure, genotype distribution, and likelihood of cross-fertilization during coinfection and also investigated the extent of naturally occurring antigenic diversity for the E. tenella AMA1 homolog. Using genome-wide Sequenom SNP-based haplotyping, targeted sequencing, and single-cell genotyping, we show that in this coccidian the functionality of EtAMA1 appears to outweigh immune evasion. This result is in direct contrast to the situation in Plasmodium and most likely is underpinned by the biology of the direct and acute coccidian life cycle in the definitive host.

Genes involved in host-parasite interactions can be revealed by their correlated expression.

Molecular interactions between a parasite and its host are key to the ability of the parasite to enter the host and persist. Our understanding of the genes and proteins involved in these interactions is limited. To better understand these processes it would be advantageous to have a range of methods to predict pairs of genes involved in such interactions. Correlated gene expression profiles can be used to identify molecular interactions within a species. Here we have extended the concept to different species, showing that genes with correlated expression are more likely to encode proteins, which directly or indirectly participate in host-parasite interaction. We go on to examine our predictions of molecular interactions between the malaria parasite and both its mammalian host and insect vector. Our approach could be applied to study any interaction between species, for example, between a host and its parasites or pathogens, but also symbiotic and commensal pairings.

Comparative genomics of the apicomplexan parasites Toxoplasma gondii and Neospora caninum: Coccidia differing in host range and transmission strategy.

Toxoplasma gondii is a zoonotic protozoan parasite which infects nearly one third of the human population and is found in an extraordinary range of vertebrate hosts. Its epidemiology depends heavily on horizontal transmission, especially between rodents and its definitive host, the cat. Neospora caninum is a recently discovered close relative of Toxoplasma, whose definitive host is the dog. Both species are tissue-dwelling Coccidia and members of the phylum Apicomplexa; they share many common features, but Neospora neither infects humans nor shares the same wide host range as Toxoplasma, rather it shows a striking preference for highly efficient vertical transmission in cattle. These species therefore provide a remarkable opportunity to investigate mechanisms of host restriction, transmission strategies, virulence and zoonotic potential. We sequenced the genome of N. caninum and transcriptomes of the invasive stage of both species, undertaking an extensive comparative genomics and transcriptomics analysis. We estimate that these organisms diverged from their common ancestor around 28 million years ago and find that both genomes and gene expression are remarkably conserved. However, in N. caninum we identified an unexpected expansion of surface antigen gene families and the divergence of secreted virulence factors, including rhoptry kinases. Specifically we show that the rhoptry kinase ROP18 is pseudogenised in N. caninum and that, as a possible consequence, Neospora is unable to phosphorylate host immunity-related GTPases, as Toxoplasma does. This defense strategy is thought to be key to virulence in Toxoplasma. We conclude that the ecological niches occupied by these species are influenced by a relatively small number of gene products which operate at the host-parasite interface and that the dominance of vertical transmission in N. caninum may be associated with the evolution of reduced virulence in this species.

Characterization and gene expression analysis of the cir multi-gene family of Plasmodium chabaudi chabaudi (AS).

BACKGROUND: The pir genes comprise the largest multi-gene family in Plasmodium, with members found in P. vivax, P. knowlesi and the rodent malaria species. Despite comprising up to 5% of the genome, little is known about the functions of the proteins encoded by pir genes. P. chabaudi causes chronic infection in mice, which may be due to antigenic variation. In this model, pir genes are called cirs and may be involved in this mechanism, allowing evasion of host immune responses. In order to fully understand the role(s) of CIR proteins during P. chabaudi infection, a detailed characterization of the cir gene family was required. RESULTS: The cir repertoire was annotated and a detailed bioinformatic characterization of the encoded CIR proteins was performed. Two major sub-families were identified, which have been named A and B. Members of each sub-family displayed different amino acid motifs, and were thus predicted to have undergone functional divergence. In addition, the expression of the entire cir repertoire was analyzed via RNA sequencing and microarray. Up to 40% of the cir gene repertoire was expressed in the parasite population during infection, and dominant cir transcripts could be identified. In addition, some differences were observed in the pattern of expression between the cir subgroups at the peak of P. chabaudi infection. Finally, specific cir genes were expressed at different time points during asexual blood stages. CONCLUSIONS: In conclusion, the large number of cir genes and their expression throughout the intraerythrocytic cycle of development indicates that CIR proteins are likely to be important for parasite survival. In particular, the detection of dominant cir transcripts at the peak of P. chabaudi infection supports the idea that CIR proteins are expressed, and could perform important functions in the biology of this parasite. Further application of the methodologies described here may allow the elucidation of CIR sub-family A and B protein functions, including their contribution to antigenic variation and immune evasion.

Methods of remote homology detection can be combined to increase coverage by 10% in the midnight zone.

MOTIVATION: A recent development in sequence-based remote homologue detection is the introduction of profile-profile comparison methods. These are more powerful than previous technologies and can detect potentially homologous relationships missed by structural classifications such as CATH and SCOP. As structural classifications traditionally act as the gold standard of homology this poses a challenge in benchmarking them. RESULTS: We present a novel approach which allows an accurate benchmark of these methods against the CATH structural classification. We then apply this approach to assess the accuracy of a range of publicly available methods for remote homology detection including several profile-profile methods (COMPASS, HHSearch, PRC) from two perspectives. First, in distinguishing homologous domains from non-homologues and second, in annotating proteomes with structural domain families. PRC is shown to be the best method for distinguishing homologues. We show that SAM is the best practical method for annotating genomes, whilst using COMPASS for the most remote homologues would increase coverage. Finally, we introduce a simple approach to increase the sensitivity of remote homologue detection by up to 10%. This is achieved by combining multiple methods with a jury vote. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.