Towards a reference map of Eimeria tenella sporozoite proteins by two-dimensional electrophoresis and mass spectrometry

Eimeria tenella is a parasite of great importance as a disease causing agent in the poultry industry. Until recently, biological studies have focused on specific proteins, some of which play an important role in the parasite life cycle. Post-genomic studies will make it possible to understand the complexity of the parasites and their interactions with host cells. Here we present a systematic reference map of the proteins from E. tenella sporozoites. The proteins expressed at the sporozoite stage were resolved between isoelectric points 3-10 and 4-7. They were systematically identified using mass spectrometry and 16 known Eimeria sporozoite proteins were identified on two-dimensional maps. Peptide fragmentation data from mass spectrometry were compared to single and consensus expression sequence tags in databases and to the E. tenella genome (not annotated). Among the set of unknown proteins analysed, 12 new assignments were proposed on the basis of similarities with Apicomplexa proteins. In order to define sporozoite proteins as potential targets for coccidiosis therapy, proteins were studied according to their relative abundance and immunogenicity in the sporozoite. Immunoblots of sporozoite 2D maps with chicken sera were performed and approximately 50 proteins were defined as antigens. It was shown that abundance and immunogenicity are not related in the sporozoite stage. Perspectives of gene prediction and completion of the genome annotation by a proteomic approach is discussed.     

BAC-HAPPY Mapping (BAP Mapping): A New and Efficient Protocol for Physical Mapping

Physical and linkage mapping underpin efforts to sequence and characterize the genomes of eukaryotic organisms by providing a skeleton framework for whole genome assembly. Hitherto, linkage and physical “contig” maps were generated independently prior to merging. Here, we develop a new and easy method, BAC HAPPY MAPPING (BAP mapping), that utilizes BAC library pools as a HAPPY mapping panel together with an Mbp-sized DNA panel to integrate the linkage and physical mapping efforts into one pipeline. Using Arabidopsis thaliana as an exemplar, a set of 40 Sequence Tagged Site (STS) markers spanning ∼10% of chromosome 4 were simultaneously assembled onto a BAP map compiled using both a series of BAC pools each comprising 0.7x genome coverage and dilute (0.7x genome) samples of sheared genomic DNA. The resultant BAP map overcomes the need for polymorphic loci to separate genetic loci by recombination and allows physical mapping in segments of suppressed recombination that are difficult to analyze using traditional mapping techniques. Even virtual “BAC-HAPPY-mapping” to convert BAC landing data into BAC linkage contigs is possible.     

Eimeria tenella: identification of secretory and surface proteins from expressed sequence tags

To identify new vaccine candidates, Eimeria tenella expressed sequence tags (ESTs) from public databases were analysed for secretory molecules with an especially developed automated in silico strategy termed DNAsignalP. A total of 12,187 ESTs were clustered into 2881 contigs followed by a blastx search, which resulted in a significant number of E. tenella contigs with homologies to entries in public databases. Amino acid sequences of appropriate homologous proteins were analysed for the occurrence of an N-terminal signal sequence using the algorithm signalP. The resulting list of 84 entries comprised 51 contigs whose deduced proteins showed homologies to proteins of apicomplexan parasites. Based on function or localisation, we selected candidate proteins classified as (i) secreted proteins of Apicomplexa parasites, (ii) secreted enzymes, and (iii) transport and signalling proteins. To verify our strategy experimentally, we used a functional complementation system in yeast. For five selected candidate proteins we found that these were indeed secreted. Our approach thus represents an efficient method to identify secretory and surface proteins out of EST databases.     

Old can be new again: HAPPY whole genome sequencing, mapping and assembly

During the last three decades, both genome mapping and sequencing methods have advanced significantly to provide a foundation for scientists to understand genome structures and functions in many species. Generally speaking, genome mapping relies on genome sequencing to provide basic materials, such as DNA probes and markers for their localizations, thus constructing the maps. On the other hand, genome sequencing often requires a high-resolution map as a skeleton for whole genome assembly. However, both genome mapping and sequencing have never come together in one pipeline. After reviewing mapping and next-generation sequencing methods, we would like to share our thoughts with the genome community on how to combine the HAPPY mapping technique with the new-generation sequencing, thus integrating two systems into one pipeline, called HAPPY pipeline. The pipeline starts with preparation of a HAPPY panel, followed by multiple displacement amplification for producing a relatively large quantity of DNA. Instead of conventional marker genotyping, the amplified panel DNA samples are subject to new-generation sequencing with barcode method, which allows us to determine the presence/absence of a sequence contig as a traditional marker in the HAPPY panel. Statistical analysis will then be performed to infer how close or how far away from each other these contigs are within a genome and order the whole genome sequence assembly as well. We believe that such a universal approach will play an important role in genome sequencing, mapping, and assembly of many species; thus advancing genome science and its applications in biomedicine and agriculture.     

Single-molecule analysis of genome rearrangements in cancer

Rearrangements of the genome can be detected by microarray methods and massively parallel sequencing, which identify copy-number alterations and breakpoint junctions, but these techniques are poorly suited to reconstructing the long-range organization of rearranged chromosomes, for example, to distinguish between translocations and insertions. The single-DNA-molecule technique HAPPY mapping is a method for mapping normal genomes that should be able to analyse genome rearrangements, i.e. deviations from a known genome map, to assemble rearrangements into a long-range map. We applied HAPPY mapping to cancer cell lines to show that it could identify rearrangement of genomic segments, even in the presence of normal copies of the genome. We could distinguish a simple interstitial deletion from a copy-number loss at an inversion junction, and detect a known translocation. We could determine whether junctions detected by sequencing were on the same chromosome, by measuring their linkage to each other, and hence map the rearrangement. Finally, we mapped an uncharacterized reciprocal translocation in the T-47D breast cancer cell line to about 2 kb and hence cloned the translocation junctions. We conclude that HAPPY mapping is a versatile tool for determining the structure of rearrangements in the human genome.     

A physical map of the highly heterozygous Populus genome: integration with the genome sequence and genetic map and analysis of haplotype variation

As part of a larger project to sequence the Populus genome and generate genomic resources for this emerging model tree, we constructed a physical map of the Populus genome, representing one of the few such maps of an undomesticated, highly heterozygous plant species. The physical map, consisting of 2802 contigs, was constructed from fingerprinted bacterial artificial chromosome (BAC) clones. The map represents approximately 9.4-fold coverage of the Populus genome, which has been estimated from the genome sequence assembly to be 485 ± 10 Mb in size. BAC ends were sequenced to assist long-range assembly of whole-genome shotgun sequence scaffolds and to anchor the physical map to the genome sequence. Simple sequence repeat-based markers were derived from the end sequences and used to initiate integration of the BAC and genetic maps. A total of 2411 physical map contigs, representing 97% of all clones assigned to contigs, were aligned to the sequence assembly (JGI Populus trichocarpa , version 1.0). These alignments represent a total coverage of 384 Mb (79%) of the entire poplar sequence assembly and 295 Mb (96%) of linkage group sequence assemblies. A striking result of the physical map contig alignments to the sequence assembly was the co-localization of multiple contigs across numerous regions of the 19 linkage groups. Targeted sequencing of BAC clones and genetic analysis in a small number of representative regions showed that these co-aligning contigs represent distinct haplotypes in the heterozygous individual sequenced, and revealed the nature of these haplotype sequence differences.     

Generation of Physical Map Contig-Specific Sequences Useful for Whole Genome Sequence Scaffolding

Along with the rapid advances of the nextgen sequencing technologies, more and more species are added to the list of organisms whose whole genomes are sequenced. However, the assembled draft genome of many organisms consists of numerous small contigs, due to the short length of the reads generated by nextgen sequencing platforms. In order to improve the assembly and bring the genome contigs together, more genome resources are needed. In this study, we developed a strategy to generate a valuable genome resource, physical map contig-specific sequences, which are randomly distributed genome sequences in each physical contig. Two-dimensional tagging method was used to create specific tags for 1,824 physical contigs, in which the cost was dramatically reduced. A total of 94,111,841 100-bp reads and 315,277 assembled contigs are identified containing physical map contig-specific tags. The physical map contig-specific sequences along with the currently available BAC end sequences were then used to anchor the catfish draft genome contigs. A total of 156,457 genome contigs (~79% of whole genome sequencing assembly) were anchored and grouped into 1,824 pools, in which 16,680 unique genes were annotated. The physical map contig-specific sequences are valuable resources to link physical map, genetic linkage map and draft whole genome sequences, consequently have the capability to improve the whole genome sequences assembly and scaffolding, and improve the genome-wide comparative analysis as well. The strategy developed in this study could also be adopted in other species whose whole genome assembly is still facing a challenge.     

Global protein expression analysis in apicomplexan parasites: current status

Members of the phylum Apicomplexa are important protozoan parasites that cause some of the most serious, and in some cases, deadly diseases in humans and animals. They include species from the genus Plasmodium, Toxoplasma, Eimeria, Neospora, Cryptosporidium, Babesia and Theileria. The medical, veterinary and economic impact of these pathogens on a global scale is enormous. Although chemo- and immuno-prophylactic strategies are available to control some of these parasites, they are inadequate. Currently, there is an urgent need to design new vaccines or chemotherapeutics for apicomplexan diseases. High-throughput global protein expression analyses using gel or non-gel based protein separation technologies coupled with mass spectrometry and bioinformatics provide a means to identify new drug and vaccine targets in these pathogens. Protein identification based proteomic projects in apicomplexan parasites is currently underway, with the most significant progress made in the malaria parasite, Plasmodium falciparum. More recently, preliminary two-dimensional gel electrophoresis maps of Toxoplasma gondii and Neospora caninum tachyzoites and Eimeria tenella sporozoites, have been produced, as well as for micronemes in E. tenella. In this review, the status of proteomics in the analysis of global protein expression in apicomplexan parasites will be compared and the challenges associated with these investigations discussed.     

Integrating genetics and genomics to identify new leads for the control of Eimeria spp

Eimerian parasites display a biologically interesting range of phenotypic variation. In addition to a wide spectrum of drug-resistance phenotypes that are expressed similarly by many other parasites, the Eimeria spp. present some unique phenotypes. For example, unique lines of Eimeria spp. include those selected for growth in the chorioallantoic membrane of the embryonating hens egg or for faster growth (precocious development) in the mature host. The many laboratory-derived egg-adapted or precocious lines also share a phenotype of a marked attenuation of virulence, the basis of which is different as a consequence of the in ovo or in vivo selection procedures used. Of current interest is the fact that some wild-type populations of Eimeria maxima are characterized by an ability to induce protective immunity that is strain-specific. The molecular basis of phenotypes that define Eimeria spp. is now increasingly amenable to investigation, both through technical improvements in genetic linkage studies and the availability of a comprehensive genome sequence for the caecal parasite E. tenella. The most exciting phenotype in the context of vaccination and the development of new vaccines is the trait of strain-specific immunity associated with E. maxima. Recent work in this laboratory has shown that infection of two inbred lines of White Leghorn chickens with the W strain of E. maxima leads to complete protection to challenge with the homologous parasite, but to complete escape of the heterologous H strain, i.e. the W strain induces an exquisitely strain-specific protective immune response with respect to the H strain. This dichotomy of survival in the face of immune-mediated killing has been examined further and, notably, mating between a drug-resistant W strain and a drug-sensitive H strain leads to recombination between the genetic loci responsible for the specificity of protective immunity and resistance to the anticoccidial drug robenidine. Such a finding opens the way forward for genetic mapping of the loci responsible for the induction of protective immunity and integration with the genome sequencing efforts.     

The genome of Eimeria spp., with special reference to Eimeria tenella--a coccidium from the chicken

Eimeria spp. contain at least four genomes. The nuclear genome is best studied in the avian species Eimeria tenella and comprises about 60 Mbp DNA contained within ca. 14 chromosomes; other avian and lupine species appear to possess a nuclear genome of similar size. In addition, sequence data and hybridisation studies have provided direct evidence for extrachromosomal mitochondrial and plastid DNA genomes, and double-stranded RNA segments have also been described. The unique phenotype of "precocious" development that characterises some selected lines of Eimeria spp. not only provides the basis for the first generation of live attenuated vaccines, but offers a significant entrée into studies on the regulation of an apicomplexan life-cycle. With a view to identifying loci implicated in the trait of precocious development, a genetic linkage map of the genome of E. tenella is being constructed in this laboratory from analyses of the inheritance of over 400 polymorphic DNA markers in the progeny of a cross between complementary drug-resistant and precocious parents. Other projects that impinge directly or indirectly on the genome and/or genetics of Eimeria spp. are currently in progress in several laboratories, and include the derivation of expressed sequence tag data and the development of ancillary technologies such as transfection techniques. No large-scale genomic DNA sequencing projects have been reported.     

Composite genome map and recombination parameters derived from three archetypal lineages of Toxoplasma gondii

Toxoplasma gondii is a highly successful protozoan parasite in the phylum Apicomplexa, which contains numerous animal and human pathogens. T.gondii is amenable to cellular, biochemical, molecular and genetic studies, making it a model for the biology of this important group of parasites. To facilitate forward genetic analysis, we have developed a high-resolution genetic linkage map for T.gondii. The genetic map was used to assemble the scaffolds from a 10X shotgun whole genome sequence, thus defining 14 chromosomes with markers spaced at ~300 kb intervals across the genome. Fourteen chromosomes were identified comprising a total genetic size of ~592 cM and an average map unit of ~104 kb/cM. Analysis of the genetic parameters in T.gondii revealed a high frequency of closely adjacent, apparent double crossover events that may represent gene conversions. In addition, we detected large regions of genetic homogeneity among the archetypal clonal lineages, reflecting the relatively few genetic outbreeding events that have occurred since their recent origin. Despite these unusual features, linkage analysis proved to be effective in mapping the loci determining several drug resistances. The resulting genome map provides a framework for analysis of complex traits such as virulence and transmission, and for comparative population genetic studies.     

Genome-wide SNP identification by high-throughput sequencing and selective mapping allows sequence assembly positioning using a framework genetic linkage map

Background  

Determining the position and order of contigs and scaffolds from a genome assembly within an organism's genome remains a technical challenge in a majority of sequencing projects. In order to exploit contemporary technologies for DNA sequencing, we developed a strategy for whole genome single nucleotide polymorphism sequencing allowing the positioning of sequence contigs onto a linkage map using the bin mapping method.     

The microneme proteins EtMIC4 and EtMIC5 of Eimeria tenella form a novel, ultra-high molecular mass protein complex that binds target host cells

Eimeria tenella, in common with other parasitic protozoa of the phylum Apicomplexa, invades host cells using an actinomyosin-powered "glideosome" complex and requires the secretion of adhesive proteins from the microneme organelles onto the parasite surface. Microneme proteins of E. tenella include EtMIC4, a transmembrane protein that has multiple thrombospondin type I domains and calcium-binding epidermal growth factor-like domains in its extracellular domain, and EtMIC5, a soluble protein composed of 11 tandemly repeated domains that belong to the plasminogen-apple-nematode superfamily. We show here that EtMIC4 and EtMIC5 interact to form an oligomeric, ultrahigh molecular mass protein complex. The complex was purified from lysed parasites by non-denaturing techniques, and the stoichiometry was shown to be [EtMIC4](2):[EtMIC5](1), with an octamer of EtMIC4 bound non-covalently to a tetramer of EtMIC5. The complex is formed within the parasite secretory pathway and is maintained after secretion onto the surface of the parasite. The purified complex binds to a number of epithelial cell lines in culture. Identification and characterization of this complex contributes to an overall understanding of the role of multimolecular protein complexes in specific interactions between pathogens and their hosts during infection.     

Identification of an apically-located antigen that is conserved in sporozoan parasites

Sporozoan parasites of the phylum Apicomplexa all possess common apical structures. The current study used a monoclonal antibody (mAb-E12) to identify a conserved antigen in the apical region of merozoites of seven species of Plasmodium (including rodent, primate and human pathogens), tachyzoites of Toxoplasma gondii, bradyzoites of Sarcocystis bovis, and sporozoites and merozoites of Eimeria tenella and E. acervulina. The antigen was also present in sporozoites of haemosporinid parasites. Immunofluorescence studies showed that the antigen was restricted to the apical 3rd of these invasive stages. Using immunoelectron microscopy, labeling was demonstrated in the region of the polar ring, below the paired inner membranes of the parasite pellicle, and near the subpellicular microtubules radiating from the polar ring of merozoites and sporozoites of E. tenella. The majority of the antigen could be extracted with 1% Triton-X 100, but a portion remained associated with the cytoskeletal elements. The molecule has a relative rate of migration (Mr) of 47,000 in Plasmodium spp. and 43-46,000 in coccidian species. Since the epitope recognized by mAb-E12 is highly conserved, restricted to motile stages, and appears to be associated with microtubules, this antigen could be involved in cellular motility and cellular invasion.     

Characterization of monoclonal antibodies that recognize the Eimeria tenella microneme protein MIC2

The apicomplexan pathogens of Eimeria cause coccidiosis, an intestinal disease of chickens, which has a major economic impact on the poultry industry. Members of the Apicomplexa share an assortment of unique secretory organelles (rhoptries, micronemes and dense granules) that mediate invasion of host cells and formation and modification of the parasitophorous vacuole. Among these, microneme protein 2 from Eimeria tenella(EtMIC2) has a putative function in parasite adhesion to the host cell to initiate the invasion process. To investigate the role of EtMIC2 in host parasite interactions, the production and characterization of 12 monoclonal antibodies (mabs) produced against recombinant EtMIC2 proteins is described. All mabs reacted with molecules belonging to the apical complex of sporozoites and merozoites of E. tenella, E. acervulina and E. maxima in an immunofluorescence assay. By Western blot analysis, the mabs identified a developmentally regulated protein of 42 kDa corresponding to EtMIC 2 and cross-reacted with proteins in developmental stages of E. acervulina. Collectively, these mabs are useful tools for the detailed investigation of the characterization of EtMIC2 related proteins in Eimeria species.     

Enhanced De Novo Assembly of High Throughput Pyrosequencing Data Using Whole Genome Mapping

Despite major advances in next-generation sequencing, assembly of sequencing data, especially data from novel microorganisms or re-emerging pathogens, remains constrained by the lack of suitable reference sequences. De novo assembly is the best approach to achieve an accurate finished sequence, but multiple sequencing platforms or paired-end libraries are often required to achieve full genome coverage. In this study, we demonstrated a method to assemble complete bacterial genome sequences by integrating shotgun Roche 454 pyrosequencing with optical whole genome mapping (WGM). The whole genome restriction map (WGRM) was used as the reference to scaffold de novo assembled sequence contigs through a stepwise process. Large de novo contigs were placed in the correct order and orientation through alignment to the WGRM. De novo contigs that were not aligned to WGRM were merged into scaffolds using contig branching structure information. These extended scaffolds were then aligned to the WGRM to identify the overlaps to be eliminated and the gaps and mismatches to be resolved with unused contigs. The process was repeated until a sequence with full coverage and alignment with the whole genome map was achieved. Using this method we were able to achieved 100% WGRM coverage without a paired-end library. We assembled complete sequences for three distinct genetic components of a clinical isolate of Providencia stuartii: a bacterial chromosome, a novel bla _NDM-1 plasmid, and a novel bacteriophage, without separately purifying them to homogeneity.     

Myosin diversity in Apicomplexa

A polymerase chain reaction (PCR) screen was used to examine the diversity of myosins in 7 Apicomplexan parasites: Toxoplasma gondii, Plasmodium falciparum, Neospora caninum, Eimeria tenella, Sarcocystis muris, Babesia bovis, and Cryptosporidium parvum. Using degenerate PCR primers compatible with the majority of known myosin classes, putative myosin sequences were obtained from all of these species. All of the sequences obtained showed greatest similarity to previously identified apicomplexan myosins, suggesting that the diversity of myosins in these parasites is limited. Myosin classes that are known to be widespread across the phylogenetic spectrum, e.g., the myosins I, II, and V, were not seen in the Apicomplexa. Thus, like the plants, the Apicomplexa may have evolved their own unique cohort of myosins that are responsible for the myosin-driven cellular functions observed in these parasites.