Genetic mapping in the human malaria parasite Plasmodium falciparum

Genetic mapping of determinants in drug resistance, virulence, disease susceptibility, and interaction of host-rodent malaria parasites

Genetic mapping has been widely employed to search for genes linked to phenotypes/traits of interest. Because of the ease of maintaining rodent malaria parasites in laboratory mice, many genetic crosses of rodent malaria parasites have been performed to map the parasite genes contributing to malaria parasite development, drug resistance, host immune response, and disease pathogenesis. Drs. Richard Carter, David Walliker, and colleagues at the University of Edinburgh, UK, were the pioneers in developing the systems for genetic mapping of malaria parasite traits, including characterization of genetic markers to follow the inheritance and recombination of parasite chromosomes and performing the first genetic cross using rodent malaria parasites. Additionally, many genetic crosses of inbred mice have been performed to link mouse chromosomal loci to the susceptibility to malaria parasite infections. In this chapter, we review and discuss past and recent advances in genetic marker development, performing genetic crosses, and genetic mapping of both parasite and host genes. Genetic mappings using models of rodent malaria parasites and inbred mice have contributed greatly to our understanding of malaria, including parasite development within their hosts, mechanism of drug resistance, and host-parasite interaction.

Sources of transcription variation in Plasmodium falciparum

Variation in transcript abundance can contribute to both short-term environmental response and long-term evolutionary adaptation. Most studies are designed to assess differences in mean transcription levels and do not consider other potentially important and confounding sources of transcriptional variation. Detailed quantification of variation sources will improve our ability to detect and identify the mechanisms that contribute to genome-wide transcription changes that underpin adaptive responses. To quantify innate levels of expression variation, we measured mRNA levels for more than 5000 genes in the malaria parasite, Plasmodium falciparum, among clones derived from two parasite strains across biologically and experimentally replicated batches. Using a mixed effects model, we partitioned the total variation among four sources - between strain, within strain, environmental batch effects, and stochastic noise. We found 646 genes with significant variation attributable to at least one of these sources. These genes were categorized by their predominant variation source and further examined using gene ontology enrichment analysis to associate function with each source of variation. Genes with environmental batch effect and within strain transcript variation may contribute to phenotypic plasticity, while genes with between strain variation may contribute to adaptive responses and processes that lead to parasite strain-specific survival under varied conditions.

A set of microsatellite markers to differentiate Plasmodium falciparum progeny of four genetic crosses

Background

Four Plasmodium falciparum genetic crosses (HB3×3D7, HB3×Dd2, 7G8×GB4, and 803×GB4) have produced sets of recombinant progeny that are widely used for malaria research, including investigations of anti-malarial drug resistance. It is critical to maintain the progeny free from cross-contamination. Microsatellite polymorphisms can be used to validate parasite identity.

Results

A set of 12 markers was developed that differentiates the parents of the four P. falciparum crosses. This typing set identified distinguishing patterns of inheritance (fingerprints) in segregant collections of 15 progeny clones from HB3×3D7, 32 from HB3×Dd2, 33 from 7G8×GB4, and 81 from 803×GB4. Stronger amplification was observed with shorter relative to longer alleles of individual microsatellites. In experiments with mixed parental DNAs, electropherograms showed that signals of cross-contamination can be missed when minor peaks less than 1/4 or 1/3 the height of the major peak are disregarded by threshold settings commonly used for population studies.

Conclusions

Microsatellite typing is an effective method to check the identity of P. falciparum lines and detect parasite cross-contamination in cultures; however, care must be taken not to ignore minor peaks that can be overlooked. The 12 microsatellite markers presented here provide a rapid and efficient means to distinguish the segregants of laboratory crosses. Fingerprint patterns from these markers are useful to maintain the integrity of diverse parasite lines in and between research laboratories.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2210-z) contains supplementary material, which is available to authorized users.

Multifunctional Involvement of a C2H2 Zinc Finger Protein (PbZfp) in Malaria Transmission, Histone Modification, and Susceptibility to DNA Damage Response

In sexually reproducing organisms, meiosis is an essential step responsible for generation of haploid gametes from diploid somatic cells. The quest for understanding regulatory mechanisms of meiotic recombination in Plasmodium led to identification of a gene encoding a protein that contains 11 copies of C₂H₂ zinc fingers (ZnF). Reverse genetic approaches were used to create Plasmodium berghei parasites either lacking expression of full-length Plasmodium berghei zinc finger protein (PbZfp) (knockout [KO]) or expressing PbZfp lacking C-terminal zinc finger region (truncated [Trunc]). Mice infected with KO parasites survived two times longer (P < 0.0001) than mice infected with wild-type (WT) parasites. In mosquito transmission experiments, the infectivity of KO and Trunc parasites was severely compromised (>95% oocyst reduction). KO parasites revealed a total lack of trimethylation of histone 3 at several lysine residues (K4, K27, and K36) without any effect on acetylation patterns (H3K9, H3K14, and H4K16). Reduced DNA damage and reduced expression of topoisomerase-like Spo11 in the KO parasites with normal Rad51 expression further suggest a functional role for PbZfp during genetic recombination that involves DNA double-strand break (DSB) formation followed by DNA repair. These finding raise the possibility of some convergent similarities of PbZfp functions to functions of mammalian PRDM9, also a C₂H₂ ZnF protein with histone 3 lysine 4 (H3K4) methyltransferase activity. These functions include the major role played by the latter in binding recombination hotspots in the genome during meiosis and trimethylation of the associated histones and subsequent chromatin recruitment of topoisomerase-like Spo11 to catalyze DNA DSB formation and DMC1/Rad51-mediated DNA repair and homologous recombination.

Malaria parasites are haploid throughout their life cycle except for a brief time period when zygotes are produced as a result of fertilization between male and female gametes during transmission through the mosquito vector. The reciprocal recombination events that follow zygote formation ensure orderly segregation of homologous chromosomes during meiosis, creating genetic diversity among offspring. Studies presented in the current manuscript identify a novel C2H2 ZnF-containing protein exhibiting multifunctional roles in parasite virulence, mosquito transmission, and homologous recombination during meiosis. Understanding the transmission biology of malaria will result in the identification of novel targets for transmission-blocking intervention approaches.

Distinct genomic architecture of Plasmodium falciparum populations from South Asia

Previous whole genome comparisons of Plasmodium falciparum populations have not included collections from the Indian subcontinent, even though two million Indians contract malaria and about 50,000 die from the disease every year. Stratification of global parasites has revealed spatial relatedness of parasite genotypes on different continents. Here, genomic analysis was further improved to obtain country-level resolution by removing var genes and intergenic regions from distance calculations. P. falciparum genomes from India were found to be most closely related to each other. Their nearest neighbors were from Bangladesh and Myanmar, followed by Thailand. Samples from the rest of Southeast Asia, Africa and South America were increasingly more distant, demonstrating a high-resolution genomic-geographic continuum. Such genome stratification approaches will help monitor variations of malaria parasites within South Asia and future changes in parasite populations that may arise from in-country and cross-border migrations.

A P450 gene associated with robust resistance to DDT in ciliated protozoan, Tetrahymena thermophila by efficient degradation

Analysis of metabolic mechanisms of dichlorodiphenyltrichloroethane (DDT) accumulation and degradation in microorganisms, which could be used to reduce its hazard to higher organisms at the higher in the food chain, have not been investigated. Robust resistance to DDT (grows well in 256 mg/L DDT) and a surprising ability to degrade DDT (more than 70% DDT within 4h) were found in the ciliated protozoan Tetrahymena thermophila. A P450 gene (CYP5013C2) was found to respond specifically to DDT treatment. In the presence of 256 mg/L DDT, cells with overexpressing CYP5013C2 (p450-OE) grew faster and degraded DDT more efficiently than wild-type (WT) cells, while cells with CYP5013C2 partially knocked down (p450-KD) grew slower and exhibited reduced ability to degrade DDT compared to WT cells. Both dichlorodiphenyldichloroethylene (DDE) and dichlorodiphenyldichloroethane (DDD) were detected in cells after exposure to DDT, and the concentration of DDD in the p450-OE strain gradually decreased from 0.5 to 4h. Thus, we argue that this P450 gene (CYP5013C2), by efficiently degrading DDT to DDD, is associated with robust resistance to DDT in Tetrahymena, and that a strain overexpressing this gene has the potential to serve as bioreactor that degrades environmental DDT.

High prevalence of pfcrt-CVIET haplotype in isolates from asymptomatic and symptomatic patients in south-central Oromia, Ethiopia

Background

As a result of extensive chloroquine resistance (CQR) in Plasmodium falciparum in late 1990s, Ethiopia replaced CQ with sulphadoxine-pyrimethamine (SP) as first-line drug, which in turn was replaced by artemisinin combination therapy in 2004. Plasmodium falciparum resistance to CQ is determined by the mutation at K76T of the P. falciparum chloroquine resistance transporter (pfcrt) gene. Understanding diversity in the P. falciparum genome is crucial since it has the potential to influence important phenotypes of the parasite such as drug resistance. Limited data is available regarding the type of pfcrt mutant allelic type, the effect of CQ withdrawal and diversity of the parasite population in south-central Oromia, Ethiopia.

Methods

Finger-pricked blood spotted on Whatman 3MM filter papers were collected from falciparum malaria patients. Parasite DNA was extracted from individual blood spots on the filter papers. The presence of K76T mutations was determined using nested PCR for all isolates. Complete sequencing of mutations in pfcrt 72-76 was done for a set of randomly selected resistant isolates. Four microsatellite (MS) markers were analysed to determine the heterozygosity.

Results

Although CQ was withdrawn for more than a decade, 100% of the parasites still carried the pfcrt K76T mutation. All isolates were mutant at the K76T polymorphism. Based on combinations of MS markers, seven different Ethiopian CQR variants (E1-E7) were identified. Heterozygosity (H_e) for MS flanking the pfcrt chloroquine resistance allele ranged from 0.00 (mscrt -29, -29.268 kb) to 0.21 (mscrt -2, -2.814 kb). H_e ranged from 0.00 (msint 3, 0 kb) to 0.19 (msint 2, 0 kb) for MS within the pfcrt gene. Both intronic and MS flanking the pfcrt gene showed low levels of diversity.

Conclusion

pfcrt CQR allele seems to be fixed in the study area. Of the different haplotypes associated with CQR, only the CVIET genotype was identified. No reversal to the wild-type has occurred in Ethiopia unlike in many Africa countries where CQR parasites declined after cessation of CQ use. Decreased diversity in CQR isolates surrounding pfcrt suggests CQ selection and homogenization among CQR parasite population. While mutation in msint 3 and mscrt -29 of the mutant pfcrt allele is being fixed, it seems that mutations in msint 2 and mscrt -2 are still evolving and may indicate the start of re-diversification of the population from a fixed 76 T population.

The landscape of inherited and de novo copy number variants in a Plasmodium falciparum genetic cross

BACKGROUND

Copy number is a major source of genome variation with important evolutionary implications. Consequently, it is essential to determine copy number variant (CNV) behavior, distributions and frequencies across genomes to understand their origins in both evolutionary and generational time frames. We use comparative genomic hybridization (CGH) microarray and the resolution provided by a segregating population of cloned progeny lines of the malaria parasite, Plasmodium falciparum, to identify and analyze the inheritance of 170 genome-wide CNVs.

RESULTS

We describe CNVs in progeny clones derived from both Mendelian (i.e. inherited) and non-Mendelian mechanisms. Forty-five CNVs were present in the parent lines and segregated in the progeny population. Furthermore, extensive variation that did not conform to strict Mendelian inheritance patterns was observed. 124 CNVs were called in one or more progeny but in neither parent: we observed CNVs in more than one progeny clone that were not identified in either parent, located more frequently in the telomeric-subtelomeric regions of chromosomes and singleton de novo CNVs distributed evenly throughout the genome. Linkage analysis of CNVs revealed dynamic copy number fluctuations and suggested mechanisms that could have generated them. Five of 12 previously identified expression quantitative trait loci (eQTL) hotspots coincide with CNVs, demonstrating the potential for broad influence of CNV on the transcriptional program and phenotypic variation.

CONCLUSIONS

CNVs are a significant source of segregating and de novo genome variation involving hundreds of genes. Examination of progeny genome segments provides a framework to assess the extent and possible origins of CNVs. This segregating genetic system reveals the breadth, distribution and dynamics of CNVs in a surprisingly plastic parasite genome, providing a new perspective on the sources of diversity in parasite populations.

Watch out for malaria: still a leading cause of child death worldwide

Background

Due to the efforts in malaria control promoted by the World Health Organization (WHO), the reported malaria burden is being reduced throughout the world. Nevertheless, malaria remains a leading cause of child death worldwide.

Aims

purpose of the paper is to summarize the main historical steps in fighting malaria, from the first descriptions to the last ones.

Results

a case of probable autochthonous malaria has been recently described in Italy, raising concern over the possibility of resurgence of malaria in countries previously interested by this disease. Moreover, both the constant threat of the parasite and vector mosquito developing resistance to medicines and insecticides, and the on-going climate change make the challenge of eradicating malaria really difficult. Therefore, malaria is still an actual disease, requiring adequate programs of surveillance, stronger health systems in poor countries, and efforts in order to develop new and effective tools in malaria control. WHO has definitely demonstrated the effects of "social determinants" on health. So, eradication strategies cannot be based only on a scientific background, because culture, politics, power, resources and wars have a profound impact on health and disease. These elements should be introduced in all the programs of malaria control.

Conclusions

malaria is still an actual disease with great public health implications, and the approaches for control and prevention should have the appropriate social and political context in addition to the science involved in order to save lives of children at risk.

Discovery, mechanisms of action and combination therapy of artemisinin

Despite great international efforts, malaria still inflicts an enormous toll on human lives, especially in Africa. Throughout history, antimalarial medicines have been one of the most powerful tools in malaria control. However, the acquisition and spread of parasite strains that are resistant to multiple antimalarial drugs have become one of the greatest challenges to malaria treatment, and are associated with the increase in morbidity and mortality in many malaria-endemic countries. To deal with this grave situation, artemisinin-based combinatory therapies (ACTs) have been introduced and widely deployed in malarious regions. Artemisinin is a new class of antimalarial compounds discovered by Chinese scientists from the sweet wormwood Artemisia annua. The potential development of resistance to artemisinins by Plasmodium falciparum threatens the usable lifespan of ACTs, and therefore is a subject of close surveillance and extensive research. Studies at the Thai-Cambodian border, a historical epicenter of multidrug resistance, have detected reduced susceptibility to artemisinins as manifested by prolonged parasite-clearance times, raising considerable concerns on resistance development. Despite this significance, there is still controversy on the mode of action of artemisinins. Although a number of potential cellular targets of artemisinins have been proposed, they remain to be verified experimentally. Here, we review the history of artemisinin discovery, discuss the mode of action and potential drug targets, and present strategies to elucidate resistance mechanisms.

Graph theoretical approach to study eQTL: a case study of Plasmodium falciparum

Motivation: Analysis of expression quantitative trait loci (eQTL) significantly contributes to the determination of gene regulation programs. However, the discovery and analysis of associations of gene expression levels and their underlying sequence polymorphisms continue to pose many challenges. Methods are limited in their ability to illuminate the full structure of the eQTL data. Most rely on an exhaustive, genome scale search that considers all possible locus–gene pairs and tests the linkage between each locus and gene.

Result: To analyze eQTLs in a more comprehensive and efficient way, we developed the Graph based eQTL Decomposition method (GeD) that allows us to model genotype and expression data using an eQTL association graph. Through graph-based heuristics, GeD identifies dense subgraphs in the eQTL association graph. By identifying eQTL association cliques that expose the hidden structure of genotype and expression data, GeD effectively filters out most locus–gene pairs that are unlikely to have significant linkage. We apply GeD on eQTL data from Plasmodium falciparum, the human malaria parasite, and show that GeD reveals the structure of the relationship between all loci and all genes on a whole genome level. Furthermore, GeD allows us to uncover additional eQTLs with lower FDR, providing an important complement to traditional eQTL analysis methods.

Contact:przytyck@ncbi.nlm.nih.gov

Large-scale genotyping and genetic mapping in Plasmodium parasites

The completion of many malaria parasite genomes provides great opportunities for genomewide characterization of gene expression and high-throughput genotyping. Substantial progress in malaria genomics and genotyping has been made recently, particularly the development of various microarray platforms for large-scale characterization of the Plasmodium falciparum genome. Microarray has been used for gene expression analysis, detection of single nucleotide polymorphism (SNP) and copy number variation (CNV), characterization of chromatin modifications, and other applications. Here we discuss some recent advances in genetic mapping and genomic studies of malaria parasites, focusing on the use of high-throughput arrays for the detection of SNP and CNV in the P. falciparum genome. Strategies for genetic mapping of malaria traits are also discussed.

Drug resistance and genetic mapping in Plasmodium falciparum

Drug resistance in malaria parasites is a serious public health burden, and resistance to most of the antimalarial drugs currently in use has been reported. A better understanding of the molecular mechanisms of drug resistance is urgently needed to slow or circumvent the spread of resistance, to allow local treatments to be deployed more effectively to prolong the life span of the current drugs, and to develop new drugs. Although mutations in genes determining resistance to drugs such as chloroquine and the antifolates have been identified, we still do not have a full understanding of the resistance mechanisms, and genes that contribute to resistance to many other drugs remain to be discovered. Genetic mapping is a powerful tool for the identification of mutations conferring drug resistance in malaria parasites because most drug-resistant phenotypes were selected within the past 60 years. High-throughput methods for genotyping large numbers of single nucleotide polymorphisms (SNPs) and microsatellites (MSs) are now available or are being developed, and genome-wide association studies for malaria traits will soon become a reality. Here we discuss strategies and issues related to mapping genes contributing to drug resistance in the human malaria parasite Plasmodium falciparum.

Accumulation of CVIET Pfcrt allele of Plasmodium falciparum in placenta of pregnant women living in an urban area of Dakar, Senegal

OBJECTIVES

Chemoprophylaxis is recommended during pregnancy to reduce the risk of placental infection. However, in areas with increasing drug resistance, it can trigger selection of resistant parasites in the placenta and increase the frequency of placental malaria. The objective of this study was to analyse the selection of drug-resistant parasites in the placenta in an area where chloroquine was still recommended as prophylaxis.

PATIENTS AND METHODS

We analysed the polymorphism of parasites from matched placental and venous blood samples at the time of delivery from women in Dakar. Polymorphism of the isolates was studied using nested PCR typing of MSA1 and MSA2 genes, and full sequence of PfCRT exon 2.

RESULTS

Of 692 women recruited at delivery, 72 had placental malaria. Two Pfcrt exon 2 genotypes were found, and 86% of the placentas had monoallelelic CVIET infection compared with 39% that had peripheral blood infection. Mixed parasite populations of CVIET/CVMNK occurred in 53% of the peripheral blood samples but only in 7% of the infected placentas. This selection of CVIET in placenta was not related to a decreased polymorphism of the parasites, as a large diversity of MSA1 and MSA2 was found in both placenta and venous blood. This diversity confirms that a multiplicity of circulation isolates can occur at low parasite transmission. msp1 and msp2 genotyping revealed mostly distinct populations of parasites in venous and placental blood.

CONCLUSIONS

These data suggest that, even in low transmission areas, diverse parasite populations can accumulate in the placenta during pregnancy despite strong selection at the PfCRT locus due to chemoprophylaxis with chloroquine.

A network to monitor antimalarial drug resistance: a plan for moving forward

The spread of resistance to antimalarial drugs has required changes in the recommended first-line treatment for falciparum malaria in almost all regions. Most drugs recommended currently are combinations of a long-acting antimalarial and an artemisinin derivative. This article presents the rationale for establishing a web-based, open-access database of antimalarial drug resistance and efficacy: the World Antimalarial Resistance Network (WARN). The goal of this network is to assemble the tools and information that will enable the malaria community to collate, analyze and share contemporary information on antimalarial-drug efficacy in all endemic regions so that decisions on antimalarial-drug use are based on solid evidence.

Using the genome to dissect the molecular basis of drug resistance

The need to understand the genetic basis of drug resistance in human pathogens has never been greater. The global incidence of drug-resistant organisms, such as those that cause malaria, continues to rise, while the repertoire of effective, inexpensive drugs is declining. Genomic technologies, such as DNA microarrays and full-genome sequencing offer new hope in advancing our understanding of the underlying genetic processes that facilitate a resistance phenotype. Importantly, evidence that drug resistance in many organisms can be a multigene, complex phenomenon implies that unbiased, genome-wide scans of diversity will be required to fully understand the molecular mechanisms of both established and novel resistance traits. While the potential application of full-genome approaches for deciphering mechanisms of drug resistance has yet to be fully realized, this review evaluates drug resistance in human malaria parasites and discusses the exciting role genome-based systems could play in monitoring drug resistance, as well as guiding the implementation of efficient therapeutic strategies for malaria. The approaches reviewed within this article will be applicable to all known or emerging microbial pathogens.

Contrasting genetic structure in Plasmodium vivax populations from Asia and South America

Populations of Plasmodium falciparum show striking differences in linkage disequilibrium, population differentiation and diversity, but only fragmentary data exists on the genetic structure of Plasmodium vivax. We genotyped nine tandem repeat loci bearing 2-8 bp motifs from 345 P. vivax infections collected from three Asian countries and from five locations in Colombia. We observed 9-37 alleles per locus and high diversity (He=0.72-0.79, mean=0.75) in all countries. Numbers of multiple clone infections varied considerably: these were rare in Colombia and India, but > 60% of isolates carried multiple alleles in at least one locus in Thailand and Laos. However, only one or two of the nine loci show >1 allele in many samples, suggesting that mutation within infections may result in overestimation of true multiple carriage rates. Identical nine-locus genotypes were frequently found in Colombian populations, contributing to strong linkage disequilibrium. These identical genotypes were strongly clustered in time, consistent with epidemic transmission of clones and subsequent breakdown of allelic associations, suggesting high rates of inbreeding and low effective recombination rates in this country. In contrast, identical genotypes were rare and loci were randomly associated in all three Asian populations, consistent with higher rates of outcrossing and recombination. We observed low but significant differentiation between different Asian countries (standardized FST = 0.13-0.45). In comparison, we see greater differentiation between collection locations within Colombia (standardized FST = 0.4-0.7), and strong differentiation between continents (standardized FST = 0.48-0.79). The observed heterogeneity in multiple clone carriage rates, linkage disequilibrium and population differentiation are similar in some, but not all, respects to those observed in P. falciparum, and have important implications for the design of association mapping studies, and interpretation of P. vivax epidemiology.

Development and characterization of microsatellite markers for the canine hookworm, Ancylostoma caninum

Microsatellites are repetitive genomic elements that show high levels of variation and therefore provide excellent tools to study the genetics of eukaryotic organisms. Hookworms are extremely common and important nematode parasites of humans and animals, causing potentially serious disease morbidity. Control of hookworms in dogs is achieved by frequent treatment with anthelmintics, and in humans, anthelmintics are frequently administered in a mass-treatment community-wide approach. Understanding the population genetics of hookworms has important implications for studies on the development and spread of drug resistance. We investigated the genome of Ancylostoma caninum for microsatellites by developing and then screening an enriched genomic library as well as by data mining published sequences of a whole genome shotgun library. Investigations revealed a high abundance of trinucleotide repeats. Dinucleotide repeats were characterized by a high number of AT, GA, and GT repeats. After testing and optimization of 68 markers, a panel of 34 polymorphic microsatellite markers were selected. Microsatellite analysis of hookworm isolates revealed a high degree of polymorphism, which was not influenced by the length of the repeats. This panel of microsatellite markers makes it possible to pursue investigations on the population genetics of A. caninum. Furthermore, a number of the markers demonstrated suitability for analysis of the human hookworm species Necator americanus and A. duodenale.

Comparative genomic and phylogenetic analyses of calcium ATPases and calcium-regulated proteins in the apicomplexa

The phylum Apicomplexa comprises a large group of early branching eukaryotes that includes a number of human and animal parasites. Calcium controls a number of vital processes in apicomplexans including protein secretion, motility, and differentiation. Despite the importance of calcium as a second messenger, very little is known about the systems that control homeostasis or that regulate calcium signaling in parasites. The recent completion of many apicomplexan genomes provides new opportunity to define calcium response pathways in this group of parasites in comparison to model organisms. Whole-genome comparison between the apicomplexans Plasmodium spp., Cryptosporidium spp., and Toxoplasma gondii revealed the presence of several P-Type Ca2+ transporting ATPases including a single endoplasmic reticulum (ER)-type sarcoplasmic-endoplasmic reticulum Ca2+ ATPase, several Golgi-like Ca2+ ATPases, and a single Ca2+/H+ exchanger. Only T. gondii showed evidence of plasma membrane-type Ca2+ ATPases or voltage-gated calcium channels. Despite pharmacological evidence for IP3 and ryanodine-mediated calcium release, animal-type calcium channels were not readily identified in parasites, indicating they are more similar to plants. Downstream of calcium release, a variety of EF-hand-containing proteins regulate calcium responses. Our analyses detected a single conserved calmodulin (CaM) homologue, 3 distinct centrin (CETN)-caltractin-like proteins, one of which is shared with ciliates, and a variety of deep-branching, CaM-CETN-like proteins. Apicomplexans were also found to contain a wide array of calcium-dependent protein kinases (CDPKs), which are commonly found in plants. Toxoplasma gondii contains more than 20 CDPK or CDPK-related kinases, which likely regulate a variety of responses including secretion, motility, and differentiation. Genomic and phylogenetic comparisons revealed that apicomplexans contain a variety of unusual calcium response pathways that are distinct from those seen in vertebrates. Notably, plant-like pathways for calcium release channels and calcium-dependent kinases are found in apicomplexans. The experimental flexibility of T. gondii should allow direct experimental manipulation of these pathways to validate their biological roles. The central importance of calcium in signaling and development, and the novel characteristics of many of these systems, indicates that parasite calcium pathways may be exploited as new therapeutic targets for intervention.

SeqDoC: rapid SNP and mutation detection by direct comparison of DNA sequence chromatograms

Background

This paper describes SeqDoC, a simple, web-based tool to carry out direct comparison of ABI sequence chromatograms. This allows the rapid identification of single nucleotide polymorphisms (SNPs) and point mutations without the need to install or learn more complicated analysis software.

Results

SeqDoC produces a subtracted trace showing differences between a reference and test chromatogram, and is optimised to emphasise those characteristic of single base changes. It automatically aligns sequences, and produces straightforward graphical output. The use of direct comparison of the sequence chromatograms means that artefacts introduced by automatic base-calling software are avoided. Homozygous and heterozygous substitutions and insertion/deletion events are all readily identified. SeqDoC successfully highlights nucleotide changes missed by the Staden package 'tracediff' program.

Conclusion

SeqDoC is ideal for small-scale SNP identification, for identification of changes in random mutagenesis screens, and for verification of PCR amplification fidelity. Differences are highlighted, not interpreted, allowing the investigator to make the ultimate decision on the nature of the change.

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

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