Artemisinin resistance in Plasmodium falciparum is associated with an altered temporal pattern of transcription
- Sachel Mok†1,
- Mallika Imwong†2, 3,
- Margaret J Mackinnon4,
- Joan Sim1,
- Ramya Ramadoss1,
- Poravuth Yi9,
- Mayfong Mayxay5, 6,
- Kesinee Chotivanich2,
- Kek-Yee Liong1,
- Bruce Russell7,
- Duong Socheat9,
- Paul N Newton5, 8,
- Nicholas PJ Day3, 8,
- Nicholas J White3, 8,
- Peter R Preiser1,
- François Nosten8, 10,
- Arjen M Dondorp3, 8Email author and
- Zbynek Bozdech1Email author
© Mok et al; licensee BioMed Central Ltd. 2011
Received: 25 March 2011
Accepted: 3 August 2011
Published: 3 August 2011
Artemisinin resistance in Plasmodium falciparum malaria has emerged in Western Cambodia. This is a major threat to global plans to control and eliminate malaria as the artemisinins are a key component of antimalarial treatment throughout the world. To identify key features associated with the delayed parasite clearance phenotype, we employed DNA microarrays to profile the physiological gene expression pattern of the resistant isolates.
In the ring and trophozoite stages, we observed reduced expression of many basic metabolic and cellular pathways which suggests a slower growth and maturation of these parasites during the first half of the asexual intraerythrocytic developmental cycle (IDC). In the schizont stage, there is an increased expression of essentially all functionalities associated with protein metabolism which indicates the prolonged and thus increased capacity of protein synthesis during the second half of the resistant parasite IDC. This modulation of the P. falciparum intraerythrocytic transcriptome may result from differential expression of regulatory proteins such as transcription factors or chromatin remodeling associated proteins. In addition, there is a unique and uniform copy number variation pattern in the Cambodian parasites which may represent an underlying genetic background that contributes to the resistance phenotype.
The decreased metabolic activities in the ring stages are consistent with previous suggestions of higher resilience of the early developmental stages to artemisinin. Moreover, the increased capacity of protein synthesis and protein turnover in the schizont stage may contribute to artemisinin resistance by counteracting the protein damage caused by the oxidative stress and/or protein alkylation effect of this drug. This study reports the first global transcriptional survey of artemisinin resistant parasites and provides insight to the complexities of the molecular basis of pathogens with drug resistance phenotypes in vivo.
KeywordsPlasmodium falciparum, in vivo artemisinin-resistance field isolates comparative genomics comparative transcriptomics
Artemisinin combination therapy (ACT) is recommended by the World Health Organization as the first-line treatment for falciparum malaria in all endemic regions [1, 2]. The excellent effectiveness and tolerability of ACTs brought new enthusiasm into world-wide efforts to eliminate human malaria which until today accounts for 243 million cases of infection and 863,000 deaths per annum . The core components of ACTs - artemisinin and its derivatives, provide an important alternative to quinoline and antifolate-based compounds. Resistance to these older compounds that emerged on the Thai-Cambodian border and subsequently spread across the world has severely compromised their use and contributed to a dramatic rise in malaria morbidity prior to introduction of the ACTs in the late 1990's [4–7]. Learning from past mistakes, much effort is being invested in proper management of ACTs in order to sustain their efficacy and prevent the spread of resistance .
In spite of these efforts, there have been sporadic reports of artemisinin resistance in-vivo and in-vitro for many years (from Yunnan Province, Southwest China , Vietnam  and French Guiana ). Although the biological and clinical significance of these reports were uncertain , these early warning signs suggested a possibility of emergence of malaria parasites resistant to artemisinin [12–14]. Recently, unequivocal evidence of reduced artemisinin susceptibility from Western Cambodia has been reported . Curiously, this was also the epicenter of chloroquine and sulfadoxine-pyrimethamine resistance. Dondorp et al. (2009) documented markedly prolonged parasite clearance times (median PCT 84 hours (interquartile range 60 to 96 hours) in Pailin, Western Cambodia. This compares with a median PCT of 48 hours (36 to 66 hours) on the Western border of Thailand . Since this study, reports of delayed parasite clearance have emerged in other parts of the region, including the Thai-Myanmar border . Although it has yet to be established whether artemisinin resistance has spread westward, the possibility of the spread of resistant parasites through Asia to Africa would be disastrous.
The mechanism of artemisinin resistance is unknown. The resistant phenotype detected in Western Cambodia does not associate with any polymorphisms in the established drug resistance markers . In vitro susceptibility testing of parasites which are cleared abnormally slowly in vivo, showed essentially no shift in IC50 (50% inhibition concentration) values in vitro . This apparent discrepancy between the experimental and the clinical data may be explained by the reduced susceptibility of Plasmodium parasites at only the ring stage (first third) of its 48 hours intraerythrocytic developmental cycle (IDC) [15, 16]. Another phenomenon that has been suggested to explain artemisinin resistance is an increased propensity for these parasites to form "dormant" (or quiescent) rings under artemisinin exposure [17, 18]. However, this process, suggested by in vitro studies, is unlikely to explain the slow first order decline in parasitemia with time that was observed in Western Cambodian patients treated with artemisinin-based drugs . The reduced artemisinin susceptibility phenotype of the resistant malaria parasites exhibits a heritable pattern suggesting that it has a genetic basis . Identifying the genetic determinants will be crucial for understanding the molecular basis of artemisinin resistance and will also provide an important molecular tool for epidemiological surveys.
Here we carry out genome-wide gene expression analyses in order to identify key elements of a transcriptional profile underlining artemisinin resistance. We show that the P. falciparum parasites with slow clearance after artemisinin treatment exhibit reduced expression levels of generic metabolic (e.g. glycolysis, nucleotide metabolism) or cellular (e.g. DNA replication) pathways in the ring and trophozoite stages but strong increased expression of essentially all functionalities associated with protein synthesis, folding and trafficking in the schizont stages. This specific "tune-up" of the transcriptional pattern in these resistant P. falciparum isolates is associated with altered expression of a number of genes involved in cell cycle regulation, transcription regulation, chromatin remodeling as well as intracellular signaling. Altogether these results provide the first set of testable genetic markers associated with artemisinin resistance.
Gene expression associated with artemisinin resistance
The main purpose of these studies was to characterize the transcriptional profile associated with artemisinin resistance in field isolates of P. falciparum . For this, we conducted DNA microarray analyses of parasites collected from patients in Pailin, Western Cambodia which had slow parasite clearance (detectable parasitemia 78-96 hours following administration of ACTs; isolates CP025, CP037 and CP040; Additional file 1). For comparison we analyzed transcriptomes of additional Southeast Asian isolates collected from Xepon, Savannakhet Province, Laos (isolates BMT061, BMT076, BMT077, XPN003), from Mae Sot, Thailand (isolates NHP2094, NHP4459, NHP4460), and one additional isolate from Pailin (CP022) with normal clearance (~54 hour). At the time of collection, all isolates exhibited high synchronicity; ~100% parasites were at the ring stage (Additional file 2). Parasites were then cultured for up to 48 hours in vitro and total RNA was isolated from samples harvested at regular intervals of 2-8 hours.
To identify genes whose expression is associated with artemisinin resistance, we compared global transcription levels between the three artemisinin resistant isolates and the remaining sensitive isolates by analyzing steady state mRNA levels in 3 developmental stages; ring (12-16 hpi: hours post invasion), trophozoite (24-28 hpi) and schizont (32-36 hpi) (Figure 1b). For that, we included the experimental time points that show highest correlation value to the in vitro reference transcriptome  within these sampling intervals (Figure 1b). Overall, we identified 160 (3.9%), 198 (4.9%) and 373 (9.2%) genes that are up-regulated and 278 (6.9%), 344 (8.5%) and 232 (5.7%) genes that are down-regulated in the resistant parasites (p-value < 0.01) at the ring, trophozoite and schizont stages, respectively (Figure 1c, Additional file 5). Interestingly, the majority of the genes over-expressed in the ring stage of the resistant parasites correspond to trophozoite- specific genes (peak mRNA abundance in the normal IDC) and the genes over-expressed in the schizont stages represent mainly ring and trophozoite specific transcripts (Figure 1c). The genes up-regulated in the trophozoite stage are evenly distributed between ring and schizont specific transcripts. This indicates that the transcriptional up-regulations associated with artemisinin resistance constitute either accelerated timing of transcription of trophozoite and schizont genes at the ring stage, or prolonged expression of ring and trophozoite transcripts in the schizont stage. Intriguingly, this pattern was not observed for down-regulated genes whose expression in artemisinin resistant parasites coincides with their expected stage specificity, albeit at lower levels than expected in the artemisinin sensitive isolates (Figure 1c). Taken together, these results suggest that artemisinin resistance is associated with specific modifications of the IDC transcriptional cascade that involve a large number of genes and presumably alter the levels and temporal distributions of biological and cellular functions.
Functional analysis of artemisinin resistance associated genes
In summary, our data suggest that the specific modulation of the IDC transcriptional cascade observed in artemisinin resistant parasites has at least two major physiological implications that include: down-regulation of metabolic and cellular pathways in the first half of the IDC (up to 28 hpi) and prolonged transcription up-regulation of functionalities associated with protein synthesis and their supporting activities in the late stages (~36 hpi). Interestingly, this transcriptional program was observed in all three Western Cambodian isolates with the delayed clearance but not in the isolate collected in the same region (CP022) that had a normal clearance time (Figure 2b; Additional file 7). This suggests that the altered transcriptional pattern is not a simple reflection of geographical or other generic differences between P. falciparum isolates but is indeed associated with artemisinin resistance.
Differential expression of regulatory proteins may contribute to artemisinin resistance
Functional enrichment analyses of the top 5% of genes from both extremes of the distribution (see Figure 3a grey boxes) did not uncover any additional functional groups compared to the stage-wise GSEA (Figure 2), however, visual inspection of these gene groups identified 41 putative transcription regulators, of which 22 are under-expressed and 19 are over-expressed. These include transcription factors, RNA binding proteins, cell cycle regulators, chromatin remodeling associated proteins, histone modification enzymes as well as histones (Figure 3b, c, for complete list see Additional file 5 Sheet "Rank product score"). The most remarkable examples are strong down-regulation of genes coding for histone deacetylase 1 protein (HDAC1) (PFI1260c), recently proposed to play a major role in the transcriptional regulation , and CDK-activating kinase assembly factor (PFE0610c) whose paralogue was shown be associated with HDAC1 in P. falciparum . Also, the artemisinin resistant parasites show a considerable up-regulation of the AP2 containing transcription factor (PF14_0079), bromo-domain containing protein (PFA0510w), putative histone acetyltransferase (hat; PF14_0350), zinc finger protein (PF14_0416), and the P. falciparum homologue of yeast histone chaperone Rttp106-like transcriptional regulator (PFE0870w). In addition to these DNA interacting regulatory proteins, we identified 10 genes encoding the RNA-recognition motif (RRM) that are involved in regulating mRNA stability or translational repression. Quantitative real-time PCR carried out with 4 of these genes (PFA0510w, PF14_0079, PF11_0061 and PF14_0057) further confirmed that these genes are significantly differentially expressed between the resistant and sensitive isolates (Figure 3d). Differential expression of these genes encoding regulatory factors may contribute to the global changes in the transcriptome observed in the artemisinin resistant parasites.
Copy number variations (CNV) and genotypes of artemisinin resistant isolates
In addition to the CGH analysis, we carried out "standard" genotyping of genes coding for merozoite surface proteins: msp1, msp2 and glutamate-rich protein: glurp together with sequencing of drug-resistant markers: multidrug resistant protein 1 (pfmdr1) and dihydropteroate synthase (pfdhps). This analysis revealed that the four Cambodian and two Lao isolates are genetically distinct from each other carrying different alleles of the tested genes (Figure 4b). Of these, CP025 and CP037 appeared highly related, exhibiting identical msp1 and msp2 nested PCR products of KO33 and IC/3d7 alleles (Figure 4b). However, these two isolates differ in 2 non-synonymous SNPs of codons K214Q and K1181T and a synonymous SNP in codon 171 in the pfmdr1 gene (Figure 4c). In addition, CP025 and CP037 exhibit subtle differences in their CNV pattern as measured by CGH (Figure 4a). Interestingly, we found that none of the 6 isolates carry the N86Y pfmdr1 mutation (Figure 4c) whose presence is frequent in many Southeast Asian regions  and which is commonly associated with chloroquine resistance and mefloquine hypersensitivity . Furthermore, we observed significant haplotype sharing amongst analyzed isolates such as CP040 and CP022 that share the "NEFSNVD" haplotype of pfmdr1, containing the non-synonymous Y184F mutation. Moreover, CP025, CP037, CP040, CP022 and XPN003 share the IAGEAA haplotype of pfdhps with non-synonymous triple mutants S436A, G437A and K540E, which were previously linked to resistance to sulfadoxine-pyrimethamine [38, 39]. The mutations in dhps may reflect the selection pressure as a result of an extensive antifolate drug use in this region over several decades. In summary, these genetic studies suggest that while the Western Cambodian isolates are not identical isogenic clones, they share a recent common ancestor that is evident from their CNV profiles. Hence, the highly unique and also uniform CNV pattern identified in these parasites may represent a genetic background that contributes to development of artemisinin resistance and possibly other drug resistant phenotypes in this region.
Mechanism of artemisinin resistance
Artemisinin resistance of P. falciparum is a major threat to malaria control. Understanding its molecular basis is thus essential for determining treatment strategies, mapping the spread of resistance and guiding elimination [20, 40–44]. Mutations or amplification of genes encoding transporters or target enzymes have been identified as resistance mechanisms to other antimalarial drugs. It is possible that resistance to artemisinin is unlike these classical mechanisms in vitro but instead results from a complex series of genetic and epigenetic events affecting multiple pathways. Here we showed that artemisinin resistant parasites are characterized by a specific modulation of the IDC transcriptome that affects a broad spectrum of genes and biological functions. First, we detected a specific down-regulation of many ring stage-specific metabolic pathways such as energy metabolism, nucleotide synthesis etc. Second, the progress to mature schizont stages was marked by significant increase in transcription of genes associated with essentially all functionalities involved in protein metabolism (Figure 2). Based on these two observations, we propose a plausible model for molecular basis of the artemisinin resistance observed in these Cambodian isolates that is founded on two (possibly independent) events: (i) Limited metabolic activity of the ring stages could lead to lower levels of hemoglobin digestion, less ferrous ions produced and thus reduced activation and conversion of artemisinin drugs to reactive intermediates . This effect may be linked to the high prevalence of hemoglobinopathies, especially Hemoglobin E and α-thalassemia in this area conferring increased oxidant stress [46, 47]. This is possibly compounded by decades of unregulated artemisinin use in Western Cambodia which may select for parasites with greater resistance to oxidant stress during the ring form development. Alternatively, the lower metabolic activities of the ring stages may be a prerequisite of the ability of the resistant parasite to become dormant under the artemisinin pressure as demonstrated by several in vitro analyses [17, 18]. Although this phenomenon is unlikely to be the major factor in artemisinin resistance, it is worth noting that the results obtained in our study are in agreement with the transcriptional analyses of the in vitro derived artemisinin tolerant parasites that include mRNA levels of the hypoxanthine phosphoribosyl transferase (PF10_0121), a cell cycle regulator (PFE1415w) and a Heat Shock Protein 70 kDa (PF08_0054) (Additional file 9) . (ii) The increased activity of protein metabolism may contribute to the resistance by withstanding the damaging effects of artemisinin on parasite's proteins by increasing rates of protein turnover and synthesis and thus compensate for loss of active proteins. This is consistent with the presumed artemisinin mode of action in which the drug inflicts substantial damage to the parasite cell either by an oxidative stress (via reactive free radicals) or by direct alkylation of a wide spectrum of cellular components such as proteins, heme, and lipids [48–53].
Transcriptional regulation and artemisinin resistance
It is important to note that the differences in mRNA levels linked to artemisinin resistance were not the most dramatic transcriptional differences detected among the field isolates. Functional enrichment analysis of genes that vary the most in their expression among the 11 field isolates irrespective of their artemisinin sensitivities, revealed a distinct set of functionalities that varied in gene expression amongst field isolates in general. These include mainly genes encoding proteins of host-parasite interactions, Maurer's clefts proteins, invasion-related proteins, as well as factors of fatty acid metabolism, glycolysis and hemoglobin digestion (Additional file 10). This indicates that the observed transcriptional differences between artemisinin resistant and sensitive strains may represent a new event in transcriptional differentiation of field isolates that does not occur in other parasite populations. Although very little is known about molecular factors of transcriptional regulation in P. falciparum, a number of previous studies [25, 54–57] have indicated the existence of broad range regulatory mechanisms that affect large groups of genes controlling the progression of the Plasmodium life cycle. Although more studies that involve large epidemiological surveys followed by extensive molecular analyses will have to be carried out to support such putative mechanisms, here we explored genes whose expression was affected in all three developmental stages in the resistant Western Cambodian parasites. One of the most remarkable examples of consistently down-regulated genes is pfhdac1 and its putative interacting proteins partner CDK-activating kinase, both of which show dramatic down-regulation in all three developmental stages (Figure 3). In our recent study, we showed that HDAC activities play a highly dynamic role in regulating the transcriptome of the Plasmodium life cycle and that their inhibition leads to widespread transcriptional changes. Intriguingly, the transcription profile of the resistant parasites is reminiscent of the effect of HDAC inhibitors on the P. falciparum gene expression in which most gene induction represent accelerated transcription in the ring and trophozoite stages and prolonged mRNA expression during the schizont stage .
Another gene whose differential expression is significant is histone 4 (PF11_0061). Histone 4 is the only nucleosomal subunit present as a single copy gene and thus its down-regulation could have severe implication on the nucleosomal assembly during the P. falciparum schizont stage. Depletion of histone 4 in both yeast and human cells causes an arrest or dramatic delay in the progression of the S-phase, which is caused by insufficient nucleosomal assembly and subsequently inhibition of DNA synthesis [58, 59]. We observed a modest but statistically significant decrease of histone 4 in the artemisinin resistant P. falciparum parasites (Figure 3d) that could analogously lead to a delayed onset of schizogony and thus prolonged expression of genes associated with the typical trophozoite functions such as protein synthesis and hemoglobin digestion (Figure 2). Since the CNV profiles were similar among the four Cambodian isolates regardless of their tolerance level to artemisinin, more studies will be needed to understand their role in regulating transcriptional levels in artemisinin resistant parasites. More importantly, the coherence of CNV profiles among the Cambodian isolates indicates that the parasites in this region share a recent common ancestor. Similar conclusion was made by an admixture study carried out in Southeast Asia that demonstrated a sub-group of recently expanded parasite populations in Western Cambodia . This could be a result of intense administration of a successive number of antimalarials from chloroquine to mefloquine to sulfodoxine-pyrimethamine to treat malaria cases in the Western Cambodian region over the last 80 years. It is plausible that a population with a restricted genotype has emerged with a genetic background that has higher propensity to give rise to resistance of the P. falciparum parasites to artemisinin.
Here we present the most comprehensive transcriptional analysis of the artemisinin tolerant malaria parasites collected at the site of its origin, Western Cambodia. Although much more extensive research on multiple levels will be needed to elucidate the complete mechanism, the results from these studies provide the first testable hypotheses regarding both the global pattern of the transcriptional cascade and candidate genes involved in artemisinin resistance.
Sample collection from patient and during ex-vivo IDC
On admission to the hospital and upon consent, 10 ml of total infected blood was taken from each adult P. falciparum-infected patient from the various field sites (Pailin, Cambodia; Xepon, Savannakhet Province, Laos; Mae Sot, Thailand). Samples were filtered by CF11 column purification to remove white blood cells and obtain infected red blood cells which were subjected to in vitro culture for one generation of the IDC over 48 hours. For every isolate, a total of 6 to 10 samples were harvested at regular intervals of 2 to 8 hours throughout the IDC. 100-500 μl of packed red blood cells with 1% - 10% parasitemia were collected per time point sampling.
RNA extraction, cDNA synthesis & Microarray Hybridizations
Total RNA was isolated from all ex-vivo samples using the trizol-chloroform-isopropanol precipitation method as described . Synthesis of target DNA for microarray hybridization was carried out as previously described [27, 61] with the following exceptions: Reverse-transcription was carried out using 200 units of SuperScript II enzyme (Invitrogen, USA) to create the first strand cDNA. Following PCR amplification with 19 cycles and purification of the product, a total of 5-7_g of double-stranded DNA was obtained for each sample. From this, 4 μg of DNA was labeled with fluorescent Cy5 dye and used for the microarray hybridization (GE Amersham, USA). For the reference pool, cDNA consisting of all IDC stages of 3d7 strain was made identically and labeled with Cy3 dye. Equal amounts of labeled sample from each time point and reference RNA pool was hybridized on the 70-mer P. falciparum cDNA microarray chip containing 10,680 oligonucleotides representing 5,343 coding genes . The microarray hybridization was carried out at 65°C in the automated MAUI hybridization system (BioMicro Systems, USA). The microarray chips were scanned using the GenePix 4000 B scanner and GenePix Pro 6.0 software (Molecular Devices, USA).
Microarray data Processing and analysis
The initial acquired microarray data was filtered to include microarray hybridization "spots" with at least 95% pixels having signal intensity above 2 standard deviations from background for both Cy3 and Cy5 fluorescence intensity. Subsequently, array Lowess normalization was applied to all arrays as implemented by Acuity 4.0 software (Molecular Devices, USA). The assembly of the IDC transcriptome for all isolates was carried out using Fast Fourier Transform as previously described  with following modifications; the expression value for each gene is represented by an average value obtained from all oligonucleotide probes in each transcript/gene. To construct the phaseogram of reference IDC transcriptome for further comparisons, genes with expression data in at least 19 out of 24 (80%) time points were extracted, any missing values were imputed and polynomial fit of order 6 was applied to generate smoothed data. This was then analyzed by Fast Fourier Transform method and all 4,634 genes with the mean-centered gene expression log2 ratios were sorted according to phase from -π to π (Additional file 2). The phaseogram of each isolate contains mean-centered gene expression data ordered according to the phase value. For each isolate's time course, only genes having at least 80% of data present were extracted and any missing data values were imputed by KNN (Kth nearest neighbor) method using R version 2.10.1 package and impute function.
Mapping parasite's sample age to the reference IDC
To map the age of each isolate time point sample relative to the progression of in vitro lab strain IDC, Spearman Rank Correlation Coefficient (SRCC) values between global mRNA profiles for each isolate time point and time points in the reference IDC transcriptome  (every 2 hour sample time point of the in vitro Dd2 lifecycle) were calculated. The stage (hpi) corresponding to the peak SRCC value was assigned as the best estimate of the age of the parasite for that sample collection (see colored box in Figure 1b). Only samples which had SRCC above a value of 0.35 were considered in the following analysis. To select time points for downstream differential analysis, we calculated the frequency distribution of isolate time points found within a window of three consecutive hpi and selected three windows based on the highest representation of isolates with at least one sampling time point found within that window. Isolates with sample time points found within the 4 hour windows of 12-16 hpi, 24-28 hpi and 32-36 hpi corresponding to mid-ring, early-trophozoite and early-schizont stages were grouped accordingly as representing these stages.
Identifying Differentially Expressed Genes
Where μ a and μ b are the mean expression log2 ratio of resistant and sensitive parasites respectively and σ a and σ b are the standard deviations of resistant and sensitive group respectively. Differentially expressed genes were defined with p-value cut-off < 0.01 from Student t- distribution of the z-scores carried out for the 3 stages. Gene set enrichment analysis  was carried out on the resistant vs. sensitive parasites using SNR. A total of 1,261 gene sets comprising of KEGG pathways , MPM and functional groups based on previous studies [21, 64–66] were included. Nominal p-value and false discovery rate (FDR) were calculated for each gene set from the observed normalized enrichment score (NES) against a null distribution using a gene set-based permutation test. Gene sets with p-value < 0.05 and FDR < 25% were considered to be statistically significant.
Calculate Timing of Peak mRNA abundance of genes
Where θ is the phase value of the gene.
Calculation of Rank Product Score
Where the rankg14hpi, rankg26hpi, rankg34hpi corresponds to the rank score of the gene at 14 hpi, 26 hpi or 34 hpi respectively.
Copy Number Difference Analysis
To identify significant CNVs among the 6 isolates, we performed comparative genomic hybridization using genomic DNA extracted from the isolates (CP022, CP025, CP037, CP040, BMT077, XPN003) against genomic DNA of the reference lab strain of 3d7 as previously described . We carried out subtraction analysis between any 2 isolates in multiple pair-wise comparisons and considered only segments that have at least 1.7-fold change in mean amplitude and T-coefficient greater than 3.5 in minimum 2 consecutive gene probes between any 2 isolates as significant. This was performed using the Genomic Alteration Detection Analysis (GADA) program in R package .
Real-time PCR (qPCR)
Real time PCR and the relative quantification of gene expression was performed in at least triplicates using the comparative CT method of calculating -ΔΔCt values and taking PFC0965w as a reference control gene, as previously described .
Genotyping of isolates and Sequencing of pfmdr1 and pfdhps genes
Genotyping of msp1, msp2 and glurp were carried out as previously described  with positive control of lab clone, 3d7, in a nested PCR reaction and products were visualized on a 2.5% agarose gel. Nested PCR reactions were performed as previously described  to generate fragments of mdr1 and dhps genes containing commonly found SNPs and products were sequenced with corresponding forward and reverse primers using ABI Big Dye Terminator cycle sequencing method.
List of abbreviations used
Intra-erythrocytic developmental cycle
hours post invasion
Parasite clearance time
Spearman rank correlation coefficient
Pearson correlation coefficient
Signal to noise ratio
Normalized enrichment score
Gene set enrichment analysis
Kyoto Encyclopedia of Genes and Genomes
Malaria metabolic pathway maps
Copy number variants
Polymerase chain reaction
Single nucleotide polymorphism
Maniphone Khanthavong, Odai Chanthongthip, Samlane Phompida, Tiengkham Pongvongsa, and Sengchanh Koumphiengmaniseng for their technical help and support. This work was supported by grant of the Singapore National Medical Council (NMRC) # IRG10may058.
- Nosten F, White NJ: Artemisinin-based combination treatment of falciparum malaria. Am J Trop Med Hyg. 2007, 77: 181-192.PubMed
- White NJ: Qinghaosu (artemisinin): the price of success. Science. 2008, 320: 330-334. 10.1126/science.1155165.PubMedView Article
- Mendis K, Rietveld A, Warsame M, Bosman A, Greenwood B, Wernsdorfer WH: From malaria control to eradication: The WHO perspective. Trop Med Int Health. 2009, 14: 802-809. 10.1111/j.1365-3156.2009.02287.x.PubMedView Article
- Bruce-Chwatt LJ: Malaria and its control: present situation and future prospects. Annu Rev Public Health. 1987, 8: 75-110. 10.1146/annurev.pu.08.050187.000451.PubMedView Article
- Wellems TE: Transporter of a malaria catastrophe. Nat Med. 2004, 10: 1169-1171. 10.1038/nm1104-1169.PubMedView Article
- Wellems TE, Hayton K, Fairhurst RM: The impact of malaria parasitism: from corpuscles to communities. J Clin Invest. 2009, 119: 2496-2505. 10.1172/JCI38307.PubMed CentralPubMedView Article
- Wongsrichanalai C, Pickard AL, Wernsdorfer WH, Meshnick SR: Epidemiology of drug-resistant malaria. Lancet Infect Dis. 2002, 2: 209-218. 10.1016/S1473-3099(02)00239-6.PubMedView Article
- Yang H, Liu D, Yang Y, Fan B, Yang P, Li X, Li C, Dong Y, Yang C: Changes in susceptibility of Plasmodium falciparum to artesunate in vitro in Yunnan Province, China. Trans R Soc Trop Med Hyg. 2003, 97: 226-228. 10.1016/S0035-9203(03)90127-1.PubMedView Article
- Huong NM, Hewitt S, Davis TM, Dao LD, Toan TQ, Kim TB, Hanh NT, Phuong VN, Nhan DH, Cong LD: Resistance of Plasmodium falciparum to antimalarial drugs in a highly endemic area of southern Viet Nam: a study in vivo and in vitro. Trans R Soc Trop Med Hyg. 2001, 95: 325-329. 10.1016/S0035-9203(01)90254-8.PubMedView Article
- Jambou R, Legrand E, Niang M, Khim N, Lim P, Volney B, Ekala MT, Bouchier C, Esterre P, Fandeur T: Resistance of Plasmodium falciparum field isolates to in-vitro artemether and point mutations of the SERCA-type PfATPase6. The Lancet. 2005, 366: 1960-1963. 10.1016/S0140-6736(05)67787-2.View Article
- Kachur SP, MacArthur JR, Slutsker L: A call to action: addressing the challenge of artemisinin-resistant malaria. Expert Rev Anti Infect Ther. 2010, 8: 365-366. 10.1586/eri.10.23.PubMedView Article
- Noedl H, Se Y, Schaecher K, Smith BL, Socheat D, Fukuda MM: Evidence of artemisinin-resistant malaria in western Cambodia. N Engl J Med. 2008, 359: 2619-2620. 10.1056/NEJMc0805011.PubMedView Article
- Carrara VI, Zwang J, Ashley EA, Price RN, Stepniewska K, Barends M, Brockman A, Anderson T, McGready R, Phaiphun L, Proux S, van Vugt M, Hutagalung R, Lwin KM, Phyo AP, Preechapornkul P, Imwong M, Pukrittayakamee S, Singhasivanon P, White NJ, Nosten F: Changes in the treatment responses to artesunate-mefloquine on the northwestern border of Thailand during 13 years of continuous deployment. PLoS One. 2009, 4: e4551-10.1371/journal.pone.0004551.PubMed CentralPubMedView Article
- Lim P, Wongsrichanalai C, Chim P, Khim N, Kim S, Chy S, Sem R, Nhem S, Yi P, Duong S, Bouth DM, Genton B, Beck HP, Gobert JG, Rogers WO, Coppee JY, Fandeur T, Mercereau-Puijalon O, Ringwald P, Le Bras J, Ariey F: Decreased in vitro susceptibility of Plasmodium falciparum isolates to artesunate, mefloquine, chloroquine, and quinine in Cambodia from 2001 to 2007. Antimicrob Agents Chemother. 2010, 54: 2135-2142. 10.1128/AAC.01304-09.PubMed CentralPubMedView Article
- Dondorp AM, Nosten F, Yi P, Das D, Phyo AP, Tarning J, Lwin KM, Ariey F, Hanpithakpong W, Lee SJ, Ringwald P, Silamut K, Imwong M, Chotivanich K, Lim P, Herdman T, An SS, Yeung S, Singhasivanon P, Day NP, Lindegardh N, Socheat D, White NJ: Artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2009, 361: 455-467. 10.1056/NEJMoa0808859.PubMed CentralPubMedView Article
- Saralamba S, Pan-Ngum W, Maude RJ, Lee SJ, Tarning J, Lindegardh N, Chotivanich K, Nosten F, Day NP, Socheat D, White NJ, Dondorp AM, White LJ: Intrahost modeling of artemisinin resistance in Plasmodium falciparum. Proc Natl Acad Sci USA. 2011, 108: 397-402. 10.1073/pnas.1006113108.PubMed CentralPubMedView Article
- Teuscher F, Gatton ML, Chen N, Peters J, Kyle DE, Cheng Q: Artemisinin-induced dormancy in plasmodium falciparum: duration, recovery rates, and implications in treatment failure. J Infect Dis. 2010, 202: 1362-1368. 10.1086/656476.PubMed CentralPubMedView Article
- Witkowski B, Lelievre J, Lopez Barragan MJ, Laurent V, Su Xz, Berry A, Benoit-Vical F: Increased Tolerance to Artemisinin in Plasmodium falciparum Is Mediated by a Quiescence Mechanism. Antimicrobial Agents and Chemotherapy. 2010, 54: 1872-1877. 10.1128/AAC.01636-09.PubMed CentralPubMedView Article
- Nosten F: Waking the sleeping beauty. J Infect Dis. 2010, 202: 1300-1301. 10.1086/656478.PubMed CentralPubMedView Article
- Anderson Tim JC, Nair S, Nkhoma S, Williams Jeff T, Imwong M, Yi P, Socheat D, Das D, Chotivanich K, Day Nicholas PJ, White Nicholas J, Dondorp Arjen M: High Heritability of Malaria Parasite Clearance Rate Indicates a Genetic Basis for Artemisinin Resistance in Western Cambodia. The Journal of Infectious Diseases. 2010, 201: 1326-1330. 10.1086/651562.PubMed CentralPubMedView Article
- Bozdech Z, Llinas M, Pulliam BL, Wong ED, Zhu J, DeRisi JL: The transcriptome of the intraerythrocytic developmental cycle of Plasmodium falciparum. PLoS Biol. 2003, 1: E5-PubMed CentralPubMedView Article
- Foth BJ, Zhang N, Chaal BK, Sze SK, Preiser PR, Bozdech Z: Quantitative time-course profiling of parasite and host cell proteins in the human malaria parasite Plasmodium falciparum. Mol Cell Proteomics. 2011
- Lemieux JE, Gomez-Escobar N, Feller A, Carret C, Amambua-Ngwa A, Pinches R, Day F, Kyes SA, Conway DJ, Holmes CC, Newbold CI: Statistical estimation of cell-cycle progression and lineage commitment in Plasmodium falciparum reveals a homogeneous pattern of transcription in ex vivo culture. Proc Natl Acad Sci USA. 2009, 106: 7559-7564. 10.1073/pnas.0811829106.PubMed CentralPubMedView Article
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102: 15545-15550. 10.1073/pnas.0506580102.PubMed CentralPubMedView Article
- Chaal BK, Gupta AP, Wastuwidyaningtyas BD, Luah YH, Bozdech Z: Histone deacetylases play a major role in the transcriptional regulation of the Plasmodium falciparum life cycle. PLoS Pathog. 2010, 6: e1000737-10.1371/journal.ppat.1000737.PubMed CentralPubMedView Article
- Halbert J, Ayong L, Equinet L, Le Roch K, Hardy M, Goldring D, Reininger L, Waters N, Chakrabarti D, Doerig C: A Plasmodium falciparum transcriptional cyclin-dependent kinase-related kinase with a crucial role in parasite proliferation associates with histone deacetylase activity. Eukaryot Cell. 2010, 9: 952-959. 10.1128/EC.00005-10.PubMed CentralPubMedView Article
- Mackinnon MJ, Li J, Mok S, Kortok MM, Marsh K, Preiser PR, Bozdech Z: Comparative transcriptional and genomic analysis of Plasmodium falciparum field isolates. PLoS Pathog. 2009, 5: e1000644-10.1371/journal.ppat.1000644.PubMed CentralPubMedView Article
- Dharia NV, Sidhu AB, Cassera MB, Westenberger SJ, Bopp SE, Eastman RT, Plouffe D, Batalov S, Park DJ, Volkman SK, Wirth DF, Zhou Y, Fidock DA, Winzeler EA: Use of high-density tiling microarrays to identify mutations globally and elucidate mechanisms of drug resistance in Plasmodium falciparum. Genome Biol. 2009, 10: R21-10.1186/gb-2009-10-2-r21.PubMed CentralPubMedView Article
- Cheeseman IH, Gomez-Escobar N, Carret CK, Ivens A, Stewart LB, Tetteh KK, Conway DJ: Gene copy number variation throughout the Plasmodium falciparum genome. BMC Genomics. 2009, 10: 353-10.1186/1471-2164-10-353.PubMed CentralPubMedView Article
- Kidgell C, Volkman SK, Daily J, Borevitz JO, Plouffe D, Zhou Y, Johnson JR, Le Roch K, Sarr O, Ndir O, Mboup S, Batalov S, Wirth DF, Winzeler EA: A systematic map of genetic variation in Plasmodium falciparum. PLoS Pathog. 2006, 2: e57-10.1371/journal.ppat.0020057.PubMed CentralPubMedView Article
- Jiang H, Yi M, Mu J, Zhang L, Ivens A, Klimczak LJ, Huyen Y, Stephens RM, Su XZ: Detection of genome-wide polymorphisms in the AT-rich Plasmodium falciparum genome using a high-density microarray. BMC Genomics. 2008, 9: 398-10.1186/1471-2164-9-398.PubMed CentralPubMedView Article
- Ribacke U, Mok BW, Wirta V, Normark J, Lundeberg J, Kironde F, Egwang TG, Nilsson P, Wahlgren M: Genome wide gene amplifications and deletions in Plasmodium falciparum. Mol Biochem Parasitol. 2007, 155: 33-44. 10.1016/j.molbiopara.2007.05.005.PubMedView Article
- Uhlemann AC, McGready R, Ashley EA, Brockman A, Singhasivanon P, Krishna S, White NJ, Nosten F, Price RN: Intrahost selection of Plasmodium falciparum pfmdr1 alleles after antimalarial treatment on the northwestern border of Thailand. J Infect Dis. 2007, 195: 134-141. 10.1086/509809.PubMed CentralPubMedView Article
- Price RN, Uhlemann AC, Brockman A, McGready R, Ashley E, Phaipun L, Patel R, Laing K, Looareesuwan S, White NJ, Nosten F, Krishna S: Mefloquine resistance in Plasmodium falciparum and increased pfmdr1 gene copy number. Lancet. 2004, 364: 438-447. 10.1016/S0140-6736(04)16767-6.PubMed CentralPubMedView Article
- Nair S, Miller B, Barends M, Jaidee A, Patel J, Mayxay M, Newton P, Nosten F, Ferdig MT, Anderson TJ: Adaptive copy number evolution in malaria parasites. PLoS Genet. 2008, 4: e1000243-10.1371/journal.pgen.1000243.PubMed CentralPubMedView Article
- Price RN, Cassar C, Brockman A, Duraisingh M, van Vugt M, White NJ, Nosten F, Krishna S: The pfmdr1 gene is associated with a multidrug-resistant phenotype in Plasmodium falciparum from the western border of Thailand. Antimicrobial Agents and Chemotherapy. 1999, 43: 2943-2949.PubMed CentralPubMed
- Mita T, Kaneko A, Hombhanje F, Hwaihwanje I, Takahashi N, Osawa H, Tsukahara T, Masta A, Lum JK, Kobayakawa T, Ishizaki T, Bjorkman A: Role of pfmdr1 mutations on chloroquine resistance in Plasmodium falciparum isolates with pfcrt K76T from Papua New Guinea. Acta Trop. 2006, 98: 137-144. 10.1016/j.actatropica.2006.03.002.PubMedView Article
- Enosse S, Magnussen P, Abacassamo F, Gomez-Olive X, Ronn AM, Thompson R, Alifrangis M: Rapid increase of Plasmodium falciparum dhfr/dhps resistant haplotypes, after the adoption of sulphadoxine-pyrimethamine as first line treatment in 2002, in southern Mozambique. Malar J. 2008, 7: 115-10.1186/1475-2875-7-115.PubMed CentralPubMedView Article
- Babiker HA, Al-Saai S, Kheir A, Abdel-Muhsin AMA, Al-Ghazali A, Nwakanma D, Swedberg G: Distinct haplotypes of dhfr and dhps among Plasmodium falciparum isolates in an area of high level of sulfadoxine-pyrimethamine (SP) resistance in eastern Sudan. Infection Genetics and Evolution. 2009, 9: 778-783. 10.1016/j.meegid.2009.04.010.View Article
- Witkowski B, Lelievre J, Lopez Barragan MJ, Laurent V, Su XZ, Berry A, Benoit-Vical F: Increased tolerance to artemisinin in Plasmodium falciparum is mediated by a quiescence mechanism. Antimicrob Agents Chemother. 2010
- Imwong M, Dondorp AM, Nosten F, Yi P, Mungthin M, Hanchana S, Das D, Phyo AP, Lwin KM, Pukrittayakamee S, Lee SJ, Saisung S, Koecharoen K, Nguon C, Day NP, Socheat D, White NJ: Exploring the contribution of candidate genes to artemisinin resistance in Plasmodium falciparum. Antimicrob Agents Chemother. 2010, 54: 2886-2892. 10.1128/AAC.00032-10.PubMed CentralPubMedView Article
- Chen N, Chavchich M, Peters JM, Kyle DE, Gatton ML, Cheng Q: Deamplification of pfmdr1-Containing Amplicon on Chromosome 5 in Plasmodium falciparum Is Associated with Reduced Resistance to Artelinic Acid In Vitro. Antimicrobial Agents and Chemotherapy. 2010, 54: 3395-3401. 10.1128/AAC.01421-09.PubMed CentralPubMedView Article
- Hunt P, Afonso A, Creasey A, Culleton R, Sidhu ABS, Logan J, Valderramos SG, McNae I, Cheesman S, Rosario Vd, Carter R, Fidock DA, Cravo P: Gene encoding a deubiquitinating enzyme is mutated in artesunate- and chloroquine-resistant rodent malaria parasites. Molecular Microbiology. 2007, 65: 27-40. 10.1111/j.1365-2958.2007.05753.x.PubMed CentralPubMedView Article
- Hunt P, Martinelli A, Modrzynska K, Borges S, Creasey A, Rodrigues L, Beraldi D, Loewe L, Fawcett R, Kumar S, Thomson M, Trivedi U, Otto TD, Pain A, Blaxter M, Cravo P: Experimental evolution, genetic analysis and genome re-sequencing reveal the mutation conferring artemisinin resistance in an isogenic lineage of malaria parasites. BMC Genomics. 2010, 11: 499-10.1186/1471-2164-11-499.PubMed CentralPubMedView Article
- Klonis N, Crespo-Ortiz MP, Bottova I, Abu-Bakar N, Kenny S, Rosenthal PJ, Tilley L: Artemisinin activity against Plasmodium falciparum requires hemoglobin uptake and digestion. Proc Natl Acad Sci USA. 2011, 108: 11405-11410. 10.1073/pnas.1104063108.PubMed CentralPubMedView Article
- Flint J, Harding RM, Boyce AJ, Clegg JB: The population genetics of the haemoglobinopathies. Baillieres Clin Haematol. 1998, 11: 1-51. 10.1016/S0950-3536(98)80069-3.PubMedView Article
- Wiwanitkit V: Genetic disorders and malaria in Indo-China region. J Vector Borne Dis. 2008, 45: 98-104.PubMed
- Meshnick SR, Taylor TE, Kamchonwongpaisan S: Artemisinin and the antimalarial endoperoxides: from herbal remedy to targeted chemotherapy. Microbiol Rev. 1996, 60: 301-315.PubMed CentralPubMed
- Robert A, Bonduelle C, Laurent SAL, Meunier B: Heme alkylation by artemisinin and trioxaquines. J Phys Org Chem. 2006, 19: 562-569. 10.1002/poc.1059.View Article
- Yang YZ, Little B, Meshnick SR: Alkylation of proteins by artemisinin. Effects of heme, pH, and drug structure. Biochem Pharmacol. 1994, 48: 569-573.PubMedView Article
- Bhisutthibhan J, Meshnick SR: Immunoprecipitation of [(3)H]dihydroartemisinin translationally controlled tumor protein (TCTP) adducts from Plasmodium falciparum-infected erythrocytes by using anti-TCTP antibodies. Antimicrob Agents Chemother. 2001, 45: 2397-2399. 10.1128/AAC.45.8.2397-2399.2001.PubMed CentralPubMedView Article
- Bhisutthibhan J, Pan XQ, Hossler PA, Walker DJ, Yowell CA, Carlton J, Dame JB, Meshnick SR: The Plasmodium falciparum translationally controlled tumor protein homolog and its reaction with the antimalarial drug artemisinin. J Biol Chem. 1998, 273: 16192-16198. 10.1074/jbc.273.26.16192.PubMedView Article
- Asawamahasakda W, Ittarat I, Pu YM, Ziffer H, Meshnick SR: Reaction of antimalarial endoperoxides with specific parasite proteins. Antimicrob Agents Chemother. 1994, 38: 1854-1858.PubMed CentralPubMedView Article
- Hu G, Cabrera A, Kono M, Mok S, Chaal BK, Haase S, Engelberg K, Cheemadan S, Spielmann T, Preiser PR, Gilberger TW, Bozdech Z: Transcriptional profiling of growth perturbations of the human malaria parasite Plasmodium falciparum. Nat Biotechnol. 2010, 28: 91-98. 10.1038/nbt.1597.PubMedView Article
- Campbell TL, De Silva EK, Olszewski KL, Elemento O, Llinas M: Identification and Genome-Wide Prediction of DNA Binding Specificities for the ApiAP2 Family of Regulators from the Malaria Parasite. PLoS Pathog. 2010, 6: e1001165-10.1371/journal.ppat.1001165.PubMed CentralPubMedView Article
- Gonzales JM, Patel JJ, Ponmee N, Jiang L, Tan A, Maher SP, Wuchty S, Rathod PK, Ferdig MT, Templeton TJ: Regulatory Hotspots in the Malaria Parasite Genome Dictate Transcriptional Variation. PLoS Biology. 2008, 6: e238-10.1371/journal.pbio.0060238.PubMed CentralPubMedView Article
- Daily JP, Scanfeld D, Pochet N, Le Roch K, Plouffe D, Kamal M, Sarr O, Mboup S, Ndir O, Wypij D, Levasseur K, Thomas E, Tamayo P, Dong C, Zhou Y, Lander ES, Ndiaye D, Wirth D, Winzeler EA, Mesirov JP, Regev A: Distinct physiological states of Plasmodium falciparum in malaria-infected patients. Nature. 2007, 450: 1091-1095. 10.1038/nature06311.PubMedView Article
- Kim UJ, Han M, Kayne P, Grunstein M: Effects of histone H4 depletion on the cell cycle and transcription of Saccharomyces cerevisiae. EMBO J. 1988, 7: 2211-2219.PubMed CentralPubMed
- Prado F, Aguilera A: Partial Depletion of Histone H4 Increases Homologous Recombination-Mediated Genetic Instability. Molecular and Cellular Biology. 2005, 25: 1526-1536. 10.1128/MCB.25.4.1526-1536.2005.PubMed CentralPubMedView Article
- Mu J, Myers RA, Jiang H, Liu S, Ricklefs S, Waisberg M, Chotivanich K, Wilairatana P, Krudsood S, White NJ, Udomsangpetch R, Cui L, Ho M, Ou F, Li H, Song J, Li G, Wang X, Seila S, Sokunthea S, Socheat D, Sturdevant DE, Porcella SF, Fairhurst RM, Wellems TE, Awadalla P, Su XZ: Plasmodium falciparum genome-wide scans for positive selection, recombination hot spots and resistance to antimalarial drugs. Nat Genet. 2010, 42: 268-271. 10.1038/ng.528.PubMed CentralPubMedView Article
- Bozdech Z, Mok S, Hu G, Imwong M, Jaidee A, Russell B, Ginsburg H, Nosten F, Day NP, White NJ, Carlton JM, Preiser PR: The transcriptome of Plasmodium vivax reveals divergence and diversity of transcriptional regulation in malaria parasites. Proc Natl Acad Sci USA. 2008, 105: 16290-16295. 10.1073/pnas.0807404105.PubMed CentralPubMedView Article
- Hu G, Llinas M, Li J, Preiser PR, Bozdech Z: Selection of long oligonucleotides for gene expression microarrays using weighted rank-sum strategy. BMC Bioinformatics. 2007, 8: 350-10.1186/1471-2105-8-350.PubMed CentralPubMedView Article
- Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M: The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004, 32: D277-280. 10.1093/nar/gkh063.PubMed CentralPubMedView Article
- Ginsburg H: Malaria Parasite Metabolic Pathways. [http://sites.huji.ac.il/malaria/]
- Le Roch KG, Zhou Y, Blair PL, Grainger M, Moch JK, Haynes JD, De La Vega P, Holder AA, Batalov S, Carucci DJ, Winzeler EA: Discovery of gene function by expression profiling of the malaria parasite life cycle. Science. 2003, 301: 1503-1508. 10.1126/science.1087025.PubMedView Article
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.PubMed CentralPubMedView Article
- Pique-Regi R, Ortega A, Asgharzadeh S: Joint estimation of copy number variation and reference intensities on multiple DNA arrays using GADA. Bioinformatics. 2009, 25: 1223-1230. 10.1093/bioinformatics/btp119.PubMed CentralPubMedView Article
- Snounou G: Genotyping of Plasmodium spp. Nested PCR. Methods Mol Med. 2002, 72: 103-116.PubMed
- Crameri A, Marfurt J, Mugittu K, Maire N, Regos A, Coppee JY, Sismeiro O, Burki R, Huber E, Laubscher D, Puijalon O, Genton B, Felger I, Beck HP: Rapid microarray-based method for monitoring of all currently known single-nucleotide polymorphisms associated with parasite resistance to antimalaria drugs. J Clin Microbiol. 2007, 45: 3685-3691. 10.1128/JCM.01178-07.PubMed CentralPubMedView Article
- Ginsburg H: Progress in in silico functional genomics: the malaria Metabolic Pathways database. Trends Parasitol. 2006, 22: 238-240. 10.1016/j.pt.2006.04.008.PubMedView Article
- Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.PubMed CentralPubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.