Japanese Honeybees (Apis cerana japonica Radoszkowski, 1877) May Be Resilient to Land Use Change

Simple Summary Pollinators are threatened globally by growing urban sprawl and agriculture. The Western Honeybee (Apis mellifera) readily adapts to whatever food is available, so people have made it the most widely distributed pollinator across the world. Previous research has suggested that the Western Honeybee may be less resilient to land use change outside of its natural range. This study examines a different honeybee species—the Japanese Honeybee (Apis cerana japonica). Unlike the Western Honeybee, this species is found almost exclusively in its natural range in Japan. Consequently, it may be better adapted to its local food sources and therefore more resilient. Working in southern Japan, in the Nagasaki and Saga prefectures, we looked at the nectar and pollen that the Japanese Honeybee feeds on. Their food intake was then examined in relation to local land use composition. We found minimal impact of increasing urban sprawl on the forage of the Japanese Honeybee. This goes against previous studies on the Western Honeybee elsewhere in the world. Though in need of a direct comparison with Western Honeybee, these preliminary results could be due to differences in urban green infrastructure in Japan, or due to an adaptation by the Japanese honeybee to its surroundings. Abstract Pollinators are being threatened globally by urbanisation and agricultural intensification, driven by a growing human population. Understanding these impacts on landscapes and pollinators is critical to ensuring a robust pollination system. Remote sensing data on land use attributes have previously linked honeybee nutrition to land use in the Western Honeybee (Apis mellifera L.). Here, we instead focus on the less commonly studied Apis cerana japonica—the Japanese Honeybee. Our study presents preliminary data comparing forage (honey and pollen) with land use across a rural-urban gradient from 22 sites in Kyushu, southern Japan. Honey samples were collected from hives between June 2018 and August 2019. Pollen were collected and biotyped from hives in urban and rural locations (n = 4). Previous studies of honey show substantial variation in monosaccharide content. Our analysis of A. cerana japonica honey found very little variation in glucose and fructose (which accounted for 97% of monosaccharides), despite substantial differences in surrounding forage composition. As expected, we observed temporal variation in pollen foraged by A. cerana japonica, likely dependent on flowering phenology. These preliminary results suggest that the forage and nutrition of A. cerana japonica may not be negatively affected by urban land use. This highlights the need for further comparative studies between A. cerana japonica and A. mellifera as it could suggest a resilience in pollinators foraging in their native range.


Introduction
Globally, land use change is driving biodiversity loss. Primarily driven by anthropogenic uses, such as agricultural intensification and urban development, these losses the Western honeybee has adapted to forage more on "mass-flowering plants" than the Japanese honeybee. Beekeepers and academics have observed that the Japanese honeybee instead forages on more diverse flower sources, collecting nectar and pollen from anything available, rather than through the "flower constancy" behaviour observed in Western honeybee [34,35]. The beekeepers believe that "if nectar is very scarce (like in urban area), cerana can more successfully find small patches on which to survive (maybe in a small garden)" [36]. Observations by this group comparing Western honeybee and Japanese honeybee when kept in the same apiary have led the beekeepers to believe that "mellifera is not good at collecting nectar in summer compared with cerana", especially in urban environments with less spatially extensive flower patches [36]. It has recently been identified that Western honeybee are more susceptible than Japanese honeybee to predation from various native hornets in Japan [29,30].
Beekeeping is growing in Japan, as more people take it up as a hobby or a business interest. Investigating the threats to pollinators and their ecosystem services are equally increasing in importance. Identifying land uses and plant species that support A. cerana japonica may be key in maintaining the success of both urban and rural commercial and hobbyist beekeepers, as their numbers continue to grow.

Aims and Scope
The sugar composition of honey is directly linked to the flowers that have been foraged by bees to make it. Phytochemicals derived from the flowers convey unique biochemical properties to honey, which has led to marketization of particular honeys for their medicinal properties (such as "Manuka") [37][38][39]. The glucose/fructose ratio of honey can impact its palatability for humans [40,41], and higher glucose levels contribute to the production of hydrogen peroxide, an important antimicrobial compound found in honey [42,43].
This study aims to investigate pollen and nectar foraging by A. cerana japonica, studying urban and rural populations maintained by hobbyists in Nagasaki and Saga prefectures, on Japan's southern island of Kyushu. Here, beekeeping is practiced by a small beekeeping community, largely producing honey part-time for personal use and sale, with a few full-time commercial producers. A. cerana japonica are almost exclusively chosen by hobbyist beekeepers in Nagasaki (though A. mellifera continue to be favoured by commercial beekeepers). The factors affecting pollinator health and honey production in Japan are not well studied, so this study focuses on answering the two following key questions: 1.
How does land use in Japan affect the honey produced by A. cerana japonica? 2.
How does the time of year and location affect the pollen collected and honey produced by A. cerana japonica?

Honey Sampling
Honey was collected from 22 hives across Nagasaki and Saga ( Figure 1) between June 2018 and August 2019. A questionnaire was given to the owner of each hive to determine information such as the species of bee, the location of their hive and the environment surrounding the hive (Supplementary Materials 1: Questionnaire).
Honey samples were analysed using high-performance liquid chromatography (HPLC) to determine sugar composition (fructose, glucose and maltose), following methods for honey analysis used in Ouchemoukh et al. [44]. These sugars were selected due to previously observed geographic variation in composition [45].
Crystallised honey samples were heated in a 90 • C water bath until clear and left to cool to room temperature. All samples were then diluted using High Performance Liquid Chromatography (HPLC)-grade water to a ratio of 1 µL honey per 100 mL water. Two repeats of 10 µL solution were run through the HPLC, testing for fructose, glucose and maltose levels. Three repeats per sample were conducted and the means taken.

Land Use Composition
To analyse the correlation between land cover and Japanese honeybee nutrition, data were sourced from the Japanese Ministry of the Environment's Biodiversity Centre (Table  S2). The composition and configuration of different land uses in the surrounding 1, 3 and 5 km radii of each hive were measured using the radius tool in ArcGIS 10.8.1 (ESRI, US).
Radii of 1, 3 and 5 km were chosen as these cover the range of foraging distances travelled by bees from their hive [47]. Land cover classes that accounted for <0.5% of total cover within a buffer zone were excluded from analysis. The dominant land use (the land use contributing the greatest percentage of land cover) and the ratio of urban-to-rural land uses were then calculated using these data. The distance from each hive to the nearest urban area was also calculated based on methods established in Clermont et al. [50].

Statistical Analysis
Honey composition data were analysed to determine inter-hive variance in sugar content. Monosaccharide composition of the honey samples was analysed by Pearson's correlation with three land use factors: distance to urban areas, dominant surrounding land use and ratio of rural-to-urban land use. Critical P-scaling was performed to control Aliquots were analysed on an Agilent Analytical 1200LC HPLC machine (Agilent Systems, UK) using a Thermo Dionex CarboPac PA20 Analytical column, 3 × 150 mm (Thermofisher, UK). A Pulsed Amperometric Detection (PAD) detector was used, and samples were transported in HPLC grade water/200 nM NaOH. Detection peaks were quantified against a dilution series of standards for fructose, glucose and maltose.

Pollen Sampling
Pollen samples were collected at least once every three weeks for nine weeks from 13.6.2019 to 8.8.2019 at two urban (Hives 15, 16) and two rural hives (Hives 11,19) in Nagasaki-ken. Pollen was sampled via pollen trapping, direct pollen basket collection and brood chamber sampling, as these methods have previously been used to sample pollen successfully [46,47]. Samples at all hives were collected from the brood chamber, as this was found to be the quickest and most effective method, as well as arguably the least disruptive. Typhoons during weeks four and seven limited sampling during these periods; one hive was abandoned by the colony in week eight (detailed information available in Table S1). Biotyping of pollen grains was used to provide a descriptive measure of the diversity of foraged pollen in samples [48]. Pollen samples were purified via acetolysis using methods based on Jones [49]; each sample was imaged five times with a digital microscope. Pollen grains were then counted and identified into biotypes based on physical properties (Supplementary Appendix 1: Pollen Identification), due to low sample site replication, no further analysis was performed on the pollen data.

Land Use Composition
To analyse the correlation between land cover and Japanese honeybee nutrition, data were sourced from the Japanese Ministry of the Environment's Biodiversity Centre (Table S2). The composition and configuration of different land uses in the surrounding 1, 3 and 5 km radii of each hive were measured using the radius tool in ArcGIS 10.8.1 (ESRI, US).
Radii of 1, 3 and 5 km were chosen as these cover the range of foraging distances travelled by bees from their hive [47]. Land cover classes that accounted for <0.5% of total cover within a buffer zone were excluded from analysis. The dominant land use (the land use contributing the greatest percentage of land cover) and the ratio of urban-to-rural land uses were then calculated using these data. The distance from each hive to the nearest urban area was also calculated based on methods established in Clermont et al. [50].

Statistical Analysis
Honey composition data were analysed to determine inter-hive variance in sugar content. Monosaccharide composition of the honey samples was analysed by Pearson's correlation with three land use factors: distance to urban areas, dominant surrounding land use and ratio of rural-to-urban land use. Critical P-scaling was performed to control false discovery rate (FDR) on these multiple land use analyses [51], for these analyses critical P was set at α = 0.0115.
Pollen count data were not analysed due to low hive level replication. The number of observed biotypes (as a rough approximation of pollen species richness), and sample date and hive are presented in Supplementary Materials. All other analyses were all carried out in R statistical software version 4.0.5 [52].

Honey Sugar Composition
Honey samples from 22 different hives across Japan, were analysed for fructose, glucose and maltose content. Glucose accounted for the majority of mono-saccharides within the honey samples, accounting for 59.3 ± 0.4% (mean ± S.E.); followed by fructose at 38.1 ± 0.5%, and maltose sugars being present in trace amounts at 2.5 ± 0.2%. The sugar ratios of honey samples did not change significantly between hives and apiary location (H = 21, df = 21, p = 0.459; Figure 2).
Although fructose and glucose content remained consistent between hives, significant inter-hive variance was observed in maltose content (F = 0.465, df = 21, p = 0.029). The distance to the nearest urban area had no effect on the proportions of sugars in honey samples Honey samples from 22 different hives across Japan, were analysed for fructose, glu cose and maltose content. Glucose accounted for the majority of mono-saccharides within the honey samples, accounting for 59.3 ± 0.4% (mean ± S.E.); followed by fructose at 38. ± 0.5%, and maltose sugars being present in trace amounts at 2.5 ± 0.2%. The sugar ratio of honey samples did not change significantly between hives and apiary location (H = 21 df = 21, p = 0.459; Figure 2).

Figure 2.
Proportion of fructose, maltose and glucose in honey samples collected from 22 Apis spe cies hives across Kyushu, Japan. Yellow dots represent hives where pie chart could not be placed directly in the correct location. Where multiple hives had the same location, pie charts are labelled with hive number. Sugar proportions were calculated using HPLC. Clustered hives located in Na gasaki, within the black rectangle, can be found in Figure S1.
Although fructose and glucose content remained consistent between hives, signifi cant inter-hive variance was observed in maltose content (F = 0.465, df = 21, p = 0.029). Th distance to the nearest urban area had no effect on the proportions of sugars in honey Figure 2. Proportion of fructose, maltose and glucose in honey samples collected from 22 Apis species hives across Kyushu, Japan. Yellow dots represent hives where pie chart could not be placed directly in the correct location. Where multiple hives had the same location, pie charts are labelled with hive number. Sugar proportions were calculated using HPLC. Clustered hives located in Nagasaki, within the black rectangle, can be found in Figure S2.

Pollen Biotype Composition
In the process of biotyping, 262 microscope images were counted across 131 samples from four A. cerana japonica hives collected between June and September 2019. Across these samples, 50 biotypes were identified (Supplementary Materials: Pollen). Biotypes 1, 3, 10, 14, 17 and 20 had a total abundance greater than 1000 grains across all sample images, and thus were deemed 'dominant' biotypes. A description of these biotypes can be found in Table 1. Statistical analyses of pollen biotypes use these groups throughout, as well as 'others' (combined totals from biotypes 2, 4-9, 11-13, 15, 16, 18, 19 and 20-52).

Land Use Composition
Within a 1 km radius of the 22 hives, eight land types were dominant (Table S2). Within a 3 km radius, this was six land types, and in a 5 km radius, this was five land types. On average, within a 1 km radius around each hive, rural land use was dominant, covering an average of 80.0 ± 6.7% of all land. This land cover varied between hives (Figure 3 and Figure S2), but this variation was not reflected in fructose or glucose content of honey. The maltose content of honey increased as rural land cover within 1 km radius of the hive increased, but not significantly under P-scaling to control FDR (Table 2). Table 1. Dominant biotype classifications and descriptions used to identify pollen grains found in samples collected from five Apis cerana japonica hives in Japan, between June and September 2019. Potential family and species names are also given where possible, based on pollen samples collected directly from flowering plants and identified using a palynomorph guide of Japanese flora by Shimakura [53]. Reference images are approximately relatively sized, as indicated by scale bar (50 µm).

Pollen Biotype Composition
In the process of biotyping, 262 microscope images were counted across 131 samples from four A. cerana japonica hives collected between June and September 2019. Across these samples, 50 biotypes were identified (Supplementary Materials: Pollen). Biotypes 1, 3, 10, 14, 17 and 20 had a total abundance greater than 1000 grains across all sample images, and thus were deemed 'dominant' biotypes. A description of these biotypes can be found in Table 1. Statistical analyses of pollen biotypes use these groups throughout, as well as 'others' (combined totals from biotypes 2, 4-9, 11-13, 15, 16, 18, 19 and 20-52). Table 1. Dominant biotype classifications and descriptions used to identify pollen grains found in samples collected from five Apis cerana japonica hives in Japan, between June and September 2019. Potential family and species names are also given where possible, based on pollen samples collected directly from flowering plants and identified using a palynomorph guide of Japanese flora by Shimakura [53]. Reference images are approximately relatively sized, as indicated by scale bar (50µm).

Biotype
Possible Species and Family Name (s) Description Reference Image

Pollen Biotype Composition
In the process of biotyping, 262 microscope images were counted across 131 samples from four A. cerana japonica hives collected between June and September 2019. Across these samples, 50 biotypes were identified (Supplementary Materials: Pollen). Biotypes 1, 3, 10, 14, 17 and 20 had a total abundance greater than 1000 grains across all sample images, and thus were deemed 'dominant' biotypes. A description of these biotypes can be found in Table 1. Statistical analyses of pollen biotypes use these groups throughout, as well as 'others' (combined totals from biotypes 2, 4-9, 11-13, 15, 16, 18, 19 and 20-52). Table 1. Dominant biotype classifications and descriptions used to identify pollen grains found in samples collected from five Apis cerana japonica hives in Japan, between June and September 2019. Potential family and species names are also given where possible, based on pollen samples collected directly from flowering plants and identified using a palynomorph guide of Japanese flora by Shimakura [53]. Reference images are approximately relatively sized, as indicated by scale bar (50µm).

Biotype
Possible Species and Family Name (s) Description Reference Image

Pollen Biotype Composition
In the process of biotyping, 262 microscope images were counted across 131 samples from four A. cerana japonica hives collected between June and September 2019. Across these samples, 50 biotypes were identified (Supplementary Materials: Pollen). Biotypes 1, 3, 10, 14, 17 and 20 had a total abundance greater than 1000 grains across all sample images, and thus were deemed 'dominant' biotypes. A description of these biotypes can be found in Table 1. Statistical analyses of pollen biotypes use these groups throughout, as well as 'others' (combined totals from biotypes 2, 4-9, 11-13, 15, 16, 18, 19 and 20-52). Table 1. Dominant biotype classifications and descriptions used to identify pollen grains found in samples collected from five Apis cerana japonica hives in Japan, between June and September 2019. Potential family and species names are also given where possible, based on pollen samples collected directly from flowering plants and identified using a palynomorph guide of Japanese flora by Shimakura [53]. Reference images are approximately relatively sized, as indicated by scale bar (50µm).

Biotype
Possible Species and Family Name (s) Description Reference Image Similar to 16, but more indents, elongated. From end on: small-medium, three rounded sides, triangle inside, with points between indents.

Land Use Composition
Within a 1 km radius of the 22 hives, eight land types were dominant (Table S2). Within a 3 km radius, this was six land types, and in a 5 km radius, this was five land types. On average, within a 1 km radius around each hive, rural land use was dominant, covering an average of 80.0 ± 6.7% of all land. This land cover varied between hives (Figures 3 and S2), but this variation was not reflected in fructose or glucose content of honey. The maltose content of honey increased as rural land cover within 1 km radius of the hive increased, but not significantly under P-scaling to control FDR (Table 2).

Discussion
In this study, fructose and glucose sugars made up the majority (an average of 97.5 ± 0.9%) of the three monosaccharide sugars analysed from the honey of the Japanese honeybee. Our study was based on sugar analysis from sites across south Japan (n = 22), as well as pollen collections from a smaller number of sites (n = 4), sampled repeatedly (n = 13 per hive). Our results show minimal variance in the sugar composition across a wide geographical area, and a broad diversity of pollen types associated with these bees.
Despite the small replication size of pollen trapping, the sampling regime we attempted was able to show temporal variation in pollen forage composition. Flowering date is a significant factor impacting pollinator foraging [9,54,55], and our study was no exception. As a result, it is important to focus on ensuring a consistent provision of pollen sources throughout the year, by maintaining heterogeneity of flowering times within vegetation cover [12,48].

The Importance of Honey Sugar Composition
The biochemical properties (i.e., antibacterial activity) of honey varies according to its sources, as do its physical and chemical properties [38,39,56]. These properties are as important to humans for their medical properties as they are to the honeybees for their potential to act as a "bee pharmacy" [57,58]. Variation in these properties has been suggested to be part of a suite of factors negatively impacting honeybee health globally [16,59].
Consequently, flower resources (availability and botanical origin) in the urban environment and subsequent variance in honey sugar composition have the potential to substantially impact bee health. Previous research in the Western honeybee from western urban-rural gradient landscapes has shown that increasing urban land use has a negative impact on the nutrition of these eusocial bees [9][10][11][12][13]. Yet, in our study measures of land use composition, for example rural-to-urban ratio or distance to urban landscapes, did not significantly impact sugar composition. Maltose levels appeared to vary substantially more between hives. Maltose is often added to man-made products, as it is not commonly naturally occurring [60]. However, studies have shown that the sugar composition of honey changes over time, due to various chemical processes enabled by the heat of hives [40,61]. This can lead to the production of maltose from glucose and fructose, meaning that the maltose may have formed during storage [61].
No link between nutrition of the Japanese honeybee and the composition of their environment suggests that Japanese honeybees in Japan may not be perturbed by the negative effects of urbanisation observed in other studies of honeybees. The question remains, as to why: is it land use being less detrimental, or are Japanese honeybee being more resilient to land use change?

Could Urban Land Use Be Less Detrimental?
As observed by beekeepers involved in our study, when comparing Western honeybee with Japanese honeybee, the latter may prefer to forage on a wider array of available forage [36]. The low variance in monosaccharide composition in our study of Japanese honeybee, compared with the significant inter-hive variance in sugar composition from studies of Western honeybee [45], suggests Japanese honeybee may be better at regulating the balance of sugars in their nectar forage intake.
One could suggest this is due to differences in the landscape composition impacting the bees. Within the context of our study (Nagasaki and Saga), the rural-urban matrix of landscapes was significantly variable between sites (Figures S1 and S2). Our findings may show that Japanese honeybee could be successfully exploiting the diverse nectar sources available to them in these environments. This may suggest that Japanese urban landscapes are not as detrimental to Japanese honeybee as western urban environments are to Western honeybee, as found in existing studies that have explored the impact of urban environments on pollen forage and pollinator health [50,62,63]. Few studies have explicitly studied the link between landscape "quality" and nectar foraging as has been attempted in our study. Homogenous landscapes (i.e., agricultural monocultures) have previously been shown to negatively impact the nectar foraging of Western honeybees [64]. The effect of urban environments on nectar foraging remains largely unexplored for Western honeybees however. Furthermore, there is a growing literature on urban green infrastructure that suggests cultural and demographic factors influence the distribution and composition of urban green spaces [65,66].
Elsewhere in the world, where urban land uses dominate, the managed floral composition and overall lack of availability of nectar has been shown to have detrimental effects on honeybee nutrition [50,63]. Previous studies directly comparing the effects of land use composition on honeybee nutrition found significant effects on pollinator foraging and nutrition [9,67,68]. The lack of impact from landscape composition on nutrition of the Japanese honeybee here suggests urban flora may not be perturbing foraging by Japanese honeybee. Cultural, historic and heritage factors that impact the composition of natural elements in urban settings may explain these differences. If urban green spaces are providing a stable and sufficient nectar source for Japanese honeybee [65,69], prefectures like Nagasaki and Saga may become important case studies in urban landscapes that support local bee populations.
The social construct of what constitutes "urban" in Japan is observably distinct to other sites in which pollinator ecology has been researched, in terms of extent and age of urban green spaces, the social expectation for preserving and cultivating these spaces, as well as social engagement with natural history. [65,[70][71][72]. In previous studies, urban environments are often dominated by invasive or alien species, supplanting native flora [73,74]. Cultural factors leading to an emphasis on native flora [65,71], as well as strict plant quarantine and biosecurity measures [75], may have contributed to urban environments in Japan being more "native" than western urbanised landscapes. Emergent analyses of urban green infrastructure suggest remote sensing data may bear out this comparison [65,66,76], but a direct comparative analysis with sufficiently controlled parameters is currently lacking in the literature.

Could Japanese Honeybees Be More Resilient?
The other driver of our observed minimal disturbance of Japanese honeybee within urban environments and agricultural rural environments is that Japanese honeybee itself may be better adapted to these environments. Due to the lack of significant correlation with land use composition on honeybee nutrition, we suggest that in Nagasaki and Saga, Japanese honeybee may in fact be "coping" with the presence and expansion of urbanisation. Previous studies by Garbuzov et al. [69] and Lowenstein et al. [8] suggest that pollinators can adapt successfully to urban living. However, notably, due to a lack of direct comparison with Western honeybee or experiments relocating Japanese honeybee to a non-localised environment, determining this link remains speculative.
Some studies suggested that honeybees may prefer key plant species because of their ability to perceive higher nutritional value in these plants [77,78]. We suggest that it is possible the Japanese honeybee may have co-evolved with their native flora to have these advantageous preferences. Whereas the Western honeybee, being both a global generalist and alien species to this location, is not connected with the local flora via eco-evolutionary interactions and therefore is lacking this foraging and nutritional advantage.
When comparing the biology of Western honeybee with Japanese honeybee, studies of the "waggle dance" properties suggest that the flight capacity of Western honeybee may be up to two times larger than Japanese honeybee [35,79]. The suggestion, that different species of honeybee might possess distinct 'dialects' of the waggle dance, remains controversial. Direct comparisons of waggle dance properties between these two species suggest that their "dialects" (i.e., the forage distance-to-dance duration) scale differently [29,80]. Foraging over larger distances is typically advantageous to generalist foragers like the Western honeybee, though it suggests that Japanese honeybee may be better adapted to smaller areas with forage sources that are more difficult to locate [80].
Denser, more homogeneous and/or larger urban areas may have a significant effect on Japanese honeybee [54]. The area studied in this paper did not have many high population density (Nagasaki: 335 people/km 2 , Saga: 342 people/km 2 ) or large urban areas, so this may be a relevant factor to incorporate in future studies. However, based on the results of this study alone, we found that the proximity of hives to urban areas is not an issue for regulation of honey composition, which means that the Japanese honeybee may be resilient in increasingly urbanised areas. This contradicts research on Western honeybee, which has highlighted the threats of urban and agricultural expansion. Here, we may be observing a difference between Western honeybee and Japanese honeybee, a resilience to land use variance in Japanese honeybee not found in other honeybee species [12,63]. One key factor when considering local environmental impacts on pollinators is activity outside of their native range. Whereas Western honeybee has a cosmopolitan distribution, well outside of its native range [81], Japanese honeybee is almost exclusively present within its native range in Japan. Any apparent resilience of this latter species may be linked with this, highlighting it as an important factor in future studies of pollinator resilience relative to its range/distribution.

Conclusions
Our study found that the sugar composition of honeys collected by Japanese honeybee does not vary significantly between hives across a wide geographical area. Furthermore, no significant correlation was found with land use composition surrounding each hive. Our understanding of the importance of sugar composition for hive health, the impact of botanical origins on sugar composition, and local knowledge on the foraging behaviour of Japanese honeybee leads us to conclude that this particular bee species may be sufficiently able to adapt to local land use conditions. By this, we mean the Japanese honeybee honey forage is not negatively affected by land use, as has been seen in studies of the Western honeybee [64].
Comparing Japanese honeybee to Western honeybee would potentially identify factors allowing the Japanese honeybee to successfully exploit both urban and rural environments in Japan, when the Western honeybee tends to be less successful across the planet. It would also be beneficial to identify local plant species abundances relative to the pollen spectrum foraged by the Japanese honeybee, to maximise benefits of any vegetation conservation efforts.
To discover more of the effects of various external factors on honey composition, more complex compositional analyses, including other sugars such as melezitose and sucrose, as well as other components such as pollen grains could be relevant [44,82,83]. In summary, our study has provided a first-look at the nutrition of A. cerana japonica in its native habitat, which suggests these bees in particular may be resilient to the effects of land use change that have negatively affected other honeybee species elsewhere in the world.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/insects12080685/s1, Supplementary Materials 1: Questionnaire. Appendix 1: Pollen biotype identification; Table S1. Pollen collection dates for Apis cerana japonica hives in Nagasaki-ken, Japan; Table S2. The land types used to classify vegetation cover across Nagasaki and Saga provinces of Japan. Figure S1. Proportion of fructose, maltose and glucose in honey samples collected from five Apis cerana japonica hives in Nagasaki, Japan; Figure S2. Proportion of land uses in a 1 km radius surrounding five Apis cerana japonica hives in Nagasaki, Japan. Supplementary materials: Questionnaire data.
Author Contributions: All authors contributed to conceiving and writing the manuscript. P.D. performed statistical analysis, L.C. and T.O. performed honey and pollen sampling, P.D. and L.C. performed chemical analyses, species identification were performed by L.C. All authors have read and agreed to the published version of the manuscript.
Funding: P.D. is supported by an ESRC Impact Acceleration Account grant awarded by Lancaster University. Fieldwork by L.C. was supported by the Japan Student Services Organisation (JASSO) exchange programme between Lancaster University and Nagasaki University. T.O. acknowledges departmental support from Nagasaki University.

Institutional Review Board Statement: Not applicable.
Data Availability Statement: All data generated or analysed during this study are included in this published article and its supplementary information files.