{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,9,30]],"date-time":"2026-09-30T23:29:00Z","timestamp":1790810940042,"version":"4.1.0"},"update-to":[{"DOI":"10.1371\/journal.pcbi.1008500","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2021,3,17]],"date-time":"2021-03-17T00:00:00Z","timestamp":1615939200000}}],"reference-count":37,"publisher":"Public Library of Science (PLoS)","issue":"3","license":[{"start":{"date-parts":[[2021,3,5]],"date-time":"2021-03-05T00:00:00Z","timestamp":1614902400000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":["www.ploscompbiol.org"],"crossmark-restriction":false},"short-container-title":["PLoS Comput Biol"],"abstract":"<jats:p><jats:italic>Mycobacterium tuberculosis<\/jats:italic> complex (MTC) CRISPR locus diversity has long been studied solely investigating the presence\/absence of a known set of spacers. Unveiling the genetic mechanisms of its evolution requires a more exhaustive reconstruction in a large amount of representative strains.<\/jats:p>\n<jats:p>In this article, we point out and resolve, with a new pipeline, the problem of CRISPR reconstruction based directly on short read sequences in <jats:italic>M<\/jats:italic>. <jats:italic>tuberculosis<\/jats:italic>. We first show that the process we set up, that we coin as \u201cCRISPRbuilder-TB\u201d (<jats:ext-link xmlns:xlink=\"http:\/\/www.w3.org\/1999\/xlink\" ext-link-type=\"uri\" xlink:href=\"https:\/\/github.com\/cguyeux\/CRISPRbuilder-TB\" xlink:type=\"simple\">https:\/\/github.com\/cguyeux\/CRISPRbuilder-TB<\/jats:ext-link>), allows an efficient reconstruction of simulated or real CRISPRs, even when including complex evolutionary steps like the insertions of mobile elements. Compared to more generalist tools, the whole process is much more precise and robust, and requires only minimal manual investigation. Second, we show that more than 1\/3 of the currently complete genomes available for this complex in the public databases contain largely erroneous CRISPR loci. Third, we highlight how both the classical experimental <jats:italic>in vitro<\/jats:italic> approach and the basic <jats:italic>in silico<\/jats:italic> spoligotyping provided by existing analytic tools miss a whole diversity of this locus in MTC, by not capturing duplications, spacer and direct repeats variants, and IS<jats:italic>6110<\/jats:italic> insertion locations. This description is extended in a second article that describes MTC-CRISPR diversity and suggests general rules for its evolution. This work opens perspectives for an in-depth exploration of <jats:italic>M<\/jats:italic>. <jats:italic>tuberculosis<\/jats:italic> CRISPR loci diversity and of mechanisms involved in its evolution and its functionality, as well as its adaptation to other CRISPR locus-harboring bacterial species.<\/jats:p>","DOI":"10.1371\/journal.pcbi.1008500","type":"journal-article","created":{"date-parts":[[2021,3,5]],"date-time":"2021-03-05T18:45:01Z","timestamp":1614969901000},"page":"e1008500","update-policy":"https:\/\/doi.org\/10.1371\/journal.pcbi.corrections_policy","source":"Crossref","is-referenced-by-count":19,"title":["CRISPRbuilder-TB: \u201cCRISPR-builder for tuberculosis\u201d. Exhaustive reconstruction of the CRISPR locus in mycobacterium tuberculosis complex using SRA"],"prefix":"10.1371","volume":"17","author":[{"ORCID":"https:\/\/orcid.org\/0000-0003-0195-4378","authenticated-orcid":true,"given":"Christophe","family":"Guyeux","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-4672-2140","authenticated-orcid":true,"given":"Christophe","family":"Sola","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Camille","family":"No\u00fbs","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-3749-7981","authenticated-orcid":true,"given":"Guislaine","family":"Refr\u00e9gier","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"340","published-online":{"date-parts":[[2021,3,5]]},"reference":[{"issue":"1","key":"pcbi.1008500.ref001","doi-asserted-by":"crossref","first-page":"841","DOI":"10.1186\/s12864-020-07178-6","article-title":"Unexpected diversity of CRISPR unveils some evolutionary patterns of repeated sequences in Mycobacterium tuberculosis","volume":"21","author":"G Refr\u00e9gier","year":"2020","journal-title":"BMC Genomics"},{"issue":"5","key":"pcbi.1008500.ref002","doi-asserted-by":"crossref","first-page":"1057","DOI":"10.1111\/j.1365-2958.1993.tb00976.x","article-title":"Nature of DNA polymorphism in the direct repeat cluster of Mycobacterium tuberculosis; application for strain differentiation by a novel typing method","volume":"10","author":"PM Groenen","year":"1993","journal-title":"Mol Microbiol."},{"issue":"4","key":"pcbi.1008500.ref003","doi-asserted-by":"crossref","first-page":"907","DOI":"10.1128\/JCM.35.4.907-914.1997","article-title":"Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology","volume":"35","author":"J Kamerbeek","year":"1997","journal-title":"J Clin Microbiol"},{"issue":"6","key":"pcbi.1008500.ref004","doi-asserted-by":"crossref","first-page":"23","DOI":"10.1186\/1471-2180-6-23","article-title":"Mycobacterium tuberculosis complex genetic diversity: mining the fourth international spoligotyping database (SpolDB4) for classification, Population Genetics, and Epidemiology","volume":"6","author":"K Brudey","year":"2006","journal-title":"BMC Microbiol"},{"issue":"1","key":"pcbi.1008500.ref005","first-page":"131","article-title":"Strain classification of Mycobacterium tuberculosis: congruence between large sequence polymorphisms and spoligotypes","volume":"15","author":"M Kato-Maeda","year":"2011","journal-title":"Int J Tuberc Lung Dis"},{"key":"pcbi.1008500.ref006","doi-asserted-by":"crossref","first-page":"4812","DOI":"10.1038\/ncomms5812","article-title":"A robust SNP barcode for typing Mycobacterium tuberculosis complex strains","volume":"5","author":"F Coll","year":"2014","journal-title":"Nat Commun."},{"issue":"12","key":"pcbi.1008500.ref007","doi-asserted-by":"crossref","first-page":"1535","DOI":"10.1038\/ng.3704","article-title":"Mycobacterium tuberculosis lineage 4 comprises globally distributed and geographically restricted sublineages","volume":"48","author":"D Stucki","year":"2016","journal-title":"Nat Genet."},{"issue":"1","key":"pcbi.1008500.ref008","doi-asserted-by":"crossref","first-page":"11597","DOI":"10.1038\/s41598-018-29986-3","article-title":"Evidence for Host-Bacterial Co-evolution via Genome Sequence Analysis of 480 Thai Mycobacterium tuberculosis Lineage 1 Isolates","volume":"8","author":"P Palittapongarnpim","year":"2018","journal-title":"Sci Rep."},{"issue":"1","key":"pcbi.1008500.ref009","doi-asserted-by":"crossref","first-page":"9227","DOI":"10.1038\/s41598-017-10018-5","article-title":"Evolutionary pathway analysis and unified classification of East Asian lineage of Mycobacterium tuberculosis","volume":"7","author":"E Shitikov","year":"2017","journal-title":"Sci Rep"},{"issue":"5","key":"pcbi.1008500.ref010","doi-asserted-by":"crossref","first-page":"325","DOI":"10.1089\/crispr.2018.0033","article-title":"Classification and Nomenclature of CRISPR-Cas Systems: Where from Here?","volume":"1","author":"KS Makarova","year":"2018","journal-title":"CRISPR J."},{"issue":"1","key":"pcbi.1008500.ref011","doi-asserted-by":"crossref","first-page":"1496","DOI":"10.1096\/fj.201800557RR","article-title":"Mycobacterium tuberculosis type III-A CRISPR\/Cas system crRNA and its maturation have atypical features","volume":"33","author":"W Wei","year":"2019","journal-title":"FASEB J."},{"issue":"1","key":"pcbi.1008500.ref012","doi-asserted-by":"crossref","first-page":"168","DOI":"10.1186\/s12864-017-3560-6","article-title":"Structure and variation of CRISPR and CRISPR-flanking regions in deleted-direct repeat region Mycobacterium tuberculosis complex strains","volume":"18","author":"PJ Freidlin","year":"2017","journal-title":"BMC Genomics"},{"issue":"14","key":"pcbi.1008500.ref013","doi-asserted-by":"crossref","first-page":"4865","DOI":"10.1073\/pnas.0305634101","article-title":"Functional and evolutionary genomics of Mycobacterium tuberculosis: insights from genomic deletions in 100 strains","volume":"101","author":"AG Tsolaki","year":"2004","journal-title":"Proc Natl Acad Sci U S A"},{"key":"pcbi.1008500.ref014","doi-asserted-by":"crossref","first-page":"2393","DOI":"10.1128\/JB.182.9.2393-2401.2000","article-title":"Genetic variation and evolutionary origin of the Direct repeat locus of Mycobacterium tuberculosis complex bacteria","volume":"182","author":"JDA van Embden","year":"2000","journal-title":"J Bacteriol."},{"issue":"1","key":"pcbi.1008500.ref015","doi-asserted-by":"crossref","first-page":"19","DOI":"10.1186\/s13073-016-0270-7","article-title":"SpoTyping: fast and accurate in silico Mycobacterium spoligotyping from sequence reads","volume":"8","author":"E Xia","year":"2016","journal-title":"Genome Med"},{"key":"pcbi.1008500.ref016","doi-asserted-by":"crossref","first-page":"209","DOI":"10.1186\/1471-2105-8-209","article-title":"CRISPR recognition tool (CRT): a tool for automatic detection of clustered regularly interspaced palindromic repeats","volume":"8","author":"C Bland","year":"2007","journal-title":"BMC Bioinformatics"},{"issue":"22","key":"pcbi.1008500.ref017","doi-asserted-by":"crossref","first-page":"2991","DOI":"10.1093\/bioinformatics\/bts544","article-title":"SpolPred: rapid and accurate prediction of Mycobacterium tuberculosis spoligotypes from short genomic sequences","volume":"28","author":"F Coll","year":"2012","journal-title":"Bioinformatics"},{"issue":"5","key":"pcbi.1008500.ref018","doi-asserted-by":"crossref","first-page":"821","DOI":"10.1101\/gr.074492.107","article-title":"Velvet: algorithms for de novo short read assembly using de Bruijn graphs","volume":"18","author":"DR Zerbino","year":"2008","journal-title":"Genome Res"},{"key":"pcbi.1008500.ref019","doi-asserted-by":"crossref","first-page":"W52","DOI":"10.1093\/nar\/gkm360","article-title":"CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats","volume":"35","author":"I Grissa","year":"2007","journal-title":"Nucleic Acids Res"},{"key":"pcbi.1008500.ref020","doi-asserted-by":"crossref","first-page":"18","DOI":"10.1186\/1471-2105-8-18","article-title":"PILER-CR: fast and accurate identification of CRISPR repeats","volume":"8","author":"RC Edgar","year":"2007","journal-title":"BMC Bioinformatics"},{"issue":"10","key":"pcbi.1008500.ref021","doi-asserted-by":"crossref","first-page":"e105","DOI":"10.1093\/nar\/gkt183","article-title":"Crass: identification and reconstruction of CRISPR from unassembled metagenomic data","volume":"41","author":"CT Skennerton","year":"2013","journal-title":"Nucleic Acids Res"},{"issue":"6","key":"pcbi.1008500.ref022","doi-asserted-by":"crossref","first-page":"461","DOI":"10.1089\/cmb.2015.0226","article-title":"CRISPR Detection From Short Reads Using Partial Overlap Graphs","volume":"23","author":"I Ben-Bassat","year":"2016","journal-title":"J Comput Biol."},{"issue":"17","key":"pcbi.1008500.ref023","doi-asserted-by":"crossref","first-page":"i520","DOI":"10.1093\/bioinformatics\/btw456","article-title":"Assemble CRISPRs from metagenomic sequencing data","volume":"32","author":"J Lei","year":"2016","journal-title":"Bioinformatics"},{"issue":"Suppl 5","key":"pcbi.1008500.ref024","doi-asserted-by":"crossref","first-page":"100","DOI":"10.1186\/s12918-018-0618-2","article-title":"On the reconstruction of the ancestral bacterial genomes in genus Mycobacterium and Brucella","volume":"12","author":"C Guyeux","year":"2018","journal-title":"BMC Syst Biol"},{"issue":"6","key":"pcbi.1008500.ref025","doi-asserted-by":"crossref","first-page":"e0216347","DOI":"10.1371\/journal.pone.0216347","article-title":"Evaluation of chloroplast genome annotation tools and application to analysis of the evolution of coffee species","volume":"14","author":"C Guyeux","year":"2019","journal-title":"PLoS One"},{"issue":"7","key":"pcbi.1008500.ref026","doi-asserted-by":"crossref","first-page":"e0130912","DOI":"10.1371\/journal.pone.0130912","article-title":"Genomics and Machine Learning for Taxonomy Consensus: The Mycobacterium tuberculosis Complex Paradigm","volume":"10","author":"J Aze","year":"2015","journal-title":"PLoS One"},{"issue":"6","key":"pcbi.1008500.ref027","doi-asserted-by":"crossref","first-page":"1908","DOI":"10.1128\/JCM.00025-15","article-title":"PhyResSE: a Web Tool Delineating Mycobacterium tuberculosis Antibiotic Resistance and Lineage from Whole-Genome Sequencing Data","volume":"53","author":"S Feuerriegel","year":"2015","journal-title":"J Clin Microbiol"},{"issue":"1","key":"pcbi.1008500.ref028","doi-asserted-by":"crossref","first-page":"613","DOI":"10.1186\/s12864-018-4988-z","article-title":"Identifying mixed Mycobacterium tuberculosis infections from whole genome sequence data","volume":"19","author":"B Sobkowiak","year":"2018","journal-title":"BMC Genomics"},{"key":"pcbi.1008500.ref029","article-title":"CRISPRCasFinder, an update of CRISRFinder, includes a portable version, enhanced performance and integrates search for Cas proteins","author":"D Couvin","year":"2018","journal-title":"Nucleic Acids Res."},{"issue":"6","key":"pcbi.1008500.ref030","doi-asserted-by":"crossref","first-page":"e1002441","DOI":"10.1371\/journal.pgen.1002441","article-title":"Diverse CRISPRs evolving in human microbiomes","volume":"8","author":"M Rho","year":"2012","journal-title":"PLoS Genet."},{"issue":"3","key":"pcbi.1008500.ref031","doi-asserted-by":"crossref","first-page":"346","DOI":"10.1016\/j.tube.2014.02.005","article-title":"PolyTB: A genomic variation map for Mycobacterium tuberculosis","volume":"94","author":"F Coll","year":"2014","journal-title":"Tuberculosis (Edinb)"},{"key":"pcbi.1008500.ref032","doi-asserted-by":"crossref","first-page":"12567","DOI":"10.1038\/srep12567","article-title":"Analysis of IS6110 insertion sites provide a glimpse into genome evolution of Mycobacterium tuberculosis","volume":"5","author":"T Roychowdhury","year":"2015","journal-title":"Sci Rep."},{"issue":"10","key":"pcbi.1008500.ref033","doi-asserted-by":"crossref","first-page":"eaat5869","DOI":"10.1126\/sciadv.aat5869","article-title":"Global expansion of Mycobacterium tuberculosis lineage 4 shaped by colonial migration and local adaptation","volume":"4","author":"OB Brynildsrud","year":"2018","journal-title":"Sci Adv."},{"issue":"1","key":"pcbi.1008500.ref034","doi-asserted-by":"crossref","first-page":"21","DOI":"10.3201\/eid2001.130652","article-title":"Progenitor \"Mycobacterium canettii\" clone responsible for lymph node tuberculosis epidemic, Djibouti","volume":"20","author":"Y Blouin","year":"2014","journal-title":"Emerg Infect Dis"},{"key":"pcbi.1008500.ref035","article-title":"Macro-geographical specificities of the prevailing tuberculosis epidemic as seen through SITVIT2, an updated version of the Mycobacterium tuberculosis genotyping database","author":"D Couvin","year":"2018","journal-title":"Infect Genet Evol"},{"key":"pcbi.1008500.ref036","doi-asserted-by":"crossref","first-page":"e5090","DOI":"10.7717\/peerj.5090","article-title":"MIRU-profiler: a rapid tool for determination of 24-loci MIRU-VNTR profiles from assembled genomes of Mycobacterium tuberculosis","volume":"6","author":"R Rajwani","year":"2018","journal-title":"PeerJ"},{"issue":"1","key":"pcbi.1008500.ref037","doi-asserted-by":"crossref","first-page":"69","DOI":"10.1093\/dnares\/dsu041","article-title":"MetaVelvet-SL: an extension of the Velvet assembler to a de novo metagenomic assembler utilizing supervised learning","volume":"22","author":"Afiahayati","year":"2015","journal-title":"DNA Res"}],"updated-by":[{"DOI":"10.1371\/journal.pcbi.1008500","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2021,3,17]],"date-time":"2021-03-17T00:00:00Z","timestamp":1615939200000}}],"container-title":["PLOS Computational Biology"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1008500","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2021,3,17]],"date-time":"2021-03-17T17:57:18Z","timestamp":1616003838000},"score":1,"resource":{"primary":{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1008500"}},"subtitle":[],"editor":[{"given":"Mihaela","family":"Pertea","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2021,3,5]]},"references-count":37,"journal-issue":{"issue":"3","published-online":{"date-parts":[[2021,3,5]]}},"URL":"https:\/\/doi.org\/10.1371\/journal.pcbi.1008500","relation":{"new_version":[{"id-type":"doi","id":"10.1371\/journal.pcbi.1008500","asserted-by":"object"}]},"ISSN":["1553-7358"],"issn-type":[{"value":"1553-7358","type":"electronic"}],"subject":[],"published":{"date-parts":[[2021,3,5]]}}}