Clinal variation and selection of MDH allozymes in honey bee populations

Clinal variation in natural populations of Drosophila melanogaster: An old debate about natural selection and neutral processes

Distinguishing between environmental adaptations and neutral processes poses a challenge in population genetics and evolutionary studies, particularly when phenomena can be explained by both processes. Clines are genotypic or phenotypic characters correlated with environmental variables, because of that correlation, they are used as examples of spatially varying selection. At the same time, many genotypic clines can be explained by demographic history, like isolation by distance or secondary contact zones. Clines have been extensively studied in Drosophila melanogaster, especially in North America and Australia, where they are attributed to both differential selection and various demographic processes. This review explores existing literature supporting this conclusion and suggests new approaches to better understand the influence of these processes on clines. These innovative approaches aim to shed light on the longstanding debate regarding the importance of natural selection versus neutral processes in maintaining variation in natural populations.

Addressing emerging issues in entomology: 2023 student debates

The Entomological Society of America (ESA) Student Debates is an annual student competition at the ESA Annual Meeting organized by Student Debates Subcommittee (SDS) members of the ESA Student Affairs Committee. In conjunction with the 2023 ESA Annual Meeting theme, 'Insects and influence: Advancing entomology's impact on people and policy', the theme of this year's student debate was 'Addressing emerging issues in entomology'. With the aid of ESA membership, the SDS selected the following debate topics: (1) Should disclosure of artificial intelligence large language models in scientific writing always be required? and (2) Is it more important to prioritize honey bee or native pollinator health for long-term food security within North America? Four student teams from across the nation, composed of 3-5 student members and a professional advisor, were assigned a topic and stance. Over the course of 5 months, all team members researched and prepared for their assigned topic before debating live with an opposing team at the 2023 ESA Annual Meeting in National Harbor, Maryland. SDS members additionally prepared and presented introductions for each debate topic to provide unbiased backgrounds to the judges and audience for context in assessing teams' arguments. The result was an engaging discussion between our teams, judges, and audience members on emerging issues facing entomology and its impact on people and policy, such as scientific communication and food security, that brought attention to the complexities involved when debating topics concerning insects and influence.

Metabolic rate shapes differences in foraging efficiency among honeybee foragers

Metabolic rate is the rate at which organisms process energy and is often considered as the fundamental driver of life history processes. The link between metabolic rate and life history is critically mediated via foraging, which shapes the energy acquisition patterns of an individual. This predicts that individuals with different metabolic rates likely vary in their foraging strategies, although such a link has rarely been empirically investigated in the context of optimal foraging theory—a powerful framework for understanding how animals maximize their foraging returns. Many central place foragers such as honeybees maximize their energetic efficiency rather than the rate of energetic gain, given the critical role of energetic costs in foraging decisions. We therefore tested if individuals with low and high metabolic rates differ in efficiency maximization, using genetic lines of honeybees with different metabolic rates. Our results show that low metabolic rate foragers visit more flowers during a foraging trip and have a higher energetic efficiency than high metabolic rate foragers in both low and high resource conditions. We discuss the significance of these results in the context of division of labor and the adaptive role of phenotypic diversity in metabolic rate in a social insect colony.

Mitochondrial DNA Variation in Peruvian Honey Bee (Apis mellifera L.) Populations Using the tRNAleu-cox2 Intergenic Region

Simple Summary

Currently, the genetic diversity of Peruvian honey bee populations is unknown. Only two studies were carried out and suggest that many regions of Peru present Africanized honey bee colonies and a varied degree of Africanization. To molecularly characterize and know more about the genetic background of Peruvian honey bees, the highly polymorphic tRNAleu-cox2 was used. This study analyzed 512 colonies in three regions of Peru: Lima, Piura, and Junín. The results indicated that 65% colonies correspond to lineage A (African), 33.8% colonies to lineage C (Eastern European), and 1.2% colonies to lineage M (Western European). A total of 24 haplotypes were identified: 16 haplotypes belong to lineage A (sub-lineage A_I (13), sub-lineage A_III (03)), lineage C (06), and lineage M (02), and 15 of them are for the first time reported and represented by A1t, A1u, A1w, A4p, A4q, A4s, A4t, A4u, A4v, A4w, 30d, A30e, A65, M7b, and M7c. Piura and Lima presented higher proportions of African haplotypes and lower proportions of haplotypes from lineage C than Lima. Very few haplotypes of lineage M were identified, whose presence could be due to accidental purchases or traces of honey bee introductions from lineage M in the 19th century. Hence, studies about the diversity and genetic structure of Peruvian honey bee populations are necessary to promote adequate, sustainable management and establish conservation and breeding programs.

Abstract

Mitochondrial DNA variations of Peruvian honey bee populations were surveyed by using the tRNA^leu-cox2 intergenic region. Only two studies have characterized these populations, indicating the presence of Africanized honey bee colonies in different regions of Peru and varied levels of Africanization, but the current status of its genetic diversity is unknown. A total of 512 honey bee colonies were sampled from three regions to characterize them. Our results revealed the presence of European and African haplotypes: the African haplotypes identified belong to sub-lineage A_I (13) and sub-lineage A_III (03), and the European haplotypes to lineages C (06) and M (02). Of 24 haplotypes identified, 15 new sequences are reported here (11 sub-lineage A_I, 2 sub-lineage A_III, and 2 lineage M). Peruvian honey bee populations presented a higher proportion from African than European haplotypes. High proportions of African haplotype were reported for Piura and Junín, unlike Lima, which showed more European haplotypes from lineage C. Few colonies belonging to lineage M would represent accidental purchase or traces of the introduction into Peru in the 19th century.

Thermal adaptation in the honeybee (Apis mellifera) via changes to the structure of malate dehydrogenase

In honeybees there are three alleles of cytosolic malate dehydrogenase gene: F, M and S. Allele frequencies are correlated with environmental temperature, suggesting that the alleles have temperature-dependent fitness benefits. We determined the enzyme activity of each allele across a range of temperatures in vitro The F and S alleles have higher activity and are less sensitive to high temperatures than the M allele, which loses activity after incubation at temperatures found in the thorax of foraging bees in hot climates. Next, we predicted the protein structure of each allele and used molecular dynamics simulations to investigate their molecular flexibility. The M allozyme is more flexible than the S and F allozymes at 50°C, suggesting a plausible explanation for its loss of activity at high temperatures, and has the greatest structural flexibility at 15°C, suggesting that it can retain some enzyme activity at cooler temperatures. MM bees recovered from 2 h of cold narcosis significantly better than all other genotypes. Combined, these results explain clinal variation in malate dehydrogenase allele frequencies in the honeybee at the molecular level.

Revisiting the Iberian honey bee (Apis mellifera iberiensis) contact zone: maternal and genome-wide nuclear variations provide support for secondary contact from historical refugia

Dissecting diversity patterns of organisms endemic to Iberia has been truly challenging for a variety of taxa, and the Iberian honey bee is no exception. Surveys of genetic variation in the Iberian honey bee are among the most extensive for any honey bee subspecies. From these, differential and complex patterns of diversity have emerged, which have yet to be fully resolved. Here, we used a genome-wide data set of 309 neutrally tested single nucleotide polymorphisms (SNPs), scattered across the 16 honey bee chromosomes, which were genotyped in 711 haploid males. These SNPs were analysed along with an intergenic locus of the mtDNA, to reveal historical patterns of population structure across the entire range of the Iberian honey bee. Overall, patterns of population structure inferred from nuclear loci by multiple clustering approaches and geographic cline analysis were consistent with two major clusters forming a well-defined cline that bisects Iberia along a northeastern-southwestern axis, a pattern that remarkably parallels that of the mtDNA. While a mechanism of primary intergradation or isolation by distance could explain the observed clinal variation, our results are more consistent with an alternative model of secondary contact between divergent populations previously isolated in glacial refugia, as proposed for a growing list of other Iberian taxa. Despite current intense honey bee management, human-mediated processes have seemingly played a minor role in shaping Iberian honey bee genetic structure. This study highlights the complexity of the Iberian honey bee patterns and reinforces the importance of Iberia as a reservoir of Apis mellifera diversity.

Signatures of selection in the Iberian honey bee (Apis mellifera iberiensis) revealed by a genome scan analysis of single nucleotide polymorphisms

Understanding the genetic mechanisms of adaptive population divergence is one of the most fundamental endeavours in evolutionary biology and is becoming increasingly important as it will allow predictions about how organisms will respond to global environmental crisis. This is particularly important for the honey bee, a species of unquestionable ecological and economical importance that has been exposed to increasing human-mediated selection pressures. Here, we conducted a single nucleotide polymorphism (SNP)-based genome scan in honey bees collected across an environmental gradient in Iberia and used four FST -based outlier tests to identify genomic regions exhibiting signatures of selection. Additionally, we analysed associations between genetic and environmental data for the identification of factors that might be correlated or act as selective pressures. With these approaches, 4.4% (17 of 383) of outlier loci were cross-validated by four FST -based methods, and 8.9% (34 of 383) were cross-validated by at least three methods. Of the 34 outliers, 15 were found to be strongly associated with one or more environmental variables. Further support for selection, provided by functional genomic information, was particularly compelling for SNP outliers mapped to different genes putatively involved in the same function such as vision, xenobiotic detoxification and innate immune response. This study enabled a more rigorous consideration of selection as the underlying cause of diversity patterns in Iberian honey bees, representing an important first step towards the identification of polymorphisms implicated in local adaptation and possibly in response to recent human-mediated environmental changes.

Ecological adaptation of diverse honey bee (Apis mellifera) populations

Background

Honey bees are complex eusocial insects that provide a critical contribution to human agricultural food production. Their natural migration has selected for traits that increase fitness within geographical areas, but in parallel their domestication has selected for traits that enhance productivity and survival under local conditions. Elucidating the biochemical mechanisms of these local adaptive processes is a key goal of evolutionary biology. Proteomics provides tools unique among the major ‘omics disciplines for identifying the mechanisms employed by an organism in adapting to environmental challenges.

Results

Through proteome profiling of adult honey bee midgut from geographically dispersed, domesticated populations combined with multiple parallel statistical treatments, the data presented here suggest some of the major cellular processes involved in adapting to different climates. These findings provide insight into the molecular underpinnings that may confer an advantage to honey bee populations. Significantly, the major energy-producing pathways of the mitochondria, the organelle most closely involved in heat production, were consistently higher in bees that had adapted to colder climates. In opposition, up-regulation of protein metabolism capacity, from biosynthesis to degradation, had been selected for in bees from warmer climates.

Conclusions

Overall, our results present a proteomic interpretation of expression polymorphisms between honey bee ecotypes and provide insight into molecular aspects of local adaptation or selection with consequences for honey bee management and breeding. The implications of our findings extend beyond apiculture as they underscore the need to consider the interdependence of animal populations and their agro-ecological context.

Allozyme markers to help define the South American origins of Microctonus hyperodae (Hymenoptera: Braconidae) established in New Zealand for biological control of Argentine stem weevil

The thelytokous parasitoid, Microctonus hyperodae Loan, was collected from eight South American locations and introduced to New Zealand in 1991 for biological control of Argentine stem weevil, Listronotus bonariensis (Kuschel) (Coleoptera: Curculionidae). Parasitoids from each population were released in equal numbers at each New Zealand site to give them the same opportunities to establish. Population markers have been sought to identify the South American geographic populations that have become most successful in New Zealand. These would assist in determining the importance of concepts such as climate matching and host-parasitoid coevolution to the establishment of natural enemies in new regions for biological control. Vertical polyacrylamide electrophoresis was used to survey 16 enzymes and ten calcium binding proteins, and this paper reports variation at three putative loci. Malate dehydrogenase, a dihydrolipoamide dehydrogenase isozyme and a calcium binding protein exhibited clear genetic variation, each with two alleles. All M. hyperodae isofemale lines from east of the Andes mountains shared one genotype, all but one from west of the Andes shared another, while a population from within the Andes contained both genotypes. This variation was highly congruent with previously described morphometric variation. At two loci, the maintenance of heterozygotes, and the absence of homozygotes, within isofemale lines suggested M. hyperodae thelytoky is apomictic.

Clinal variation and selection on MDH allozymes in honeybees in Chile

Evidence of clinal variation and selection on Mdh-1 locus was observed in 27 samples from 22 sites in a 2800 km north-south transect across Chile. We found a negative correlation among F allele normalized frequency and mean temperature, and minimum temperature of January and July, as well as a positive correlation among S allele normalized frequency and annual mean, and minimum January temperatures. Our results lend weight to the idea that Chilean honeybee populations of colder areas have higher F allele frequencies, supporting previous claims that Mdh-1 allozymes of Apis mellifera are subject to temperature-mediated selection.

Environmental and genetic influences on flight metabolic rate in the honey bee, Apis mellifera

Flying honey bees demonstrate highly variable metabolic rates. The lowest reported values (approximately 0.3 Wg(-1)) occur in tethered bees generating the minimum lift to support their body weight, free-flying 2-day old bees, winter bees, or bees flying at high air temperatures (45 degrees C). The highest values (approximately 0.8 Wg(-1)) occur in foragers that are heavily loaded or flying in low-density air. In different studies, flight metabolic rate has increased, decreased, or remained constant with air temperature. Current research collectively suggests that this variation occurs because flight metabolic rates decrease at thorax temperatures above or below 38 degrees C. At 30 degrees C, approximately 30% of colonial energy is spent during typical foraging, so variation in flight metabolic rate can strongly affect colony-level energy balance. Higher air temperatures tend to increase colonial net gain rates, efficiencies and honey storage rates due to lower metabolic rates during flight and in the hive. Variation in flight metabolism has a clear genetic basis. Different genetic strains of honey bees often differ in flight metabolic rate, and these differences in flight physiology can be correlated with foraging effort, suggesting a possible pathway for selection effects on flight metabolism.

An ancient transpecific polymorphism shows extreme divergence in a multitrait cline in an intertidal snail (Nucella lapillus (L.))

Clines in intraspecific genetic variation are frequently associated with an environmental transition. Here, divergence among nucleotide sequences of two nuclear loci, cytosolic and mitochondrial malate dehydrogenase (cMDH and mMDH, respectively), is described, in a multitrait cline over a distance of ca. 3 km where shell phenotype, allozyme, mitochondrial DNA haplotype, and centric fusion (Robertsonian translocations) frequencies covary with temperature and humidity and change abruptly in a continuous population of the dog-whelk (Nucella lapillus), a common intertidal snail of the north temperate Atlantic. Protein electrophoresis has already shown two alleles of mMDH varying from fixation of one allele to near fixation of the other, whereas cMDH appears to be monomorphic. The results of this study show a striking disparity in nucleotide sequence divergence among alleles at the two loci, with extreme molecular differentiation in one of them. Four alleles of cMDH were found to have nucleotide and amino acid sequence divergences of 0.4% and 0.3%, respectively. In contrast, the two mMDH cDNA alleles differed by 23% and 20% at the nucleotide and amino acid levels, respectively. Analysis of a 91-bp partial nucleotide sequence of mMDH from Nucella freycineti, the closest relative of N. lapillus, revealed two similar alleles and indicated that the divergence in mMDH in N. lapillus represents an ancient transpecific polymorphism in these Nucella. Together with earlier studies on variation in N. lapillus, it is argued that the polymorphism in mMDH and the clines in N. lapillus represent the presence of two persistent coadapted gene complexes, multitrait coevolving genetic solutions to environmental variation, which may presently enable this snail to exploit a diverse environment successfully.

Genetic and morphological variation in three snow pool Aedes mosquito species of the subgenus Ochlerotatus (Diptera: Culicidae)

Thirteen Aedes hexodontus Dyar populations from throughout the western United States were examined for genetic and morphological variation. Analysis of allozyme frequencies at 16 loci revealed a pattern of genetic variation that formed a north-south cline across Washington, Oregon, and California in the number of alleles per locus, the percent of polymorphic loci, and the frequency of one allele of aconitase-1. Comparison of the genetic profile of Ae. hexodontus populations to two other widely distributed members of the punctor subgroup, Aedes punctor (Kirby) and Aedes abserratus (Felt & Young), revealed only one diagnostic locus for all three species. Seven loci exhibited frequency differences among species but were not diagnostic. Morphological characters also exhibited little variation within and among the three species. The adult females differed only in the scaling of the probasisternum. This area was extensively scaled in 91% of Ae. hexodontus specimens but bare or only lightly scaled in Ae. puntor and Ae. abserratus. No other differences were observed in the adult females or the male genitalia in any of the three species. The larvae of Ae. abserratus could be separated by the single-branched seta 2-X. Six larval characters differed between Ae. hexodontus and Ae. punctor but the ranges of each character overlapped and none were diagnostic. These comparisons indicated that Ae. hexodontus is a single species, at least in the southern part of its range. Also, genetic and morphological comparison of the three species within the punctor subgroup attested to the close relationship hypothesized for these mosquitoes.

Flight respiration and energetics

We use a comparative approach to examine some of the physiological traits that make flight possible. Comparisons of related fliers and runners suggest that fliers generally have higher aerobic metabolic capacities than runners but that the difference is highly dependent on the taxa studied. The high metabolic rates of fliers relative to runners, especially in insects, are correlated with high locomotory muscle cycle frequencies and low efficiencies of conversion of metabolic power to mechanical power. We examine some factors that produce variation in flight respiration and energetics. Air temperature strongly affects the flight metabolic rate of some insects and birds. Flight speed interacts with flier mass, so that small fliers tend to exhibit a J-shaped power curve and larger fliers a U-shaped power curve. As body size increases, mass-specific aerobic flight metabolism decreases in most studies, but mass-specific power output is constant or increases, leading to an increase in efficiency with size. Intraspecific studies have revealed specific genetically based effects on flight metabolism and power output and multiple ecological correlates of flight capabilities.

Cloning and primary structure of putative cytosolic and mitochondrial malate dehydrogenase from the mollusc Nucella lapillus (L.)

The evolutionary history of the malate dehydrogenase (MDH) gene family [NAD-dependent MDH; EC 1.1.1.37 and NAD(P)-dependent MDH; EC 1.1.1.82] has received much attention. MDHs have also featured extensively as electrophoretic markers in population genetics and evolutionary ecology, and in many cases, intraspecific variation in MDH has been correlated with environmental variables. However, while the amino acid residues essential for MDH function are known, no studies have examined intraspecific nucleotide variation despite evidence indicating that natural selection may be operating on this locus. This study presents two sets of degenerate oligonucleotide PCR primers to facilitate the cloning of cytosolic MDH (cMDH) and mitochondrial MDH (mMDH) from a broad range of animals (cMDH) and animals and plants (mMDH). These primers were used to obtain putative cMDH and mMDH cDNAs from the mollusc Nucella lapillus. The N. lapillus cMDH cDNA was found to encode a putative cMDH protein of 334aa and 36kDa, while the mMDH cDNA encoded a putative mature mMDH protein of 315aa and 33kDa. The putative amino acid sequences of the two compartmentalised N. lapillus MDHs are presented and compared to other known MDH sequences.