In order to compare phylogenetic methods and to reconstruct the evolutionary history of the tubercle bacilli, a set of macro-array-based genotyping data of Mycobacterium tuberculosis clinical isolates (called spoligotyping for spacer oligonucleotide typing, which assays the variability of the Direct Repeat -DR- locus), was analyzed in four settings of the Caribbean region (Guadeloupe, Martinique, Cuba and Haiti). A set of 47 alleles, split into 26 shared and 21 unique alleles) representative of 321 individual M. tuberculosis clinical isolates from patients residing in the above regions was studied. The following methods (and software in brackets) were investigated: numerical taxonomy distance methods (TAXOTRON), maximum parsimony procedure (PAUP), median-joining networks (NETWORK), and nested clade analysis (GEODIS). Results using these methods were analyzed, compared and discussed. The latter method (GEODIS) was investigated in detail by introducing geographical data together with genetic variability results to detect a link between population structure and population history, and to test the null hypothesis of no association between geography and genotypes. Irrespective of the methods used, our findings demonstrate that a core structure of four families (or clades) of M. tuberculosis strains is highly prevalent within the islands studied, indirectly reflecting passed colonization history of these different settings. Specificity of M. tuberculosis genotypes in each of the islands is discussed in the light of their respective colonial and contemporary histories.