We analyzed 17 microsatellites (AmpFLSTR® Yfiler® Life Technologies) using Genescan v. 3.7 (Applied Biosystems). Electrophoresis was performed in an ABBI 3130 analyzer. Our Laboratory has been granted the YHRD Quality Test Certificate for Yfiler analysis.
We built median-joining networks (Bandelt, Forster & Röhl 1999) using the NETWORK software (fluxus-engineering.com) with weights assigned following Muzzio et al. (2010) . We used Rho dating methods (Forster et al. 1996 , Saillard et al. 2000) , considering the rate based on population comparison 6.9 104/ STR/generation (Zhivotovsky et al. 2004 ). The geographic distribution of different sub-haplogroups was analyzed by building isofrequency graphs with Surfer12 software, applying the Kriging algorithm each time that the number of samples was higher than 20 (Golden Software 2002) .