Free Access
Volume 41, Number 6, November-December 2010
Page(s) 589 - 601
Published online 15 March 2010

© INRA/DIB-AGIB/EDP Sciences, 2010


Unraveling the structural complexity of Apis cerana (Fabricius, 1793) has been a historically continuous process. Having largely abandoned the typological approach of Maa (1953) and adopted the multivariate techniques of DuPraw (1964, 1965), Ruttner (1988) published the first comprehensive, multivariate statistical analyses of A. cerana, which proved seminal in stimulating further research. Subsequent authors used Ruttner’s interpretations of A. cerana as a new baseline and concentrated on morphoclusters derived from multivariate analyses on a microscale level (Muzaffar and Ahmad, 1989; Pesenko et al., 1989; Rinderer et al., 1989; Limbipichai, 1990; Otis and Hadisoesilo, 1990; Singh et al., 1990; Sulistianto, 1990; Szabo, 1990; Verma, 1990, 1992; Verma et al., 1989, 1994; Kwon and Huh, 1992; Ono, 1992; Damus, 1995; Hadisoesilo et al., 1995; Fuchs et al., 1996; Sylvester et al., 1998) as well as a more regional, mesoscale level (Peng et al., 1989; Yang, 1986, 2001; Diniz-Filho et al., 1993; Damus, 1995; Damus and Otis, 1997; Hepburn et al., 2001b; Tilde et al., 2000; Kuang et al., 2002; Radloff and Hepburn, 2002; Tan et al., 2002; Radloff et al., 2003, 2005a,b,c). Hepburn and Radloff (1996, 1998) designated statistically defined populations of honey bees under the new coinage of “morphoclusters”. This obviated the use of infraspecific names previously applied to the A. cerana group (Hepburn et al., 2001a), because the typological approach is fundamentally incompatible with a statistical one.

Here we report the results of a multivariate morphometric analysis of A. cerana across its full geographical range and identify the statistically defined morphoclusters and subcluster populations within them. We consider the relationships of morphoclusters to each other, their geophysical and ecological surroundings and present a new distribution map for A. cerana based on 2687 references (Hepburn and Hepburn, 2006). Comments on gene flow or degree of temporal reproductive isolation are suggested from swarming and migration data. Finally, we apply geographically based common names to the morphoclusters and subclusters as is the general practice in economic entomology and ornithology.


2.1. Honey bees

The worker honey bee specimens, used for the combined mainland and oceanic A. cerana analysis, derive from several different sources: (1) new material was collected in Japan, Korea and Russia (25 colonies); (2) raw databases of the Institut für Bienenkunde at Oberursel (Germany), Grahamstown (South Africa), Bogor (Indonesia) and Kunming (China); (3) raw databases of previously published results for India and Nepal (Singh et al., 1990; Singh and Verma, 1993; Verma et al., 1989, 1994; Hepburn et al., 2001b), for Afghanistan and Pakistan (Radloff et al., 2005a), for Philippines and Indonesia (Rinderer et al., 1989; Sulistianto, 1999; Tilde et al., 2000) (n = 964 colonies; Tab. I, supplementary material). Previous classification results of processed data for Japan from Akahira and Sakagami (1959a,b) and for Sri Lanka from Szabo (1990) and additional information based on statistically processed data for Thailand, Indonesia and Malaysia from Limbipichai (1990); Damus (1995); Damus and Otis (1997) and for Thailand from Sylvester et al. (1998) were incorporated into the final population structure and classification of A. cerana table (Radloff et al., 2005c). The geographical origins of A. cerana worker honey bees from 1067 colonies at 354 localities are listed in Table I (supplementary material) and shown in Figure 1.

thumbnail Figure 1

Geographical distribution of the six Apis cerana morphocluster groups and identifiable subspecific populations obtained in the current analysis.

2.2. Measurements

Twelve morphological characters of worker honey bees common to the Ruttner/Oberursel, Kunming and Grahamstown databases as well as those in the present study were utilized (n = 964 colonies). Their Ruttner (1988) and Verma et al. (1989, 1994) character numbers are given in round and square brackets, respectively, as follows: length of metafemur (5) [30], length of metatibia (6) [31], metatarsus length (7) [32], tergum 3, longitudinal (9) [42], tergum 4, longitudinal (10) [45], sternum 3, longitudinal (11) [49], wax plate of sternum 3, longitudinal (12) [48], sternum 6, longitudinal (15) [50], forewing length (17) [5], wing angle B4 (22) [10], wing angle D7 (23) [11] and wing angle G18 (25) [13].

2.3. Data analysis

Multivariate statistical analysis of the data (using 12 morphometric characters) included principal components analysis to identify possible morphoclusters and subclusters for A. cerana, hierarchical cluster analysis to identify homogeneous morphoclusters, linear discriminant analysis to determine the percentages of correct classification of colonies in each morphocluster, Wilks’ lambda statistic for testing significant differences between the multiple means of the characters entered into the discriminant functions and Levene’s F statistic procedure for testing heterogeneity of variances (Johnson and Wichern, 2002). All statistical analyses were performed using Statistica (StatSoft, 2007).

2.4. Definitions

The sympatric occurrence of A. cerana with other medium-sized bees, A. koschevnikovi, A. nigrocincta and A. nuluensis in Southeastern Asia unfortunately means that an undeterminable amount of some previous “A. cerana”literature may inadvertently include data derived from species other than A. cerana (Hepburn et al., 2001a). To assist in overcoming this problem, we list metric characters that in combination separate these four species of Apis even if they are pinned museum specimens as follows. The cubital index of the forewings are A.cerana = 3.9, A.koschevnikovi = 7.2, A.nigrocincta = 3.7, and A.nuluensis = 2.4, which quickly separates paired comparisons for all but an A. cerana and A. nigrocincta option. To further separate A. cerana from A. nigrocincta, the length of the basal portion of the radial cell of the forewing is 1.2 mm in A.cerana, 1.8 mm in A.nigrocincta; length of the apical portion of the radial cell is 1.8 mm in A.cerana, 1.1 mm in A. nigrocincta; the length of the labial palp is 1.8 mm in A. cerana, 3.7 mm in A. nigrocincta.


3.1. Morphoclusters

Principal components analysis using twelve morphometric characters of worker honey bees ((5), (6), (7), (9), (10), (11), (12), (15), (17), B4, D7, G18) was carried out on the morphometric mainland and oceanic databases of Oberursel, Kunming, Bogor, Grahamstown, Rinderer et al. (1989), Sulistianto (1990), and Tilde et al. (2000) comprising 964 colonies. Three principal components with eigenvalues greater than one were isolated. PC 1: size-related characters (5), (6), (7), (9), (10), (11), (12), (15) and (17) with component loadings between 0.76 and 0.95 accounted for 63.1% of the variation; PC 2: angles of venation (22) and (23) with component loadings 0.66 and 0.75, respectively, accounted for 10.4% of the variation; PC 3: angle of venation (25) with component loading 0.91 accounted for 8.4% of the variation. The three principal components accounted for a total of 81.9% of the variation in the data.

thumbnail Figure 2

Hierarchical clustering dendrogram for A. cerana, derived from complete linkage clustering on 12 morphometric characters (Ruttner numbers, (5), (6), (7), (9), (10), (11), (12), (15), (17), B4 (22), D7 (23) and G18 (25)) averaged for regional groups and countries.

PC plots using both the first and second PC scores and the first and third PC scores did not reveal distinct morphoclusters. Further substructuring of the PC plots by introducing local labelling and running a hierarchical cluster analysis using the mean scores for PC 1 to 3 to identify homogeneous morphoclusters revealed six main overlapping morphoclusters and are defined as follows (Fig. 2):

(1) Morphocluster I, here named “Northern cerana”, extends from northern Afghanistan and Pakistan through northwest India, across southern Tibet, northern Myanmar, China and then northeasterly into Korea, far eastern Russia and Japan (Fig. 2). Six subclusters or populations are morphometrically discernible within this morphocluster: (a) an “Indus”group in Afghanistan, Pakistan and Kashmir; (b) a “Himachali”group in Himachal Pradesh, India; (c) an “Aba”group in southern Ganshu and central and northern Sichuan provinces in China, northern China and Russia (larger bees); (d) a subcluster in central and eastern China; (e) a “southern”cerana subcluster in southern Yunnan, Guangdong, Guangxi and Hainan in China; and (f) a “Japonica” group in Japan and Korea (Fig. 1). Morphocluster I bees have been previously named as follows: A. skorikovi, A. c. abansis, A. c. abanensis, A. c. bijjieca, A. c. cathayca, A. c. cerana, A. c. fantsun, A. c. hainana, A. c. hainanensis, A. c. heimifeng, A. c. indica, A. c. japonica, A. c. javana, A. c. kweiyanga, A. c. maerkang, A. c. pekinga, A. c. peroni, A. c. skorikovi, A. c. shankianga and A. c. twolareca.

(2) Morphocluster II, here named “Himalayan cerana”includes the bees of northern India: (a) northwest, (c) northeast; and some of southern (c) Tibet and Nepal. Two subclusters are discernible within this morphocluster, which we term the bees of the northwest the “Hills”group, and those of the northeast, the “Ganges”group (Figs. 1 and 2). Morphocluster II bees have previously been named A. c. skorikovi, A. c. indica, A. c. himalayana, and A. c. himalaya.

(3) Morphocluster III, here named “Indian Plains cerana”occurs across the plains of central and southern India and Sri Lanka as a fairly uniform population, long known as “plains cerana”for this subcontinent (Figs. 1 and 2). Morphocluster III bees have only previously been termed A. c. indica.

(4) Morphocluster IV, here named “Indo-chinese cerana”form a compact group in Myanmar, northern Thailand, Laos, Cambodia and more southern Vietnam (Figs. 1 and 2). Morphocluster IV bees have been previously named A. c. indica and A. c. javana.

(5) Morphocluster V, here named “Philippinecerana”is restricted to the Philippines, but excluding most of Palawan Island which group with Morphocluster VI. Within these islands there are subclusters and we term these bees respectively after the major island groups there: “Luzon”bees, “Mindanao”bees, and “Visayas”bees. The latter two subclusters showed closer morphometric similarity than the former (Figs. 1 and 2). Morphocluster V bees have been previously named A. philippina, A. c. philippina and A. c. samarensis.

(6) Morphocluster VI, here named “Indo-Malayan cerana”, extend from southern Thailand, through Malaysia and Indonesia. This large area consists of a morphometrically rather uniform bee below the South China Sea. Three subclusters are discernible within this morphocluster: (a) Philippines (Palawan), Malaysia (North Borneo), Indonesia (Kalimantan) bees; (b) Malay Peninsula, Sumatera, and some Sulawesi bees; (c) Indonesia (Java, Bali, Irian Jaya, some Sulawesi and Sumatera) bees (Figs. 1 and 2). Morphocluster VI bees have been previously termed A. cerana, A. indica, A. javana, A. c. johni, A. lieftincki, A. peroni, A. vechti linda and A. v. vechti.

thumbnail Figure 3

Discriminant analysis plot of six morphoclusters of A. cerana using colony means of 12 common morphometric characters. Numbers indicate colonies from morphoclusters: 1 = colonies from northern Afghanistan, Pakistan, Kashmir, Himachal Pradesh, northern Myanmar and northern Vietnam, China, Russia, Japan and Korea; 2 = colonies from northwestern and northeastern India, some of southern Tibet and Nepal; 3 = colonies from central and southern India and Sri Lanka; 4 = colonies from Myanmar, Thailand, Cambodia and central and southern Vietnam; 5 = colonies from Philippines (Luzon, Mindanao and Visayas); 6 = colonies from Philippines (Palawan), Malaysia and Indonesia.

The results of a linear discriminant analysis using 12 morphometric characters classified: 95.3% of colonies (n = 342) from Morphocluster I correctly and misclassified 17 colonies; 94.5% of colonies (n = 155) from Morphocluster II correctly and misclassified 9 colonies; 95.1% of colonies (n = 77) from Morphocluster III correctly and misclassified 4 colonies; 95.3% of colonies (n = 121) from Morphocluster IV correctly and misclassified 6 colonies; 89.5% of colonies (n = 68) from Morphocluster V correctly and misclassified 8 colonies; and 80.4% of colonies (n = 123) from Morphocluster VI correctly and misclassified 30 colonies. Overall 92.3% of the colonies were correctly classified into the morphclusters. A significant difference between the means of the six morphoclusters was established (Wilks’ test: Λ = 0.0089, F60,4428, P < 0.0001). Canonical 1 and 2 scores plot showed the morphoclusters are not disjoint (Fig. 3, confidence ellipses at 95%).

3.2. Taxonomy of A. cerana

The ICZN summary of the taxonomy for A. cerana is summarized here, with the formal synonymy of subspecific and other names noted as required nomenclaturally:

Apis cerana Fabricius

Apis cerana Fabricius, 1793: 327.

Apis indica Fabricius, 1798: 274.

Apis socialis Latreille, 1804a: 390.

Apis peroni Latreille, 1804b: 173.

Apis gronovii Guillou, 1841: 323.

Apis perrottetii Guérin-Méneville, 1844: 461.

Apis delessertii Guérin-Méneville, 1844: 461.

Apis sinensis Smith, 1865: 380.

Apis mellifica variety japonica Radoszkowski, 1887: 436.

Apis delesserti Buttel-Reepen, 1906: 168. Nomen vanum (unjustified emendation) [ICZN Article 33.2.3].

Apis indica variety javana Enderlein, 1906: 337.

Apis indrea Baldensperger, 1928: 173. Lapsus calami.

Apis johni Skorikov, 1929: 251.

Apis indica philippina Skorikov, 1929: 252.

Apis indica skorikovi Maa, 1944: 4. Nomen nudum (vide skorikovi Engel, 1999). [ICZN Article 13.1].

Apis mellifera gandhiana Muttoo, 1951: 153.

Apis (Sigmatapis) lieftincki Maa, 1953: 572.

Apis (Sigmatapis) samarensis Maa, 1953: 580.

Apis indica sinensis ussuriensis Goetze, 1964: 26. Nomen nudum et invalidum [ICZN Articles 10.2, 13.1, 15.2].

Apis cerana abansis Yang & Kuang, In Yang, Xu, Kuang, & Liu, 1986: 90. Syn. nov.

Apis cerana hainana Yang & Kuang, In Yang, Xu, Kuang, & Liu, 1986: 91. Syn. nov.

Apis cerana fantsun mountaina Zhuang, 1989: 61. Nomen invalidum [ICZN Articles 10.2, 15.2]. Syn. nov.

Apis cerana shenp countica Zhuang, 1993 [sic]: 61. Lapsus calami.

Apis ceran [sic] cathayca Zhuang, 1989: 62. Syn. nov.

Apis ceran [sic] fantsan mountaina Zhuang, 1993 [sic]: 62. Lapsus calami.

Apis ceran [sic] sheng countica Zhuang, 1989: 62. Nomen invalidum [ICZN Articles 10.2, 15.2]. Syn. nov.

Apis cerana pekinga Zhuang, 1989: 62. Syn. nov.

Apis cerana bijieca Zhuang, 1989: 62. Syn. nov.

Apis cerana kweiyanga Zhuang, 1989: 62. Syn. nov.

Apis cerana twolakeca Zhuang, 1989: 62. Syn. nov.

Apis cerana shankianga Zhuang, 1989: 62. Syn. nov.

Apis cerana himalaya Smith, 1991: 154. Nomen nudum (erroneous and invalid Latinization of informal term “Himalaya race”as applied by Ruttner, 1988: 121) [ICZN Article 13.1].

Apis cerena Willis, Winston, & Honda, 1992: 169. Lapsus calami.

Apis cerana himalayana Ruttner, 1992: 296. Nomen nudum (misspelling of already invalid name himalaya Smith, 1991). [ICZN Article 13.1]

Apis cerana twolareca Chen, 1993: 21. Lapsus calami.

Apis cerana heimifeng Engel, 1999: 179. Syn. nov.

Apis cerana skorikovi Engel, 1999: 180. Syn. nov.

Apis cerana abaensis Hepburn, Smith, Radloff & Otis, 2001: 6. Nomen nudum et invalidum, et nomen vanum [ICZN Articles 13.1, 16.1, 16.4, 33.2.3]. Syn. nov.

Apis cerana hainanensis Hepburn, Smith, Radloff & Otis, 2001: 6. Nomen nudum et invalidum, et nomen vanum [ICZN Articles 13.1, 16.1, 16.4, 33.2.3]. Syn. nov.


4.1. Trinomial names

Trinomial epithets (e.g. Apis cerana cerana) represent a historical artefact of earlier nomenclatural systems and can create considerable confusion when applied in modern systems. While prior to 1969 such names were usually applied for extreme forms of variation, more modern applications of trinomial names have been for ranges of variation (and not necessarily a single variant) over a specified geographic region. Thus, while the names can be used to refer to a circumscribed cluster within a variable species [as they are meant under rules of the ICZN, 4th Edition, 1999], often they are misapplied and collapse back into older systems of use. Rather than revising and updating the system of names, subsequent authors perpetuated this historical artefact and by the end of the 20th century, some 31 putative “taxa”of A. cerana had been established (Hepburn et al., 2001a), although many were not valid under modern ICZN rules (e.g., briefly summarized by Engel, 1999). The question must be asked: how did these “names”come about and are any of them actually biologically meaningful? The short answer is that all of these trinomials arose as artefacts of sampling coupled with the use by honey bee researchers of long out-dated taxonomic procedures.

According to modern taxonomic practice none of the “formal”Latin names historically applied to variations in Apis would be considered suitable, and many are not even legitimate. As such, those that were validly established under the nomenclatural rules of their day should be accordingly recognized as synonyms. Similarly, the remaining, largely post-1970, names also should be considered synonyms, with the caveat for several that they were invalidly established under any taxonomic concept. Detailed synonymic summaries in standard taxonomic format for A. cerana are given in Engel (1999). As a consequence of these published synonymies, trinomials such as A. c. japonica or A. c. cerana no longer have any official, nomenclatural standing in Apis classification and under the rules of the ICZN. All former subspecies simply fall into disuse. By the same token, the morphocluster names such as “Himalayan cerana”or “Indo-Malayan cerana”, as we propose are not binding, but offer Apis biologists an alternative solution in which the contiguous populations of A. cerana can be sensibly separated in what is still a biologically meaningful way.

4.2. Constraints on A. cerana morphoclusters

The interactive effects of limited sampling coupled with the Linnean system of nomenclature have shaped the historical process of unravelling the A. cerana complex. The importance of sample size, sample characters and inter-locality sampling distances has now been resolved (Radloff and Hepburn, 1998). In a series of regional parallel studies on the A. cerana of China, Tan Ken et al. (2002, 2003) showed that bees from the northern high-altitude areas of Yunnan Province were clearly larger and darker, and showed similarities with samples from Beijing, Nepal and north India, whereas bees from southern Yunnan clustered with the bees from Thailand and Vietnam. These results are completely consistent with those of Radloff et al. (2005b) for the bees of southern Yunnan. Morphometric analyses of Apis cerana from oceanic Asia yielded two distinct morphoclusters: (1) the bees of Japan and (2) the bees of all the other south Asian islands bringing the then total number of morphoclusters to seven (Radloff et al., 2005c). On completion of the above series of regional mesoscale studies, the newly formed comprehensive data set for all A. cerana was subjected to multivariate morphometric analysis. The final result is that six morphoclusters of A. cerana were obtained (Fig. 1).

When all of the mesoscale morphoclusters of the present authors are compared with the new macroscale results, the only discrepancies are that in the former (1) the bees of the Philippines were included with those of Indonesia and Borneo and (2) the bees of Japan are now placed in the Northern Asia morphocluster of the latter. However there are differences between the mapped morphocluster results of Ruttner (1988), Damus and Otis (1997) and our present ones. These discrepancies are best explained by sampling differences in each study, which affected the degree of morphometric discrimination of the honey bees of Japan. Ruttner (1988) only had access to a very small sample of large A. cerana from China and none from Russia. The only Morphocluster I bees available to him were from the far northwest of the A. cerana range (Afghanistan and Pakistan) and some 6000 km distant from Japan, the bees of which form a subcluster in a continuum of A.cerana Morphocluster I. Gaps in sampling inevitably resulted in the differences between Afghani and Japanese A. cerana being statistically magnified. The dataset of Damus and Otis (1997) was based on the much smaller bees of the more southerly oceanic islands (Philippines, Indonesia, Borneo, etc.) with the same effects. The data of Tan et al. (2008) is entirely consistent with that presented here.

4.3. Biogeography

A. cerana naturally occurs on an Asian landscape of some 30 M km2 encompassing a series of climatic zones from tropical moist rainforest, wet-dry tropical savanna, mid-latitude steppe, dry mid-latitude grasslands, moist continental deciduous forest and taiga (Köppen Climate Classification System – cf. Müller, 1982). It is of interest to note how the distributions of the different morphoclusters and, possibly of subclusters within them, may be related to these various physiographic and climatic factors. The Philippine V and Indo-Malayan VI Morphoclusters both lay below the South China Sea in the tropical wet climatic zone, yet are morphometrically differentiated. Reproductive isolation among them can be interpreted in a geological context (Voris, 2000). Palawan bees are in morphocluster VI and link with bees from Borneo which is not surprising because the sea depth between Borneo and Palawan is much shallower than between Palawan and the other 3 Philippine islands (Voris, 2000).

During the last 250 KBP of the Pleistocene, the Philippines were physically separated from the Indo-Malayan region about 75% of the time, clearly allowing differentiation through genetic drift. Differentiation into biometric subclusters within the Philippine group among the major islands systems reflects that they too were physically separated for 70% of the 250 KBP. In stark contrast, the bees of the Indo-Malayan morphocluster occur in the Indonesian chain of islands on islands that have been connected 50% of the last 250 KBP and nearly 70% of the last 17 K (Voris, 2000), which clearly forms a genetically more panmictic population. These geological data are completely consistent with the Euclidean linkage distances between these two oceanic morphoclusters.

Proceeding to the continental mainland, land links between the Indo-Malayan (VI) and Southeast Asian (IV) Morphoclusters resulted from lowered sea levels nearly 50% of the last 250 KBP and nearly 70% of the last 17 KBP, but there were never land connections between them and the Indian subcontinent. Here the Euclidean linkage distances provide a trend in the right direction. Nonetheless, there is a sharp transition between Indo-Malayan VI and the Southeast Asian Morphocluster IV, from a tropical dipterocarp rain forest to mixed moist deciduous forest with a wet and dry tropical climate to the northeast of the Isthmus of Kra in the northern Thai Malay Peninsula (Hughes et al., 2003). The rain forests actually extend up the western side of the isthmus along the Chin Hills-Arakan Yoma forests running between Myanmar and Bangladesh, the Naga Hills of India, the highlands of northern Laos and Chuor Phnum of Cambodia. These ecological discontinuities probably provide physiographic and ecological boundaries between the Southeast Asian IV and the adjoining eastern edge of the sub-Himalayan Morphocluster II.

The sub-Himalayan (II) and Southeast Asian (IV) Morphoclusters are adjacent in terms of their Euclidean linkage distances and subtend northern Morphocluster I along the Himalayas and then give way to a complex range of mountains (Hengduanshan, Wuliangshan, Taniantawengshan, Nushan and Qingshuilangshan), which turn southwards along the present Myanmar-China border. The sub-Himalayan morphocluster contains two recognizable subclusters, the larger Indus and smaller Ganges populations. The morphometric distinctness of these two populations arises from the fact that they differ significantly in the seasonal timing of reproductive swarming (Hepburn et al., 2001b).

Northern Morphocluster I extends over the largest land areas where A. cerana occurs (Fig. 1) and consists of six biometrically definable subclusters. The Indus population (Morphocluster I) occurs in a tropical wet and dry climate while its northwestern neighbour, the Himachali population occurs in the subtropical steppe. These climatic differences provide a degree of reproductive isolation between the populations because of seasonal differences in swarming so reducing the possibility of gene flow between them (Muzaffar and Ahmad, 1989; Saraf, 1985; Sharma, 1960). The Himachali population is also largely isolated from the eastern (north and south) Chinese subclusters by the intervening subarctic climate of the Tibetan highlands. Movingeastwards, the next populations of Morphocluster I are the northerly Aba and central and eastern China and Japonica biometric subclusters. The Aba area is in a subtropical mountain system while that of the central and eastern China and Japonica subclusters are situated in subtropical humid forest. The Aba and central and eastern subclusters have only partial overlap in seasonal swarming (Gong, 1983; Liu and Lai, 1990; Zhu, 1996). The South China subcluster (S Yunnan, Hainan, Guangdong, Guangxi) is essentially tropical.

As a final observation, it can be noted that A. cerana has been introduced to areas where it did not occur. In the 19th century it was introduced into the northeastern provinces of Heilongjiang and Jilin, China from maritime Russia (Chen, 1993). Likewise it has been introduced to Ambon Island to the east of Sulawesi (Maa, 1953) as well as Iran (Gassparian, 1977). That A. cerana is an invasive species is evident from the fact that it was introduced from Java into Papua New Guinea in the late 1970’s and now occurs throughout New Guinea (Anderson, 1989). From there it has reached islands in the Torres Strait (Dunn, 1992) and even more recently New Britain and the Solomon Islands (Anderson, 2005 – unpublish. data). A. cerana from Papua New Guinea are morphometrically indistinguishable from those of Java (Radloff et al., 2005c). It has also been intercepted several times at Darwin, Brisbane and near Perth, Australia.

Online material

Download PDF file.

Table I

Geographical origins, altitude and colony sample sizes (n) of A. cerana used in the morphometric analysis of honeybees of Asia.


  • Akahira Y., Sakagami S.F. (1959) A biometrical study on the Japanese honey bee, observations upon some populations of Kyûshû (studies on the Japanese honey bee, Apis cerana cerana Fabr. IX), J. Hokkaido Gakugei Univ. 2, 353–362. [Google Scholar]
  • Akahira Y., Sakagami S.F. (1959) Observations on the variability of wing venation in the honeybees, J. Fac. Sci. Hokkaido Univ. VI 14, 175–184. [Google Scholar]
  • Anderson D.L. (1989) Mites and pathogens of the European honey bee Apis mellifera and the eastern hive bee Apis cerana in Papua New Guinea, DSIR, Entomology Division, Auckland. [Google Scholar]
  • Baldensperger P.J. (1928) Bees in their natural state and with beekeepers, Bee World 9, 173–174. [Google Scholar]
  • Buttel-Reepen H. von (1906) Apistica, Beiträge zur Systematik, Biologie, sowie zur geschichtlichen und geographischen Verbreitung der Honigbiene (Apis mellifica L.), ihrer Varietäten und der überigen Apis-Arten, Mitt. Zool. Mus. Berlin 3, 117–201. [Google Scholar]
  • Chen Y.-C. (1993) Apiculture in China, China Agricultural Publication, Beijing [in Chinese]. [Google Scholar]
  • Damus M.S. (1995) A morphometric and genetic analysis of honey bee (Apis cerana F.) samples from Malaysia: population discrimination and relationships, Thesis, University of Guelph, Guelph, Ontario. [Google Scholar]
  • Damus M.S., Otis G.W. (1997) A morphometric analysis of Apis cerana F and Apis nigrocincta Smith populations from Southeast Asia, Apidologie 28, 309–323. [CrossRef] [EDP Sciences] [Google Scholar]
  • Diniz-Filho J.A.F., Malapsina O., Pignata M.I.B. (1993) Geographic variation in Apis cerana indica F.: a spatial autocorrelation analysis of morphometric patterns, J. Apic. Res. 32, 65–72. [Google Scholar]
  • Dunn K.J. (1992) Exotic Asian bee detected in Torres Strait, Bee Briefs 9, 18–19. [Google Scholar]
  • DuPraw E.J. (1964) Non-Linnean taxonomy, Nature 202, 849-852. [CrossRef] [Google Scholar]
  • DuPraw E.J. (1965) Non-Linnean taxonomy and the systematics of honeybees, Syst. Zool. 14, 1–24. [CrossRef] [Google Scholar]
  • Enderlein G. (1906) Neue Honigbienen und Beiträge zur Kentnnis der Verbreitung der Gattung Apis, Stett. Entomol. Ztg. 67, 331–344. [Google Scholar]
  • Engel M.S. (1999) The taxonomy of recent and fossil honey bees (Hymenoptera: Apidae, Apis), J. Hymenoptera Res. 8, 165–196. [Google Scholar]
  • Fabricius J.C. (1793) Entomologia Systematica Emendata et Aucta, Secundum Classes, Ordines, Genera, Species adiectis Synonymis, Locis, Observationibus, Descriptionibus [Vol. 2], Proft, Hafniae [Copenhagen]. [Google Scholar]
  • Fabricius J.C. (1798) Supplementum Entomologiae Systematicae, Proft, Hafniae [Copenhagen]. [Google Scholar]
  • Fuchs S., Koeniger N., Tingek S. (1996) The morphometric position of Apis nuluensis Tingek, Koeniger and Koeniger, 1996 within cavity-nesting honey bees, Apidologie 27, 397–405. [CrossRef] [EDP Sciences] [Google Scholar]
  • Gassparian S. (1977) Studies on Apis indica cerana in eastern part of Iran, in: Proc. XXVI Int. Apic. Congr., Adelaide, Australia, pp. 93–96. [Google Scholar]
  • Goetze G.K.L. (1964) Die Honigbiene in natürlicher und künstlicher Zuchtauslese, Teil 1: Systematik, Zeugung und Vererbung, Monogr. Angew. Entomol. 19, 1–120. [Google Scholar]
  • Gong Y.F. (1983) The natural beekeeping conditions and honeybee races in China, J. Fujian Agric. Coll. 12, 241–249. [Google Scholar]
  • Guérin-Méneville F.E. (1844) Iconographie du Règne Animal de G. Cuvier, ou Représentation d’après Nature de l’une des Espèces les plus Remarquables et Souvent non Encore Figurées, de Chaque Genre d’Animaux [Vol. 3], Baillière, Paris. [Google Scholar]
  • Guillou E.J.F. le (1841) Des Insectes Hyménoptères recueillis dans le Voyage de circumnavigation des Corvettes l’Astrolabe et la Zélée, Ann. Soc. Entomol. F. 10, 311–324. [Google Scholar]
  • Hadisoesilo S., Otis G.W., Meixner M. (1995) Two distinct populations of cavity-nesting honeybees (Hymenoptera: Apidae) in south Sulawesi, Indonesia, J. Kans. Entomol. Soc. 68, 399–407. [Google Scholar]
  • Hepburn H.R., Hepburn C. (2006) Bibliography of Apis cerana, Apidologie 37, 651–652. [CrossRef] [EDP Sciences] [Google Scholar]
  • Hepburn H.R., Radloff S.E. (1996) Morphometric and pheromonal analyses of Apis mellifera L. along a transect from the Sahara to the Pyrenees, Apidologie 27, 35–45. [CrossRef] [EDP Sciences] [Google Scholar]
  • Hepburn H.R., Radloff S.E. (1998) Honeybees of Africa, Springer-Verlag, Berlin. [Google Scholar]
  • Hepburn H.R., Smith D.R., Radloff S.E., Otis G.W. (2001) Infraspecific categories of Apis cerana: morphometric, allozymal and mtDNA diversity, Apidologie 32, 3–23. [CrossRef] [EDP Sciences] [Google Scholar]
  • Hepburn H.R., Radloff S.E., Verma S., Verma L.R. (2001) Morphometric analysis of Apis cerana populations in the southern Himalayan region, Apidologie 32, 435–447. [CrossRef] [EDP Sciences] [Google Scholar]
  • Hughes J.B., Round P.D., Woodruff D.S. (2003) The Indochinese-Sundaic faunal transition at the Isthmus of Kra: an analysis of resident forest bird species distributions, J. Biogeogr. 30, 569–580. [CrossRef] [Google Scholar]
  • ICZN (1999) International Code of Zoological Nomenclature [4th Edition], International Trust for Zoological Nomenclature, London. [Google Scholar]
  • Johnson R.A., Wichern D.W. (2002) Applied Multivariate Statistical Analysis, 5th ed., Prentice Hall, Upper Saddle River, New Jersey. [Google Scholar]
  • Kuang B.Y., Tan K., Liu Y.Q. (2002) Honeybee biology, Yunnan Science Publishing House, Kunming [in Chinese]. [Google Scholar]
  • Kwon Y.J., Huh E.Y. (1992) Electronmorphometric classification of the native honey bees from Korea. Part I. Discriminant analysis on the wing characters, Korean J. Apic. 7, 1–31. [Google Scholar]
  • Latreille P.A. (1804) Mémoire sur un gâteau de ruche d’une abeille des Grandes-Indes, et sur différences des abeilles proprement dites, vivant en grande société, de l’ancien continent et du nouveau, Ann. Mus. Hist. Nat. 4, 383–394. [Google Scholar]
  • Latreille P.A. (1804) Notice des espèces d’abeilles vivant en grande société, ou abeilles proprement dites, et description d’espèces nouvelles, Ann. Mus. Hist. Nat. 5, 161–178. [Google Scholar]
  • Limbipichai K. (1990) Morphometric studies on the eastern honey bee (Apis cerana Fabricius) in Thailand and the Malaysian peninsula, Thesis, Chulalongkorn University, Bangkok. [Google Scholar]
  • Liu C.S., Lai Y.S. (1990) Swarming in Apis cerana cerana, Apic. Sci. Tech. 3, 2–4. [Google Scholar]
  • Maa T.C. (1944) On the classification and phylogeny of the Chinese honeybees, Entomol. Shaowuana 1, 4–5. [Google Scholar]
  • Maa T.C. (1953) An inquiry into the systematics of the tribus Apidini or honeybees (Hymenoptera), Treubia 21, 525–640. [Google Scholar]
  • Müller M.J. (1982) Selected climatic data for a global set of standard stations for vegetation science, Junk, The Hague. [Google Scholar]
  • Muttoo R.N. (1951) The correct scientific nomenclature for our Indian hive bees, Indian Bee J. 13, 150–153. [Google Scholar]
  • Muzaffar N., Ahmad R. (1989) Distribution and competition of Apis spp. in Pakistan, in: Proc. 4th Int. Conf. Apic. Trop. Clim., Cairo, Egypt, pp. 449–452. [Google Scholar]
  • Ono M. (1992) The Asian honey bees (Apis spp.), Honeybee Sci. 13, 19–22. [Google Scholar]
  • Otis G.W., Hadisoesilo S. (1990) Honey bee survey of south Sulawesi, J. Penelitian Kehutanan 4, 1–3. [Google Scholar]
  • Peng Y.S., Nasr M.E., Locke S.J. (1989) Geographical races of Apis cerana Fabricius in China and their distribution, Review of recent Chinese publications and a preliminary statistical analysis, Apidologie 20, 9–20. [CrossRef] [EDP Sciences] [Google Scholar]
  • Pesenko Y.A., Lelei A.S., Radchenko V.G., Filatkin G.N. (1989) Chinese wax-bee, Apis cerana cerana F. (Hymenoptera, Apidae) in the Far East of the USSR, Ent. Obozr. 68, 527–548. [Google Scholar]
  • Radloff S.E., Hepburn H.R. (1998) The matter of sampling distance and confidence levels in the subspecific classification of honeybees, Apis mellifera L., Apidologie 29, 491–501. [CrossRef] [EDP Sciences] [Google Scholar]
  • Radloff S.E., Hepburn H.R. (2002) Multivariate morphometric analysis of the Apis cerana populations of southern Asia, in: Proc. 14th Int. Congr., IUSSI, Sapporo, Japan, p. 209. [Google Scholar]
  • Radloff S.E., Hepburn H.R., Bangay L. (2003) Quantitative analysis of intercolonial and intracolonial morphometric variance in honeybees, Apis mellifera and Apis cerana, Apidologie 34, 339–351. [CrossRef] [EDP Sciences] [Google Scholar]
  • Radloff S.E., Hepburn H.R., Fuchs S. (2005) The morphometric affinities of Apis cerana of the Hindu Kush and Himalayan regions of western Asia, Apidologie 36, 25–30. [CrossRef] [EDP Sciences] [Google Scholar]
  • Radloff S.E., Hepburn H.R., Hepburn C., Fuchs S., Otis G.W., Sein M.M., Aung H.L., Pham H.T., Tam D.Q., Nuru A.M., Tan K. (2005) Multivariate morphometric analysis of Apis cerana of southern mainland Asia, Apidologie 36, 127–139. [CrossRef] [EDP Sciences] [Google Scholar]
  • Radloff S.E., Hepburn H.R., Fuchs S., Otis G.W., Hadisoesilo S., Hepburn C., Tan K. (2005) Multivariate morphometric analysis of the Apis cerana populations of oceanic Asia, Apidologie 36, 475–492. [CrossRef] [EDP Sciences] [Google Scholar]
  • Radoszkowski O. (1887) Hyménoptères de Korée, Horae Soc. Entomol. Ross. 21, 428–436. [Google Scholar]
  • Rinderer T.E., Koeniger G., Tingek S., Mardan M.B. (1989) A morphological comparison of the cavity dwelling honeybees of Borneo Apis koschevnikovi (ButtelReepen, 1906) and Apis cerana (Fabricius, 1793), Apidologie 20, 405–411. [CrossRef] [EDP Sciences] [Google Scholar]
  • Ruttner F. (1988) Biogeography and taxonomy of honeybees, Springer-Verlag, Berlin. [Google Scholar]
  • Ruttner F. (1992) Naturgeschichte der Honigbienen, Ehrenwirth, Munich. [Google Scholar]
  • Saraf S.K. (1985) Beekeeping in Gurez, Jammu and Kashmir, India, Indian Bee J. 47, 36–37. [Google Scholar]
  • Sharma P.L. (1960) Observations on the swarming and mating habits of the Indian honeybee, Bee World 41, 121–125. [Google Scholar]
  • Singh M.P., Verma L.R. (1993) Morphometric comparison of three geographic populations of the northeast Himalayan Apis cerana, Asian apiculture, in: Connor L.J., Rinderer T.E., Sylvester H.A., Wongsiri S. (Eds.), Wicwas Press, Cheshire, Connecticut, pp. 67–80. [Google Scholar]
  • Singh M.P., Verma L.R., Daly H.V. (1990) Morphometric analysis of the Indian honeybee in the northeast Himalayan Region, J. Apic. Res. 29, 3–14. [Google Scholar]
  • Skorikov A.S. (1929) Eine neue Basis für eine Revision der Gattung Apis L., Rep. Appl. Entomol., Leningrad 4, 249–270. [Google Scholar]
  • Smith D.R. (1991) Mitochondrial DNA and honey bee biogeography, in: Smith D.R. (Ed.), Diversity in the Genus Apis, Westview Press, Boulder, pp. 131–176. [Google Scholar]
  • Smith F. (1865) On the species and varieties of the honey-bees belonging to the genus Apis, Ann. Mag. Nat. Hist., Ser. 3, 15, 372–380. [Google Scholar]
  • StatSoft Inc. (2007) STATISTICA, version 8.1, [Google Scholar]
  • Sulistianto A. (1990) Morphometric analysis of Indonesian honeybees, Thesis, University of Wales College of Cardiff, Cardiff. [Google Scholar]
  • Sylvester H.A., Limbipichai K., Wongsiri S., Rinderer T.E., Mardan M.B. (1998) Morphometric studies of Apis cerana in Thailand and the Malaysian peninsula, J. Apic. Res. 37, 137–145. [Google Scholar]
  • Szabo T.I. (1990) Morphometric characteristics of Apis cerana from Sri Lanka, Apidologie 21, 505–509. [CrossRef] [EDP Sciences] [Google Scholar]
  • Tan K., Fuchs S., Koeniger N., Ruiguang Z. (2002) Morphological characterization of Apis cerana (Fabr.) in the Province Yunnan of China, Apidologie 34, 553–561. [Google Scholar]
  • Tan K., Hepburn H.R., Radloff S.E., Fuchs S., Fan X., Zhang L., Yang M. (2008) Multivariate morphometric analysis of the Apis cerana of China, Apidologie 39, 343–353. [CrossRef] [EDP Sciences] [Google Scholar]
  • Tan K., Zhang X., He S., Zhou D.Y. (2003) Morphometrics and biogeography of Apis cerana in Yunnan, Apic. China 53, 4–6. [Google Scholar]
  • Tilde A.C., Fuchs S., Koeniger N., Cervancia C.R. (2000) Morphometric diversity of Apis cerana Fabr. within the Philippines, Apidologie 31, 249–264. [CrossRef] [EDP Sciences] [Google Scholar]
  • Verma L.R. (1990) Genetic diversity and conservation of Apis cerana in Hindu Kush-Himalayas, in: Proc. 11th Int. Congr., IUSSI, Bangalore, India, pp. 118–119. [Google Scholar]
  • Verma L.R. (1992) Species and genetic diversity in Himalayan honeybees, Honeybees in mountain agriculture, in: Verma L.R. (Ed.), Oxford and IBH Publishing Co. Pvt. Ltd., New Delhi, India, pp. 39–49. [Google Scholar]
  • Verma L.R., Kafle G.P., Sharma A., Mattu V.K. (1989) Biometry of Apis cerana of Nepal Himalayas, in: Proc. 4th Int. Conf. Apic. Trop. Clim., Cairo, Egypt, pp. 458–465. [Google Scholar]
  • Verma L.R., Mattu V.K., Daly H.V. (1994) Multivariate morphometrics of the Indian honeybee in the northwest Himalayan Region, Apidologie 25, 203–223. [CrossRef] [EDP Sciences] [Google Scholar]
  • Voris H.K. (2000) Maps of Pleistocene sea levels in Southeast Asia: shorelines, river systems and time durations, J. Biogeogr. 27, 1153–1167. [CrossRef] [Google Scholar]
  • Willis L.G., Winston M.L., Honda B.M. (1992) Phylogenetic relationships in the honeybee (genus Apis) as determined by the sequence of the cytochrome oxidase II region of mitochondrial DNA, Mol. Phylogenet. Evol. 1, 169–178. [CrossRef] [PubMed] [Google Scholar]
  • Yang G.H. (1986) The distribution of the eastern honeybee (Apis cerana Fab.) in China and an analysis of its subspecies, J. Yunnan Agric. Univ. 12, 202–206. [Google Scholar]
  • Yang G.H. (2001) The biology of Apis cerana, G.H. Yang, Beijing [in Chinese]. [Google Scholar]
  • Yang G., Xu S., Kuang B., Liu Y. (1986) The distribution of the genus [sic] Apis cerana Fab. in China and some of its subspecies, J. Yunnan Agric. Univ. 1, 89–92. [Google Scholar]
  • Zhu Y.Y. (1996) The causes of swarming in Apis cerana cerana, J. Bee (China) 10, 23 [in Chinese]. [Google Scholar]
  • Zhuang D. (1989) New subspecies of Apis cerana, Southwest China J. Agric. Sci. 2, 61–65. [Google Scholar]

All Tables

Table I

Geographical origins, altitude and colony sample sizes (n) of A. cerana used in the morphometric analysis of honeybees of Asia.

All Figures

thumbnail Figure 1

Geographical distribution of the six Apis cerana morphocluster groups and identifiable subspecific populations obtained in the current analysis.

In the text
thumbnail Figure 2

Hierarchical clustering dendrogram for A. cerana, derived from complete linkage clustering on 12 morphometric characters (Ruttner numbers, (5), (6), (7), (9), (10), (11), (12), (15), (17), B4 (22), D7 (23) and G18 (25)) averaged for regional groups and countries.

In the text
thumbnail Figure 3

Discriminant analysis plot of six morphoclusters of A. cerana using colony means of 12 common morphometric characters. Numbers indicate colonies from morphoclusters: 1 = colonies from northern Afghanistan, Pakistan, Kashmir, Himachal Pradesh, northern Myanmar and northern Vietnam, China, Russia, Japan and Korea; 2 = colonies from northwestern and northeastern India, some of southern Tibet and Nepal; 3 = colonies from central and southern India and Sri Lanka; 4 = colonies from Myanmar, Thailand, Cambodia and central and southern Vietnam; 5 = colonies from Philippines (Luzon, Mindanao and Visayas); 6 = colonies from Philippines (Palawan), Malaysia and Indonesia.

In the text