Pathogen Species Identification from Metagenomes in Ancient Remains: The Challenge of Identifying Human Pathogenic Species of Trypanosomatidae via Bioinformatic Tools
Identifieur interne : 000C94 ( Pmc/Corpus ); précédent : 000C93; suivant : 000C95Pathogen Species Identification from Metagenomes in Ancient Remains: The Challenge of Identifying Human Pathogenic Species of Trypanosomatidae via Bioinformatic Tools
Auteurs : Denis Sereno ; Franck Dorkeld ; Mohammad Akhoundi ; Pascale PerrinSource :
- Genes [ 2073-4425 ] ; 2018.
Abstract
Accurate species identification from ancient DNA samples is a difficult task that would shed light on the evolutionary history of pathogenic microorganisms. The field of palaeomicrobiology has undoubtedly benefited from the advent of untargeted metagenomic approaches that use next-generation sequencing methodologies. Nevertheless, assigning ancient DNA at the species level is a challenging process. Recently, the gut microbiome analysis of three pre-Columbian Andean mummies (Santiago-Rodriguez et al., 2016) has called into question the identification of
Url:
DOI: 10.3390/genes9080418
PubMed: 30127280
PubMed Central: 6115896
Links to Exploration step
PMC:6115896Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Pathogen Species Identification from Metagenomes in Ancient Remains: The Challenge of Identifying Human Pathogenic Species of <italic>Trypanosomatidae</italic>
via Bioinformatic Tools</title>
<author><name sortKey="Sereno, Denis" sort="Sereno, Denis" uniqKey="Sereno D" first="Denis" last="Sereno">Denis Sereno</name>
<affiliation><nlm:aff id="af1-genes-09-00418">IRD, Montpellier University, InterTryp, 34394 Montpellier, France</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Dorkeld, Franck" sort="Dorkeld, Franck" uniqKey="Dorkeld F" first="Franck" last="Dorkeld">Franck Dorkeld</name>
<affiliation><nlm:aff id="af2-genes-09-00418">INRA-UMR 1062 CBGP (INRA, IRD, CIRAD), Montpellier SupAgro, Montferrier-sur-Lez, 34988 Languedoc Roussillon, France;<email>franck.dorkeld@inra.fr</email>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Akhoundi, Mohammad" sort="Akhoundi, Mohammad" uniqKey="Akhoundi M" first="Mohammad" last="Akhoundi">Mohammad Akhoundi</name>
<affiliation><nlm:aff id="af3-genes-09-00418">Parasitology-Mycology Department, Avicenne Hospital, AP-HP, 93000 Bobigny, France;<email>m.akhoundi@yahoo.com</email>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Perrin, Pascale" sort="Perrin, Pascale" uniqKey="Perrin P" first="Pascale" last="Perrin">Pascale Perrin</name>
<affiliation><nlm:aff id="af4-genes-09-00418">Montpellier University, IRD, CNRS, MIVEGEC, 34394 Montpellier, France</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">30127280</idno>
<idno type="pmc">6115896</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC6115896</idno>
<idno type="RBID">PMC:6115896</idno>
<idno type="doi">10.3390/genes9080418</idno>
<date when="2018">2018</date>
<idno type="wicri:Area/Pmc/Corpus">000C94</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000C94</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Pathogen Species Identification from Metagenomes in Ancient Remains: The Challenge of Identifying Human Pathogenic Species of <italic>Trypanosomatidae</italic>
via Bioinformatic Tools</title>
<author><name sortKey="Sereno, Denis" sort="Sereno, Denis" uniqKey="Sereno D" first="Denis" last="Sereno">Denis Sereno</name>
<affiliation><nlm:aff id="af1-genes-09-00418">IRD, Montpellier University, InterTryp, 34394 Montpellier, France</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Dorkeld, Franck" sort="Dorkeld, Franck" uniqKey="Dorkeld F" first="Franck" last="Dorkeld">Franck Dorkeld</name>
<affiliation><nlm:aff id="af2-genes-09-00418">INRA-UMR 1062 CBGP (INRA, IRD, CIRAD), Montpellier SupAgro, Montferrier-sur-Lez, 34988 Languedoc Roussillon, France;<email>franck.dorkeld@inra.fr</email>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Akhoundi, Mohammad" sort="Akhoundi, Mohammad" uniqKey="Akhoundi M" first="Mohammad" last="Akhoundi">Mohammad Akhoundi</name>
<affiliation><nlm:aff id="af3-genes-09-00418">Parasitology-Mycology Department, Avicenne Hospital, AP-HP, 93000 Bobigny, France;<email>m.akhoundi@yahoo.com</email>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Perrin, Pascale" sort="Perrin, Pascale" uniqKey="Perrin P" first="Pascale" last="Perrin">Pascale Perrin</name>
<affiliation><nlm:aff id="af4-genes-09-00418">Montpellier University, IRD, CNRS, MIVEGEC, 34394 Montpellier, France</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">Genes</title>
<idno type="eISSN">2073-4425</idno>
<imprint><date when="2018">2018</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><p>Accurate species identification from ancient DNA samples is a difficult task that would shed light on the evolutionary history of pathogenic microorganisms. The field of palaeomicrobiology has undoubtedly benefited from the advent of untargeted metagenomic approaches that use next-generation sequencing methodologies. Nevertheless, assigning ancient DNA at the species level is a challenging process. Recently, the gut microbiome analysis of three pre-Columbian Andean mummies (Santiago-Rodriguez et al., 2016) has called into question the identification of <italic>Leishmania</italic>
in South America. The accurate assignment would be important because it will provide some key elements that are linked to the evolutionary scenario for visceral leishmaniasis agents in South America. Here, we recovered the metagenomic data filed in the metagenomics RAST server (MG-RAST) to identify the different members of the <italic>Trypanosomatidae</italic>
family that have infected these ancient remains. For this purpose, we used the ultrafast metagenomic sequence classifier, based on an exact alignment of k-mers (Kraken) and Bowtie2, an ultrafast and memory-efficient tool for aligning sequencing reads to long reference sequences. The analyses, which have been conducted on the most exhaustive genomic database possible on <italic>Trypanosomatidae</italic>
, show that species assignments could be biased by a lack of some genomic sequences of <italic>Trypanosomatidae</italic>
species (strains). Nevertheless, our work raises the issue of possible co-infections by multiple members of the <italic>Trypanosomatidae</italic>
family in these three pre-Columbian mummies. In the three mummies, we show the presence of DNA that is reminiscent of a probable co-infection with <italic>Leptomonas seymouri</italic>
, a parasite of insect’s gut, and <italic>Lotmaria</italic>
.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Santiago Rodriguez, T M" uniqKey="Santiago Rodriguez T">T.M. Santiago-Rodriguez</name>
</author>
<author><name sortKey="Fornaciari, G" uniqKey="Fornaciari G">G. Fornaciari</name>
</author>
<author><name sortKey="Luciani, S" uniqKey="Luciani S">S. Luciani</name>
</author>
<author><name sortKey="Dowd, S E" uniqKey="Dowd S">S.E. Dowd</name>
</author>
<author><name sortKey="Toranzos, G A" uniqKey="Toranzos G">G.A. Toranzos</name>
</author>
<author><name sortKey="Marota, I" uniqKey="Marota I">I. Marota</name>
</author>
<author><name sortKey="Cano, R J" uniqKey="Cano R">R.J. Cano</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sereno, D" uniqKey="Sereno D">D. Sereno</name>
</author>
<author><name sortKey="Akhoundi, M" uniqKey="Akhoundi M">M. Akhoundi</name>
</author>
<author><name sortKey="Dorkeld, F" uniqKey="Dorkeld F">F. Dorkeld</name>
</author>
<author><name sortKey="Oury, B" uniqKey="Oury B">B. Oury</name>
</author>
<author><name sortKey="Momen, H" uniqKey="Momen H">H. Momen</name>
</author>
<author><name sortKey="Perrin, P" uniqKey="Perrin P">P. Perrin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gerszten, E" uniqKey="Gerszten E">E. Gerszten</name>
</author>
<author><name sortKey="Allison, M J" uniqKey="Allison M">M.J. Allison</name>
</author>
<author><name sortKey="Maguire, B" uniqKey="Maguire B">B. Maguire</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Costa, M A" uniqKey="Costa M">M.A. Costa</name>
</author>
<author><name sortKey="Matheson, C" uniqKey="Matheson C">C. Matheson</name>
</author>
<author><name sortKey="Lachetta, L" uniqKey="Lachetta L">L. Lachetta</name>
</author>
<author><name sortKey="Llagostera, A" uniqKey="Llagostera A">A. Llagostera</name>
</author>
<author><name sortKey="Appenzeller, O" uniqKey="Appenzeller O">O. Appenzeller</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Salter, S J" uniqKey="Salter S">S.J. Salter</name>
</author>
<author><name sortKey="Cox, M J" uniqKey="Cox M">M.J. Cox</name>
</author>
<author><name sortKey="Turek, E M" uniqKey="Turek E">E.M. Turek</name>
</author>
<author><name sortKey="Calus, S T" uniqKey="Calus S">S.T. Calus</name>
</author>
<author><name sortKey="Cookson, W O" uniqKey="Cookson W">W.O. Cookson</name>
</author>
<author><name sortKey="Moffatt, M F" uniqKey="Moffatt M">M.F. Moffatt</name>
</author>
<author><name sortKey="Turner, P" uniqKey="Turner P">P. Turner</name>
</author>
<author><name sortKey="Parkhill, J" uniqKey="Parkhill J">J. Parkhill</name>
</author>
<author><name sortKey="Loman, N J" uniqKey="Loman N">N.J. Loman</name>
</author>
<author><name sortKey="Walker, A W" uniqKey="Walker A">A.W. Walker</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Warinner, C" uniqKey="Warinner C">C. Warinner</name>
</author>
<author><name sortKey="Herbig, A" uniqKey="Herbig A">A. Herbig</name>
</author>
<author><name sortKey="Mann, A" uniqKey="Mann A">A. Mann</name>
</author>
<author><name sortKey="Fellows Yates, J A" uniqKey="Fellows Yates J">J.A. Fellows Yates</name>
</author>
<author><name sortKey="Weir, C L" uniqKey="Weir C">C.L. Weir</name>
</author>
<author><name sortKey="Burbano, H A" uniqKey="Burbano H">H.A. Burbano</name>
</author>
<author><name sortKey="Orlando, L" uniqKey="Orlando L">L. Orlando</name>
</author>
<author><name sortKey="Krause, J" uniqKey="Krause J">J. Krause</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Eisenhofer, R" uniqKey="Eisenhofer R">R. Eisenhofer</name>
</author>
<author><name sortKey="Cooper, A" uniqKey="Cooper A">A. Cooper</name>
</author>
<author><name sortKey="Weyrich, L S" uniqKey="Weyrich L">L.S. Weyrich</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Eisenhofer, R" uniqKey="Eisenhofer R">R. Eisenhofer</name>
</author>
<author><name sortKey="Weyrich, L S" uniqKey="Weyrich L">L.S. Weyrich</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Toranzos, G A" uniqKey="Toranzos G">G.A. Toranzos</name>
</author>
<author><name sortKey="Santiago Rodriguez, T M" uniqKey="Santiago Rodriguez T">T.M. Santiago-Rodriguez</name>
</author>
<author><name sortKey="Cano, R J" uniqKey="Cano R">R.J. Cano</name>
</author>
<author><name sortKey="Fornaciari, G" uniqKey="Fornaciari G">G. Fornaciari</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wood, D E" uniqKey="Wood D">D.E. Wood</name>
</author>
<author><name sortKey="Salzberg, S L" uniqKey="Salzberg S">S.L. Salzberg</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Langmead, B" uniqKey="Langmead B">B. Langmead</name>
</author>
<author><name sortKey="Salzberg, S L" uniqKey="Salzberg S">S.L. Salzberg</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Aslett, M" uniqKey="Aslett M">M. Aslett</name>
</author>
<author><name sortKey="Aurrecoechea, C" uniqKey="Aurrecoechea C">C. Aurrecoechea</name>
</author>
<author><name sortKey="Berriman, M" uniqKey="Berriman M">M. Berriman</name>
</author>
<author><name sortKey="Brestelli, J" uniqKey="Brestelli J">J. Brestelli</name>
</author>
<author><name sortKey="Brunk, B P" uniqKey="Brunk B">B.P. Brunk</name>
</author>
<author><name sortKey="Carrington, M" uniqKey="Carrington M">M. Carrington</name>
</author>
<author><name sortKey="Depledge, D P" uniqKey="Depledge D">D.P. Depledge</name>
</author>
<author><name sortKey="Fischer, F" uniqKey="Fischer F">F. Fischer</name>
</author>
<author><name sortKey="Gajria, B" uniqKey="Gajria B">B. Gajria</name>
</author>
<author><name sortKey="Gao, X" uniqKey="Gao X">X. Gao</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ondov, B D" uniqKey="Ondov B">B.D. Ondov</name>
</author>
<author><name sortKey="Bergman, N H" uniqKey="Bergman N">N.H. Bergman</name>
</author>
<author><name sortKey="Phillippy, A M" uniqKey="Phillippy A">A.M. Phillippy</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Li, H" uniqKey="Li H">H. Li</name>
</author>
<author><name sortKey="Handsaker, B" uniqKey="Handsaker B">B. Handsaker</name>
</author>
<author><name sortKey="Wysoker, A" uniqKey="Wysoker A">A. Wysoker</name>
</author>
<author><name sortKey="Fennell, T" uniqKey="Fennell T">T. Fennell</name>
</author>
<author><name sortKey="Ruan, J" uniqKey="Ruan J">J. Ruan</name>
</author>
<author><name sortKey="Homer, N" uniqKey="Homer N">N. Homer</name>
</author>
<author><name sortKey="Marth, G" uniqKey="Marth G">G. Marth</name>
</author>
<author><name sortKey="Abecasis, G" uniqKey="Abecasis G">G. Abecasis</name>
</author>
<author><name sortKey="Durbin, R" uniqKey="Durbin R">R. Durbin</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Breniere, S F" uniqKey="Breniere S">S.F. Brenière</name>
</author>
<author><name sortKey="Waleckx, E" uniqKey="Waleckx E">E. Waleckx</name>
</author>
<author><name sortKey="Barnabe, C" uniqKey="Barnabe C">C. Barnabé</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Aufderheide, A C" uniqKey="Aufderheide A">A.C. Aufderheide</name>
</author>
<author><name sortKey="Salo, W" uniqKey="Salo W">W. Salo</name>
</author>
<author><name sortKey="Madden, M" uniqKey="Madden M">M. Madden</name>
</author>
<author><name sortKey="Streitz, J" uniqKey="Streitz J">J. Streitz</name>
</author>
<author><name sortKey="Buikstra, J" uniqKey="Buikstra J">J. Buikstra</name>
</author>
<author><name sortKey="Guhl, F" uniqKey="Guhl F">F. Guhl</name>
</author>
<author><name sortKey="Arriaza, B" uniqKey="Arriaza B">B. Arriaza</name>
</author>
<author><name sortKey="Renier, C" uniqKey="Renier C">C. Renier</name>
</author>
<author><name sortKey="Wittmers, L E" uniqKey="Wittmers L">L.E. Wittmers</name>
</author>
<author><name sortKey="Fornaciari, G" uniqKey="Fornaciari G">G. Fornaciari</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tomasini, N" uniqKey="Tomasini N">N. Tomasini</name>
</author>
<author><name sortKey="Diosque, P" uniqKey="Diosque P">P. Diosque</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Akhoundi, M" uniqKey="Akhoundi M">M. Akhoundi</name>
</author>
<author><name sortKey="Downing, T" uniqKey="Downing T">T. Downing</name>
</author>
<author><name sortKey="Votypka, J" uniqKey="Votypka J">J. Votypka</name>
</author>
<author><name sortKey="Kuhls, K" uniqKey="Kuhls K">K. Kuhls</name>
</author>
<author><name sortKey="Lukes, J" uniqKey="Lukes J">J. Lukes</name>
</author>
<author><name sortKey="Cannet, A" uniqKey="Cannet A">A. Cannet</name>
</author>
<author><name sortKey="Ravel, C" uniqKey="Ravel C">C. Ravel</name>
</author>
<author><name sortKey="Marty, P" uniqKey="Marty P">P. Marty</name>
</author>
<author><name sortKey="Delaunay, P" uniqKey="Delaunay P">P. Delaunay</name>
</author>
<author><name sortKey="Kasbari, M" uniqKey="Kasbari M">M. Kasbari</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Barral, A" uniqKey="Barral A">A. Barral</name>
</author>
<author><name sortKey="Badaro, R" uniqKey="Badaro R">R. Badaro</name>
</author>
<author><name sortKey="Barral Netto, M" uniqKey="Barral Netto M">M. Barral-Netto</name>
</author>
<author><name sortKey="Grimaldi, G" uniqKey="Grimaldi G">G. Grimaldi</name>
</author>
<author><name sortKey="Momen, H" uniqKey="Momen H">H. Momen</name>
</author>
<author><name sortKey="Carvalho, E M" uniqKey="Carvalho E">E.M. Carvalho</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Barral, A" uniqKey="Barral A">A. Barral</name>
</author>
<author><name sortKey="Pedral Sampaio, D" uniqKey="Pedral Sampaio D">D. Pedral-Sampaio</name>
</author>
<author><name sortKey="Grimaldi, G" uniqKey="Grimaldi G">G. Grimaldi</name>
</author>
<author><name sortKey="Momen, H" uniqKey="Momen H">H. Momen</name>
</author>
<author><name sortKey="Macmahon Pratt, D" uniqKey="Macmahon Pratt D">D. MacMahon-Pratt</name>
</author>
<author><name sortKey="Ribeiro De Jesus, A" uniqKey="Ribeiro De Jesus A">A. Ribeiro de Jesus</name>
</author>
<author><name sortKey="Almeida, R" uniqKey="Almeida R">R. Almeida</name>
</author>
<author><name sortKey="Badaro, R" uniqKey="Badaro R">R. Badaro</name>
</author>
<author><name sortKey="Barral Netto, M" uniqKey="Barral Netto M">M. Barral-Netto</name>
</author>
<author><name sortKey="Carvalho, E M" uniqKey="Carvalho E">E.M. Carvalho</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kraeva, N" uniqKey="Kraeva N">N. Kraeva</name>
</author>
<author><name sortKey="Butenko, A" uniqKey="Butenko A">A. Butenko</name>
</author>
<author><name sortKey="Hlavacova, J" uniqKey="Hlavacova J">J. Hlavacova</name>
</author>
<author><name sortKey="Kostygov, A" uniqKey="Kostygov A">A. Kostygov</name>
</author>
<author><name sortKey="Myskova, J" uniqKey="Myskova J">J. Myskova</name>
</author>
<author><name sortKey="Grybchuk, D" uniqKey="Grybchuk D">D. Grybchuk</name>
</author>
<author><name sortKey="Lestinova, T" uniqKey="Lestinova T">T. Lestinova</name>
</author>
<author><name sortKey="Votypka, J" uniqKey="Votypka J">J. Votypka</name>
</author>
<author><name sortKey="Volf, P" uniqKey="Volf P">P. Volf</name>
</author>
<author><name sortKey="Opperdoes, F" uniqKey="Opperdoes F">F. Opperdoes</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ghosh, S" uniqKey="Ghosh S">S. Ghosh</name>
</author>
<author><name sortKey="Banerjee, P" uniqKey="Banerjee P">P. Banerjee</name>
</author>
<author><name sortKey="Sarkar, A" uniqKey="Sarkar A">A. Sarkar</name>
</author>
<author><name sortKey="Datta, S" uniqKey="Datta S">S. Datta</name>
</author>
<author><name sortKey="Chatterjee, M" uniqKey="Chatterjee M">M. Chatterjee</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tritschler, M" uniqKey="Tritschler M">M. Tritschler</name>
</author>
<author><name sortKey="Retschnig, G" uniqKey="Retschnig G">G. Retschnig</name>
</author>
<author><name sortKey="Ya Ez, O" uniqKey="Ya Ez O">O. Yañez</name>
</author>
<author><name sortKey="Williams, G R" uniqKey="Williams G">G.R. Williams</name>
</author>
<author><name sortKey="Neumann, P" uniqKey="Neumann P">P. Neumann</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Da Costa, A P" uniqKey="Da Costa A">A.P. da Costa</name>
</author>
<author><name sortKey="Costa, F B" uniqKey="Costa F">F.B. Costa</name>
</author>
<author><name sortKey="Soares, H S" uniqKey="Soares H">H.S. Soares</name>
</author>
<author><name sortKey="Ramirez, D G" uniqKey="Ramirez D">D.G. Ramirez</name>
</author>
<author><name sortKey="Mesquita, E T" uniqKey="Mesquita E">E.T. Mesquita</name>
</author>
<author><name sortKey="Gennari, S M" uniqKey="Gennari S">S.M. Gennari</name>
</author>
<author><name sortKey="Marcili, A" uniqKey="Marcili A">A. Marcili</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Martinez, D Y" uniqKey="Martinez D">D.Y. Martínez</name>
</author>
<author><name sortKey="Verdonck, K" uniqKey="Verdonck K">K. Verdonck</name>
</author>
<author><name sortKey="Kaye, P M" uniqKey="Kaye P">P.M. Kaye</name>
</author>
<author><name sortKey="Adaui, V" uniqKey="Adaui V">V. Adaui</name>
</author>
<author><name sortKey="Polman, K" uniqKey="Polman K">K. Polman</name>
</author>
<author><name sortKey="Llanos Cuentas, A" uniqKey="Llanos Cuentas A">A. Llanos-Cuentas</name>
</author>
<author><name sortKey="Dujardin, J C" uniqKey="Dujardin J">J.C. Dujardin</name>
</author>
<author><name sortKey="Boelaert, M" uniqKey="Boelaert M">M. Boelaert</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Dedet, J P" uniqKey="Dedet J">J.P. Dedet</name>
</author>
<author><name sortKey="Pratlong, F" uniqKey="Pratlong F">F. Pratlong</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">Genes (Basel)</journal-id>
<journal-id journal-id-type="iso-abbrev">Genes (Basel)</journal-id>
<journal-id journal-id-type="publisher-id">genes</journal-id>
<journal-title-group><journal-title>Genes</journal-title>
</journal-title-group>
<issn pub-type="epub">2073-4425</issn>
<publisher><publisher-name>MDPI</publisher-name>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">30127280</article-id>
<article-id pub-id-type="pmc">6115896</article-id>
<article-id pub-id-type="doi">10.3390/genes9080418</article-id>
<article-id pub-id-type="publisher-id">genes-09-00418</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>Pathogen Species Identification from Metagenomes in Ancient Remains: The Challenge of Identifying Human Pathogenic Species of <italic>Trypanosomatidae</italic>
via Bioinformatic Tools</article-title>
</title-group>
<contrib-group><contrib contrib-type="author"><contrib-id contrib-id-type="orcid" authenticated="true">https://orcid.org/0000-0002-0034-9291</contrib-id>
<name><surname>Sereno</surname>
<given-names>Denis</given-names>
</name>
<xref ref-type="aff" rid="af1-genes-09-00418">1</xref>
<xref rid="c1-genes-09-00418" ref-type="corresp">*</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Dorkeld</surname>
<given-names>Franck</given-names>
</name>
<xref ref-type="aff" rid="af2-genes-09-00418">2</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Akhoundi</surname>
<given-names>Mohammad</given-names>
</name>
<xref ref-type="aff" rid="af3-genes-09-00418">3</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Perrin</surname>
<given-names>Pascale</given-names>
</name>
<xref ref-type="aff" rid="af4-genes-09-00418">4</xref>
<xref rid="c1-genes-09-00418" ref-type="corresp">*</xref>
</contrib>
</contrib-group>
<aff id="af1-genes-09-00418"><label>1</label>
IRD, Montpellier University, InterTryp, 34394 Montpellier, France</aff>
<aff id="af2-genes-09-00418"><label>2</label>
INRA-UMR 1062 CBGP (INRA, IRD, CIRAD), Montpellier SupAgro, Montferrier-sur-Lez, 34988 Languedoc Roussillon, France;<email>franck.dorkeld@inra.fr</email>
</aff>
<aff id="af3-genes-09-00418"><label>3</label>
Parasitology-Mycology Department, Avicenne Hospital, AP-HP, 93000 Bobigny, France;<email>m.akhoundi@yahoo.com</email>
</aff>
<aff id="af4-genes-09-00418"><label>4</label>
Montpellier University, IRD, CNRS, MIVEGEC, 34394 Montpellier, France</aff>
<author-notes><corresp id="c1-genes-09-00418"><label>*</label>
Correspondence: <email>denis.sereno@ird.fr</email>
(D.S.); <email>pascale.perrin@umontpellier.fr</email>
(P.P.)</corresp>
</author-notes>
<pub-date pub-type="epub"><day>20</day>
<month>8</month>
<year>2018</year>
</pub-date>
<pub-date pub-type="collection"><month>8</month>
<year>2018</year>
</pub-date>
<volume>9</volume>
<issue>8</issue>
<elocation-id>418</elocation-id>
<history><date date-type="received"><day>06</day>
<month>7</month>
<year>2018</year>
</date>
<date date-type="accepted"><day>16</day>
<month>8</month>
<year>2018</year>
</date>
</history>
<permissions><copyright-statement>© 2018 by the authors.</copyright-statement>
<copyright-year>2018</copyright-year>
<license license-type="open-access"><license-p>Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">http://creativecommons.org/licenses/by/4.0/</ext-link>
).</license-p>
</license>
</permissions>
<abstract><p>Accurate species identification from ancient DNA samples is a difficult task that would shed light on the evolutionary history of pathogenic microorganisms. The field of palaeomicrobiology has undoubtedly benefited from the advent of untargeted metagenomic approaches that use next-generation sequencing methodologies. Nevertheless, assigning ancient DNA at the species level is a challenging process. Recently, the gut microbiome analysis of three pre-Columbian Andean mummies (Santiago-Rodriguez et al., 2016) has called into question the identification of <italic>Leishmania</italic>
in South America. The accurate assignment would be important because it will provide some key elements that are linked to the evolutionary scenario for visceral leishmaniasis agents in South America. Here, we recovered the metagenomic data filed in the metagenomics RAST server (MG-RAST) to identify the different members of the <italic>Trypanosomatidae</italic>
family that have infected these ancient remains. For this purpose, we used the ultrafast metagenomic sequence classifier, based on an exact alignment of k-mers (Kraken) and Bowtie2, an ultrafast and memory-efficient tool for aligning sequencing reads to long reference sequences. The analyses, which have been conducted on the most exhaustive genomic database possible on <italic>Trypanosomatidae</italic>
, show that species assignments could be biased by a lack of some genomic sequences of <italic>Trypanosomatidae</italic>
species (strains). Nevertheless, our work raises the issue of possible co-infections by multiple members of the <italic>Trypanosomatidae</italic>
family in these three pre-Columbian mummies. In the three mummies, we show the presence of DNA that is reminiscent of a probable co-infection with <italic>Leptomonas seymouri</italic>
, a parasite of insect’s gut, and <italic>Lotmaria</italic>
.</p>
</abstract>
<kwd-group><kwd><italic>Trypanosomatidae</italic>
</kwd>
<kwd>kraken taxonomic assignment tool</kwd>
<kwd>bowtie2 fast short reads aligner</kwd>
<kwd>ancient DNA</kwd>
<kwd>parasitome</kwd>
<kwd>co-infection</kwd>
</kwd-group>
</article-meta>
</front>
<body><sec id="sec1-genes-09-00418"><title>1. Introduction</title>
<p>Santiago-Rodriguez et al. reported for the first time some evidence of the occurrence of <italic>Leishmania</italic>
DNA in the guts of Andean mummies dating to pre-Columbian times, and they proposed an assignment to <italic>L. donovani</italic>
[<xref rid="B1-genes-09-00418" ref-type="bibr">1</xref>
]. In South America, the circulation of such <italic>Leishmania</italic>
species is currently unknown, and has never before been documented in human remains. Therefore, the proper identification of the <italic>Leishmania</italic>
species that would have infected mummies before Iberian colonization remains a major concern, and would bring new elements to the puzzle of the possible evolutionary scenarios [<xref rid="B2-genes-09-00418" ref-type="bibr">2</xref>
]. <italic>L. donovani</italic>
and <italic>L. infantum</italic>
are by far the most common <italic>Leishmania</italic>
species responsible for the visceral form (visceral leishmaniasis, VL) of the disease in both the Old World and the New World. So far, all of the cases of leishmaniasis described on pre-Columbian mummies are reminiscent of cutaneous (CL) or mucocutaneous (MCL) lesions: CL was observed on a mummy dating to a cultural group from 700–800 AD and found in a cemetery in Peru [<xref rid="B3-genes-09-00418" ref-type="bibr">3</xref>
], and MCL was found on four samples in the archaeological cemetery of Coyo Oriente (skulls approximately 1000 years old) located in the desert of San Pedro de Atacama, Northern Chile [<xref rid="B4-genes-09-00418" ref-type="bibr">4</xref>
]. These sites come from a time period predating European contacts. A <italic>Leishmania</italic>
infection was confirmed in these remains using a Polymerase chain reaction (PCR) approach [<xref rid="B4-genes-09-00418" ref-type="bibr">4</xref>
] (amplification of fragments of the inositol monophosphate dehydrogenase gene, the kinetoplast minicircle, the amino acid permease AAP13LD, and the adenylate kinase gene). Furthermore, the amplified sequences differed from those of <italic>L. donovani</italic>
.</p>
<p>The rise of NGS (next-generation sequencing) technologies has opened a new field of systematic investigations in metagenomics. When these new technologies are applied to samples of ancient human remains, they provide valuable information to scientists working on the evolutionary history of infectious diseases. To avoid environmental contamination and accurately authenticate ancient DNA, a number of precautions and rules have been enacted [<xref rid="B5-genes-09-00418" ref-type="bibr">5</xref>
,<xref rid="B6-genes-09-00418" ref-type="bibr">6</xref>
]. Among them: (i) the requirement of a dedicated ancient DNA laboratory to minimize and manage contamination, particularly during sample collection, and the risk of exogenous DNA contaminations from the laboratory; and (ii) microbial ancient DNA damage should exhibit patterns of DNA damage and fragmentation; nevertheless, these patterns vary according to source context and species, but also according to the workflow and enzymes used during library preparation. The characterization of gut microbiomes of pre-Columbian mummies by Santiago-Rodriguez et al. using unbiased metagenomic approaches gave the opportunity to detect traces of specific pathogens [<xref rid="B1-genes-09-00418" ref-type="bibr">1</xref>
]. Nevertheless, several issues have been raised regarding the age of sequenced DNA and the possibility of environmental contaminations [<xref rid="B7-genes-09-00418" ref-type="bibr">7</xref>
,<xref rid="B8-genes-09-00418" ref-type="bibr">8</xref>
]. The answer to these questions has been published [<xref rid="B9-genes-09-00418" ref-type="bibr">9</xref>
]. However, if environmental contamination is a crucial problem when working on bacteria, for human pathogenic <italic>Trypanosomatidae</italic>
, environmental contamination is unlikely to occur. Indeed, these parasitic eukaryotes require an adequate environment to survive and proliferate. They are not free-living microorganisms, and have not had the ability to survive very long outside their hosts. Thus, contamination with such organisms is unlikely to occur during manipulation and the laboratory processing of samples. DNA contamination of the mummies by contact with some infected insects that transmit them (sandfly, <italic>Psychodidae</italic>
-<italic>Phlebotominae</italic>
, for <italic>Leishmania</italic>
or <italic>Triatoma</italic>
, <italic>Reduviidae</italic>
-<italic>Triatominae</italic>
, for <italic>Trypanosoma cruzi)</italic>
cannot be ruled out. Thus, because of the very low probability of such contamination, the available metagenomic data [<xref rid="B1-genes-09-00418" ref-type="bibr">1</xref>
] are of great interest to detect the presence of parasites belonging to the <italic>Trypanosomatidae</italic>
familly in these ancient remains.</p>
<p>Basic Local Alignment Search Tool (BLAST) alignments, which rely on finding the best alignment to a panel of genomic sequences, were often the traditional approach to assign a taxonomic label to an unknown sequence. However, unambiguous assignment at the species level is very hard, and this tool is very expensive in central processing unit (CPU) time for NGS data analysis on local computers, even if new facilities are now available via the use of Cloud public infrastructure or computer cluster. To shed light on the causative agent of leishmaniasis infection (and more broadly on parasites belonging to the <italic>Trypanosomatidae</italic>
family) in these ancient remains, we used new software dedicated to metagenomic data analysis. Kraken, a bioinformatic program [<xref rid="B10-genes-09-00418" ref-type="bibr">10</xref>
], presents numerous advantages over other programs, including its speed of performing analysis on metagenomes. The identification at the species level by Kraken is based on the use of exact-match database queries of k-mers, rather than on alignment similarity. This new approach was applied to the metagenomic rapid annotation using subsystems technology (MG-RAST) pre-processed metagenomic data [<xref rid="B1-genes-09-00418" ref-type="bibr">1</xref>
], and the results were compared to those obtained with Bowtie2 [<xref rid="B11-genes-09-00418" ref-type="bibr">11</xref>
].</p>
</sec>
<sec id="sec2-genes-09-00418"><title>2. Materials and Methods</title>
<sec id="sec2dot1-genes-09-00418"><title>2.1. Data</title>
<p>The metagenomic data that we used of the gut microbiomes from three pre-Columbian Andean mummies (FI3, FI9, and FI12) were available in the MG-RAST server (<uri xlink:href="http://blog.mg-rast.org/">http://blog.mg-rast.org/</uri>
) (MGRAST IDs 4629033.3, 4630170.3 and 4626489.3, respectively).</p>
<p>We built a reference genome database composed of all of the complete <italic>Trypanosomatidae</italic>
genomes collected from the NCBI (<uri xlink:href="http://www.ncbi.nlm.nih.gov/genome/">http://www.ncbi.nlm.nih.gov/genome/</uri>
) and complemented by genomes from TriTrypDB release 37 (25 April 2018) (Kinetoplastid Genomics Resource <uri xlink:href="http://tritrypdb.org/tritrypdb/">http://tritrypdb.org/tritrypdb/</uri>
) [<xref rid="B12-genes-09-00418" ref-type="bibr">12</xref>
]. The total number of <italic>Trypanosomatidae</italic>
genomes in the database is 79 (<xref ref-type="app" rid="app1-genes-09-00418">Table S1</xref>
).</p>
</sec>
<sec id="sec2dot2-genes-09-00418"><title>2.2. Methods</title>
<p>We analyzed the data with Kraken (version 1.0), which is a system for ultrafast metagenomic sequence classification using exact alignment [<xref rid="B10-genes-09-00418" ref-type="bibr">10</xref>
] (<uri xlink:href="http://ccb.jhu.edu/software/kraken/">http://ccb.jhu.edu/software/kraken/</uri>
). It relies on the development of a database that contains records consisting of a k-mer and the lowest common ancestor (LCA) of all of the organisms whose genomes contain that k-mer [<xref rid="B10-genes-09-00418" ref-type="bibr">10</xref>
]. We built a reference database with the default parameters of Kraken (k = 31). We downloaded the pre-processed sequences (this step filters sequences based on length, the number of ambiguous bases, and quality value) from the MG-RAST server, and then compared them with the non-redundant, custom-built <italic>Trypanosomatidae</italic>
database. The results were visualized with a metagenomic visualization tool, Krona [<xref rid="B13-genes-09-00418" ref-type="bibr">13</xref>
]. To complete the analysis, the same set of metagenome data was analyzed with Bowtie2—a program for the rapid alignment of gapped reads—using the sensitive option [<xref rid="B11-genes-09-00418" ref-type="bibr">11</xref>
]. To extract the coverage from the Binary Alignment Map (BAM) alignment files, we used the samtools program version 1.08, command idxstats [<xref rid="B14-genes-09-00418" ref-type="bibr">14</xref>
]. The reference genome database was used for analyses conducted with the two methodologies.</p>
</sec>
</sec>
<sec id="sec3-genes-09-00418"><title>3. Results and Discussion</title>
<p>In a first attempt to compare the respective limits of the two tools in the species assignment, we focused on microorganisms belonging to the <italic>Trypanosomatidae</italic>
family. These were chosen for two reasons: first, because reported <italic>Leishmania</italic>
infection is already reported in some mummies [<xref rid="B1-genes-09-00418" ref-type="bibr">1</xref>
], and second, because of the endemicity of Chagas disease [<xref rid="B15-genes-09-00418" ref-type="bibr">15</xref>
] and leishmaniasis in this South American region [<xref rid="B16-genes-09-00418" ref-type="bibr">16</xref>
].</p>
<p>We conducted analyses for the three mummies on pre-processed data, as indicated in the methods section. To confirm that the detected traces of <italic>Trypanosomatidae</italic>
DNA are not the result of environmental contaminations, we used the data passed through the screening step (i.e., following the dereplication and the duplicate read inferred sequencing error estimation (DRISEE) steps) in MG-RAST. The results confirmed clearly a co-infection with at least two pathogenic Trypanosomatids and the fact that is not the result of a contamination</p>
<sec id="sec3dot1-genes-09-00418"><title>3.1. Trypanosomes and Chagas Disease</title>
<p>The genome representativity of our database is 55% (11 genomes out of approximately 20 <italic>Trypanosoma</italic>
species currently described) (<xref ref-type="fig" rid="genes-09-00418-f001">Figure 1</xref>
A). The reference database includes trypanosomes responsible for American trypanosomiasis (<italic>T. cruzi</italic>
) and trypanosomes responsible for human African trypanosomiasis and animal trypanosomiasis.</p>
<p>The analysis performed with Kraken reveals that a non-negligible proportion of reads (5%, 9%, and 10% for FI3, FI9, and FI12) is attributed to the <italic>Trypanosomatidae</italic>
family (<xref ref-type="fig" rid="genes-09-00418-f002">Figure 2</xref>
A). The genus <italic>Trypanosoma</italic>
is clearly identified in ancient DNA from the guts of the three mummies FI3, FI9, and FI12, with 63%, 69%, and 52% of the detected <italic>Trypanosomatidae</italic>
community, respectively (<xref ref-type="fig" rid="genes-09-00418-f002">Figure 2</xref>
B). As expected, almost all of the reads belonging to the genus <italic>Trypanosoma</italic>
can be attributed to <italic>T. cruzi</italic>
(<xref ref-type="fig" rid="genes-09-00418-f003">Figure 3</xref>
A–C) for the three mummies. These results corroborate the conclusion of Santiago-Rodriguez et al. [<xref rid="B1-genes-09-00418" ref-type="bibr">1</xref>
]. Currently, <italic>T. cruzi</italic>
is split into six genetic lineages or discrete typing units (DTUs) named TcI, TcIV, TcII, TcIII, TcV, and TcVI, respectively, and a seventh one called TcBat. For both FI3 and FI9, the highest number of reads matches the <italic>T. cruzi</italic>
strain Tula cl2 (<xref ref-type="app" rid="app1-genes-09-00418">Table S2</xref>
), which belongs to TcI (DTU I). This DTU is widely represented in the genomic database that we gathered (five of the 11 genomes currently available). For FI12, reads that match Tula cl2 are scarce (<xref ref-type="app" rid="app1-genes-09-00418">Table S2</xref>
). Nevertheless, because of the relatively low proportion (<10%) of reads that match with a genome filed in our database, it is probable that the DTU of the infecting <italic>T. cruzi</italic>
strain is not yet represented.</p>
<p>Analysis with Bowtie2 clearly confirmed a <italic>T. cruzi</italic>
infection (53%, 49%, and 57% for FI3, FI9, and FI12, respectively) in these mummies (<xref ref-type="fig" rid="genes-09-00418-f002">Figure 2</xref>
C). The majority of reads matched <italic>T. cruzi</italic>
strains Y and Tula cl2 (<xref ref-type="app" rid="app1-genes-09-00418">Table S3</xref>
) rather than other genomes filed in the database. Concerning the alignment with <italic>T. cruzi</italic>
Y, the reads matched 178 and 203 contigs in FI3 and FI9, respectively, out of 9821 (<xref ref-type="app" rid="app1-genes-09-00418">Table S4</xref>
). A large majority of reads (80% in FI3 and 90% in FI9) matched 52 identical contigs for both mummies. For <italic>T. cruzi</italic>
Tula cl2, reads of mummy FI3 and mummy FI9 matched 42 and 39 contigs, respectively, out of 5300, and almost 100% of the reads for both mummies matched 28 identical contigs (<xref ref-type="app" rid="app1-genes-09-00418">Table S4</xref>
). Overall, this set of results indicates that the infecting <italic>T. cruzi</italic>
strain is probably not in our reference database.</p>
<p>Most of the genomes of <italic>T. cruzi</italic>
strains are lacking: currently, only 11 complete genomes are available in the NCBI and/or TriTrypDB databases out of more than 1902 deposited strains to date [<xref rid="B17-genes-09-00418" ref-type="bibr">17</xref>
]. Likewise, the database we used contains no representative genomes belonging to DTUs III, IV, or VII (TcBat). Interestingly, ancient <italic>T. cruzi</italic>
DNA has also been identified in human mummies dating from the same period (Chinchorro culture) and in the same geographical region, Southern Peru [<xref rid="B18-genes-09-00418" ref-type="bibr">18</xref>
]. In Bolivia and Peru, strains belonging to the DTU I and DTU V clades are the main <italic>T. cruzi</italic>
strains isolated, followed by strains of DTUs IV, III, and VI [<xref rid="B15-genes-09-00418" ref-type="bibr">15</xref>
]. The knowledge of the complex evolutionary history of <italic>T. cruzi</italic>
, which involved genetic exchanges [<xref rid="B19-genes-09-00418" ref-type="bibr">19</xref>
], and the existence of hybrid DTUs will certainly benefit from the identification of the infecting strains at the DTU level from this ancient DNA material. Insight into the circulating strains in these mummies should generate important elements for the calibration of such a dynamically evolving scenario. Nevertheless, our analysis highlights that the ancient <italic>T. cruzi</italic>
DNA present in the three mummies cannot be assigned to <italic>T. cruzi</italic>
CL Brener (DTU VI).</p>
</sec>
<sec id="sec3dot2-genes-09-00418"><title>3.2. Leishmania and Leishmaniasis</title>
<p>All of the <italic>Leishmania</italic>
genomes that are currently available in NCBI or in TriTrypDB are included in the reference database. The database represents only 38% of genomes from 21 of the 54 <italic>Leishmania</italic>
species that have been currently identified [<xref rid="B20-genes-09-00418" ref-type="bibr">20</xref>
]. Microorganisms from the subgenera <italic>Leishmania</italic>
and <italic>Viannia</italic>
are well represented (<xref ref-type="fig" rid="genes-09-00418-f001">Figure 1</xref>
B). Human pathogenic <italic>Leishmania</italic>
spp. are far better represented in the database (14/20), with the subgenera <italic>Leishmania</italic>
and <italic>Viannia</italic>
well represented (80% and 83%, respectively). No genome from the subgenus <italic>Porcisia</italic>
(<italic>Paraleishmania</italic>
section) is available (<xref rid="genes-09-00418-t001" ref-type="table">Table 1</xref>
). This subgenus includes some human pathogenic <italic>Leishmania</italic>
sp., such as <italic>L. colombiensis</italic>
.</p>
<p>The NGS reads analyzed with the metagenomic sequence classifier Kraken highlight a probable co-infection with a <italic>Leishmania</italic>
parasite for the three mummies FI3, FI9, and FI12 (<xref ref-type="fig" rid="genes-09-00418-f002">Figure 2</xref>
B), with 3%, 3%, and 10% of the reads attributed to <italic>Trypanosomatidae</italic>
that belong to the genus <italic>Leishmania</italic>
, respectively. From <xref ref-type="fig" rid="genes-09-00418-f003">Figure 3</xref>
, it is clear that Kraken assigns reads to a large number of <italic>Leishmania</italic>
species, each with a low percentage. Clearly, no predominant <italic>Leishmania</italic>
spp. is detected, and <italic>L. donovani</italic>
does not appear as a probable infecting <italic>Leishmania</italic>
species in all of these mummies. Therefore, these results prompt us to question the <italic>Leishmania</italic>
species infecting the mummies. Since DNA originates from internal tissues, we have to look for the <italic>Leishmania</italic>
spp. That are currently known to affect mainly internal organs and cause VL, namely, <italic>L. donovani</italic>
, <italic>L. infantum</italic>
, <italic>L. tropica</italic>
, <italic>L. martiniquensis</italic>
, and <italic>L. colombiensis</italic>
. At this point, we can exclude an infection by the first four L. species cited above. The sequence of the <italic>L. colombiensis</italic>
genome (<italic>Paraleishmania</italic>
section) is currently not available, and may be the <italic>Leishmania</italic>
agent that we are seeking. Currently, in Peru, the country from which mummies originate, <italic>Leishmania</italic>
pathogens for humans include <italic>L. peruviana</italic>
, <italic>L. guyanensis</italic>
, <italic>L. amazonensis</italic>
, <italic>L. lainsoni</italic>
, and <italic>L. braziliensis</italic>
[<xref rid="B20-genes-09-00418" ref-type="bibr">20</xref>
]. Some <italic>Leishmania</italic>
species that usually cause cutaneous forms have, under some circumstances, the capacity to disseminate into internal organs [<xref rid="B21-genes-09-00418" ref-type="bibr">21</xref>
,<xref rid="B22-genes-09-00418" ref-type="bibr">22</xref>
]. Such unusual clinical presentation is frequently associated with an immunocompromised state or dysfunction of the T helper-mediated immune response [<xref rid="B2-genes-09-00418" ref-type="bibr">2</xref>
]. Our study clearly shows that a low proportion of reads are attributed to <italic>L. peruviana</italic>
, <italic>L. amazonensis</italic>
, <italic>L. braziliensis</italic>
, and <italic>L. guyanensis</italic>
(species causing CL). Such low frequency in read assignment calls into question the validity of the identification (<xref ref-type="fig" rid="genes-09-00418-f002">Figure 2</xref>
) (<xref ref-type="app" rid="app1-genes-09-00418">Table S2</xref>
). Nevertheless, the occurrence of an infection by multiple <italic>Leishmania</italic>
sp. cannot be ruled out as well as infection by hybrid strains of <italic>Leishmania</italic>
[<xref rid="B20-genes-09-00418" ref-type="bibr">20</xref>
]. From these results, it is clear that to be accurate in its species assignment, the reference database used by Kraken needs to be highly representative of the described species and strains. The more exhaustive the database, the better the assignment accuracy will be.</p>
<p>The alignment program Bowtie2 detects the presence of <italic>Leishmania</italic>
in three mummies (FI3, FI9, and FI12) (<xref ref-type="fig" rid="genes-09-00418-f002">Figure 2</xref>
C). However, unlike the mummy FI12, in the mummies FI3 and FI9, the reads match a large panel of <italic>Leishmania</italic>
sp., especially in the mummy FI3 (<xref ref-type="app" rid="app1-genes-09-00418">Table S3</xref>
). Then, it is not possible to unequivocally assign DNA to a specific <italic>Leishmania</italic>
species and <italic>L. donovani</italic>
in particular due to the number of genomes of different <italic>Leishmania</italic>
spp. matched.</p>
<p>Our analysis reveals that the genome of the infecting <italic>Leishmania</italic>
species, whose DNA is detected in these ancient remains, is probably not present in our reference database. We may consider the occurrence of an ancient Leishmania species. Nevertheless, since our reference database is not exhaustive, we cannot rule out the presence of <italic>Leishmania</italic>
sp. causing CL or MCL or the co-infection of multiple <italic>Leishmania</italic>
spp. Still, the detection of <italic>Leishmania</italic>
DNA is clearly not the result of an artefact or a contamination.</p>
</sec>
<sec id="sec3dot3-genes-09-00418"><title>3.3. Other Leishmaniinae (Crithidiatae)</title>
<p><italic>Trypanosoma</italic>
and <italic>Leishmania</italic>
are dixenous parasites, meaning that their life cycle includes an invertebrate as a first host and a vertebrate as a second host. Several genomes of monoxenous (one-host) <italic>Trypanosomatidae</italic>
are available and belong to the genera <italic>Angomonas</italic>
, <italic>Crithidia</italic>
, <italic>Leptomonas</italic>
, <italic>Lotmaria,</italic>
and <italic>Strigomonas</italic>
(<xref ref-type="app" rid="app1-genes-09-00418">Table S1</xref>
). The database genome representativity for <italic>Leptomonas</italic>
is 5%, with only two genomes available out of the 39 <italic>Leptomonas</italic>
species currently described (<xref ref-type="fig" rid="genes-09-00418-f001">Figure 1</xref>
C). The analysis performed on the three mummies revealed that a relatively high proportion of reads attributed to the family <italic>Trypanosomatidae</italic>
is reminiscent of the genus <italic>Leptomonas</italic>
(13%, 18%, and 13% for F13, FI9, and FI12, respectively) (<xref ref-type="fig" rid="genes-09-00418-f002">Figure 2</xref>
B). A large majority of these reads are assigned to <italic>L. seymouri</italic>
and not <italic>L. pyrrochoris</italic>
after an analysis with Kraken (<xref ref-type="fig" rid="genes-09-00418-f003">Figure 3</xref>
A–C) (<xref ref-type="app" rid="app1-genes-09-00418">Table S2</xref>
). Nevertheless, because of the lack of available genomes for 37 out of 39 <italic>Leptomonas</italic>
sp., the assignment to <italic>L. seymouri</italic>
has to be taken with caution, and an assignment at the strain level cannot be made.</p>
<p>Bowtie2 also confirms the recurrent occurrence of <italic>L. seymouri</italic>
in all of the mummies (<xref ref-type="fig" rid="genes-09-00418-f002">Figure 2</xref>
C), and a high number of reads match the <italic>L. seymouri</italic>
strain BHU-1095 (<xref ref-type="app" rid="app1-genes-09-00418">Table S3</xref>
) in mummies FI3 and FI9 (56,245 and 125,577, respectively). In both mummies, the reads match with eight and nine scaffolds (out of 1216), respectively, and a large majority of reads match on only five scaffolds (<xref ref-type="app" rid="app1-genes-09-00418">Table S4</xref>
). <italic>Leptomonas</italic>
species are usually found in the gut of insects, but they have the potential to infect mammals as an opportunistic parasite. Nevertheless, their infective capacity in mammals seems to be limited to immunocompromised hosts [<xref rid="B23-genes-09-00418" ref-type="bibr">23</xref>
]. Interestingly, <italic>L. seymouri</italic>
has been repeatedly isolated from VL patients infected by <italic>L. donovani</italic>
in India [<xref rid="B24-genes-09-00418" ref-type="bibr">24</xref>
].</p>
<p>To our surprise, the analysis performed with Kraken reveals that an equivalent proportion of reads (2%) is assigned to <italic>L. passim</italic>
in the three mummies (<xref ref-type="fig" rid="genes-09-00418-f002">Figure 2</xref>
B and <xref ref-type="fig" rid="genes-09-00418-f003">Figure 3</xref>
). Via Bowtie2 analysis, a very low number of reads matches this genus in the mummies (<xref ref-type="fig" rid="genes-09-00418-f002">Figure 2</xref>
C and <xref ref-type="app" rid="app1-genes-09-00418">Table S3</xref>
). <italic>L. passim</italic>
is the founding member of this genus (<xref ref-type="fig" rid="genes-09-00418-f001">Figure 1</xref>
). This member of Crithidiatae is described as a common parasite of the honey bee <italic>Apis mellifera</italic>
[<xref rid="B25-genes-09-00418" ref-type="bibr">25</xref>
], which calls into question its presence in a human gut microbiome.</p>
</sec>
</sec>
<sec id="sec4-genes-09-00418"><title>4. Conclusions</title>
<p>The strength of our approach is to work with the most exhaustive reference database of <italic>Trypanosomatidae</italic>
genomic sequences that are possible, combined with two programs that are complementary in their metagenomic analysis approaches (k-mer searching versus local alignment). These allowed us to gain insight into the identity of the infecting <italic>Trypanosomatidae</italic>
agent. Analysis using the sequence classifier (Kraken) unambiguously confirmed <italic>T. cruzi</italic>
infection and undoubtedly <italic>Leishmania</italic>
infection in the three Andean mummies. In addition, our analysis provides new information on co-infection by at least two human pathogenic trypanosomatids, <italic>Leishmania</italic>
spp. and <italic>T. cruzi</italic>
, in all of the mummies with available metagenomes. This type of co-infection is known to occur in humans as well as in some wild mammals [<xref rid="B26-genes-09-00418" ref-type="bibr">26</xref>
,<xref rid="B27-genes-09-00418" ref-type="bibr">27</xref>
]. Unfortunately, it is not possible to go further in the assignment to the species level for <italic>Leishmania</italic>
and DTU level for <italic>T. cruzi</italic>
. This may be due to the lack of a number of <italic>Leishmania</italic>
genomes, particularly those of the <italic>Paraleishmania</italic>
section, and of members of <italic>T. cruzi</italic>
belonging to some DTUs. It also highlights a pattern of polyinfection coupled with an opportunistic trypanosomatid, i.e., <italic>L. seymouri</italic>
[<xref rid="B28-genes-09-00418" ref-type="bibr">28</xref>
]. Therefore, future studies with an exhaustive reference database are necessary to better understand the interrelationships that shape the microbial community and play a role in the evolution of the parasitome.</p>
</sec>
</body>
<back><app-group><app id="app1-genes-09-00418"><title>Supplementary Materials</title>
<p>The following are available online at <uri xlink:href="http://www.mdpi.com/2073-4425/9/8/418/s1">http://www.mdpi.com/2073-4425/9/8/418/s1</uri>
, Table S1: List of <italic>Trypanosomatidae</italic>
genome sequences available and included in the reference database. Table S2. (A) Results of Kraken analysis on preprocessed metagenomic data for the mummy FI3 (mgm_46.29033.3); (B) results of Kraken analysis on preprocessed metagenomic data for the mummy FI9 (mgm_46.30170.3); (C) results of Kraken analysis on preprocessed metagenomic data for the mummy FI12 (mgm_46.26489.3); Table S3. Results of Bowtie2 analysis concerning the three mummies FI3, FI9 and FI12. Table S4. Distribution of reads matching (a) <italic>Leptomonas seymouri</italic>
, (b) <italic>Trypanosoma cruzi</italic>
Y and (c) <italic>Trypanosoma cruzi</italic>
Tula cl2 on the contigs for the three mummies.</p>
<supplementary-material content-type="local-data" id="genes-09-00418-s001"><media xlink:href="genes-09-00418-s001.zip"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
</app>
</app-group>
<notes><title>Author Contributions</title>
<p>Writing-Original Draft Preparation, P.P. and D.S.; Bioinformatic analyses, F.D.; Writing-Review & Editing, P.P., D.S., F.D. and M.A.</p>
</notes>
<notes><title>Funding</title>
<p>The authors received no specific funding for this work.</p>
</notes>
<notes notes-type="COI-statement"><title>Conflicts of Interest</title>
<p>The authors declare no conflict of interest.</p>
</notes>
<ref-list><title>References</title>
<ref id="B1-genes-09-00418"><label>1.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Santiago-Rodriguez</surname>
<given-names>T.M.</given-names>
</name>
<name><surname>Fornaciari</surname>
<given-names>G.</given-names>
</name>
<name><surname>Luciani</surname>
<given-names>S.</given-names>
</name>
<name><surname>Dowd</surname>
<given-names>S.E.</given-names>
</name>
<name><surname>Toranzos</surname>
<given-names>G.A.</given-names>
</name>
<name><surname>Marota</surname>
<given-names>I.</given-names>
</name>
<name><surname>Cano</surname>
<given-names>R.J.</given-names>
</name>
</person-group>
<article-title>Taxonomic and predicted metabolic profiles of the human gut microbiome in pre-Columbian mummies</article-title>
<source>FEMS Microbiol. Ecol.</source>
<year>2016</year>
<volume>92</volume>
<pub-id pub-id-type="doi">10.1093/femsec/fiw182</pub-id>
<pub-id pub-id-type="pmid">27559027</pub-id>
</element-citation>
</ref>
<ref id="B2-genes-09-00418"><label>2.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sereno</surname>
<given-names>D.</given-names>
</name>
<name><surname>Akhoundi</surname>
<given-names>M.</given-names>
</name>
<name><surname>Dorkeld</surname>
<given-names>F.</given-names>
</name>
<name><surname>Oury</surname>
<given-names>B.</given-names>
</name>
<name><surname>Momen</surname>
<given-names>H.</given-names>
</name>
<name><surname>Perrin</surname>
<given-names>P.</given-names>
</name>
</person-group>
<article-title>What pre-Columbian mummies could teach us about South American leishmaniases?</article-title>
<source>Pathog. Dis.</source>
<year>2017</year>
<volume>75</volume>
<pub-id pub-id-type="doi">10.1093/femspd/ftx019</pub-id>
<pub-id pub-id-type="pmid">28423167</pub-id>
</element-citation>
</ref>
<ref id="B3-genes-09-00418"><label>3.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gerszten</surname>
<given-names>E.</given-names>
</name>
<name><surname>Allison</surname>
<given-names>M.J.</given-names>
</name>
<name><surname>Maguire</surname>
<given-names>B.</given-names>
</name>
</person-group>
<article-title>Paleopathology in South American mummies: A review and new findings</article-title>
<source>Pathobiology</source>
<year>2012</year>
<volume>79</volume>
<fpage>247</fpage>
<lpage>256</lpage>
<pub-id pub-id-type="doi">10.1159/000334087</pub-id>
<pub-id pub-id-type="pmid">22722564</pub-id>
</element-citation>
</ref>
<ref id="B4-genes-09-00418"><label>4.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Costa</surname>
<given-names>M.A.</given-names>
</name>
<name><surname>Matheson</surname>
<given-names>C.</given-names>
</name>
<name><surname>Lachetta</surname>
<given-names>L.</given-names>
</name>
<name><surname>Llagostera</surname>
<given-names>A.</given-names>
</name>
<name><surname>Appenzeller</surname>
<given-names>O.</given-names>
</name>
</person-group>
<article-title>Ancient leishmaniasis in a highland desert of Northern Chile</article-title>
<source>PLoS ONE</source>
<year>2009</year>
<volume>4</volume>
<pub-id pub-id-type="doi">10.1371/journal.pone.0006983</pub-id>
<pub-id pub-id-type="pmid">19746163</pub-id>
</element-citation>
</ref>
<ref id="B5-genes-09-00418"><label>5.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Salter</surname>
<given-names>S.J.</given-names>
</name>
<name><surname>Cox</surname>
<given-names>M.J.</given-names>
</name>
<name><surname>Turek</surname>
<given-names>E.M.</given-names>
</name>
<name><surname>Calus</surname>
<given-names>S.T.</given-names>
</name>
<name><surname>Cookson</surname>
<given-names>W.O.</given-names>
</name>
<name><surname>Moffatt</surname>
<given-names>M.F.</given-names>
</name>
<name><surname>Turner</surname>
<given-names>P.</given-names>
</name>
<name><surname>Parkhill</surname>
<given-names>J.</given-names>
</name>
<name><surname>Loman</surname>
<given-names>N.J.</given-names>
</name>
<name><surname>Walker</surname>
<given-names>A.W.</given-names>
</name>
</person-group>
<article-title>Reagent and laboratory contamination can critically impact sequence-based microbiome analyses</article-title>
<source>BMC Biol.</source>
<year>2014</year>
<volume>12</volume>
<elocation-id>87</elocation-id>
<pub-id pub-id-type="doi">10.1186/s12915-014-0087-z</pub-id>
<pub-id pub-id-type="pmid">25387460</pub-id>
</element-citation>
</ref>
<ref id="B6-genes-09-00418"><label>6.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Warinner</surname>
<given-names>C.</given-names>
</name>
<name><surname>Herbig</surname>
<given-names>A.</given-names>
</name>
<name><surname>Mann</surname>
<given-names>A.</given-names>
</name>
<name><surname>Fellows Yates</surname>
<given-names>J.A.</given-names>
</name>
<name><surname>Weir</surname>
<given-names>C.L.</given-names>
</name>
<name><surname>Burbano</surname>
<given-names>H.A.</given-names>
</name>
<name><surname>Orlando</surname>
<given-names>L.</given-names>
</name>
<name><surname>Krause</surname>
<given-names>J.</given-names>
</name>
</person-group>
<article-title>A robust framework for microbial archaeology</article-title>
<source>Annu. Rev. Genom. Hum. Genet.</source>
<year>2017</year>
<volume>18</volume>
<fpage>321</fpage>
<lpage>356</lpage>
<pub-id pub-id-type="doi">10.1146/annurev-genom-091416-035526</pub-id>
<pub-id pub-id-type="pmid">28460196</pub-id>
</element-citation>
</ref>
<ref id="B7-genes-09-00418"><label>7.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Eisenhofer</surname>
<given-names>R.</given-names>
</name>
<name><surname>Cooper</surname>
<given-names>A.</given-names>
</name>
<name><surname>Weyrich</surname>
<given-names>L.S.</given-names>
</name>
</person-group>
<article-title>Reply to Santiago-Rodriguez et al.: Proper authentication of ancient DNA is essential</article-title>
<source>FEMS Microbiol. Ecol.</source>
<year>2017</year>
<volume>93</volume>
<pub-id pub-id-type="doi">10.1093/femsec/fix042</pub-id>
<pub-id pub-id-type="pmid">28369385</pub-id>
</element-citation>
</ref>
<ref id="B8-genes-09-00418"><label>8.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Eisenhofer</surname>
<given-names>R.</given-names>
</name>
<name><surname>Weyrich</surname>
<given-names>L.S.</given-names>
</name>
</person-group>
<article-title>Proper authentication of ancient DNA is still essential</article-title>
<source>Genes</source>
<year>2018</year>
<elocation-id>9</elocation-id>
<pub-id pub-id-type="doi">10.3390/genes9030122</pub-id>
<pub-id pub-id-type="pmid">29495401</pub-id>
</element-citation>
</ref>
<ref id="B9-genes-09-00418"><label>9.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Toranzos</surname>
<given-names>G.A.</given-names>
</name>
<name><surname>Santiago-Rodriguez</surname>
<given-names>T.M.</given-names>
</name>
<name><surname>Cano</surname>
<given-names>R.J.</given-names>
</name>
<name><surname>Fornaciari</surname>
<given-names>G.</given-names>
</name>
</person-group>
<article-title>Proper authentication of ancient DNA is essential, yes; but so are undogmatic approaches</article-title>
<source>FEMS Microbiol. Ecol.</source>
<year>2017</year>
<volume>93</volume>
<pub-id pub-id-type="doi">10.1093/femsec/fix043</pub-id>
<pub-id pub-id-type="pmid">28402470</pub-id>
</element-citation>
</ref>
<ref id="B10-genes-09-00418"><label>10.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wood</surname>
<given-names>D.E.</given-names>
</name>
<name><surname>Salzberg</surname>
<given-names>S.L.</given-names>
</name>
</person-group>
<article-title>Kraken: Ultrafast metagenomic sequence classification using exact alignment</article-title>
<source>Genome Biol.</source>
<year>2014</year>
<volume>15</volume>
<pub-id pub-id-type="doi">10.1186/gb-2014-15-3-r46</pub-id>
<pub-id pub-id-type="pmid">24580807</pub-id>
</element-citation>
</ref>
<ref id="B11-genes-09-00418"><label>11.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Langmead</surname>
<given-names>B.</given-names>
</name>
<name><surname>Salzberg</surname>
<given-names>S.L.</given-names>
</name>
</person-group>
<article-title>Fast gapped-read alignment with Bowtie2</article-title>
<source>Nat. Methods</source>
<year>2012</year>
<volume>9</volume>
<fpage>357</fpage>
<lpage>359</lpage>
<pub-id pub-id-type="doi">10.1038/nmeth.1923</pub-id>
<pub-id pub-id-type="pmid">22388286</pub-id>
</element-citation>
</ref>
<ref id="B12-genes-09-00418"><label>12.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Aslett</surname>
<given-names>M.</given-names>
</name>
<name><surname>Aurrecoechea</surname>
<given-names>C.</given-names>
</name>
<name><surname>Berriman</surname>
<given-names>M.</given-names>
</name>
<name><surname>Brestelli</surname>
<given-names>J.</given-names>
</name>
<name><surname>Brunk</surname>
<given-names>B.P.</given-names>
</name>
<name><surname>Carrington</surname>
<given-names>M.</given-names>
</name>
<name><surname>Depledge</surname>
<given-names>D.P.</given-names>
</name>
<name><surname>Fischer</surname>
<given-names>F.</given-names>
</name>
<name><surname>Gajria</surname>
<given-names>B.</given-names>
</name>
<name><surname>Gao</surname>
<given-names>X.</given-names>
</name>
<etal></etal>
</person-group>
<article-title>TriTrypDB: A functional genomic resource for the <italic>Trypanosomatidae</italic>
</article-title>
<source>Nucl. Acids Res.</source>
<year>2010</year>
<volume>38</volume>
<fpage>D457</fpage>
<lpage>D462</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gkp851</pub-id>
<pub-id pub-id-type="pmid">19843604</pub-id>
</element-citation>
</ref>
<ref id="B13-genes-09-00418"><label>13.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ondov</surname>
<given-names>B.D.</given-names>
</name>
<name><surname>Bergman</surname>
<given-names>N.H.</given-names>
</name>
<name><surname>Phillippy</surname>
<given-names>A.M.</given-names>
</name>
</person-group>
<article-title>Interactive metagenomic visualization in a Web browser</article-title>
<source>BMC Bioinform.</source>
<year>2011</year>
<volume>12</volume>
<elocation-id>385</elocation-id>
<pub-id pub-id-type="doi">10.1186/1471-2105-12-385</pub-id>
<pub-id pub-id-type="pmid">21961884</pub-id>
</element-citation>
</ref>
<ref id="B14-genes-09-00418"><label>14.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname>
<given-names>H.</given-names>
</name>
<name><surname>Handsaker</surname>
<given-names>B.</given-names>
</name>
<name><surname>Wysoker</surname>
<given-names>A.</given-names>
</name>
<name><surname>Fennell</surname>
<given-names>T.</given-names>
</name>
<name><surname>Ruan</surname>
<given-names>J.</given-names>
</name>
<name><surname>Homer</surname>
<given-names>N.</given-names>
</name>
<name><surname>Marth</surname>
<given-names>G.</given-names>
</name>
<name><surname>Abecasis</surname>
<given-names>G.</given-names>
</name>
<name><surname>Durbin</surname>
<given-names>R.</given-names>
</name>
<collab>1000 Genome Project Data Processing Subgroup</collab>
</person-group>
<article-title>The sequence alignment/map (SAM) format and SAMtools</article-title>
<source>Bioinformatics</source>
<year>2009</year>
<volume>25</volume>
<fpage>2078</fpage>
<lpage>2079</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/btp352</pub-id>
<pub-id pub-id-type="pmid">19505943</pub-id>
</element-citation>
</ref>
<ref id="B15-genes-09-00418"><label>15.</label>
<element-citation publication-type="web"><person-group person-group-type="author"><collab>World Health Organization</collab>
</person-group>
<article-title>Fact Sheet on Chagas Disease</article-title>
<comment>Available online: <ext-link ext-link-type="uri" xlink:href="http://www.who.int/chagas/resources/en/">http://www.who.int/chagas/resources/en/</ext-link>
</comment>
<date-in-citation content-type="access-date" iso-8601-date="2018-04-15">(accessed on 15 April 2018)</date-in-citation>
</element-citation>
</ref>
<ref id="B16-genes-09-00418"><label>16.</label>
<element-citation publication-type="web"><person-group person-group-type="author"><collab>World Health Organization</collab>
</person-group>
<article-title>Fact Sheet on Leishmaniasis</article-title>
<comment>Available online: <ext-link ext-link-type="uri" xlink:href="http://www.who.int/leishmaniasis/resources/en/">http://www.who.int/leishmaniasis/resources/en/</ext-link>
</comment>
<date-in-citation content-type="access-date" iso-8601-date="2018-04-15">(accessed on 15 April 2018)</date-in-citation>
</element-citation>
</ref>
<ref id="B17-genes-09-00418"><label>17.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Brenière</surname>
<given-names>S.F.</given-names>
</name>
<name><surname>Waleckx</surname>
<given-names>E.</given-names>
</name>
<name><surname>Barnabé</surname>
<given-names>C.</given-names>
</name>
</person-group>
<article-title>Over six thousand <italic>Trypanosoma cruzi</italic>
strains classified into discrete typing units (DTUs): Attempt at an inventory</article-title>
<source>PLoS Negl. Trop. Dis.</source>
<year>2016</year>
<volume>10</volume>
<pub-id pub-id-type="doi">10.1371/journal.pntd.0004792</pub-id>
<pub-id pub-id-type="pmid">27571035</pub-id>
</element-citation>
</ref>
<ref id="B18-genes-09-00418"><label>18.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Aufderheide</surname>
<given-names>A.C.</given-names>
</name>
<name><surname>Salo</surname>
<given-names>W.</given-names>
</name>
<name><surname>Madden</surname>
<given-names>M.</given-names>
</name>
<name><surname>Streitz</surname>
<given-names>J.</given-names>
</name>
<name><surname>Buikstra</surname>
<given-names>J.</given-names>
</name>
<name><surname>Guhl</surname>
<given-names>F.</given-names>
</name>
<name><surname>Arriaza</surname>
<given-names>B.</given-names>
</name>
<name><surname>Renier</surname>
<given-names>C.</given-names>
</name>
<name><surname>Wittmers</surname>
<given-names>L.E.</given-names>
<suffix>Jr.</suffix>
</name>
<name><surname>Fornaciari</surname>
<given-names>G.</given-names>
</name>
<etal></etal>
</person-group>
<article-title>A 9000-year record of chagas’ disease</article-title>
<source>Proc. Natl. Acad. Sci. USA</source>
<year>2004</year>
<volume>101</volume>
<fpage>2034</fpage>
<lpage>2039</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.0307312101</pub-id>
<pub-id pub-id-type="pmid">14766963</pub-id>
</element-citation>
</ref>
<ref id="B19-genes-09-00418"><label>19.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tomasini</surname>
<given-names>N.</given-names>
</name>
<name><surname>Diosque</surname>
<given-names>P.</given-names>
</name>
</person-group>
<article-title>Evolution of <italic>Trypanosoma cruzi</italic>
: Clarifying hybridisations, mitochondrial introgressions and phylogenetic relationships between major lineages</article-title>
<source>Mem. Inst. Oswaldo Cruz</source>
<year>2015</year>
<volume>110</volume>
<fpage>403</fpage>
<lpage>413</lpage>
<pub-id pub-id-type="doi">10.1590/0074-02760140401</pub-id>
<pub-id pub-id-type="pmid">25807469</pub-id>
</element-citation>
</ref>
<ref id="B20-genes-09-00418"><label>20.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Akhoundi</surname>
<given-names>M.</given-names>
</name>
<name><surname>Downing</surname>
<given-names>T.</given-names>
</name>
<name><surname>Votypka</surname>
<given-names>J.</given-names>
</name>
<name><surname>Kuhls</surname>
<given-names>K.</given-names>
</name>
<name><surname>Lukes</surname>
<given-names>J.</given-names>
</name>
<name><surname>Cannet</surname>
<given-names>A.</given-names>
</name>
<name><surname>Ravel</surname>
<given-names>C.</given-names>
</name>
<name><surname>Marty</surname>
<given-names>P.</given-names>
</name>
<name><surname>Delaunay</surname>
<given-names>P.</given-names>
</name>
<name><surname>Kasbari</surname>
<given-names>M.</given-names>
</name>
<etal></etal>
</person-group>
<article-title><italic>Leishmania</italic>
infections: Molecular targets and diagnosis</article-title>
<source>Mol. Asp. Med.</source>
<year>2017</year>
<volume>57</volume>
<fpage>1</fpage>
<lpage>29</lpage>
<pub-id pub-id-type="doi">10.1016/j.mam.2016.11.012</pub-id>
<pub-id pub-id-type="pmid">28159546</pub-id>
</element-citation>
</ref>
<ref id="B21-genes-09-00418"><label>21.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Barral</surname>
<given-names>A.</given-names>
</name>
<name><surname>Badaro</surname>
<given-names>R.</given-names>
</name>
<name><surname>Barral-Netto</surname>
<given-names>M.</given-names>
<suffix>Jr.</suffix>
</name>
<name><surname>Grimaldi</surname>
<given-names>G.</given-names>
</name>
<name><surname>Momen</surname>
<given-names>H.</given-names>
</name>
<name><surname>Carvalho</surname>
<given-names>E.M.</given-names>
</name>
</person-group>
<article-title>Isolation of <italic>Leishmania mexicana amazonensis</italic>
from the bone marrow in a case of American visceral leishmaniasis</article-title>
<source>Am. J. Trop. Med. Hyg.</source>
<year>1986</year>
<volume>35</volume>
<fpage>732</fpage>
<lpage>734</lpage>
<pub-id pub-id-type="doi">10.4269/ajtmh.1986.35.732</pub-id>
<pub-id pub-id-type="pmid">3728794</pub-id>
</element-citation>
</ref>
<ref id="B22-genes-09-00418"><label>22.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Barral</surname>
<given-names>A.</given-names>
</name>
<name><surname>Pedral-Sampaio</surname>
<given-names>D.</given-names>
<suffix>Jr.</suffix>
</name>
<name><surname>Grimaldi</surname>
<given-names>G.</given-names>
</name>
<name><surname>Momen</surname>
<given-names>H.</given-names>
</name>
<name><surname>MacMahon-Pratt</surname>
<given-names>D.</given-names>
</name>
<name><surname>Ribeiro de Jesus</surname>
<given-names>A.</given-names>
</name>
<name><surname>Almeida</surname>
<given-names>R.</given-names>
</name>
<name><surname>Badaro</surname>
<given-names>R.</given-names>
</name>
<name><surname>Barral-Netto</surname>
<given-names>M.</given-names>
</name>
<name><surname>Carvalho</surname>
<given-names>E.M.</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Leishmaniasis in Bahia, Brazil: Evidence that <italic>Leishmania amazonensis</italic>
produces a wide spectrum of clinical diseases</article-title>
<source>Am. J. Trop. Med. Hyg.</source>
<year>1991</year>
<volume>44</volume>
<fpage>536</fpage>
<lpage>546</lpage>
<pub-id pub-id-type="doi">10.4269/ajtmh.1991.44.536</pub-id>
<pub-id pub-id-type="pmid">2063957</pub-id>
</element-citation>
</ref>
<ref id="B23-genes-09-00418"><label>23.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kraeva</surname>
<given-names>N.</given-names>
</name>
<name><surname>Butenko</surname>
<given-names>A.</given-names>
</name>
<name><surname>Hlavacova</surname>
<given-names>J.</given-names>
</name>
<name><surname>Kostygov</surname>
<given-names>A.</given-names>
</name>
<name><surname>Myskova</surname>
<given-names>J.</given-names>
</name>
<name><surname>Grybchuk</surname>
<given-names>D.</given-names>
</name>
<name><surname>Lestinova</surname>
<given-names>T.</given-names>
</name>
<name><surname>Votypka</surname>
<given-names>J.</given-names>
</name>
<name><surname>Volf</surname>
<given-names>P.</given-names>
</name>
<name><surname>Opperdoes</surname>
<given-names>F.</given-names>
</name>
<etal></etal>
</person-group>
<article-title><italic>Leptomonas seymouri</italic>
: Adaptations to the dixenous life cycle analyzed by genome sequencing, transcriptome profiling and co-infection with <italic>Leishmania donovani</italic>
</article-title>
<source>PLoS Pathog.</source>
<year>2015</year>
<volume>11</volume>
<pub-id pub-id-type="doi">10.1371/journal.ppat.1005127</pub-id>
<pub-id pub-id-type="pmid">26317207</pub-id>
</element-citation>
</ref>
<ref id="B24-genes-09-00418"><label>24.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ghosh</surname>
<given-names>S.</given-names>
</name>
<name><surname>Banerjee</surname>
<given-names>P.</given-names>
</name>
<name><surname>Sarkar</surname>
<given-names>A.</given-names>
</name>
<name><surname>Datta</surname>
<given-names>S.</given-names>
</name>
<name><surname>Chatterjee</surname>
<given-names>M.</given-names>
</name>
</person-group>
<article-title>Coinfection of <italic>Leptomonas seymouri</italic>
and <italic>Leishmania donovani</italic>
in Indian leishmaniasis</article-title>
<source>J. Clin. Microbiol.</source>
<year>2012</year>
<volume>50</volume>
<fpage>2774</fpage>
<lpage>2778</lpage>
<pub-id pub-id-type="doi">10.1128/JCM.00966-12</pub-id>
<pub-id pub-id-type="pmid">22622439</pub-id>
</element-citation>
</ref>
<ref id="B25-genes-09-00418"><label>25.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tritschler</surname>
<given-names>M.</given-names>
</name>
<name><surname>Retschnig</surname>
<given-names>G.</given-names>
</name>
<name><surname>Yañez</surname>
<given-names>O.</given-names>
</name>
<name><surname>Williams</surname>
<given-names>G.R.</given-names>
</name>
<name><surname>Neumann</surname>
<given-names>P.</given-names>
</name>
</person-group>
<article-title>Host sharing by the honey bee parasites <italic>Lotmaria passim</italic>
and <italic>Nosema ceranae</italic>
</article-title>
<source>Ecol. Evol.</source>
<year>2017</year>
<volume>7</volume>
<fpage>1850</fpage>
<lpage>1857</lpage>
<pub-id pub-id-type="doi">10.1002/ece3.2796</pub-id>
<pub-id pub-id-type="pmid">28331592</pub-id>
</element-citation>
</ref>
<ref id="B26-genes-09-00418"><label>26.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>da Costa</surname>
<given-names>A.P.</given-names>
</name>
<name><surname>Costa</surname>
<given-names>F.B.</given-names>
</name>
<name><surname>Soares</surname>
<given-names>H.S.</given-names>
</name>
<name><surname>Ramirez</surname>
<given-names>D.G.</given-names>
</name>
<name><surname>Mesquita</surname>
<given-names>E.T.</given-names>
</name>
<name><surname>Gennari</surname>
<given-names>S.M.</given-names>
</name>
<name><surname>Marcili</surname>
<given-names>A.</given-names>
</name>
</person-group>
<article-title><italic>Trypanosoma cruzi</italic>
and <italic>Leishmania infantum chagasi</italic>
infection in wild mammals from Maranhão state, Brazil</article-title>
<source>Vector Borne Zoonotic Dis.</source>
<year>2015</year>
<volume>15</volume>
<fpage>656</fpage>
<lpage>666</lpage>
<pub-id pub-id-type="doi">10.1089/vbz.2015.1771</pub-id>
<pub-id pub-id-type="pmid">26501369</pub-id>
</element-citation>
</ref>
<ref id="B27-genes-09-00418"><label>27.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Martínez</surname>
<given-names>D.Y.</given-names>
</name>
<name><surname>Verdonck</surname>
<given-names>K.</given-names>
</name>
<name><surname>Kaye</surname>
<given-names>P.M.</given-names>
</name>
<name><surname>Adaui</surname>
<given-names>V.</given-names>
</name>
<name><surname>Polman</surname>
<given-names>K.</given-names>
</name>
<name><surname>Llanos-Cuentas</surname>
<given-names>A.</given-names>
</name>
<name><surname>Dujardin</surname>
<given-names>J.C.</given-names>
</name>
<name><surname>Boelaert</surname>
<given-names>M.</given-names>
</name>
</person-group>
<article-title>Tegumentary leishmaniasis and coinfections other than HIV</article-title>
<source>PLoS Negl. Trop. Dis.</source>
<year>2018</year>
<volume>12</volume>
<pub-id pub-id-type="doi">10.1371/journal.pntd.0006125</pub-id>
<pub-id pub-id-type="pmid">29494584</pub-id>
</element-citation>
</ref>
<ref id="B28-genes-09-00418"><label>28.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dedet</surname>
<given-names>J.P.</given-names>
</name>
<name><surname>Pratlong</surname>
<given-names>F.</given-names>
</name>
</person-group>
<article-title><italic>Leishmania</italic>
, <italic>Trypanosoma</italic>
and monoxenous Trypanosomatids as emerging opportunistic agents</article-title>
<source>J. Eukaryot. Microbiol.</source>
<year>2000</year>
<volume>47</volume>
<fpage>37</fpage>
<lpage>39</lpage>
<pub-id pub-id-type="doi">10.1111/j.1550-7408.2000.tb00008.x</pub-id>
<pub-id pub-id-type="pmid">10651294</pub-id>
</element-citation>
</ref>
</ref-list>
</back>
<floats-group><fig id="genes-09-00418-f001" orientation="portrait" position="float"><label>Figure 1</label>
<caption><p>Proportion of available genomes (in blue) on the number of currently named species (in red) for the various genera belonging to the <italic>Trypanosomatidae</italic>
family. (<bold>A</bold>
) genus <italic>Trypanosoma</italic>
; (<bold>B</bold>
) genus <italic>Leishmania</italic>
; (<bold>C</bold>
) other <italic>Trypanosomatidae</italic>
genera.</p>
</caption>
<graphic xlink:href="genes-09-00418-g001"></graphic>
</fig>
<fig id="genes-09-00418-f002" orientation="portrait" position="float"><label>Figure 2</label>
<caption><p>Composition of gut microbiota for the three Andean mummies studied by Santiago-Rodriguez et al. [<xref rid="B1-genes-09-00418" ref-type="bibr">1</xref>
]: FI3, FI9, and FI12, respectively. (<bold>A</bold>
) Proportion of reads assigned to the <italic>Trypanosomatidae</italic>
family for each mummy after analysis with Kraken. Reads not matching <italic>Trypanosomatidae</italic>
sequences appear in grey. Reads matching <italic>Trypanosomatidae</italic>
sequences are flesh colored; (<bold>B</bold>
) In <italic>Trypanosomatidae</italic>
: numbers and percentages of reads assigned to the genera <italic>Trypanosoma</italic>
(purple), <italic>Leishmania</italic>
(blue), <italic>Leptomonas</italic>
(red), <italic>Lotmaria</italic>
(green). Reads from mummies FI3, FI9, and FI12 matching to any other <italic>Trypanosomatidae</italic>
are shown in orange; (<bold>C</bold>
) Results obtained by the Bowtie2 alignment sequence for mummies FI3, FI9, and FI12, respectively. For facility in comparison of the results gathered with both software the same color reference was used.</p>
</caption>
<graphic xlink:href="genes-09-00418-g002"></graphic>
</fig>
<fig id="genes-09-00418-f003" orientation="portrait" position="float"><label>Figure 3</label>
<caption><p>Profiles of microbiomes from three Andean mummies [<xref rid="B1-genes-09-00418" ref-type="bibr">1</xref>
]: FI3, FI9, and FI12, for <italic>Trypanosomatidae</italic>
. Circles represent taxonomic classification in ascending order up to the species level (outermost circle) with their relative abundance. These graphs were generated using the program Krona. Less abundant taxa are listed outside the charts.</p>
</caption>
<graphic xlink:href="genes-09-00418-g003a"></graphic>
<graphic xlink:href="genes-09-00418-g003b"></graphic>
<graphic xlink:href="genes-09-00418-g003c"></graphic>
</fig>
<table-wrap id="genes-09-00418-t001" orientation="portrait" position="float"><object-id pub-id-type="pii">genes-09-00418-t001_Table 1</object-id>
<label>Table 1</label>
<caption><p>Representativity of sequenced genomes available in relation to the number of known <italic>Leishmania</italic>
species that are pathogenic for humans, according the updated classification of Akhoundi et al. (2017) [<xref rid="B18-genes-09-00418" ref-type="bibr">18</xref>
].</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th align="center" valign="middle" style="border-top:solid thin;border-bottom:solid thin" rowspan="1" colspan="1"></th>
<th align="center" valign="middle" style="border-top:solid thin;border-bottom:solid thin" rowspan="1" colspan="1"><italic>Leishmania</italic>
</th>
<th align="center" valign="middle" style="border-top:solid thin;border-bottom:solid thin" rowspan="1" colspan="1"><italic>Viannia</italic>
</th>
<th align="center" valign="middle" style="border-top:solid thin;border-bottom:solid thin" rowspan="1" colspan="1"><italic>Sauroleishmania</italic>
</th>
<th align="center" valign="middle" style="border-top:solid thin;border-bottom:solid thin" rowspan="1" colspan="1"><italic>Mundinia</italic>
</th>
<th align="center" valign="middle" style="border-top:solid thin;border-bottom:solid thin" rowspan="1" colspan="1"><italic>L. australiensis</italic>
<break></break>
<italic>(L. macropodum)</italic>
</th>
<th align="center" valign="middle" style="border-top:solid thin;border-bottom:solid thin" rowspan="1" colspan="1"><italic>Porcisia</italic>
</th>
<th align="center" valign="middle" style="border-top:solid thin;border-bottom:solid thin" rowspan="1" colspan="1"><italic>Endotrypanum</italic>
</th>
</tr>
</thead>
<tbody><tr><td align="center" valign="middle" rowspan="1" colspan="1">Number of species pathogenic for humans</td>
<td align="center" valign="middle" rowspan="1" colspan="1">10</td>
<td align="center" valign="middle" rowspan="1" colspan="1">6</td>
<td align="center" valign="middle" rowspan="1" colspan="1">0</td>
<td align="center" valign="middle" rowspan="1" colspan="1">1</td>
<td align="center" valign="middle" rowspan="1" colspan="1">1</td>
<td align="center" valign="middle" rowspan="1" colspan="1">2</td>
<td align="center" valign="middle" rowspan="1" colspan="1">0</td>
</tr>
<tr><td align="center" valign="middle" rowspan="1" colspan="1">Number of genomes available</td>
<td align="center" valign="middle" rowspan="1" colspan="1">8</td>
<td align="center" valign="middle" rowspan="1" colspan="1">5</td>
<td align="center" valign="middle" rowspan="1" colspan="1">0</td>
<td align="center" valign="middle" rowspan="1" colspan="1">1</td>
<td align="center" valign="middle" rowspan="1" colspan="1">0</td>
<td align="center" valign="middle" rowspan="1" colspan="1">0</td>
<td align="center" valign="middle" rowspan="1" colspan="1">0</td>
</tr>
<tr><td align="center" valign="middle" style="border-bottom:solid thin" rowspan="1" colspan="1">Representativity in %</td>
<td align="center" valign="middle" style="border-bottom:solid thin" rowspan="1" colspan="1">80</td>
<td align="center" valign="middle" style="border-bottom:solid thin" rowspan="1" colspan="1">83</td>
<td align="center" valign="middle" style="border-bottom:solid thin" rowspan="1" colspan="1">0</td>
<td align="center" valign="middle" style="border-bottom:solid thin" rowspan="1" colspan="1">100</td>
<td align="center" valign="middle" style="border-bottom:solid thin" rowspan="1" colspan="1">0</td>
<td align="center" valign="middle" style="border-bottom:solid thin" rowspan="1" colspan="1">0</td>
<td align="center" valign="middle" style="border-bottom:solid thin" rowspan="1" colspan="1">0</td>
</tr>
</tbody>
</table>
</table-wrap>
</floats-group>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/MersV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000C94 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000C94 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Sante |area= MersV1 |flux= Pmc |étape= Corpus |type= RBID |clé= PMC:6115896 |texte= Pathogen Species Identification from Metagenomes in Ancient Remains: The Challenge of Identifying Human Pathogenic Species of Trypanosomatidae via Bioinformatic Tools }}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:30127280" \ | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \ | NlmPubMed2Wicri -a MersV1
![]() | This area was generated with Dilib version V0.6.33. | ![]() |