Whole genome profiling of spontaneous and chemically induced mutations in Toxoplasma gondii

BMC Genomics

Marc-Jan Gubbels³  Gabor T Marth³  Ira J Blader¹  Kevin M Brown²  Brian Benenati³  Bradley I Coleman³  Andrew Farrell³ 

[1] Present address: Department of Microbiology and Immunology, State University of New York, 138 Farber Hall, 3435 Main Street, Buffalo, NY 14214, USA;Present address: Department of Molecular Microbiology, Washington University School of Medicine, MPRB 9230, 4940 Parkview Place, St. Louis, MO 63110, USA;Department of Biology, Boston College, Higgins Hall 355, 140 Commonwealth Avenue, Chestnut Hill, MA 02467, USA

关键词: Apicomplexa;    SNV calling;    In vitro adaptation;    Chemical mutagenesis;    Whole genome sequencing;   

Others  :  1217244
DOI  :  10.1186/1471-2164-15-354

 received in 2013-12-23, accepted in 2014-05-02,  发布年份 2014

【 摘 要 】

Background

Next generation sequencing is helping to overcome limitations in organisms less accessible to classical or reverse genetic methods by facilitating whole genome mutational analysis studies. One traditionally intractable group, the Apicomplexa, contains several important pathogenic protozoan parasites, including the Plasmodium species that cause malaria.

Here we apply whole genome analysis methods to the relatively accessible model apicomplexan, Toxoplasma gondii, to optimize forward genetic methods for chemical mutagenesis using N-ethyl-N-nitrosourea (ENU) and ethylmethane sulfonate (EMS) at varying dosages.

Results

By comparing three different lab-strains we show that spontaneously generated mutations reflect genome composition, without nucleotide bias. However, the single nucleotide variations (SNVs) are not distributed randomly over the genome; most of these mutations reside either in non-coding sequence or are silent with respect to protein coding. This is in contrast to the random genomic distribution of mutations induced by chemical mutagenesis. Additionally, we report a genome wide transition vs transversion ratio (ti/tv) of 0.91 for spontaneous mutations in Toxoplasma, with a slightly higher rate of 1.20 and 1.06 for variants induced by ENU and EMS respectively. We also show that in the Toxoplasma system, surprisingly, both ENU and EMS have a proclivity for inducing mutations at A/T base pairs (78.6% and 69.6%, respectively).

Conclusions

The number of SNVs between related laboratory strains is relatively low and managed by purifying selection away from changes to amino acid sequence. From an experimental mutagenesis point of view, both ENU (24.7%) and EMS (29.1%) are more likely to generate variation within exons than would naturally accumulate over time in culture (19.1%), demonstrating the utility of these approaches for yielding proportionally greater changes to the amino acid sequence. These results will not only direct the methods of future chemical mutagenesis in Toxoplasma, but also aid in designing forward genetic approaches in less accessible pathogenic protozoa as well.

【 授权许可】

   
2014 Farrell et al.; licensee BioMed Central Ltd.

【 图 表 】

Figure 1.

Figure 2.

Figure 3.

Figure 4.

Figure 5.

Figure 6.

Figure 7.

【 参考文献 】

• [1]White NJ, Pukrittayakamee S, Hien TT, Faiz MA, Mokuolu OA, Dondorp AM: Malaria. Lancet 2013, 383:723-735.
• [2]Montoya JG, Liesenfeld O: Toxoplasmosis. Lancet 2004, 363:1965-1976.
• [3]Kim K, Weiss LM: Toxoplasma gondii: the model apicomplexan. Int J Parasitol 2004, 34:423-432.
• [4]Pfefferkorn ER, Pfefferkorn LC: Arabinosyl nucleosides inhibit Toxoplasma gondii and allow the selection of resistant mutants. J Parasitol 1976, 62:993-999.
• [5]Pfefferkorn ER, Pfefferkorn LC: Toxoplasma gondii: isolation and preliminary characterization of temperature-sensitive mutants. Exp Parasitol 1976, 39:365-376.
• [6]Pfefferkorn ER, Schwartzman JD, Kasper LH: Toxoplasma gondii: use of mutants to study the host-parasite relationship. Ciba Found Symp 1983, 99:74-91.
• [7]Farrell A, Thirugnanam S, Lorestani A, Dvorin JD, Eidell KP, Ferguson DJ, Anderson-White BR, Duraisingh MT, Marth GT, Gubbels MJ: A DOC2 protein identified by mutational profiling is essential for apicomplexan parasite exocytosis. Science 2012, 335:218-221.
• [8]Gubbels MJ, Lehmann M, Muthalagi M, Jerome ME, Brooks CF, Szatanek T, Flynn J, Parrot B, Radke J, Striepen B, White MW: Forward genetic analysis of the apicomplexan cell division cycle in Toxoplasma gondii. PLoS Pathog 2008, 4:e36.
• [9]White MW, Jerome ME, Vaishnava S, Guerini M, Behnke M, Striepen B: Genetic rescue of a Toxoplasma gondii conditional cell cycle mutant. Mol Microbiol 2005, 55:1060-1071.
• [10]Pfefferkorn ER, Pfefferkorn LC: Toxoplasma gondii: characterization of a mutant resistant to 5-fluorodeoxyuridine. Exp Parasitol 1977, 42:44-55.
• [11]Black MW, Arrizabalaga G, Boothroyd JC: Ionophore-resistant mutants of Toxoplasma gondii reveal host cell permeabilization as an early event in egress. Mol Cell Biol 2000, 20:9399-9408.
• [12]Arrizabalaga G, Ruiz F, Moreno S, Boothroyd JC: Ionophore-resistant mutant of Toxoplasma gondii reveals involvement of a sodium/hydrogen exchanger in calcium regulation. J Cell Biol 2004, 165:653-662.
• [13]Garrison E, Treeck M, Ehret E, Butz H, Garbuz T, Oswald BP, Settles M, Boothroyd J, Arrizabalaga G: A forward genetic screen reveals that calcium-dependent protein kinase 3 regulates egress in Toxoplasma. PLoS Pathog 2012, 8:e1003049.
• [14]Uyetake L, Ortega-Barria E, Boothroyd JC: Isolation and characterization of a cold-sensitive attachment/invasion mutant of Toxoplasma gondii. Exp Parasitol 2001, 97:55-59.
• [15]Pfefferkorn ER, Pfefferkorn LC: Quantitative studies of the mutagenesis of Toxoplasma gondii. J Parasitol 1979, 65:364-370.
• [16]Sikora P, Chawade A, Larsson M, Olsson J, Olsson O: Mutagenesis as a Tool in Plant Genetics, Functional Genomics, and Breeding. Int J Plant Genomics 2011, 13.
• [17]Anderson P: Mutagenesis. Methods Cell Biol 1995, 48:31-58.
• [18]Flibotte S, Edgley ML, Chaudhry I, Taylor J, Neil SE, Rogula A, Zapf R, Hirst M, Butterfield Y, Jones SJ, Marra MA, Barstead RJ, Moerman DG: Whole-genome profiling of mutagenesis in Caenorhabditis elegans. Genetics 2010, 185:431-441.
• [19]Acevedo-Arozena A, Wells S, Potter P, Kelly M, Cox RD, Brown SD: ENU mutagenesis, a way forward to understand gene function. Annu Rev Genomics Hum Genet 2008, 9:49-69.
• [20]Amsterdam A, Hopkins N: Mutagenesis strategies in zebrafish for identifying genes involved in development and disease. Trends Genet 2006, 22:473-478.
• [21]Barbaric I, Wells S, Russ A, Dear TN: Spectrum of ENU-induced mutations in phenotype-driven and gene-driven screens in the mouse. Environ Mol Mutagen 2007, 48:124-142.
• [22]Greene EA, Codomo CA, Taylor NE, Henikoff JG, Till BJ, Reynolds SH, Enns LC, Burtner C, Johnson JE, Odden AR, Comai L, Henikoff S: Spectrum of chemically induced mutations from a large-scale reverse-genetic screen in Arabidopsis. Genetics 2003, 164:731-740.
• [23]Hanash SM, Boehnke M, Chu EH, Neel JV, Kuick RD: Nonrandom distribution of structural mutants in ethylnitrosourea-treated cultured human lymphoblastoid cells. Proc Natl Acad Sci U S A 1988, 85:165-169.
• [24]Hillier LW, Marth GT, Quinlan AR, Dooling D, Fewell G, Barnett D, Fox P, Glasscock JI, Hickenbotham M, Huang W, Magrini VJ, Richt RJ, Sander SN, Stewart DA, Stromberg M, Tsung EF, Wylie T, Schedl T, Wilson RK, Mardis ER: Whole-genome sequencing and variant discovery in C. elegans. Nat Methods 2008, 5:183-188.
• [25]Smith DR, Quinlan AR, Peckham HE, Makowsky K, Tao W, Woolf B, Shen L, Donahue WF, Tusneem N, Stromberg MP, Stewart DA, Zhang L, Ranade SS, Warner JB, Lee CC, Coleman BE, Zhang Z, McLaughlin SF, Malek JA, Sorenson JM, Blanchard AP, Chapman J, Hillman D, Chen F, Rokhsar DS, McKernan KJ, Jeffries TW, Marth GT, Richardson PM: Rapid whole-genome mutational profiling using next-generation sequencing technologies. Genome Res 2008, 18:1638-1642.
• [26]Sarin S, Prabhu S, O'Meara MM, Pe'er I, Hobert O: Caenorhabditis elegans mutant allele identification by whole-genome sequencing. Nat Methods 2008, 5:865-867.
• [27]Shen Y, Sarin S, Liu Y, Hobert O, Pe'er I: Comparing platforms for C. elegans mutant identification using high-throughput whole-genome sequencing. PLoS One 2008, 3:e4012.
• [28]Sarin S, Bertrand V, Bigelow H, Boyanov A, Doitsidou M, Poole RJ, Narula S, Hobert O: Analysis of multiple ethyl methanesulfonate-mutagenized Caenorhabditis elegans strains by whole-genome sequencing. Genetics 2010, 185:417-430.
• [29]Zuryn S, Le Gras S, Jamet K, Jarriault S: A strategy for direct mapping and identification of mutations by whole-genome sequencing. Genetics 2010, 186:427-430.
• [30]Voz ML, Coppieters W, Manfroid I, Baudhuin A, Von Berg V, Charlier C, Meyer D, Driever W, Martial JA, Peers B: Fast homozygosity mapping and identification of a zebrafish ENU-induced mutation by whole-genome sequencing. PLoS One 2012, 7:e34671.
• [31]Shiwa Y, Fukushima-Tanaka S, Kasahara K, Horiuchi T, Yoshikawa H: Whole-Genome Profiling of a Novel Mutagenesis Technique Using Proofreading-Deficient DNA Polymerase delta. Int J Evol Biol 2012, 2012:860797.
• [32]Hobert O: The impact of whole genome sequencing on model system genetics: get ready for the ride. Genetics 2010, 184:317-319.
• [33]Brown KM, Suvorova ES, Farrell A, Wiley GB, Gubbels MJ, Marth GT, Gaffney PM, White M, Blader IJ: Chemical genetic screening identifies a small molecule that blocks Toxoplasma growth by inhibiting both host- and parasite-encoded kinases. PLoS Path 2014. accepted for publication
• [34]Sugi T, Kobayashi K, Takamea H, Gong H, Ishiwa A: Identification of mutations in TgMAPK1 of Toxoplasma gondii conferring the resistance to 1NM-PP1. Int J Parasitol Drugs Drug Resist 2013, 3:93-101.
• [35]Coleman BI, Gubbels MJ: A genetic screen to isolate Toxoplasma gondii host cell egress mutants. J Vis Exp 2012, 60:e3807.
• [36]Roos DS, Donald RG, Morrissette NS, Moulton AL: Molecular tools for genetic dissection of the protozoan parasite Toxoplasma gondii. Methods Cell Biol 1994, 45:27-63.
• [37]Sabin AB: Toxoplasmic encephalitis in children. J Am Med Assoc 1941, 116:801-807.
• [38]Donald RG, Carter D, Ullman B, Roos DS: Insertional tagging, cloning, and expression of the Toxoplasma gondii hypoxanthine-xanthine-guanine phosphoribosyltransferase gene. Use as a selectable marker for stable transformation. J Biol Chem 1996, 271:14010-14019.
• [39]Dobrowolski JM, Sibley LD: Toxoplasma invasion of mammalian cells is powered by the actin cytoskeleton of the parasite. Cell 1996, 84:933-939.
• [40]Teo CF, Zhou XW, Bogyo M, Carruthers VB: Cysteine protease inhibitors block Toxoplasma gondii microneme secretion and cell invasion. Antimicrob Agents Chemother 2007, 51:679-688.
• [41]Gubbels MJ, Li C, Striepen B: High-throughput growth assay for Toxoplasma gondii using yellow fluorescent protein. Antimicrob Agents Chemother 2003, 47:309-316.
• [42]Eidell KP, Burke T, Gubbels MJ: Development of a screen to dissect Toxoplasma gondii egress. Mol Biochem Parasitol 2010, 171:97-103.
• [43]Carey KL, Westwood NJ, Mitchison TJ, Ward GE: A small-molecule approach to studying invasive mechanisms of Toxoplasma gondii. Proc Natl Acad Sci U S A 2004, 101:7433-7438.
• [44]Lee W-P, Stromberg M, Ward A, Stewart C, Garrison E, Marth GT: MOSAIK: A hash-based algorithm for accurate next-generation sequencing short-read mapping. PLoS One 2014, 9:e90581.
• [45]Marth GT, Korf I, Yandell MD, Yeh RT, Gu Z, Zakeri H, Stitziel NO, Hillier L, Kwok PY, Gish WR: A general approach to single-nucleotide polymorphism discovery. Nat Genet 1999, 23:452-456.
• [46]Garrison E, Marth GT: Haplotype-Based Variant Detection from Short-Read Sequencing. 2012. http://arxiv.org/abs/1207.3907v2 webcite
• [47]Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA: Circos: an information aesthetic for comparative genomics. Genome Res 2009, 19:1639-1645.
• [48]Gajria B, Bahl A, Brestelli J, Dommer J, Fischer S, Gao X, Heiges M, Iodice J, Kissinger JC, Mackey AJ, Pinney DF, Roos DS, Stoeckert CJ Jr, Wang H, Brunk BP: ToxoDB: an integrated Toxoplasma gondii database resource. Nucleic Acids Res 2008, 36:D553-D556.
• [49]Sibley LD, Boothroyd JC: Virulent strains of Toxoplasma gondii comprise a single clonal lineage. Nature 1992, 359:82-85.
• [50]Martinelli A, Henriques G, Cravo P, Hunt P: Whole genome re-sequencing identifies a mutation in an ABC transporter (mdr2) in a Plasmodium chabaudi clone with altered susceptibility to antifolate drugs. Int J Parasitol 2010, 41:165-171.
• [51]Araya CL, Payen C, Dunham MJ, Fields S: Whole-genome sequencing of a laboratory-evolved yeast strain. BMC Genomics 2010, 11:88. BioMed Central Full Text
• [52]Le Crom S, Schackwitz W, Pennacchio L, Magnuson JK, Culley DE, Collett JR, Martin J, Druzhinina IS, Mathis H, Monot F, Seiboth B, Cherry B, Rey M, Berka R, Kubicek CP, Baker SE, Margeot A: Tracking the roots of cellulase hyperproduction by the fungus Trichoderma reesei using massively parallel DNA sequencing. Proc Natl Acad Sci U S A 2009, 106:16151-16156.
• [53]Yang Z, Yoder AD: Estimation of the transition/transversion rate bias and species sampling. J Mol Evol 1999, 48:274-283.
• [54]Zhang J: Rates of conservative and radical nonsynonymous nucleotide substitutions in mammalian nuclear genes. J Mol Evol 2000, 50:56-68.
• [55]Khan A, Taylor S, Su C, Sibley LD, Paulsen I, Ajioka JW: Genetics and Genome Organization of Toxoplasma Gondii. Norwich, UK: Horizon Scientific Press; 2007.
• [56]Pfefferkorn ER: Toxoplasma gondii: the enzymic defect of a mutant resistant to 5-fluorodeoxyuridine. Exp Parasitol 1978, 44:26-35.
• [57]Schwartzman JD, Pfefferkorn ER: Pyrimidine synthesis by intracellular Toxoplasma gondii. J Parasitol 1981, 67:150-158.
• [58]Donald RG, Roos DS: Insertional mutagenesis and marker rescue in a protozoan parasite: cloning of the uracil phosphoribosyltransferase locus from Toxoplasma gondii. Proc Natl Acad Sci U S A 1995, 92:5749-5753.
• [59]Goecks J, Eberhard C, Too T, Nekrutenko A, Taylor J: Web-based visual analysis for high-throughput genomics. BMC Genomics 2013, 14:397. BioMed Central Full Text
• [60]Sun YV, Levin AM, Boerwinkle E, Robertson H, Kardia SL: A scan statistic for identifying chromosomal patterns of SNP association. Genet Epidemiol 2006, 30:627-635.
• [61]Yang N, Farrell A, Niedelman W, Melo M, Lu D, Julien L, Marth GT, Gubbels MJ, Saeij JP: Genetic basis for phenotypic differences between different Toxoplasma gondii type I strains. BMC Genomics 2013, 14:467. BioMed Central Full Text
• [62]Martincorena I, Seshasayee AS, Luscombe NM: Evidence of non-random mutation rates suggests an evolutionary risk management strategy. Nature 2012, 485:95-98.
• [63]Fox BA, Ristuccia JG, Gigley JP, Bzik DJ: Efficient gene replacements in Toxoplasma gondii strains deficient for nonhomologous end-joining. Eukaryot Cell 2009, 8:520-529.
• [64]Nordstrom KJ, Albani MC, James GV, Gutjahr C, Hartwig B, Turck F, Paszkowski U, Coupland G, Schneeberger K: Mutation identification by direct comparison of whole-genome sequencing data from mutant and wild-type individuals using k-mers. Nat Biotechnol 2013, 31:325-330.
• [65]Kohler S, Delwiche CF, Denny PW, Tilney LG, Webster P, Wilson RJ, Palmer JD, Roos DS: A plastid of probable green algal origin in Apicomplexan parasites. Science 1997, 275:1485-1489.
• [66]Reiff SB, Vaishnava S, Striepen B: The HU protein is important for apicoplast genome maintenance and inheritance in Toxoplasma gondii. Eukaryot Cell 2012, 11:905-915.
• [67]Matsuzaki M, Kikuchi T, Kita K, Kojima S, Kuroiwa T: Large amounts of apicoplast nucleoid DNA and its segregation in Toxoplasma gondii. Protoplasma 2001, 218:180-191.
• [68]Pastink A, Heemskerk E, Nivard MJ, Van Vliet CJ, Vogel EW: Mutational specificity of ethyl methanesulfonate in excision-repair-proficient and -deficient strains of Drosophila melanogaster. Mol Gen Genet 1991, 229:213-218.
• [69]Klungland A, Laake K, Hoff E, Seeberg E: Spectrum of mutations induced by methyl and ethyl methanesulfonate at the hprt locus of normal and tag expressing Chinese hamster fibroblasts. Carcinogenesis 1995, 16:1281-1285.
• [70]Onyango DO, Naguleswaran A, Delaplane S, Reed A, Kelley MR, Georgiadis MM, Sullivan WJ Jr: Base excision repair apurinic/apyrimidinic endonucleases in apicomplexan parasite Toxoplasma gondii. DNA Repair 2011, 10:466-475.
• [71]Quwailid MM, Hugill A, Dear N, Vizor L, Wells S, Horner E, Fuller S, Weedon J, McMath H, Woodman P, Edwards D, Campbell D, Rodger S, Carey J, Roberts A, Glenister P, Lalanne Z, Parkinson N, Coghill EL, McKeone R, Cox S, Willan J, Greenfield A, Keays D, Brady S, Spurr N, Gray I, Hunter J, Brown SD, Cox RD: A gene-driven ENU-based approach to generating an allelic series in any gene. Mamm Genome 2004, 15:585-591.
• [72]Timmermann B, Jarolim S, Russmayer H, Kerick M, Michel S, Kruger A, Bluemlein K, Laun P, Grillari J, Lehrach H, Breitenbach M, Ralser M: A new dominant peroxiredoxin allele identified by whole-genome re-sequencing of random mutagenized yeast causes oxidant-resistance and premature aging. Aging (Albany NY) 2010, 2:475-486.
• [73]Schumacher MA, Carter D, Scott DM, Roos DS, Ullman B, Brennan RG: Crystal structures of Toxoplasma gondii uracil phosphoribosyltransferase reveal the atomic basis of pyrimidine discrimination and prodrug binding. EMBO J 1998, 17:3219-3232.
• [74]Longley DB, Harkin DP, Johnston PG: 5-fluorouracil: mechanisms of action and clinical strategies. Nat Rev Cancer 2003, 3:330-338.