Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Heteroplasmy in the Mitochondrial Genomes of Human Lice and Ticks Revealed by High Throughput Sequencing</title><author><name sortKey="Xiong, Haoyu" sort="Xiong, Haoyu" uniqKey="Xiong H" first="Haoyu" last="Xiong">Haoyu Xiong</name><affiliation><nlm:aff id="aff1"><addr-line>Parasitology Section, School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia</addr-line></nlm:aff></affiliation></author><author><name sortKey="Barker, Stephen C" sort="Barker, Stephen C" uniqKey="Barker S" first="Stephen C." last="Barker">Stephen C. Barker</name><affiliation><nlm:aff id="aff1"><addr-line>Parasitology Section, School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia</addr-line></nlm:aff></affiliation></author><author><name sortKey="Burger, Thomas D" sort="Burger, Thomas D" uniqKey="Burger T" first="Thomas D." last="Burger">Thomas D. Burger</name><affiliation><nlm:aff id="aff1"><addr-line>Parasitology Section, School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia</addr-line></nlm:aff></affiliation></author><author><name sortKey="Raoult, Didier" sort="Raoult, Didier" uniqKey="Raoult D" first="Didier" last="Raoult">Didier Raoult</name><affiliation><nlm:aff id="aff2"><addr-line>Unité de Recherche sur les Maladies Infectieuses et Tropicales Emergentes (URMITE), UMR CNRS 6236 IRD 198, Faculté de Mé decine, Mé diterranée Infection, Aix-Marseille Université, Marseille, France</addr-line></nlm:aff></affiliation></author><author><name sortKey="Shao, Renfu" sort="Shao, Renfu" uniqKey="Shao R" first="Renfu" last="Shao">Renfu Shao</name><affiliation><nlm:aff id="aff3"><addr-line>GeneCology Research Centre, Faculty of Science, Health, Education and Engineering, University of the Sunshine Coast, Maroochydore, Queensland, Australia</addr-line></nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">24058467</idno><idno type="pmc">3772822</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3772822</idno><idno type="RBID">PMC:3772822</idno><idno type="doi">10.1371/journal.pone.0073329</idno><date when="2013">2013</date><idno type="wicri:Area/Pmc/Corpus">002B15</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">002B15</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Heteroplasmy in the Mitochondrial Genomes of Human Lice and Ticks Revealed by High Throughput Sequencing</title><author><name sortKey="Xiong, Haoyu" sort="Xiong, Haoyu" uniqKey="Xiong H" first="Haoyu" last="Xiong">Haoyu Xiong</name><affiliation><nlm:aff id="aff1"><addr-line>Parasitology Section, School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia</addr-line></nlm:aff></affiliation></author><author><name sortKey="Barker, Stephen C" sort="Barker, Stephen C" uniqKey="Barker S" first="Stephen C." last="Barker">Stephen C. Barker</name><affiliation><nlm:aff id="aff1"><addr-line>Parasitology Section, School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia</addr-line></nlm:aff></affiliation></author><author><name sortKey="Burger, Thomas D" sort="Burger, Thomas D" uniqKey="Burger T" first="Thomas D." last="Burger">Thomas D. Burger</name><affiliation><nlm:aff id="aff1"><addr-line>Parasitology Section, School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia</addr-line></nlm:aff></affiliation></author><author><name sortKey="Raoult, Didier" sort="Raoult, Didier" uniqKey="Raoult D" first="Didier" last="Raoult">Didier Raoult</name><affiliation><nlm:aff id="aff2"><addr-line>Unité de Recherche sur les Maladies Infectieuses et Tropicales Emergentes (URMITE), UMR CNRS 6236 IRD 198, Faculté de Mé decine, Mé diterranée Infection, Aix-Marseille Université, Marseille, France</addr-line></nlm:aff></affiliation></author><author><name sortKey="Shao, Renfu" sort="Shao, Renfu" uniqKey="Shao R" first="Renfu" last="Shao">Renfu Shao</name><affiliation><nlm:aff id="aff3"><addr-line>GeneCology Research Centre, Faculty of Science, Health, Education and Engineering, University of the Sunshine Coast, Maroochydore, Queensland, Australia</addr-line></nlm:aff></affiliation></author></analytic><series><title level="j">PLoS ONE</title><idno type="eISSN">1932-6203</idno><imprint><date when="2013">2013</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>The typical mitochondrial (mt) genomes of bilateral animals consist of 37 genes on a single circular chromosome. The mt genomes of the human body louse, <italic>Pediculus humanus</italic>, and the human head louse, <italic>Pediculus capitis</italic>, however, are extensively fragmented and contain 20 minichromosomes, with one to three genes on each minichromosome. Heteroplasmy, i.e. nucleotide polymorphisms in the mt genome within individuals, has been shown to be significantly higher in the mt <italic>cox1</italic> gene of human lice than in humans and other animals that have the typical mt genomes. To understand whether the extent of heteroplasmy in human lice is associated with mt genome fragmentation, we sequenced the entire coding regions of all of the mt minichromosomes of six human body lice and six human head lice from Ethiopia, China and France with an Illumina HiSeq platform. For comparison, we also sequenced the entire coding regions of the mt genomes of seven species of ticks, which have the typical mitochondrial genome organization of bilateral animals. We found that the level of heteroplasmy varies significantly both among the human lice and among the ticks. The human lice from Ethiopia have significantly higher level of heteroplasmy than those from China and France (P<sub>t</sub><0.05). The tick, <italic>Amblyomma cajennense</italic>, has significantly higher level of heteroplasmy than other ticks (P<sub>t</sub><0.05). Our results indicate that heteroplasmy level can be substantially variable within a species and among closely related species, and does not appear to be determined by single factors such as genome fragmentation.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Rantanen, A" uniqKey="Rantanen A">A Rantanen</name></author><author><name sortKey="Larsson, N" uniqKey="Larsson N">N Larsson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lightowlers, Rn" uniqKey="Lightowlers R">RN Lightowlers</name></author><author><name sortKey="Chinnery, Pf" uniqKey="Chinnery P">PF Chinnery</name></author><author><name sortKey="Turnbull, Dm" uniqKey="Turnbull D">DM Turnbull</name></author><author><name sortKey="Howell, N" uniqKey="Howell N">N Howell</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Boore, Jl" uniqKey="Boore J">JL Boore</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Song, Gj" uniqKey="Song G">GJ Song</name></author><author><name sortKey="Lewis, V" uniqKey="Lewis V">V Lewis</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chiaratti, Mr" uniqKey="Chiaratti M">MR Chiaratti</name></author><author><name sortKey="Meirelles, Fvv" uniqKey="Meirelles F">FvV Meirelles</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kmiec, B" uniqKey="Kmiec B">B Kmiec</name></author><author><name sortKey="Woloszynska, M" uniqKey="Woloszynska M">M Woloszynska</name></author><author><name sortKey="Janska, H" uniqKey="Janska H">H Janska</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Light, Je" uniqKey="Light J">JE Light</name></author><author><name sortKey="Allen, Jm" uniqKey="Allen J">JM Allen</name></author><author><name sortKey="Long, Lm" uniqKey="Long L">LM Long</name></author><author><name sortKey="Carter, Te" uniqKey="Carter T">TE Carter</name></author><author><name sortKey="Barrow, L" uniqKey="Barrow L">L Barrow</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="He, Y" uniqKey="He Y">Y He</name></author><author><name sortKey="Wu, J" uniqKey="Wu J">J Wu</name></author><author><name sortKey="Dressman, Dc" uniqKey="Dressman D">DC Dressman</name></author><author><name sortKey="Iacobuzio Donahue, C" uniqKey="Iacobuzio Donahue C">C Iacobuzio-Donahue</name></author><author><name sortKey="Markowitz, Sd" uniqKey="Markowitz S">SD Markowitz</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Goto, H" uniqKey="Goto H">H Goto</name></author><author><name sortKey="Dickins, B" uniqKey="Dickins B">B Dickins</name></author><author><name sortKey="Afgan, E" uniqKey="Afgan E">E Afgan</name></author><author><name sortKey="Paul, Im" uniqKey="Paul I">IM Paul</name></author><author><name sortKey="Taylor, J" uniqKey="Taylor J">J Taylor</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Payne, Ba" uniqKey="Payne B">BA Payne</name></author><author><name sortKey="Wilson, Ij" uniqKey="Wilson I">IJ Wilson</name></author><author><name sortKey="Yu Wai Man, P" uniqKey="Yu Wai Man P">P Yu-Wai-Man</name></author><author><name sortKey="Coxhead, J" uniqKey="Coxhead J">J Coxhead</name></author><author><name sortKey="Deehan, D" uniqKey="Deehan D">D Deehan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chinnery, Pf" uniqKey="Chinnery P">PF Chinnery</name></author><author><name sortKey="Samuels, Dc" uniqKey="Samuels D">DC Samuels</name></author><author><name sortKey="Elso, J" uniqKey="Elso J">J Elso</name></author><author><name sortKey="Turnbull, Dm" uniqKey="Turnbull D">DM Turnbull</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Santos, C" uniqKey="Santos C">C Santos</name></author><author><name sortKey="Sierra, B" uniqKey="Sierra B">B Sierra</name></author><author><name sortKey="Alvarez, L" uniqKey="Alvarez L">L Alvarez</name></author><author><name sortKey="Ramos, A" uniqKey="Ramos A">A Ramos</name></author><author><name sortKey="Fernandez, E" uniqKey="Fernandez E">E Fernández</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Legros, F" uniqKey="Legros F">F Legros</name></author><author><name sortKey="Malka, F" uniqKey="Malka F">F Malka</name></author><author><name sortKey="Frachon, P" uniqKey="Frachon P">P Frachon</name></author><author><name sortKey="Manuel Rojo, Al" uniqKey="Manuel Rojo A">AL Manuel Rojo</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Boyce, Tm" uniqKey="Boyce T">TM Boyce</name></author><author><name sortKey="Zwick, Me" uniqKey="Zwick M">ME Zwick</name></author><author><name sortKey="Aquadro, Cf" uniqKey="Aquadro C">CF Aquadro</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nardi, F" uniqKey="Nardi F">F Nardi</name></author><author><name sortKey="Caraoelli, A" uniqKey="Caraoelli A">A Caraoelli</name></author><author><name sortKey="Fanciulli, Pp" uniqKey="Fanciulli P">PP Fanciulli</name></author><author><name sortKey="Dallai, R" uniqKey="Dallai R">R Dallai</name></author><author><name sortKey="Frati, F" uniqKey="Frati F">F Frati</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Marcade, I" uniqKey="Marcade I">I Marcadé</name></author><author><name sortKey="Cordaux, R" uniqKey="Cordaux R">R Cordaux</name></author><author><name sortKey="Doublet, V" uniqKey="Doublet V">V Doublet</name></author><author><name sortKey="Debenest, C" uniqKey="Debenest C">C Debenest</name></author><author><name sortKey="Bouchon, D" uniqKey="Bouchon D">D Bouchon</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kann, Lm" uniqKey="Kann L">LM Kann</name></author><author><name sortKey="Rosenblum, Eb" uniqKey="Rosenblum E">EB Rosenblum</name></author><author><name sortKey="Rand, Dm" uniqKey="Rand D">DM Rand</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Volz Lingenhohl, A" uniqKey="Volz Lingenhohl A">A Volz-Lingenhöhl</name></author><author><name sortKey="Solignac, M" uniqKey="Solignac M">M Solignac</name></author><author><name sortKey="Sperlich, D" uniqKey="Sperlich D">D Sperlich</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Leeuwen, Tv" uniqKey="Leeuwen T">TV Leeuwen</name></author><author><name sortKey="Bsartel, V" uniqKey="Bsartel V">V Bsartel</name></author><author><name sortKey="Pottelberge, Sv" uniqKey="Pottelberge S">SV Pottelberge</name></author><author><name sortKey="Nieuwenhjuyse, Pv" uniqKey="Nieuwenhjuyse P">PV Nieuwenhjuyse</name></author><author><name sortKey="Nauen, R" uniqKey="Nauen R">R Nauen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lavrov, Dv" uniqKey="Lavrov D">DV Lavrov</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shao, R" uniqKey="Shao R">R Shao</name></author><author><name sortKey="Kirkness, Ef" uniqKey="Kirkness E">EF Kirkness</name></author><author><name sortKey="Barker, Sc" uniqKey="Barker S">SC Barker</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shao, R" uniqKey="Shao R">R Shao</name></author><author><name sortKey="Zhu, Xq" uniqKey="Zhu X">XQ Zhu</name></author><author><name sortKey="Barker, Sc" uniqKey="Barker S">SC Barker</name></author><author><name sortKey="Herd, K" uniqKey="Herd K">K Herd</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Herd, K" uniqKey="Herd K">K Herd</name></author><author><name sortKey="Barker, Sc" uniqKey="Barker S">SC Barker</name></author><author><name sortKey="Shao, R" uniqKey="Shao R">R Shao</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Black, Wc" uniqKey="Black W">WC Black</name></author><author><name sortKey="Roehrdanz, Rl" uniqKey="Roehrdanz R">RL Roehrdanz</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shao, R" uniqKey="Shao R">R Shao</name></author><author><name sortKey="Aoki, Y" uniqKey="Aoki Y">Y Aoki</name></author><author><name sortKey="Mitani, H" uniqKey="Mitani H">H Mitani</name></author><author><name sortKey="Tabuchi, N" uniqKey="Tabuchi N">N Tabuchi</name></author><author><name sortKey="Barker, Sc" uniqKey="Barker S">SC Barker</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shao, R" uniqKey="Shao R">R Shao</name></author><author><name sortKey="Barker, Sc" uniqKey="Barker S">SC Barker</name></author><author><name sortKey="Mitani, H" uniqKey="Mitani H">H Mitani</name></author><author><name sortKey="Aoki, Y" uniqKey="Aoki Y">Y Aoki</name></author><author><name sortKey="Fukunaga, M" uniqKey="Fukunaga M">M Fukunaga</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shao, R" uniqKey="Shao R">R Shao</name></author><author><name sortKey="Barker, Sc" uniqKey="Barker S">SC Barker</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Burger, Td" uniqKey="Burger T">TD Burger</name></author><author><name sortKey="Shao, R" uniqKey="Shao R">R Shao</name></author><author><name sortKey="Beati, L" uniqKey="Beati L">L Beati</name></author><author><name sortKey="Miller, H" uniqKey="Miller H">H Miller</name></author><author><name sortKey="Barker, Sc" uniqKey="Barker S">SC Barker</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kambhampati, S" uniqKey="Kambhampati S">S Kambhampati</name></author><author><name sortKey="Smith, Pt" uniqKey="Smith P">PT Smith</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Simon, C" uniqKey="Simon C">C Simon</name></author><author><name sortKey="Frati, F" uniqKey="Frati F">F Frati</name></author><author><name sortKey="Beckenbach, A" uniqKey="Beckenbach A">A Beckenbach</name></author><author><name sortKey="Crespi, B" uniqKey="Crespi B">B Crespi</name></author><author><name sortKey="Liu, H" uniqKey="Liu H">H Liu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Altschul, Sf" uniqKey="Altschul S">SF Altschul</name></author><author><name sortKey="Gish, W" uniqKey="Gish W">W Gish</name></author><author><name sortKey="Miller, W" uniqKey="Miller W">W Miller</name></author><author><name sortKey="Myers, Ew" uniqKey="Myers E">EW Myers</name></author><author><name sortKey="Lipman, Dj" uniqKey="Lipman D">DJ Lipman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Mcginnis S Amp Madden, Tl" uniqKey="Mcginnis S Amp Madden T">TL McGinnis S & Madden</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lowe Tm Amp, Eddysr" uniqKey="Lowe Tm Amp E">EddySR Lowe TM &</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Laslett D Amp Canback, B" uniqKey="Laslett D Amp Canback B">B Laslett D & Canback</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Storey, Jd" uniqKey="Storey J">JD Storey</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Storey, Jd" uniqKey="Storey J">JD Storey</name></author><author><name sortKey="Tibshirani, R" uniqKey="Tibshirani R">R Tibshirani</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nakamura, K" uniqKey="Nakamura K">K Nakamura</name></author><author><name sortKey="Oshima, T" uniqKey="Oshima T">T Oshima</name></author><author><name sortKey="Morimoto, T" uniqKey="Morimoto T">T Morimoto</name></author><author><name sortKey="Ikeda, S" uniqKey="Ikeda S">S Ikeda</name></author><author><name sortKey="Yoshikawa, H" uniqKey="Yoshikawa H">H Yoshikawa</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Minoche, A" uniqKey="Minoche A">A Minoche</name></author><author><name sortKey="Dohm, Jc" uniqKey="Dohm J">JC Dohm</name></author><author><name sortKey="Himmelbauer, H" uniqKey="Himmelbauer H">H Himmelbauer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sachadyn, P" uniqKey="Sachadyn P">P Sachadyn</name></author><author><name sortKey="Zhang, Xm" uniqKey="Zhang X">XM Zhang</name></author><author><name sortKey="Clark, Ld" uniqKey="Clark L">LD Clark</name></author><author><name sortKey="Naviaux, Rk" uniqKey="Naviaux R">RK Naviaux</name></author><author><name sortKey="Heber Katz, E" uniqKey="Heber Katz E">E Heber-Katz</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Messenger, La" uniqKey="Messenger L">LA Messenger</name></author><author><name sortKey="Llewellyn, Ms" uniqKey="Llewellyn M">MS Llewellyn</name></author><author><name sortKey="Bhattacharyya, T" uniqKey="Bhattacharyya T">T Bhattacharyya</name></author><author><name sortKey="Franze, O" uniqKey="Franze O">O Franze</name></author><author><name sortKey="Lewis, Md" uniqKey="Lewis M">MD Lewis</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Quinlan, Ar" uniqKey="Quinlan A">AR Quinlan</name></author><author><name sortKey="Stewart, Da" uniqKey="Stewart D">DA Stewart</name></author><author><name sortKey="Stromberg, Mp" uniqKey="Stromberg M">MP Stromberg</name></author><author><name sortKey="Marth, Gt" uniqKey="Marth G">GT Marth</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sosa, Mx" uniqKey="Sosa M">MX Sosa</name></author><author><name sortKey="Sivakumar, Ka" uniqKey="Sivakumar K">KA Sivakumar</name></author><author><name sortKey="Maragh, S" uniqKey="Maragh S">S Maragh</name></author><author><name sortKey="Veeramachaneni, V" uniqKey="Veeramachaneni V">V Veeramachaneni</name></author><author><name sortKey="Hariharan, R" uniqKey="Hariharan R">R Hariharan</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Li, Mk" uniqKey="Li M">MK Li</name></author><author><name sortKey="Schonberg, A" uniqKey="Schonberg A">A Schönberg</name></author><author><name sortKey="Schaefer, M" uniqKey="Schaefer M">M Schaefer</name></author><author><name sortKey="Schroeder, R" uniqKey="Schroeder R">R Schroeder</name></author><author><name sortKey="Nasidze, I" uniqKey="Nasidze I">I Nasidze</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">PLoS One</journal-id><journal-id journal-id-type="iso-abbrev">PLoS ONE</journal-id><journal-id journal-id-type="publisher-id">plos</journal-id><journal-id journal-id-type="pmc">plosone</journal-id><journal-title-group><journal-title>PLoS ONE</journal-title></journal-title-group><issn pub-type="epub">1932-6203</issn><publisher><publisher-name>Public Library of Science</publisher-name><publisher-loc>San Francisco, USA</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">24058467</article-id><article-id pub-id-type="pmc">3772822</article-id><article-id pub-id-type="publisher-id">PONE-D-13-12541</article-id><article-id pub-id-type="doi">10.1371/journal.pone.0073329</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group></article-categories><title-group><article-title>Heteroplasmy in the Mitochondrial Genomes of Human Lice and Ticks Revealed by High Throughput Sequencing</article-title><alt-title alt-title-type="running-head">Mitochondrial Heteroplasmy in Human Lice and Ticks</alt-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Xiong</surname><given-names>Haoyu</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author"><name><surname>Barker</surname><given-names>Stephen C.</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author"><name><surname>Burger</surname><given-names>Thomas D.</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author"><name><surname>Raoult</surname><given-names>Didier</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author"><name><surname>Shao</surname><given-names>Renfu</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Parasitology Section, School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia</addr-line></aff><aff id="aff2"><label>2</label><addr-line>Unité de Recherche sur les Maladies Infectieuses et Tropicales Emergentes (URMITE), UMR CNRS 6236 IRD 198, Faculté de Mé decine, Mé diterranée Infection, Aix-Marseille Université, Marseille, France</addr-line></aff><aff id="aff3"><label>3</label><addr-line>GeneCology Research Centre, Faculty of Science, Health, Education and Engineering, University of the Sunshine Coast, Maroochydore, Queensland, Australia</addr-line></aff><contrib-group><contrib contrib-type="editor"><name><surname>Ouzounis</surname><given-names>Christos A.</given-names></name><role>Editor</role><xref ref-type="aff" rid="edit1"></xref></contrib></contrib-group><aff id="edit1"><addr-line>The Centre for Research and Technology, Hellas, Greece</addr-line></aff><author-notes><corresp id="cor1">* E-mail: <email>rshao@usc.edu.au</email></corresp><fn fn-type="conflict"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="con"><p>Conceived and designed the experiments: HX RS. Performed the experiments: HX TB RS. Contributed reagents/materials/analysis tools: DR SB RS. Wrote the paper: HX RS.</p></fn></author-notes><pub-date pub-type="collection"><year>2013</year></pub-date><pub-date pub-type="epub"><day>13</day><month>9</month><year>2013</year></pub-date><volume>8</volume><issue>9</issue><elocation-id>e73329</elocation-id><history><date date-type="received"><day>20</day><month>3</month><year>2013</year></date><date date-type="accepted"><day>10</day><month>7</month><year>2013</year></date></history><permissions><copyright-year>2013</copyright-year><copyright-holder>Shao et al</copyright-holder><license><license-p>This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p></license></permissions><abstract><p>The typical mitochondrial (mt) genomes of bilateral animals consist of 37 genes on a single circular chromosome. The mt genomes of the human body louse, <italic>Pediculus humanus</italic>, and the human head louse, <italic>Pediculus capitis</italic>, however, are extensively fragmented and contain 20 minichromosomes, with one to three genes on each minichromosome. Heteroplasmy, i.e. nucleotide polymorphisms in the mt genome within individuals, has been shown to be significantly higher in the mt <italic>cox1</italic> gene of human lice than in humans and other animals that have the typical mt genomes. To understand whether the extent of heteroplasmy in human lice is associated with mt genome fragmentation, we sequenced the entire coding regions of all of the mt minichromosomes of six human body lice and six human head lice from Ethiopia, China and France with an Illumina HiSeq platform. For comparison, we also sequenced the entire coding regions of the mt genomes of seven species of ticks, which have the typical mitochondrial genome organization of bilateral animals. We found that the level of heteroplasmy varies significantly both among the human lice and among the ticks. The human lice from Ethiopia have significantly higher level of heteroplasmy than those from China and France (P<sub>t</sub><0.05). The tick, <italic>Amblyomma cajennense</italic>, has significantly higher level of heteroplasmy than other ticks (P<sub>t</sub><0.05). Our results indicate that heteroplasmy level can be substantially variable within a species and among closely related species, and does not appear to be determined by single factors such as genome fragmentation.</p></abstract><funding-group><funding-statement>This work was supported by the Australian Research Council (DP0662755 to R.S. and DP120100240 to R.S. and S.C.B) (<ext-link ext-link-type="uri" xlink:href="http://www.arc.gov.au/">http://www.arc.gov.au/</ext-link>). R.S. also acknowledges the funding support from the Australian Government for an Australia-China Science & Research Fund Group Mission visit to China (ACSRF00980) (<ext-link ext-link-type="uri" xlink:href="http://www.innovation.gov.au/SCIENCE/INTERNATIONALCOLLABORATION/ACSRF/Pages/default.aspx">http://www.innovation.gov.au/SCIENCE/INTERNATIONALCOLLABORATION/ACSRF/Pages/default.aspx</ext-link>). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement></funding-group><counts><page-count count="9"></page-count></counts></article-meta></front><body><sec id="s1"><title>Introduction</title><p>Mitochondrial (mt) genomes exist usually in multiple copies per cell and vary in copy number among different cells and organisms. For instance, there are 30 to 40 copies of mt DNA molecules in the sperm of <italic>Caenorhabditis elegans</italic> whereas ∼250 copies of mt DNA molecules in the oocyte <xref ref-type="bibr" rid="pone.0073329-Lemire1">[1]</xref>. A mouse sperm has ∼100 copies of mt DNA molecules whereas a mouse oocyte can have up to 100,000 copies of mt DNA molecules <xref ref-type="bibr" rid="pone.0073329-Rantanen1">[2]</xref>. In humans, a sperm has 100 to 1,500 copies of mt DNA molecules whereas an oocyte can have up to 400,000 copies <xref ref-type="bibr" rid="pone.0073329-Lightowlers1">[3]</xref>, <xref ref-type="bibr" rid="pone.0073329-Boore1">[4]</xref>, <xref ref-type="bibr" rid="pone.0073329-Song1">[5]</xref>, <xref ref-type="bibr" rid="pone.0073329-Chiaratti1">[6]</xref>. The multiple copies of mt DNA molecules in a cell or an organism often are not identical in sequence. Heteroplasmy, i.e. variation in mt genome sequence within a cell, an organism or an individual, is common in bilateral animals but usually occurs at very low levels <xref ref-type="bibr" rid="pone.0073329-Kmiec1">[7]</xref>, <xref ref-type="bibr" rid="pone.0073329-Light1">[8]</xref>, <xref ref-type="bibr" rid="pone.0073329-He1">[9]</xref>, <xref ref-type="bibr" rid="pone.0073329-Goto1">[10]</xref>, <xref ref-type="bibr" rid="pone.0073329-Payne1">[11]</xref>.</p><p>Heteroplasmy has been most studied in humans in association with mt diseases, aging and cancers; nevertheless, heteroplasmy is present in health individuals too <xref ref-type="bibr" rid="pone.0073329-Goto1">[10]</xref>, <xref ref-type="bibr" rid="pone.0073329-Chinnery1">[12]</xref>, <xref ref-type="bibr" rid="pone.0073329-Santos1">[13]</xref>, <xref ref-type="bibr" rid="pone.0073329-Legros1">[14]</xref>. In arthropods, heteroplasmy in the mitochondrial A+T-rich regions (also called control regions) has been most studied <xref ref-type="bibr" rid="pone.0073329-Boyce1">[15]</xref>, <xref ref-type="bibr" rid="pone.0073329-Nardi1">[16]</xref>, <xref ref-type="bibr" rid="pone.0073329-Marcad1">[17]</xref>; heteroplasmy in mt tRNA genes and protein-coding genes has also been reported <xref ref-type="bibr" rid="pone.0073329-Kann1">[18]</xref>, <xref ref-type="bibr" rid="pone.0073329-VolzLingenhhl1">[19]</xref>, <xref ref-type="bibr" rid="pone.0073329-Leeuwen1">[20]</xref>. Andrea <italic>et al</italic>. (1992) showed deletion heteroplasmies in six protein-coding genes and four tRNA genes in <italic>Drosophila subobscura</italic><xref ref-type="bibr" rid="pone.0073329-VolzLingenhhl1">[19]</xref>. Leeuwen <italic>et al</italic>. (2008) showed that heteroplasmy in <italic>cytb</italic> gene was linked to resistance to insecticide bifenazate in mites; such resistance was inherited in a non-Mendelian manner <xref ref-type="bibr" rid="pone.0073329-Leeuwen1">[20]</xref>.</p><p>For bilateral animals, the mt genomes typically consist of a single circular chromosome that is ∼16 kb in size and has 37 genes <xref ref-type="bibr" rid="pone.0073329-Boore1">[4]</xref>, <xref ref-type="bibr" rid="pone.0073329-Lavrov1">[21]</xref>. The mt genomes of the human body louse, <italic>Pediculus humanus</italic>, and the human head louse, <italic>Pediculus capitis</italic>, however, contain 20 minichromosomes; each minichromosome is 3–4 kb in size and has one to three genes <xref ref-type="bibr" rid="pone.0073329-Shao1">[22]</xref>, <xref ref-type="bibr" rid="pone.0073329-Shao2">[23]</xref>. Herd <italic>et al</italic>. (2012) showed that the heteroplasmy level in the mt <italic>cox1</italic> gene of the human body lice and the human head lice is much higher than that in humans and other animals <xref ref-type="bibr" rid="pone.0073329-Herd1">[24]</xref>. At a lower coverage (∼5x) with Sanger sequencing, Herd <italic>et al</italic>. (2012) found 13 heteroplasmic sites in <italic>cox1</italic> gene in human body lice and four heteroplasmic sites in this gene in human head lice, in comparison to only one heteroplasmic site in this gene in humans revealed by high throughput sequencing <xref ref-type="bibr" rid="pone.0073329-He1">[9]</xref>. Herd <italic>et al</italic>. (2012) suggested that the high-level heteroplasmy might be linked to the fragmented mt genome organization in the human lice, possibly through increased recombination activities between minichromosomes. In the present study, we investigated further the link between heteroplasmy level and mt genome organization by sequencing in high coverage the entire coding regions (15.4 kb) of all of the mt minichromosomes of six human body lice and six human head lice collected in Ethiopia, China and France with an Illumina HiSeq platform. For comparison, we also sequenced the entire coding regions of the mt genomes of seven species of ticks. Ticks are arthropods, so are human lice. Ticks, however, have the typical mt genome organization of bilateral animals <xref ref-type="bibr" rid="pone.0073329-Black1">[25]</xref>, <xref ref-type="bibr" rid="pone.0073329-Shao3">[26]</xref>, <xref ref-type="bibr" rid="pone.0073329-Shao4">[27]</xref>, <xref ref-type="bibr" rid="pone.0073329-Shao5">[28]</xref>, <xref ref-type="bibr" rid="pone.0073329-Burger1">[29]</xref>. Comparison of human lice with ticks, thus, will allow us to test whether or not heteroplasmy level is indeed linked to changes in mt genome organization.</p></sec><sec sec-type="materials|methods" id="s2"><title>Materials and Methods</title><sec id="s2a"><title>Ethical statement</title><p>The research ethics committees of The University of Queensland and the University of the Sunshine Coast were consulted for the experimental work undertaken. No ethical approval or permits were required for conducting the research reported in this manuscript according the laws and regulations in Australia and Queensland.</p></sec><sec id="s2b"><title>Sample collection, DNA extraction, PCR amplification and Illumina sequencing</title><p>Total DNA was extracted individually from six human body lice, <italic>Pediculus humanus</italic>, six human head lice, <italic>Pediculus capitis</italic>, from Ethiopia, China and France, and seven species of ticks with DNeasy Tissue Kit (QIAGEN) (<xref ref-type="table" rid="pone-0073329-t001">Table 1</xref>). For the human lice, a forward primer (USF1PH1: <named-content content-type="gene">5′-GAAATTAAAATTTCAACAAATCTCAACTCG-3′</named-content>) and a reverse primer (PHR10: <named-content content-type="gene">5′-CCCCCCAAGCTATTTATAGCTTGGAGTATTAACGG-3′</named-content>) were designed from two conserved sequence blocks adjacent to the coding region of each mt minichromosome of the human body louse and the human head louse <xref ref-type="bibr" rid="pone.0073329-VolzLingenhhl1">[19]</xref>, <xref ref-type="bibr" rid="pone.0073329-Shao2">[23]</xref>. These primers were used to amplify the coding regions of the 20 mt minichromosomes of the human body lice and the human head lice. The PCR conditions were: 96°C for 2 minutes; 39 cycles of 98°C for 10 seconds, 62°C for 30 seconds, 72°C for 2.5 minutes; followed by 72°C for 5 minutes and then 25°C for 2 minutes. Each 25-µl PCR reaction contains 0.5 µl La Taq, 2 µl forward primer (10 µM), 2 µl reverse primer (10 µM), 4 µl dNTP (2.5 mM), 2.5 µl La Taq Buffer (10x), 1 µl DNA template, and 13 µl deionized water. Negative controls were always included for contamination check. PCR amplicons were purified and concentrated with the Wizard® SV Gel and PCR Clean-Up System (Promega) and then submitted to the Beijing Genome Institute (BGI) for deep sequencing. Sequencing libraries were constructed for each louse and tick sample at the BGI. For library construction, PCR amplicons were sheared into 180-bp fragments, purified and then end-repaired. Adapters were ligated to both ends of sheared fragments; DNA cluster preparation were then performed and the products were sequenced with an Illumina Hiseq 2000 platform at the BGI (details available at: <ext-link ext-link-type="uri" xlink:href="http://bgitechsolutions.com/service-solutions/services/genomics/de-novo-sequencing">http://bgitechsolutions.com/service-solutions/services/genomics/de-novo-sequencing</ext-link>; <ext-link ext-link-type="uri" xlink:href="http://support.illumina.com/sequencing/documentation.ilmn">http://support.illumina.com/sequencing/documentation.ilmn</ext-link>).</p><table-wrap id="pone-0073329-t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0073329.t001</object-id><label>Table 1</label><caption><title>Samples of human lice and ticks used in this study.</title></caption><alternatives><graphic id="pone-0073329-t001-1" xlink:href="pone.0073329.t001"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Sample</td><td align="left" rowspan="1" colspan="1">Country</td><td align="left" rowspan="1" colspan="1">Host</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Body louse, <italic>Pediculus humanus</italic> (B2470B1)</td><td align="left" rowspan="1" colspan="1">China</td><td align="left" rowspan="1" colspan="1">Human</td></tr><tr><td align="left" rowspan="1" colspan="1">Head louse, <italic>Pediculus capitis</italic> (B2471H5)</td><td align="left" rowspan="1" colspan="1">China</td><td align="left" rowspan="1" colspan="1">Human</td></tr><tr><td align="left" rowspan="1" colspan="1">Body louse, <italic>Pediculus humanus</italic> (B2557B)</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">Human</td></tr><tr><td align="left" rowspan="1" colspan="1">Head louse, <italic>Pediculus capitis</italic> (B2557H)</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">Human</td></tr><tr><td align="left" rowspan="1" colspan="1">Body louse, <italic>Pediculus humanus</italic> (B2558B)</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">Human</td></tr><tr><td align="left" rowspan="1" colspan="1">Head louse, <italic>Pediculus capitis</italic> (B2558H)</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">Human</td></tr><tr><td align="left" rowspan="1" colspan="1">Body louse, <italic>Pediculus humanus</italic> (B2516B2)</td><td align="left" rowspan="1" colspan="1">Ethiopia</td><td align="left" rowspan="1" colspan="1">Human</td></tr><tr><td align="left" rowspan="1" colspan="1">Head louse, <italic>Pediculus capitis</italic> (B2516H2)</td><td align="left" rowspan="1" colspan="1">Ethiopia</td><td align="left" rowspan="1" colspan="1">Human</td></tr><tr><td align="left" rowspan="1" colspan="1">Body louse, <italic>Pediculus humanus</italic> (B2560B)</td><td align="left" rowspan="1" colspan="1">Ethiopia</td><td align="left" rowspan="1" colspan="1">Human</td></tr><tr><td align="left" rowspan="1" colspan="1">Head louse, <italic>Pediculus capitis</italic> (B2560H)</td><td align="left" rowspan="1" colspan="1">Ethiopia</td><td align="left" rowspan="1" colspan="1">Human</td></tr><tr><td align="left" rowspan="1" colspan="1">Body louse, <italic>Pediculus humanus</italic> (B2563B2)</td><td align="left" rowspan="1" colspan="1">Ethiopia</td><td align="left" rowspan="1" colspan="1">Human</td></tr><tr><td align="left" rowspan="1" colspan="1">Head louse, <italic>Pediculus capitis</italic> (B2563H)</td><td align="left" rowspan="1" colspan="1">Ethiopia</td><td align="left" rowspan="1" colspan="1">Human</td></tr><tr><td align="left" rowspan="1" colspan="1">Hard tick, <italic>Haemaphysalis formosensis</italic> (B2340)</td><td align="left" rowspan="1" colspan="1">Japan</td><td align="left" rowspan="1" colspan="1">Unknown</td></tr><tr><td align="left" rowspan="1" colspan="1">Hard tick, <italic>Haemaphysalis parva</italic> (B2345)</td><td align="left" rowspan="1" colspan="1">Romania</td><td align="left" rowspan="1" colspan="1">Unknown</td></tr><tr><td align="left" rowspan="1" colspan="1">Hard tick, <italic>Rhipicephalus microplus</italic> (B2487)</td><td align="left" rowspan="1" colspan="1">Cambodia</td><td align="left" rowspan="1" colspan="1">Cattle</td></tr><tr><td align="left" rowspan="1" colspan="1">Hard tick, <italic>Amblyomma cajennense</italic> (B2479)</td><td align="left" rowspan="1" colspan="1">Brazil</td><td align="left" rowspan="1" colspan="1">Monkey</td></tr><tr><td align="left" rowspan="1" colspan="1">Hard tick, <italic>Rhipicephalus geigyi</italic> (B2538)</td><td align="left" rowspan="1" colspan="1">Burkina Faso</td><td align="left" rowspan="1" colspan="1">Unknown</td></tr><tr><td align="left" rowspan="1" colspan="1">Soft tick, <italic>Argas</italic> sp. (B2384)</td><td align="left" rowspan="1" colspan="1">Springbok</td><td align="left" rowspan="1" colspan="1">Bird</td></tr><tr><td align="left" rowspan="1" colspan="1">Soft tick, <italic>Otobius megnini</italic> (B2512)</td><td align="left" rowspan="1" colspan="1">Madagascar</td><td align="left" rowspan="1" colspan="1">Horse</td></tr></tbody></table></alternatives></table-wrap><p>For the tick species, short regions of <italic>cox1, cytb</italic> and <italic>rrnS</italic> genes were amplified first and sequenced using Sanger method with arthropod-conserved primers <xref ref-type="bibr" rid="pone.0073329-Kambhampati1">[30]</xref>, <xref ref-type="bibr" rid="pone.0073329-Simon1">[31]</xref>. Specific primers were then designed for each species from these regions, and used in most species with conserved tick primers in <italic>cytb</italic> and <italic>rrnS</italic>, designed from known sequences (<xref ref-type="supplementary-material" rid="pone.0073329.s001">Table S1</xref>). Entire mitochondrial genomes were then amplified in two overlapping fragments for the soft ticks: one fragment from <italic>cytb</italic> to <italic>cox1</italic>, and the other from <italic>cox1</italic> to <italic>cytb</italic>. The mt genomes of hard ticks were amplified in three overlapping fragments that span from <italic>cytb</italic> to <italic>c</italic>o<italic>x1</italic>, from <italic>cox1</italic> to <italic>rrnS</italic>, and from <italic>rrnS</italic> to <italic>cytb</italic>. The long amplicons from the soft and hard ticks were subject to Illumina sequencing. TaKaRa Ex Taq DNA polymerase kit was used to amplify the short <italic>cox1</italic>, <italic>cytb</italic> and <italic>rrnS</italic> fragments (<1 kb). Expand Long Range dNTPack kit (Roche) was used to amplify the long fragments (>1 kb). PCR conditions were optimized for each reaction, with the annealing temperature adjusted to suit the primers used, and the extension time set to one minute per kb of expected amplicon. General PCR conditions for Ex Taq were: 94°C for 60 seconds, followed by 40 cycles of 98°C for 10 seconds, 60°C for 30 seconds, 72°C for 1 minute and a final extension of 72°C for 2 minutes. General PCR conditions for Expand dNTPack were: 92°C for 2 minutes, 10 cycles of 92°C for 10 seconds, 55°C for 15 seconds, 60°C for 8 minutes, followed by 25 cycles of 92°C for 10 seconds 55°C for 15 seconds, 60°C for 8 minutes (increasing by 20 seconds per cycle), and a final extension of 68°C for 7 minutes. PCR products were examined on 1% agarose gel stained with ethidium bromide. DNA Molecular Weight Marker VII (Roche Diagnostics) and Low DNA Mass Ladder (Invitrogen) were used, respectively, to estimate the length and concentration of PCR products. Wizard SV Gel and PCR Clean-up System (Promega) was used to purify PCR products for sequencing. Sanger sequencing was with the ABI Prism BigDye v3.1 Terminator kit (Applied Biosystems) and the Applied Biosystems 3730<italic>xl</italic> DNA Analyzer at the Australian Genome Research Facility (AGRF). Illumina sequencing was with the HiSeq 2000 platform at the BGI, as stated above for the human lice. The sequence reads from the mitochondrial genomes of the human lice and ticks have been deposited in NCBI Sequence Read Archive (SRA) under the accession number SRR616254 and SRR616964.</p></sec><sec id="s2c"><title>Sequence and statistical analyses</title><p>Raw Illumina sequence-reads were processed before assembly. The adapter sequences were removed from each sequence read; the reads that contain more than 50% low quality bases (quality value < = 5) were excluded from further analysis. The Q20 bases rate of each pool is above 99.26% (Q20 refers to less than 1% error rate in the sequencing base reads). Processed sequence reads were assembled into contigs with Geneious (Biomatters), using the mitochondrial genome sequences of the human lice and ticks available as references <xref ref-type="bibr" rid="pone.0073329-Shao1">[22]</xref>, <xref ref-type="bibr" rid="pone.0073329-Shao2">[23]</xref>, <xref ref-type="bibr" rid="pone.0073329-Burger1">[29]</xref>. The assembly parameters were: minimum overlap length 80 bp, and minimum identity 95%. Consensus sequences were made for each contig with the “Generate consensus sequence” function in Geneious with the settings: 1) Threshold: Highest Quality – use chromatogram quality to call the best base; and 2) Assign Quality: Total – Sum quality of contributing bases minus non-contributing bases. Protein-coding genes and rRNA genes were identified with BLAST searches of GenBank <xref ref-type="bibr" rid="pone.0073329-Altschul1">[32]</xref>, <xref ref-type="bibr" rid="pone.0073329-McGinnisSMadden1">[33]</xref>. tRNA genes were identified with tRNA-Scan <xref ref-type="bibr" rid="pone.0073329-LoweTM1">[34]</xref> and ARWEN <xref ref-type="bibr" rid="pone.0073329-LaslettDCanback1">[35]</xref>. Heteroplasmy was measured with the “Find SNPs/variations” function in Geneious under the threshold “minimum variant frequency” >1.5%, and “find SNPs only”.</p><p>Statistic analyses were conducted to estimate the false positive rate and the false discovery rate for the heteroplasmy called above, according to Storey and Tibshirani (2003) <xref ref-type="bibr" rid="pone.0073329-Storey1">[36]</xref>. The false positive rate, indicated as p-value in Storey and Tibshirani (2003), was measured as the probability of a heteroplasmy to be false positive (P<sub>f</sub>) under a binomial distribution given the sequence-read coverage of that heteroplasmic site and the sequencing error of the Illumina platform. P<sub>f</sub> was calculated with the formula P<sub>f</sub> = C (n,k) * p∧k * (1-p)∧(n-k), where n is the sequence-read coverage of a heteroplasmy called at a given site of a contig; k is the times of that heteroplasmy seen at that site; C(n, k) is n!/(k! * (n-k)!); and p is the sequencing error, 0.5%, of the Illumina platform. The false discovery rate, indicated as q-value in Storey and Tibshirani (2003), was measured as the probability of a statistically significant heteroplasmy to be false positive (Q<sub>f</sub>). Q<sub>f</sub> was calculated from the P<sub>f</sub> values obtained above with the QValue program <xref ref-type="bibr" rid="pone.0073329-Storey1">[36]</xref>, <xref ref-type="bibr" rid="pone.0073329-Storey2">[37]</xref>. The statistical significance threshold for the heteroplasmy called was set at P<sub>f</sub><1% and Q<sub>f</sub><0.1%. Measurement of heteroplasmy was further corrected by excluding those at the sequencing-error hotspots revealed in Nakamura <italic>et al</italic>. (2011) and Minoche <italic>et al</italic>. (2011): 1) after four or more identical base calls; and 2) two adjacent SNP variants <xref ref-type="bibr" rid="pone.0073329-Nakamura1">[38]</xref>, <xref ref-type="bibr" rid="pone.0073329-Minoche1">[39]</xref>. One sample T test and Independent samples T test (P<sub>t</sub>) were used to compare the mean heteroplasmy levels between any two samples or genes.</p></sec></sec><sec id="s3"><title>Results</title><sec id="s3a"><title>Heteroplasmy in the human lice</title><p>Twenty contigs were obtained for each of the six human body lice and six human head lice, corresponding to the 20 mt minichromosomes of these lice. The gene content and gene arrangement in each minichromosome of these lice are the same as that reported in Shao et al. (2009, 2012) <xref ref-type="bibr" rid="pone.0073329-Shao1">[22]</xref>, <xref ref-type="bibr" rid="pone.0073329-Shao2">[23]</xref>. At the “minimum variant frequency” >1.5%, P<sub>f</sub><1% and Q<sub>f</sub><0.1%, a total of 494 heteroplasmic sites (excluding those at sequencing-error hotspots here and throughout this manuscript) were found in the mt coding regions (15.4 kb) of six human body lice and six human head lice, i.e. an average of 41.1±20.3 (mean ± standard deviation, n = 12) heteroplasmic sites per louse (<xref ref-type="table" rid="pone-0073329-t002">Table 2</xref>). These heteroplasmic sites have an average coverage of ∼6,000x. Of the 494 heteroplasmic sites, 374 sites are in the protein-coding genes, 21 heteroplasmic sites in the rRNA genes, and 99 heteroplasmic sites in the tRNA genes (<xref ref-type="table" rid="pone-0073329-t002">Table 2</xref>, <xref ref-type="supplementary-material" rid="pone.0073329.s002">Table S2</xref>). Of the 374 heteroplasmic sites in the protein-coding genes, 252 sites are synonymous, i.e. resulting in no amino acid changes, whereas the other 122 sites are nonsynonymous, i.e. resulting in amino acid changes (<xref ref-type="table" rid="pone-0073329-t002">Table 2</xref>, <xref ref-type="supplementary-material" rid="pone.0073329.s006">Table S6</xref>). There is no significant variation in the level of heteroplasmy among different mt genes in the human lice with the exception of <italic>cox1</italic> gene, which has significantly higher level of heteroplasmy (13.3±15.7 sites per kb, n = 12) than other genes (P<sub>t</sub><0.05) (<xref ref-type="fig" rid="pone-0073329-g001">Figure 1</xref>; <xref ref-type="table" rid="pone-0073329-t003">Table 3</xref>). Each of the 12 human lice share 15 to 63 heteroplasmic sites with one or more other lice (<xref ref-type="fig" rid="pone-0073329-g002">Figure 2</xref>; <xref ref-type="table" rid="pone-0073329-t004">Table 4</xref>; <xref ref-type="supplementary-material" rid="pone.0073329.s003">Table S3</xref>), indicating part of the heteroplasmic sites are inheritable and have limited conservation within species.</p><fig id="pone-0073329-g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0073329.g001</object-id><label>Figure 1</label><caption><title>Level of heteroplasmy in the mitochondrial genes of human lice and ticks.</title><p>Error bars stand for the standard error of the mean value. <italic>Amblyomma cajennense</italic> was excluded due to its extremely high level of heteroplasmy relative to other ticks.</p></caption><graphic xlink:href="pone.0073329.g001"></graphic></fig><fig id="pone-0073329-g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0073329.g002</object-id><label>Figure 2</label><caption><title>Numbers of shared and unshared heteroplasmic sites in human lice.</title><p>Error bars stand for the standard error of the mean value.</p></caption><graphic xlink:href="pone.0073329.g002"></graphic></fig><table-wrap id="pone-0073329-t002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0073329.t002</object-id><label>Table 2</label><caption><title>Heteroplasmy in the mitochondrial coding regions of human lice and ticks.</title></caption><alternatives><graphic id="pone-0073329-t002-2" xlink:href="pone.0073329.t002"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Number of heteroplasmic sites</td><td align="left" rowspan="1" colspan="1">Total of 12 lice</td><td align="left" rowspan="1" colspan="1">Average per louse</td><td align="left" rowspan="1" colspan="1">Total of seven ticks</td><td align="left" rowspan="1" colspan="1">Average per tick</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">All mitochondrial genes</td><td align="left" rowspan="1" colspan="1">494</td><td align="left" rowspan="1" colspan="1">41.1</td><td align="left" rowspan="1" colspan="1">261</td><td align="left" rowspan="1" colspan="1">37.3</td></tr><tr><td align="left" rowspan="1" colspan="1">Protein-coding genes</td><td align="left" rowspan="1" colspan="1">374</td><td align="left" rowspan="1" colspan="1">31.2</td><td align="left" rowspan="1" colspan="1">221</td><td align="left" rowspan="1" colspan="1">31.6</td></tr><tr><td align="left" rowspan="1" colspan="1">Synonymous</td><td align="left" rowspan="1" colspan="1">252</td><td align="left" rowspan="1" colspan="1">21.0</td><td align="left" rowspan="1" colspan="1">158</td><td align="left" rowspan="1" colspan="1">22.6</td></tr><tr><td align="left" rowspan="1" colspan="1">Codon 1<sup>st</sup></td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1">1.2</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">1.4</td></tr><tr><td align="left" rowspan="1" colspan="1">Codon 2<sup>nd</sup></td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">0.5</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">Codon 3<sup>rd</sup></td><td align="left" rowspan="1" colspan="1">232</td><td align="left" rowspan="1" colspan="1">19.3</td><td align="left" rowspan="1" colspan="1">148</td><td align="left" rowspan="1" colspan="1">21.2</td></tr><tr><td align="left" rowspan="1" colspan="1">Nonsynonymous</td><td align="left" rowspan="1" colspan="1">122*</td><td align="left" rowspan="1" colspan="1">10.1</td><td align="left" rowspan="1" colspan="1">63</td><td align="left" rowspan="1" colspan="1">9.0</td></tr><tr><td align="left" rowspan="1" colspan="1">Codon 1<sup>st</sup></td><td align="left" rowspan="1" colspan="1">60</td><td align="left" rowspan="1" colspan="1">5.0</td><td align="left" rowspan="1" colspan="1">26</td><td align="left" rowspan="1" colspan="1">3.7</td></tr><tr><td align="left" rowspan="1" colspan="1">Codon 2<sup>nd</sup></td><td align="left" rowspan="1" colspan="1">54</td><td align="left" rowspan="1" colspan="1">4.5</td><td align="left" rowspan="1" colspan="1">24</td><td align="left" rowspan="1" colspan="1">3.4</td></tr><tr><td align="left" rowspan="1" colspan="1">Codon 3<sup>rd</sup></td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">0.6</td><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1">1.9</td></tr><tr><td align="left" rowspan="1" colspan="1">rRNA genes</td><td align="left" rowspan="1" colspan="1">21</td><td align="left" rowspan="1" colspan="1">1.8</td><td align="left" rowspan="1" colspan="1">24</td><td align="left" rowspan="1" colspan="1">3.4</td></tr><tr><td align="left" rowspan="1" colspan="1">tRNA genes</td><td align="left" rowspan="1" colspan="1">99</td><td align="left" rowspan="1" colspan="1">8.3</td><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1">2.3</td></tr><tr><td align="left" rowspan="1" colspan="1">Anti-codon 3<sup>rd</sup></td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">0.3</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt101"><label></label><p>Note: Heteroplasmy was called at minimum variant frequency >1.5%, false positive rate P<sub>f</sub><1% and false discovery rate Q<sub>f</sub><0.1%, excluding those at sequencing error hotpots. There are 13 protein coding genes (<italic>atp8-atp6</italic>, <italic>cytb</italic>, <italic>cox1-3</italic>, <italic>nad1-6</italic> and <italic>4L</italic>), two rRNA genes (<italic>rrnL</italic> and r<italic>rnS</italic>), 22 tRNA genes. *46 of the 122 sites have amino acid hydrophobic or hydrophilic changes; 25 of the 122 sites have amino acid polar or non-polar changes.</p></fn></table-wrap-foot></table-wrap><table-wrap id="pone-0073329-t003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0073329.t003</object-id><label>Table 3</label><caption><title>Number of heteroplasmic sites in the mitochondrial genes of the human lice.</title></caption><alternatives><graphic id="pone-0073329-t003-3" xlink:href="pone.0073329.t003"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Human lice</td><td align="left" rowspan="1" colspan="1"><italic>atp8-atp6</italic></td><td align="left" rowspan="1" colspan="1"><italic>cytb</italic></td><td align="left" rowspan="1" colspan="1"><italic>cox1</italic></td><td align="left" rowspan="1" colspan="1"><italic>cox3</italic></td><td align="left" rowspan="1" colspan="1"><italic>nad1</italic></td><td align="left" rowspan="1" colspan="1"><italic>nad5</italic></td><td align="left" rowspan="1" colspan="1"><italic>nad6</italic></td><td align="left" rowspan="1" colspan="1"><italic>nad4L</italic></td><td align="left" rowspan="1" colspan="1"><italic>nad4</italic></td><td align="left" rowspan="1" colspan="1"><italic>nad2</italic></td><td align="left" rowspan="1" colspan="1"><italic>nad3</italic></td><td align="left" rowspan="1" colspan="1"><italic>cox2</italic></td><td align="left" rowspan="1" colspan="1"><italic>rrnL</italic></td><td align="left" rowspan="1" colspan="1"><italic>rrnS</italic></td><td align="left" rowspan="1" colspan="1"><italic>tRNAs</italic></td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">B2470B1</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">3(2)</td><td align="left" rowspan="1" colspan="1">15(9)</td><td align="left" rowspan="1" colspan="1">4(3)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">2(2)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">3(2)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">1(0)</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">4</td></tr><tr><td align="left" rowspan="1" colspan="1">B2471H5</td><td align="left" rowspan="1" colspan="1">2(2)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">8(2)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">7(4)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">2(0)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">10</td></tr><tr><td align="left" rowspan="1" colspan="1">B2557B</td><td align="left" rowspan="1" colspan="1">2(2)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">2(2)</td><td align="left" rowspan="1" colspan="1">3(2)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">5(3)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">1(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">1(0)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">14</td></tr><tr><td align="left" rowspan="1" colspan="1">B2557H</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">2(1)</td><td align="left" rowspan="1" colspan="1">3(2)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">5(3)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">9</td></tr><tr><td align="left" rowspan="1" colspan="1">B2558B</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">3(2)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">2(1)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">9</td></tr><tr><td align="left" rowspan="1" colspan="1">B2558H</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">3(3)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">5(2)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">10</td></tr><tr><td align="left" rowspan="1" colspan="1">B2516B2</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">52(4)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">3(3)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">6</td></tr><tr><td align="left" rowspan="1" colspan="1">B2516H2</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">45(1)</td><td align="left" rowspan="1" colspan="1">2(2)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">3</td></tr><tr><td align="left" rowspan="1" colspan="1">B2560B</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">4(1)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">3(3)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">6</td></tr><tr><td align="left" rowspan="1" colspan="1">B2560H</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">54(9)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">2(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">8</td></tr><tr><td align="left" rowspan="1" colspan="1">B2563B2</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">61(8)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">4(3)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">1(0)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">6</td></tr><tr><td align="left" rowspan="1" colspan="1">B2563H</td><td align="left" rowspan="1" colspan="1">7(3)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">5(2)</td><td align="left" rowspan="1" colspan="1">6(4)</td><td align="left" rowspan="1" colspan="1">6(3)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">8(4)</td><td align="left" rowspan="1" colspan="1">4(2)</td><td align="left" rowspan="1" colspan="1">1(0)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">14</td></tr><tr><td align="left" rowspan="1" colspan="1">Average</td><td align="left" rowspan="1" colspan="1">1.3</td><td align="left" rowspan="1" colspan="1">0.3</td><td align="left" rowspan="1" colspan="1">20.8</td><td align="left" rowspan="1" colspan="1">2.0</td><td align="left" rowspan="1" colspan="1">0.6</td><td align="left" rowspan="1" colspan="1">3.0</td><td align="left" rowspan="1" colspan="1">0.5</td><td align="left" rowspan="1" colspan="1">0.1</td><td align="left" rowspan="1" colspan="1">1.2</td><td align="left" rowspan="1" colspan="1">0.3</td><td align="left" rowspan="1" colspan="1">0.6</td><td align="left" rowspan="1" colspan="1">0.6</td><td align="left" rowspan="1" colspan="1">1.1</td><td align="left" rowspan="1" colspan="1">0.7</td><td align="left" rowspan="1" colspan="1">8.3</td></tr><tr><td align="left" rowspan="1" colspan="1">A/L</td><td align="left" rowspan="1" colspan="1">1.54</td><td align="left" rowspan="1" colspan="1">0.28</td><td align="left" rowspan="1" colspan="1">13.36</td><td align="left" rowspan="1" colspan="1">2.47</td><td align="left" rowspan="1" colspan="1">0.71</td><td align="left" rowspan="1" colspan="1">1.81</td><td align="left" rowspan="1" colspan="1">1.00</td><td align="left" rowspan="1" colspan="1">0.39</td><td align="left" rowspan="1" colspan="1">0.91</td><td align="left" rowspan="1" colspan="1">0.31</td><td align="left" rowspan="1" colspan="1">1.65</td><td align="left" rowspan="1" colspan="1">0.84</td><td align="left" rowspan="1" colspan="1">1.05</td><td align="left" rowspan="1" colspan="1">0.91</td><td align="left" rowspan="1" colspan="1">5.41</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt102"><label></label><p>Note: Heteroplasmy was called at minimum variant frequency >1.5%, false positive rate P<sub>f</sub><1% and false discovery rate Q<sub>f</sub><0.1%, excluding those at sequencing error hotpots. “H” is for human head lice and “B” is for human body lice. Numbers in bracket refers to the number of nonsynonymous sites of heteroplasmy that result in amino acid changes. A/L stands for Average/Length (kb).</p></fn></table-wrap-foot></table-wrap><table-wrap id="pone-0073329-t004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0073329.t004</object-id><label>Table 4</label><caption><title>Number of heteroplasmic sites in the mitochondrial coding regions shared between human body lice and head lice.</title></caption><alternatives><graphic id="pone-0073329-t004-4" xlink:href="pone.0073329.t004"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">B2471H5</td><td align="left" rowspan="1" colspan="1">B2557B</td><td align="left" rowspan="1" colspan="1">B2557H</td><td align="left" rowspan="1" colspan="1">B2558B</td><td align="left" rowspan="1" colspan="1">B2558H</td><td align="left" rowspan="1" colspan="1">B2516B2</td><td align="left" rowspan="1" colspan="1">B2516H2</td><td align="left" rowspan="1" colspan="1">B2560B</td><td align="left" rowspan="1" colspan="1">B2560H</td><td align="left" rowspan="1" colspan="1">B2563B2</td><td align="left" rowspan="1" colspan="1">B2563H</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">B2470B1</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">6</td></tr><tr><td align="left" rowspan="1" colspan="1">B2471H5</td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">6</td></tr><tr><td align="left" rowspan="1" colspan="1">B2557B</td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">11</td></tr><tr><td align="left" rowspan="1" colspan="1">B2557H</td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">7</td></tr><tr><td align="left" rowspan="1" colspan="1">B2558B</td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">10</td></tr><tr><td align="left" rowspan="1" colspan="1">B2558H</td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">9</td></tr><tr><td align="left" rowspan="1" colspan="1">B2516B2</td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1">55</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">B2516H2</td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">21</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">1</td></tr><tr><td align="left" rowspan="1" colspan="1">B2560B</td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">2</td></tr><tr><td align="left" rowspan="1" colspan="1">B2560H</td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">1</td></tr><tr><td align="left" rowspan="1" colspan="1">B2563B2</td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">2</td></tr></tbody></table></alternatives></table-wrap><p>The three body lice and three head lice from Ethiopia have an average of 62.2±10.6 heteroplasmic sites per louse (n = 6), which is significantly higher than that of the lice from China (an average of 33.5±2.1 sites per louse, n = 2) (P<sub>t</sub><0.05) and France (an average of 24.8±5.9 sites per louse, n = 4) (P<sub>t</sub><0.05) (<xref ref-type="table" rid="pone-0073329-t005">Table 5</xref>). Four heteroplasmic sites in four tRNA genes, <italic>trnN</italic>, <italic>trnD</italic>, <italic>trnT</italic> and <italic>trnM</italic>, are present only in the lice from Ethiopia (n = 6) but not in those from China (n = 2) and France (n = 4) (<xref ref-type="supplementary-material" rid="pone.0073329.s002">Table S2</xref>). There is no significant difference in the heteroplasmy level among the lice from Ethiopia, nor among those from China and France (P<sub>t</sub>>0.05). There is no significant difference either between the lice from China and those from France in the level of heteroplasmy (P<sub>t</sub>>0.05).</p><table-wrap id="pone-0073329-t005" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0073329.t005</object-id><label>Table 5</label><caption><title>Heteroplasmy in the mitochondrial coding regions of human lice from China, France and Ethiopia.</title></caption><alternatives><graphic id="pone-0073329-t005-5" xlink:href="pone.0073329.t005"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Louse sample#</td><td align="left" rowspan="1" colspan="1">Country</td><td align="left" rowspan="1" colspan="1">Total heteroplasmic sites</td><td align="left" rowspan="1" colspan="1">Protein-coding genes</td><td align="left" rowspan="1" colspan="1">rRNA genes</td><td align="left" rowspan="1" colspan="1">tRNA genes</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">B2470B1</td><td align="left" rowspan="1" colspan="1">China</td><td align="left" rowspan="1" colspan="1">35</td><td align="left" rowspan="1" colspan="1">30 (20)</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">4</td></tr><tr><td align="left" rowspan="1" colspan="1">B2471H5</td><td align="left" rowspan="1" colspan="1">China</td><td align="left" rowspan="1" colspan="1">32</td><td align="left" rowspan="1" colspan="1">22 (11)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">10</td></tr><tr><td align="left" rowspan="1" colspan="1">B2557B</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">32</td><td align="left" rowspan="1" colspan="1">15 (10)</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">14</td></tr><tr><td align="left" rowspan="1" colspan="1">B2557H</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">27</td><td align="left" rowspan="1" colspan="1">12 (8)</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">9</td></tr><tr><td align="left" rowspan="1" colspan="1">B2558B</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">21</td><td align="left" rowspan="1" colspan="1">8 (6)</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">9</td></tr><tr><td align="left" rowspan="1" colspan="1">B2558H</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">19</td><td align="left" rowspan="1" colspan="1">9 (6)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">10</td></tr><tr><td align="left" rowspan="1" colspan="1">B2516B2</td><td align="left" rowspan="1" colspan="1">Ethiopia</td><td align="left" rowspan="1" colspan="1">67</td><td align="left" rowspan="1" colspan="1">58 (10)</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">6</td></tr><tr><td align="left" rowspan="1" colspan="1">B2516H2</td><td align="left" rowspan="1" colspan="1">Ethiopia</td><td align="left" rowspan="1" colspan="1">50</td><td align="left" rowspan="1" colspan="1">47 (3)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">3</td></tr><tr><td align="left" rowspan="1" colspan="1">B2560B</td><td align="left" rowspan="1" colspan="1">Ethiopia</td><td align="left" rowspan="1" colspan="1">17</td><td align="left" rowspan="1" colspan="1">9 (6)</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">6</td></tr><tr><td align="left" rowspan="1" colspan="1">B2560H</td><td align="left" rowspan="1" colspan="1">Ethiopia</td><td align="left" rowspan="1" colspan="1">65</td><td align="left" rowspan="1" colspan="1">56 (9)</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">8</td></tr><tr><td align="left" rowspan="1" colspan="1">B2563B2</td><td align="left" rowspan="1" colspan="1">Ethiopia</td><td align="left" rowspan="1" colspan="1">76</td><td align="left" rowspan="1" colspan="1">69 (13)</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">6</td></tr><tr><td align="left" rowspan="1" colspan="1">B2563H</td><td align="left" rowspan="1" colspan="1">Ethiopia</td><td align="left" rowspan="1" colspan="1">53</td><td align="left" rowspan="1" colspan="1">39 (20)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">14</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt103"><label></label><p>Note: Heteroplasmy was called at minimum variant frequency >1.5%, false positive rate P<sub>f</sub><1% and false discovery rate Q<sub>f</sub><0.1%, excluding those at sequencing error hotpots. “H” is for human head lice and “B” is for human body lice. Numbers in bracket refers to the number of nonsynonymous sites of heteroplasmy that result in amino acid changes.</p></fn></table-wrap-foot></table-wrap></sec><sec id="s3b"><title>Heteroplasmy in the ticks</title><p>A single contig was obtained for each of the seven species of ticks, corresponding to the mt chromosomes of these ticks. The mt genomic organization of these ticks is the same as that reported for hard ticks and soft ticks in a number of previous studies <xref ref-type="bibr" rid="pone.0073329-Shao4">[27]</xref>, <xref ref-type="bibr" rid="pone.0073329-Shao5">[28]</xref>, <xref ref-type="bibr" rid="pone.0073329-Burger1">[29]</xref>. At the “minimum variant frequency” >1.5%, P<sub>f</sub><1% and Q<sub>f</sub><0.1%, a total of 261 heteroplasmic sites were found in the seven tick samples, an average of 37.3±22.0 sites per tick (n = 7) (<xref ref-type="table" rid="pone-0073329-t002">Table 2</xref>). These heteroplasmic sites have an average coverage of ∼1,700x. Of the 261 heteroplasmic sites, 221 sites are in the protein-coding genes, 24 sites are in the rRNA genes, and the other 16 sites are in the tRNA genes (<xref ref-type="table" rid="pone-0073329-t002">Table 2</xref>; <xref ref-type="supplementary-material" rid="pone.0073329.s004">Table S4</xref>). There is no significant variation in the level of heteroplasmy among different mt genes in the ticks (P<sub>t</sub>>0.05) (<xref ref-type="fig" rid="pone-0073329-g001">Figure 1</xref>; <xref ref-type="table" rid="pone-0073329-t006">Table 6</xref>). Of the 221 heteroplasmic sites in the protein-coding genes, 63 sites are synonymous whereas the rest are nonsynonymous. No heteroplasmic sites were shared between any two species of ticks except one in <italic>nad3</italic> shared between <italic>Amblyomma cajennense</italic> and <italic>Otobius megnini</italic>, indicating nearly no conservation of heteroplasmy between species. The level of heteroplasmy varies among different species of ticks (<xref ref-type="table" rid="pone-0073329-t007">Table 7</xref>). <italic>Amblyomma cajennense</italic> has 166 heteroplasmic sites (11.5 sites per kb, n = 1) and is significantly higher in heteroplasmy level than the other six ticks (1.1±0.9 sites per tick per kb, n = 6) (P<sub>t</sub><0.01), and even higher than the human lice (2.9±1.4 sites per louse per kb, n = 12; P<sub>t</sub><0.01).</p><table-wrap id="pone-0073329-t006" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0073329.t006</object-id><label>Table 6</label><caption><title>Number of heteropalsmic sites in the mitochondrial genes of ticks.</title></caption><alternatives><graphic id="pone-0073329-t006-6" xlink:href="pone.0073329.t006"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Ticks</td><td align="left" rowspan="1" colspan="1"><italic>atp8-atp6</italic></td><td align="left" rowspan="1" colspan="1"><italic>cytb</italic></td><td align="left" rowspan="1" colspan="1"><italic>cox1</italic></td><td align="left" rowspan="1" colspan="1"><italic>cox2</italic></td><td align="left" rowspan="1" colspan="1"><italic>cox3</italic></td><td align="left" rowspan="1" colspan="1"><italic>nad1</italic></td><td align="left" rowspan="1" colspan="1"><italic>nad2</italic></td><td align="left" rowspan="1" colspan="1"><italic>nad3</italic></td><td align="left" rowspan="1" colspan="1"><italic>nad4</italic></td><td align="left" rowspan="1" colspan="1"><italic>nad4L</italic></td><td align="left" rowspan="1" colspan="1"><italic>nad5</italic></td><td align="left" rowspan="1" colspan="1"><italic>nad6</italic></td><td align="left" rowspan="1" colspan="1"><italic>rrnL</italic></td><td align="left" rowspan="1" colspan="1"><italic>rrnS</italic></td><td align="left" rowspan="1" colspan="1"><italic>tRNAs</italic></td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1"><italic>Haemaphysalis formosensis</italic></td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">5(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>Haemaphysalis parva</italic></td><td align="left" rowspan="1" colspan="1">2(1)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">2(2)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>Rhipicephalus microplus</italic></td><td align="left" rowspan="1" colspan="1">2(1)</td><td align="left" rowspan="1" colspan="1">3(0)</td><td align="left" rowspan="1" colspan="1">6(0)</td><td align="left" rowspan="1" colspan="1">2(0)</td><td align="left" rowspan="1" colspan="1">3(0)</td><td align="left" rowspan="1" colspan="1">3(0)</td><td align="left" rowspan="1" colspan="1">2(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">6(2)</td><td align="left" rowspan="1" colspan="1">1(0)</td><td align="left" rowspan="1" colspan="1">6(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">6</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>Amblyomma cajennense</italic></td><td align="left" rowspan="1" colspan="1">13(4)</td><td align="left" rowspan="1" colspan="1">11(10)</td><td align="left" rowspan="1" colspan="1">22(0)</td><td align="left" rowspan="1" colspan="1">8(0)</td><td align="left" rowspan="1" colspan="1">7(3)</td><td align="left" rowspan="1" colspan="1">13(1)</td><td align="left" rowspan="1" colspan="1">6(3)</td><td align="left" rowspan="1" colspan="1">6(1)</td><td align="left" rowspan="1" colspan="1">19(1)</td><td align="left" rowspan="1" colspan="1">3(0)</td><td align="left" rowspan="1" colspan="1">27(10)</td><td align="left" rowspan="1" colspan="1">1(0)</td><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">8</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>Argas</italic> sp.</td><td align="left" rowspan="1" colspan="1">1(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">2(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">1(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">3(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">3(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>Rhipicephalus geigyi</italic></td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">1(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">5(4)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>Otobius megnini</italic></td><td align="left" rowspan="1" colspan="1">5(5)</td><td align="left" rowspan="1" colspan="1">1(1)</td><td align="left" rowspan="1" colspan="1">6(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">1(0)</td><td align="left" rowspan="1" colspan="1">3(3)</td><td align="left" rowspan="1" colspan="1">3(3)</td><td align="left" rowspan="1" colspan="1">5(5)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0(0)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">Average</td><td align="left" rowspan="1" colspan="1">3.4</td><td align="left" rowspan="1" colspan="1">3.0</td><td align="left" rowspan="1" colspan="1">5.3</td><td align="left" rowspan="1" colspan="1">1.4</td><td align="left" rowspan="1" colspan="1">1.7</td><td align="left" rowspan="1" colspan="1">3.0</td><td align="left" rowspan="1" colspan="1">1.6</td><td align="left" rowspan="1" colspan="1">1.6</td><td align="left" rowspan="1" colspan="1">4.7</td><td align="left" rowspan="1" colspan="1">0.6</td><td align="left" rowspan="1" colspan="1">5.1</td><td align="left" rowspan="1" colspan="1">0.1</td><td align="left" rowspan="1" colspan="1">2.1</td><td align="left" rowspan="1" colspan="1">1.3</td><td align="left" rowspan="1" colspan="1">2.3</td></tr><tr><td align="left" rowspan="1" colspan="1">A/L</td><td align="left" rowspan="1" colspan="1">4.16</td><td align="left" rowspan="1" colspan="1">2.72</td><td align="left" rowspan="1" colspan="1">3.44</td><td align="left" rowspan="1" colspan="1">2.08</td><td align="left" rowspan="1" colspan="1">2.19</td><td align="left" rowspan="1" colspan="1">3.03</td><td align="left" rowspan="1" colspan="1">1.68</td><td align="left" rowspan="1" colspan="1">4.75</td><td align="left" rowspan="1" colspan="1">3.56</td><td align="left" rowspan="1" colspan="1">2.15</td><td align="left" rowspan="1" colspan="1">3.08</td><td align="left" rowspan="1" colspan="1">0.23</td><td align="left" rowspan="1" colspan="1">2.19</td><td align="left" rowspan="1" colspan="1">1.08</td><td align="left" rowspan="1" colspan="1">1.67</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt104"><label></label><p>Notes: Heteroplasmy was called at minimum variant frequency >1.5%, false positive rate P<sub>f</sub><1% and false discovery rate Q<sub>f</sub><0.1%, excluding those at sequencing error hotpots. Numbers in bracket refers to nonsynonymous sites of heteroplasmy that cause amino acid changes.</p></fn></table-wrap-foot></table-wrap><table-wrap id="pone-0073329-t007" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0073329.t007</object-id><label>Table 7</label><caption><title>Heteroplasmy in the mitochondrial coding regions of ticks.</title></caption><alternatives><graphic id="pone-0073329-t007-7" xlink:href="pone.0073329.t007"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Ticks</td><td align="left" rowspan="1" colspan="1">Total heteroplasmic sites</td><td align="left" rowspan="1" colspan="1">Protein-coding genes</td><td align="left" rowspan="1" colspan="1">rRNA genes</td><td align="left" rowspan="1" colspan="1">tRNA genes</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1"><italic>Haemaphysalis formosensis</italic></td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">6(1)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>Haemaphysalis parva</italic></td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">5(4)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>Rhipicephalus microplus</italic></td><td align="left" rowspan="1" colspan="1">40</td><td align="left" rowspan="1" colspan="1">34(3)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">6</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>Amblyomma cajennense</italic></td><td align="left" rowspan="1" colspan="1">166</td><td align="left" rowspan="1" colspan="1">136(33)</td><td align="left" rowspan="1" colspan="1">22</td><td align="left" rowspan="1" colspan="1">8</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>Argas</italic> sp.</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">10(0)</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>Rhipicephalus geigyi</italic></td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">7(5)</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>Otobius megnini</italic></td><td align="left" rowspan="1" colspan="1">24</td><td align="left" rowspan="1" colspan="1">23(17)</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">0</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt105"><label></label><p>Note: Heteroplasmy was called at minimum variant frequency >1.5%, false positive rate P<sub>f</sub><1% and false discovery rate Q<sub>f</sub><0.1%, excluding those at sequencing error hotpots. Numbers in bracket refers to the number of nonsynonymous sites of heteroplasmy that result in amino acid changes.</p></fn></table-wrap-foot></table-wrap></sec></sec><sec id="s4"><title>Discussion</title><p>Heteroplasmy, i.e. the presence of more than one type of mitochondrial DNA sequence within an individual, is common in humans, other bilateral animals and plants <xref ref-type="bibr" rid="pone.0073329-Kmiec1">[7]</xref>. Heteroplasmy occurs at a relatively low level in humans, with 32 heteroplasmic sites revealed in a human individual by high throughput sequencing at 16,700x coverage <xref ref-type="bibr" rid="pone.0073329-He1">[9]</xref>. Studies of other animals and plants with low-coverage sequencing also showed low levels of heteroplasmy <xref ref-type="bibr" rid="pone.0073329-Chinnery1">[12]</xref>, <xref ref-type="bibr" rid="pone.0073329-Kann1">[18]</xref>, <xref ref-type="bibr" rid="pone.0073329-Leeuwen1">[20]</xref>, <xref ref-type="bibr" rid="pone.0073329-Sachadyn1">[40]</xref>, <xref ref-type="bibr" rid="pone.0073329-Messenger1">[41]</xref>. Thus, low level of heteroplasmy has been assumed to be a general feature for eukaryotes.</p><p>Herd <italic>et al</italic>. (2012) identified, however, substantially more heteroplasmic sites in mt <italic>cox1</italic> gene in the human lice than in humans, with 13 heteroplasmic sites in this gene in human body lice and four heteroplasmic sites in human head lice, in comparison to one heteroplasmic site in this gene in humans revealed by high throughput sequencing <xref ref-type="bibr" rid="pone.0073329-He1">[9]</xref>. Herd <italic>et al</italic>. (2012) suggested that the high level of heteroplasmy in the human lice might be linked to the fragmentation of mt genomes in these lice, possibly through increased recombination activities between minichromosomes <xref ref-type="bibr" rid="pone.0073329-Shao1">[22]</xref>, <xref ref-type="bibr" rid="pone.0073329-Shao2">[23]</xref>. In the present study, we investigated further whether or not changes in mt genome organization is associated with the level of heteroplasmy by sequencing the entire mt coding region (13 protein-coding genes, two rRNA genes and 22 tRNA genes) of six human body lice, six head lice and seven species of ticks with an Illumina HiSeq platform.</p><p>To ensure the accuracy of the measurement of the heteroplasmy level from high throughput sequencing data, we set the minimum heteroplasmy-detection threshold to 1.5%, i.e. only the heteroplasmic sites that have a frequency above 1.5% are counted; sites with a frequency below 1.5% are discarded to exclude sequencing errors. As a general guideline, the sequencing-error rate of the Illumina HiSeq platform is ∼0.5% or less <xref ref-type="bibr" rid="pone.0073329-He1">[9]</xref>, <xref ref-type="bibr" rid="pone.0073329-Quinlan1">[42]</xref>, <xref ref-type="bibr" rid="pone.0073329-Sosa1">[43]</xref>. Our threshold-setting to 1.5% is consistent with those in He <italic>et al</italic>. (2010) and Doan <italic>et al</italic>. (2012), which set the thresholds to 1.6% and 1.5% respectively (i.e. triple the sequencing-error rate) <xref ref-type="bibr" rid="pone.0073329-He1">[9]</xref>, <xref ref-type="bibr" rid="pone.0073329-Doan1">[44]</xref>. Furthermore, we set the statistical significance threshold at P<sub>f</sub><1.0% and Q<sub>f</sub><0.1%, i.e. the probability of a heteroplasmy we called to be false positive is smaller than 1.0%, and the probability of a statistically significant heteroplasmy we called to be false positive is smaller than 0.1% (<xref ref-type="supplementary-material" rid="pone.0073329.s005">Table S5</xref>). We also excluded the heteroplasmy that fall into sequencing-error hotspots for Ilumina Hiseq platform revealed by Nakamura <italic>et al</italic>. (2011) and Minoche <italic>et al</italic>. (2011) <xref ref-type="bibr" rid="pone.0073329-Nakamura1">[38]</xref>, <xref ref-type="bibr" rid="pone.0073329-Minoche1">[39]</xref>. By all these measures, we were able to catch the full spectrum of the heteroplasmy while minimize the influence of high throughput sequencing errors.</p><p>We should point out that the parameters we used above to call heteroplasmy from high throughput sequence data, i.e. minimum detection threshold >1.5%, P<sub>f</sub><1.0% and Q<sub>f</sub><0.1%, are arbitrary. Indeed, both more relaxed and more stringent parameters can be used, which may increase or decrease the total number of heteroplasmic sites called from the data. The choice of these parameters, however, will not affect the conclusion for the comparison of the heteroplasmy levels between any two samples as long as the same set of parameters are applied. For instance, we also tested minimum detection threshold >1.0%, >3.0% and >5.0%. These three parameters changed substantially the cut-off sequence-read coverage at which heteroplasmy is called (<xref ref-type="fig" rid="pone-0073329-g003">Figure 3</xref>) but did not change the conclusion made when minimum detection threshold >1.5% was used.</p><fig id="pone-0073329-g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0073329.g003</object-id><label>Figure 3</label><caption><title>Plots of false positive rate (P<sub>f</sub>) against heteroplasmic-site coverage at minimum variant frequency 1.0%, 1.5%, 3.0% and 5.0% respectively.</title></caption><graphic xlink:href="pone.0073329.g003"></graphic></fig><p>Despite all having fragmented mt genomes, the human lice from Ethiopia have significantly higher level of heteroplasmy (62.2±10.6 per louse, n = 6) than those from China (33.5±2.1 sites per louse, n = 2; P<sub>t</sub><0.05) and France (24.8±5.9 sites per louse, n = 4; P<sub>t</sub><0.05), indicating substantial variation in heteroplasmy level among different populations of human lice. The same is true for the seven species of ticks. Although all of these ticks have the typical mt genomes of animals with all genes on a single chromosome, <italic>Amblyomma cajennense</italic> has significantly more heteroplasmic sites (166, n = 1) than the other six species of ticks (15.8±13.6 sites per tick, n = 1; P<sub>t</sub><0.05) (<xref ref-type="table" rid="pone-0073329-t007">Table 7</xref>). Fragmentation of mt genome in the human lice does not appear to have any bearing on the level of heteroplasmy as the two human lice from China have similar level of heteroplasmy (33.5±2.1 per louse) as humans, 32 heteroplasmic sites reported in He <italic>et al</italic>. (2010) <xref ref-type="bibr" rid="pone.0073329-He1">[9]</xref> and 34 sites in Li <italic>et al</italic>. (2012) <xref ref-type="bibr" rid="pone.0073329-Li1">[45]</xref>, whereas the four human lice from France have slightly lower level of heteroplasmy (24.8±5.9 sites per louse) than humans (P<sub>t</sub>>0.05). On the other hand, the tick, <italic>Amblyomma cajennense</italic>, has significantly higher level of heteroplasmy (166, n = 1) even than the human lice from Ethiopia (62.2±10.6 per louse, n = 6; P<sub>t</sub><0.05). With the exception of <italic>Amblyomma cajennense</italic>, the other six species of ticks have significant lower level of heteroplasmy (15.8±13.6 sites per tick, n = 1; P<sub>t</sub><0.05) than both the human lice and humans, reconfirming the view that mt genome organization has no direct influence on the level of heteroplasmy.</p><p>In conclusion, our results showed substantial variation in heteroplasmy level in mt genomes among human lice from different countries and among different species of ticks. Thus, we can reject the possible link between the elevated level of heteroplasmy and the fragmented mt genome organization in the human lice. Many questions, however, remain to be answered. Why does <italic>cox1</italic> have more heteroplasmy than other genes in human lice (but not in ticks)? Why do human lice from Ethiopia have higher levels of heteroplasmy than those from China and France? Why does <italic>Amblyomma cajennense</italic> have more heteroplasmy than other species of ticks? Are there any genomic, biochemical or life history factors that may contribute to the level of heteroplasmy? And more importantly, what does the varying level of heteroplasy mean to the fitness and evolution of individuals, populations and species? With the data generated in the current study, we cannot answer these questions yet. Both lice and ticks are ectoparasites of vertebrates. Parasitic lifestyle, however, does not appear to have any role in the variation of heteroplasmy level as there is clear difference in heteroplasmy level between lice and ticks and within each of them although all of them have the parasitic lifestyle. Recombination between mt minichromosomes observed in human lice does not appear to contribute much to heteroplasmy variation either as heteroplasmic sites are not limited to recombinant regions but spread across minichromosomes. More studies are needed to look into the issue of heteroplasmy in eukaryotes with deep sequencing technologies and to address the questions raised above in the current study.</p></sec><sec sec-type="supplementary-material" id="s5"><title>Supporting Information</title><supplementary-material content-type="local-data" id="pone.0073329.s001"><label>Table S1</label><caption><p><bold>Conserved and specific primers used in this study.</bold></p><p>(DOC)</p></caption><media xlink:href="pone.0073329.s001.doc"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0073329.s002"><label>Table S2</label><caption><p><bold>Heteroplasmic sites in mitochondrial tRNA genes of human body lice and human head lice.</bold></p><p>(DOC)</p></caption><media xlink:href="pone.0073329.s002.doc"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0073329.s003"><label>Table S3</label><caption><p><bold>Shared heteroplasmic sites in human lice.</bold></p><p>(DOC)</p></caption><media xlink:href="pone.0073329.s003.doc"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0073329.s004"><label>Table S4</label><caption><p><bold>Heteroplasmic sites in mitochondrial tRNA genes of ticks.</bold></p><p>(DOC)</p></caption><media xlink:href="pone.0073329.s004.doc"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0073329.s005"><label>Table S5</label><caption><p><bold>P<sub>f</sub> and Q<sub>f</sub> values for heteroplasmic sites at different sequence-read coverage (minimum variant 1.5%, sequencing-error rate 0.5%, π0 = 0.05).</bold></p><p>(DOC)</p></caption><media xlink:href="pone.0073329.s005.doc"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0073329.s006"><label>Table S6</label><caption><p><bold>Putative amino acid changes caused by heteroplasmy in the 12 human lice and 7 ticks.</bold> AA stands for amino acid. “+” stands for positive charge. “−” stands for negative charge.</p><p>(DOC)</p></caption><media xlink:href="pone.0073329.s006.doc"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material></sec></body><back><ack><p>We thank Kate Herd for assistance with PCR experiments and Tianfang Wang for assistance with QValue analysis. We thank the anonymous reviewers for comments that have greatly improved this paper. We thank Lidia Chitimia for the <italic>H. parva</italic>; Matsuo Kobayashi and Kyoko Sawabe for the <italic>H. formosensis</italic>; Dennis Poppi and Wayne Jovgenson for the <italic>R. microplus</italic> from Cambadia; Paulo Batista and Renato Andreaffi for the <italic>R. geigyi</italic> and <italic>O. megnini</italic>. We thank Tusunjiang for helping R.S. collect human lice in China and Kassahun Desalegn Bilcha, Jemal Ali and Richard Pollack for helping S.C.B collect human lice in Ethiopia.</p></ack><ref-list><title>References</title><ref id="pone.0073329-Lemire1"><label>1</label><mixed-citation publication-type="other">Lemire B (2005) Mitochondrial genetics. WormBook, ed. T<italic>he C. elegans Research Community</italic>, WormBook. /10.1895/wormbook.1.25.1,<ext-link ext-link-type="uri" xlink:href="http://www.wormbook.org">http://www.wormbook.org</ext-link></mixed-citation></ref><ref id="pone.0073329-Rantanen1"><label>2</label><mixed-citation publication-type="journal"><name><surname>Rantanen</surname><given-names>A</given-names></name>, <name><surname>Larsson</surname><given-names>N</given-names></name> (<year>2000</year>) <article-title>Regulation of mitochondrial DNA copy number during spermatogenesis</article-title>. <source>Hum Reprod</source><volume>15</volume>: <fpage>86</fpage>–<lpage>91</lpage><pub-id pub-id-type="pmid">11041516</pub-id></mixed-citation></ref><ref id="pone.0073329-Lightowlers1"><label>3</label><mixed-citation publication-type="journal"><name><surname>Lightowlers</surname><given-names>RN</given-names></name>, <name><surname>Chinnery</surname><given-names>PF</given-names></name>, <name><surname>Turnbull</surname><given-names>DM</given-names></name>, <name><surname>Howell</surname><given-names>N</given-names></name> (<year>1997</year>) <article-title>Mammalian mitochondrial genetics: heredity, heteroplasmy and disease</article-title>. <source>Trends Genet</source><volume>13</volume>: <fpage>450</fpage>–<lpage>455</lpage><pub-id pub-id-type="pmid">9385842</pub-id></mixed-citation></ref><ref id="pone.0073329-Boore1"><label>4</label><mixed-citation publication-type="journal"><name><surname>Boore</surname><given-names>JL</given-names></name> (<year>1999</year>) <article-title>Animal mitochondrial genomes</article-title>. <source>Nucleic Acids Res</source><volume>27</volume>: <fpage>1767</fpage>–<lpage>1780</lpage><pub-id pub-id-type="pmid">10101183</pub-id></mixed-citation></ref><ref id="pone.0073329-Song1"><label>5</label><mixed-citation publication-type="journal"><name><surname>Song</surname><given-names>GJ</given-names></name>, <name><surname>Lewis</surname><given-names>V</given-names></name> (<year>2008</year>) <article-title>Mitochondrial DNA integrity and copy number in sperm from infertile men</article-title>. <source>Fertil Steril</source><volume>90</volume>: <fpage>2238</fpage>–<lpage>2244</lpage><pub-id pub-id-type="pmid">18249386</pub-id></mixed-citation></ref><ref id="pone.0073329-Chiaratti1"><label>6</label><mixed-citation publication-type="journal"><name><surname>Chiaratti</surname><given-names>MR</given-names></name>, <name><surname>Meirelles</surname><given-names>FvV</given-names></name> (<year>2010</year>) <article-title>Mitochondrial DNA copy number, a marker of viability for oocytes</article-title>. <source>Biol Reprod</source><volume>83</volume>: <fpage>1</fpage>–<lpage>2</lpage><pub-id pub-id-type="pmid">20220127</pub-id></mixed-citation></ref><ref id="pone.0073329-Kmiec1"><label>7</label><mixed-citation publication-type="journal"><name><surname>Kmiec</surname><given-names>B</given-names></name>, <name><surname>Woloszynska</surname><given-names>M</given-names></name>, <name><surname>Janska</surname><given-names>H</given-names></name> (<year>2006</year>) <article-title>Heteroplasmy as a common state of mitochondrial genetic information in plants and animals</article-title>. <source>Curr Genet</source><volume>50</volume>: <fpage>149</fpage>–<lpage>159</lpage><pub-id pub-id-type="pmid">16763846</pub-id></mixed-citation></ref><ref id="pone.0073329-Light1"><label>8</label><mixed-citation publication-type="journal"><name><surname>Light</surname><given-names>JE</given-names></name>, <name><surname>Allen</surname><given-names>JM</given-names></name>, <name><surname>Long</surname><given-names>LM</given-names></name>, <name><surname>Carter</surname><given-names>TE</given-names></name>, <name><surname>Barrow</surname><given-names>L</given-names></name>, <etal>et al</etal> (<year>2008</year>) <article-title>Geographic distributions and origins of human head lice (<italic>Pediculus humanus capitis</italic>) based on mitochondrial data</article-title>. <source>J Parasitol</source><volume>94</volume>: <fpage>1275</fpage>–<lpage>1281</lpage><pub-id pub-id-type="pmid">18576877</pub-id></mixed-citation></ref><ref id="pone.0073329-He1"><label>9</label><mixed-citation publication-type="journal"><name><surname>He</surname><given-names>Y</given-names></name>, <name><surname>Wu</surname><given-names>J</given-names></name>, <name><surname>Dressman</surname><given-names>DC</given-names></name>, <name><surname>Iacobuzio-Donahue</surname><given-names>C</given-names></name>, <name><surname>Markowitz</surname><given-names>SD</given-names></name>, <etal>et al</etal> (<year>2010</year>) <article-title>Heteroplasmic mitochondrial DNA mutations in normal and tumor cells</article-title>. <source>Nature</source><volume>25</volume>: <fpage>610</fpage>–<lpage>614</lpage><pub-id pub-id-type="pmid">20200521</pub-id></mixed-citation></ref><ref id="pone.0073329-Goto1"><label>10</label><mixed-citation publication-type="journal"><name><surname>Goto</surname><given-names>H</given-names></name>, <name><surname>Dickins</surname><given-names>B</given-names></name>, <name><surname>Afgan</surname><given-names>E</given-names></name>, <name><surname>Paul</surname><given-names>IM</given-names></name>, <name><surname>Taylor</surname><given-names>J</given-names></name>, <etal>et al</etal> (<year>2011</year>) <article-title>Dynamics of mitochondrial heteroplasmy in three families investigated via a repeatable re-sequencing study</article-title>. <source>Genome Biol</source><volume>12</volume>: <fpage>R59</fpage><pub-id pub-id-type="pmid">21699709</pub-id></mixed-citation></ref><ref id="pone.0073329-Payne1"><label>11</label><mixed-citation publication-type="journal"><name><surname>Payne</surname><given-names>BA</given-names></name>, <name><surname>Wilson</surname><given-names>IJ</given-names></name>, <name><surname>Yu-Wai-Man</surname><given-names>P</given-names></name>, <name><surname>Coxhead</surname><given-names>J</given-names></name>, <name><surname>Deehan</surname><given-names>D</given-names></name>, <etal>et al</etal> (<year>2012</year>) <article-title>Universal heteroplasmy of human mitochondrial DNA</article-title>. <source>Hum Mol Genet</source><volume>22</volume>: <fpage>384</fpage>–<lpage>390</lpage><pub-id pub-id-type="pmid">23077218</pub-id></mixed-citation></ref><ref id="pone.0073329-Chinnery1"><label>12</label><mixed-citation publication-type="journal"><name><surname>Chinnery</surname><given-names>PF</given-names></name>, <name><surname>Samuels</surname><given-names>DC</given-names></name>, <name><surname>Elso</surname><given-names>J</given-names></name>, <name><surname>Turnbull</surname><given-names>DM</given-names></name> (<year>2002</year>) <article-title>Accumulation of mitochondrial DNA mutations in ageing, cancer, and mitochondrial disease: is there a common mechanism?</article-title><source>Lancet</source><volume>360</volume>: <fpage>1323</fpage>–<lpage>1325</lpage><pub-id pub-id-type="pmid">12414225</pub-id></mixed-citation></ref><ref id="pone.0073329-Santos1"><label>13</label><mixed-citation publication-type="journal"><name><surname>Santos</surname><given-names>C</given-names></name>, <name><surname>Sierra</surname><given-names>B</given-names></name>, <name><surname>Alvarez</surname><given-names>L</given-names></name>, <name><surname>Ramos</surname><given-names>A</given-names></name>, <name><surname>Fernández</surname><given-names>E</given-names></name>, <etal>et al</etal> (<year>2008</year>) <article-title>Frequency and pattern of heteroplasmy in the control region of human mitochondrial DNA</article-title>. <source>J Mol Evol</source><volume>67</volume>: <fpage>191</fpage>–<lpage>200</lpage><pub-id pub-id-type="pmid">18618067</pub-id></mixed-citation></ref><ref id="pone.0073329-Legros1"><label>14</label><mixed-citation publication-type="journal"><name><surname>Legros</surname><given-names>F</given-names></name>, <name><surname>Malka</surname><given-names>F</given-names></name>, <name><surname>Frachon</surname><given-names>P</given-names></name>, <name><surname>Manuel Rojo</surname><given-names>AL</given-names></name> (<year>2012</year>) <article-title>Organization and dynamics of human mitochondrial DNA</article-title>. <source>J Cell Sci</source><volume>125</volume>: <fpage>2653</fpage>–<lpage>2663</lpage></mixed-citation></ref><ref id="pone.0073329-Boyce1"><label>15</label><mixed-citation publication-type="journal"><name><surname>Boyce</surname><given-names>TM</given-names></name>, <name><surname>Zwick</surname><given-names>ME</given-names></name>, <name><surname>Aquadro</surname><given-names>CF</given-names></name> (<year>1989</year>) <article-title>Mitochondrial DNA in the bark weevils: Size, structure and heteroplasmy</article-title>. <source>Genetics</source><volume>123</volume>: <fpage>825</fpage>–<lpage>836</lpage><pub-id pub-id-type="pmid">2612897</pub-id></mixed-citation></ref><ref id="pone.0073329-Nardi1"><label>16</label><mixed-citation publication-type="journal"><name><surname>Nardi</surname><given-names>F</given-names></name>, <name><surname>Caraoelli</surname><given-names>A</given-names></name>, <name><surname>Fanciulli</surname><given-names>PP</given-names></name>, <name><surname>Dallai</surname><given-names>R</given-names></name>, <name><surname>Frati</surname><given-names>F</given-names></name> (<year>2001</year>) <article-title>The Complete mitochondrial DNA sequence of the basal hexapod <italic>Tetrodontophora bielanensis</italic>: Evidence for heteroplasmy and tRNA translocations</article-title>. <source>Mol Biol Evol</source><volume>18</volume>: <fpage>1293</fpage>–<lpage>1304</lpage><pub-id pub-id-type="pmid">11420368</pub-id></mixed-citation></ref><ref id="pone.0073329-Marcad1"><label>17</label><mixed-citation publication-type="journal"><name><surname>Marcadé</surname><given-names>I</given-names></name>, <name><surname>Cordaux</surname><given-names>R</given-names></name>, <name><surname>Doublet</surname><given-names>V</given-names></name>, <name><surname>Debenest</surname><given-names>C</given-names></name>, <name><surname>Bouchon</surname><given-names>D</given-names></name>, <etal>et al</etal> (<year>2007</year>) <article-title>Structure and evolution of the atypical mitochondrial genome of <italic>Armadillidium vulgare</italic> (Isopoda, Crustacea)</article-title>. <source>J Mol Evol</source><volume>65</volume>: <fpage>651</fpage>–<lpage>659</lpage><pub-id pub-id-type="pmid">17906827</pub-id></mixed-citation></ref><ref id="pone.0073329-Kann1"><label>18</label><mixed-citation publication-type="journal"><name><surname>Kann</surname><given-names>LM</given-names></name>, <name><surname>Rosenblum</surname><given-names>EB</given-names></name>, <name><surname>Rand</surname><given-names>DM</given-names></name> (<year>1997</year>) <article-title>Aging, mating, and the evolution of mt DNA heteroplasmy in <italic>Drosophila melanogaster</italic></article-title>. <source>Proc Nati Acad Sci U S A</source><volume>95</volume>: <fpage>2372</fpage>–<lpage>2377</lpage></mixed-citation></ref><ref id="pone.0073329-VolzLingenhhl1"><label>19</label><mixed-citation publication-type="journal"><name><surname>Volz-Lingenhöhl</surname><given-names>A</given-names></name>, <name><surname>Solignac</surname><given-names>M</given-names></name>, <name><surname>Sperlich</surname><given-names>D</given-names></name> (<year>1992</year>) <article-title>Stable heteroplasmy for a large-scale deletion in the coding region of <italic>Drosophila subobscura</italic> mitochondrial DNA</article-title>. <source>Proc Nati Acad Sci U S A</source><volume>89</volume>: <fpage>11528</fpage>–<lpage>32</lpage></mixed-citation></ref><ref id="pone.0073329-Leeuwen1"><label>20</label><mixed-citation publication-type="journal"><name><surname>Leeuwen</surname><given-names>TV</given-names></name>, <name><surname>Bsartel</surname><given-names>V</given-names></name>, <name><surname>Pottelberge</surname><given-names>SV</given-names></name>, <name><surname>Nieuwenhjuyse</surname><given-names>PV</given-names></name>, <name><surname>Nauen</surname><given-names>R</given-names></name>, <etal>et al</etal> (<year>2008</year>) <article-title>Mitochondrial heteroplasmy and the evolution of insecticide resistance: Non-Mendelian inheritance in action</article-title>. <source>Proc Nati Acad Sci U S A</source><volume>105</volume>: <fpage>5980</fpage>–<lpage>5985</lpage></mixed-citation></ref><ref id="pone.0073329-Lavrov1"><label>21</label><mixed-citation publication-type="journal"><name><surname>Lavrov</surname><given-names>DV</given-names></name> (<year>2007</year>) <article-title>Key transitions in animal evolution: a mitochondrial DNA perspective</article-title>. <source>Integr Comp Biol</source><volume>47</volume>: <fpage>734</fpage>–<lpage>743</lpage><pub-id pub-id-type="pmid">21669754</pub-id></mixed-citation></ref><ref id="pone.0073329-Shao1"><label>22</label><mixed-citation publication-type="journal"><name><surname>Shao</surname><given-names>R</given-names></name>, <name><surname>Kirkness</surname><given-names>EF</given-names></name>, <name><surname>Barker</surname><given-names>SC</given-names></name> (<year>2009</year>) <article-title>The single mitochondrial chromosome typical of animals has evolved into 18 minichromosomes in the human body louse, <italic>Pediculus humanus</italic></article-title>. <source>Genome Res</source><volume>19</volume>: <fpage>904</fpage>–<lpage>912</lpage><pub-id pub-id-type="pmid">19336451</pub-id></mixed-citation></ref><ref id="pone.0073329-Shao2"><label>23</label><mixed-citation publication-type="journal"><name><surname>Shao</surname><given-names>R</given-names></name>, <name><surname>Zhu</surname><given-names>XQ</given-names></name>, <name><surname>Barker</surname><given-names>SC</given-names></name>, <name><surname>Herd</surname><given-names>K</given-names></name> (<year>2012</year>) <article-title>Evolution of extensively fragmented mitochondrial genomes in the lice of humans</article-title>. <source>Genome Biol Evol</source><volume>4</volume>: <fpage>1088</fpage>–<lpage>1101</lpage><pub-id pub-id-type="pmid">23042553</pub-id></mixed-citation></ref><ref id="pone.0073329-Herd1"><label>24</label><mixed-citation publication-type="journal"><name><surname>Herd</surname><given-names>K</given-names></name>, <name><surname>Barker</surname><given-names>SC</given-names></name>, <name><surname>Shao</surname><given-names>R</given-names></name> (<year>2012</year>) <article-title>High-level of heteroplasmy in the mitochondrial cox1-minichromosome of the human body louse, <italic>Pediculus humanus</italic>, and the human head louse, <italic>Pediculus capitis</italic></article-title>. <source>Open Genomics J</source><volume>5</volume>: <fpage>14</fpage>–<lpage>17</lpage></mixed-citation></ref><ref id="pone.0073329-Black1"><label>25</label><mixed-citation publication-type="journal"><name><surname>Black</surname><given-names>WC</given-names></name>, <name><surname>Roehrdanz</surname><given-names>RL</given-names></name> (<year>1998</year>) <article-title>Mitochondrial gene order is not conserved in arthropods: prostriate and metastriate tick mitochondrial genomes</article-title>. <source>Mol Biol Evol</source><volume>15</volume>: <fpage>1772</fpage>–<lpage>1785</lpage><pub-id pub-id-type="pmid">9866211</pub-id></mixed-citation></ref><ref id="pone.0073329-Shao3"><label>26</label><mixed-citation publication-type="journal"><name><surname>Shao</surname><given-names>R</given-names></name>, <name><surname>Aoki</surname><given-names>Y</given-names></name>, <name><surname>Mitani</surname><given-names>H</given-names></name>, <name><surname>Tabuchi</surname><given-names>N</given-names></name>, <name><surname>Barker</surname><given-names>SC</given-names></name>, <etal>et al</etal> (<year>2004</year>) <article-title>The mitochondrial genomes of soft ticks have an arrangement of genes that has remained unchanged for over 400 million years</article-title>. <source>Insect Mol Biol</source><volume>13</volume>: <fpage>219</fpage>–<lpage>224</lpage><pub-id pub-id-type="pmid">15157222</pub-id></mixed-citation></ref><ref id="pone.0073329-Shao4"><label>27</label><mixed-citation publication-type="journal"><name><surname>Shao</surname><given-names>R</given-names></name>, <name><surname>Barker</surname><given-names>SC</given-names></name>, <name><surname>Mitani</surname><given-names>H</given-names></name>, <name><surname>Aoki</surname><given-names>Y</given-names></name>, <name><surname>Fukunaga</surname><given-names>M</given-names></name> (<year>2005</year>) <article-title>Evolution of duplicate control regions in the mitochondrial genomes of metazoa: a case study with Australasian Ixodes ticks</article-title>. <source>Mol Biol Evol</source><volume>22</volume>: <fpage>620</fpage>–<lpage>629</lpage><pub-id pub-id-type="pmid">15537802</pub-id></mixed-citation></ref><ref id="pone.0073329-Shao5"><label>28</label><mixed-citation publication-type="journal"><name><surname>Shao</surname><given-names>R</given-names></name>, <name><surname>Barker</surname><given-names>SC</given-names></name> (<year>2007</year>) <article-title>Mitochondrial genomes of parasitic arthropods: implications for studies of population genetics and evolution</article-title>. <source>Parasitology</source><volume>134</volume>: <fpage>153</fpage>–<lpage>167</lpage><pub-id pub-id-type="pmid">17032475</pub-id></mixed-citation></ref><ref id="pone.0073329-Burger1"><label>29</label><mixed-citation publication-type="journal"><name><surname>Burger</surname><given-names>TD</given-names></name>, <name><surname>Shao</surname><given-names>R</given-names></name>, <name><surname>Beati</surname><given-names>L</given-names></name>, <name><surname>Miller</surname><given-names>H</given-names></name>, <name><surname>Barker</surname><given-names>SC</given-names></name> (<year>2012</year>) <article-title>Phylogenetic analysis of ticks (Acari: Ixodida) using mitochondrial genomes and nuclear rRNA genes indicates that the genus <italic>Amblyomma</italic> is polyphyletic</article-title>. <source>Mol Phylogenet Evol</source><volume>64</volume>: <fpage>45</fpage>–<lpage>55</lpage><pub-id pub-id-type="pmid">22465402</pub-id></mixed-citation></ref><ref id="pone.0073329-Kambhampati1"><label>30</label><mixed-citation publication-type="journal"><name><surname>Kambhampati</surname><given-names>S</given-names></name>, <name><surname>Smith</surname><given-names>PT</given-names></name> (<year>1995</year>) <article-title>PCR primers for the amplification of four insect mitochondrial gene fragments</article-title>. <source>Insect Mol Biol</source><volume>4</volume>: <fpage>233</fpage>–<lpage>236</lpage><pub-id pub-id-type="pmid">8825760</pub-id></mixed-citation></ref><ref id="pone.0073329-Simon1"><label>31</label><mixed-citation publication-type="journal"><name><surname>Simon</surname><given-names>C</given-names></name>, <name><surname>Frati</surname><given-names>F</given-names></name>, <name><surname>Beckenbach</surname><given-names>A</given-names></name>, <name><surname>Crespi</surname><given-names>B</given-names></name>, <name><surname>Liu</surname><given-names>H</given-names></name>, <etal>et al</etal> (<year>1994</year>) <article-title>Evolution, weighting, and phylogenetic utility of mitochondrial gene-sequences and a compilation of conserved polymerases chain reaction primers</article-title>. <source>Ann Entomol Soc Am</source><volume>87</volume>: <fpage>651</fpage>–<lpage>701</lpage></mixed-citation></ref><ref id="pone.0073329-Altschul1"><label>32</label><mixed-citation publication-type="journal"><name><surname>Altschul</surname><given-names>SF</given-names></name>, <name><surname>Gish</surname><given-names>W</given-names></name>, <name><surname>Miller</surname><given-names>W</given-names></name>, <name><surname>Myers</surname><given-names>EW</given-names></name>, <name><surname>Lipman</surname><given-names>DJ</given-names></name> (<year>1990</year>) <article-title>Basic local alignment search tool</article-title>. <source>J Mol Biol</source><volume>215</volume>: <fpage>403</fpage>–<lpage>410</lpage><pub-id pub-id-type="pmid">2231712</pub-id></mixed-citation></ref><ref id="pone.0073329-McGinnisSMadden1"><label>33</label><mixed-citation publication-type="journal"><name><surname>McGinnis S & Madden</surname><given-names>TL</given-names></name> (<year>2004</year>) <article-title>BLAST: at the core of a powerful and diverse set of sequence analysis tools</article-title>. <source>Nucleic Acids Res</source><volume>32</volume>: <fpage>20</fpage>–<lpage>25</lpage></mixed-citation></ref><ref id="pone.0073329-LoweTM1"><label>34</label><mixed-citation publication-type="journal"><name><surname>Lowe TM &</surname><given-names>EddySR</given-names></name> (<year>1997</year>) <article-title>tRNAscan-SE: A Program for Improved Detection of Transfer RNA Genes in Genomic Sequence</article-title>. <source>Nucleic Acids Res</source><volume>25</volume>: <fpage>955</fpage>–<lpage>964</lpage><pub-id pub-id-type="pmid">9023104</pub-id></mixed-citation></ref><ref id="pone.0073329-LaslettDCanback1"><label>35</label><mixed-citation publication-type="journal"><name><surname>Laslett D & Canback</surname><given-names>B</given-names></name> (<year>2008</year>) <article-title>ARWEN: a program to detect tRNA genes in metazoan mitochondrial nucleotide sequences</article-title>. <source>Bioinformatics</source><volume>24(2)</volume>: <fpage>172</fpage>–<lpage>175</lpage><pub-id pub-id-type="pmid">18033792</pub-id></mixed-citation></ref><ref id="pone.0073329-Storey1"><label>36</label><mixed-citation publication-type="journal"><name><surname>Storey</surname><given-names>JD</given-names></name> (<year>2002</year>) <article-title>A direct approach to false discovery rates</article-title>. <source>J R Stat Soc Series B</source><volume>64</volume>: <fpage>479</fpage>–<lpage>498</lpage></mixed-citation></ref><ref id="pone.0073329-Storey2"><label>37</label><mixed-citation publication-type="journal"><name><surname>Storey</surname><given-names>JD</given-names></name>, <name><surname>Tibshirani</surname><given-names>R</given-names></name> (<year>2003</year>) <article-title>Statistical significance for genome-wide studies</article-title>. <source>Proc Nati Acad Sci U S A</source><volume>100</volume>: <fpage>9440</fpage>–<lpage>9445</lpage></mixed-citation></ref><ref id="pone.0073329-Nakamura1"><label>38</label><mixed-citation publication-type="journal"><name><surname>Nakamura</surname><given-names>K</given-names></name>, <name><surname>Oshima</surname><given-names>T</given-names></name>, <name><surname>Morimoto</surname><given-names>T</given-names></name>, <name><surname>Ikeda</surname><given-names>S</given-names></name>, <name><surname>Yoshikawa</surname><given-names>H</given-names></name>, <etal>et al</etal> (<year>2011</year>) <article-title>Sequence-specific error profile of Illumina sequencers</article-title>. <source>Nucleic Acids Res</source><volume>39</volume>: <fpage>e90</fpage><pub-id pub-id-type="pmid">21576222</pub-id></mixed-citation></ref><ref id="pone.0073329-Minoche1"><label>39</label><mixed-citation publication-type="journal"><name><surname>Minoche</surname><given-names>A</given-names></name>, <name><surname>Dohm</surname><given-names>JC</given-names></name>, <name><surname>Himmelbauer</surname><given-names>H</given-names></name> (<year>2011</year>) <article-title>Evaluation of genomic high-throughput sequencing data generated on Illumina HiSeq and Genome Analyzer systems</article-title>. <source>Genome Biol</source><volume>12</volume>: <fpage>R112</fpage><pub-id pub-id-type="pmid">22067484</pub-id></mixed-citation></ref><ref id="pone.0073329-Sachadyn1"><label>40</label><mixed-citation publication-type="journal"><name><surname>Sachadyn</surname><given-names>P</given-names></name>, <name><surname>Zhang</surname><given-names>XM</given-names></name>, <name><surname>Clark</surname><given-names>LD</given-names></name>, <name><surname>Naviaux</surname><given-names>RK</given-names></name>, <name><surname>Heber-Katz</surname><given-names>E</given-names></name> (<year>2008</year>) <article-title>Naturally-Occurring Mitochondrial DNA heteroplasmy in the MRL mouse</article-title>. <source>Mitochondrion</source><volume>8</volume>: <fpage>358</fpage>–<lpage>366</lpage><pub-id pub-id-type="pmid">18761428</pub-id></mixed-citation></ref><ref id="pone.0073329-Messenger1"><label>41</label><mixed-citation publication-type="journal"><name><surname>Messenger</surname><given-names>LA</given-names></name>, <name><surname>Llewellyn</surname><given-names>MS</given-names></name>, <name><surname>Bhattacharyya</surname><given-names>T</given-names></name>, <name><surname>Franze</surname><given-names>O</given-names></name>, <name><surname>Lewis</surname><given-names>MD</given-names></name>, <etal>et al</etal> (<year>2012</year>) <article-title>Multiple mitochondrial introgression events and heteroplasmy in <italic>Trypanosoma cruzi</italic> revealed by maxicircle MLST and next generation sequencing</article-title>. <source>PLoS Negl Trop Dis</source><volume>6</volume>: <fpage>1584</fpage>–<lpage>1596</lpage></mixed-citation></ref><ref id="pone.0073329-Quinlan1"><label>42</label><mixed-citation publication-type="journal"><name><surname>Quinlan</surname><given-names>AR</given-names></name>, <name><surname>Stewart</surname><given-names>DA</given-names></name>, <name><surname>Stromberg</surname><given-names>MP</given-names></name>, <name><surname>Marth</surname><given-names>GT</given-names></name> (<year>2008</year>) <article-title>Pyrobayes: Animproved base caller for SNP discovery in pyrosequences</article-title>. <source>Nat Methods</source><volume>5</volume>: <fpage>179</fpage>–<lpage>181</lpage><pub-id pub-id-type="pmid">18193056</pub-id></mixed-citation></ref><ref id="pone.0073329-Sosa1"><label>43</label><mixed-citation publication-type="journal"><name><surname>Sosa</surname><given-names>MX</given-names></name>, <name><surname>Sivakumar</surname><given-names>KA</given-names></name>, <name><surname>Maragh</surname><given-names>S</given-names></name>, <name><surname>Veeramachaneni</surname><given-names>V</given-names></name>, <name><surname>Hariharan</surname><given-names>R</given-names></name>, <etal>et al</etal> (<year>2012</year>) <article-title>Next-generation sequencing of human mitochondrial reference genomes uncovers high heteroplasmy frequency</article-title>. <source>PLoS Comput Biol</source><volume>8</volume>: <fpage>e1002737</fpage><pub-id pub-id-type="pmid">23133345</pub-id></mixed-citation></ref><ref id="pone.0073329-Doan1"><label>44</label><mixed-citation publication-type="other">Doan R, Cohen ND, Sawyer J, Ghaffari N, Johnson CD, <etal>et al</etal>.. (2012) Whole-genome sequencing and genetic variant analysis of a quarter Horse mare. BMC Genomics 13 10.1186/1471-2164-13-78.</mixed-citation></ref><ref id="pone.0073329-Li1"><label>45</label><mixed-citation publication-type="journal"><name><surname>Li</surname><given-names>MK</given-names></name>, <name><surname>Schönberg</surname><given-names>A</given-names></name>, <name><surname>Schaefer</surname><given-names>M</given-names></name>, <name><surname>Schroeder</surname><given-names>R</given-names></name>, <name><surname>Nasidze</surname><given-names>I</given-names></name>, <etal>et al</etal> (<year>2010</year>) <article-title>Detecting heteroplasmy from high-throughput sequencing of complete human mitochondrial DNA genomes</article-title>. <source>Am J Hum Genet</source><volume>87</volume>: <fpage>237</fpage>–<lpage>249</lpage><pub-id pub-id-type="pmid">20696290</pub-id></mixed-citation></ref></ref-list></back></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Asie/explor/AustralieFrV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 002B15 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 002B15 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Wicri/Asie
|area= AustralieFrV1
|flux= Pmc
|étape= Corpus
|type= RBID
|clé=
|texte=
}}
| This area was generated with Dilib version V0.6.33. |  |