Haplotype network estimation: The intraspecific haplotype network was reconstructed using the algorithm described by Templeton et al (1992). The method was designed to estimate phylogenies at low levels of divergence and begins by calculating the overall limits of parsimony to reconstruct intraspecific relationships of haplotypes using a statistic from neutral coalescent theory (Hudson, 1989). The parameter (Watterson, 1975) was estimated for our sample with DnaSP version 3.14 (Rozas and Rozas, 1999) and used to solve H in equation 10 of Hudson (1989). Once the use of parsimony was validated, the haplotype network was estimated by finding haplotypeconnections with probabilities above the 0.95 limit. TCS version 1.0 (Clement et al, 2000) implements the algorithm in Templeton et al (1992). It was used to estimate those probabilities and for network construction.