Integrated mapping, chromosomal sequencing and sequence analysis of Cryptosporidium parvum

Paving the Way: Contributions of Big Data to Apicomplexan and Kinetoplastid Research

In the age of big data an important question is how to ensure we make the most out of the resources we generate. In this review, we discuss the major methods used in Apicomplexan and Kinetoplastid research to produce big datasets and advance our understanding of Plasmodium, Toxoplasma, Cryptosporidium, Trypanosoma and Leishmania biology. We debate the benefits and limitations of the current technologies, and propose future advancements that may be key to improving our use of these techniques. Finally, we consider the difficulties the field faces when trying to make the most of the abundance of data that has already been, and will continue to be, generated.

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.

A review of Cryptosporidium spp. and their detection in water

Cryptosporidium spp. are one of the most important waterborne pathogens worldwide and a leading cause of mortality from waterborne gastrointestinal diseases. Detection of Cryptosporidium spp. in water can be very challenging due to their low numbers and the complexity of the water matrix. This review describes the biology of Cryptosporidium spp. and current methods used in their detection with a focus on C. parvum and C. hominis. Among the methods discussed and compared are microscopy, immunology-based methods using monoclonal antibodies, molecular methods including PCR (polymerase chain reaction)-based assays, and emerging aptamer-based methods. These methods have different capabilities and limitations, but one common challenge is the need for better sensitivity and specificity, particularly in the presence of contaminants. The application of DNA aptamers in the detection of Cryptosporidium spp. oocysts shows promise in overcoming these challenges, and there will likely be significant developments in aptamer-based sensors in the near future.

Development of novel methodology for the molecular differentiation of Cryptosporidium parvum gp60 subtypes via high resolution melting analysis

Cryptosporidium species subtypes are generally identified via DNA sequencing of the gp60 gene tandem repeat motif region. Due to the immunogenic nature of its glycoprotein products, gp60 is subject to host selective pressures, genetic recombination and evolutionary processes that drive extensive polymorphism at this locus. The elucidation of the polymorphic nature of this gene has led to the current mainstay in Cryptosporidium subtyping nomenclature.

This study aimed to develop a real-time polymerase chain reaction based method utilising a post-PCR application, high resolution melting (HRM) analysis, in conjunction with the abovementioned gp60 nomenclature system, in order to differentiate between Cryptosporidium parvum gp60 subtypes. Subtype differentiation is based on the difference between the melting temperatures of individual subtypes conferred by variations in the polymorphic region of gp60.

• Nested gp60 primers were designed to amplify a target region of <200 base pairs for effective HRM analysis

• This method presents a rapid, sensitive, cost effective alternative to conventional sequencing.

• This method is highly flexible and may be applied to other loci in order to facilitate multi-locus analysis and improve the discriminative abilities of the method.

Comparative analyses and structural insights of new class glutathione transferases in Cryptosporidium species

Cryptosporidiosis, caused by protozoan parasites of the genus Cryptosporidium, is estimated to rank as a leading cause in the global burden of neglected zoonotic parasitic diseases. This diarrheal disease is the second leading cause of death in children under 5 years of age. Based on the C. parvum transcriptome data, glutathione transferase (GST) has been suggested as a drug target against this pathogen. GSTs are diverse multifunctional proteins involved in cellular defense and detoxification in organisms and help pathogens to alleviate chemical and environmental stress. In this study, we performed genome-wide data mining, identification, classification and in silico structural analysis of GSTs in fifteen Cryptosporidium species. The study revealed the presence three GSTs in each of the Cryptosporidium species analyzed in the study. Based on the percentage identity and comprehensive comparative phylogenetic analysis, we assigned Cryptosporidium species GSTs to three new GST classes, named Vega (ϑ), Gamma (γ) and Psi (ψ). The study also revealed an atypical thioredoxin-like fold in the C. parvum GST1 of the Vega class, whereas C. parvum GST2 of the Gamma class and C. melagridis GST3 of the Psi class has a typical thioredoxin-like fold in the N-terminal region. This study reports the first comparative analysis of GSTs in Cryptosporidium species.

Divergent Copies of a Cryptosporidium parvum-Specific Subtelomeric Gene

Subtype families of Cryptosporidium parvum differ in host range, with IIa and IId being found in a broad range of animals, IIc in humans, and IIo and IIp in some rodents. Previous studies indicated that the subtelomeric cgd6_5520-5510 gene in C. parvum is lost in many Cryptosporidium species, and could potentially contribute to the broad host range of the former. In this study, we identified the presence of a second copy of the gene in some C. parvum subtype families with a broad host range, and showed sequence differences among them. The sequence differences in the cgd6_5520-5510 gene were not segregated by the sequence type of the 60 kDa glycoprotein gene. Genetic recombination appeared to have played a role in generating divergent nucleotide sequences between copies and among subtype families. These data support the previous conclusion on the potential involvement of the insulinase-like protease encoded by the subtelomeric cgd6_5520-5510 gene in the broad host range of C. parvum IIa and IId subtypes.

ApiAP2 Transcription Factors in Apicomplexan Parasites

Apicomplexan parasites are protozoan organisms that are characterised by complex life cycles and they include medically important species, such as the malaria parasite Plasmodium and the causative agents of toxoplasmosis (Toxoplasma gondii) and cryptosporidiosis (Cryptosporidium spp.). Apicomplexan parasites can infect one or more hosts, in which they differentiate into several morphologically and metabolically distinct life cycle stages. These developmental transitions rely on changes in gene expression. In the last few years, the important roles of different members of the ApiAP2 transcription factor family in regulating life cycle transitions and other aspects of parasite biology have become apparent. Here, we review recent progress in our understanding of the different members of the ApiAP2 transcription factor family in apicomplexan parasites.

Comparison of current methods used to detect Cryptosporidium oocysts in stools

In this review all of the methods that are currently in use for the investigation of Cryptosporidium in stool material are highlighted and critically discussed. It appears that more qualifications and background knowledge in this field regarding the diagnosis of the Cryptosporidium parasite is required. Furthermore, there is no standardization for the protocols that are commonly used to either detect oocysts in faeces or to diagnose the Cryptosporidium infection. It is therefore necessary to initiate further education and research that will assist in improving the accuracy of the diagnosis of Cryptosporidium oocysts in the faecal micro-cosmos. Where ambient concentrations of oocysts are low in stool material, detection becomes a formidable task. Procedures for ring tests and the standardization of multi-laboratory testing are recommended. It is also necessary to enhance the routine surveillance capacity of cryptosporidiosis and to improve the safety against it, considering the fact that this disease is under diagnosed and under reported. This review is intended to stimulate research that could lead to future improvements and further developments in monitoring the diagnostic methodologies that will assist in harmonizing Cryptosporidium oocysts in stool diagnosis.

An ApiAP2 member regulates expression of clonally variant genes of the human malaria parasite Plasmodium falciparum

Variegated surface antigen expression is key to chronic infection and pathogenesis of the human malaria parasite Plasmodium falciparum. This protozoan parasite expresses distinct surface molecules that are encoded by clonally variant gene families such as var, rif and stevor. The molecular mechanisms governing activation of individual members remain ill-defined. To investigate the molecular events of the initial transcriptional activation process we focused on a member of the apicomplexan ApiAP2 transcription factor family predicted to bind to the 5′ upstream regions of the var gene family, AP2-exp (PF3D7_1466400). Viable AP2-exp mutant parasites rely on expressing no less than a short truncated protein including the N-terminal AP2 DNA-binding domain. RNA-seq analysis in mutant parasites revealed transcriptional changes in a subset of exported proteins encoded by clonally variant gene families. Upregulation of RIFINs and STEVORs was validated at the protein levels. In addition, morphological alterations were observed on the surface of the host cells infected by the mutants. This work points to a complex regulatory network of clonally variant gene families in which transcription of a subset of members is regulated by the same transcription factor. In addition, we highlight the importance of the non-DNA binding AP2 domain in functional gene regulation.

Genome-wide upstream motif analysis of Cryptosporidium parvum genes clustered by expression profile

Background

There are very few molecular genetic tools available to study the apicomplexan parasite Cryptosporidium parvum. The organism is not amenable to continuous in vitro cultivation or transfection, and purification of intracellular developmental stages in sufficient numbers for most downstream molecular applications is difficult and expensive since animal hosts are required. As such, very little is known about gene regulation in C. parvum.

Results

We have clustered whole-genome gene expression profiles generated from a previous study of seven post-infection time points of 3,281 genes to identify genes that show similar expression patterns throughout the first 72 hours of in vitro epithelial cell culture. We used the algorithms MEME, AlignACE and FIRE to identify conserved, overrepresented DNA motifs in the upstream promoter region of genes with similar expression profiles. The most overrepresented motifs were E2F (5′-TGGCGCCA-3′); G-box (5′-G.GGGG-3′); a well-documented ApiAP2 binding motif (5′-TGCAT-3′), and an unknown motif (5′-[A/C] AACTA-3′). We generated a recombinant C. parvum DNA-binding protein domain from a putative ApiAP2 transcription factor [CryptoDB: cgd8_810] and determined its binding specificity using protein-binding microarrays. We demonstrate that cgd8_810 can putatively bind the overrepresented G-box motif, implicating this ApiAP2 in the regulation of many gene clusters.

Conclusion

Several DNA motifs were identified in the upstream sequences of gene clusters that might serve as potential cis-regulatory elements. These motifs, in concert with protein DNA binding site data, establish for the first time the beginnings of a global C. parvum gene regulatory map that will contribute to our understanding of the development of this zoonotic parasite.

Library of Apicomplexan Metabolic Pathways: a manually curated database for metabolic pathways of apicomplexan parasites

The Library of Apicomplexan Metabolic Pathways (LAMP, http://www.llamp.net) is a web database that provides near complete mapping from genes to the central metabolic functions for some of the prominent intracellular parasites of the phylum Apicomplexa. This phylum includes the causative agents of malaria, toxoplasmosis and theileriosis-diseases with a huge economic and social impact. A number of apicomplexan genomes have been sequenced, but the accurate annotation of gene function remains challenging. We have adopted an approach called metabolic reconstruction, in which genes are systematically assigned to functions within pathways/networks for Toxoplasma gondii, Neospora caninum, Cryptosporidium and Theileria species, and Babesia bovis. Several functions missing from pathways have been identified, where the corresponding gene for an essential process appears to be absent from the current genome annotation. For each species, LAMP contains interactive diagrams of each pathway, hyperlinked to external resources and annotated with detailed information, including the sources of evidence used. We have also developed a section to highlight the overall metabolic capabilities of each species, such as the ability to synthesize or the dependence on the host for a particular metabolite. We expect this new database will become a valuable resource for fundamental and applied research on the Apicomplexa.

Insights into the genome structure and copy-number variation of Eimeria tenella

BACKGROUND

Eimeria is a genus of parasites in the same phylum (Apicomplexa) as human parasites such as Toxoplasma, Cryptosporidium and the malaria parasite Plasmodium. As an apicomplexan whose life-cycle involves a single host, Eimeria is a convenient model for understanding this group of organisms. Although the genomes of the Apicomplexa are diverse, that of Eimeria is unique in being composed of large alternating blocks of sequence with very different characteristics - an arrangement seen in no other organism. This arrangement has impeded efforts to fully sequence the genome of Eimeria, which remains the last of the major apicomplexans to be fully analyzed. In order to increase the value of the genome sequence data and aid in the effort to gain a better understanding of the Eimeria tenella genome, we constructed a whole genome map for the parasite.

RESULTS

A total of 1245 contigs representing 70.0% of the whole genome assembly sequences (Wellcome Trust Sanger Institute) were selected and subjected to marker selection. Subsequently, 2482 HAPPY markers were developed and typed. Of these, 795 were considered as usable markers, and utilized in the construction of a HAPPY map. Markers developed from chromosomally-assigned genes were then integrated into the HAPPY map and this aided the assignment of a number of linkage groups to their respective chromosomes. BAC-end sequences and contigs from whole genome sequencing were also integrated to improve and validate the HAPPY map. This resulted in an integrated HAPPY map consisting of 60 linkage groups that covers approximately half of the estimated 60 Mb genome. Further analysis suggests that the segmental organization first seen in Chromosome 1 is present throughout the genome, with repeat-poor (P) regions alternating with repeat-rich (R) regions. Evidence of copy-number variation between strains was also uncovered.

CONCLUSIONS

This paper describes the application of a whole genome mapping method to improve the assembly of the genome of E. tenella from shotgun data, and to help reveal its overall structure. A preliminary assessment of copy-number variation (extra or missing copies of genomic segments) between strains of E. tenella was also carried out. The emerging picture is of a very unusual genome architecture displaying inter-strain copy-number variation. We suggest that these features may be related to the known ability of this parasite to rapidly develop drug resistance.

Assessment of polymorphic genetic markers for multi-locus typing of Cryptosporidium parvum and Cryptosporidium hominis

The use of high resolution molecular tools to study Cryptosporidium parvum and Cryptosporidium hominis intra-species variation is becoming common practice, but there is currently no consensus in the methods used. The most commonly applied tool is partial gp60 gene sequence analysis. However, multi-locus schemes are acknowledged to improve resolution over analysis of a single locus, which neglects potential re-assortment of genes during the sexual phase of the Cryptosporidium life-cycle. Multi-locus markers have been investigated in isolates from a variety of sampling frames, in varying combinations and using different assays and methods of analysis. To identify the most informative markers as candidates for the development of a standardised multi-locus fragment size-based typing (MLFT) scheme to integrate with epidemiological analyses, we examined the published literature. A total of 31 MLFT studies were found, employing 55 markers of which 45 were applied to both C. parvum and C. hominis. Of the studies, 11 had sufficient raw data, from three or more markers, and a sampling frame containing at least 50 samples, for meaningful in-depth analysis using assessment criteria based on the sampling frame, study size, number of markers investigated in each study, marker characteristics (>2 nucleotide repeats) and the combinations of markers generating all possible multi-locus genotypes. Markers investigated differed between C. hominis and C. parvum. When each scheme was analysed for the fewest markers required to identify 95% of all MLFTs, some redundancy was identified in all schemes; an average redundancy of 40% for C. hominis and 27% for C. parvum. Ranking markers, based on the most productive combinations, identified two different sets of potentially most informative candidate markers, one for each species. These will be subjected to technical evaluation including typability (percentage of samples generating a complete multi-locus type) and discriminatory power by direct fragment size analysis and analysed for correlation with epidemiological data in suitable sampling frames. The establishment of a group of users and agreed subtyping scheme for improved epidemiological and public health investigations of C. parvum and C. hominis will facilitate further developments and consideration of technological advances in a harmonised manner.

Identification of putative potassium channel homologues in pathogenic protozoa

K⁺ channels play a vital homeostatic role in cells and abnormal activity of these channels can dramatically alter cell function and survival, suggesting that they might be attractive drug targets in pathogenic organisms. Pathogenic protozoa lead to diseases such as malaria, leishmaniasis, trypanosomiasis and dysentery that are responsible for millions of deaths each year worldwide. The genomes of many protozoan parasites have recently been sequenced, allowing rational design of targeted therapies. We analyzed the genomes of pathogenic protozoa and show the existence within them of genes encoding putative homologues of K⁺ channels. These protozoan K⁺ channel homologues represent novel targets for anti-parasitic drugs. Differences in the sequences and diversity of human and parasite proteins may allow pathogen-specific targeting of these K⁺ channel homologues.

Genomics and population biology of Cryptosporidium species

We describe recent advances in the genomics and population biology of Cryptosporidium parvum and C. hominis, the causative agents of cryptosporidiosis in humans and animals. Many basic aspects of the biology of Cryptosporidium species remain to be investigated and effective drugs to control cryptosporidiosis are not available. Sequencing and annotation of the genome of C. parvum and C. hominis has uncovered unique features of the metabolism of these species. The recently sequenced genome of the gastric species C. muris is providing new insights into the evolution of the genus. Cryptosporidian sequence information has facilitated the identification of polymorphic genetic markers. Genotyping of oocysts excreted by human and animal hosts using such markers has revealed many new species and genotypes, and is leading to a better understanding of the epidemiology of cryptosporidiosis.

Identification of intracellular and plasma membrane calcium channel homologues in pathogenic parasites

Ca²⁺ channels regulate many crucial processes within cells and their abnormal activity can be damaging to cell survival, suggesting that they might represent attractive therapeutic targets in pathogenic organisms. Parasitic diseases such as malaria, leishmaniasis, trypanosomiasis and schistosomiasis are responsible for millions of deaths each year worldwide. The genomes of many pathogenic parasites have recently been sequenced, opening the way for rational design of targeted therapies. We analyzed genomes of pathogenic protozoan parasites as well as the genome of Schistosoma mansoni, and show the existence within them of genes encoding homologues of mammalian intracellular Ca²⁺ release channels: inositol 1,4,5-trisphosphate receptors (IP₃Rs), ryanodine receptors (RyRs), two-pore Ca²⁺ channels (TPCs) and intracellular transient receptor potential (Trp) channels. The genomes of Trypanosoma, Leishmania and S. mansoni parasites encode IP₃R/RyR and Trp channel homologues, and that of S. mansoni additionally encodes a TPC homologue. In contrast, apicomplexan parasites lack genes encoding IP₃R/RyR homologues and possess only genes encoding TPC and Trp channel homologues (Toxoplasma gondii) or Trp channel homologues alone. The genomes of parasites also encode homologues of mammalian Ca²⁺influx channels, including voltage-gated Ca²⁺ channels and plasma membrane Trp channels. The genome of S. mansoni also encodes Orai Ca²⁺ channel and STIM Ca²⁺ sensor homologues, suggesting that store-operated Ca²⁺ entry may occur in this parasite. Many anti-parasitic agents alter parasite Ca²⁺ homeostasis and some are known modulators of mammalian Ca²⁺ channels, suggesting that parasite Ca²⁺ channel homologues might be the targets of some current anti-parasitic drugs. Differences between human and parasite Ca²⁺ channels suggest that pathogen-specific targeting of these channels may be an attractive therapeutic prospect.

Structure and RNA recognition by the snRNA and snoRNA transport factor PHAX

Small nuclear and small nucleolar RNAs (snRNAs and snoRNAs) are critical components of snRNPs and snoRNPs and play an essential role in the maturation of, respectively, mRNAs and rRNAs within the nucleus of eukaryotic cells. Complex and specific pathways exist for the assembly of snRNPs and snoRNPs, involving, for instance, nucleocytoplasmic transport of snRNAs and intranuclear transport between compartments of snoRNAs. The phosphorylated adaptor for nuclear export (PHAX) is required for nuclear export of snRNAs in metazoans and also involved in the intranuclear transport of snoRNAs to Cajal bodies. PHAX contains a conserved single-stranded nucleic acid binding domain (RNA_GG_bind domain) with no sequence homology with any other known RNA-binding module. Here, we report NMR and X-ray crystallography studies that elucidate the structural basis for RNA recognition by the PHAX RNA-binding domain (PHAX-RBD). The crystal structure of the RNA_GG_bind domain from the parasite Cryptosporidium parvum (Cp RBD) forms well-folded dimers in solution in the absence of any ligand. The human PHAX-RBD is monomeric and only adopts a tertiary fold upon RNA binding. The PHAX-RBD represents a novel helical fold and binds single-stranded RNA with micromolar affinity without sequence specificity. RNA recognition by human PHAX-RBD is consistent with mutational analysis that affects RNA binding and PHAX-mediated nuclear export. Our data suggest that the PHAX-RBD mediates auxiliary RNA contacts with the snRNA and snoRNA substrates that are required for transport and/or substrate release.

Conservation and divergence of known apicomplexan transcriptional regulons

BACKGROUND

The apicomplexans are a diverse phylum of parasites causing an assortment of diseases including malaria in a wide variety of animals and lymphoproliferation in cattle. Little is known about how these varied parasites regulate their transcriptional regulons. Even less is known about how regulon systems, consisting of transcription factors and target genes together with their associated biological process, evolve in these diverse parasites.

RESULTS

In order to obtain insights into the differences in transcriptional regulation between these parasites we compared the orthology profiles of putative malaria transcription factors across species and examined the enrichment patterns of four binding sites across eleven apicomplexans. About three-fifths of the factors are broadly conserved in several phylogenetic orders of sequenced apicomplexans. This observation suggests the existence of regulons whose regulation is conserved across this ancient phylum. Transcription factors not broadly conserved across the phylum are possibly involved in regulon systems that have diverged between species. Examining binding site enrichment patterns in light of transcription factor conservation patterns suggests a second mode via which regulon systems may diverge - rewiring of existing transcription factors and their associated binding sites in specific ways. Integrating binding sites with transcription factor conservation patterns also facilitated prediction of putative regulators for one of the binding sites.

CONCLUSIONS

Even though transcription factors are underrepresented in apicomplexans, the distribution of these factors and their associated regulons reflect common and family-specific transcriptional regulatory processes.

Identification and functional characterization of cis-regulatory elements in the apicomplexan parasite Toxoplasma gondii

BACKGROUND

Toxoplasma gondii is a member of the phylum Apicomplexa, which consists entirely of parasitic organisms that cause several diseases of veterinary and human importance. Fundamental mechanisms of gene regulation in this group of protistan parasites remain largely uncharacterized. Owing to their medical and veterinary importance, genome sequences are available for several apicomplexan parasites. Their genome sequences reveal an apparent paucity of known transcription factors and the absence of canonical cis-regulatory elements. We have approached the question of gene regulation from a sequence perspective by mining the genomic sequence data to identify putative cis-regulatory elements using a de novo approach.

RESULTS

We have identified putative cis-regulatory elements present upstream of functionally related groups of genes and subsequently characterized the function of some of these conserved elements using reporter assays in the parasite. We show a sequence-specific role in gene-expression for seven out of eight identified elements.

CONCLUSIONS

This work demonstrates the power of pure sequence analysis in the absence of expression data or a priori knowledge of regulatory elements in eukaryotic organisms with compact genomes.

Role of CpSUB1, a subtilisin-like protease, in Cryptosporidium parvum infection in vitro

The apicomplexan parasite Cryptosporidium is a significant cause of diarrheal disease worldwide. Previously, we reported that a Cryptosporidium parvum subtilisin-like serine protease activity with furin-type specificity cleaves gp40/15, a glycoprotein that is proteolytically processed into gp40 and gp15, which are implicated in mediating infection of host cells. Neither the enzyme(s) responsible for the protease activity in C. parvum lysates nor those that process gp40/15 are known. There are no furin or other proprotein convertase genes in the C. parvum genome. However, a gene encoding CpSUB1, a subtilisin-like serine protease, is present. In this study, we cloned the CpSUB1 genomic sequence and expressed and purified the recombinant prodomain. Reverse transcriptase PCR analysis of RNA from C. parvum-infected HCT-8 cells revealed that CpSUB1 is expressed throughout infection in vitro. In immunoblots, antiserum to the recombinant CpSUB1 prodomain revealed two major bands, of approximately 64 kDa and approximately 48 kDa, for C. parvum lysates and proteins "shed" during excystation. In immunofluorescence assays, the antiserum reacted with the apical region of sporozoites and merozoites. The recombinant prodomain inhibited protease activity and processing of recombinant gp40/15 by C. parvum lysates but not by furin. Since prodomains are often selective inhibitors of their cognate enzymes, these results suggest that CpSUB1 may be a likely candidate for the protease activity in C. parvum and for processing of gp40/15. Importantly, the recombinant prodomain inhibited C. parvum infection of HCT-8 cells. These studies indicate that CpSUB1 plays a significant role in infection of host cells by the parasite and suggest that this enzyme may serve as a target for intervention.

New eukaryotic systematics: a phylogenetic perspective of developmental gene expression in the Apicomplexa

The phylum Apicomplexa consists of obligate intracellular protistan parasites, some of which are responsible for global disease causing serious morbidity and mortality in humans and animals. Understanding the mechanisms of gene expression that drive the cellular changes required to complete their life cycles will be critical in combating infection and disease. Plasmodium spp. and Toxoplasma gondii have served as good models for growth and development in the Apicomplexa. Elucidating developmental gene expression relies on comparisons with known mechanisms and their DNA, RNA and protein components. Transcriptional profiling across asexual development suggests a model where a cascade of gene expression results in a "just-in-time" production process that makes products only when needed. Some mechanisms that control transcription such as chromatin/histone modification are highly conserved in the phylum compared with the traditional model organisms, yeast, worms, flies and mammals. Studies exploiting this phenomenon show great potential for both investigating the effects of chromatin structure on developmental gene expression, and helping to identify genes that are expressed in a stage-specific manner. Transcription factors and their cognate cis-acting binding sites have been difficult to identify. This may be because the DNA binding motifs that have evolved to act as transcription factors in the Apicomplexa, e.g. the AP2 family, may be more like plants than the traditional model organisms. A new eukaryotic phylogenetic model comprised of six super-groups divides the traditional model organisms, plants and the Apicomplexa into separate super-groups. This phylogenetic model helps explain why basic functions such as transcriptional regulation appear be a composite of mechanisms in the Apicomplexa compared with what is known from other eukaryotes.

Sequencing and analysis of chromosome 1 of Eimeria tenella reveals a unique segmental organization

Eimeria tenella is an intracellular protozoan parasite that infects the intestinal tracts of domestic fowl and causes coccidiosis, a serious and sometimes lethal enteritis. Eimeria falls in the same phylum (Apicomplexa) as several human and animal parasites such as Cryptosporidium, Toxoplasma, and the malaria parasite, Plasmodium. Here we report the sequencing and analysis of the first chromosome of E. tenella, a chromosome believed to carry loci associated with drug resistance and known to differ between virulent and attenuated strains of the parasite. The chromosome--which appears to be representative of the genome--is gene-dense and rich in simple-sequence repeats, many of which appear to give rise to repetitive amino acid tracts in the predicted proteins. Most striking is the segmentation of the chromosome into repeat-rich regions peppered with transposon-like elements and telomere-like repeats, alternating with repeat-free regions. Predicted genes differ in character between the two types of segment, and the repeat-rich regions appear to be associated with strain-to-strain variation.

An efficient method for multi-locus molecular haplotyping

Many methods exist for genotyping--revealing which alleles an individual carries at different genetic loci. A harder problem is haplotyping--determining which alleles lie on each of the two homologous chromosomes in a diploid individual. Conventional approaches to haplotyping require the use of several generations to reconstruct haplotypes within a pedigree, or use statistical methods to estimate the prevalence of different haplotypes in a population. Several molecular haplotyping methods have been proposed, but have been limited to small numbers of loci, usually over short distances. Here we demonstrate a method which allows rapid molecular haplotyping of many loci over long distances. The method requires no more genotypings than pedigree methods, but requires no family material. It relies on a procedure to identify and genotype single DNA molecules, and reconstruction of long haplotypes by a 'tiling' approach. We demonstrate this by resolving haplotypes in two regions of the human genome, harbouring 20 and 105 single-nucleotide polymorphisms, respectively. The method can be extended to reconstruct haplotypes of arbitrary complexity and length, and can make use of a variety of genotyping platforms. We also argue that this method is applicable in situations which are intractable to conventional approaches.

Use of HAPPY mapping for the higher order assembly of the Tetrahymena genome

Tetrahymena thermophila is the best studied of the ciliates, a diversified and successful lineage of eukaryotic protists. Mirroring the way in which many metazoans partition their germ line and soma into distinct cell types, ciliates separate germ line and soma into two distinct nuclei in a single cell. The diploid, transcriptionally silent micronucleus undergoes meiosis and fertilization during sexual reproduction and determines the genotype of the progeny; in contrast, the expressed macronucleus contains many copies of hundreds of small chromosomes, determines the cell's phenotype, and is inherited only through vegetative reproduction. Here we demonstrate the power of HAPPY physical mapping to aid the complete assembly of T. thermophila macronuclear chromosomes from shotgun sequence scaffolds. The finished genome, one of only two ciliate genomes shotgun sequenced, will shed valuable additional light upon the biology of this extraordinary, diverse, and, from a genomics standpoint, as yet largely unexplored evolutionary branch of eukaryotes.

Differential evolution of repetitive sequences in Cryptosporidium parvum and Cryptosporidium hominis

Cryptosporidium parvum and Cryptosporidium hominis are two morphologically identical species of Apicomplexan protozoa infecting humans. Although the genomes of these species are 97% identical, their host range is strikingly different. C. parvum infects humans and animals and is primarily a zoonotic infection, whereas C. hominis is typically not detected in animals. The extent of genetic polymorphism in both species has been surveyed locally, but not on a larger geographical scale. Herein, a collection of unrelated C. parvum and C. hominis isolates was genotyped using multiple, randomly distributed micro- and minisatellites. In average, minisatellites, consisting of tandemly repeated sequence motifs of 6-24 basepair, were more polymorphic than microsatellites. When the average number of micro- and minisatellite alleles per locus was used as a measure of heterogeneity, no difference between C. parvum and C. hominis was found. However, the frequency distribution of alleles in both species was significantly different and in 6 of the 14 loci the size of the C. parvum and C. hominis repeats did not overlap. Assuming that C. parvum and C. hominis evolved from a common ancestor, these observations suggest a differential evolution of repeat length at these loci.

CryptoDB: a Cryptosporidium bioinformatics resource update

The database, CryptoDB (), is a community bioinformatics resource for the AIDS-related apicomplexan-parasite, Cryptosporidium. CryptoDB integrates whole genome sequence and annotation with expressed sequence tag and genome survey sequence data and provides supplemental bioinformatics analyses and data-mining tools. A simple, yet comprehensive web interface is available for mining and visualizing the data. CryptoDB is allied with the databases PlasmoDB and ToxoDB via ApiDB, an NIH/NIAID-fundedBioinformatics Resource Center. Recent updates to CryptoDB include the deposition of annotated genome sequences for Cryptosporidium parvum and Cryptosporidium hominis, migration to a relational database (GUS), a new query and visualization interface and the introduction of Web services.

Evolutionary, structural and functional relationships revealed by comparative analysis of syntenic genes in Rhizobiales

BACKGROUND

Comparative genomics has provided valuable insights into the nature of gene sequence variation and chromosomal organization of closely related bacterial species. However, questions about the biological significance of gene order conservation, or synteny, remain open. Moreover, few comprehensive studies have been reported for rhizobial genomes.

RESULTS

We analyzed the genomic sequences of four fast growing Rhizobiales (Sinorhizobium meliloti, Agrobacterium tumefaciens, Mesorhizobium loti and Brucella melitensis). We made a comprehensive gene classification to define chromosomal orthologs, genes with homologs in other replicons such as plasmids, and those which were species-specific. About two thousand genes were predicted to be orthologs in each chromosome and about 80% of these were syntenic. A striking gene colinearity was found in pairs of organisms and a large fraction of the microsyntenic regions and operons were similar. Syntenic products showed higher identity levels than non-syntenic ones, suggesting a resistance to sequence variation due to functional constraints; also, an unusually high fraction of syntenic products contained membranal segments. Syntenic genes encode a high proportion of essential cell functions, presented a high level of functional relationships and a very low horizontal gene transfer rate. The sequence variability of the proteins can be considered the species signature in response to specific niche adaptation. Comparatively, an analysis with genomes of Enterobacteriales showed a different gene organization but gave similar results in the synteny conservation, essential role of syntenic genes and higher functional linkage among the genes of the microsyntenic regions.

CONCLUSION

Syntenic bacterial genes represent a commonly evolved group. They not only reveal the core chromosomal segments present in the last common ancestor and determine the metabolic characteristics shared by these microorganisms, but also show resistance to sequence variation and rearrangement, possibly due to their essential character. In Rhizobiales and Enterobacteriales, syntenic genes encode a high proportion of essential cell functions and presented a high level of functional relationships.

Gene analysis of Cryptosporidium parvum HNJ-1 strain isolated in Japan

We analyzed genetically Cryptosporidium parvum HNJ-1 strain, which is the Japanese reference strain isolated from human in Japan. DNA sequences of genes for thrombospondin-related adhesive protein of Cryptosporidium-1 and Cryptosporidium-2 (TRAP-C1, TRAP-C2), heat shock protein 70 (HSP70), oocyst wall protein (COWP), beta-tubulin, alpha-tubulin, polythreonine-region (Poly-T), elongation factor 1 alpha (EF-1 alpha), and 18S rRNA of this strain were determined. They showed high rate of homology to published sequences of genotype 2 strains, which were considered to be infective to both humans and animals. However, HNJ-1 had synonymous and non-synonymous substitutions in the nucleotide sequence of TRAP-C1 and beta-tubulin among HNJ-1 and published sequences of genotype 2 strains. These results implied that HNJ-1 strain was an unique subpopulation of genotype 2 strain of C. parvum.

Cryptosporidium parvum mitochondrial-type HSP70 targets homologous and heterologous mitochondria

A mitochondrial HSP70 gene (Cp-mtHSP70) is described for the apicomplexan Cryptosporidium parvum, an agent of diarrhea in humans and animals. Mitochondrial HSP70 is known to have been acquired from the proto-mitochondrial endosymbiont. The amino acid sequence of Cp-mtHSP70 shares common domains with mitochondrial and proteobacterial homologues, including 34 amino acids of an NH2-terminal mitochondrion-like targeting presequence. Phylogenetic reconstruction places Cp-mtHSP70 within the mitochondrial clade of HSP70 homologues. Using reverse transcription-PCR, Cp-mtHSP70 mRNA was observed in C. parvum intracellular stages cultured in HCT-8 cells. Polyclonal antibodies to Cp-mtHSP70 recognize a approximately 70-kDa protein in Western blot analysis of sporozoite extracts. Both fluorescein- and immunogold-labeled anti-Cp-mtHSP70 localize to a single mitochondrial compartment in close apposition to the nucleus. Furthermore, the NH2-terminal presequence of Cp-mtHSP70 can correctly target green fluorescent protein to the single mitochondrion of the apicomplexan Toxoplasma gondii and the mitochondrial network of the yeast Saccharomyces cerevisiae. When this presequence was truncated, the predicted amphiphilic alpha-helix was shown to be essential for import into the yeast mitochondrion. These data further support the presence of a secondarily reduced relict mitochondrion in C. parvum.

Phylogenomic evidence supports past endosymbiosis, intracellular and horizontal gene transfer in Cryptosporidium parvum

BACKGROUND

The apicomplexan parasite Cryptosporidium parvum is an emerging pathogen capable of causing illness in humans and other animals and death in immunocompromised individuals. No effective treatment is available and the genome sequence has recently been completed. This parasite differs from other apicomplexans in its lack of a plastid organelle, the apicoplast. Gene transfer, either intracellular from an endosymbiont/donor organelle or horizontal from another organism, can provide evidence of a previous endosymbiotic relationship and/or alter the genetic repertoire of the host organism. Given the importance of gene transfers in eukaryotic evolution and the potential implications for chemotherapy, it is important to identify the complement of transferred genes in Cryptosporidium.

RESULTS

We have identified 31 genes of likely plastid/endosymbiont (n = 7) or prokaryotic (n = 24) origin using a phylogenomic approach. The findings support the hypothesis that Cryptosporidium evolved from a plastid-containing lineage and subsequently lost its apicoplast during evolution. Expression analyses of candidate genes of algal and eubacterial origin show that these genes are expressed and developmentally regulated during the life cycle of C. parvum.

CONCLUSIONS

Cryptosporidium is the recipient of a large number of transferred genes, many of which are not shared by other apicomplexan parasites. Genes transferred from distant phylogenetic sources, such as eubacteria, may be potential targets for therapeutic drugs owing to their phylogenetic distance or the lack of homologs in the host. The successful integration and expression of the transferred genes in this genome has changed the genetic and metabolic repertoire of the parasite.

Complete genome sequence of the apicomplexan, Cryptosporidium parvum

The apicomplexan Cryptosporidium parvum is an intestinal parasite that affects healthy humans and animals, and causes an unrelenting infection in immunocompromised individuals such as AIDS patients. We report the complete genome sequence of C. parvum, type II isolate. Genome analysis identifies extremely streamlined metabolic pathways and a reliance on the host for nutrients. In contrast to Plasmodium and Toxoplasma, the parasite lacks an apicoplast and its genome, and possesses a degenerate mitochondrion that has lost its genome. Several novel classes of cell-surface and secreted proteins with a potential role in host interactions and pathogenesis were also detected. Elucidation of the core metabolism, including enzymes with high similarities to bacterial and plant counterparts, opens new avenues for drug development.

The iron-sulfur cluster assembly genes iscS and iscU of Entamoeba histolytica were acquired by horizontal gene transfer

Background

Iron-sulfur (FeS) proteins are present in all living organisms and play important roles in electron transport and metalloenzyme catalysis. The maturation of FeS proteins in eukaryotes is an essential function of mitochondria, but little is known about this process in amitochondriate eukaryotes. Here we report on the identification and analysis of two genes encoding critical FeS cluster (Isc) biosynthetic proteins from the amitochondriate human pathogen Entamoeba histolytica.

Results

E. histolytica IscU and IscS were found to contain all features considered essential for their biological activity, including amino acid residues involved in substrate and/or co-factor binding. The IscU protein differs significantly from other eukaryotic homologs and resembles the long type isoforms encountered in some bacteria. Phylogenetic analyses of E. histolytica IscS and IscU showed a close relationship with homologs from Helicobacter pylori and Campylobacter jejuni, to the exclusion of mitochondrial isoforms.

Conclusions

The bacterial-type FeS cluster assembly genes of E. histolytica suggest their lateral acquisition from epsilon proteobacteria. This is a clear example of horizontal gene transfer (HGT) from eubacteria to unicellular eukaryotic organisms, a phenomenon known to contribute significantly to the evolution of eukaryotic genomes.

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.