Agricultural Science Digest

  • Chief EditorArvind kumar

  • Print ISSN 0253-150X

  • Online ISSN 0976-0547

  • NAAS Rating 5.52

  • SJR 0.156

Frequency :
Bi-monthly (February, April, June, August, October and December)
Indexing Services :
BIOSIS Preview, Biological Abstracts, Elsevier (Scopus and Embase), AGRICOLA, Google Scholar, CrossRef, CAB Abstracting Journals, Chemical Abstracts, Indian Science Abstracts, EBSCO Indexing Services, Index Copernicus
Agricultural Science Digest, volume 41 issue 3 (september 2021) : 397-404

Detection of Heterozygosity and Allelic Frequency of Indian and Southeast Asian Types of Oecophylla smaragdina (Fabricius) (Hymenoptera, Formicidae) in Bangladesh using Microsatellite Markers

Md Mamunur Rahman1,2,*, Shingo Hosoishi1, Kazuo Ogata1
1Institute of Tropical Agriculture, Kyushu University, Fukuoka 819-0395, Japan.
2Department of Entomology, Faculty of Agriculture, Bangabandhu Sheikh Mujibur Rahman Agricultural University, Gazipur 1706, Bangladesh.
Cite article:- Rahman Mamunur Md, Hosoishi Shingo, Ogata Kazuo (2021). Detection of Heterozygosity and Allelic Frequency of Indian and Southeast Asian Types of Oecophylla smaragdina (Fabricius) (Hymenoptera, Formicidae) in Bangladesh using Microsatellite Markers . Agricultural Science Digest. 41(3): 397-404. doi: 10.18805/ag.D-281.
Background: The two types of Asian weaver ant species, Oecophylla smaragdina, Indian types and Southeast (SE) Asian types can co-occur in the same habitat. Previous study on mtDNA based phylogenesis of O. smaragdina in Bangladesh has revealed the presence of the mixture of both Indian and SE Asian type of O. smaragdina. However, the nuclear DNA based study showed some inconsistency leading the chance of hybridization as colonies encountered from two localities contained O. smaragdina workers from more than one type as identified via longwave length rhodopsin (LWRh) and mitochondrial COI and Cytb gene analysis. The purpose of present study was to detect hybridization status and to find out the identical allele for Indian and SE Asian type of O. smaragdina from different colonies.  

Method: A total of seven O. smaragdina colonies had been collected from Bangladesh during 2013 to December 2016. Microsatellite analysis were performed for Gene mapping and heterozygosity study with 11 microsatellite loci through electropherogrum study generated by peak scanner software.

Result: Among 11 microsatellite loci, two loci, MS 8.24 and MS 6.45 provided the identical allele position for Indian and SE Asian types while no allelic position were identified with nine microsatellite loci. The results inferred that weaver ant of both Indian and SE Asian type engage in bidirectional interspecific gene flow. The strong possibility of hybridization between these two types is inferred, which indicates that barriers to hybridization do not exist or can be completely overcome. This is unexpected, given that they are non-sister species and broadly sympatric in nature. These findings also indicate the possible occurrence of hybrid-hybrid crosses in some localities of weaver ant, a phenomenon rarely observed in ants. 
The weaver ant species, Oecophylla smaragdina, is widely distributed from India through Southeast Asia to northern Australia including many tropical western Pacific islands (Weterrer, 2017). The ant is arboreal, observed in natural forests, fruit orchards and homestead woodlands and is one of the important biological control agents (Dlusky, 2008). The phylogenetic study of O. smaragdina based on mitochondrial Cytb and COI genes identified two major types where Indian types occurred mainly in India and Sri Lanka while the Southeast Asian (SE Asian) types have been observed in most of the SE Asian countries including Bangladesh (Azuma et al., 2006, 2002). However, recent phylogenetic study revealed the occurrence and dominancy of Indian type in the western part (Rahman et al., 2017a) and SE Asian types in the eastern part and overlapping population of O. smaragdina in the central part of Bangladesh (Rahman et al., 2017b) that can lead the chance of hybridization. Recently, for inferring the evidence of hybridization a comprehensive view of evolutionary history by analyzing nuclear and mitochondrial DNA between species level phylogeny in Hymenoptera was found effective and has been used extensively (Roos et al., 2011; (Ascher et al., 2001; Cameron and Williams, 2003). Azuma et al., (2006) through long wavelength rhodopsin (LWRh) analysis categorized the Indian and SE Asian type of O. smaragdina population as Smaragdina B and Smaragdina A and suggested that SE Asian are derived and monophyletic while   the Indian type is the ancestral. Discordant genetic relationships usually occurred from mitochondrial introgression or the incomplete lineage sorting (Eto et al., 2013) which was observed from the two types of weaver ants in Bangladesh (Rahman et al., 2020). The occured discordancy can partially reveal the hybrid phenomenon but cannot fully detect the hybridization status of the colony by revealing the identical alleles.

Microsatellite markers, also known as simple sequence repeats (SSRs) is very useful for revealing the heterozygosity and finding the identical alleles with combination of mitochondrial and nuclear DNA results. It’s a nuclear DNA marker and combines all the alleles from each parental species for F1 detection and can be used to identify the polymorphism by detecting heterozygosity and with high rate of mutation (Loxdale and Lushai, 1998; Hancock, 1999).  The microsatellite markers consist of short but repeated units of around two to six base pairs in length and include both coding and non-coding regions of up to 200 bp (Beukeboom and Zwaan, 2005). By combining the alleles from each parental species in all loci, microsatellites can indicate F1 hybrids as microsatellite is considered as nuclear DNA marker, although the event of backcrossing has the significant effect on detection of hybridization (Goodman et al., 1999).  Through mtDNA, only the maternal identity can be traced without the detection of hybrid. Therefore, Hybridization can be revealed through the inconsistency result of mt DNA and nuclear DNA analysis along with backcrossing with the parental specie (Mardulyn and Cameron, 1999). However, by only microsatellite study, hybridization can be undetected if backcrossing occurred repeatedly because back crossing sometimes causes the complete elimination of the nuclear genome of the maternal species from hybrid (Schilck-Steiner et al., 2011).
 
Mallet (2005) reported that hybridization often takes place among the closely related species because of insufficient evolvement of reproductive barriers. This may lead to inconsistency between nDNA and mtDNA types within the cell of a hybrid and thus, the occurrence of both Indian and SE Asia types, created the possibility of a hybrid colony of O. smaragdina, especially in the overlapping zone in Bangladesh. In the present context, four possible combinations can happened within the cell if the crossing between the male and female Indian and SE Asian types occurred. The occurrence of either Indian or SE Asian in both mtDNA and nDNA is consistent, i.e., both the mtDNA and nDNA are Indian type or both are the SE Asian type. The inconsistency might be found out if the nDNA are Indian type and mtDNA are SE Asian type or the vice-versa (Fig 1a). In the case, where mtDNA are SE Asian type, the crossing might to be expected between the male of Indian type and female of SE Asian type and mtDNA are maternally inherited. Through, comparing the results of the nucleotide sequences obtained from nDNA analysis by LW Rh and mtDNA and nDNA inconsistency was observed (Rahman et al., 2020) and it partially reveals the hybrid phenomenon but cannot fully detect the hybridization status of the colony. In addition, in the diploid cell, the nucleus can contain both the Indian and SE Asian types. However,  the crossing pattern of occurring either the male or female from Indian or SE Asian type can be detected through mtDNA analysis (Fig 1b). Microsatellite markers could identify the heterozygosity in the above-mentioned condition. The detailed explanation about the justification of choosing nuclear DNA and microsatellite analysis for revealing the hybridization status were described in Fig 1 Therefore, the purpose of this study was to facilitate further confirmation of heterozygosity along with the actual hybrid condition within the diploid nucleus by detecting identical allele between Indian and SE Asian types within the colony.
 

Fig 1: mtDNA and nDNA inconsistency pattern and the heterozygous status within the cell due to crossing of male and female of different Indian and SE Asian types.

Localities of O. smaragdina sampling in Bangladesh
 
A total of seven O. smaragdina colonies had been collected from seven localities belonging to eight divisions of Bangladesh during 2013 to December 2016. The detailed locality information is presented in Fig 2 and Table 1.
 

Fig 2: The sampling sites of Oecophylla smaragdina used in microsatellite study. The locality information is corresponding to that of Table 1.


 

Table 1: Detailed locality information of weaver ant sampling in Bangladesh for microsatellite study.


 
Sample preparation
 
The samples of O. smaragdina collected from Bangladesh were preserved in 99% ethanol prior to DNA extraction. All the molecular analysis and sequencing were done in the laboratory Institute of Tropical Agriculture of Kyushu University, Japan.
 
DNA extraction
 
For molecular analyses, specimens were randomly chosen from seven colonies of seven different localities. Ten individuals from each colony were selected randomly for extraction of DNA. For the convenient discussion, we grouped sample no. 133 to 142 from L3, 143 to 152 from L4, 153 to 160 from L2, 163 to 170 from L7, 173 to 180 from L1, 181 to 188 from L5 and 189 to 195 from L6. Genomic DNA was extracted from the fore, middle and hind legs of specimens that were preserved in alcohol by using QIAGEN DNeasy Blood and Tissue kit (Qiagen, Maryland, USA) following manufacturer’s instruction. Samples were vortexed by adding 180 µl tissue lysis buffer (Buffer-ATL) and 20 µl proteinase K. Buffer-AL (lysis buffer) of 200 µl were added to the samples and were incubated at 55°C for 48 hours. The contaminant of DNA was washed out using two wash buffers (Buffer- AW1 and Buffer- AW2). The unbound proteins were removed by the elution buffer (Buffer- AE). Extracted DNA were purified using DNA purification buffer (Buffer- AL). All centrifugation steps were completed at room temperature as mentioned by Crozier et al., (1993) and Crozier et al., (1994). The colony mates of the specimens used for DNA analysis were preserved in the laboratory of Institute of Tropical Agriculture, Kyushu University after DNA extraction.
 
Primers and PCR status followed for microsatellite study
 
The primers used in PCR and with its corresponding loci information are given in Table 2. PCR was performed in a 15 μl reaction volume containing 1 μl × reaction buffer (10XEx7Eq buffer), 1 μl × dNTP with 1.8 mM MgCl2, 1 μl of each primer, 0.1 μl of Ex Taq DNA polymerase (TAKARA, USA) and about 1 μl template DNA. Forward primers were 5'-end labeled with 6-FAM fluorescent dye. Cycling parameters were followed as mentioned by Azuma et al., (2004). Fragments were run on an ABI PRISM 310 genetic analyzer (Applied Biosystems). For gene mapping by fragment analysis and data processing, Peak Scanner software vr 1.0 (Applied Biosystems) was used. For identifying the heterozygosity by the observation of number of peaks by sharing the locus number were determined by electropherogrum study given in Fig 3. This figure is the representative method of the total study.
 

Table 2: List of loci used in microsatellite analysis of O. smaragdina in Bangladesh modified from Azuma et al. (2004).


 

Fig 3: Electropherogrum study to find out the sharing allele position of each individuals per colony.

In this study, 11 microsatellite primers (MS 2.2.2, MS 2.3, MS 5.10, MS 2.14, MS 3.2, MS 6.45, MS 6.29, MS 5.2, MS 7.1, MS 8.24 and MS 6.47) (Azuma et al., 2004) were used. However, we failed to analyze the samples using MS 5.2 microsatellite primer. The results of heterozygosity study using each microsatellite primer were presented in Table 3.
 

Table 3: Allele record for 10 microsatellite loci in seven localities.


 
Use of different microsatellite primers for detecting heterozygosity and allelic frequency
 
Microsatellite primer MS 5.10
 
With this primer, we analyzed a total of 68 samples from seven localities. Among them, all the samples from the locality L4, L46, L27, L67 L50 and L66 were sharing both the loci 98 and 99 while the samples from locality 13 were sharing both the loci 96 and 97 (Table 3). L13 was identified as Indian type of Oecophylla smaragdina by both mitochondrial and nuclear DNA analysis. However, L34, L46, L67, L50 was determined as SE Asian type and L27 and L66 was characterized as Indian type. So, in this case for the primer MS 5.10, we failed to distinguish identical alleles for Indian and SE Asian types although we found heterozygosity as it showed the two peaks in almost all the samples.
 
Microsatellite primer MS 2.14
 
This primer  detected several allele positions for Indian and SE Asian types. The samples from L46, L27 shared the loci positions 231 and 234 whereas sample from L34 shared the locus 221, in addition, 1 sample from L50 shared the locus 262 with loci 231 and 234.  The samples from locality L13 and L66 showed a quite different pattern of sharing allele frequencies. The samples from L13, shared the loci 245, 247, 262, 270 and 273. The samples from L66 shared the loci 238, 252 and 270 (Table 3). So, no distinguished locus position sharing by Indian and SE Asian types of samples were observed.
 
Microsatellite primer MS 3.2
 
For identifying the identical allele, I tested 32 samples from 4 localities. The samples from L27 and L67 shared the locus 81 while the samples from L13 shared the locus 86. The samples from locality 50 shared both locus 81 and locus 86 (Table 3).
 
Microsatellite primer MS 7.1, MS 6.29 and MS 6.47
 
Eight samples from each of two localities (L67 and L13) were tested by microsatellite primer MS 7.1, MS 6.29 and MS 6.47. All 08 samples shared the locus 100, locus 128 and locus 172 for MS 7.1, MS 6.47 and MS 6.29, respectively (Table 3). So, no locus was identified that can separate Indian and SE Asian type of Oecophylla in Bangladesh.
 
Microsatellite primer MS 2.3 and MS 2.2.2
 
Sixty samples from seven locaities were tested each for MS 2.3 and MS 2.2.2. All the samples from MS 2.3 shared the locus 144, while two samples from L34 and L46 shared both the loci 100 and 144. In the case of MS 2.2.2, all the samples sharing the locus 103 except the samples from L13 and L66 that were sharing the locus 97. The samples from other locality that were identified as Indian type also shared the locus position 97 along with some samples from L34 which was identified as SE Asian type. So, with these two primers the identical allele position was not possible to distinguish between Indian and SE Asian types (Table 3).
 
Successful detection of identical allele
 
The identical allele between Indian and SE Asian types of O. smaragdina population in Bangladesh were possibly detected with the microsatellite primer MS 6.45 and MS 8.24. For MS 6.45, four samples from the locality L67 were analyzed and all 04 samples shared the locus 238. However, the samples from locality L13 shared the locus 244 and locus 250. As the colony of L13 locality was identified as Indian type and L67 was SE Asian type so, these locus positions was identical. Similarly, the samples with 8.24, among 08 samples, 4 samples from L67 shared the locus 277 while the rest 4 samples from locality L13 shared the locus 262 (Table 3).
 
There was not too many evidence of such heterozygous condition within the colony of O. smragdina in India or any other SE Asian country and this is the first report of such mixed colony in Bangladesh as well. Similarly,  Roos et al., (2011) studied tracing the evolution and hybridization of colobine monkey in the Asian continent, found several hybridization patterns by testing the mitochondrial and nuclear DNA which was found to be supported by Blaimer (2012). The observed phylogenetic incongruences  resulted by hybridization among ancestral lineages most likely occured, due to the presence of potential contact zones like today’s Bangladesh, Myanmar and the northeast of India, which is suggested as hybridization area (Karanth et al., 2008). However, several big mountains and big rivers in the border region of today’s Myanmar, India and China might have been a possible diversification hotspot (Chakraborty et al., 2007) which leads to develop such hybridization pattern.
 
Divergence time is also an important factor for many ant genera. Oecophylla is thought to be a significant factor in such distribution (Bolton, 1995). Diversification within groups in this continent was recorded from the middle Pliocene to early Pleistocene (Azuma et al., 2002). After this period, the world has encountered a significant climatic change. It might also affect the distribution of Oecophylla in different parts of the world. During Last Glacial Maximum (LGM), the tropic region shifted southward and it retained northward after glaciation. This study suggested that the Indian and Southeast Asian clades of O. smaragdina expanded their distribution northward along suitable regions with high temperature and humidity and then the two clades supposedly encountered and overlapped in central Bangladesh. Similar trends were also observed in the case of study of the origin of Asian elephants (Vidya et al., 2009). They suggested a contraction-expansion scenario during climatic oscillation leads to geographical overlaps of two mtDNA clades created the allopatric population of Asian elephants in India, Sri Lanka and Myanmar. In the case of weaver ant in Bangladesh, as both the Indian and SE Asian types were dominated  in Western and Eastern part of the country, respectively and there was not such a big border of separating those two populations (Rahman et al., 2017a). The probability of the contact of both the types was considered very usual. Pusch et al., (2006) reported that ant colony with heterozygous produce hybrid workers. This type of gene introgression may increase the genetic diversity of the hybrid relative to its parental species and can lead to hybrid vigor (Helms and Keller, 2003). In this case, the O. smaragdina colony from the locality 11 and 15 in Bangladesh might have the influence of the arising of new evolutionary lineages. This hybridization therefore, can be considered as a process of the evolutionary significance of O. smaragdina in Bangladesh.
We are thankful to the authority of Bangabandhu Sheikh Mujibur Rahman Agricultural University (BSMRAU), Authority of Bangladesh Tea Research Institute, Sreemangal, Sylhet for helping us to collect samples. We express our sincere gratitude to Dr. Masaru Matsumoto, Institute of Tropical Agriculture of Kyushu University, for his technical support. We are grateful to Dr. Noriko Azuma and Ms. Yukiko Asaka, Sapporo, Japan for their valuable advice and suggestions. We are also thankful to Dr. Akinori Ozaki, Institute of Tropical Agriculture for providing facilities during sampling in Bangladesh. This work was supported in part by JSPS KAKENHI [Grant-in-Aid for Scientific Research (B)] Grant Number 26304014, MEXT, Japan.

  1. Ascher, J.S., Danforth, B.N. and Ji, S. (2001). Phylogenetic utility of the major opsin in bees (Hymenoptera: Apoidea): A reassessment. Molecular Phylogenetic Evolution. 19: 76-93. doi:10.1006/mpev.2001.0911.

  2. Azuma, N., Kikuchi, T., Ogata, K. and Higashi, S. (2002). Molecular phylogeny among local populations of weaver ant Oecophylla smaragdina. Zoological Science. 19: 1321-1328. doi:10.2108/zsj.19.1321.

  3. Azuma, N., Takahashi, J., Higashi, S., Sasaki, M. (2004). Microsatellite loci for the weaver ant Oecophylla smaragdina. Mol. Ecol. Notes. 4. 608-610. doi:10.1111/j.1471-8286.2004.00752.x.

  4. Azuma, N., Ogata, K., Kikuchi, T. and Higashi, S. (2006). Phylogeography of Asian weaver ants, Oecophylla smaragdina. Ecological Research. 21: 126-136. doi:10.1007/s11284-005-0101-6.

  5. Beukeboom, L.W., Zwaan, B.J. (2005). Genetics. In: Insects as Natural Enemies: a Practical Perspective. [Jervis M.A. (eds)]. Springer. pp. 167-218.

  6. Blaimer, B.B. (2012). Taxonomy and species-groups of the subgenus Crematogaster (Orthocrema) in the Malagasy region (Hymenoptera, Formicidae). Zookeys. 199: 23-70. doi:10.3897/zookeys.199.2631.

  7. Bolton, B. (1995). A new general catalogue of the ants of the world. Harvard University Press, London.

  8. Cameron, S.A. and Williams, P.H. (2003). Phylogeny of bumble bees in the New World subgenus Fervidobombus (Hymenoptera: Apidae): Congruence of molecular and morphological data. Molecular Phylogenetic Evolution. 28: 552-563. doi:10.1016/S1055-7903(03)00056-3.

  9. Chakraborty, D., Ramakrishnan, U., Panor, J., Mishra, C. and Sinha, A. (2007). Phylogenetic relationships and morphometric affinities of the Arunachal macaque Macaca munzala, a newly described primate from Arunachal Pradesh, northeastern India. Molecular Phylogenetic Evolution. 44: 838-849. doi:10.1016/j.ympev.2007.04.007.

  10. Crozier, R.H. and Crozier, Y.C. (1993). The mitochondrial genome of the honeybee Apis mellifera: Complete sequence and genome organization. Genetics. 133: 97-117. doi:10.1111/j.1365-2583.1993.tb00131.x.

  11. Crozier, R.H., Dobric, H.T., Graur, D., Corneut, J.M. and Taylor, R.W.(1994). Mitochondrial DNA sequence evidence on the phylogeny of Australian Jack-jumper ants of Myrmecia pilosula complex. Molecular Phylogenetic Evolution. 4(1): 20-30.

  12. Dlussky, G.M., Wappler, T., Wedmann, S. (2008). New middle Eocene formicid species from Germany and the evolution of weaver ants. Acta Palaeontol. Pol. 53: 615-626. doi:10.4202/app.2008.0406.

  13. Eto, K., Matsui, M. and Sugahar, T. (2013). Discordance between mitochondrial DNA genealogy and nuclear DNA genetic structure in the two morphotypes of Rana tagoi tagoi (Amphibia: Anura: Ranidae) in the Kinki Region, Japan. Zoological Science. 30(7): 553-558. doi:10.2108/    zsj.30.553.

  14. Goodman, S.J., Barton, N.H., Swanson, G., Abernethy, K., and Pemberton, J.M. (1999). Introgression through rare hybridization: a genetic study of a hybrid zone between red and sika deer (genus Cervus) in Argyll, Scotland. Genetics. 152: 355-371.

  15. Hancock, J.M., Simon, M. (2005). Simple sequence repeats in proteins and their significance for network evolution. Genetics. 345: 113-118.

  16. Helms Cahan, S. and Keller, L. (2003). Complex hybrid origin of genetic caste determination in harvester ants. Nature. 424: 306-309.

  17. Karanth, K.P., Singh, L., Collura, R.V. and Stewart, C.B. (2008). Molecular phylogeny and biogeography of langurs and leaf monkeys of South Asia (Primates: Colobinae). Molecular Phylogenetic Evolution. 46: 683-694. doi: 10.1016/j.ympev.2007.11.026.

  18. Loxdale, H.D. and Lushai, G. (1998). Molecular markers in entomology. Bulletin of Entomological Research. 88: 577-600.

  19. Mallet, J. (2005). Hybridization as an invasion of the genome. Trends in Ecology and Evolution. 20: 229-237.

  20. Mardulyn, P. and Cameron, S.A. (1999). The Major Opsin in Bees (Insecta: Hymenoptera): A promising nuclear gene for higher level phylogenetics. Molecular Phylogenetic Evolution. 12: 168-176. doi:10.1006/mpev.1998.0606.

  21. Pusch, K., Heinze, J., Foitzik, S. (2006). The influence of hybridization on colony structure in the ant species Temnothorax nylanderi and T. crassispinus. Insectes Soc. 53: 439-445. doi:10.1007/s00040-005-0891-8.

  22. Rahman, M.M., Hosoishi, S. and Ogata, K. (2017b). Phylogenetic analysis reveals the overlapping distribution of the Indian and southeast Asian clades of Oecophylla smaragdina (Fabricius). Journal of the Faculty of Agriculture, Kyushu University, 62: 429-434.

  23. Rahman, M.M., Hosoishi, S. and Ogata, K., (2017a). Phylogenetic position of the western Bangladesh populations of weaver ant, Oecophylla smaragdina (Fabricius) (Hymenoptera, Formicidae). Sociobiology. 64: 437. doi:10.13102/sociobiology.v64i4.1153.

  24. Rahman, M.M., Hosoishi, S. and Ogata, K. (2020). Discordance between mitochondrial and nuclear DNA genes suggests the possibility of hybridization of Indian and southeast Asian types of Oecophylla smaragdina (Fabricius) (Hymenoptera, Formicidae) in Bangladesh. Indian Journal of Agricultural Research (In press).

  25. Roos, C., Zinner, D., Kubatko, L.S., Schwarz, C., Yang, M., Meyer, D., Nash, S.D., Xing, J., Batzer, M.A., Bremeier, M., Leendertz, F.H., Ziegler, T., Perwitasari-Farajallah, D., Nadler, T., Walter, L. and Oshtorholz, M. (2011). Nuclear versus mitochondrial DNA: evidence for hybridization in colobine monkeys. BMC Evolutionary Biology. 11: 77. doi:10.1186/1471-2148-11-77.

  26. Schlick-Steiner, B.C., Steiner, F.M., Seifert, B., Stauffer, C., Christian, E. and Crozier, R.H. (2010). Integrative taxonomy: a multi-source approach to exploring biodiversity. Annual Review of Entomology. 55: 421-438.

  27. Vidya, T.N., Sukumar, R., Melnick, D.J.(2009). Range-wide mtDNA phylogeography yields insights into the origins of Asian elephants. Proc Biol Sci. 276: 893-902. doi:10.1098/rspb.2008.1494.

  28. Wetterer, J.K. (2017). Geographic distribution of the weaver ant Oecophylla smaragdina. Asian Myrmecology. 9: 1-12. doi:10.20362/am.009004. 

Editorial Board

View all (0)