Next Article in Journal
Fossils Shed a New Light on the Diversity and Disparity of the Family Limoniidae (Diptera, Nematocera)
Next Article in Special Issue
Landscape and Local Drivers Affecting Flying Insects along Fennel Crops (Foeniculum vulgare, Apiaceae) and Implications for Its Yield
Previous Article in Journal
Chance or Necessity—The Fungi Co−Occurring with Formica polyctena Ants
Previous Article in Special Issue
Effects of Urbanization on Plant–Pollinator Interactions in the Tropics: An Experimental Approach Using Exotic Plants
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metabarcoding in Diet Assessment of Heterotrigona itama Based on trnL Marker towards Domestication Program

1
Agrobiodiversity and Environmental Research Centre, MARDI Headquarter, 43400 Serdang, Selangor, Malaysia
2
Centre for Insect Systematics, Department of Biological Sciences and Biotechnology, Faculty of Science and Technology, Universiti Kebangsaan Malaysia, 43600 Bangi, Selangor, Malaysia
3
Department of Biological Sciences and Biotechnology, Faculty of Science and Technology, Universiti Kebangsaan Malaysia, 43600 Bangi, Selangor, Malaysia
*
Author to whom correspondence should be addressed.
Submission received: 30 January 2021 / Revised: 15 February 2021 / Accepted: 17 February 2021 / Published: 28 February 2021
(This article belongs to the Collection Insects on Pollination Ecology)

Abstract

:

Simple Summary

Heterotrigona itama is the most abundant stingless bee species commonly found in many ecosystems and is actively domesticated in the meliponiculture industry. In this research, the honey quality is the main criterion to be evaluated and is correlated with the plant species that the bees consume as their main diet. The objective of this study was to investigate the diet of H. itama derived from 12 populations in Malaysia. As a result, five plant phyla and 70 plant families with 262 species were obtained, of which four species were the abundant species consumed by H. itama. These findings are very valuable for strategizing the management of the H. itama domestication process specifically in a mono-cropping system and in a netted structure.

Abstract

Honey quality is the main criterion used for evaluating honey production in the stingless bee Heterotrigona itama, and it is correlated with the plant species consumed as its main diet. The objective of this study was to obtain the metabarcode data from 12 populations of H. itama species throughout Malaysia (Borneo and Peninsular Malaysia) using the trnL marker. A total of 262 species under 70 families and five phyla of plants were foraged by H. itama in the studied populations. Spermatophyta and Magnoliophyta were recorded as the two most abundant phyla foraged, at 55.95% and 32.39%, respectively. Four species, Garcinia oblongifolia, Muntingia calabura, Mallotus pellatus, and Pinus squamata, occurred abundantly and were consumed by H. itama in all the populations. These data are considered as a fundamental finding that is specific to the diet of H. itama for strategizing the management of the domestication process specifically in a mono-cropping system and in a netted structure. Thus, based on these findings, we recommend Momordica charantia, Melastoma sp., and Cucumis sativa as the best choices of food plant species to be planted and utilized by H. itama in meliponiculture.

1. Introduction

The stingless bee is known as a pollinator of many plant species in forests and agricultural ecosystems. For example, stingless bee species have become a pollinating agent of the mangosteen (Garcinia mangostana), liposu (Bacaurea lanceolate), and star fruit (Averrhoa bilimbi) [1]. From over 500 species of stingless bees recorded worldwide [2], about 45 species are found in Malaysia [3], with Heterotrigona itama Cockerell, 1918 being the most common and abundant species in most ecosystems. This is most likely due to the availability of its food sources [4]. The ubiquitous H. itama has become the domesticated species of choice in our local meliponiculture industry [5], as well as in neighboring countries [6]. To date, 50,000 colonies of H. itama have been reared for honey production in Malaysia and this trend is increasing yearly.
According to Vit et al. (2013) [7], H. itama is able to continue its honey production due to the exclusive interactions and coevolution between the stingless bees and specific angiosperm species of flowering trees [8]. In the meliponiculture industry, the breeders and the market players are keen to evaluate the honey quality, which is determined by the plant nectar and/or pollens that are consumed by the bees. In this regard, the best flowering plant species in relation to H. itama should be planted as the main food source, for in situ reproduction, as well as for ex situ conservation [9].
According to Afik et al. (2014) [10], Absy et al. (2018) [11], and Jaapar et al. (2019) [12], H. itama normally collects pollen, nectar, and resin when it forages away from the nest. However, in some cases, the bees also collect and obtain nutrients or salts from animal carcasses [13]. Therefore, in order to utilize this species within the setting of a netting structure, it is crucial that the food plant species is prioritized. This is because the plants depend solely on the stingless bees to act as their pollinating agent to perform pollination and, simultaneously, for the operators to increase the honey yield [14]. Furthermore, the stingless bee rearing program within the orchard, and in the natural habitat, was determined as the best practice to support the ecosystem services and to protect the bee species from population decline and extinction [15].
The application of DNA barcoding has become the latest approach or tool in biotechnology, showing benefits for human health and marine sustainability [16]. The application of DNA molecules has also solved many problems and issues related to the meliponiculture industry worldwide, such as accessing the honey compounds in several honey varieties [17], stingless bee identification through DNA barcoding [18], and plant species identification as food sources for the stingless bees [19].
Information on the plant species foraged or used as the diet for H. itama has never been widely investigated despite the great success of H. itama in the meliponiculture industry, specifically in honey production. However, a lot of information is available on the pollen collection of H. itama, although the question as to whether the bees consume pollen or only store it in pollen pots has yet to be clarified. In this regard, metabarcoding analysis by using a trnL marker on the H. itama diet is highly crucial and is the main objective for this study. It is important that the flowering plants consumed in the natural or wild habitat by H. itama are nurtured on the farms or orchards located together with H. itama colonies for sustainable reproduction; this will also ensure the best quality of honey production.

2. Materials and Methods

2.1. Sampling of the Stingless Bees

A total of 36 individuals of Heterotrigona itama were sampled and collected from 12 meliponiculture sites throughout Malaysia (Table 1 and Figure 1), with three individuals collected from three colonies from each locality.

2.2. DNA Extraction

Three individuals from different colonies in each locality were dissected and the DNA was extracted from the whole body using NucleoSpin® DNA Insect protocols (Macherey-Nagel, Germany). DNA samples from three individuals (25 µL × 3) were pooled from each locality prior to amplification and totaling up to only 12 DNA samples for the amplification process.

2.3. DNA Sample QC

The quality of the pooled purified DNAs samples was monitored on 1% TAE agarose gel. The DNA concentration was measured using a spectrophotometer (Implen NanoPhotometer® N60/N50) and fluorometric quantification using an iQuant™ Broad Range dsDNA Quantification kit. Preparation of the library was done through two rounds of polymerase chain reaction (PCR). The first PCR was carried out using trnL primers, PCR forward primer, g-A49425 (5′-GGGCAATCCTGAGCCAA-3′) and PCR reverse primers, H-B49466 (5′- CCATTGAGTCTCTGCACCTATC-3′) [20]. All the PCR reactions were carried out with Q5® High-Fidelity DNA Polymerase (NEB). A total of 25 µL first PCR MasterMix consisted of 12.5 µL of 2× KAPA HiFi HotStart ReadyMix, 5 µL of forward and reverse primers, and 2.5 µL DNA. The amplification was performed under the following protocol: initial denaturation of 95 °C for 3 min, followed by 25 cycles of denaturation at 95 °C for 30 s, primer annealing at 55 °C for 30 s and extension at 72 °C for 30 s, with final elongation at 72 °C for 5 min, and stored at 4 °C until PCR clean-up.
The DNA was purified using AMPure XP beads to purify the trnL amplicon away from free primers and primer dimers in the PCR clean-up procedure. The AMPure Xp beads were taken out and left at room temperature. Amplicon PCR products were centrifuged at 1000× g at 20 °C for 1 min to collect condensation. AMPure XP beads were vortexed for 30 s to make sure that the beads were evenly dispersed. A total of 20 µL of AMPure XP beads were added to each well of the amplicon PCR plate and entire volumes were pipetted up and down 10 times gently. Then the sample was incubated for 5 min at room temperature. The plate was placed on a magnetic stand for 2 min and then the supernatant was removed. Then, the beads were washed with 80% alcohol. Each sample was added to 200 µL 80% ethanol and was incubated for 30 s and the supernatant was removed afterwards. This procedure was repeated twice and ethanol was removed. The amplicon PCR plate was placed on the magnetic stand for 10 min to allow the beads to air-dry. Then, the amplicon PCR plate was removed from the magnetic stand and to each well was added 52.5 µL of 10 mM Tris pH 8.5 followed by mixing to ensure that the beads were fully resuspended. Amplicon PCR plates were incubated at room temperature for 2 min, then placed on the magnetic stand for 2 min until the supernatant was clear. Only 50 µL of supernatant from the amplicon PCR plate were transferred to a new 96-well PCR plate.

2.4. Second PCR (Index PCR)

This involved a dual indices Illumina sequencing adapter using a Nextera XT index kit. The PCR mixture consisted of 25 µL of 2× KAPA HiFi HotStart ReadyMix, 5 µL of Index 1 primers (N7XX) and Index 2 primers (S5XX), and 10 µL PCR grade water, then the plate was centrifuged at 1000× g at 20 °C for 1 min. PCR amplification for DNA templates with indexes was performed using the following profiles: polymerase activation at 95 °C for 3 min, followed by 8 cycles of denaturation at 95 °C for 30 s, annealing at 55 °C for 30 s, and extension at 72 °C for 30 s, and final extension at 72 °C for 5 min. The quality of the libraries was measured using Agilent Bioanalyzer 2100 System by Agilent DNA 1000 Kit and fluorometric quantification by Helixyte GreenTM Quantifying Reagent.

2.5. Next-Generation Sequencing

The libraries were normalized and pooled according to the protocol recommended by Illumina and proceed to sequencing using MiSeq platform using 250 PE.

2.6. Data Analysis

Paired-end reads were first had sequence adaptors and low-quality reads removed using BBDuk of the BBTools package. After this, the forward and reverse reads were merged using USEARCH v11.0.667 [21]. All sequences that were shorter than 150 bp or longer than 600 bp were removed from the downstream processing. Reads were then aligned with trnL sequences based on NCBI database. After these quality assessment steps, reads were clustered de novo into operational taxonomic units (OTUs) at 97% similarity using UPARSE v11.0.667 [22]; rare OTUs with fewer than 2 reads (doubleton), which were often spurious, were deleted from downstream processing. A single representative sequence from each OTU was randomly chosen, and Pynast was used to align and construct a phylogenetic tree to compare with NCBI database. Then, taxonomic assignment of OTUs was achieved using QIIME V1.9.1 [23] A phylogenetic dendrogram was constructed at the species level using Bray–Curtis distance with 1000 bootstraps to define the relationship between localities of H. itama using Paleontological Statistics Software Package for Education (PAST 3) software. Other statistical analysis was done in R V3.6.1 [24].

3. Results

A total of 527,398 plant trnL sequences were generated from 12 samples of H. itama and ranged from 7448 to 56,612. Sabah TE showed the most reads (70,068), followed by Perlis (56,612), Sarawak Beladin (55,714), Kedah (54,937), and Perak (52,216), while the fewest reads were from Pool Sarawak-Bintulu (7448). The Shannon–Wiener index (H’) indicated that Sabah TE had the highest diversity with H’ of 2.73 and 63 OTU, and the lowest was recorded for Selangor at 0.644 and 19 OTU. The same was found for Simpson 1-D and Chao1, in which Sabah TE showed the most OTUs (1 − D = 0.899, Chao1 = 83), while the fewest was for H. itama from Selangor (1 − D = 0.243, Chao1 = 22.33) (Table 2). Furthermore, the OTU number from Borneo (Sabah TE, Sabah Telipok, Sarawak Kuching, Pool Sarawak Bintulu, Sarawak Beladin, Sabah Tuaran, and Sarawak Miri) was 328 and higher than that of Peninsular Malaysia (159). The rarefaction curve was not asymptote since the methodology yielded up to 7500 reads only. H. itama from Sabah TE and Sabah Telipok showed the highest sequencing depth compared to the other localities (Figure 2).
A total of 262 species from 70 families and five phyla were assigned at 97% similarity (Figure 3). The most common family recorded was Euphorbiaceae (27), followed by Fabaceae (20), Cucurbitaceae (16) and Fagaceae (13). Spermatophyta was the most abundant phylum (55.95%) recorded in samples from all the study localities, followed by Magnoliophyta (32.39%), Tracheophyta (4.72%), and Anthophyta (1.97%), while the least was Embryophyta (0.13%). There were also 4.85% unknown phyla due to the blast not hitting any plant DNA (Table 3). The Sabah Telipok data showed the highest abundance of Spermatophyta (18.24%), followed by Perak (17.07%), Selangor (15.20%), and Kedah (11.76%), while Pool Sarawak Bintulu showed the lowest (0.34%). Furthermore, Embryophyta was only recorded from Sarawak Beladin, indicating that this locality represented all the bee food plant phyla in this study (Table 4).
Based on the analyzed data, we found that samples from Sarawak Beladin presented all the six phyla that were collected by H. itama, while Magnoliophyta represented the most common phylum collected, at 78.7% (Table 3). Most of the Tracheophyta were found in Sabah Tuaran (18%), Anthophyta in Sabah Telipok (10.56%), Magnoliophyta in Sarawak Beladin (87.7%) Spermatophyta in Perak (95.8%), Embryophyta in Sarawak Beladin (1.16%), and the unknown phyla mostly found in Pool Sarawak Bintulu (80.8%) (Table 4). Most of the species in Tracheophyta were found in Sabah Tuaran (18%), Anthophyta in Sabah Telipok (10.56%), Magnoliophyta in Sarawak Beladin (87.7%) Spermatophyta in Perak (95.8%), and Embryophyta in Sarawak Beladin (1.16%), while an unknown phylum was mainly found in Pool Sarawak Bintulu (80.8%) (Table 4).
The heat map constructed showed the 30 most abundant plant species for 12 localities foraged by H. itama (Figure 4). The most dominant species found in all locations were Garcinia oblongifolia, Muntingia calabura, Mallotus pellatus, Pinus squamata, and Phoebe puwenensis. Meanwhile, Umtiza listeriana, Spathodea campanulate, Ischaemum aristatum, and Dillenia indica could only be found in several localities. Samples from Perak (P) and Perlis (Pr) showed that H. itama had consumed more from Garcinia oblongifolia (p = 56.71%, Pr = 69.63%), followed by Muntingia calabura (p = 32.67%, Pr = 21.34%).
Meanwhile, Sabah TE (STE), and Sabah Telipok (SBT) samples showed that the most abundant food plant utilized by H. itama was Phoebe puwenensis (STE = 25.82%, SBT = 31.56%). Kedah (K) and Selangor (S) showed the highest totals for Phoenix canariensis (K = 37.55%, S = 87.01%), while other samples showed that H. itama foraged mainly for Macarangga tanarius (Sarawak Beladin = 67.04%), Garcinia mangostana (Terengganu = 41.08%), Annona reticulata (Sabah Tuaran = 30.58%), Cucumis melo (Sarawak Miri = 49.93%), and Pinus squamata (Sarawak Bintulu = 4.89%).
The lowest totals of plant DNA found from all localities were: Typha angustifolia (Sabah Telipok = 0.167), Aclepias tuberosa (Perak = 0.16), Heliamphora pulchella (Perlis = 0.004%), Triptilion cordifolium (Sabah TE = 0.15%), Xanthocercis rabiensis (Sarawak Beladin = 0.076%), Chamaecrista polita (Kedah = 0.006%), Macarangga bancana (Sarawak Kuching = 0.01%), Iryanthera sagotiana (Terengganu = 0.058%), Mallotus mollisimus (Sabah Tuaran = 0.071%), Manihot pruinose (Sarawak Miri = 0.003%), Brassica oleraceae (Selangor = 1.680%), and Laguncularia racemose (Pool Sarawak Bintulu = 0.582%) (Figure 5). In addition, there were overlaps in plant species between the sampling sites (Figure 6). Samples from Sabah and Sarawak showed a larger overlap compared to between Sarawak and Peninsular Malaysia. There were close relationships of plant communities between results from Terengganu and Sarawak, Miri, and Sarawak, Kuching, and Sabah, Tuaran based on the phylogenetic dendrogram constructed based on Bray–Curtis distance (Figure 7).

4. Discussion

The information on the plant species that were foraged by the stingless bee, H. itama has been summarized in this study based on the data obtained from metabarcoding analysis. The metabarcode has been presented and proven to yield more thorough and rigorous findings compared to the melissopalynology technique [25]. Chloroplast and nuclear barcoding could be amplified from the pollen, and therefore, the metabarcode would give the best and most rapid means of identifying the plant species. Selvaraju et al. (2019) [26] noted that there were only 60 species under 34 families of trees revealed by the melissapalynology analysis from the west coast Malaysia. Likewise, Ghazi (2015) [27] also found only 59 plant species in an island ecosystem, and most of the species came from the underutilized fruits. However, more than 140 species of pollen from the melissapalynology analysis remained unidentified [1], which indicated that the method is still lacking compared to metabarcoding. A total of 262 plant species under 70 families was recorded in this study, therefore, we concluded that a broad diversity of plant species had been foraged by H. itama. Furthermore, in terms of the total plant species detected, our totals were distinctly higher than the melissapalynology method. Interestingly, an understory level plant family, i.e., Vitaceae, and a top canopy level family, i.e., Bombaceae, have been recorded to be consumed and foraged by H. itama using the metabarcode analysis. In addition, the species Pterisanthes stonei, considered as rare in the secondary forest, had also been recorded in the diet of H. itama, which is strongly supported by the results of the metabarcode in this study.
The samples of H. itama were collected from various localities throughout Malaysia (Peninsular Malaysia and Borneo) to represent a holistic view of H. itama’s diet, in consideration of its limited range of foraging behavior of up to 300 m distances only [28], and at a specific time and for a specific duration [29]. The data which were obtained from 12 localities at different settings of the meliponiculture sites presented a novel finding for understanding the diet of H. itama. Flower pollens are basically collected by H. itama, and then transferred to its hind legs before being stored in the pollen pots [30]. This behavior highlighted that H. itama either consumed the pollen or unintentionally attached the pollen onto its body. Therefore, the metagenomic analysis using trnL can explore the plant species really selected by H. itama as food for its energy sources [31]. The results obtained can be used as a management tool to domesticate H. itama species in the netting structure for pollination purposes in honey production. The data can be applied at meliponicultural sites to enhance the quality of the honey [8].
Generally, all H. itama samples collected from Borneo (Sabah and Sarawak) showed a large number of OTUs, ranging from 35 to 64, as compared to Peninsular Malaysia (19–42) (Table 2). This data supported the species richness of the food plants, which is much higher on Borneo compared to Peninsular Malaysia, as similarly found in the study of [32]. Our results also showed that in Borneo, the highest diversity index was from Sabah as compared to Sarawak. In observing the high richness and diversity of plants in Borneo, Saw and Chung (2015) [33] and Pereira et al. (2019) [34] reported that over 12,000 taxa of seed plants occur in Sabah and Sarawak, while only 8200 species occur in Peninsular Malaysia. According to Wan Mohd Jaafar et al. (2020) [35], a large coverage of forests would support a high diversity of plants species, as represented in Terengganu, Perak, and Kedah, with 654,625 ha, 1,019,052 ha, and 330,585 ha of forest area, respectively, which correspondingly serve as the main habitats for 80% of the stingless bee species in Malaysia [3].
Spermatophyta was the most abundant phylum found in this study, followed by Magnoliophyta. Spermatophyta is a plant group that produces seeds [36], while Magnoliophyta produces flowers. Both phyla are very active in producing a lot of pollen. As seed producers, members of Spermatophyta normally depend on pollinators, almost 80% of which are bee species [37,38]. Magnoliophyta, likewise, constantly produces flowers which attract the bees to collect pollen and nectar [39]. However, some of the Magnoliophyta species have a specifically pointed corolla structure to enable the bees to enter the flowers and collect the nectar. In such a situation, the bees will sometimes bore a hole in the corolla to steal the nectar without collecting the pollen from the flowers [40]. Many studies on the food source of stingless bees confirmed that Magnoliophyta and Spermatophyta were the two phyla most frequently visited by stingless bees [41].
Based on the analyzed data, we found that nectar from all the six phyla collected by H. itama were represented in the samples from Sarawak Beladin, in which Magnoliophyta was the most common at 78.7% (Table 4). However, the number of OTUs, which was only 39, was not as high as in the other localities (Table 2). Beladin is in the Betong division in Sarawak, where Meludam National Park is located. This national park has its own uniqueness due to some of the plant species recorded here that are poorly recognized [42], and most probably consist of endemic species. The national park still has 3.29% undisturbed area [43] and is identified as a refuge for Sarawak’s most endemic and rarest species [44]. Beladin, Sarawak also presented a very low percentage of Embryophyta (1.19%), but none of the species in this phylum could be found in other localities in Peninsular Malaysia and Borneo.
The family composition of our study results was dominated by Euphorbiaceae, Fabaceae, Fagaceae, and Cucurbitaceae (Figure 3). Normally, the stingless bees would forage, for their diet, for three types of honey sources, namely, nectar from flowers for nectar honey, nectar from insects for honeydew honey, and nectar from leaf shoots for extrafloral honey [45,46]. The Euphorbiaceae family produces unlimited extrafloral honey, and some of the species are associated with ants [47]. In Malaysia, the stingless bees harvest pollen from the flowers of Macarangga sp. and Hevea sp. to produce extrafloral honey. Studies by Bahri (2018) [48] highlighted that stingless bees were domesticated under rubber trees to increase yield and production. Although the rubber tree does not produce pollen constantly and only once a year normally, our data have revealed that Hevea sp. DNA was found inside the stingless bee bodies, indicating that besides the pollen, DNA could also be derived from other plant parts, as revealed by using metabarcoding.
From our study results, it can be concluded that in some situations, stingless bees will forage for nectar without collecting pollen. Referring to Jaapar et al. (2018) [49] and Asma et al. (2019) [50], stingless bees collect pollen in the morning from 9 a.m. to 11 a.m., and then again from 3 p.m. to 4 p.m. The worker bees also normally collect resin to produce extrafloral honey before returning to the hive [49,51]. Interestingly, members of Fabaceae, Fagaceae, and Cucurbitaceae normally produce yellowish-colored flowers that are preferably collected by the stingless bees [1], and Fabaceae has been identified as the most frequently visited family for collecting pollen by the stingless bees [52]. Thus, although Fagaceae had been determined as the dominant family in the mixed orchards or secondary forests wherein the meliponiculture activities were sited [53], H. itama, however, preferably selected plant species belonging to Fabaceae in their foraging behavior.
Based on the heat map (Figure 4), four species, namely, Garcinia oblongifolia, Mallotus pellatus, Pinus squamata, and Phoebe puwenensis, were found in all the localities to be foraged by H. itama [52,54,55]. According to Azuan et al. (2008) [56], Perak, Perlis, and Kedah showed high densities of Garcinia oblongifolia communities. Garcinia is a genus that produces fruits and resin for stingless bees, from which the nectar and resins are collected by the bees during the flowering season [57]. We estimated that the DNA of Garcinia oblongifolia and Pinus squamata was sourced from the plant resins, collected by the bees when foraging for structural materials for their hives [12]. Meanwhile, H. itama was attracted to both Mallatus pellatus and Phobe puwenensis due to the abundance of their nectar and pollen.
According to Nagamitsu et al. (1999) [4], several plants species from the ground to the canopy tops are foraged by stingless bees in their quest for nectar and pollen, which, consequently, would result in their pollination of the flowers. The stingless bee has also been identified as an important pollinator of dipterocarp species in the forests of Malaysia [58], as well as in agricultural areas [50]. However, a very good management strategy is required for the domestication process of stingless bees, especially in mono-cropping systems, as these insects play a major role as the pollinating agent for crops as well as for honey production.
Based on these findings, the data generated in this study can be utilized in many broad-scale applications. For example, DNA of the durian species Durio zibethinus has been detected in the gut of H. itama, even though it has not been recognized as a highly abundant species based on the metabarcode data. Therefore, there is a high potential of utilizing H. itama as an alternative pollinator to be domesticated in large acreage durian orchards in Malaysia and other tropical countries and is thus highly recommended. A study by Wayo et al. (2018) [59] has confirmed that stingless bee species are also pollinators of durian, besides the well-known bat species. The dual role of stingless bees in pollination and honey production can be fully exploited by farmers and local communities to generate more income. Thus, our study data and findings can be utilized in many areas of application to support our socio-economic activities in the near future. Likewise, in sericulture, the domestication process of stingless bees can be successfully developed under mono-cropping systems and in netted structures by planting all the possible plant species foraged by these bees to produce higher value and better quality honey as well as improved brood production as a promising way forward.

5. Conclusion

The taxonomic information (phylum, family, genus, species, OTU, diversity) on the plant species foraged by H. itama has been presented in this study through the metabarcoding analysis obtained from the DNA of adult H. itama. The fast-growing plant species under the phyla Spermatophyta and Magnoliophyta, such as Momordica charantia, Melastoma sp., and Cucumis sativa, which produce abundant flowers and nectar, are recommended to be planted in mono-cropping systems for an effective and successful domestication process within the netted structure. Therefore, the food plant species obtained from the metabarcoding analysis in this study are considered as key species in the domestication process of H. itama for high-quality honey production and for greater success in the meliponiculture industry.

Author Contributions

Methodology (sampling), M.S.Z., N.F.R., N.M.; software (analysis), J.F., A.S.B.; writing—original draft preparation, J.F.; writing—review and editing, S.Y.; J.F., A.S.B.; supervision, R.J., B.M.M.-Z., S.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Agriculture, Malaysia (KG0020910) and GP-2019-K013317.

Institutional Review Board Statement

Not applicable.

Acknowledgments

The authors would like to thank Maimon Abdullah for her kind editing and critical comments of the final draft copy of this publication.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ghazi, R.; Zulqurnain, N.S.; Azmi, W.A. Melittopalynological Studies of Stingless Bees from the East Coast of Peninsular Malaysia. In Pot-Pollen in Stingless Bee Melittology; Springer: Berlin/Heidelberg, Germany, 2018; pp. 77–88. [Google Scholar]
  2. Michener, C.D. The Bees of the World, 2nd ed.; JHU Press: Baltimore, MA, USA, 2007. [Google Scholar]
  3. Jaapar, M.F.; Halim, M.; Mispan, M.R.; Jajuli, R.; Saranum, M.M.; Zainuddin, M.Y.; Ghazi, R.; Abd Ghani, I. The diversity and abundance of stingless bee (Hymenoptera: Meliponini) in Peninsular Malaysia. Adv. Environ. Biol. 2016, 10, 1–8. [Google Scholar]
  4. Roubik, D.W.; Momose, K.; Inoue, T.; Nagamitsu, T. Preference in flower visits and partitioning in pollen diets of stingless bees in an Asian tropical rain forest. Popul. Ecol. 1999, 41, 195–202. [Google Scholar] [CrossRef]
  5. Jaapar, M.F.; Jajuli, R. Lebah Kelulut Malaysia: Biologi dan Penternakan; Penerbit MARDI: Serdang, Malaysia, 2015. [Google Scholar]
  6. Rattanawannee, A.; Duangphakdee, O. Southeast Asian meliponiculture for sustainable livelihood. In Modern Beekeeping: Bases for Sustainable Production; Ranz, R.E.R., Ed.; IntechOpen: London, UK, 2020; pp. 1–17. [Google Scholar]
  7. Vit, P.; Pedro, S.R.; Roubik, D. (Eds.) Pot-Honey: A Legacy of Stingless Bees; Springer Science & Business Media: New York, NY, USA, 2013. [Google Scholar]
  8. Badrulhisham, N.S.R.; Ab Hamid, S.N.P.; Ismail, M.A.H.; Yong, Y.K.; Zakuan, N.M.; Harith, H.H.; Saidi, H.I.; Nurdin, A. Harvested locations influence the total phenolic content, antioxidant levels, cytotoxic, and anti-inflammatory activities of stingless bee honey. J. Asia-Pac. Èntomol. 2020, 23, 950–956. [Google Scholar] [CrossRef]
  9. Shankar, U.; Abrol, D.P.; Singh, A.K. Plants for bees Pongamia pinnata (L.). Pierre. J. Palynol. 2017, 53, 133–137. [Google Scholar]
  10. Afik, O.; Delaplane, K.S.; Shafir, S.; Moo-Valle, H.; Javier, J.; Quezada-Euán, J.J. Nectar minerals as regulators of flower vis-itation in stingless bees and nectar hoarding wasps. J. Chem. Ecol. 2014, 40, 476–483. [Google Scholar] [CrossRef] [PubMed]
  11. Absy, M.L.; Rech, A.R.; Ferreira, M.G. Pollen collected by stingless bees: A contribution to understanding Amazonian bio-diversity. In Pot-Pollen in Stingless Bee Melittology; Vit, P., Pedro, S., Roubik, D., Eds.; Springer: Berlin/Heidelberg, Germany, 2018; pp. 29–46. [Google Scholar]
  12. Jaapar, M.F.; Nasarodin, N.S.; Reward, N.F.; Jajuli, R.; Abd Ghani, I. Notes on resin collected by stingless bees in Taman Tropika Tasik Kenyir, Terengganu, Malaysia. Serangga 2019, 24, 81–89. [Google Scholar]
  13. Dos Santos, W.E.; Carneiro, L.T.; Alves, A.C.F.; Creão-Duarte, A.J.; Martins, C.F. Stingless Bees (Hymenoptera: Apidae: Meliponini) Attracted to Animal Carcasses in the Brazilian Dry Forest and Implications for Forensic Entomology. Sociobiology 2014, 61, 490–493. [Google Scholar] [CrossRef] [Green Version]
  14. Greco, M.K.; Spooner-Hart, R.N.; Beattie, A.G.A.C.; Barchia, I.; Holford, P. Australian stingless bees improve greenhouse Capsicum production. J. Apic. Res. 2011, 50, 102–115. [Google Scholar] [CrossRef]
  15. Jaffé, R.; Pope, N.; Carvalho, A.T.; Maia, U.M.; Blochtein, B.; De Carvalho, C.A.L.; Carvalho-Zilse, G.A.; Freitas, B.M.; Menezes, C.; Ribeiro, M.D.F.; et al. Bees for Development: Brazilian Survey Reveals How to Optimize Stingless Beekeeping. PLoS ONE 2015, 10, e0121157. [Google Scholar] [CrossRef] [PubMed]
  16. Xiong, X.; Yao, L.; Ying, X.; Lu, L.; Guardone, L.; Armani, A.; Guidi, A.; Xiong, X. Multiple fish species identified from China’s roasted Xue Yu fillet products using DNA and mini-DNA barcoding: Implications on human health and marine sustainability. Food Control. 2018, 88, 123–130. [Google Scholar] [CrossRef]
  17. Talebi, M.; Talebi, M.; Farkhondeh, T.; Samarghandian, S. Molecular mechanism-based therapeutic properties of honey. Biomed. Pharmacother. 2020, 130, 110590. [Google Scholar] [CrossRef] [PubMed]
  18. Ramírez, S.R.; Nieh, J.C.; Quental, T.B.; Roubik, D.W.; Imperatriz-Fonseca, V.L.; Pierce, N.E. A molecular phylogeny of the stingless bee genus Melipona (Hymenoptera: Apidae). Mol. Phylogenet. Evol. 2010, 56, 519–525. [Google Scholar] [CrossRef]
  19. Joseph, S.; Murphy, D.J.; Miller, J.T.; Bhave, M. Application of Molecular Markers for Identification of Potential Salt Tolerant Plant Species for use in Agroforestry and Saline Land Reclamation. APCBEE Procedia 2013, 5, 514–519. [Google Scholar] [CrossRef]
  20. Taberlet, P.; Coissac, E.; Pompanon, F.; Gielly, L.; Miquel, C.; Valentini, A.; Vermat, T.; Corthier, G.; Brochmann, C.; Willerslev, E. Power and limitations of the chloroplast trn L (UAA) intron for plant DNA barcoding. Nucleic Acids Res. 2007, 35, e14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Edgar, R.C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010, 26, 2460–2461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Caporaso, J.G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F.D.; Costello, E.K.; Fierer, N.; Gonzalez, P.A.; Goodrich, J.K.; Gordon, J.I.; et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods. 2010, 7, 335–336. [Google Scholar] [CrossRef] [Green Version]
  23. Cattley, S.; Arthur, J.W. BioManager: The use of a bioinformatics web application as a teaching tool in undergraduate bioinformatics training. Brief Bioinform. 2007, 8, 457–465. [Google Scholar] [CrossRef] [Green Version]
  24. R Core Team. R: A Language and Environment for Statistical Computing. 2018. Available online: https://www.R-project.org/ (accessed on 1 January 2021).
  25. Hawkins, J.; De Vere, N.; Griffith, A.; Ford, C.R.; Allainguillaume, J.; Hegarty, M.J.; Baillie, L.; Adams-Groom, B. Using DNA metabarcoding to identify the floral composition of honey: A new tool for investigating honey bee foraging preferences. PLoS ONE 2015, 10, e0134735. [Google Scholar] [CrossRef] [Green Version]
  26. Selvaraju, K.; Vikram, P.; Soon, J.M.; Krishnan, K.T.; Mohammed, A. Melissopalynological, physicochemical and antioxidant properties of honey from West Coast of Malaysia. J. Food Sci. Technol. 2019, 56, 2508–2521. [Google Scholar] [CrossRef] [Green Version]
  27. Ghazi, R. Melissopalynology and Foraging Activity of the Stingless Bees, Heterotrigona itama (Hymenoptera: Apidae) at Taman Tropika Kenyir, Terengganu. Master’s Thesis, Universiti Malaysia Terengganu, Terengganu, Malaysia, 2015. [Google Scholar]
  28. Kuhn-Neto, B.; Contrera, F.A.L.; Castro, M.S.; Nieh, J.C. Long distance foraging and recruitment by a stingless bee, Melipona mandacaia. Apidologie 2009, 40, 472–480. [Google Scholar] [CrossRef] [Green Version]
  29. Leonhardt, S.D.; Dworschak, K.; Eltz, T.; Blüthgen, N. Foraging loads of stingless bees and utilisation of stored nectar for pollen harvesting. Apidologie 2007, 38, 125–135. [Google Scholar] [CrossRef] [Green Version]
  30. Boongird, S.; Michener, C.D. Pollen and Propolis Collecting by Male Stingless Bees (Hymenoptera: Apidae). J. Kans. Èntomol. Soc. 2010, 83, 47–50. [Google Scholar] [CrossRef]
  31. Tanaka, K.; Nozaki, A.; Nakadai, H.; Shiwa, Y.; Shimizu-Kadota, M. Using pollen DNA metabarcoding to profile nectar sources of urban beekeeping in Kōtō-ku, Tokyo. BMC Res. Notes 2020, 13, 1–6. [Google Scholar] [CrossRef] [PubMed]
  32. Beaman, J. Mount Kinabalu: Hotspot of plant diversity in Borneo. In Plant Diversity and Complexity Patterns: Local, Regional and Global Dimensions. Proceedings of an International Symposium Held at the Royal Danish Academy of Sciences and Letters in Copenhagen, Denmark, 25–28 May 2003; Det Kongelige Danske Videnskabernes Selskab: Copenhagen, Denmark, 2005; Volume 55, pp. 103–127. [Google Scholar]
  33. Saw, L.; Chung, R. The flora of Malaysia projects. Rodriguésia 2015, 66, 947–960. [Google Scholar] [CrossRef] [Green Version]
  34. Pereira, J.T.; Sugau, J.B.; Nilus, R.; Khoo, E.; Tze Lui, S.T.; Leung, L.Y.; Damit, A.; Mustapeng, A.M.A.; Kiat, H.P.; Tanggaraju, S.; et al. Conservation efforts towards contributing to the plant conservation strategies in Sabah, Malaysia. In Proceedings of the 56th Annual Meeting of the Association of Tropical Biology and Conservation, Antananarivo, MA, USA, 30 July–3 August 2019. [Google Scholar]
  35. Jaafar, W.S.W.M.; Maulud, K.N.A.; Kamarulzaman, A.M.M.; Raihan, A.; Sah, S.M.; Ahmad, A.; Saad, S.N.M.; Azmi, A.T.M.; Syukri, N.K.A.J.; Khan, W.R. The Influence of Deforestation on Land Surface Temperature—A Case Study of Perak and Kedah, Malaysia. Forests 2020, 11, 670. [Google Scholar] [CrossRef]
  36. Jenking, C.M.; Boyce, A. Phylum Spermatophyta—Introduction. In The Diversity of Life; Springer: Berlin/Heidelberg, Germany, 1979; pp. 129–131. [Google Scholar]
  37. Van der Sluijs, J.P.; Vaage, N. Pollinators and global food security: The need for holistic global stewardship. Food Ethics 2016, 1, 75–91. [Google Scholar] [CrossRef] [Green Version]
  38. Hung, K.-L.J.; Kingston, J.M.; Albrecht, M.; Holway, D.A.; Kohn, J.R. The worldwide importance of honey bees as pollinators in natural habitats. Proc. R. Soc. B Boil. Sci. 2018, 285, 20172140. [Google Scholar] [CrossRef] [Green Version]
  39. Takhtajan, A. Phylum Magnoliophyta (flowering plants). In Flowering Plants; Springer: Berlin/Heidelberg, Germany, 2009; pp. 1–6. [Google Scholar]
  40. Amolin, A.V.; Ogol, I.N. Trophic relations of wasps (Hymenoptera: Vespidae) with flowering plants (Magnoliophyta) in Donbass. Euroasian Èntomol. J. 2019, 18, 365–376. [Google Scholar] [CrossRef] [Green Version]
  41. Lucek, K.; Galli, A.; Gurten, S.; Hohmann, N.; Maccagni, A.; Patsiou, T.; Willi, Y. Metabarcoding of honey to assess differences in plant-pollinator interactions between urban and non-urban sites. Apidologie 2019, 50, 317–329. [Google Scholar] [CrossRef]
  42. Dow, R.A.; Ngiam, R.W.J.; Ahmad, R. Odonata of Maludam National Park, Sarawak, Malaysia. J. Threat. Taxa 2015, 7, 6764–6773. [Google Scholar] [CrossRef]
  43. Wong, J. Status of Peat Swamp Forest of Sarawak; Sarawak Forestry Corporation/Alterra Wageningen UR: Kuching, Malaysia, 2003. [Google Scholar]
  44. Hazebroek, H.P. National Parks of Sarawak; Natural History Pub. (Borneo): Kota Kinabalu, Malaysia, 2000. [Google Scholar]
  45. Sanz, M.; Gonzalez, M.; de Lorenzo, C.; Sanz, J.; Martínez-CastroI, I. A contribution to the differentiation between nectar honey and honeydew honey. Food Chem. 2005, 91, 313–317. [Google Scholar] [CrossRef]
  46. Fiala, B.; Maschwitz, U. Extrafloral nectaries in the genus Macaranga (Euphorbiaceae) in Malaysia: Comparative studies of their possible significance as predispositions for myrmecophytism. Biol. J. Linn. Soc. 1991, 44, 287–305. [Google Scholar] [CrossRef] [Green Version]
  47. Heckroth, H.-P.; Fiala, B.; Gullan, P.J.; Idris, A.H.; Maschwitz, U. The soft scale (Coccidae) associates of Malaysian ant-plants. J. Trop. Ecol. 1998, 14, 427–443. [Google Scholar] [CrossRef]
  48. Bahri, S. Pollen profile by stingless bee (Heterotrigona itama) reared in rubber smallholding environment at Tepoh, Terengganu. Malays. J. Microsc. 2018, 14, 38–54. [Google Scholar]
  49. Jaapar, M.F.; Jajuli, R.; Mispan, M.R.; Ghani, I.A. Foraging behavior of stingless bee Heterotrigona itama (Cockerell, 1918) (Hymenoptera: Apidae: Meliponini). AIP Conf. Proc. 2018, 1940, 020037. [Google Scholar] [CrossRef]
  50. Asma, S.S.; Adam, N.A.; Najwa, M.Y.S.; Asilah, N.; Syukri, T.S.; Asiah, W.N. Foraging behavior of stingless bee (Geniotrigona thoracica and Heterotrigona itama) on star fruit trees (Averrhoa carambola L.). IOP Conf. Ser. Earth Environ. Sci. 2019, 383, 012011. [Google Scholar] [CrossRef]
  51. Wallace, H.M.; Lee, D.J. Resin-foraging by colonies of Trigona sapiens and T. hockingsi (Hymenoptera: Apidae, Meliponini) and consequent seed dispersal of Corymbia torelliana (Myrtaceae). Apidologie 2010, 41, 428–435. [Google Scholar] [CrossRef] [Green Version]
  52. Bill, N.M. Palynology of Stingless Bee, Tetragonula laeviceps (Hymenoptera: Meliponinae). Master’s Thesis, Faculty of Agro-Based Industry, Universiti Malaysia Kelantan, Kelantan, Malaysia, 2016. [Google Scholar]
  53. Zhi-Shu, X.I.; Yu-Shan, W.A.; Zhi-Bin, Z.H. Seed bank and the factors influencing it for three Fagaceae species in Dujiangyan Region, Sichuan. Biodivers. Sci. 2001, 9, 373–381. [Google Scholar]
  54. Engel, M.S.; Dingemans-Bakels, F. Nectar and pollen resources for stingless bees (Meliponinae, Hymenoptera) in Surinam (South America). Apidologie 1980, 11, 341–350. [Google Scholar] [CrossRef] [Green Version]
  55. Vongsak, B.; Kongkiatpaiboon, S.; Jaisamut, S.; Machana, S.; Pattarapanich, C. In vitro alpha glucosidase inhibition and free-radical scavenging activity of propolis from Thai stingless bees in mangosteen orchard. Rev. Bras. Farm. 2015, 25, 445–450. [Google Scholar] [CrossRef] [Green Version]
  56. Azuan, A.; Maya Izar, K.; Salma, I.; Raziah, M.L. Distribution and diversity of Garcinia species in Malaysia. In Proceedings of the Biodiversity and National Development Achievements, Opportunities and Challenges Conference, Kuala Lumpur, Malaysia, 28–30 May 2008. [Google Scholar]
  57. Slaa, E.J.; Tack, A.J.; Sommeijer, M.J. The effect of intrinsic and extrinsic factors on flower constancy in stingless bees. Apidologie 2003, 34, 457–468. [Google Scholar] [CrossRef] [Green Version]
  58. Momose, K.; Yumoto, T.; Nagamitsu, T.; Kato, M.; Nagamasu, H.; Sakai, S.; Harrison, R.D.; Itioka, T.; Hamid, A.A.; Inoue, T. Pollination biology in a lowland dipterocarp forest in Sarawak, Malaysia. I. Characteristics of the plant-pollinator community in a lowland dipterocarp forest. Am. J. Bot. 1998, 85, 1477–1501. [Google Scholar] [CrossRef] [PubMed]
  59. Wayo, K.; Phankaew, C.; Stewart, A.B.; Bumrungsri, S. Bees are supplementary pollinators of self-compatible chiropterophilous durian. J. Trop. Ecol. 2018, 34, 41–52. [Google Scholar] [CrossRef]
Figure 1. Collection sites of Heterotrigona itama samples throughout Malaysia.
Figure 1. Collection sites of Heterotrigona itama samples throughout Malaysia.
Insects 12 00205 g001
Figure 2. The rarefaction curve of the trnL gene sequence for different localities of H. itama calculated for operational taxonomic units (OTUs) at 97% similarity.
Figure 2. The rarefaction curve of the trnL gene sequence for different localities of H. itama calculated for operational taxonomic units (OTUs) at 97% similarity.
Insects 12 00205 g002
Figure 3. Families of plant species based on trnL result.
Figure 3. Families of plant species based on trnL result.
Insects 12 00205 g003
Figure 4. Heatmap at species level. The 30 most abundant genera were used in hierarchical clustering. A darker color indicates a more dominant genus.
Figure 4. Heatmap at species level. The 30 most abundant genera were used in hierarchical clustering. A darker color indicates a more dominant genus.
Insects 12 00205 g004
Figure 5. Relative frequency of plant species from different localities.
Figure 5. Relative frequency of plant species from different localities.
Insects 12 00205 g005
Figure 6. Venn diagrams summarizing the number of shared plant species among different locations. (a) Shared plant species between Sarawak and Peninsular Malaysia (Terengganu, Kedah, Perak, Perlis, Selangor). (b) Shared plant species between Sabah (Sabah Telipok, Sabah TE, and Sabah Tuaran) and Sarawak (Sarawak Beladin, Pool Sarawak Bintulu, and Sarawak Kuching).
Figure 6. Venn diagrams summarizing the number of shared plant species among different locations. (a) Shared plant species between Sarawak and Peninsular Malaysia (Terengganu, Kedah, Perak, Perlis, Selangor). (b) Shared plant species between Sabah (Sabah Telipok, Sabah TE, and Sabah Tuaran) and Sarawak (Sarawak Beladin, Pool Sarawak Bintulu, and Sarawak Kuching).
Insects 12 00205 g006
Figure 7. Phylogenetic tree dendrogram based on trnL using Bray–Curtis distance (at family level).
Figure 7. Phylogenetic tree dendrogram based on trnL using Bray–Curtis distance (at family level).
Insects 12 00205 g007
Table 1. List of samples used for metabarcoding analysis.
Table 1. List of samples used for metabarcoding analysis.
Sample CodeLocationGridEcosystem
Sabah TEBorneo. Sabah: Lahad Datu5.01248, 118.0970218Crop field adjacent to secondary forest
TerengganuPeninsular Malaysia. Terengganu: Jerangau5.0655177, 102.8652274Crop field adjacent to secondary forest
Sabah TelipokBorneo. Sabah: Telipok6.1732126, 116.2984958Crop field adjacent to hill forest
PerakPeninsular Malaysia. Perak: Lenggong5.1662481, 100.8801384Crop field
KedahPeninsular Malaysia. Kedah: Pokok Sena6.126510, 100.559629Crop field
Sarawak KuchingBorneo. Sarawak: Lundu1.758078, 109.806310Crop field adjacent to coastal area
Pool Sarawak BintuluBorneo. Sarawak: MARDI Bintulu3.358613, 113.430816Crop field adjacent to oil palm estate
PerlisPeninsular Malaysia. Perlis: Chuping6.427755, 100.304074Crop field
Sarawak BeladinBorneo. Sarawak: Beladin1.6273364, 111.1958089Crop field
SelangorPeninsular Malaysia. Selangor: MARDI Serdang2.990892, 101.702434Crop field
Sabah TuaranBornoe. Sabah: Tuaran6.0523621, 116.1788793Crop field adjacent to river
Sarawak MiriBorneo. Sarawak: Miri4.102966, 113.853467Crop field adjacent to coastal area
Table 2. The numbers of effective trnL gene sequences. Number of observed OTUs, alpha diversity indices for the plant DNA from 12 localities of H. itama.
Table 2. The numbers of effective trnL gene sequences. Number of observed OTUs, alpha diversity indices for the plant DNA from 12 localities of H. itama.
LocalitiesSequencesOTUsChao1Shannon–WeinerSimpsonEvenness
Sabah TE70,0686383.0002.7320.90023.056
Terengganu34,2104253.3752.1260.82619.769
Sabah Telipok50,1736287.5002.3720.84126.144
Perak52,2163259.5001.1470.56427.901
Kedah54,9373742.6001.6720.73922.126
Sarawak Kuching20,6424554.752.0230.85122.244
Pool Sarawak Bintulu74483542.0001.9390.79418.047
Perlis56,6122934.0001.0440.47627.786
Sarawak Beladin55,7143948.0001.4940.54126.098
Selangor47,5321922.3300.6450.24329.470
Sabah Tuaran44,0603653.5002.1300.84016.899
Sarawak Miri33,7864886.5001.8570.71525.844
Total527,398487
Table 3. Relative abundance at the phylum level of plant communities from different localities.
Table 3. Relative abundance at the phylum level of plant communities from different localities.
PhylumLocalities %
Tracheophyta4.72
Anthophyta1.97
Magnoliophyta32.39
Spermatophyta55.95
Embryophyta0.13
Unknown4.85
Table 4. Relative abundance at the phylum level of plant communities in Heterotrigona itama from different localities.
Table 4. Relative abundance at the phylum level of plant communities in Heterotrigona itama from different localities.
PhylumBorneo (%)Peninsular Malaysia (%)
Sabah TESabah TelipokSabah TuaranSarawak KuchingPool Sarawak BintuluSarawak BeladinSarawak MiriTerengganu PerakKedahPerlisSelangor
Tracheophyta6.662.811815.71.60.119.398.91.1140.90.080
Anthophyta6.3210.560.010.0101.130.010.010.010.30.170.05
Magnoliophyta41.9257.3738.3394.0978.722.432.52.934.75.63.6
Spermatophyta42.727.943.527.4813.417.152.343.295.863.294.193
Embryophyta000001.16000000
Unknown2.31.2530.01417.780.81.6415.715.20.140.70.013.2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fahimee, J.; Badrulisham, A.S.; Zulidzham, M.S.; Reward, N.F.; Muzammil, N.; Jajuli, R.; Md-Zain, B.M.; Yaakop, S. Metabarcoding in Diet Assessment of Heterotrigona itama Based on trnL Marker towards Domestication Program. Insects 2021, 12, 205. https://0-doi-org.brum.beds.ac.uk/10.3390/insects12030205

AMA Style

Fahimee J, Badrulisham AS, Zulidzham MS, Reward NF, Muzammil N, Jajuli R, Md-Zain BM, Yaakop S. Metabarcoding in Diet Assessment of Heterotrigona itama Based on trnL Marker towards Domestication Program. Insects. 2021; 12(3):205. https://0-doi-org.brum.beds.ac.uk/10.3390/insects12030205

Chicago/Turabian Style

Fahimee, Jaapar, Aqilah Sakinah Badrulisham, Mohd Sani Zulidzham, Nurul Farisa Reward, Nizar Muzammil, Rosliza Jajuli, Badrul Munir Md-Zain, and Salmah Yaakop. 2021. "Metabarcoding in Diet Assessment of Heterotrigona itama Based on trnL Marker towards Domestication Program" Insects 12, no. 3: 205. https://0-doi-org.brum.beds.ac.uk/10.3390/insects12030205

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop