The Apicomplexan whole-genome phylogeny: an analysis of incongruence among gene trees

The Tubulin Superfamily in Apicomplexan Parasites

Microtubules and specialized microtubule-containing structures are assembled from tubulins, an ancient superfamily of essential eukaryotic proteins. Here, we use bioinformatic approaches to analyze features of tubulins in organisms from the phylum Apicomplexa. Apicomplexans are protozoan parasites that cause a variety of human and animal infectious diseases. Individual species harbor one to four genes each for α- and β-tubulin isotypes. These may specify highly similar proteins, suggesting functional redundancy, or exhibit key differences, consistent with specialized roles. Some, but not all apicomplexans harbor genes for δ- and ε-tubulins, which are found in organisms that construct appendage-containing basal bodies. Critical roles for apicomplexan δ- and ε-tubulin are likely to be limited to microgametes, consistent with a restricted requirement for flagella in a single developmental stage. Sequence divergence or the loss of δ- and ε-tubulin genes in other apicomplexans appears to be associated with diminished requirements for centrioles, basal bodies, and axonemes. Finally, because spindle microtubules and flagellar structures have been proposed as targets for anti-parasitic therapies and transmission-blocking strategies, we discuss these ideas in the context of tubulin-based structures and tubulin superfamily properties.

Tukey's Depth for Object Data

We develop a novel exploratory tool for non-Euclidean object data based on data depth, extending celebrated Tukey's depth for Euclidean data. The proposed metric halfspace depth, applicable to data objects in a general metric space, assigns to data points depth values that characterize the centrality of these points with respect to the distribution and provides an interpretable center-outward ranking. Desirable theoretical properties that generalize standard depth properties postulated for Euclidean data are established for the metric halfspace depth. The depth median, defined as the deepest point, is shown to have high robustness as a location descriptor both in theory and in simulation. We propose an efficient algorithm to approximate the metric halfspace depth and illustrate its ability to adapt to the intrinsic data geometry. The metric halfspace depth was applied to an Alzheimer's disease study, revealing group differences in the brain connectivity, modeled as covariance matrices, for subjects in different stages of dementia. Based on phylogenetic trees of 7 pathogenic parasites, our proposed metric halfspace depth was also used to construct a meaningful consensus estimate of the evolutionary history and to identify potential outlier trees.

Tropical principal component analysis on the space of phylogenetic trees

MOTIVATION

Due to new technology for efficiently generating genome data, machine learning methods are urgently needed to analyze large sets of gene trees over the space of phylogenetic trees. However, the space of phylogenetic trees is not Euclidean, so ordinary machine learning methods cannot be directly applied. In 2019, Yoshida et al. introduced the notion of tropical principal component analysis (PCA), a statistical method for visualization and dimensionality reduction using a tropical polytope with a fixed number of vertices that minimizes the sum of tropical distances between each data point and its tropical projection. However, their work focused on the tropical projective space rather than the space of phylogenetic trees. We focus here on tropical PCA for dimension reduction and visualization over the space of phylogenetic trees.

RESULTS

Our main results are 2-fold: (i) theoretical interpretations of the tropical principal components over the space of phylogenetic trees, namely, the existence of a tropical cell decomposition into regions of fixed tree topology; and (ii) the development of a stochastic optimization method to estimate tropical PCs over the space of phylogenetic trees using a Markov Chain Monte Carlo approach. This method performs well with simulation studies, and it is applied to three empirical datasets: Apicomplexa and African coelacanth genomes as well as sequences of hemagglutinin for influenza from New York.

AVAILABILITY AND IMPLEMENTATION

Dataset: http://polytopes.net/Data.tar.gz. Code: http://polytopes.net/tropica_MCMC_codes.tar.gz.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Computational characterisation of Toxoplasma gondii FabG (3-oxoacyl-[acyl-carrier-protein] reductase): a combined virtual screening and all-atom molecular dynamics simulation study

Toxoplasma gondii is an opportunistic obligate parasite, ubiquitous around the globe with seropositivity rates that range from 10% to 90% and infection by the parasite of pregnant women causes pre-natal death of the foetus in most cases and severe neurodegenerative syndromes in some. No vaccine is currently available, and since drug-resistance is common among T. gondii strains, discovering lead compounds for drug design using diverse tactics is necessary. In this study, the sole constituent isoform of an enzymatic 3-oxoacyl-[acyl-carrier-protein] reduction step in an apicoplast-located fatty acid biosynthesis pathway was chosen as a possible drug target. FASII is prokaryotic therefore, targeting it would pose fewer side-effects to human hosts. After a homology 3D modelling of TgFabG, a high-throughput virtual screening of 9867 compounds, the elimination of ligands was carried out by a flexible ligand molecular docking and 200 ns molecular dynamics simulations, with additional DCCM and PC plot analyses. Molecular Dynamics and related post-MD analyses of the top 3 TgFabG binders selected for optimal free binding energies, showed that L2 maintained strong H-bonds with TgFabG and facilitated structural reorientation expected of FabGs, namely an expansion of the Rossmann Fold and a flexible lid capping. The most flexible TgFabG sites were the α7 helix (the flexible lid region) and the β4-α4 and β5-α6 loops. For TgFabG-L2, the movements of these regions toward the active site enabled greater ligand stability. Thus, L2 ("Skimmine"; PubChem ID: 320361), was ultimately selected as the optimal candidate for the discovery of lead compounds for rational drug design.Communicated by Ramaswamy H. Sarma.

imPhy: Imputing Phylogenetic Trees with Missing Information Using Mathematical Programming

Advances in modern genomics have allowed researchers to apply phylogenetic analyses on a genome-wide scale. While large volumes of genomic data can be generated cheaply and quickly, data missingness is a non-trivial and somewhat expected problem. Since the available information is often incomplete for a given set of genetic loci and individual organisms, a large proportion of trees that depict the evolutionary history of a single genetic locus, called gene trees, fail to contain all individuals. Data incompleteness causes difficulties in data collection, information extraction, and gene tree inference. Furthermore, identifying outlying gene trees, which can represent horizontal gene transfers, gene duplications, or hybridizations, is difficult when data is missing from the gene trees. The typical approach is to remove all individuals with missing data from the gene trees, and focus the analysis on individuals whose information is fully available - a huge loss of information. In this work, we propose and design an optimization-based imputation approach to infer the missing distances between leaves in a set of gene trees via a mixed integer non-linear programming model. We also present a new research pipeline, imPhy, that can (i) simulate a set of gene trees with leaves randomly missing in each tree, (ii) impute the missing pairwise distances in each gene tree, (iii) reconstruct the gene trees using the Neighbor Joining (NJ) and Unweighted Pair Group Method with Arithmetic Mean (UPGMA) methods, and (iv) analyze and report the efficiency of the reconstruction. To impute the missing leaves, we employ our newly proposed non-linear programming framework, and demonstrate its capability in reconstructing gene trees with incomplete information in both simulated and empirical datasets. In the empirical datasets apicomplexa and lungfish, our imputation has very small normalized mean square errors, even in the extreme case where 50 percent of the individuals in each gene tree are missing. Data, software, and user manuals can be found at https://github.com/yasuiniko/imPhy.

The tyrosine transporter of Toxoplasma gondii is a member of the newly defined apicomplexan amino acid transporter (ApiAT) family

Apicomplexan parasites are auxotrophic for a range of amino acids which must be salvaged from their host cells, either through direct uptake or degradation of host proteins. Here, we describe a family of plasma membrane-localized amino acid transporters, termed the Apicomplexan Amino acid Transporters (ApiATs), that are ubiquitous in apicomplexan parasites. Functional characterization of the ApiATs of Toxoplasma gondii indicate that several of these transporters are important for intracellular growth of the tachyzoite stage of the parasite, which is responsible for acute infections. We demonstrate that the ApiAT protein TgApiAT5-3 is an exchanger for aromatic and large neutral amino acids, with particular importance for L-tyrosine scavenging and amino acid homeostasis, and that TgApiAT5-3 is critical for parasite virulence. Our data indicate that T. gondii expresses additional proteins involved in the uptake of aromatic amino acids, and we present a model for the uptake and homeostasis of these amino acids. Our findings identify a family of amino acid transporters in apicomplexans, and highlight the importance of amino acid scavenging for the biology of this important phylum of intracellular parasites.

The Apicomplexa comprise a large number of parasitic protozoa that have obligate intracellular lifestyles and cause significant human and animal diseases, including malaria, cryptosporidiosis, toxoplasmosis, coccidiosis in poultry, and various cattle fevers. Apicomplexans must scavenge essential nutrients from their hosts in order to proliferate and cause disease, including a range of amino acids. The direct uptake of these nutrients is presumed to be mediated by transporter proteins located in the plasma membrane of intracellular stages, although the identities of these proteins are poorly defined. Using a combination of bioinformatic, genetic, cell biological, and physiological approaches, we have characterized an apicomplexan-specific family of plasma membrane-localized transporter proteins that we have called the Apicomplexan Amino acid Transporters (ApiATs). We show that TgApiAT5-3, a member of the family in the apicomplexan Toxoplasma gondii, is an exchanger for aromatic and large neutral amino acids. In particular, it is critical for uptake of tyrosine, and for parasite virulence in a mouse infection model. We conclude that ApiATs are a family of plasma membrane transporters that play crucial roles in amino acid scavenging by apicomplexan parasites.

Tropical Principal Component Analysis and Its Application to Phylogenetics

Principal component analysis is a widely used method for the dimensionality reduction of a given data set in a high-dimensional Euclidean space. Here we define and analyze two analogues of principal component analysis in the setting of tropical geometry. In one approach, we study the Stiefel tropical linear space of fixed dimension closest to the data points in the tropical projective torus; in the other approach, we consider the tropical polytope with a fixed number of vertices closest to the data points. We then give approximative algorithms for both approaches and apply them to phylogenetics, testing the methods on simulated phylogenetic data and on an empirical dataset of Apicomplexa genomes.

A Gene-Based Positive Selection Detection Approach to Identify Vaccine Candidates Using Toxoplasma gondii as a Test Case Protozoan Pathogen

Over the last two decades, various in silico approaches have been developed and refined that attempt to identify protein and/or peptide vaccines candidates from informative signals encoded in protein sequences of a target pathogen. As to date, no signal has been identified that clearly indicates a protein will effectively contribute to a protective immune response in a host. The premise for this study is that proteins under positive selection from the immune system are more likely suitable vaccine candidates than proteins exposed to other selection pressures. Furthermore, our expectation is that protein sequence regions encoding major histocompatibility complexes (MHC) binding peptides will contain consecutive positive selection sites. Using freely available data and bioinformatic tools, we present a high-throughput approach through a pipeline that predicts positive selection sites, protein subcellular locations, and sequence locations of medium to high T-Cell MHC class I binding peptides. Positive selection sites are estimated from a sequence alignment by comparing rates of synonymous (dS) and non-synonymous (dN) substitutions among protein coding sequences of orthologous genes in a phylogeny. The main pipeline output is a list of protein vaccine candidates predicted to be naturally exposed to the immune system and containing sites under positive selection. Candidates are ranked with respect to the number of consecutive sites located on protein sequence regions encoding MHCI-binding peptides. Results are constrained by the reliability of prediction programs and quality of input data. Protein sequences from Toxoplasma gondii ME49 strain (TGME49) were used as a case study. Surface antigen (SAG), dense granules (GRA), microneme (MIC), and rhoptry (ROP) proteins are considered worthy T. gondii candidates. Given 8263 TGME49 protein sequences processed anonymously, the top 10 predicted candidates were all worthy candidates. In particular, the top ten included ROP5 and ROP18, which are T. gondii virulence determinants. The chance of randomly selecting a ROP protein was 0.2% given 8263 sequences. We conclude that the approach described is a valuable addition to other in silico approaches to identify vaccines candidates worthy of laboratory validation and could be adapted for other apicomplexan parasite species (with appropriate data).

The Apicomplexa-specific glucosamine-6-phosphate N-acetyltransferase gene family encodes a key enzyme for glycoconjugate synthesis with potential as therapeutic target

Apicomplexa form a phylum of obligate parasitic protozoa of great clinical and veterinary importance. These parasites synthesize glycoconjugates for their survival and infectivity, but the enzymatic steps required to generate the glycosylation precursors are not completely characterized. In particular, glucosamine-phosphate N-acetyltransferase (GNA1) activity, needed to produce the essential UDP-N-acetylglucosamine (UDP-GlcNAc) donor, has not been identified in any Apicomplexa. We scanned the genomes of Plasmodium falciparum and representatives from six additional main lineages of the phylum for proteins containing the Gcn5-related N-acetyltransferase (GNAT) domain. One family of GNAT-domain containing proteins, composed by a P. falciparum sequence and its six apicomplexan orthologs, rescued the growth of a yeast temperature-sensitive GNA1 mutant. Heterologous expression and in vitro assays confirmed the GNA1 enzymatic activity in all lineages. Sequence, phylogenetic and synteny analyses suggest an independent origin of the Apicomplexa-specific GNA1 family, parallel to the evolution of a different GNA1 family in other eukaryotes. The inability to disrupt an otherwise modifiable gene target suggests that the enzyme is essential for P. falciparum growth. The relevance of UDP-GlcNAc for parasite viability, together with the independent evolution and unique sequence features of Apicomplexa GNA1, highlights the potential of this enzyme as a selective therapeutic target against apicomplexans.

PhyloChromoMap, a Tool for Mapping Phylogenomic History along Chromosomes, Reveals the Dynamic Nature of Karyotype Evolution in Plasmodium falciparum

The genome of Plasmodium falciparum, the causative agent of malaria in Africa, has been extensively studied since it was first fully sequenced in 2002. However, many open questions remain, including understanding the chromosomal context of molecular evolutionary changes (e.g., relationship between chromosome map and phylogenetic conservation, patterns of gene duplication, and patterns of selection). Here, we present PhyloChromoMap, a method that generates a phylogenomic map of chromosomes from a custom-built bioinformatics pipeline. Using P. falciparum 3D7 as a model, we analyze 2,116 genes with homologs in up to 941 diverse eukaryotic, bacterial and archaeal lineages. We estimate the level of conservation along chromosomes based on conservation across clades, and identify “young” regions (i.e., those with recent or fast evolving genes) that are enriched in subtelomeric regions as compared with internal regions. We also demonstrate that patterns of molecular evolution for paralogous genes differ significantly depending on their location as younger paralogs tend to be found in subtelomeric regions whereas older paralogs are enriched in internal regions. Combining these observations with analyses of synteny, we demonstrate that subtelomeric regions are actively shuffled among chromosome ends, which is consistent with the hypothesis that these regions are prone to ectopic recombination. We also assess patterns of selection by comparing dN/dS ratios of gene family members in subtelomeric versus internal regions, and we include the important antigenic gene family var. These analyses illustrate the highly dynamic nature of the karyotype of P. falciparum, and provide a method for exploring genome dynamics in other lineages.

Comparative analysis of low complexity regions in Plasmodia

Low complexity regions (LCRs) are a common feature shared by many genomes, but their evolutionary and functional significance remains mostly unknown. At the core of the uncertainty is a poor understanding of the mechanisms that regulate their retention in genomes, whether driven by natural selection or neutral evolution. Applying a comparative approach of LCRs to multiple strains and species is a powerful approach to identify patterns of conservation in these regions. Using this method, we investigate the evolutionary history of LCRs in the genus Plasmodium based on orthologous protein coding genes shared by 11 species and strains from primate and rodent-infecting pathogens. We find multiple lines of evidence in support of natural selection as a major evolutionary force shaping the composition and conservation of LCRs through time and signatures that their evolutionary paths are species specific. Our findings add a comparative analysis perspective to the debate on the evolution of LCRs and harness the power of sequence comparisons to identify potential functionally important LCR candidates.

Maximizing Power in Phylogenetics and Phylogenomics: A Perspective Illuminated by Fungal Big Data

Since its original inception over 150 years ago by Darwin, we have made tremendous progress toward the reconstruction of the Tree of Life. In particular, the transition from analyzing datasets comprised of small numbers of loci to those comprised of hundreds of loci, if not entire genomes, has aided in resolving some of the most vexing of evolutionary problems while giving us a new perspective on biodiversity. Correspondingly, phylogenetic trees have taken a central role in fields that span ecology, conservation, and medicine. However, the rise of big data has also presented phylogenomicists with a new set of challenges to experimental design, quantitative analyses, and computation. The sequencing of a number of very first genomes presented significant challenges to phylogenetic inference, leading fungal phylogenomicists to begin addressing pitfalls and postulating solutions to the issues that arise from genome-scale analyses relevant to any lineage across the Tree of Life. Here we highlight insights from fungal phylogenomics for topics including systematics and species delimitation, ecological and phenotypic diversification, and biogeography while providing an overview of progress made on the reconstruction of the fungal Tree of Life. Finally, we provide a review of considerations to phylogenomic experimental design for robust tree inference. We hope that this special issue of Advances in Genetics not only excites the continued progress of fungal evolutionary biology but also motivates the interdisciplinary development of new theory and methods designed to maximize the power of genomic scale data in phylogenetic analyses.

Normalizing Kernels in the Billera-Holmes-Vogtmann Treespace

As costs of genome sequencing have dropped precipitously, development of efficient bioinformatic methods to analyze genome structure and evolution have become ever more urgent. For example, most published phylogenomic studies involve either massive concatenation of sequences, or informal comparisons of phylogenies inferred on a small subset of orthologous genes, neither of which provides a comprehensive overview of evolution or systematic identification of genes with unusual and interesting evolution (e.g., horizontal gene transfers, gene duplication, and subsequent neofunctionalization). We are interested in identifying such "outlying" gene trees from the set of gene trees and estimating the distribution of trees over the "tree space". This paper describes an improvement to the kdetrees algorithm, an adaptation of classical kernel density estimation to the metric space of phylogenetic trees (Billera-Holmes-Vogtman treespace), whereby the kernel normalizing constants, are estimated through the use of the novel holonomic gradient methods. As in the original kdetrees paper, we have applied kdetrees to a set of Apicomplexa genes. The analysis identified several unreliable sequence alignments that had escaped previous detection, as well as a gene independently reported as a possible case of horizontal gene transfer. The updated version of the kdetrees software package is available both from CRAN (the official R package system), as well as from the official development repository on Github. ( github.com/grady/kdetrees).

Principal component analysis and the locus of the Fréchet mean in the space of phylogenetic trees

Evolutionary relationships are represented by phylogenetic trees, and a phylogenetic analysis of gene sequences typically produces a collection of these trees, one for each gene in the analysis. Analysis of samples of trees is difficult due to the multi-dimensionality of the space of possible trees. In Euclidean spaces, principal component analysis is a popular method of reducing high-dimensional data to a low-dimensional representation that preserves much of the sample’s structure. However, the space of all phylogenetic trees on a fixed set of species does not form a Euclidean vector space, and methods adapted to tree space are needed. Previous work introduced the notion of a principal geodesic in this space, analogous to the first principal component. Here we propose a geometric object for tree space similar to the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$k$\end{document}th principal component in Euclidean space: the locus of the weighted Fréchet mean of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$k+1$\end{document} vertex trees when the weights vary over the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$k$\end{document}-simplex. We establish some basic properties of these objects, in particular showing that they have dimension \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$k$\end{document}, and propose algorithms for projection onto these surfaces and for finding the principal locus associated with a sample of trees. Simulation studies demonstrate that these algorithms perform well, and analyses of two datasets, containing Apicomplexa and African coelacanth genomes respectively, reveal important structure from the second principal components.

Comparative transcriptomics of female and male gametocytes in Plasmodium berghei and the evolution of sex in alveolates

BACKGROUND

The clinical symptoms of malaria are caused by the asexual replication of Plasmodium parasites in the blood of the vertebrate host. To spread to new hosts, however, the malaria parasite must differentiate into sexual forms, termed gametocytes, which are ingested by a mosquito vector. Sexual differentiation produces either female or male gametocytes, and involves significant morphological and biochemical changes. These transformations prepare gametocytes for the rapid progression to gamete formation and fertilisation, which occur within 20 min of ingestion. Here we present the transcriptomes of asexual, female, and male gametocytes in P. berghei, and a comprehensive statistically-based differential-expression analysis of the transcriptional changes that underpin this sexual differentiation.

RESULTS

RNA-seq analysis revealed numerous differences in the transcriptomes of female and male gametocytes compared to asexual stages. Overall, there is net downregulation of transcripts in gametocytes compared to asexual stages, with this trend more marked in female gametocytes. Our analysis identified transcriptional changes in previously-characterised gametocyte-specific pathways, which validated our approach. We also detected many previously-unreported female- and male-specific pathways and genes. Transcriptional biases in stage and gender were then used to investigate sex-specificity and sexual dimorphism of Plasmodium in an evolutionary context. Sex-related gene expression is well conserved between Plasmodium species, but relatively poorly conserved in related organisms outside this genus. This pattern of conservation is most evident in genes necessary for both male and female gametocyte formation. However, this trend is less pronounced for male-specific genes, which are more highly conserved outside the genus than genes specific to female development.

CONCLUSIONS

We characterised the transcriptional changes that are integral to the development of the female and male sexual forms of Plasmodium. These differential-expression patterns provide a vital insight into understanding the gender-specific characteristics of this essential stage that is the primary target for treatments that block parasite transmission. Our results also offer insight into the evolution of sex genes through Alveolata, and suggest that many Plasmodium sex genes evolved within the genus. We further hypothesise that male gametocytes co-opted pre-existing cellular machinery in their evolutionary history, whereas female gametocytes evolved more through the development of novel, parasite-specific pathways.

Assessment and optimization of Theileria parva sporozoite full-length p67 antigen expression in mammalian cells

Delivery of various forms of recombinant Theileria parva sporozoite antigen (p67) has been shown to elicit antibody responses in cattle capable of providing protection against East Coast fever, the clinical disease caused by T. parva. Previous formulations of full-length and shorter recombinant versions of p67 derived from bacteria, insect, and mammalian cell systems are expressed in non-native and highly unstable forms. The stable expression of full-length recombinant p67 in mammalian cells has never been described and has remained especially elusive. In this study, p67 was expressed in human-derived cells as a full-length, membrane-linked protein and as a secreted form by omission of the putative transmembrane domain. The recombinant protein expressed in this system yielded primarily two products based on Western immunoblot analysis, including one at the expected size of 67 kDa, and one with a higher than expected molecular weight. Through treatment with PNGase F, our data indicate that the larger product of this mammalian cell-expressed recombinant p67 cannot be attributed to glycosylation. By increasing the denaturing conditions, we determined that the larger sized mammalian cell-expressed recombinant p67 product is likely a dimeric aggregate of the protein. Both forms of this recombinant p67 reacted with a monoclonal antibody to the p67 molecule, which reacts with the native sporozoite. Additionally, through this work we developed multiple mammalian cell lines, including both human and bovine-derived cell lines, transduced by a lentiviral vector, that are constitutively able to express a stable, secreted form of p67 for use in immunization, diagnostics, or in vitro assays. The recombinant p67 developed in this system is immunogenic in goats and cattle based on ELISA and flow cytometric analysis. The development of a mammalian cell system that expresses full-length p67 in a stable form as described here is expected to optimize p67-based immunization.

East Coast fever, caused by the tick-borne protozoan parasite Theileria parva, is a disease that results in significant bovine morbidity, mortality, and production losses in regions of sub-Saharan Africa. Susceptible cattle develop clinical signs within a 7–14 days of exposure, which often progress to severe pulmonary edema and death. Control of East Coast fever in affected regions of Africa is largely prohibited by the lack of an affordable and efficacious vaccine. Furthermore, pastoralist farmers in affected regions of Africa often lack resources to prevent losses due to East Coast fever, so these production losses play a significant role in food security and protein availability. Experimental immunization of cattle with a recombinant T. parva-derived antigen, p67, has shown promise in preventing East Coast fever, but this antigen is extremely difficult to produce in full-length in sufficient quantities, and results of immunization studies using truncated recombinant p67 products are highly inconsistent. In this study, p67 antigen production was optimized and produced for use in future immunization studies. Optimization of p67-based immunization strategies is an important step forward in the development of a sustainable, next-generation vaccine against T. parva, which is urgently needed to minimize losses associated with East Coast fever.

Genome microsatellite diversity within the Apicomplexa phylum

The Apicomplexa phylum groups include unicellular and obligate intracellular protozoan parasites with an apical complex used for attachment and invasion to host cells. In this study, we analyze single sequence repeats (SSRs) in the whole genome of 20 apicomplexan organisms that represent four different lineages within the phylum. Only perfect SSRs with at least 12 nucleotides and composed of 2-6 mers were included. To better understand the association of SSR types with the genomic regions, the SSRs were classified accordingly with the genomic location into exon, intron and intergenic categories. Our results showed heterogeneous SSRs density within the studied genomes. However, the most frequent SSRs types were di- and tri-nucleotide repeats. The former was associated with intergenic regions, while the latter was associated with exon regions.

Origin of a major infectious disease in vertebrates: The timing of Cryptosporidium evolution and its hosts

Protozoan parasites of the genus Cryptosporidium infect all vertebrate groups and display some host specificity in their infections. It is therefore possible to assume that Cryptosporidium parasites evolved intimately aside with vertebrate lineages. Here we propose a scenario of Cryptosporidium–Vertebrata coevolution testing the hypothesis that the origin of Cryptosporidium parasites follows that of the origin of modern vertebrates. We use calibrated molecular clocks and cophylogeny analyses to provide and compare age estimates and patterns of association between these clades. Our study provides strong support for the evolution of parasitism of Cryptosporidium with the rise of the vertebrates about 600 million years ago (Mya). Interestingly, periods of increased diversification in Cryptosporidium coincides with diversification of crown mammalian and avian orders after the Cretaceous-Palaeogene (K-Pg) boundary, suggesting that adaptive radiation to new mammalian and avian hosts triggered the diversification of this parasite lineage. Despite evidence for ongoing host shifts we also found significant correlation between protozoan parasites and vertebrate hosts trees in the cophylogenetic analysis. These results help us to understand the underlying macroevolutionary mechanisms driving evolution in Cryptosporidium and may have important implications for the ecology, dynamics and epidemiology of cryptosporidiosis disease in humans and other animals.

Extensive Shared Chemosensitivity between Malaria and Babesiosis Blood-Stage Parasites

The apicomplexan parasites that cause malaria and babesiosis invade and proliferate within erythrocytes. To assess the potential for common antiparasitic treatments, we measured the sensitivities of multiple species of Plasmodium and Babesia parasites to the chemically diverse collection of antimalarial compounds in the Malaria Box library. We observed that these parasites share sensitivities to a large fraction of the same inhibitors and we identified compounds with strong babesiacidal activity.

The Genome of Haemoproteus tartakovskyi and Its Relationship to Human Malaria Parasites

The phylogenetic relationships among hemosporidian parasites, including the origin of Plasmodium falciparum, the most virulent malaria parasite of humans, have been heavily debated for decades. Studies based on multiple-gene sequences have helped settle many of these controversial phylogenetic issues. However, denser taxon sampling and genome-wide analyses are needed to confidently resolve the evolutionay relationships among hemosporidian parasites. Genome sequences of several Plasmodium parasites are available but only for species infecting primates and rodents. To root the phylogenetic tree of Plasmodium, genomic data from related parasites of birds or reptiles are required. Here, we use a novel approach to isolate parasite DNA from microgametes and describe the first genome of a bird parasite in the sister genus to Plasmodium, Haemoproteus tartakovskyi. Similar to Plasmodium parasites, H. tartakovskyi has a small genome (23.2 Mb, 5,990 genes) and a GC content (25.4%) closer to P. falciparum (19.3%) than to Plasmodium vivax (42.3%). Combined with novel transcriptome sequences of the bird parasite Plasmodium ashfordi, our phylogenomic analyses of 1,302 orthologous genes demonstrate that mammalian-infecting malaria parasites are monophyletic, thus rejecting the repeatedly proposed hypothesis that the ancestor of Laverania parasites originated from a secondary host shift from birds to humans. Genes and genomic features previously found to be shared between P. falciparum and bird malaria parasites, but absent in other mammal malaria parasites, are therefore signatures of maintained ancestral states. We foresee that the genome of H. tartakovskyi will open new directions for comparative evolutionary analyses of malarial adaptive traits.

Profiles of low complexity regions in Apicomplexa

BACKGROUND

Low complexity regions (LCRs) are a ubiquitous feature in genomes and yet their evolutionary history and functional roles are unclear. Previous studies have shown contrasting evidence in favor of both neutral and selective mechanisms of evolution for different sets of LCRs suggesting that modes of identification of these regions may play a role in our ability to discern their evolutionary history. To further investigate this issue, we used a multiple threshold approach to identify species-specific profiles of proteome complexity and, by comparing properties of these sets, determine the influence that starting parameters have on evolutionary inferences.

RESULTS

We find that, although qualitatively similar, quantitatively each species has a unique LCR profile which represents the frequency of these regions within each genome. Inferences based on these profiles are more accurate in comparative analyses of genome complexity as they allow to determine the relative complexity of multiple genomes as well as the type of repetitiveness that is most common in each. Based on the multiple threshold LCR sets obtained, we identified predominant evolutionary mechanisms at different complexity levels, which show neutral mechanisms acting on highly repetitive LCRs (e.g., homopolymers) and selective forces becoming more important as heterogeneity of the LCRs increases.

CONCLUSIONS

Our results show how inferences based on LCRs are influenced by the parameters used to identify these regions. Sets of LCRs are heterogeneous aggregates of regions that include homo- and heteropolymers and, as such, evolve according to different mechanisms. LCR profiles provide a new way to investigate genome complexity across species and to determine the driving mechanism of their evolution.

An ensemble of specifically targeted proteins stabilizes cortical microtubules in the human parasite Toxoplasma gondii

The human parasite Toxoplasma gondii has 22 regularly spaced microtubules associated with the cortex. This work defines the differential localization of associated proteins, explores the biophysical constraints on specific targeting along the cortical microtubules, and investigates the function of these proteins in stabilizing the polymers.

Although all microtubules within a single cell are polymerized from virtually identical subunits, different microtubule populations carry out specialized and diverse functions, including directional transport, force generation, and cellular morphogenesis. Functional differentiation requires specific targeting of associated proteins to subsets or even subregions of these polymers. The cytoskeleton of Toxoplasma gondii, an important human parasite, contains at least five distinct tubulin-based structures. In this work, we define the differential localization of proteins along the cortical microtubules of T. gondii, established during daughter biogenesis and regulated by protein expression and exchange. These proteins distinguish cortical from mitotic spindle microtubules, even though the assembly of these subsets is contemporaneous during cell division. Finally, proteins associated with cortical microtubules collectively protect the stability of the polymers with a remarkable degree of functional redundancy.

A maximum pseudo-likelihood approach for phylogenetic networks

Background

Several phylogenomic analyses have recently demonstrated the need to account simultaneously for incomplete lineage sorting (ILS) and hybridization when inferring a species phylogeny. A maximum likelihood approach was introduced recently for inferring species phylogenies in the presence of both processes, and showed very good results. However, computing the likelihood of a model in this case is computationally infeasible except for very small data sets.

Results

Inspired by recent work on the pseudo-likelihood of species trees based on rooted triples, we introduce the pseudo-likelihood of a phylogenetic network, which, when combined with a search heuristic, provides a statistical method for phylogenetic network inference in the presence of ILS. Unlike trees, networks are not always uniquely encoded by a set of rooted triples. Therefore, even when given sufficient data, the method might converge to a network that is equivalent under rooted triples to the true one, but not the true one itself. The method is computationally efficient and has produced very good results on the data sets we analyzed. The method is implemented in PhyloNet, which is publicly available in open source.

Conclusions

Maximum pseudo-likelihood allows for inferring species phylogenies in the presence of hybridization and ILS, while scaling to much larger data sets than is currently feasible under full maximum likelihood. The nonuniqueness of phylogenetic networks encoded by a system of rooted triples notwithstanding, the proposed method infers the correct network under certain scenarios, and provides candidates for further exploration under other criteria and/or data in other scenarios.

Phylogeny, Morphology, and Metabolic and Invasive Capabilities of Epicellular Fish Coccidium Goussia janae

To fill the knowledge gap on the biology of the fish coccidian Goussia janae, RNA extracted from exogenously sporulated oocysts was sequenced. Analysis by Trinity and Trinotate pipelines showed that 84.6% of assembled transcripts share the highest similarity with Toxoplasma gondii and Neospora caninum. Phylogenetic and interpretive analyses from RNA-seq data provide novel insight into the metabolic capabilities, composition of the invasive machinery and the phylogenetic relationships of this parasite of cold-blooded vertebrates with other coccidians. This allows re-evaluation of the phylogenetic position of G. janae and sheds light on the emergence of the highly successful obligatory intracellularity of apicomplexan parasites. G. janae possesses a partial glideosome and along with it, the metabolic capabilities and adaptions of G. janae might provide cues as to how apicomplexans adjusted to extra- or intra-cytoplasmic niches and also to become obligate intracellular parasites. Unlike the similarly localized epicellular Cryptosporidium spp., G. janae lacks the feeder organelle necessary for directly scavenging nutrients from the host. Transcriptome analysis indicates that G. janae possesses metabolic capabilities comparable to T. gondii. Additionally, this enteric coccidium might also access host cell nutrients given the presence of a recently identified gene encoding the molecular sieve at the parasitophorous vacuole membrane.

Shared elements of host-targeting pathways among apicomplexan parasites of differing lifestyles

Apicomplexans are a diverse group of obligate parasites occupying different intracellular niches that require modification to meet the needs of the parasite. To efficiently manipulate their environment, apicomplexans translocate numerous parasite proteins into the host cell. Whereas some parasites remain contained within a parasitophorous vacuole membrane (PVM) throughout their developmental cycle, others do not, a difference that affects the machinery needed for protein export. A signal-mediated pathway for protein export into the host cell has been characterized in Plasmodium parasites, which maintain the PVM. Here, we functionally demonstrate an analogous host-targeting pathway involving organellar staging prior to secretion in the related bovine parasite, Babesia bovis, a parasite that destroys the PVM shortly after invasion. Taking into account recent identification of a similar signal-mediated pathway in the coccidian parasite Toxoplasma gondii, we suggest a model in which this conserved pathway has evolved in multiple steps from signal-mediated trafficking to specific secretory organelles for controlled secretion to a complex protein translocation process across the PVM.

Parasite Calcineurin Regulates Host Cell Recognition and Attachment by Apicomplexans

Apicomplexans invade a variety of metazoan host cells through mechanisms involving host cell receptor engagement and secretion of parasite factors to facilitate cellular attachment. We find that the parasite homolog of calcineurin, a calcium-regulated phosphatase complex central to signal transduction in eukaryotes, also contributes to host cell invasion by the malaria parasite Plasmodium falciparum and related Toxoplasma gondii. Using reverse-genetic and chemical-genetic approaches, we determine that calcineurin critically regulates and stabilizes attachment of extracellular P. falciparum to host erythrocytes before intracellular entry and has similar functions in host cell engagement by T. gondii. Calcineurin-mediated Plasmodium invasion is strongly associated with host receptors required for host cell recognition, and calcineurin function distinguishes this form of receptor-mediated attachment from a second mode of host-parasite adhesion independent of host receptors. This specific role of calcineurin in coordinating physical interactions with host cells highlights an ancestral mechanism for parasitism used by apicomplexans.

The ghost plastid of Choreocolax polysiphoniae

Parasitism has evolved innumerable times among eukaryotes. Red algal parasites alone have independently evolved over 100 times. The accepted evolutionary paradigm proposes that red algal parasites arise by first infecting a close relative and over time diversifying and infecting more distantly related species. This provides a natural evolutionary gradient of relationships between hosts and parasites that share a photosynthetic common ancestor. Upon infection, the parasite deposits its organelles into the host cell and takes over, spreading through cell-cell connections. Microscopy and molecular studies have demonstrated that the parasites do not maintain their own plastid, but rather abscond with a dedifferentiated host plastid as they pack up spores for dispersal. We sequenced a ~90 kb plastid genome from the parasite Choreocolax polysiphoniae, which has lost genes for light harvesting and photosynthesis. Furthermore, the presence of a native C. polysiphoniae plastid indicates that not all red algal parasites follow the same evolutionary pathway to parasitism. Along with the 167 kb plastid genome of its host, Vertebrata lanosa, these plastids are the first to be sequenced from the Ceramiales.

A lactate and formate transporter in the intraerythrocytic malaria parasite, Plasmodium falciparum

The intraerythrocytic malaria parasite relies primarily on glycolysis to fuel its rapid growth and reproduction. The major byproduct of this metabolism, lactic acid, is extruded into the external medium. In this study, we show that the human malaria parasite Plasmodium falciparum expresses at its surface a member of the microbial formate-nitrite transporter family (PfFNT), which, when expressed in Xenopus laevis oocytes, transports both formate and lactate. The transport characteristics of PfFNT in oocytes (pH-dependence, inhibitor-sensitivity and kinetics) are similar to those of the transport of lactate and formate across the plasma membrane of mature asexual-stage P. falciparum trophozoites, consistent with PfFNT playing a major role in the efflux of lactate and hence in the energy metabolism of the intraerythrocytic parasite.

Ciliate species diversity and host-parasitoid codiversification in Pseudocollinia infecting krill, with description of Pseudocollinia similis sp. nov

All parasitoid apostome ciliates infecting krill in the northeastern Pacific are currently assigned to the genus Pseudocollinia. Each krill specimen is apparently infected by only 1 Pseudocollinia species. We describe Pseudocollinia similis sp. nov., discovered infecting the krill Thysanoessa spinifera off Oregon, USA. Its protomite-tomite stage resembles that of P. beringensis, which infects T. inermis (type host species), T. longipes, and T. raschii females in the Bering Sea. These ciliates have similar numbers of somatic kineties (18-21 vs. 16-20) and typically have 3 oral kineties. Furthermore, these 2 apostomes are sister species on gene trees based on sequences of small subunit rRNA (0.06% difference) and cytochrome c oxidase subunit 1 (cox1; 30% difference). P. brintoni and P. oregonensis are closely related as a separate group from P. similis and P. beringensis. The similar tree topologies based on the cox1 sequences of 21 host krill individuals representing 6 krill species (Euphausia pacifica, Nyctiphanes simplex, T. inermis, T. longipes, T. raschii, and T. spinifera) and the apostomes isolated from these krill suggest host-parasitoid codiversification. However, this hypothesis was statistically rejected by an approximately unbiased test in which the host tree topology was used to model parasitoid evolution (p ≤ 0.05).

kdetrees: Non-parametric estimation of phylogenetic tree distributions

MOTIVATION

Although the majority of gene histories found in a clade of organisms are expected to be generated by a common process (e.g. the coalescent process), it is well known that numerous other coexisting processes (e.g. horizontal gene transfers, gene duplication and subsequent neofunctionalization) will cause some genes to exhibit a history distinct from those of the majority of genes. Such 'outlying' gene trees are considered to be biologically interesting, and identifying these genes has become an important problem in phylogenetics.

RESULTS

We propose and implement kdetrees, a non-parametric method for estimating distributions of phylogenetic trees, with the goal of identifying trees that are significantly different from the rest of the trees in the sample. Our method compares favorably with a similar recently published method, featuring an improvement of one polynomial order of computational complexity (to quadratic in the number of trees analyzed), with simulation studies suggesting only a small penalty to classification accuracy. Application of kdetrees to a set of Apicomplexa genes identified several unreliable sequence alignments that had escaped previous detection, as well as a gene independently reported as a possible case of horizontal gene transfer. We also analyze a set of Epichloë genes, fungi symbiotic with grasses, successfully identifying a contrived instance of paralogy.

AVAILABILITY AND IMPLEMENTATION

Our method for estimating tree distributions and identifying outlying trees is implemented as the R package kdetrees and is available for download from CRAN.

A survey of innovation through duplication in the reduced genomes of twelve parasites

We characterize the prevalence, distribution, divergence, and putative functions of detectable two-copy paralogs and segmental duplications in the Apicomplexa, a phylum of parasitic protists. Apicomplexans are mostly obligate intracellular parasites responsible for human and animal diseases (e.g. malaria and toxoplasmosis). Gene loss is a major force in the phylum. Genomes are small and protein-encoding gene repertoires are reduced. Despite this genomic streamlining, duplications and gene family amplifications are present. The potential for innovation introduced by duplications is of particular interest. We compared genomes of twelve apicomplexans across four lineages and used orthology and genome cartography to map distributions of duplications against genome architectures. Segmental duplications appear limited to five species. Where present, they correspond to regions enriched for multi-copy and species-specific genes, pointing toward roles in adaptation and innovation. We found a phylum-wide association of duplications with dynamic chromosome regions and syntenic breakpoints. Trends in the distribution of duplicated genes indicate that recent, species-specific duplicates are often tandem while most others have been dispersed by genome rearrangements. These trends show a relationship between genome architecture and gene duplication. Functional analysis reveals: proteases, which are vital to a parasitic lifecycle, to be prominent in putative recent duplications; a pair of paralogous genes in Toxoplasma gondii previously shown to produce the rate-limiting step in dopamine synthesis in mammalian cells, a possible link to the modification of host behavior; and phylum-wide differences in expression and subcellular localization, indicative of modes of divergence. We have uncovered trends in multiple modes of duplicate divergence including sequence, intron content, expression, subcellular localization, and functions of putative recent duplicates that highlight the role of duplications in the continuum of forces that have shaped these genomes.

Fast algorithms and heuristics for phylogenomics under ILS and hybridization

BACKGROUND

Phylogenomic analyses involving whole-genome or multi-locus data often entail dealing with incongruent gene trees. In this paper, we consider two causes of such incongruence, namely, incomplete lineage sorting (ILS) and hybridization, and consider both parsimony and probabilistic criteria for dealing with them.

RESULTS

Under the assumption of ILS, computing the probability of a gene tree given a species tree is a very hard problem. We present a heuristic for speeding up the computation, and demonstrate how it scales up computations to data sizes that are not feasible to analyze using current techniques, while achieving very good accuracy. Further, under the assumption of both ILS and hybridization, computing the probability of a gene tree and parsimoniously reconciling it with a phylogenetic network are both very hard problems. We present two exact algorithms for these two problems that speed up existing techniques significantly and enable analyses of much larger data sets than is currently feasible.

CONCLUSION

Our heuristics and algorithms enable phylogenomic analyses of larger (in terms of numbers of taxa) data sets than is currently feasible. Further, our methods account for ILS and hybridization, thus allowing analyses of reticulate evolutionary histories.

Malaria's many mates: past, present, and future of the systematics of the order Haemosporida

Malaria has been one of the most important diseases of humans throughout history and continues to be a major public health concern. The 5 species of Plasmodium that cause the disease in humans are part of the order Haemosporida, a diverse group of parasites that all have heteroxenous life cycles, alternating between a vertebrate host and a free-flying, blood-feeding dipteran vector. Traditionally, the identification and taxonomy of these parasites relied heavily on life-history characteristics, basic morphological features, and the host species infected. However, molecular approaches to resolving the phylogeny of the group have sometimes challenged many of these traditional hypotheses. One of the greatest debates has concerned the origin of the most virulent of the human-infecting parasites, Plasmodium falciparum, with early results suggesting a close relationship with an avian parasite. Subsequent phylogenetic studies placed it firmly within the mammalian clade instead, but the avian origin hypothesis has been revived with recent genome-based analyses. The rooting of the tree of Haemosporida has also been inconsistent, and the various topologies that result certainly affect our interpretation of the history of the group. There is clearly a pressing need to obtain a much more complete degree of taxon sampling of haemosporidians, as well as a greater number of characters before confidence can be placed in any hypothesis regarding the evolutionary history of the order. There are numerous challenges moving forward, particularly for generating complete genome sequences of avian and saurian parasites.

A novel clade of unique eukaryotic ribonucleotide reductase R2 subunits is exclusive to apicomplexan parasites

Apicomplexa are protist parasites of tremendous medical and economic importance, causing millions of deaths and billions of dollars in losses each year. Apicomplexan-related diseases may be controlled via inhibition of essential enzymes. Ribonucleotide reductase (RNR) provides the only de novo means of synthesizing deoxyribonucleotides, essential precursors for DNA replication and repair. RNR has long been the target of antibacterial and antiviral therapeutics. However, targeting this ubiquitous protein in eukaryotic pathogens may be problematic unless these proteins differ significantly from that of their respective host. The typical eukaryotic RNR enzymes belong to class Ia, and the holoenzyme consists minimally of two R1 and two R2 subunits (α₂β₂). We generated a comparative, annotated, structure-based, multiple-sequence alignment of R2 subunits, identified a clade of R2 subunits unique to Apicomplexa, and determined its phylogenetic position. Our analyses revealed that the apicomplexan-specific sequences share characteristics with both class I R2 and R2lox proteins. The putative radical-harboring residue, essential for the reduction reaction by class Ia R2-containing holoenzymes, was not conserved within this group. Phylogenetic analyses suggest that class Ia subunits are not monophyletic and consistently placed the apicomplexan-specific clade sister to the remaining class Ia eukaryote R2 subunits. Our research suggests that the novel apicomplexan R2 subunit may be a promising candidate for chemotherapeutic-induced inhibition as it differs greatly from known eukaryotic host RNRs and may be specifically targeted.

Parsimonious inference of hybridization in the presence of incomplete lineage sorting

Hybridization plays an important evolutionary role in several groups of organisms. A phylogenetic approach to detect hybridization entails sequencing multiple loci across the genomes of a group of species of interest, reconstructing their gene trees, and taking their differences as indicators of hybridization. However, methods that follow this approach mostly ignore population effects, such as incomplete lineage sorting (ILS). Given that hybridization occurs between closely related organisms, ILS may very well be at play and, hence, must be accounted for in the analysis framework. To address this issue, we present a parsimony criterion for reconciling gene trees within the branches of a phylogenetic network, and a local search heuristic for inferring phylogenetic networks from collections of gene-tree topologies under this criterion. This framework enables phylogenetic analyses while accounting for both hybridization and ILS. Further, we propose two techniques for incorporating information about uncertainty in gene-tree estimates. Our simulation studies demonstrate the good performance of our framework in terms of identifying the location of hybridization events, as well as estimating the proportions of genes that underwent hybridization. Also, our framework shows good performance in terms of efficiency on handling large data sets in our experiments. Further, in analysing a yeast data set, we demonstrate issues that arise when analysing real data sets. Although a probabilistic approach was recently introduced for this problem, and although parsimonious reconciliations have accuracy issues under certain settings, our parsimony framework provides a much more computationally efficient technique for this type of analysis. Our framework now allows for genome-wide scans for hybridization, while also accounting for ILS.

Targeted enrichment: maximizing orthologous gene comparisons across deep evolutionary time

Estimated phylogenies of evolutionarily diverse taxa will be well supported and more likely to be historically accurate when the analysis contains large amounts of data–many genes sequenced across many taxa. Inferring such phylogenies for non-model organisms is challenging given limited resources for whole-genome sequencing. We take advantage of genomic data from a single species to test the limits of hybridization-based enrichment of hundreds of exons across frog species that diverged up to 250 million years ago. Enrichment success for a given species depends greatly on the divergence time between it and the reference species, and the resulting alignment contains a significant proportion of missing data. However, our alignment generates a well-supported phylogeny of frogs, suggesting that this technique is a practical solution towards resolving relationships across deep evolutionary time.

Algorithms for genome-scale phylogenetics using gene tree parsimony

The use of genomic data sets for phylogenetics is complicated by the fact that evolutionary processes such as gene duplication and loss, or incomplete lineage sorting (deep coalescence) cause incongruence among gene trees. One well-known approach that deals with this complication is gene tree parsimony, which, given a collection of gene trees, seeks a species tree that requires the smallest number of evolutionary events to explain the incongruence of the gene trees. However, a lack of efficient algorithms has limited the use of this approach. Here, we present efficient algorithms for SPR and TBR-based local search heuristics for gene tree parsimony under the 1) duplication, 2) loss, 3) duplication-loss, and 4) deep coalescence reconciliation costs. These novel algorithms improve upon the time complexities of previous algorithms for these problems by a factor of n, where n is the number of species in the collection of gene trees. Our algorithms provide a substantial improvement in runtime and scalability compared to previous implementations and enable large-scale gene tree parsimony analyses using any of the four reconciliation costs. Our algorithms have been implemented in the software packages DupTree and iGTP, and have already been used to perform several compelling phylogenetic studies.

High-throughput RNA sequencing profiles and transcriptional evidence of aerobic respiratory enzymes in sporulating oocysts and sporozoites of Eimeria tenella

Seven species of Eimeria are responsible for coccidiosis in chickens. Eimeria tenella is one of the most pathogenic parasites since it is associated with high mortality and great economic impact. The life cycle of the parasite includes development in the environment and in the intestinal tract. We conducted RNA sequencing using a next generation sequencer to obtain transcriptome information from the sporulating oocysts, and sporozoites. We collected 2.8 million 75 bp reads of a short-tag sequence, and 25,880 contigs were generated by the Oases assembler. A Blastx search of GenBank databases revealed that 7780 contigs (30.1%) had significant homology with deposited sequence data (E-value <1e-6); among these contigs, 6051 contigs were similar to those of Toxoplasma gondii while only 513 contigs (6.6%) were similar to those of E. tenella. After an orthological analysis conducted with the UniProt database of T. gondii, 6661 contigs were distributed within the categories of cellular components (1528 gene categories), biological processes (861 gene categories), and molecular functions (241 gene categories). The significantly matched contigs contained high numbers of enzymes associated with glycolysis, TCA, and the pentose-phosphate pathway. Most of the enzymes, measured by quantitative reverse transcription-PCR, were up-regulated in sporulating stage. These results suggest that the intracellular carbohydrate amylopectin could be used as an energy source for ATP production including glycolysis and the pentose-phosphate pathway, which generates NADPH and pentoses. Our data also suggest that Eimeria might possess a partial or similar pathway to the TCA cycle essential for aerobic respiration. Furthermore, the newly annotated and non-annotated contigs might contain E. tenella-specific or novel sequences.

Cryptosporidium pathogenicity and virulence

Cryptosporidium is a protozoan parasite of medical and veterinary importance that causes gastroenteritis in a variety of vertebrate hosts. Several studies have reported different degrees of pathogenicity and virulence among Cryptosporidium species and isolates of the same species as well as evidence of variation in host susceptibility to infection. The identification and validation of Cryptosporidium virulence factors have been hindered by the renowned difficulties pertaining to the in vitro culture and genetic manipulation of this parasite. Nevertheless, substantial progress has been made in identifying putative virulence factors for Cryptosporidium. This progress has been accelerated since the publication of the Cryptosporidium parvum and C. hominis genomes, with the characterization of over 25 putative virulence factors identified by using a variety of immunological and molecular techniques and which are proposed to be involved in aspects of host-pathogen interactions from adhesion and locomotion to invasion and proliferation. Progress has also been made in the contribution of host factors that are associated with variations in both the severity and risk of infection. Here we provide a review comprised of the current state of knowledge on Cryptosporidium infectivity, pathogenesis, and transmissibility in light of our contemporary understanding of microbial virulence.

Aurora at the pole and equator: overlapping functions of Aurora kinases in the mitotic spindle

The correct assembly and timely disassembly of the mitotic spindle is crucial for the propagation of the genome during cell division. Aurora kinases play a central role in orchestrating bipolar spindle establishment, chromosome alignment and segregation. In most eukaryotes, ranging from amoebas to humans, Aurora activity appears to be required both at the spindle pole and the kinetochore, and these activities are often split between two different Aurora paralogues, termed Aurora A and B. Polar and equatorial functions of Aurora kinases have generally been considered separately, with Aurora A being mostly involved in centrosome dynamics, whereas Aurora B coordinates kinetochore attachment and cytokinesis. However, double inactivation of both Aurora A and B results in a dramatic synergy that abolishes chromosome segregation. This suggests that these two activities jointly coordinate mitotic progression. Accordingly, recent evidence suggests that Aurora A and B work together in both spindle assembly in metaphase and disassembly in anaphase. Here, we provide an outlook on these shared functions of the Auroras, discuss the evolution of this family of mitotic kinases and speculate why Aurora kinase activity may be required at both ends of the spindle microtubules.

Inference of reticulate evolutionary histories by maximum likelihood: the performance of information criteria

BACKGROUND

Maximum likelihood has been widely used for over three decades to infer phylogenetic trees from molecular data. When reticulate evolutionary events occur, several genomic regions may have conflicting evolutionary histories, and a phylogenetic network may provide a more adequate model for representing the evolutionary history of the genomes or species. A maximum likelihood (ML) model has been proposed for this case and accounts for both mutation within a genomic region and reticulation across the regions. However, the performance of this model in terms of inferring information about reticulate evolution and properties that affect this performance have not been studied.

RESULTS

In this paper, we study the effect of the evolutionary diameter and height of a reticulation event on its identifiability under ML. We find both of them, particularly the diameter, have a significant effect. Further, we find that the number of genes (which can be generalized to the concept of "non-recombining genomic regions") that are transferred across a reticulation edge affects its detectability. Last but not least, a fundamental challenge with phylogenetic networks is that they allow an arbitrary level of complexity, giving rise to the model selection problem. We investigate the performance of two information criteria, the Akaike Information Criterion (AIC) and the Bayesian Information Criterion (BIC), for addressing this problem. We find that BIC performs well in general for controlling the model complexity and preventing ML from grossly overestimating the number of reticulation events.

CONCLUSION

Our results demonstrate that BIC provides a good framework for inferring reticulate evolutionary histories. Nevertheless, the results call for caution when interpreting the accuracy of the inference particularly for data sets with particular evolutionary features.

An evolutionary perspective on the kinome of malaria parasites

Malaria parasites belong to an ancient lineage that diverged very early from the main branch of eukaryotes. The approximately 90-member plasmodial kinome includes a majority of eukaryotic protein kinases that clearly cluster within the AGC, CMGC, TKL, CaMK and CK1 groups found in yeast, plants and mammals, testifying to the ancient ancestry of these families. However, several hundred millions years of independent evolution, and the specific pressures brought about by first a photosynthetic and then a parasitic lifestyle, led to the emergence of unique features in the plasmodial kinome. These include taxon-restricted kinase families, and unique peculiarities of individual enzymes even when they have homologues in other eukaryotes. Here, we merge essential aspects of all three malaria-related communications that were presented at the Evolution of Protein Phosphorylation meeting, and propose an integrated discussion of the specific features of the parasite's kinome and phosphoproteome.

Sequencing of the smallest Apicomplexan genome from the human pathogen Babesia microti

We have sequenced the genome of the emerging human pathogen Babesia microti and compared it with that of other protozoa. B. microti has the smallest nuclear genome among all Apicomplexan parasites sequenced to date with three chromosomes encoding ∼3500 polypeptides, several of which are species specific. Genome-wide phylogenetic analyses indicate that B. microti is significantly distant from all species of Babesidae and Theileridae and defines a new clade in the phylum Apicomplexa. Furthermore, unlike all other Apicomplexa, its mitochondrial genome is circular. Genome-scale reconstruction of functional networks revealed that B. microti has the minimal metabolic requirement for intraerythrocytic protozoan parasitism. B. microti multigene families differ from those of other protozoa in both the copy number and organization. Two lateral transfer events with significant metabolic implications occurred during the evolution of this parasite. The genomic sequencing of B. microti identified several targets suitable for the development of diagnostic assays and novel therapies for human babesiosis.