a b s t r a c t honey bees collected from 32 different localities in greece were studied based on the geometric morphometrics approach using the coordinates of 19 landmarks located at wing vein intersections. procrustes analysis, principal component analysis, and canonical variate analysis (cVa) detected population variability among the studied samples. according to the principal component analysis (pca) of pooled data from each locality, the most differentiated populations were the populations from the aegean island localities astypalaia, chios, and kythira. however, the populations with the most distant according to the canonical variate analysis performed on all measurements were the populations from heraklion and chania (both from crete island). these results can be used as a starting point for the use of geometric morphometrics in the discrimination of honey bee populations in greece and the establishment of conservation areas for local honey bee populations.keywords: apis mellifera, geometric morphometrics, greece, morphological discrimination. Traditionally, intraspecific taxonomy of the honey bee Apis mellifera L. has been based on morphology. Twenty-seven subspecies of A. mellifera are currently recognized on the basis of morphometric characteristics (Ruttner, 1988(Ruttner, , 1992Sheppard et al., 1997;Engel, 1999;Sheppard and Meixner, 2003). These subspecies are also described as 'geographic races' because their distributions correspond to distinct geographic areas. (Garnery et al., 1993; Arias and Sheppard, 1996;Franck et al., 2000Franck et al., , 2001Miguel et al., 2007;Cánovas et al., 2008). Traditional morphometrics, so-called classical morphometry, was the only method to describe the diversity of honey bees for a long time (Ruttner, 1988), but biochemical methods have been developed (i.e., allozymes, isoenzymes) (Nunamaker and Wilson, 1982; Badino et al., 1988;Bouga et al., 2005b;Ivanova et al., 2012)
76Morphological discrimination of Greek honey bees Tarès et al., 1993), mitochondrial DNA (mtDNA) (Moritz et al., 1986;Smith et al., 1989;Smith, 1991;Hunt and Page, 1992;Garnery et al., 1993;Oldroyd et al., 1995; Arias and Sheppard, 1996;Pedersen, 1996;De la Rúa et al., 2000), and microsatellites (Estoup et al., 1993;Garnery et al., 1998), are also used to study genetic variability in honey bees. Several attempts have been made to simplify morphological determination, such as reducing the number of characters (Dedej and Nazzi, 1994;Cermak and Kaspar, 2000), using a single wing cell (Francoy et al., 2006), or using computer software designed specifically for morphometric measurements (Daly et al., 1982; Batra, 1988;Tofilski, 2004). Dupraw (1964Dupraw ( , 1965a was the first to use a set of quantitative characteristics for wing venation and correctly classified many subspecies. In 1978, Daly and Balling successfully applied this method to differentiate Africanized and European honey bees in South America. Later, Daly et al. (1982) began to use digital measurements to investigate the morphometrics of ...