A HAPPY map of Cryptosporidium parvum

Long-read assembly and comparative evidence-based reanalysis of Cryptosporidium genome sequences reveal expanded transporter repertoire and duplication of entire chromosome ends including subtelomeric regions

Cryptosporidiosis is a leading cause of waterborne diarrheal disease globally and an important contributor to mortality in infants and the immunosuppressed. Despite its importance, the Cryptosporidium community has only had access to a good, but incomplete, Cryptosporidium parvum IOWA reference genome sequence. Incomplete reference sequences hamper annotation, experimental design, and interpretation. We have generated a new C. parvum IOWA genome assembly supported by Pacific Biosciences (PacBio) and Oxford Nanopore long-read technologies and a new comparative and consistent genome annotation for three closely related species: C. parvum, Cryptosporidium hominis, and Cryptosporidium tyzzeri We made 1926 C. parvum annotation updates based on experimental evidence. They include new transporters, ncRNAs, introns, and altered gene structures. The new assembly and annotation revealed a complete Dnmt2 methylase ortholog. Comparative annotation between C. parvum, C. hominis, and C. tyzzeri revealed that most "missing" orthologs are found, suggesting that the biological differences between the species must result from gene copy number variation, differences in gene regulation, and single-nucleotide variants (SNVs). Using the new assembly and annotation as reference, 190 genes are identified as evolving under positive selection, including many not detected previously. The new C. parvum IOWA reference genome assembly is larger, gap free, and lacks ambiguous bases. This chromosomal assembly recovers all 16 chromosome ends, 13 of which are contiguously assembled. The three remaining chromosome ends are provisionally placed. These ends represent duplication of entire chromosome ends including subtelomeric regions revealing a new level of genome plasticity that will both inform and impact future research.

Challenges for Cryptosporidium Population Studies

Cryptosporidiosis is ranked sixth in the list of the most important food-borne parasites globally, and it is an important contributor to mortality in infants and the immunosuppressed. Recently, the number of genome sequences available for this parasite has increased drastically. The majority of the sequences are derived from population studies of Cryptosporidium parvum and Cryptosporidium hominis, the most important species causing disease in humans. Work with this parasite is challenging since it lacks an optimal, prolonged, in vitro culture system, which accurately reproduces the in vivo life cycle. This obstacle makes the cloning of isolates nearly impossible. Thus, patient isolates that are sequenced represent a population or, at times, mixed infections. Oocysts, the lifecycle stage currently used for sequencing, must be considered a population even if the sequence is derived from single-cell sequencing of a single oocyst because each oocyst contains four haploid meiotic progeny (sporozoites). Additionally, the community does not yet have a set of universal markers for strain typing that are distributed across all chromosomes. These variables pose challenges for population studies and require careful analyses to avoid biased interpretation. This review presents an overview of existing population studies, challenges, and potential solutions to facilitate future population analyses.

Development of Molecular Diagnosis Using Multiplex Real-Time PCR and T4 Phage Internal Control to Simultaneously Detect Cryptosporidium parvum, Giardia lamblia, and Cyclospora cayetanensis from Human Stool Samples

This study aimed to develop a new multiplex real-time PCR detection method for 3 species of waterborne protozoan parasites (Cryptosporidium parvum, Giardia lamblia, and Cyclospora cayetanensis) identified as major causes of traveler's diarrhea. Three target genes were specifically and simultaneously detected by the TaqMan probe method for multiple parasitic infection cases, including Cryptosporidium oocyst wall protein for C. parvum, glutamate dehydrogenase for G. lamblia, and internal transcribed spacer 1 for C. cayetanensis. Gene product 21 for bacteriophage T4 was used as an internal control DNA target for monitoring human stool DNA amplification. TaqMan probes were prepared using 4 fluorescent dyes, FAMTM, HEXTM, Cy5TM, and CAL Fluor Red® 610 on C. parvum, G. lamblia, C. cayetanensis, and bacteriophage T4, respectively. We developed a novel primer-probe set for each parasite, a primer-probe cocktail (a mixture of primers and probes for the parasites and the internal control) for multiplex real-time PCR analysis, and a protocol for this detection method. Multiplex real-time PCR with the primer-probe cocktail successfully and specifically detected the target genes of C. parvum, G. lamblia, and C. cayetanensis in the mixed spiked human stool sample. The limit of detection for our assay was 2×10 copies for C. parvum and for C. cayetanensis, while it was 2×103 copies for G. lamblia. We propose that the multiplex real-time PCR detection method developed here is a useful method for simultaneously diagnosing the most common causative protozoa in traveler's diarrhea.

Whole-Genome Restriction Mapping by "Subhaploid"-Based RAD Sequencing: An Efficient and Flexible Approach for Physical Mapping and Genome Scaffolding

Assembly of complex genomes using short reads remains a major challenge, which usually yields highly fragmented assemblies. Generation of ultradense linkage maps is promising for anchoring such assemblies, but traditional linkage mapping methods are hindered by the infrequency and unevenness of meiotic recombination that limit attainable map resolution. Here we develop a sequencing-based "in vitro" linkage mapping approach (called RadMap), where chromosome breakage and segregation are realized by generating hundreds of "subhaploid" fosmid/bacterial-artificial-chromosome clone pools, and by restriction site-associated DNA sequencing of these clone pools to produce an ultradense whole-genome restriction map to facilitate genome scaffolding. A bootstrap-based minimum spanning tree algorithm is developed for grouping and ordering of genome-wide markers and is implemented in a user-friendly, integrated software package (AMMO). We perform extensive analyses to validate the power and accuracy of our approach in the model plant Arabidopsis thaliana and human. We also demonstrate the utility of RadMap for enhancing the contiguity of a variety of whole-genome shotgun assemblies generated using either short Illumina reads (300 bp) or long PacBio reads (6-14 kb), with up to 15-fold improvement of N50 (∼816 kb-3.7 Mb) and high scaffolding accuracy (98.1-98.5%). RadMap outperforms BioNano and Hi-C when input assembly is highly fragmented (contig N50 = 54 kb). RadMap can capture wide-range contiguity information and provide an efficient and flexible tool for high-resolution physical mapping and scaffolding of highly fragmented assemblies.

Identification and discrimination of Toxoplasma gondii, Sarcocystis spp., Neospora spp., and Cryptosporidium spp. by righ-resolution melting analysis

The objective of this study was to standardize the high-resolution melting method for identification and discrimination of Toxoplasma gondii, Sarcocystis spp., Neospora spp., and Cryptosporidium spp. by amplification of 18S ribosomal DNA (rDNA) using a single primer pair. The analyses were performed on individual reactions (containing DNA from a single species of a protozoan), on duplex reactions (containing DNA from two species of protozoa in each reaction), and on a multiplex reaction (containing DNA of four parasites in a single reaction). The proposed method allowed us to identify and discriminate the four species by analyzing the derivative, normalized, and difference melting curves, with high reproducibility among and within the experiments, as demonstrated by low coefficients of variation (less than 2.2% and 2.0%, respectively). This is the first study where this method is used for discrimination of these four species of protozoa in a single reaction.

Comparative genomic analysis reveals occurrence of genetic recombination in virulent Cryptosporidium hominis subtypes and telomeric gene duplications in Cryptosporidium parvum

Background

Cryptosporidium hominis is a dominant species for human cryptosporidiosis. Within the species, IbA10G2 is the most virulent subtype responsible for all C. hominis–associated outbreaks in Europe and Australia, and is a dominant outbreak subtype in the United States. In recent yearsIaA28R4 is becoming a major new subtype in the United States. In this study, we sequenced the genomes of two field specimens from each of the two subtypes and conducted a comparative genomic analysis of the obtained sequences with those from the only fully sequenced Cryptosporidium parvum genome.

Results

Altogether, 8.59-9.05 Mb of Cryptosporidium sequences in 45–767 assembled contigs were obtained from the four specimens, representing 94.36-99.47% coverage of the expected genome. These genomes had complete synteny in gene organization and 96.86-97.0% and 99.72-99.83% nucleotide sequence similarities to the published genomes of C. parvum and C. hominis, respectively. Several major insertions and deletions were seen between C. hominis and C. parvum genomes, involving mostly members of multicopy gene families near telomeres. The four C. hominis genomes were highly similar to each other and divergent from the reference IaA25R3 genome in some highly polymorphic regions. Major sequence differences among the four specimens sequenced in this study were in the 5′ and 3′ ends of chromosome 6 and the gp60 region, largely the result of genetic recombination.

Conclusions

The sequence similarity among specimens of the two dominant outbreak subtypes and genetic recombination in chromosome 6, especially around the putative virulence determinant gp60 region, suggest that genetic recombination plays a potential role in the emergence of hyper-transmissible C. hominis subtypes. The high sequence conservation between C. parvum and C. hominis genomes and significant differences in copy numbers of MEDLE family secreted proteins and insulinase-like proteases indicate that telomeric gene duplications could potentially contribute to host expansion in C. parvum.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1517-1) contains supplementary material, which is available to authorized users.

Genetic mapping and coccidial parasites: past achievements and future prospects

Coccidial parasites including Cryptosporidium parvum, Cyclospora cayetanensis, Neospora caninum, Toxoplasma gondii and the Eimeria species can cause severe disease of medical and veterinary importance. As many as one-third of the human population may carry T. gondii infection, and Eimeria are thought to cost the global poultry production industry in excess of US$2 billion per annum. Despite their significance, effective vaccines are scarce and have been confined to the veterinary field. As sequencing and genotyping technologies continue to develop, genetic mapping remains a valuable tool for the identification of genes that underlie phenotypic traits of interest and the assembly of contiguous genome sequences. For the coccidian, cross-fertilization still requires in vivo infection, a feature of their life cycle which limits the use of genetic mapping strategies. Importantly, the development of population-based approaches has now removed the need to isolate clonal lines for genetic mapping of selectable traits, complementing the classical clone-based techniques. To date, four coccidial species, representing three genera, have been investigated using genetic mapping. In this review we will discuss recent progress with these species and examine the prospects for future initiatives.

The challenges of studying complex and dynamic regions of the human genome

Recent work has emphasised that the human genome is not simple and static, but complex and dynamic. This review focuses on the regions that are particularly hard to dissect and analyse, yet hold clues to how the genome changes during evolution and disease. I begin by summarising recent key advances in the understanding of the variable structure of our genome, and then I discuss a medley of methods that may allow us to analyse this structure in fine detail. In the final part, I describe potential future developments in this field, and make an argument that, just as we routinely genotype single-nucleotide polymorphisms now and will routinely re-sequence genomes in the near future, we should be aiming to physically re-map the individual human genome for each individual we study.

Development and initial characterization of a HAPPY panel for mapping the X. tropicalis genome

HAPPY mapping was designed to pursue the analysis of approximately random HAPloid DNA breakage samples using the PolYmerase chain reaction for mapping genomes. In the present study, we improved the method and integrated two other molecular techniques into the process: whole genome amplification and the Sequenom SNP (single nucleotide polymorphism) genotyping assay in order to facilitate whole genome mapping of X. tropicalis. The former technique amplified enough DNA materials to genotype a large number of markers, while the latter allowed for relatively high throughput marker genotyping with multiplex assays on the HAPPY lines. A total of 58 X. tropicalis genes were genotyped on an initial panel of 383 HAPPY lines, which contributed to formation of a working panel of 146 lines. Further genotyping of 29 markers on the working panel led to construction of a HAPPY map for the X. tropicalis genome. We believe that our improved HAPPY method described in the present study has paved the way for the community to map different genomes with a simple, but powerful approach.

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.

Comparative efficacy of conventional primer sets in detection of Cryptosporidium parvum for diagnostic use

In this study, the sensitivity and specificity of different previously described primer sets for Cryptosporidium parvum detection by polymerase chain reaction (PCR) was evaluated. For this purpose, the primer sets defined by Cacciò et al. (FEMS Microbiol Lett 170(1):173-179, 1999) (tub), Widmer et al. (Appl Environ Microbiol 64(11):4477-4481, 1998) (btub) and Rochelle et al. (Appl Environ Microbiol 63:2029-2037, 1997) (cphsp), respectively, were used. Deoxyribonucleic acid (DNA) was isolated from three different sample materials: (1) from the faeces of an experimentally C. parvum-infected calf, (2) from purified C. parvum oocysts, and (3) from C. parvum-infected HCT-8 cell cultures. The DNA samples were subjected to PCR reactions with each of the three given primer sets to investigate sensitivity and suitability for routine use. The primers described by Cacciò et al. (FEMS Microbiol Lett 170(1):173-179, 1999) (TUB) were superior regarding sensitivity and specificity in terms of detection of C. parvum in faeces, in purified oocysts and also in cell culture, and may thus be applied for routine diagnostic use in common sample materials.

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.

Cryptosporidium: genomic and biochemical features

Recent progress in understanding the unique biochemistry of the two closely related human enteric pathogens Cryptosporidium parvum and Cryptosporidium hominis has been stimulated by the elucidation of the complete genome sequences for both pathogens. Much of the work that has occurred since that time has been focused on understanding the metabolic pathways encoded by the genome in hopes of providing increased understanding of the parasite biology, and in the identification of novel targets for pharmacological interventions. However, despite identifying the genes encoding enzymes that participate in many of the major metabolic pathways, only a hand full of proteins have actually been the subjects of detailed scrutiny. Thus, much of the biochemistry of these parasites remains a true mystery.

A hundred-year retrospective on cryptosporidiosis

Tyzzer discovered the genus Cryptosporidium a century ago, and for almost 70 years cryptosporidiosis was regarded as an infrequent and insignificant infection that occurred in the intestines of vertebrates and caused little or no disease. Its association with gastrointestinal illness in humans and animals was recognized only in the early 1980s. Over the next 25 years, information was generated on the disease's epidemiology, biology, cultivation, taxonomy and development of molecular tools. Milestones include: (i) recognition in 1980 of cryptosporidiosis as an acute enteric disease; (ii) its emergence as a chronic opportunistic infection that complicates AIDS; (iii) acknowledgement of impact on the water industry once it was shown to be waterborne; and (iv) study of Cryptosporidium genomics.

Genetic characterization of Cryptosporidium species from humans in Spain

Several species of Cryptosporidium have been associated with infection. Cryptosporidium parvum and Cryptosporidium hominis are the main agents of cryptosporidiosis in humans. Stool samples from 108 Cryptosporidium-infected patients were submitted to PCR-RFLP analysis for a 553-bp fragment of Cryptosporidium oocyst wall protein (COWP) gene and an 826-864 bp fragment of the small-subunit ribosomal RNA (SSU-rRNA) gene. Ninety-two patients were immunocompetent children and 16 were HIV-infected adults. C. hominis was detected in 69 patients (59 immunocompetent and 10 HIV-infected); C. parvum, in 34 patients (28 immunocompetent and 6 HIV-infected); and C. meleagridis and C. felis in one patient each (both immunocompetent children). Three samples yielded negative results. C. parvum was significantly more frequent in children from rural areas than in those of urban residence (p=0.010). As far as we know, this is the first surveillance study about the molecular characterization of Cryptosporidium in humans performed in Spain. The finding of zoonotic species infecting humans calls for further research on this subject.

Understanding anthelmintic resistance: The need for genomics and genetics

Anthelmintic resistance is a major problem for the control of many parasitic nematode species and has become a major constraint to livestock production in many parts of the world. In spite of its increasing importance, there is still a poor understanding of the molecular and genetic basis of resistance. It is unclear which mutations contribute most to the resistance phenotype and how resistance alleles arise, are selected and spread in parasite populations. The main strategy used to identify mutations responsible for anthelmintic resistance has been to undertake experimental studies on candidate genes. These genes have been chosen predominantly on the basis of our knowledge of drug mode-of-action and the identification of mutations that can confer resistance in model organisms. The application of these approaches to the analysis of benzimidazole and ivermectin resistance is reviewed and the reasons for their relative success or failure are discussed. The inherent limitation of candidate gene studies is that they rely on very specific and narrow assumptions about the likely identity of resistance-associated genes. In contrast, forward genetic and functional genomic approaches do not make such assumptions, as illustrated by the successful application of these techniques in the study of insecticide resistance. Although there is an urgent need to apply these powerful approaches to anthelmintic resistance research, the basic methodologies and resources are still lacking. However, these are now being developed for the trichostrongylid nematode Haemonchus contortus and the current progress and research priorities in this area are discussed.

New genotypes and factors associated with Cryptosporidium detection in mussels (Mytilus spp.) along the California coast

A 3 year study was conducted to evaluate mussels as bioindicators of faecal contamination in coastal ecosystems of California. Haemolymph samples from 4680 mussels (Mytilus spp.) were tested for Cryptosporidium genotypes using PCR amplification and DNA sequence analysis. Our hypotheses were that mussels collected from sites near livestock runoff or human sewage outflow would be more likely to contain the faecal pathogen Cryptosporidium than mussels collected distant to these sites, and that the prevalence would be greatest during the wet season when runoff into the nearshore marine environment was highest. To test these hypotheses, 156 batches of sentinel mussels were collected quarterly at nearshore marine sites considered at higher risk for exposure to livestock runoff, higher risk for exposure to human sewage, or lower risk for exposure to both faecal sources. Cryptosporidium genotypes detected in Haemolymph samples from individual mussels included Cryptosporidium parvum, Cryptosporidium felis, Cryptosporidium andersoni, and two novel Cryptosporidium spp. Factors significantly associated with detection of Cryptosporidium spp. in mussel batches were exposure to freshwater outflow and mussel collection within a week following a precipitation event. Detection of Cryptosporidium spp. was not associated with higher or lower risk status for exposure to livestock faeces or human sewage sources. This study showed that mussels can be used to monitor water quality in California and suggests that humans and animals ingesting faecal-contaminated water and shellfish may be exposed to both host-specific and anthropozoonotic Cryptosporidium genotypes of public health significance.

Integrating genetics and genomics to identify new leads for the control ofEimeriaspp

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, theEimeriaspp. present some unique phenotypes. For example, unique lines ofEimeriaspp. 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 thein ovoorin vivoselection procedures used. Of current interest is the fact that some wild-type populations ofEimeria maximaare characterized by an ability to induce protective immunity that is strain-specific. The molecular basis of phenotypes that defineEimeriaspp. 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 parasiteE. tenella. The most exciting phenotype in the context of vaccination and the development of new vaccines is the trait of strain-specific immunity associated withE. maxima. Recent work in this laboratory has shown that infection of two inbred lines of White Leghorn chickens with the W strain ofE. maximaleads 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.

An AFLP-based genetic linkage map of Plasmodium chabaudi chabaudi

Background

Plasmodium chabaudi chabaudi can be considered as a rodent model of human malaria parasites in the genetic analysis of important characters such as drug resistance and immunity. Despite the availability of some genome sequence data, an extensive genetic linkage map is needed for mapping the genes involved in certain traits.

Methods

The inheritance of 672 Amplified Fragment Length Polymorphism (AFLP) markers from two parental clones (AS and AJ) of P. c. chabaudi was determined in 28 independent recombinant progeny clones. These, AFLP markers and 42 previously mapped Restriction Fragment Length Polymorphism (RFLP) markers (used as chromosomal anchors) were organized into linkage groups using Map Manager software.

Results

614 AFLP markers formed linkage groups assigned to 10 of 14 chromosomes, and 12 other linkage groups not assigned to known chromosomes. The genetic length of the genome was estimated to be about 1676 centiMorgans (cM). The mean map unit size was estimated to be 13.7 kb/cM. This was slightly less then previous estimates for the human malaria parasite, Plasmodium falciparum

Conclusion

The P. c. chabaudi genetic linkage map presented here is the most extensive and highly resolved so far available for this species. It can be used in conjunction with the genome databases of P. c chabaudi, P. falciparum and Plasmodium yoelii to identify genes underlying important phenotypes such as drug resistance and strain-specific immunity.

CryptoDB: the Cryptosporidium genome resource

CryptoDB (http://CryptoDB.org) represents a collaborative effort to locate all genome data for the apicomplexan parasite Cryptosporidium parvum in a single user-friendly database. CryptoDB currently houses the genomic sequence data for both the human type 1 H strain and the bovine type 2 IOWA strain in addition to all other available EST and GSS sequences obtained from public repositories. All data are available for data mining via BLAST, keyword searches of pre-computed BLASTX results and user-defined or PROSITE motif pattern searches. Release 1.0 of CryptoDB contains approximately 19 million bases of genome sequence for the H and IOWA strains and an additional approximately 24 million bases of GSS and EST sequence obtained from other sources. Open reading frames greater than 50 and 100 amino acids have been generated for all sequences and all data are available for bulk download. This database, like other apicomplexan parasite databases, has been built utilizing the PlasmoDB model.

Genome Sequencing and Comparative Genomics of Tropical Disease Pathogens

The sequencing of eukaryotic genomes has lagged behind sequencing of organisms in the other domains of life, archae and bacteria, primarily due to their greater size and complexity. With recent advances in high-throughput technologies such as robotics and improved computational resources, the number of eukaryotic genome sequencing projects has increased significantly. Among these are a number of sequencing projects of tropical pathogens of medical and veterinary importance, many of which are responsible for causing widespread morbidity and mortality in peoples of developing countries. Uncovering the complete gene complement of these organisms is proving to be of immense value in the development of novel methods of parasite control, such as antiparasitic drugs and vaccines, as well as the development of new diagnostic tools. Combining pathogen genome sequences with the host and vector genome sequences is promising to be a robust method for the identification of host-pathogen interactions. Finally, comparative sequencing of related species, especially of organisms used as model systems in the study of the disease, is beginning to realize its potential in the identification of genes, and the evolutionary forces that shape the genes, that are involved in evasion of the host immune response.

Integrated mapping, chromosomal sequencing and sequence analysis of Cryptosporidium parvum

The apicomplexan Cryptosporidium parvum is one of the most prevalent protozoan parasites of humans. We report the physical mapping of the genome of the Iowa isolate, sequencing and analysis of chromosome 6, and approximately 0.9 Mbp of sequence sampled from the remainder of the genome. To construct a robust physical map, we devised a novel and general strategy, enabling accurate placement of clones regardless of clone artefacts. Analysis reveals a compact genome, unusually rich in membrane proteins. As in Plasmodium falciparum, the mean size of the predicted proteins is larger than that in other sequenced eukaryotes. We find several predicted proteins of interest as potential therapeutic targets, including one exhibiting similarity to the chloroquine resistance protein of Plasmodium. Coding sequence analysis argues against the conventional phylogenetic position of Cryptosporidium and supports an earlier suggestion that this genus arose from an early branching within the Apicomplexa. In agreement with this, we find no significant synteny and surprisingly little protein similarity with Plasmodium. Finally, we find two unusual and abundant repeats throughout the genome. Among sequenced genomes, one motif is abundant only in C. parvum, whereas the other is shared with (but has previously gone unnoticed in) all known genomes of the Coccidia and Haemosporida. These motifs appear to be unique in their structure, distribution and sequences.

Genomics and its impact on parasitology and the potential for development of new parasite control methods

Parasitic organisms remain the scourge of the developed and underdeveloped worlds. Malaria, schistosomiasis, leishmaniasis, and trypanosomiasis, for example, still result in a large number of human deaths each year worldwide, while drug resistance among nematodes still poses a major problem to the livestock industries. Genome projects involving parasitic organisms are now abundant, and technologies for the investigations of the parasite transcriptome and proteome are well established. There is no doubt the era of the "omics" is with parasitology, and current trends in the discipline are addressing fundamental biological questions that can make best use of the new technologies, as well as the vast amount of new data being generated. Will this become the "golden age of molecular parasitology," leading to the control of parasitic diseases that have plagued mankind for hundreds of years? The primary aim of this paper is to review advances in the general area of parasite genomics, and to outline where the application of "omics" technologies can and have impacted on the development of new control methods for parasitic organisms.

Polymerase chain reaction-based diagnosis of infection with Cryptosporidium in children with primary immunodeficiencies

BACKGROUND

Patients with deficient cell-mediated immunity are prone to chronic biliary tract infection with Cryptosporidium, which can lead to the development of sclerosing cholangitis and acute cryptosporidiosis after bone marrow transplantation (BMT). The organism is very difficult to detect during asymptomatic periods.

METHODS

PCR techniques were compared with standard microscopy for detecting the organism in such patients. Amplification targets were two fragments of the 18S ribosomal RNA gene (unnested) and part of the Cryptosporidium oocyst wall protein gene (nested and unnested). Twenty eight-patients with primary immunodeficiencies were studied including: CD40 ligand deficiency (13); undefined combined immunodeficiency (10); major histocompatibility complex II deficiency (2); and other defects (3). Samples analyzed included stool, bile and liver tissue.

RESULTS

Of 25 patients tested prospectively, Cryptosporidium could be detected by PCR but not by microscopy in 12, only 3 of whom had a known history of infection. Five of this group had sclerosing cholangitis. Nine of the PCR-positive patients subsequently underwent BMT and 5 developed acute posttransplant diarrhea and cholangiopathy associated with Cryptosporidium excretion. Of the 13 PCR-negative patients, 3 had cholangiopathy (sclerosing cholangitis in 1 and minor changes in 2). Four of these underwent BMT and none developed cryptosporidiosis. In 3 patients, studied only after developing post-BMT cholangiopathy and diarrhea, Cryptosporidium was detected by PCR but not by microscopy. Genotyping and sequencing showed multiple types of Cryptosporidium in approximately one-third of positive cases.

CONCLUSIONS

These results indicate that PCR-based procedures are more sensitive than microscopy for detecting Cryptosporidium in patients with immunodeficiencies.

HAPPY mapping in a plant genome: reconstruction and analysis of a high-resolution physical map of a 1.9 Mbp region of Arabidopsis thaliana chromosome 4

HAPPY mapping is an in vitro approach for defining the order and spacing of DNA markers directly on native genomic DNA. This cloning-free technique is based on analysing the segregation of markers amplified from high molecular weight genomic DNA which has been broken randomly and 'segregated' by limiting dilution into subhaploid samples. It is a uniquely versatile tool, allowing for the construction of genome maps with flexible ranges and resolutions. Moreover, it is applicable to plant genomes, for which many of the techniques pioneered in animal genomes are inapplicable or inappropriate. We report here its demonstration in a plant genome by reconstructing the physical map of a 1.9 Mbp region around the FCA locus of Arabidopsis thaliana. The resulting map, spanning around 10% of chromosome 4, is in excellent agreement with the DNA sequence and has a mean marker spacing of 16 kbp. We argue that HAPPY maps of any required resolution can be made immediately and with relatively little effort for most plant species and, furthermore, that such maps can greatly aid the construction of regional or genome-wide physical maps.

Sequence of Plasmodium falciparum chromosomes 1, 3-9 and 13

Since the sequencing of the first two chromosomes of the malaria parasite, Plasmodium falciparum, there has been a concerted effort to sequence and assemble the entire genome of this organism. Here we report the sequence of chromosomes 1, 3-9 and 13 of P. falciparum clone 3D7--these chromosomes account for approximately 55% of the total genome. We describe the methods used to map, sequence and annotate these chromosomes. By comparing our assemblies with the optical map, we indicate the completeness of the resulting sequence. During annotation, we assign Gene Ontology terms to the predicted gene products, and observe clustering of some malaria-specific terms to specific chromosomes. We identify a highly conserved sequence element found in the intergenic region of internal var genes that is not associated with their telomeric counterparts.

Synteny comparison between apes and human using fine-mapping of the genome

Comparing the genomes of the great apes and human should provide novel information concerning the origins of humankind. Relative to the great apes, the human karyotype has one fewer chromosome pair, as human chromosome 2 derived from the telomeric fusion of two ancestral primate chromosomes. To identify the genomic rearrangements that accompanied human speciation, we initiated a comparative study between human, chimpanzee, and gorilla. Using the HAPPY mapping method, an acellular adaptation of the radiation hybrid method, we mapped a few hundred markers on the human, chimpanzee, and gorilla genomes. This allowed us to identify several chromosome rearrangements, in particular a pericentric inversion and a translocation. We precisely localized the synteny breakpoint that led to the formation of human chromosome 2. This breakpoint was confirmed by FISH mapping.

Detection and genotyping of oocysts of Cryptosporidium parvum by real-time PCR and melting curve analysis

Several real-time PCR procedures for the detection and genotyping of oocysts of Cryptosporidium parvum were evaluated. A 40-cycle amplification of a 157-bp fragment from the C. parvum beta-tubulin gene detected individual oocysts which were introduced into the reaction mixture by micromanipulation. SYBR Green I melting curve analysis was used to confirm the specificity of the method when DNA extracted from fecal samples spiked with oocysts was analyzed. Because C. parvum isolates infecting humans comprise two distinct genotypes, designated type 1 and type 2, real-time PCR methods for discriminating C. parvum genotypes were developed. The first method used the same beta-tubulin amplification primers and two fluorescently labeled antisense oligonucleotide probes spanning a 49-bp polymorphic sequence diagnostic for C. parvum type 1 and type 2. The second genotyping method used SYBR Green I fluorescence and targeted a polymorphic coding region within the GP900/poly(T) gene. Both methods discriminated between type 1 and type 2 C. parvum on the basis of melting curve analysis. To our knowledge, this is the first report describing the application of melting curve analysis for genotyping of C. parvum oocysts.

Comparison of fluorescence, antigen and PCR assays to detect Cryptosporidium parvum in fecal specimens

To optimize routine screening for cryptosporidiosis, 198 stool samples from patients at risk and from calves were examined by enzyme immunoassay (EIA), a direct fluorescent-antibody (DFA) and a modified immunofluorescence assay. Ninety-nine samples were positive in at least one assay, whereas 99 were negative in all three assays. Sensitivity of antigen EIA and DFA were similar (94%, 95% CI: 88-98%, and 91%, 95% CI: 84-95%). The modified immunofluorescence was significantly less sensitive (64%, 95% CI: 55-74%). 149 samples were also examined by two nested PCR assays targeting either the 18S rRNA or Cryptosporidium outer wall protein (COWP) gene. A PCR product was amplified from 86 out of 89 samples being positive in at least one other assay (sensitivity 97%, 95% CI: 91-99%). None was obtained from 60 samples negative in the three other assays. PCR assays did not increase the detection rate. Antigen EIA or DFA appear sufficient for routine Cryptosporidium screening of fecal samples.

Genetic complementation in apicomplexan parasites

A robust forward genetic model for Apicomplexa could greatly enhance functional analysis of genes in these important protozoan pathogens. We have developed and successfully tested a genetic complementation strategy based on genomic insertion in Toxoplasma gondii. Adapting recombination cloning to genomic DNA, we show that complementing sequences can be shuttled between parasite genome and bacterial plasmid, providing an efficient tool for the recovery and functional assessment of candidate genes. We show complementation, gene cloning, and biological verification with a mutant parasite lacking hypoxanthine-xanthine-guanine phosphoribosyltransferase and a T. gondii cDNA library. We also explored the utility of this approach to clone genes based on function from other apicomplexan parasites using Toxoplasma as a surrogate. A heterologous library containing Cryptosporidium parvum genomic DNA was generated, and we identified a C. parvum gene coding for inosine 5-monophosphate-dehydrogenase (IMPDH). Interestingly, phylogenetic analysis demonstrates a clear eubacterial origin of this gene and strongly suggests its lateral transfer from epsilon-proteobacteria. The prokaryotic origin of this enzyme might make it a promising target for therapeutics directed against Cryptosporidium.

Characterisation of global protein expression by two-dimensional electrophoresis and mass spectrometry: proteomics of Toxoplasma gondii

The development of tools for the analysis of global gene expression is vital for the optimal exploitation of the data on parasite genomes that are now being generated in abundance. Recent advances in two-dimensional electrophoresis (2-DE), mass spectrometry and bioinformatics have greatly enhanced the possibilities for mapping and characterisation of protein populations. We have employed these developments in a proteomics approach for the analysis of proteins expressed in the tachyzoite stage of Toxoplasma gondii. Over 1000 polypeptides were reproducibly separated by high-resolution 2-DE using the pH ranges 4-7 and 6-11. Further separations using narrow range gels suggest that at least 3000-4000 polypeptides should be resolvable by 2-DE using multiple single pH unit gels. Mass spectrometry was used to characterise a variety of protein spots on the 2-DE gels. Peptide mass fingerprints, acquired by matrix-assisted laser desorption/ionisation-(MALDI) mass spectrometry, enabled unambiguous protein identifications to be made where full gene sequence information was available. However, interpretation of peptide mass fingerprint data using the T. gondii expressed sequence tag (EST) database was less reliable. Peptide fragmentation data, acquired by post-source decay mass spectrometry, proved a more successful strategy for the putative identification of proteins using the T. gondii EST database and protein databases from other organisms. In some instances, several protein spots appeared to be encoded by the same gene, indicating that post-translational modification and/or alternative splicing events may be a common feature of functional gene expression in T. gondii. The data demonstrate that proteomic analyses are now viable for T. gondii and other protozoa for which there are good EST databases, even in the absence of complete genome sequence. Moreover, proteomics is of great value in interpreting and annotating EST databases.

A status report on the sequencing and annotation of the P. falciparum genome

Almost 5 years ago, an international consortium of sequencing centers and funding agencies was formed to sequence the genome of the human malaria parasite Plasmodium falciparum. A novel chromosome by chromosome shotgun strategy was devised to sequence this very AT-rich genome. Two of the 14 chromosomes have been completed and the remaining chromosomes are in the final stages of gap closure. The consortium recently developed plans for the annotation and analysis of the complete genome sequence and its publication in 2002.

A high-resolution HAPPY map of Dictyostelium discoideum chromosome 6

We have made a high-resolution HAPPY map of chromosome 6 of Dictyostelium discoideum consisting of 300 sequence-tagged sites with an average spacing of 14 kb along the approximately 4-Mb chromosome. The majority of the marker sequences were derived from randomly chosen clones from four different chromosome 6-enriched plasmid libraries or from subclones of YACs previously mapped to chromosome 6. The map appears to span the entire chromosome, although marker density is greater in some regions than in others and is lowest within the telomeric region. Our map largely supports previous gene-based maps of this chromosome but reveals a number of errors in the physical map. In addition, we find that a high proportion of the plasmid sequences derived from gel-enriched chromosome 6 (and that form the basis of a chromosome-specific sequencing project) originates from other chromosomes.

The initiation translation factor eIF-4A of Cryptosporidium parvum is encoded by two distinct mRNA forms and shows DNA sequence polymorphism distinguishing genotype 1 and 2 isolates

The eukaryotic translation initiation factor eIF-4A is an ATP-dependent RNA helicase involved in ribosome attachment to the 5' end of mRNAs. Employing as a probe a Cryptosporidium parvum genomic amplicon encoding a partial polypeptide related to eIF-4A, we screened a C. parvum sporozoite cDNA library to clone the full length of the gene. Two complete cDNAs were characterized, Cp.F6 and Cp.F10, which consisted of 1,900 and 1,418 bp, respectively. The overlapping portions of the sequences shared 100% identity and encoded a polypeptide of 405 amino acids whose identity to known eIF-4A molecules ranged between 77 and 39%. The 2 cDNAs differed in the length of their respective 3' untranslated regions, of 577 bp in Cp.F6 and 72 bp in Cp.F10, in both of which a putative polyadenylation signal was identified. The structure of the cloned cDNAs, along with genomic Southern blot data indicating that eIF-4A is encoded by a single copy gene, strongly suggested that Cp.F6 and Cp.F10 reflect a differential 3' end processing of mRNA precursors, not observed so far in C. parvum. Northern blot analysis confirmed that the sporozoites express 2 eIF-4A mRNAs and showed that the lower molecular weight transcript is 10- to 20-fold more abundant. We also investigated the polymorphism of the eIF-4A gene and defined a novel polymerase chain reaction-restriction fragment length polymorphism marker discriminating between C. parvum isolates of genotypes 1 and 2.

Cryptosporidium parvum: the many secrets of a small genome

The coccidium Cryptosporidium parvum is an obligate intracellular parasite of the phylum Apicomplexa. It infects the gastrointestinal tract of humans and livestock, and represents the third major cause of diarrhoeal disease worldwide. Scarcely considered for decades due to its apparently non-pathogenic nature, C. parvum has been studied very actively over the last 15 years, after its medical relevance as a dangerous opportunistic parasite and widespread water contaminant was fully recognised. Despite the lack of an efficient in vitro culture system and appropriate animal models, significant advances have been made in this relatively short period of time towards understanding C. parvum biology, immunology, genetics and epidemiology. Until recently, very little was known about the genome of C. parvum, with even basic issues, such as the number and size of chromosomes, being the object of a certain controversy. With the advent of pulsed field gradient electrophoresis and the introduction of molecular biology techniques, the overall structure and fine organisation of the genome of C. parvum have started to be disclosed. Organised into eight chromosomes distributed in a very narrow range of molecular masses, the genome of C. parvum is one of the smallest so far described among unicellular eukaryotic organisms. Although fewer than 30 C. parvum genes have been cloned so far, information about the overall structure of the parasite genome has increased exponentially over the last 2 years. From the first karyotypic analyses to the recent development of physical maps for individual chromosomes, this review will try to describe the state-of-the-art of our knowledge on the nuclear genome of C. parvum and will discuss the available experimental evidence concerning the presence of extra-chromosomal elements.

Genetic characterization of Cryptosporidium strains from 218 patients with diarrhea diagnosed as having sporadic cryptosporidiosis

Samples of whole feces in which Cryptosporidium oocysts were recognized by hospital laboratories were collected from 218 patients with diarrhea. All samples were reexamined by light microscopy, and oocysts were detected in 211 samples. A simple and rapid procedure for the extraction of DNA from whole feces was developed, and this was used to amplify fragments of the Cryptosporidium outer wall protein (COWP), the thrombospondin-related adhesive protein C1 (TRAP-C1), and the 18S rRNA genes by PCR. For seven samples oocysts were not detected by microscopy and DNA failed to be amplified by the three PCR procedures. Among the 211 samples "positive" by microscopy, the sensitivities of PCRs for the 18S rRNA, COWP, and TRAP-C1 gene fragments were 97, 91, and 66%, respectively. The sensitivities of all three PCR procedures increased with increasing numbers of oocysts as observed by microscopy. Two genotypes of the COWP and TRAP-C1 genes can be detected by PCR-restriction fragment length polymorphism analysis. With this series of samples, the same genotypes of the COWP and TRAP-C1 genes always segregated together. A combined genotyping data set was produced for isolates from 194 samples: 74 (38%) were genotype 1 and 120 (62%) were genotype 2. Genotype 2 was detected in a significantly greater proportion of the samples with small numbers of oocysts, and genotype 1 was detected in a significantly greater proportion of the samples with larger numbers of oocysts. There were no significant differences in the distribution of the genotypes by patient sex and age. The distribution of the genotypes was significantly different both in patients with a history of foreign travel and in those from different regions in England.

A random survey of the Cryptosporidium parvum genome

Cryptosporidium parvum is an obligate intracellular pathogen responsible for widespread infections in humans and animals. The inability to obtain purified samples of this organism's various developmental stages has limited the understanding of the biochemical mechanisms important for C. parvum development or host-parasite interaction. To identify C. parvum genes independent of their developmental expression, a random sequence analysis of the 10.4-megabase genome of C. parvum was undertaken. Total genomic DNA was sheared by nebulization, and fragments between 800 and 1,500 bp were gel purified and cloned into a plasmid vector. A total of 442 clones were randomly selected and subjected to automated sequencing by using one or two primers flanking the cloning site. In this way, 654 genomic survey sequences (GSSs) were generated, corresponding to >320 kb of genomic sequence. These sequences were assembled into 408 contigs containing >250 kb of unique sequence, representing approximately 2.5% of the C. parvum genome. Comparison of the GSSs with sequences in the public DNA and protein databases revealed that 107 contigs (26%) displayed similarity to previously identified proteins and rRNA and tRNA genes. These included putative genes involved in the glycolytic pathway, DNA, RNA, and protein metabolism, and signal transduction pathways. The repetitive sequence elements identified included a telomere-like sequence containing hexamer repeats, 57 microsatellite-like elements composed of dinucleotide or trinucleotide repeats, and a direct repeat sequence. This study demonstrates that large-scale genomic sequencing is an efficient approach to analyze the organizational characteristics and information content of the C. parvum genome.

Speculation on whether a vaccine against cryptosporidiosis is a reality or fantasy

In this paper the authors question whether the development of a vaccine against cryptosporidiosis could be taken into consideration. The necessity and feasibility of such a vaccine for human and veterinary application is discussed. Developmental stages within the life cycle of the parasite that might act as possible targets for vaccine development are summarised, as well as the target antigens offered by molecular biology and immunology studies. Vaccination trials against cryptosporidiosis carried out so far, including the active and passive immunisation approach, are also overviewed. It seems that with respect to a Cryptosporidium vaccine two target groups can be considered: children of the developing world and neonatal ruminants. Antigens representing possible candidates for a subunit vaccine were identified based on their function, location and/or the immune response they evoke. While the active vaccination of newborn calves, lambs and goat kids has to face a number of important limitations, the passive immunisation approach, where dams were immunised to protect their progeny by colostral transfer, was proven to be a valuable alternative. Finally, a number of points of action for the near future are put forward.

beta-tubulin mRNA as a marker of Cryptosporidium parvum oocyst viability

Determining the viability of waterborne Cryptosporidium parvum oocysts remains a technical challenge. rRNA and mRNA were evaluated in a reverse transcription (RT)-PCR assay as potential markers of oocyst viability. The rationale for this approach is the rapid turnover and postmortem decay of cellular RNA. The beta-tubulin mRNA and an anonymous mRNA transcript were chosen as potential markers because they are the only mRNA species in C. parvum known to possess introns. This feature facilitated the distinction between genuine RT-PCR products and PCR products originating from copurifying DNA. Prolonged incubation at room temperature of initially viable oocysts resulted in a gradual decrease in mRNA levels, which correlated with the loss of oocyst infectivity to neonatal mice. In contrast, oocysts stored at 4 degrees C for over 39 weeks maintained their infectivity and displayed no decrease in the level of beta-tubulin RT-PCR product. The postmortem decay of two mRNA species demonstrates that RT-PCR analysis can provide information on the viability of C. parvum oocysts. The methodological similarity between PCR detection and RT-PCR viability analysis could facilitate the development of a combined detection and viability assay.