Modulation of gene expression in drug resistant Leishmania is associated with gene amplification, gene deletion and chromosome aneuploidy
© Ubeda et al.; licensee BioMed Central Ltd. 2008
Received: 25 February 2008
Accepted: 18 July 2008
Published: 18 July 2008
Drug resistance can be complex, and several mutations responsible for it can co-exist in a resistant cell. Transcriptional profiling is ideally suited for studying complex resistance genotypes and has the potential to lead to novel discoveries. We generated full genome 70-mer oligonucleotide microarrays for all protein coding genes of the human protozoan parasites Leishmania major and Leishmania infantum. These arrays were used to monitor gene expression in methotrexate resistant parasites.
Leishmania is a eukaryotic organism with minimal control at the level of transcription initiation and few genes were differentially expressed without concomitant changes in DNA copy number. One exception was found in Leishmania major, where the expression of whole chromosomes was down-regulated. The microarrays highlighted several mechanisms by which the copy number of genes involved in resistance was altered; these include gene deletion, formation of extrachromosomal circular or linear amplicons, and the presence of supernumerary chromosomes. In the case of gene deletion or gene amplification, the rearrangements have occurred at the sites of repeated (direct or inverted) sequences. These repeats appear highly conserved in both species to facilitate the amplification of key genes during environmental changes. When direct or inverted repeats are absent in the vicinity of a gene conferring a selective advantage, Leishmania will resort to supernumerary chromosomes to increase the levels of a gene product.
Aneuploidy has been suggested as an important cause of drug resistance in several organisms and additional studies should reveal the potential importance of this phenomenon in drug resistance in Leishmania.
The protozoan parasite Leishmania is distributed worldwide and is responsible for a wide spectrum of diseases, including cutaneous, mucocutaneous and visceral leishmaniasis. No vaccines are presently available against Leishmania infections  and treatments rely primarily on chemotherapy. The chemotherapeutic arsenal is limited and resistance to the mainstay of pentavalent antimonials has reached epidemic proportions in parts of India . Several studies dealing with drug resistance in Leishmania have highlighted the plasticity of the Leishmania genome [3, 4]. The antifolate methotrexate (MTX) has been one of the first and most widely used drugs for understanding drug-induced plasticity and resistance mechanisms [5–8]. While Leishmania is sensitive to MTX, the drug is not used clinically to treat leishmaniasis. However, Leishmania is a folic acid auxotroph and studies of MTX resistance mechanisms have highlighted several novel aspects of folate metabolism in this parasite that could be exploited for drug interventions [9, 10]. Indeed, the development of novel antifolate molecules for Leishmania and related parasites has been ongoing in several laboratories [11–13].
Leishmania resists MTX by a number of mechanisms. Leishmania has the capacity to transport folic acid, but this activity is often impaired in MTX resistant cells [8, 14–17]. The main Leishmania folate transporter FT1 has been isolated [18, 19] and is part of a large family of folate biopterin transporter (FBT) proteins with 14 members in Leishmania (AA Ouameur et al., unpublished data). Rearrangements of FBT genes are correlated with MTX resistance [19–21]. A frequent mechanism of drug resistance in Leishmania is gene amplification . Small chromosomal regions of 20-70 kb that are part of one of the 36 Leishmania chromosomes are amplified as part of extrachromosomal elements . These elements are usually formed by recombination between repeated homologous sequences [22–24]. Amplification of the gene coding for the target dihydrofolate reductase-thymidylate synthase (DHFR-TS) has been described in MTX resistant parasites [5, 6, 25–29]. Work on MTX resistance also led to the characterization of the pteridine reductase PTR1, whose main function is to reduce pterins. However, when overexpressed it can also reduce folic acid and lead to MTX resistance by by-passing DHFR-TS activity [30–33]. The PTR1 gene is frequently amplified as part of extrachromosomal circular or linear amplicons [6, 16, 22, 34–38]. In addition to these three main mechanisms of resistance, perturbation in folate metabolism [39, 40], in one carbon metabolism  or in DNA metabolism  have also been associated with MTX resistance. Several of these mutations can co-exist in the same cell, demonstrating that resistance can be a complex multi-gene phenomenon. Genome wide expression profiling scans represent a useful tool for understanding complex resistance mechanisms and may lead either to the discovery of novel resistance mechanisms and/or could provide clues about mechanisms of gene rearrangements.
Indeed, DNA microarrays have been useful for investigating the mode of action of drugs  and mechanisms of resistance (reviewed in [44–46]). DNA microarrays for Leishmania have evolved from random genomic DNA clones [47–50], cDNA clones [51, 52], targeted PCR fragments , selected 70-mer oligonucleotides [53, 54] to full genome microarrays [55, 56]. Targeted microarrays have been used previously for the study of drug resistance in Leishmania [29, 52, 54, 57]. We present here the generation of full genome DNA microarrays for both L. major and L. infantum and their use in the study of one L. major and one L. infantum MTX resistant mutant. These genome wide expression profiling experiments illustrate the complexity of resistance mechanisms present in the same cell. They allowed the definition of the precise mechanisms leading to the formation of extrachromosomal circular and linear amplicons, the definition of gene deletion events and revealed the involvement of aneuploidy in the complex genotype of MTX resistance.
RNA expression profiling in methotrexate resistant Leishmaniacells
The DNA microarray data were supported by selected quantitative real-time reverse transcription PCR (qRT-PCR) assays in both the L. major and L. infantum mutants (Figure 3). In only two cases we found a discrepancy between the two techniques. LmjF04.0160 and its orthologue LinJ04_V3.0160 were found down-regulated in both mutants using DNA microarrays, but this was confirmed only in the L. major mutant by qRT-PCR (Figure 3). The other discrepancy between microarray and qRT-PCR data was for FT1, but this is explained by a gene deletion event (see below). The only other gene that was modulated similarly in the two mutants was the ABC protein gene ABCA2 and this was confirmed by qRT-PCR (Figure 3). Other genes were modulated in both mutants but in different ways. While LmjF31.0720 was down-regulated in L. major MTX60.4, its orthologue LinJ31_V3.0750 in L. infantum MTX20.5 was overexpressed (Figure 3). Otherwise, genes differentially expressed were specific to individual mutants.
The differential gene expression of the MTX resistant mutants was also represented in a chromosome by chromosome fashion (Figure 2). This has permitted us to visualize regions that are differently expressed (red/orange, corresponding to overexpressed genes in the mutants). Two regions were clearly overexpressed in the L. infantum MTX20.5 mutant. One region was on chromosome 6 (DHFR-TS loci) and the second was in the left portion of chromosome 23 (Figure 2a). For the L. major MTX60.4 mutant, we also saw an increase in expression of selected genes present on chromosome 6 (DHFR-TS loci), but we also observed a number of whole chromosomes (for example, chromosome 22; colored predominantly red in Figure 2b).
Extrachromosomal circular amplification of DHFR-TS
Linear amplification of PTR1
Decrease in gene expression due to deletion of folate transporter genes
Selection for MTX resistance and chromosome aneuploidy
A number of hypotheses can explain this whole chromosome-specific gene regulation and we tested whether the copy number of specific chromosomes changed upon MTX selection in L. major MTX60.4. Quantitative Southern blot analyses with two distinct probes derived from chromosome 22 revealed that if the wild-type cells contain two homologous copies of chromosome 22 (Leishmania is a diploid organism), L. major MTX60.4 had four copies (Figure 7a, lanes 1 and 2). Similarly, L. major MTX60.4 had three copies of chromosome 28 compared to wild-type cells (Figure 7b, lanes 1 and 2). The probes used are physically far apart, indicating a change in ploidy of the whole chromosome. However, this change in chromosome copy number was not observed for chromosomes 11 and 12 (Figure 7c,d). Aneuploidy of specific chromosomes and drug resistance has been described in cancer cells (reviewed in ) and fungi [62, 63]. To test this possibility, we generated a revertant line of L. major MTX60.4 by successive passages in the absence of MTX; under these conditions, resistance to the drug decreased (Figure 1b). Revertant cells were not as sensitive as wild-type cells to MTX but this is expected as a deletion of FT1 (Figure 6) will lead to resistant parasites . The aneuploidy of chromosomes 22 and 28 regressed to diploidy (similar to wild-type diploidy) after 30 passages, thus circumstantially linking resistance levels (Figure 1b) and copy number of these chromosomes (Figure 7a,b, lanes 2-6). With the cells now diploid, additional passages (for example, passage 42) did not decrease resistance further.
Comparative genomic hybridization
The use of DNA microarrays is now useful to understand both the mode of action of drugs and the mechanisms of drug resistance (reviewed in [44–46]). Since Leishmania has no control at the level of transcription initiation , it is unlikely that drug response profiling using microarrays will be helpful to understand the mode of action of drugs in Leishmania. Results using MTX as a lead drug and qRT-PCR to monitor key genes, such as DHFR-TS, PTR1, and FT1, appeared to confirm this lack of RNA modulation of target genes upon drug exposure (unpublished observations). This is unfortunate, as the mode of action of most anti-Leishmania drugs is unknown. Nonetheless, microarrays are likely to be useful for studying resistance in Leishmania since it is often mediated by gene amplification [3, 4] and we show here that DNA arrays hybridized to cDNAs were most valuable for detecting gene amplification events (Figures 2, 4, and 5). Since resistance is mostly correlated with gene amplification, we also used CGH and found a good qualitative correlation between RNA expression profiling and CGH (Figure 8). The technique of CGH was found to be technically simpler, but since there are clear examples of modulation in RNA level (for example, increased RNA stability) without changes in copy number of DNA in drug resistant Leishmania [64–66] (Figure 3, and Figure 7 for chromosomes 11 and 12), hybridization with cDNAs is likely to be more comprehensive. Nonetheless, modulation in RNA levels without changes in copy number of a gene is an infrequent event in drug resistant Leishmania. The use of both L. infantum and L. major MTX resistant mutants validated the design of our multi-species array but has also illustrated that the cellular resistance genotype can be complex and differ considerably between different mutants selected for resistance to the same drug. The modulation in expression of a few genes was common to both mutants, and only ABCA2 and DHFR-TS could be confirmed by qRT-PCR (Figure 3). Down-regulation of the ABC protein gene ABCA2 has never been described in MTX resistant Leishmania cells and additional investigations would be required to test whether it has any role in MTX resistance.
DHFR-TS was the first amplified gene studied in a protozoan parasite  but its exact mechanism of amplification has never been reported. In addition to detecting gene amplification events, microarray data, whether derived from RNA expression profiling or CGH, were also useful in mapping the exact regions that were amplified. We show that DHFR-TS is amplified in L. major MTX60.4 as an extrachromosomal circle through homologous recombination between non-coding repeated sequences (Figure 4). This is consistent with other loci that were also found to be amplified by homologous recombination between relatively long repeated sequences [22–24]. Blast searches have shown that these exact repeated sequences are found only on chromosome 6. Remarkably, the same similar repeated sequences (albeit with different sizes) have also been conserved in L. infantum (Figure S1 in Additional data file 2). The same observation was made for the inverted repeats close to PTR1 that were conserved between L. major and L. infantum. L. major and L. infantum are thought to have diverged 0.5 million years ago  and it thus seems that there is considerable selective pressure to keep these repeated sequences intact. Since folates and pterins are important for Leishmania growth, it is possible that the presence of these repeats may allow a strategy to rapidly increase DHFR-TS or PTR1 levels in conditions of limited substrates. With its lack of transcription initiation control, Leishmania may utilize this alternative strategy of flanking key metabolic genes by repeated sequences to amplify these genes when required. Consistent with this proposal, DNA amplification has been observed in Leishmania cells subjected to nutrient shocks .
PTR1 is a well established MTX resistance gene product [30, 31] and the amplification of its gene was first reported as part of extrachromosomal circles [6, 34–36]. Linear amplification of PTR1 with inverted duplications was described later [16, 24, 37] and linear amplicons could be precursors of circular amplicons . Linear amplicons derived from other loci than the PTR1 region with inverted duplications have also been described in Leishmania [69–73]. The microarray hybridization data have enabled the elaboration of a plausible model for the generation of a linear amplicon that contained large inverted duplications formed at the site of inverted repeats (Figure 5). This is consistent with other models of gene amplification in Leishmania [16, 37] where inverted repeats seem to be a major pathway to generate amplified large DNA palindromes (inverted duplications), as described in Tetrahymena , yeast  and mammalian cancer cells [76, 77]. One of the large inverted duplications extends from the inverted repeats, where rearrangement has occurred, to the telomeric sequences (Figure 5). These data exclude the necessity of chromosomal breaks/rearrangements at two independent positions, but it remains to be determined whether a double-stranded break, a single-stranded break or blocks in replication are facilitating inverted repeat annealing.
Gene deletions were thought to be associated with MTX resistance in Leishmania [19, 20] but had not yet been characterized at the molecular level. The microarray data, either derived from RNA expression profiling or CGH, has led to the observation that a diploid non-conservative deletion occurred by homologous recombination between two members of the large FBT gene family (Figure 6). The mechanism of gene deletion thus resembles the mechanism of amplification. Usually, amplification in Leishmania is conservative, and only a few instances of non-conservative amplification (loss of one allele) have been described in it [3, 22, 23]. In the L. major MTX60.4 mutant, we observed a diploid deletion of the FT1 gene (Figure 6). It is not known whether the second allele is deleted by homologous recombination or by a gene conversion event such as a loss of heterozygosity, but there is a strong selection pressure to delete FT1, the main folate (and MTX) transporter in Leishmania. Without FT1, cells can become resistant to MTX but folates or related molecules will still need to be transported. It will be of interest to determine whether the fusion FBT protein produced by the recombination event (Figure 6) is active or not.
The microarray approach has shown that modulation of gene expression could (rarely) be due to differential RNA expression without changes in copy number (Figure 3) ; it could be more frequently due to gene amplification (Figures 4 and 5) and, as determined now, to gene deletion (Figure 6). Two novel strategies were highlighted through the use of microarrays. In the L. major MTX60.4 mutant, the entire set of genes of chromosomes 11 and 12 is down-regulated while all the genes present on chromosomes 22 and 28 and possibly a few other chromosomes are overexpressed. The mechanism underlying an upregulation in gene expression results from a change in chromosome ploidy (Figure 7). Changes in ploidy have been observed when attempting to inactivate essential genes in Leishmania , but not in resistant parasites. We recently observed a similar phenomenon with other resistant Leishmania cells (P Leprohon et al., unpublished data), suggesting that chromosome aneuploidy is part of the Leishmania arsenal for responding to drug pressure. There was a good correlation between resistance levels and the copy number of these supernumerary chromosomes (Figures 1 and 7), linking this genetic event to the resistance phenotype. Obviously, additional studies will be required to determine which gene(s) is (are) responsible for resistance. A putative mechanism for increasing the levels of a gene product in Leishmania would thus be to generate supernumerary chromosomes. This may occur when direct or inverted repeats are absent in the vicinity of a gene conferring a selective advantage. While this is plausible, especially for an organism lacking control at the level of transcription initiation, this drug induced aneuploidy has been well documented in cells with transcriptional control, such as cancer cells (reviewed in ) or fungi [62, 63]. The mechanism of down-regulation of whole chromosome expression does not seem to involve a change in chromosome number (Figures 7 and 8) and may involve epigenetic factors that will need to be investigated.
The microarray approach was useful in highlighting several mechanisms used by resistant cells to modulate the copy number of genes by: gene deletion or extrachromosomal circular or linear amplicons; through supernumerary chromosomes; and by decreasing the expression of whole chromosomes by a mechanism that remains to be identified. In the case of the first two events, the rearrangements have occurred at the site of repeated (direct or inverted) sequences. It is possible that these repeats are not randomly distributed to allow the amplification of specific chromosomal regions. Using DNA microarrays it was shown that inverted duplications are frequent in cancer cells; these are not randomly distributed, and a subset are associated with gene amplification . The availability of DNA microarrays for Leishmania has highlighted the role of repeated sequences and of chromosome ploidy in responding to environmental changes. Aneuploidy has been suggested as an important cause of cancer specific drug resistance  and further work should reveal the potential importance of this phenomenon in drug resistance in Leishmania.
Materials and methods
The wild-type strain L. major LV39 and the mutants L. major MTX60.4 have been described previously . The L. infantum strain (MHOM/MA/67/ITMAP-263) was selected in vitro in a stepwise fashion starting with its EC50 (0.5 μM) with doubling concentrations of MTX when cells were adapted to yield L. infantum MTX20.5 growing at 20 μM of MTX. All cells were grown in M199 medium supplemented with 10% heat-inactivated fetal bovine serum and 5 μg/ml hemin at 25°C.
Chromosomes in agarose blocks were prepared and separated by PFGE as described previously . For Southern blot and PCR, genomic DNA was isolated using the DNAzol technique (Invitrogen, Carlsbad, CA, USA) as recommended by the manufacturer. Southern blots, hybridization, and washing conditions were done following standard protocols . For chromosome copy number investigation, Southern spots were quantified using ImageQuant 5.2 (GE Healthcare, Upsala, Sweden) and the reference gene α-tubulin was used for normalization.
L. infantum and L. majorDNA oligonucleotides full genome microarray design
The recent completion of the sequence of the L. major  and L. infantum  genomes, allowed the generation of multispecies high-density oligonucleotide microarrays. Our analysis of open reading frame sequence conservation between L. major and L. infantum revealed that these two species share 91-96% nucleotide identity, suggesting that interspecies microarray probes can be designed. Therefore, 70-mer oligonucleotides were designed for each open reading frame of L. infantum and L. major using automated bioinformatic procedures. The genomes of both species were first compared using BLAST and homologous genes were grouped together. Probes were designed with consistent thermodynamic properties. Probes were initially designed for L. infantum with the added requirement that the region targeted by the probes had perfect homology between both species. For common probes, up to 2 mismatches (out of 70 nucleotides) were tolerated. In the case that more than two mismatches were present in a given gene between L. infantum and L. major, a new probe was designed specifically for L. major (956 probes). The microarray included a total of 8,978 70-mer probes that recognized with no mismatches all L. infantum genes (8,184, GeneDB version 3) and also all L. major genes (8,370, GeneDB version 5.1) with a small percentage of the probes having at most 2 mismatches. Also, 372 control probes were included in the microarray for assessing synthesis variability, and location of the probe within a given open reading frame and of mismatches on hybridization. The probes were synthesized in 384-well plates by Invitrogen. The microarrays were printed on SuperChip (Erie Scientific, Portsmouth, NH, USA) using a BioRobotics MicroGrid (Genomic solutions Inc, Ann Arbor, MI, USA). Each probe was printed in duplicate. Our microarray platform is described in the Gene Expression Omnibus (GEO) with accession number GPL6855.
Total RNA preparation and labeling
Total RNA was isolated from 108 Leishmania cells during the mid-log phase using RNeasy Plus Mini Kit (QIAGEN, Hilden, Germany). The RNA preparation was treated with TURBO DNase (Ambion, Austin, TX, USA) to avoid any genomic contamination. The purity, integrity and quantity of the RNA were assessed on the Agilent 2100 bioanalyzer with the RNA 6000 Nano LabChip reagent set (Agilent Technologies, Santa Clara, CA, USA). For each probe, 10 μg of RNA were converted to aminoallyl-dUTP incorporated cDNA using random hexamers (Roche, Basel, Switzerland) and the SuperScript III RNase H Reverse Transcriptase (Invitrogen). Probes were thereafter coupled to the fluorescent dye Alexa Fluor555 or Alexa Fluor647 (Invitrogen) following the manufacturer's recommendations. Fluorescent probes were then purified with MinElute Spin Columns (QIAGEN) and quantified spectrophotometrically.
Genomic DNA preparation and labeling
Genomic DNA from 108 cells was isolated using the DNAzol technique (Invitrogen) as recommended by the manufacturer. Total DNA was then fragmented by successive passages through 22G1" and 27G 1/2" needles (Becton Dickinson Franklin Lakes, NJ, USA). Fragmented DNA was then double digested with PvuII and MseI restriction enzymes. Digested DNA was purified by phenol-chloroform, followed by an ethanol precipitation. For each probe, 4 μg of purified fragmented and digested genomic DNA were converted to fluorescently labeled DNA using Cy5- or Cy3-dCTP (Amersham, Piscataway, NJ, USA), random hexamers (Roche) and the exo- Klenow DNA polymerase (NEB, Ipswich, MA, USA). Fluorescent probes were then purified with ArrayIt columns (TeleChem International, Sunnyvale, CA, USA) and quantified spectrophotometrically.
Prehybridization and hybridization were performed at 42°C under immersion (Corning chambers, Corning, NY, USA). Slides were prehybridized for 90 minutes in PreHYB Solution (5× Denhardt, 30% formamide, 6× SSPE, 0.5% SDS, 100 μg/ml salmon sperm DNA). Then, slides were first washed 2 times at 42°C for 5 minutes in 2× SSC, 0.1% SDS with gentle agitation. Subsequent washes were at room temperature, 3 minutes each, in 1× SSC, 0.2× SSC and 0.05× SSC. Slides were then dipped in 100% isopropanol and dried by centrifugation. For hybridization, Alexa Fluor555 and 647 cDNA probes were dried and resuspended in the HYB solution (2.5× Denhardt, 30% formamide, 6× SSPE, 0.5% SDS, 100 μg/ml salmon sperm DNA, 750 μg/ml yeast tRNA), then mixed, denatured 5 minutes at 95°C and cooled slowly to 42°C. Mixed probes were applied on the array under a lifterslip. Hybridization was performed for 16 h. Washes after hybridization were the same as those described for the prehybridization.
Fluorescence detection, data processing and statistical analysis
The Perkin Elmer ScanArray 4000XL Scanner was used for image acquisition (Perkin Elmer, Waltham, MA, USA). GenePix Pro 6.0 image analysis software (Axon Instruments, Union City, CA, USA) was used to quantify the fluorescence signal intensities of the array features. Four different RNA preparations of each mutant and their respective wild-type strain were analyzed, including dye-swaps. Raw data from GenePix were imported in R 2.2.1 for normalization and statistical analyses were performed using the LIMMA (version 2.7.3) package [83–85]. Before processing, probes were flagged according to the hybridization signal quality . Weights were assigned to each array in order to give less weight to arrays of lesser quality . Data were corrected using background subtraction based on convolution of normal and exponential distributions . Intra-array normalization was carried out using the 'print-tip loess' statistical method and inter-array normalization was done by using the 'quantiles of A' method for each array . Statistical analysis was done using linear model fitting and standard errors were moderated using a simple empirical Bayes . Multiple testing corrections were done using the FDR method with a threshold p-value of 0.05. Only genes statistically significant with an absolute log ratio greater than 0.58 (log2 1.5) were considered as differentially expressed. Species comparison was performed only on probes that had less than two mismatches when hybridized to either Leishmania species. GeneSpring GX 3.1 was used for the generation of scatter plots and for chromosome by chromosome analysis. The entire data set has been deposited in GEO under the accession number series GSE9949. The comparative genomic hybridization data are deposited under reference number GSE11623.
Three independent RNA preparations were conducted for each condition. First-strand cDNA was synthesized from 2 μg of total RNA using the Superscript III RNase H Reverse Transcriptase enzyme and random hexamers (Roche) according to the manufacturer's instructions. The resulting cDNA samples were stored at -20°C until use. Control PCR amplification was carried out using primers from different internal controls (GAPDH and actin) to evaluate the uniformity of cDNA synthesis in different samples. Primers, TaqMan probes, experimental procedures and quantification for qRT-PCR of the folate transporter genes was as described (AA Ouameur et al., unpublished data) using the glyceraldehyde-3-phosphate dehydrogenase gene (GAPDH) for normalization. For all other genes, equal amounts of cDNA were run in triplicate and amplified in a 15 μl reaction containing 7.5 μl of 2× Universal PCR Master Mix (Applied Biosystems, Foster City, CA, USA), 10 nM of Z-tailed forward primer, 100 nM of reverse primer, 250 nM of Amplifluor Uniprimer probe (Chemicon Int., Temecula, CA, USA), and 1 μl of cDNA target. Reactions were performed at the Gene Quantification core laboratory of the Centre de Génomique de Québec using the Applied Biosystems Prism 7900 Sequence Detector . Amplification was normalized to two genes showing a highly stable expression in wild-type and resistant strains: LinJ18_V3.0630/LmjF18.0620 encoding a putative 60S ribosomal protein L10a, and LinJ36_V3.0850/LmjF36.2500 encoding a chromatin assembly factor 1 subunit b-like protein.
Additional data files
The following additional data are available with the online version of this paper. Additional data file 1 contains Table S1, which lists the differential expression measured by the full-genome microarray analysis. Additional data file 2 contains supplementary Figures S1-S3. Additional data file 3 contains supplementary Figure S4, which shows the results of the comparative genomic hybridization analyses of L. major MTX60.4 versus the respective wild-type cells.
comparative genomic hybridization
folate biopterin transporter
Gene Expression Omnibus
pulsed-field gel electrophoresis
quantitative real-time reverse transcription PCR.
We are grateful to Dr Eric Madore from the Centre Génomique du Centre de Recherche en Infectiologie for help during the optimization process of the microarray hybridizations. This work was funded in part by a CIHR group grant to JC, MOl, MOu, BP and MJT and operating grants to MOu. JMU is a Strategic Training Fellow of the Strategic Training Program in Microbial Resistance, a partnership of the CIHR Institute of Infection and Immunity and the Fonds de Recherche en Santé du Québec. AAO and FR are recipients of CIHR studentships. JC holds the Canada Research Chair in Medical Genomics, MJT holds the Canada Research Chair in Human Immuno-Retrovirology. BP and MOl are Burroughs Wellcome Fund New Investigator in Molecular Parasitology and the holders of FRSQ senior scholarships. MOu is a Burroughs Wellcome Fund Scholar in Molecular Parasitology and holds the Canada Research Chair in Antimicrobial Resistance.
- Handman E: Leishmania sis: current status of vaccine development. Clin Microbiol Rev. 2001, 14: 229-243. 10.1128/CMR.14.2.229-243.2001.PubMedPubMed CentralView ArticleGoogle Scholar
- Sundar S, More DK, Singh MK, Singh VP, Sharma S, Makharia A, Kumar PC, Murray HW: Failure of pentavalent antimony in visceral leishmaniasis in India: report from the center of the Indian epidemic. Clin Infect Dis. 2000, 31: 1104-1107. 10.1086/318121.PubMedView ArticleGoogle Scholar
- Beverley SM: Gene amplification in Leishmania. Annu Rev Microbiol. 1991, 45: 417-444. 10.1146/annurev.mi.45.100191.002221.PubMedView ArticleGoogle Scholar
- Borst P, Ouellette M: New mechanisms of drug resistance in parasitic protozoa. Annu Rev Microbiol. 1995, 49: 427-460. 10.1146/annurev.mi.49.100195.002235.PubMedView ArticleGoogle Scholar
- Coderre JA, Beverley SM, Schimke RT, Santi DV: Overproduction of a bifunctional thymidylate synthetase-dihydrofolate reductase and DNA amplification in methotrexate-resistant Leishmania tropica. Proc Natl Acad Sci USA. 1983, 80: 2132-2136. 10.1073/pnas.80.8.2132.PubMedPubMed CentralView ArticleGoogle Scholar
- Beverley SM, Coderre JA, Santi DV, Schimke RT: Unstable DNA amplifications in methotrexate-resistant Leishmania consist of extrachromosomal circles which relocalize during stabilization. Cell. 1984, 38: 431-439. 10.1016/0092-8674(84)90498-7.PubMedView ArticleGoogle Scholar
- Garvey EP, Santi DV: Stable amplified DNA in drug-resistant Leishmania exists as extrachromosomal circles. Science. 1986, 233: 535-540. 10.1126/science.3726545.PubMedView ArticleGoogle Scholar
- Dewes H, Ostergaard HL, Simpson L: Impaired drug uptake in methotrexate resistant Crithidia fasciculata without changes in dihydrofolate reductase activity or gene amplification. Mol Biochem Parasitol. 1986, 19: 149-161. 10.1016/0166-6851(86)90120-9.PubMedView ArticleGoogle Scholar
- Nare B, Luba J, Hardy LW, Beverley S: New approaches to Leishmania chemotherapy: pteridine reductase 1 (PTR1) as a target and modulator of antifolate sensitivity [In Process Citation]. Parasitology. 1997, 114 (Suppl): S101-110.PubMedGoogle Scholar
- Ouellette M, Drummelsmith J, El Fadili A, Kundig C, Richard D, Roy G: Pterin transport and metabolism in Leishmania and related trypanosomatid parasites. Int J Parasitol. 2002, 32: 385-398. 10.1016/S0020-7519(01)00346-0.PubMedView ArticleGoogle Scholar
- Hardy LW, Matthews W, Nare B, Beverley SM: Biochemical and genetic tests for inhibitors of Leishmania pteridine pathways. Exp Parasitol. 1997, 87: 157-169. 10.1006/expr.1997.4207.PubMedView ArticleGoogle Scholar
- Chowdhury SF, Di Lucrezia R, Guerrero RH, Brun R, Goodman J, Ruiz-Perez LM, Pacanowska DG, Gilbert IH: Novel inhibitors of Leishmania l dihydrofolate reductase. Bioorg Med Chem Lett. 2001, 11: 977-980. 10.1016/S0960-894X(01)00089-0.PubMedView ArticleGoogle Scholar
- Khabnadideh S, Pez D, Musso A, Brun R, Perez LM, Gonzalez-Pacanowska D, Gilbert IH: Design, synthesis and evaluation of 2,4-diaminoquinazolines as inhibitors of trypanosomal and leishmanial dihydrofolate reductase. Bioorg Med Chem. 2005, 13: 2637-2649. 10.1016/j.bmc.2005.01.025.PubMedView ArticleGoogle Scholar
- Ellenberger TE, Beverley SM: Reductions in methotrexate and folate influx in methotrexate-resistant lines of Leishmania major are independent of R or H region amplification. J Biol Chem. 1987, 262: 13501-13506.PubMedGoogle Scholar
- Kaur K, Coons T, Emmett K, Ullman B: Methotrexate-resistant Leishmania donovani genetically deficient in the folate-methotrexate transporter. J Biol Chem. 1988, 263: 7020-7028.PubMedGoogle Scholar
- Papadopoulou B, Roy G, Ouellette M: Frequent amplification of a short chain dehydrogenase gene as part of circular and linear amplicons in methotrexate resistant Leishmania. Nucleic Acids Res. 1993, 21: 4305-4312. 10.1093/nar/21.18.4305.PubMedPubMed CentralView ArticleGoogle Scholar
- Gamarro F, Chiquero MJ, Amador MV, Legare D, Ouellette M, Castanys S: P-glycoprotein overexpression in methotrexate-resistant Leishmania tropica. Biochem Pharmacol. 1994, 47: 1939-1947. 10.1016/0006-2952(94)90067-1.PubMedView ArticleGoogle Scholar
- Cunningham ML, Beverley SM: Pteridine salvage throughout the Leishmania infectious cycle: implications for antifolate chemotherapy. Mol Biochem Parasitol. 2001, 113: 199-213. 10.1016/S0166-6851(01)00213-4.PubMedView ArticleGoogle Scholar
- Richard D, Leprohon P, Drummelsmith J, Ouellette M: Growth phase regulation of the main folate transporter of Leishmania infantum and its role in methotrexate resistance. J Biol Chem. 2004, 279: 54494-54501. 10.1074/jbc.M409264200.PubMedView ArticleGoogle Scholar
- Richard D, Kundig C, Ouellette M: A new type of high affinity folic acid transporter in the protozoan parasite Leishmania and deletion of its gene in methotrexate-resistant cells. J Biol Chem. 2002, 277: 29460-29467. 10.1074/jbc.M204796200.PubMedView ArticleGoogle Scholar
- El Fadili A, Kundig C, Roy G, Ouellette M: Inactivation of the Leishmania tarentolae pterin transporter (BT1) and reductase (PTR1) genes leads to viable parasites with changes in folate metabolism and hypersensitivity to the antifolate methotrexate. J Biol Chem. 2004, 279: 18575-18582. 10.1074/jbc.M400652200.PubMedView ArticleGoogle Scholar
- Ouellette M, Hettema E, Wust D, Fase-Fowler F, Borst P: Direct and inverted DNA repeats associated with P-glycoprotein gene amplification in drug resistant Leishmania. EMBO J. 1991, 10: 1009-1016.PubMedPubMed CentralGoogle Scholar
- Grondin K, Papadopoulou B, Ouellette M: Homologous recombination between direct repeat sequences yields P-glycoprotein containing amplicons in arsenite resistant Leishmania. Nucleic Acids Res. 1993, 21: 1895-1901. 10.1093/nar/21.8.1895.PubMedPubMed CentralView ArticleGoogle Scholar
- Grondin K, Roy G, Ouellette M: Formation of extrachromosomal circular amplicons with direct or inverted duplications in drug-resistant Leishmania tarentolae. Mol Cell Biol. 1996, 16: 3587-3595.PubMedPubMed CentralView ArticleGoogle Scholar
- Hightower RC, Wong ML, Ruiz-Perez L, Santi DV: Electron microscopy of amplified DNA forms in antifolate-resistant Leishmania. J Biol Chem. 1987, 262: 14618-14624.PubMedGoogle Scholar
- Kapler GM, Beverley SM: Transcriptional mapping of the amplified region encoding the dihydrofolate reductase-thymidylate synthase of Leishmania major reveals a high density of transcripts, including overlapping and antisense RNAs. Mol Cell Biol. 1989, 9: 3959-3972.PubMedPubMed CentralView ArticleGoogle Scholar
- Arrebola R, Olmo A, Reche P, Garvey EP, Santi DV, Ruiz-Perez LM, Gonzalez-Pacanowska D: Isolation and characterization of a mutant dihydrofolate reductase- thymidylate synthase from methotrexate-resistant Leishmania cells. J Biol Chem. 1994, 269: 10590-10596.PubMedGoogle Scholar
- Kündig C, Leblanc E, Papadopoulou B, Ouellette M: Role of the locus and of the resistance gene on gene amplification frequency in methotrexate resistant Leishmania tarentolae. Nucleic Acids Res. 1999, 27: 3653-3659. 10.1093/nar/27.18.3653.PubMedPubMed CentralView ArticleGoogle Scholar
- Guimond C, Trudel N, Brochu C, Marquis N, El Fadili A, Peytavi R, Briand G, Richard D, Messier N, Papadopoulou B, Corbeil J, Bergeron MG, Légaré D, Ouellette M: Modulation of gene expression in Leishmania drug resistant mutants as determined by targeted DNA microarrays. Nucleic Acids Res. 2003, 31: 5886-5896. 10.1093/nar/gkg806.PubMedPubMed CentralView ArticleGoogle Scholar
- Callahan HL, Beverley SM: A member of the aldoketo reductase family confers methotrexate resistance in Leishmania. J Biol Chem. 1992, 267: 24165-24168.PubMedGoogle Scholar
- Papadopoulou B, Roy G, Ouellette M: A novel antifolate resistance gene on the amplified H circle of Leishmania. EMBO J. 1992, 11: 3601-3608.PubMedPubMed CentralGoogle Scholar
- Nare B, Hardy LW, Beverley SM: The roles of pteridine reductase 1 and dihydrofolate reductase- thymidylate synthase in pteridine metabolism in the protozoan parasite Leishmania major. J Biol Chem. 1997, 272: 13883-13891. 10.1074/jbc.272.21.13883.PubMedView ArticleGoogle Scholar
- Wang J, Leblanc E, Chang CF, Papadopoulou B, Bray T, Whiteley JM, Lin SX, Ouellette M: Pterin and folate reduction by the Leishmania tarentolae H locus short- chain dehydrogenase/reductase PTR1. Arch Biochem Biophys. 1997, 342: 197-202. 10.1006/abbi.1997.0126.PubMedView ArticleGoogle Scholar
- White TC, Fase-Fowler F, van Luenen H, Calafat J, Borst P: The H circles of Leishmania tarentolae are a unique amplifiable system of oligomeric DNAs associated with drug resistance. J Biol Chem. 1988, 263: 16977-16983.PubMedGoogle Scholar
- Petrillo-Peixoto ML, Beverley SM: Amplified DNAs in laboratory stocks of Leishmania tarentolae: extrachromosomal circles structurally and functionally similar to the inverted-H-region amplification of methotrexate-resistant Leishmania major. Mol Cell Biol. 1988, 8: 5188-5199.PubMedPubMed CentralView ArticleGoogle Scholar
- Hightower RC, Ruiz-Perez LM, Wong ML, Santi DV: Extrachromosomal elements in the lower eukaryote Leishmania. J Biol Chem. 1988, 263: 16970-16976.PubMedGoogle Scholar
- Olmo A, Arrebola R, Bernier V, Gonzalez-Pacanowska D, Ruiz-Perez LM: Co-existence of circular and multiple linear amplicons in methotrexate-resistant Leishmania. Nucleic Acids Res. 1995, 23: 2856-2864. 10.1093/nar/23.15.2856.PubMedPubMed CentralView ArticleGoogle Scholar
- Grondin K, Kundig C, Roy G, Ouellette M: Linear amplicons as precursors of amplified circles in methotrexate-resistant Leishmania tarentolae. Nucleic Acids Res. 1998, 26: 3372-3378. 10.1093/nar/26.14.3372.PubMedPubMed CentralView ArticleGoogle Scholar
- El Fadili A, Richard D, Kundig C, Ouellette M: Effect of polyglutamylation of methotrexate on its accumulation and the development of resistance in the protozoan parasite Leishmania. Biochem Pharmacol. 2003, 66: 999-1008. 10.1016/S0006-2952(03)00417-9.PubMedView ArticleGoogle Scholar
- Gagnon D, Foucher A, Girard I, Ouellette M: Stage specific gene expression and cellular localization of two isoforms of the serine hydroxymethyltransferase in the protozoan parasite Leishmania. Mol Biochem Parasitol. 2006, 150: 63-71. 10.1016/j.molbiopara.2006.06.009.PubMedView ArticleGoogle Scholar
- Drummelsmith J, Brochu V, Girard I, Messier N, Ouellette M: Proteome mapping of the protozoan parasite Leishmania and application to the study of drug targets and resistance mechanisms. Mol Cell Proteomics. 2003, 2: 146-155. 10.1074/mcp.M200085-MCP200.PubMedView ArticleGoogle Scholar
- Gallego C, Estevez AM, Farez E, Ruiz-Perez LM, Gonzalez-Pacanowska D: Overexpression of AP endonuclease protects Leishmania major cells against methotrexate induced DNA fragmentation and hydrogen peroxide. Mol Biochem Parasitol. 2005, 141: 191-197. 10.1016/j.molbiopara.2005.03.002.PubMedView ArticleGoogle Scholar
- Wilson M, DeRisi J, Kristensen HH, Imboden P, Rane S, Brown PO, Schoolnik GK: Exploring drug-induced alterations in gene expression in Mycobacterium tuberculosis by microarray hybridization. Proc Natl Acad Sci USA. 1999, 96: 12833-12838. 10.1073/pnas.96.22.12833.PubMedPubMed CentralView ArticleGoogle Scholar
- Lee CH, Macgregor PF: Using microarrays to predict resistance to chemotherapy in cancer patients. Pharmacogenomics. 2004, 5: 611-625. 10.1517/14622418.104.22.1681.PubMedView ArticleGoogle Scholar
- Brazas MD, Hancock RE: Using microarray gene signatures to elucidate mechanisms of antibiotic action and resistance. Drug Discov Today. 2005, 10: 1245-1252. 10.1016/S1359-6446(05)03566-X.PubMedView ArticleGoogle Scholar
- Ouellette M, Drummelsmith J, Leprohon P, El Fadili K, Foucher A, Vergnes B, Légaré D: Drug Resistance in Leishmania. 2007, Norwich, UK: Horizon PressGoogle Scholar
- Saxena A, Worthey EA, Yan S, Leland A, Stuart KD, Myler PJ: Evaluation of differential gene expression in Leishmania major Friedlin procyclics and metacyclics using DNA microarray analysis. Mol Biochem Parasitol. 2003, 129: 103-114. 10.1016/S0166-6851(03)00100-2.PubMedView ArticleGoogle Scholar
- Akopyants NS, Matlib RS, Bukanova EN, Smeds MR, Brownstein BH, Stormo GD, Beverley SM: Expression profiling using random genomic DNA microarrays identifies differentially expressed genes associated with three major developmental stages of the protozoan parasite Leishmania major. Mol Biochem Parasitol. 2004, 136: 71-86. 10.1016/j.molbiopara.2004.03.002.PubMedView ArticleGoogle Scholar
- Saxena A, Lahav T, Holland N, Aggarwal G, Anupama A, Huang Y, Volpin H, Myler PJ, Zilberstein D: Analysis of the Leishmania donovani transcriptome reveals an ordered progression of transient and permanent changes in gene expression during differentiation. Mol Biochem Parasitol. 2007, 152: 53-65. 10.1016/j.molbiopara.2006.11.011.PubMedPubMed CentralView ArticleGoogle Scholar
- Srividya G, Duncan R, Sharma P, Raju BV, Nakhasi HL, Salotra P: Transcriptome analysis during the process of in vitro differentiation of Leishmania donovani using genomic microarrays. Parasitology. 2007, 134: 1527-1539. 10.1017/S003118200700296X.PubMedView ArticleGoogle Scholar
- Almeida R, Gilmartin BJ, McCann SH, Norrish A, Ivens AC, Lawson D, Levick MP, Smith DF, Dyall SD, Vetrie D, Freeman TC, Coulson RM, Sampaio I, Schneider H, Blackwell JM: Expression profiling of the Leishmania life cycle: cDNA arrays identify developmentally regulated genes present but not annotated in the genome. Mol Biochem Parasitol. 2004, 136: 87-100. 10.1016/j.molbiopara.2004.03.004.PubMedView ArticleGoogle Scholar
- Singh N, Almeida R, Kothari H, Kumar P, Mandal G, Chatterjee M, Venkatachalam S, Govind MK, Mandal SK, Sundar S: Differential gene expression analysis in antimony-unresponsive Indian kala azar (visceral leishmaniasis) clinical isolates by DNA microarray. Parasitology. 2007, 134: 777-787. 10.1017/S0031182007002284.PubMedView ArticleGoogle Scholar
- McNicoll F, Drummelsmith J, Muller M, Madore E, Boilard N, Ouellette M, Papadopoulou B: A combined proteomic and transcriptomic approach to the study of stage differentiation in Leishmania infantum. Proteomics. 2006, 6: 3567-3581. 10.1002/pmic.200500853.PubMedView ArticleGoogle Scholar
- Leprohon P, Legare D, Girard I, Papadopoulou B, Ouellette M: Modulation of Leishmania ABC protein gene expression through life stages and among drug-resistant parasites. Eukaryot Cell. 2006, 5: 1713-1725. 10.1128/EC.00152-06.PubMedPubMed CentralView ArticleGoogle Scholar
- Holzer TR, McMaster WR, Forney JD: Expression profiling by whole-genome interspecies microarray hybridization reveals differential gene expression in procyclic promastigotes, lesion-derived amastigotes, and axenic amastigotes in Leishmania mexicana. Mol Biochem Parasitol. 2006, 146: 198-218. 10.1016/j.molbiopara.2005.12.009.PubMedView ArticleGoogle Scholar
- Leifso K, Cohen-Freue G, Dogra N, Murray A, McMaster WR: Genomic and proteomic expression analysis of Leishmania promastigote and amastigote life stages: the Leishmania genome is constitutively expressed. Mol Biochem Parasitol. 2007, 152: 35-46. 10.1016/j.molbiopara.2006.11.009.PubMedView ArticleGoogle Scholar
- El Fadili K, Messier N, Leprohon P, Roy G, Guimond C, Trudel N, Saravia NG, Papadopoulou B, Legare D, Ouellette M: Role of the ABC transporter MRPA (PGPA) in antimony resistance in Leishmania infantum axenic and intracellular amastigotes. Antimicrob Agents Chemother. 2005, 49: 1988-1993. 10.1128/AAC.49.5.1988-1993.2005.PubMedPubMed CentralView ArticleGoogle Scholar
- Clayton CE: Life without transcriptional control? From fly to man and back again. EMBO J. 2002, 21: 1881-1888. 10.1093/emboj/21.8.1881.PubMedPubMed CentralView ArticleGoogle Scholar
- Blaineau C, Bastien P, Rioux JA, Roizes G, Pages M: Long-range restriction maps of size-variable homologous chromosomes in Leishmania infantum. Mol Biochem Parasitol. 1991, 46: 292-302. 10.1016/0166-6851(91)90053-9.PubMedView ArticleGoogle Scholar
- Sanger Institute Pathogen Sequencing Unit. [http://www.genedb.org/]
- Duesberg P, Li R, Sachs R, Fabarius A, Upender MB, Hehlmann R: Cancer drug resistance: the central role of the karyotype. Drug Resist Updat. 2007, 10: 51-58. 10.1016/j.drup.2007.02.003.PubMedView ArticleGoogle Scholar
- Selmecki A, Forche A, Berman J: Aneuploidy and isochromosome formation in drug-resistant Candida albicans. Science. 2006, 313: 367-370. 10.1126/science.1128242.PubMedPubMed CentralView ArticleGoogle Scholar
- Coste A, Selmecki A, Forche A, Diogo D, Bougnoux ME, d'Enfert C, Berman J, Sanglard D: Genotypic evolution of azole resistance mechanisms in sequential Candida albicans isolates. Eukaryot Cell. 2007, 6: 1889-1904. 10.1128/EC.00151-07.PubMedPubMed CentralView ArticleGoogle Scholar
- Haimeur A, Guimond C, Pilote S, Mukhopadhyay R, Rosen BP, Poulin R, Ouellette M: Elevated levels of polyamines and trypanothione resulting from overexpression of the ornithine decarboxylase gene in arsenite-resistant Leishmania. Mol Microbiol. 1999, 34: 726-735. 10.1046/j.1365-2958.1999.01634.x.PubMedView ArticleGoogle Scholar
- Drummelsmith J, Girard I, Trudel N, Ouellette M: Differential protein expression analysis of Leishmania major reveals novel roles for methionine adenosyltransferase and S-adenosylmethionine in methotrexate resistance. J Biol Chem. 2004, 279: 33273-33280. 10.1074/jbc.M405183200.PubMedView ArticleGoogle Scholar
- Marquis N, Gourbal B, Rosen BP, Mukhopadhyay R, Ouellette M: Modulation in aquaglyceroporin AQP1 gene transcript levels in drug-resistant Leishmania. Mol Microbiol. 2005, 57: 1690-1699. 10.1111/j.1365-2958.2005.04782.x.PubMedView ArticleGoogle Scholar
- Ibrahim ME, Barker DC: The origin and evolution of the Leishmania donovani complex as inferred from a mitochondrial cytochrome oxidase II gene sequence. Infect Genet Evol. 2001, 1: 61-68. 10.1016/S1567-1348(01)00009-0.PubMedView ArticleGoogle Scholar
- Rovai L, Tripp C, Stuart K, Simpson L: Recurrent polymorphisms in small chromosomes of Leishmania tarentolae after nutrient stress or subcloning. Mol Biochem Parasitol. 1992, 50: 115-125. 10.1016/0166-6851(92)90249-J.PubMedView ArticleGoogle Scholar
- Tripp CA, Myler PJ, Stuart K: A DNA sequence (LD1) which occurs in several genomic organizations in Leishmania. Mol Biochem Parasitol. 1991, 47: 151-156. 10.1016/0166-6851(91)90174-5.PubMedView ArticleGoogle Scholar
- Navarro M, Liu J, Muthui D, Ortiz G, Segovia M, Hamers R: Inverted repeat structure and homologous sequences in the LD1 amplicons of Leishmania spp. Mol Biochem Parasitol. 1994, 68: 69-80. 10.1016/0166-6851(94)00147-2.PubMedView ArticleGoogle Scholar
- Ortiz G, Segovia M: Characterisation of the novel junctions of two minichromosomes of Leishmania major. Mol Biochem Parasitol. 1996, 82: 137-144. 10.1016/0166-6851(96)02724-7.PubMedView ArticleGoogle Scholar
- Dubessay P, Ravel C, Bastien P, Lignon MF, Ullman B, Pages M, Blaineau C: Effect of large targeted deletions on the mitotic stability of an extra chromosome mediating drug resistance in Leishmania. Nucleic Acids Res. 2001, 29: 3231-3240. 10.1093/nar/29.15.3231.PubMedPubMed CentralView ArticleGoogle Scholar
- Genest PA, ter Riet B, Dumas C, Papadopoulou B, van Luenen HG, Borst P: Formation of linear inverted repeat amplicons following targeting of an essential gene in Leishmania. Nucleic Acids Res. 2005, 33: 1699-1709. 10.1093/nar/gki304.PubMedPubMed CentralView ArticleGoogle Scholar
- Butler DK, Yasuda LE, Yao MC: An intramolecular recombination mechanism for the formation of the rRNA gene palindrome of Tetrahymena thermophila. Mol Cell Biol. 1995, 15: 7117-7126.PubMedPubMed CentralView ArticleGoogle Scholar
- Albrecht EB, Hunyady AB, Stark GR, Patterson TE: Mechanisms of sod2 gene amplification in Schizosaccharomyces pombe. Mol Biol Cell. 2000, 11: 873-886.PubMedPubMed CentralView ArticleGoogle Scholar
- Tanaka H, Tapscott SJ, Trask BJ, Yao MC: Short inverted repeats initiate gene amplification through the formation of a large DNA palindrome in mammalian cells. Proc Natl Acad Sci USA. 2002, 99: 8772-8777.PubMedPubMed CentralView ArticleGoogle Scholar
- Tanaka H, Cao Y, Bergstrom DA, Kooperberg C, Tapscott SJ, Yao MC: Intrastrand annealing leads to the formation of a large DNA palindrome and determines the boundaries of genomic amplification in human cancer. Mol Cell Biol. 2007, 27: 1993-2002. 10.1128/MCB.01313-06.PubMedPubMed CentralView ArticleGoogle Scholar
- Cruz AK, Titus R, Beverley SM: Plasticity in chromosome number and testing of essential genes in Leishmania by targeting. Proc Natl Acad Sci USA. 1993, 90: 1599-1603. 10.1073/pnas.90.4.1599.PubMedPubMed CentralView ArticleGoogle Scholar
- Tanaka H, Bergstrom DA, Yao MC, Tapscott SJ: Widespread and nonrandom distribution of DNA palindromes in cancer cells provides a structural platform for subsequent gene amplification. Nat Genet. 2005, 37: 320-327. 10.1038/ng1515.PubMedView ArticleGoogle Scholar
- Sambrook J, Fritsch EF, Maniatis T: Molecular Cloning. 1989, New York: Cold Spring Harbour Laboratory PressGoogle Scholar
- Ivens AC, Peacock CS, Worthey EA, Murphy L, Aggarwal G, Berriman M, Sisk E, Rajandream MA, Adlem E, Aert R, Anupama A, Apostolou Z, Attipoe P, Bason N, Bauser C, Beck A, Beverley SM, Bianchettin G, Borzym K, Bothe G, Bruschi CV, Collins M, Cadag E, Ciarloni L, Clayton C, Coulson RM, Cronin A, Cruz AK, Davies RM, De Gaudenzi J, et al: The genome of the kinetoplastid parasite, Leishmania major. Science. 2005, 309: 436-442. 10.1126/science.1112680.PubMedPubMed CentralView ArticleGoogle Scholar
- Peacock CS, Seeger K, Harris D, Murphy L, Ruiz JC, Quail MA, Peters N, Adlem E, Tivey A, Aslett M, Kerhornou A, Ivens A, Fraser A, Rajandream MA, Carver T, Norbertczak H, Chillingworth T, Hance Z, Jagels K, Moule S, Ormond D, Rutter S, Squares R, Whitehead S, Rabbinowitsch E, Arrowsmith C, White B, Thurston S, Bringaud F, Baldauf SL, et al: Comparative genomic analysis of three Leishmania species that cause diverse human disease. Nat Genet. 2007, 39: 839-847. 10.1038/ng2053.PubMedPubMed CentralView ArticleGoogle Scholar
- Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article3-PubMedGoogle Scholar
- Smyth GK, Speed T: Normalization of cDNA microarray data. Methods. 2003, 31: 265-273. 10.1016/S1046-2023(03)00155-5.PubMedView ArticleGoogle Scholar
- Smyth GK, Michaud J, Scott HS: Use of within-array replicate spots for assessing differential expression in microarray experiments. Bioinformatics. 2005, 21: 2067-2075. 10.1093/bioinformatics/bti270.PubMedView ArticleGoogle Scholar
- Oshlack A, Chabot AE, Smyth GK, Gilad Y: Using DNA microarrays to study gene expression in closely related species. Bioinformatics. 2007, 23: 1235-1242. 10.1093/bioinformatics/btm111.PubMedView ArticleGoogle Scholar
- Ritchie ME, Diyagama D, Neilson J, van Laar R, Dobrovic A, Holloway A, Smyth GK: Empirical array quality weights in the analysis of microarray data. BMC Bioinformatics. 2006, 7: 261-10.1186/1471-2105-7-261.PubMedPubMed CentralView ArticleGoogle Scholar
- Ritchie ME, Silver J, Oshlack A, Holmes M, Diyagama D, Holloway A, Smyth GK: A comparison of background correction methods for two-colour microarrays. Bioinformatics. 2007, 23: 2700-2707. 10.1093/bioinformatics/btm412.PubMedView ArticleGoogle Scholar
- Yang Y, Hoh J, Broger C, Neeb M, Edington J, Lindpaintner K, Ott J: Statistical methods for analyzing microarray feature data with replications. J Comput Biol. 2003, 10: 157-169. 10.1089/106652703321825946.PubMedView ArticleGoogle Scholar
- Quebec Genomics Center. [https://genome.ulaval.ca/qrtpcr]
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.