Reproductive interference between honeybee species in artificial sympatry

Reproductive isolation between closely related species is often incomplete. The Western honeybee, Apis mellifera, and the Eastern hive bee, Apis cerana, have been allopatric for millions of years, but are nonetheless similar in morphology and behaviour. During the last century, the two species were brought into contact anthropogenically, providing potential opportunities for interspecific matings. Hybrids between A. mellifera and A. cerana are inviable, so natural interspecific matings are of concern because they may reduce the viability of A. cerana and A. mellifera populations – two of the world's most important pollinators. We examined the mating behaviour of A. mellifera and A. cerana queens and drones from Caoba Basin, China and Cairns, Australia. Drone mating flight times overlap in both areas. Analysis of the spermathecal contents of queens with species‐specific genetic markers indicated that in Caoba Basin, 14% of A. mellifera queens mated with at least one A. cerana male, but we detected no A. cerana queens that had mated with A. mellifera males. Similarly, in Cairns, no A. cerana queens carried A. mellifera sperm, but one‐third of A. mellifera queens had mated with at least one A. cerana male. No hybrid embryos were detected in eggs laid by interspecifically mated A. mellifera queens in either location. However, A. mellifera queens artificially inseminated with A. cerana sperm produced inviable hybrid eggs or unfertilized drones. This suggests that reproductive interference will impact the viability of honeybee populations wherever A. cerana and A. mellifera are in contact.     

Multiple host shifts by the emerging honeybee parasite,Varroa jacobsoni

Host shifts are a key mechanism of parasite evolution and responsible for the emergence of many economically important pathogens. Varroa destructor has been a major factor in global honeybee (Apis mellifera) declines since shifting hosts from the Asian honeybee (Apis cerana) > 50 years ago. Until recently, only two haplotypes of V. destructor (Korea and Japan) had successfully host shifted to A. mellifera. In 2008, the sister species V. jacobsoni was found for the first time parasitizing A. mellifera in Papua New Guinea (PNG). This recent host shift presents a serious threat to world apiculture but also provides the opportunity to examine host shifting in this system. We used 12 microsatellites to compare genetic variation of V. jacobsoni on A. mellifera in PNG with mites on A. cerana in both PNG and surrounding regions. We identified two distinct lineages of V. jacobsoni reproducing on A. mellifera in PNG. Our analysis indicated independent host shift events have occurred through small numbers of mites shifting from local A. cerana populations. Additional lineages were found in the neighbouring Papua and Solomon Islands that had partially host shifted to A. mellifera, that is producing immature offspring on drone brood only. These mites were likely in transition to full colonization of A. mellifera. Significant population structure between mites on the different hosts suggested host shifted V. jacobsoni populations may not still reproduce on A. cerana, although limited gene flow may exist. Our studies provide further insight into parasite host shift evolution and help characterize this new Varroa mite threat to A. mellifera worldwide.     

European honeybee defense against Japanese yellow hornet using heat generation by bee‐balling behavior

The Japanese honeybee, Apis cerana japonica Radoszkowski, uses unique generation of heat by bee‐balling to defend against, overheat and kill predacious Japanese hornets. We have now observed the European honeybee, Apis mellifera Linnaeus, using similar bee‐balling behavior and heat generation against the Japanese yellow hornet, Vespa simillima xanthoptera Cameron. We monitored temperatures in the center of the bee‐ball and inside thoracic muscles of the captured hornet and found that the thoracic internal temperature (45.8 ± 2.32°C) was higher than that of the bee‐ball (44.0 ± 0.96°C). Although the thoracic temperature of captured hornets rose to the upper lethal level, defending European bees also showed some stinging attempts against the hornet, unlike the sympatric Japanese honeybee, which never stings during bee‐balling. The European honeybee bee‐balling behavior consists of three phases: (i) heating; (ii) heat‐retaining; and (iii) break up. Our results suggest that European honeybees kill hornets by raising the body temperature of hornets rapidly without stinging. The tactics of bee‐balling against hornets are complex and may be performed by extended division of labor.     

Heritable variation in circulating glucocorticoids and endocrine flexibility in a free‐living songbird

Phenotypic flexibility is a central way that organisms cope with challenging and changing environments. As endocrine signals mediate many phenotypic traits, heritable variation in hormone levels, or their context‐dependent flexibility, could present an important target for selection. Several studies have estimated the heritability of circulating glucocorticoid levels under acute stress conditions, but little is known about the potential for either baseline hormone levels or rapid endocrine flexibility to evolve. Here, we assessed the potential for selection to operate on the elevation (circulating hormone levels) and flexibility of glucocorticoid reaction norms to acute restraint stress. Multivariate animal models revealed low but significant heritability in baseline (h² = 0.13–0.14) and stress‐induced glucocorticoids (h² = 0.18), and moderate heritability in glucocorticoid flexibility in response to acute stress (h² = 0.38) in free‐living juvenile tree swallows (Tachycineta bicolor; n = 408). Baseline glucocorticoids were not genetically correlated with either stress‐induced glucocorticoids or glucocorticoid flexibility. These findings indicate that baseline glucocorticoids and the acute stress response are distinct traits that can be independently shaped by selection. Microevolutionary changes that influence the expression or flexibility of these endocrine mediators of phenotype may be an important way that populations adapt to changing environments and novel threats.     

The pheromones of laying workers in two honeybee sister species: <Emphasis Type="Italic">Apis cerana</Emphasis> and <Emphasis Type="Italic">Apis mellifera</Emphasis>

When a honeybee colony loses its queen, workers activate their ovaries and begin to lay eggs. This is accompanied by a shift in their pheromonal bouquet, which becomes more queen like. Workers of the Asian hive bee Apis cerana show unusually high levels of ovary activation and this can be interpreted as evidence for a recent evolutionary arms race between queens and workers over worker reproduction in this species. To further explore this, we compared the rate of pheromonal bouquet change between two honeybee sister species of Apis cerana and Apis mellifera under queenright and queenless conditions. We show that in both species, the pheromonal components HOB, 9-ODA, HVA, 9-HDA, 10-HDAA and 10-HDA have significantly higher amounts in laying workers than in non-laying workers. In the queenright colonies of A. mellifera and A. cerana, the ratios (9-ODA)/(9-ODA + 9-HDA + 10-HDAA + 10-HDA) are not significantly different between the two species, but in queenless A. cerana colonies the ratio is significant higher than in A. mellifera, suggesting that in A. cerana, the workers’ pheromonal bouquet is dominated by the queen compound, 9-ODA. The amount of 9-ODA in laying A. cerana workers increased by over 585% compared with the non-laying workers, that is 6.75 times higher than in A. mellifera where laying workers only had 86% more 9-ODA compared with non-laying workers.     

Herd prevalence of Enterobacteriaceae producing CTX‐M‐type and CMY‐2 β‐lactamases among Japanese dairy farms

Aims

To determine the herd prevalence of Enterobacteriaceae producing CTX‐M‐type extended‐spectrum β‐lactamases (ESBLs) among 381 dairy farms in Japan.

Methods and Results

Between 2007 and 2009, we screened 897 faecal samples using BTB lactose agar plates containing cefotaxime (2 μg ml^?1). Positive isolates were tested using ESBL confirmatory tests, PCR and sequencing for CTX‐M, AmpC, TEM and SHV. The incidence of Enterobacteriaceae producing CTX‐M‐15 (n = 7), CTX‐M‐2 (n = 12), CTX‐M‐14 (n = 3), CMY‐2 (n = 2) or CTX‐M‐15/2/14 and CMY‐2 (n = 4) in bovine faeces was 28/897 (3·1%) faecal samples. These genes had spread to Escherichia coli (n = 23) and three genera of Enterobacteriaceae (n = 5). Herd prevalence was found to be 20/381 (5·2%) dairy farms. The 23 E. coli isolates showed clonal diversity, as assessed by multilocus sequence typing and pulsed‐field gel electrophoresis. The pandemic E. coli strain ST131 producing CTX‐M‐15 or CTX‐M‐27 was not detected.

Conclusions

Three clusters of CTX‐M (CTX‐M‐15, CTX‐M‐2, CTX‐M‐14) had spread among Japanese dairy farms.

Significance and Impact of the Study

This is the first report on the prevalence of multidrug‐resistant CTX‐M‐15–producing E. coli among Japanese dairy farms.     

Development and assessment of sensitive immuno‐PCR assays for the quantification of cerebrospinal fluid three‐ and four‐repeat tau isoforms in tauopathies

Characteristic tau isoform composition of the insoluble fibrillar tau inclusions define tauopathies, including Alzheimer's disease (AD), progressive supranuclear palsy (PSP) and frontotemporal dementia with parkinsonism linked to chromosome 17/frontotemporal lobar degeneration‐tau (FTDP‐17/FTLD‐tau). Exon 10 splicing mutations in the tau gene, MAPT, in familial FTDP‐17 cause elevation of tau isoforms with four microtubule‐binding repeat domains (4R‐tau) compared to those with three repeats (3R‐tau). On the basis of two well‐characterised monoclonal antibodies against 3R‐ and 4R‐tau, we developed novel, sensitive immuno‐PCR assays for measuring the trace amounts of these isoforms in CSF. This was with the aim of assessing if CSF tau isoform changes reflect the pathological changes in tau isoform homeostasis in the degenerative brain and if these would be relevant for differential clinical diagnosis. Initial analysis of clinical CSF samples of PSP (n = 46), corticobasal syndrome (CBS; n = 22), AD (n = 11), Parkinson's disease with dementia (PDD; n = 16) and 35 controls revealed selective decreases of immunoreactive 4R‐tau in CSF of PSP and AD patients compared with controls, and lower 4R‐tau levels in AD compared with PDD. These decreases could be related to the disease‐specific conformational masking of the RD4‐binding epitope because of abnormal folding and/or aggregation of the 4R‐tau isoforms in tauopathies or increased sequestration of the 4R‐tau isoforms in brain tau pathology.     

A SNP test to identify Africanized honeybees via proportion of ‘African’ ancestry

The honeybee, Apis mellifera, is the world's most important pollinator and is ubiquitous in most agricultural ecosystems. Four major evolutionary lineages and at least 24 subspecies are recognized. Commercial populations are mainly derived from subspecies originating in Europe (75–95%). The Africanized honeybee is a New World hybrid of A. m. scutellata from Africa and European subspecies, with the African component making up 50–90% of the genome. Africanized honeybees are considered undesirable for bee‐keeping in most countries, due to their extreme defensiveness and poor honey production. The international trade in honeybees is restricted, due in part to bans on the importation of queens (and semen) from countries where Africanized honeybees are extant. Some desirable strains from the United States of America that have been bred for traits such as resistance to the mite Varroa destructor are unfortunately excluded from export to countries such as Australia due to the presence of Africanized honeybees in the USA. This study shows that a panel of 95 single nucleotide polymorphisms, chosen to differentiate between the African, Eastern European and Western European lineages, can detect Africanized honeybees with a high degree of confidence via ancestry assignment. Our panel therefore offers a valuable tool to mitigate the risks of spreading Africanized honeybees across the globe and may enable the resumption of queen and bee semen imports from the Americas.     

Locus‐specific ancestry to detect recent response to selection in admixed Swiss Fleckvieh cattle

Identification of selection signatures is one of the current endeavors of evolutionary genetics. Admixed populations may be used to infer post‐admixture selection. We calculated local ancestry for Swiss Fleckvieh, a composite of Simmental (SI) and Red Holstein Friesian (RHF), to infer such signals. Illumina Bovine SNP50 BeadChip data for 300 admixed, 88 SI and 97 RHF bulls were used. The average RHF ancestry across the whole genome was 0.70. To identify regions with high deviation from average, we considered two significance thresholds, based on a permutation test and extreme deviation from normal distribution. Regions on chromosomes 13 (46.3–47.3 Mb) and 18 (18.7–25.9 Mb) passed both thresholds in the direction of increased SI. Extended haplotype homozygosity within (iHS) and between (Rsb) populations was calculated to explore additional patterns of pre‐ and post‐admixture selection signals. The Rsb score of admixed and SI was significant in a wide region of chromosome 18 (6.6–24.6 Mb) overlapped with one area of strong local ancestry deviation. FTO, with pleiotropic effect on milk and fertility, NOD2 on dairy and NKD1 and SALL1 on fertility traits are located there. Genetic differentiation of RHF and SI (F_st), an alternative indicator of pre‐admixture selection in pure populations, was calculated. No considerable overlap of peaks of local ancestry deviations and F_st was observed. We found two regions with significant signatures of post‐admixture selection in this very young composite, applying comparatively stringent significance thresholds. The signals cover relatively large genomic areas and did not allow pinpointing of the gene(s) responsible for the apparent shift in ancestry proportions.     

Identification of contemporary selection signatures using composite log likelihood and their associations with marbling score in Korean cattle

Positive selection not only increases beneficial allele frequency but also causes augmentation of allele frequencies of sequence variants in close proximity. Signals for positive selection were detected by the statistical differences in subsequent allele frequencies. To identify selection signatures in Korean cattle, we applied a composite log‐likelihood (CLL)‐based method, which calculates a composite likelihood of the allelic frequencies observed across sliding windows of five adjacent loci and compares the value with the critical statistic estimated by 50 000 permutations. Data for a total of 11 799 nucleotide polymorphisms were used with 71 Korean cattle and 209 foreign beef cattle. As a result, 147 signals were identified for Korean cattle based on CLL estimates (P < 0.01). The signals might be candidate genetic factors for meat quality by which the Korean cattle have been selected. Further genetic association analysis with 41 intragenic variants in the selection signatures with the greatest CLL for each chromosome revealed that marbling score was associated with five variants. Intensive association studies with all the selection signatures identified in this study are required to exclude signals associated with other phenotypes or signals falsely detected and thus to identify genetic markers for meat quality.     

Targeted high‐throughput growth hormone 1 gene sequencing reveals high within‐breed genetic diversity in South African goats

This study assessed the genetic diversity in the growth hormone 1 gene (GH1) within and between South African goat breeds. Polymerase chain reaction‐targeted gene amplification together with Illumina MiSeq next‐generation sequencing (NGS) was used to generate the full length (2.54 kb) of the growth hormone 1 gene and screen for SNPs in the South African Boer (SAB) (n = 17), Tankwa (n = 15) and South African village (n = 35) goat populations. A range of 27–58 SNPs per population were observed. Mutations resulting in amino acid changes were observed at exons 2 and 5. Higher within‐breed diversity of 97.37% was observed within the population category consisting of SA village ecotypes and the Tankwa goats. Highest pairwise F_ST values ranging from 0.148 to 0.356 were observed between the SAB and both the South African village and Tankwa feral goat populations. Phylogenetic analysis indicated nine genetic clusters, which reflected close relationships between the South African populations and the other international breeds with the exception of the Italian Sarda breeds. Results imply greater potential for within‐population selection programs, particularly with SA village goats.     

Seasonal dynamics and co‐occurrence patterns of honey bee pathogens revealed by high‐throughput RT‐qPCR analysis

The health of the honey bee Apis mellifera is challenged by introduced parasites that interact with its inherent pathogens and cause elevated rates of colony losses. To elucidate co‐occurrence, population dynamics, and synergistic interactions of honey bee pathogens, we established an array of diagnostic assays for a high‐throughput qPCR platform. Assuming that interaction of pathogens requires co‐occurrence within the same individual, single worker bees were analyzed instead of collective samples. Eleven viruses, four parasites, and three pathogenic bacteria were quantified in more than one thousand single bees sampled from sixteen disease‐free apiaries in Southwest Germany. The most abundant viruses were black queen cell virus (84%), Lake Sinai virus 1 (42%), and deformed wing virus B (35%). Forager bees from asymptomatic colonies were infected with two different viruses in average, and simultaneous infection with four to six viruses was common (14%). Also, the intestinal parasites Nosema ceranae (96%) and Crithidia mellificae/Lotmaria passim (52%) occurred very frequently. These results indicate that low‐level infections in honey bees are more common than previously assumed. All viruses showed seasonal variation, while N. ceranae did not. The foulbrood bacteria Paenibacillus larvae and Melissococcus plutonius were regionally distributed. Spearman's correlations and multiple regression analysis indicated possible synergistic interactions between the common pathogens, particularly for black queen cell virus. Beyond its suitability for further studies on honeybees, this targeted approach may be, due to its precision, capacity, and flexibility, a viable alternative to more expensive, sequencing‐based approaches in nonmodel systems.     

Genome‐wide analyses suggest parallel selection for universal traits may eclipse local environmental selection in a highly mobile carnivore

Ecological and environmental heterogeneity can produce genetic differentiation in highly mobile species. Accordingly, local adaptation may be expected across comparatively short distances in the presence of marked environmental gradients. Within the European continent, wolves (Canis lupus) exhibit distinct north–south population differentiation. We investigated more than 67‐K single nucleotide polymorphism (SNP) loci for signatures of local adaptation in 59 unrelated wolves from four previously identified population clusters (northcentral Europe n = 32, Carpathian Mountains n = 7, Dinaric‐Balkan n = 9, Ukrainian Steppe n = 11). Our analyses combined identification of outlier loci with findings from genome‐wide association study of individual genomic profiles and 12 environmental variables. We identified 353 candidate SNP loci. We examined the SNP position and neighboring megabase (1 Mb, one million bases) regions in the dog (C. lupus familiaris) genome for genes potentially under selection, including homologue genes in other vertebrates. These regions included functional genes for, for example, temperature regulation that may indicate local adaptation and genes controlling for functions universally important for wolves, including olfaction, hearing, vision, and cognitive functions. We also observed strong outliers not associated with any of the investigated variables, which could suggest selective pressures associated with other unmeasured environmental variables and/or demographic factors. These patterns are further supported by the examination of spatial distributions of the SNPs associated with universally important traits, which typically show marked differences in allele frequencies among population clusters. Accordingly, parallel selection for features important to all wolves may eclipse local environmental selection and implies long‐term separation among population clusters.     

Evolution and selection of Rhg1, a copy‐number variant nematode‐resistance locus

The soybean cyst nematode (SCN) resistance locus Rhg1 is a tandem repeat of a 31.2 kb unit of the soybean genome. Each 31.2‐kb unit contains four genes. One allele of Rhg1, Rhg1‐b, is responsible for protecting most US soybean production from SCN. Whole‐genome sequencing was performed, and PCR assays were developed to investigate allelic variation in sequence and copy number of the Rhg1 locus across a population of soybean germplasm accessions. Four distinct sequences of the 31.2‐kb repeat unit were identified, and some Rhg1 alleles carry up to three different types of repeat unit. The total number of copies of the repeat varies from 1 to 10 per haploid genome. Both copy number and sequence of the repeat correlate with the resistance phenotype, and the Rhg1 locus shows strong signatures of selection. Significant linkage disequilibrium in the genome outside the boundaries of the repeat allowed the Rhg1 genotype to be inferred using high‐density single nucleotide polymorphism genotyping of 15 996 accessions. Over 860 germplasm accessions were found likely to possess Rhg1 alleles. The regions surrounding the repeat show indications of non‐neutral evolution and high genetic variability in populations from different geographic locations, but without evidence of fixation of the resistant genotype. A compelling explanation of these results is that balancing selection is in operation at Rhg1.     

SNPs selected by information content outperform randomly selected microsatellite loci for delineating genetic identification and introgression in the endangered dark European honeybee (Apis mellifera mellifera)

The honeybee (Apis mellifera) has been threatened by multiple factors including pests and pathogens, pesticides and loss of locally adapted gene complexes due to replacement and introgression. In western Europe, the genetic integrity of the native A. m. mellifera (M‐lineage) is endangered due to trading and intensive queen breeding with commercial subspecies of eastern European ancestry (C‐lineage). Effective conservation actions require reliable molecular tools to identify pure‐bred A. m. mellifera colonies. Microsatellites have been preferred for identification of A. m. mellifera stocks across conservation centres. However, owing to high throughput, easy transferability between laboratories and low genotyping error, SNPs promise to become popular. Here, we compared the resolving power of a widely utilized microsatellite set to detect structure and introgression with that of different sets that combine a variable number of SNPs selected for their information content and genomic proximity to the microsatellite loci. Contrary to every SNP data set, microsatellites did not discriminate between the two lineages in the PCA space. Mean introgression proportions were identical across the two marker types, although at the individual level, microsatellites' performance was relatively poor at the upper range of Q‐values, a result reflected by their lower precision. Our results suggest that SNPs are more accurate and powerful than microsatellites for identification of A. m. mellifera colonies, especially when they are selected by information content.     

Haplotype‐based genotyping‐by‐sequencing in oat genome research

In a de novo genotyping‐by‐sequencing (GBS) analysis of short, 64‐base tag‐level haplotypes in 4657 accessions of cultivated oat, we discovered 164741 tag‐level (TL) genetic variants containing 241224 SNPs. From this, the marker density of an oat consensus map was increased by the addition of more than 70000 loci. The mapped TL genotypes of a 635‐line diversity panel were used to infer chromosome‐level (CL) haplotype maps. These maps revealed differences in the number and size of haplotype blocks, as well as differences in haplotype diversity between chromosomes and subsets of the diversity panel. We then explored potential benefits of SNP vs. TL vs. CL GBS variants for mapping, high‐resolution genome analysis and genomic selection in oats. A combined genome‐wide association study (GWAS) of heading date from multiple locations using both TL haplotypes and individual SNP markers identified 184 significant associations. A comparative GWAS using TL haplotypes, CL haplotype blocks and their combinations demonstrated the superiority of using TL haplotype markers. Using a principal component‐based genome‐wide scan, genomic regions containing signatures of selection were identified. These regions may contain genes that are responsible for the local adaptation of oats to Northern American conditions. Genomic selection for heading date using TL haplotypes or SNP markers gave comparable and promising prediction accuracies of up to r = 0.74. Genomic selection carried out in an independent calibration and test population for heading date gave promising prediction accuracies that ranged between r = 0.42 and 0.67. In conclusion, TL haplotype GBS‐derived markers facilitate genome analysis and genomic selection in oat.     

Signatures of local adaptation along environmental gradients in a range‐expanding damselfly (Ischnura elegans)

Insect distributions are shifting rapidly in response to climate change and are undergoing rapid evolutionary change. We investigate the molecular signatures underlying local adaptation in the range‐expanding damselfly, Ischnura elegans. Using a landscape genomic approach combined with generalized dissimilarity modelling (GDM), we detect selection signatures on loci via allelic frequency change along environmental gradients. We analyse 13,612 single nucleotide polymorphisms (SNPs), derived from restriction site‐associated DNA sequencing (RADseq), in 426 individuals from 25 sites spanning the I. elegans distribution in Sweden, including its expanding northern range edge. Environmental association analysis (EAA) and the magnitude of allele frequency change along the range expansion gradient revealed significant signatures of selection in relation to high maximum summer temperature, high mean annual precipitation and low wind speeds at the range edge. SNP annotations with significant signatures of selection revealed gene functions associated with ongoing range expansion, including heat shock proteins (HSP40 and HSP70), ion transport (V‐ATPase) and visual processes (long‐wavelength‐sensitive opsin), which have implications for thermal stress response, salinity tolerance and mate discrimination, respectively. We also identified environmental thresholds where climate‐mediated selection is likely to be strong, and indicate that I. elegans is rapidly adapting to the climatic environment during its ongoing range expansion. Our findings empirically validate an integrative approach for detecting spatially explicit signatures of local adaptation along environmental gradients.     

Reproduction of parasitic mites Varroa destructor in original and new honeybee hosts

The ectoparasitic mite, Varroa destructor, shifted host from the eastern honeybee, Apis cerana, to the western honeybee, Apis mellifera. Whereas the original host survives infestations by this parasite, they are lethal to colonies of its new host. Here, we investigated a population of A. cerana naturally infested by the V. destructor Korea haplotype that gave rise to the globally invasive mite lineage. Our aim was to better characterize traits that allow for the survival of the original host to infestations by this particular mite haplotype. A known major trait of resistance is the lack of mite reproduction on worker brood in A. cerana. We show that this trait is neither due to a lack of host attractiveness nor of reproduction initiation by the parasite. However, successful mite reproduction was prevented by abnormal host development. Adult A. cerana workers recognized this state and removed hosts and parasites, which greatly affected the fitness of the parasite. These results confirm and complete previous observations of brood susceptibility to infestation in other honeybee host populations, provide new insights into the coevolution between hosts and parasites in this system, and may contribute to mitigating the large‐scale colony losses of A. mellifera due to V. destructor.     

Geographically widespread honeybee‐gut symbiont subgroups show locally distinct antibiotic‐resistant patterns

How long‐term antibiotic treatment affects host bacterial associations is still largely unknown. The honeybee‐gut microbiota has a simple composition, so we used this gut community to investigate how long‐term antibiotic treatment affects host‐associated microbiota. We investigated the phylogenetic relatedness, genomic content (GC percentage, genome size, number of genes and CRISPR) and antibiotic‐resistant genes (ARG) for strains from two abundant members of the honeybee core gut microbiota (Gilliamella apicola and Snodgrassella alvi). Domesticated honeybees are subjected to geographically different management policies, so we used two research apiaries, representing different antibiotic treatment regimens in their apiculture: low antibiotic usage (Norway) and high antibiotic usage (Arizona, USA). We applied whole‐genome shotgun sequencing on 48 G. apicola and 22 S. alvi. We identified three predominating subgroups of G. apicola in honeybees from both Norway and Arizona. For G. apicola, genetic content substantially varied between subgroups and distance similarity calculations showed similarity discrepancy between subgroups. Functional differences between subgroups, such as pectin‐degrading enzymes (G. apicola), were also identified. In addition, we identified horizontal gene transfer (HGT) of transposon (Tn10)‐associated tetracycline resistance (Tet B) across the G. apicola subgroups in the Arizonan honeybees, using interspace polymorphisms in the Tet B determinant. Our results support that honeybee‐gut symbiont subgroups can resist long‐term antibiotic treatment and maintain functionality through acquisition of geographically distinct antibiotic‐resistant genes by HGT.