Mycobacterium tuberculosis lineage 4 comprises globally distributed and geographically restricted sublineages

David Stucki, Daniela Brites, Leïla Jeljeli, Mireia Coscolla, Qingyun Liu, Andrej Trauner, Lukas Fenner, Liliana Rutaihwa, Sonia Borrell, Tao Luo, Qian Gao, Midori Kato-Maeda, Marie Ballif, Matthias Egger, Rita Macedo, Helmi Mardassi, Milagros Moreno, Griselda Tudo Vilanova, Janet Fyfe, Maria Globan, Jackson Thomas, Frances Jamieson, Jennifer L Guthrie, Adwoa Asante-Poku, Dorothy Yeboah-Manu, Eddie Wampande, Willy Ssengooba, Moses Joloba, W Henry Boom, Indira Basu, James Bower, Margarida Saraiva, Sidra E G Vaconcellos, Philip Suffys, Anastasia Koch, Robert Wilkinson, Linda Gail-Bekker, Bijaya Malla, Serej D Ley, Hans-Peter Beck, Bouke C de Jong, Kadri Toit, Elisabeth Sanchez-Padilla, Maryline Bonnet, Ana Gil-Brusola, Matthias Frank, Veronique N Penlap Beng, Kathleen Eisenach, Issam Alani, Perpetual Wangui Ndung'u, Gunturu Revathi, Florian Gehre, Suriya Akter, Francine Ntoumi, Lynsey Stewart-Isherwood, Nyanda E Ntinginya, Andrea Rachow, Michael Hoelscher, Daniela Maria Cirillo, Girts Skenders, Sven Hoffner, Daiva Bakonyte, Petras Stakenas, Roland Diel, Valeriu Crudu, Olga Moldovan, Sahal Al-Hajoj, Larissa Otero, Francesca Barletta, E Jane Carter, Lameck Diero, Philip Supply, Iñaki Comas, Stefan Niemann, Sebastien Gagneux, David Stucki, Daniela Brites, Leïla Jeljeli, Mireia Coscolla, Qingyun Liu, Andrej Trauner, Lukas Fenner, Liliana Rutaihwa, Sonia Borrell, Tao Luo, Qian Gao, Midori Kato-Maeda, Marie Ballif, Matthias Egger, Rita Macedo, Helmi Mardassi, Milagros Moreno, Griselda Tudo Vilanova, Janet Fyfe, Maria Globan, Jackson Thomas, Frances Jamieson, Jennifer L Guthrie, Adwoa Asante-Poku, Dorothy Yeboah-Manu, Eddie Wampande, Willy Ssengooba, Moses Joloba, W Henry Boom, Indira Basu, James Bower, Margarida Saraiva, Sidra E G Vaconcellos, Philip Suffys, Anastasia Koch, Robert Wilkinson, Linda Gail-Bekker, Bijaya Malla, Serej D Ley, Hans-Peter Beck, Bouke C de Jong, Kadri Toit, Elisabeth Sanchez-Padilla, Maryline Bonnet, Ana Gil-Brusola, Matthias Frank, Veronique N Penlap Beng, Kathleen Eisenach, Issam Alani, Perpetual Wangui Ndung'u, Gunturu Revathi, Florian Gehre, Suriya Akter, Francine Ntoumi, Lynsey Stewart-Isherwood, Nyanda E Ntinginya, Andrea Rachow, Michael Hoelscher, Daniela Maria Cirillo, Girts Skenders, Sven Hoffner, Daiva Bakonyte, Petras Stakenas, Roland Diel, Valeriu Crudu, Olga Moldovan, Sahal Al-Hajoj, Larissa Otero, Francesca Barletta, E Jane Carter, Lameck Diero, Philip Supply, Iñaki Comas, Stefan Niemann, Sebastien Gagneux

Abstract

Generalist and specialist species differ in the breadth of their ecological niches. Little is known about the niche width of obligate human pathogens. Here we analyzed a global collection of Mycobacterium tuberculosis lineage 4 clinical isolates, the most geographically widespread cause of human tuberculosis. We show that lineage 4 comprises globally distributed and geographically restricted sublineages, suggesting a distinction between generalists and specialists. Population genomic analyses showed that, whereas the majority of human T cell epitopes were conserved in all sublineages, the proportion of variable epitopes was higher in generalists. Our data further support a European origin for the most common generalist sublineage. Hence, the global success of lineage 4 reflects distinct strategies adopted by different sublineages and the influence of human migration.

Conflict of interest statement

Statement The authors have no competing interests as defined by Springer Nature, or other interests that might be perceived to influence the results and/or discussion reported in this paper.

Figures

Figure 1. Definition and global frequency of…
Figure 1. Definition and global frequency of Lineage 4 sublineages.
(a) We defined 10 sublineages based on the analysis of 72 MTBC Lineage 4 genome sequences published previously,. Sublineages were labeled according to Coll et al. (whenever possible) and previous designations based on spoligotyping (see Supplementary Fig. 1). Black triangles indicate sublineages identified as specialists, black circles indicate generalists. Filled shapes indicate sublineages, for which we performed deep genomic analyses. (b) Global proportion of each sublineage. A total of 3,366 MTBC Lineage 4 isolates were screened for sublineage-specific SNPs. L4.3/LAM was the most frequent sublineage globally.
Figure 2. Global distribution of Lineage 4…
Figure 2. Global distribution of Lineage 4 sublineages.
Pie charts showing proportions of the 10 Lineage 4 sublineages among all MTBC Lineage 4 isolates in each country. Circle sizes correspond to the number of isolates analyzed per country. A total of 3,366 MTBC Lineage 4 isolates were included. Color codes are as in Fig. 1.
Figure 3. Country-specific proportions of sublineages reveal…
Figure 3. Country-specific proportions of sublineages reveal generalists and specialists.
(a) The generalist sublineages L4.1.2/Haarlem, L4.3/LAM and L4.10/PGG3 were found globally at high proportions. (b) The locally restricted specialist sublineages L4.1.3/Ghana, L4.5, L4.6.1/Uganda and L4.6.2/Cameroon occurred at high frequencies in only a few countries and were restricted to certain geographical regions. Intensity of red indicates proportion of the sublineage among all Lineage 4 isolates in each country. Countries with fewer than three isolates in total are shown as “no data” and are filled white. A total of 3,366 Lineage 4 isolates were included in this analysis. The color scale for all sublineages is as indicated in Panel a, except for sublineage L4.1.3/Ghana (separate scale shown).
Figure 4. Pair-wise ratios of rates of…
Figure 4. Pair-wise ratios of rates of nonsynonymous to synonymous substitutions (dN/dS) in generalist and specialist sublineages for different gene categories.
Abbreviations: Epi experimentally confirmed human T cell epitopes; nEpi – non-epitope regions of T-cell antigens, both obtained from the Immune Epitope Database; Ess – essential genes; nEss – non-essential genes. Wilcoxon rank sum tests: L4.6.1/Uganda (N=203) Epi vs nEpi, W=4952, p<0.001; L4.6.1/Uganda (N=203) Ess vs nEss, W=1415, p<0.001; L4.3/LAM (N=293) Epi vs nEpi, W=74540, p<0.001, L4.3/LAM (n=293) Ess vs nEss W=45067, p-value=0.29; L4.1.2/Haarlem (N=228) Epi vs nEpi, W=6561, p<0.001, L4.1.2/Haarlem (N=228) Ess vs nEss W=13369, p<0.001; L4.10/PGG3 (N=301) Epi vs nEpi, W= 27335, p<0.001, L4.10/PGG3 (N=301) Ess vs nEss W= 3103, p<0.001.
Figure 5. Frequency distribution of the number…
Figure 5. Frequency distribution of the number of epitopes with nonsynonymous variants in generalist and specialist sublineages.
A total of 1,226 T cell epitopes were included in the analysis. The number above each bar corresponds to epitope counts. Generalist sublineages L4.3/LAM, L4.1.2/Haarlem and (L4.10/PGG3. Specialist sublineage L4.6.1/Uganda. Tests: L4.6.1/Uganda vs L4.3/LAM Χ2= 27.04, p<0.001; L4.6.1/Uganda vs L4.1.2/Haarlem Χ2=15.75, p<0.001; L4.6.1/Uganda vs L4.1.2/PGG3 Χ2= 68.24, p<0.001.
Figure 6. Genome-based phylogeny and diversity by…
Figure 6. Genome-based phylogeny and diversity by continent of 293 strains of the L4.3/LAM sublineage.
(a) Bayesian phylogeny with label colors indicating continent of strain origin: blue, Europe/Mediterranean; red, Sub-Saharan Africa; yellow, America; pink, Asia. Numbers on nodes indicate posterior probabilities. Pie charts indicate reconstructed ancestral geographical regions of the internal nodes. The hypothetical L4.3/LAM-ancestor is labeled and a European origin for this ancestor was supported using a Bayesian Method (shown) and a Maximum Parsimony method (Supplementary Fig. 14). The pie colors correspond to the colors of the taxa labels. (b) Boxplot of pairwise genetic distances (number of polymorphisms) of L4.3/LAM strains by continent (p-values from Wilcoxon rank sum test). (c) Nucleotide diversity per site (π), measured by continent. Error bars indicate 95% confidence intervals. MTBC isolates from countries of the continent group “Oceania“ (UN category; including Australia and New Zealand, Melanesia, Micronesia and Polynesia) were excluded for the genetic diversity analysis in panels B and C due the low number of samples.

References

    1. Futuyma DJ, Moreno G. The Evolution of Ecological Specialization. Annual Review of Ecology and Systematics. 1988;19:207–233.
    1. Woolhouse ME, Webster JP, Domingo E, Charlesworth B, Levin BR. Biological and biomedical implications of the co-evolution of pathogens and their hosts. Nat Genet. 2002;32:569–77.
    1. Woolhouse ME, Taylor LH, Haydon DT. Population biology of multihost pathogens. Science. 2001;292:1109–12.
    1. Kirzinger MW, Stavrinides J. Host specificity determinants as a genetic continuum. Trends Microbiol. 2012;20:88–93.
    1. Vouga M, Greub G. Emerging bacterial pathogens: past and beyond. Clin Microbiol Infect. 2016;22:12–21.
    1. Brites D, Gagneux S. Co-evolution of Mycobacterium tuberculosis and Homo sapiens. Immunol Rev. 2015;264:6–24.
    1. Borrell S, Gagneux S. Infectiousness, reproductive fitness and evolution of drug-resistant Mycobacterium tuberculosis. Int J Tuberc Lung Dis. 2009;13:1456–1466.
    1. Merker M, et al. Evolutionary history and global spread of the Mycobacterium tuberculosis Beijing lineage. Nat Genet. 2015;47:242–9.
    1. Luo T, et al. Southern East Asian origin and coexpansion of Mycobacterium tuberculosis Beijing family with Han Chinese. Proc Natl Acad Sci U S A. 2015;112:8136–41.
    1. Cowley D, et al. Recent and rapid emergence of W-Beijing strains of Mycobacterium tuberculosis in Cape Town, South Africa. Clin Infect Dis. 2008;47:1252–9.
    1. de Jong BC, Antonio M, Gagneux S. Mycobacterium africanum--review of an important cause of human tuberculosis in West Africa. PLoS Negl Trop Dis. 2010;4:e744.
    1. Firdessa R, et al. Mycobacterial lineages causing pulmonary and extrapulmonary tuberculosis, Ethiopia. Emerg Infect Dis. 2013;19:460–463.
    1. Gagneux S. Host-pathogen coevolution in human tuberculosis. Philos Trans R Soc Lond B Biol Sci. 2012;367:850–9.
    1. Fenner L, et al. HIV infection disrupts the sympatric host-pathogen relationship in human tuberculosis. PLoS Genet. 2013;9:e1003318.
    1. Gagneux S, et al. Variable host-pathogen compatibility in Mycobacterium tuberculosis. Proc Natl Acad Sci U S A. 2006;103:2869–2873.
    1. Baker L, Brown T, Maiden MC, Drobniewski F. Silent nucleotide polymorphisms and a phylogeny for Mycobacterium tuberculosis. Emerg Infect Dis. 2004;10:1568–77.
    1. Reed MB, et al. Major Mycobacterium tuberculosis lineages associate with patient country of origin. J Clin Microbiol. 2009;47:1119–28.
    1. Demay C, et al. SITVITWEB - A publicly available international multimarker database for studying Mycobacterium tuberculosis genetic diversity and molecular epidemiology. Infect Genet Evol. 2012;12:755–66.
    1. Coscolla M, Gagneux S. Consequences of genomic diversity in Mycobacterium tuberculosis. Semin Immunol. 2014;26:431–44.
    1. Coscolla M, Gagneux S. Does M. tuberculosis genomic diversity explain disease diversity? Drug Discov Today Dis Mech. 2010;7:e43–e59.
    1. Comas I, et al. Out-of-Africa migration and Neolithic coexpansion of Mycobacterium tuberculosis with modern humans. Nat Genet. 2013;45:1176–82.
    1. Coscolla M, et al. Novel Mycobacterium tuberculosis complex isolate from a wild chimpanzee. Emerg Infect Dis. 2013;19:969–76.
    1. Abadia E, et al. Resolving lineage assignation on Mycobacterium tuberculosis clinical isolates classified by spoligotyping with a new high-throughput 3R SNPs-based method. Infect Genet Evol. 2010;10:1066–74.
    1. Filliol I, et al. Global Phylogeny of Mycobacterium tuberculosis based on single nucleotide polymorphism (SNP) analysis: insights into tuberculosis evolution, phylogenetic accuracy of other DNA fingerprinting systems, and recommendations for a minimal standard SNP set. J Bacteriol. 2006;188:759–72.
    1. Sreevatsan S, et al. Restricted structural gene polymorphism in the Mycobacterium tuberculosis complex indicates evolutionarily recent global dissemination. Proc Natl Acad Sci U S A. 1997;94:9869–74.
    1. Homolka S, et al. High resolution discrimination of clinical Mycobacterium tuberculosis complex strains based on single nucleotide polymorphisms. PLoS One. 2012;7:e39855.
    1. Coll F, et al. A robust SNP barcode for typing Mycobacterium tuberculosis complex strains. Nat Commun. 2014;5:4812.
    1. Tsolaki AG, et al. Functional and evolutionary genomics of Mycobacterium tuberculosis: insights from genomic deletions in 100 strains. Proc Natl Acad Sci U S A. 2004;101:4865–70.
    1. Comas I, Homolka S, Niemann S, Gagneux S. Genotyping of genetically monomorphic bacteria: DNA sequencing in Mycobacterium tuberculosis highlights the limitations of current methodologies. PLoS ONE. 2009;4:e7815.
    1. Yeboah-Manu D, et al. Genotypic diversity and drug susceptibility patterns among M. tuberculosis complex isolates from South-Western Ghana. PLoS One. 2011;6:e21906.
    1. Malla B, et al. First insights into the phylogenetic diversity of Mycobacterium tuberculosis in Nepal. PLoS One. 2012;7:e52297.
    1. Fenner L, et al. Mycobacterium tuberculosis transmission in a country with low tuberculosis incidence: role of immigration and HIV infection. J Clin Microbiol. 2012;50:388–95.
    1. Ballif M, et al. Genetic diversity of Mycobacterium tuberculosis in Madang, Papua New Guinea. Int J Tuberc Lung Dis. 2012;16:1100–7.
    1. Wampande E, et al. Long-term dominance of Mycobacterium tuberculosis Uganda family in peri-urban Kampala-Uganda is not associated with cavitary disease. BMC Infect Dis. 2013;13:484.
    1. Ley SD, et al. Diversity of Mycobacterium tuberculosis and drug resistance in different provinces of Papua New Guinea. BMC Microbiol. 2014;14:307.
    1. Stucki D, et al. Two new rapid SNP-typing methods for classifying Mycobacterium tuberculosis complex into the main phylogenetic lineages. PLoS ONE. 2012;7:e41253.
    1. Bouakaze C, et al. Matrix-assisted laser desorption ionization-time of flight mass spectrometry-based single nucleotide polymorphism genotyping assay using iPLEX gold technology for identification of Mycobacterium tuberculosis complex species and lineages. J Clin Microbiol. 2011;49:3292–9.
    1. World Health Organization. Global tuberculosis control - surveillance, planning, financing. WHO; Geneva, Switzerland: 2015.
    1. Moller M, de Wit E, Hoal EG. Past, present and future directions in human genetic susceptibility to tuberculosis. FEMS Immunol Med Microbiol. 2010;58:3–26.
    1. Caws M, et al. The influence of host and bacterial genotype on the development of disseminated disease with Mycobacterium tuberculosis. PLoS Pathog. 2008;4:e1000034.
    1. Asante-Poku A, et al. Mycobacterium africanum is associated with patient ethnicity in Ghana. PLoS Negl Trop Dis. 2015;9:e3370.
    1. Ng PC, Henikoff S. SIFT: Predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003;31:3812–4.
    1. Ellstrand NC, Elam DR. Population genetic consequences of small population size: implications for plant conservation. Annual Review of Ecology and Systematics. 1993;24:217–242.
    1. Lanzas F, Karakousis PC, Sacchettini JC, Ioerger TR. Multidrug-resistant tuberculosis in panama is driven by clonal expansion of a multidrug-resistant Mycobacterium tuberculosis strain related to the KZN extensively drug-resistant M. tuberculosis strain from South Africa. J Clin Microbiol. 2013;51:3277–85.
    1. Bryant JM, et al. Inferring patient to patient transmission of Mycobacterium tuberculosis from whole genome sequencing data. BMC Infect Dis. 2013;13:110.
    1. Bryant JM, et al. Whole-genome sequencing to establish relapse or re-infection with Mycobacterium tuberculosis: a retrospective observational study. Lancet Respir Med. 2013;1:786–92.
    1. Gardy JL, et al. Whole-genome sequencing and social-network analysis of a tuberculosis outbreak. N Engl J Med. 2011;364:730–9.
    1. Casali N, et al. Evolution and transmission of drug-resistant tuberculosis in a Russian population. Nat Genet. 2014;46:279–86.
    1. Walker TM, et al. Whole-genome sequencing to delineate Mycobacterium tuberculosis outbreaks: a retrospective observational study. Lancet Infect Dis. 2012
    1. Roetzer A, et al. Whole genome sequencing versus traditional genotyping for investigation of a Mycobacterium tuberculosis outbreak: a longitudinal molecular epidemiological study. PLoS Med. 2013;10:e1001387.
    1. Farhat MR, et al. Genomic analysis identifies targets of convergent positive selection in drug-resistant Mycobacterium tuberculosis. Nat Genet. 2013;45:1183–9.
    1. Clark TG, et al. Elucidating emergence and transmission of multidrug-resistant tuberculosis in treatment experienced patients by whole genome sequencing. PLoS One. 2013;8:e83012.
    1. Jamieson FB, et al. Whole-genome sequencing of the Mycobacterium tuberculosis Manila sublineage results in less clustering and better resolution than mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) typing and spoligotyping. J Clin Microbiol. 2014;52:3795–8.
    1. Perez-Lago L, et al. Whole genome sequencing analysis of intrapatient microevolution in Mycobacterium tuberculosis: potential impact on the inference of tuberculosis transmission. J Infect Dis. 2014;209:98–108.
    1. Guerra-Assuncao JA, et al. Recurrence due to relapse or reinfection with Mycobacterium tuberculosis: a whole-genome sequencing approach in a large, population-based cohort with a high HIV infection prevalence and active follow-up. J Infect Dis. 2015;211:1154–63.
    1. Guerra-Assuncao JA, et al. Large-scale whole genome sequencing of M. tuberculosis provides insights into transmission in a high prevalence area. Elife. 2015;4
    1. Comas I, et al. Human T cell epitopes of Mycobacterium tuberculosis are evolutionarily hyperconserved. Nat Genet. 2010;42:498–503.
    1. Coscolla M, et al. M. tuberculosis T cell epitope analysis reveals paucity of antigenic variation and identifies rare variable TB antigens. Cell Host Microbe. 2015;18:538–48.
    1. Deitsch KW, Lukehart SA, Stringer JR. Common strategies for antigenic variation by bacterial, fungal and protozoan pathogens. Nat Rev Microbiol. 2009;7:493–503.
    1. Ernst JD, et al. Meeting report: NIH Workshop on the Tuberculosis Immune Epitope Database. Tuberculosis (Edinb) 2008;88:366–70.
    1. Pepperell CS, et al. The role of selection in shaping diversity of natural M. tuberculosis populations. PLoS Pathog. 2013;9:e1003543.
    1. Zhang YJ, et al. Global assessment of genomic regions required for growth in Mycobacterium tuberculosis. PLoS Pathog. 2012;8:e1002946.
    1. Esparza M, et al. PstS-1, the 38-kDa Mycobacterium tuberculosis glycoprotein, is an adhesin, which binds the macrophage mannose receptor and promotes phagocytosis. Scand J Immunol. 2015;81:46–55.
    1. Bekmurzayeva A, Sypabekova M, Kanayeva D. Tuberculosis diagnosis using immunodominant, secreted antigens of Mycobacterium tuberculosis. Tuberculosis (Edinb) 2013;93:381–8.
    1. Nagai S, Wiker HG, Harboe M, Kinomoto M. Isolation and partial characterization of major protein antigens in the culture fluid of Mycobacterium tuberculosis. Infect Immun. 1991;59:372–82.
    1. Juarez MD, Torres A, Espitia C. Characterization of the Mycobacterium tuberculosis region containing the mpt83 and mpt70 genes. FEMS Microbiol Lett. 2001;203:95–102.
    1. Coppola M, et al. Synthetic long peptide derived from Mycobacterium tuberculosis latency antigen Rv1733c protects against tuberculosis. Clin Vaccine Immunol. 2015;22:1060–9.
    1. Araujo LS, et al. Profile of interferon-gamma response to latency-associated and novel in vivo expressed antigens in a cohort of subjects recently exposed to Mycobacterium tuberculosis. Tuberculosis (Edinb) 2015;95:751–7.
    1. Yu Y, Harris AJ, Blair C, He X. RASP (Reconstruct Ancestral State in Phylogenies): a tool for historical biogeography. Mol Phylogenet Evol. 2015;87:46–9.
    1. Brudey K, et al. Mycobacterium tuberculosis complex genetic diversity : mining the fourth international spoligotyping database (SpolDB4) for classification, population genetics and epidemiology. BMC Microbiol. 2006;6:23.
    1. Comas I, et al. Population Genomics of Mycobacterium tuberculosis in Ethiopia contradicts the Virgin Soil hypothesis for human tuberculosis in Sub-Saharan Africa. Curr Biol. 2015;25:3260–6.
    1. Bos KI, et al. Pre-Columbian mycobacterial genomes reveal seals as a source of New World human tuberculosis. Nature. 2014;514:494–7.
    1. Kay GL, et al. Eighteenth-century genomes show that mixed infections were common at time of peak tuberculosis in Europe. Nat Commun. 2015;6:6717.
    1. Lazzarini LC, et al. Discovery of a novel Mycobacterium tuberculosis lineage that is a major cause of tuberculosis in Rio de Janeiro, Brazil. J Clin Microbiol. 2007;45:3891–902.
    1. Perdigao J, et al. Unraveling Mycobacterium tuberculosis genomic diversity and evolution in Lisbon, Portugal, a highly drug resistant setting. BMC Genomics. 2014;15:991.
    1. Mokrousov I, et al. Latin-American-Mediterranean lineage of Mycobacterium tuberculosis: Human traces across pathogen's phylogeography. Mol Phylogenet Evol. 2016;99:133–43.
    1. .
    1. Bates JH, Stead WW. The history of tuberculosis as a global epidemic. Med Clin North Am. 1993;77:1205–17.
    1. Gagneux S, Small PM. Global phylogeography of Mycobacterium tuberculosis and implications for tuberculosis product development. Lancet Infect Dis. 2007;7:328–37.
    1. Steiner A, Stucki D, Coscolla M, Borrell S, Gagneux S. KvarQ: targeted and direct variant calling from fastq reads of bacterial genomes. BMC Genomics. 2014;15:881.
    1. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.
    1. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38:e164.
    1. Tamura K, et al. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.
    1. Paradis E, Claude J, Strimmer K. APE: Analyses of Phylogenetics and Evolution in R language. Bioinformatics. 2004;20:289–90.
    1. Excoffier L, Laval G, Schneider S. Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinform Online. 2005;1:47–50.
    1. Waterhouse AM, Procter JB, Martin DM, Clamp M, Barton GJ. Jalview Version 2--a multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009;25:1189–91.
    1. Maddison W, Maddison D. Mesquite: a modular system for evolutionary analysis. 2001
    1. Lew JM, Kapopoulou A, Jones LM, Cole ST. TubercuList--10 years after. Tuberculosis (Edinb) 2011;91:1–7.
    1. Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572–4.
    1. Hartl D, Clarck AG. Principles of population genetics. Sinauer; 2007.
    1. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.

Source: PubMed

3
Suscribir