Common signatures of selection reveal target loci for breeding across soybean populations.

Understanding the underlying genetic bases of yield-related selection and distinguishing these changes from genetic drift are critical for both improved understanding and future success of plant breeding. Soybean [Glycine max (L.) Merr.] is a key species for world food security, yet knowledge of the mechanism of selective breeding in soybean, such as the century-long program of artificial selection in U.S. soybean germplasm, is currently limited to certain genes and loci. Here, we identify genome-wide signatures of selection in separate populations of soybean subjected to artificial selection for increased yield by multiple breeding programs in the United States. We compared the alternative soybean breeding population (AGP) created by USDA-ARS to the conventional public soybean lines (CGP) developed at three different stages of breeding (ancestral, intermediate, and elite) to identify shared signatures of selection and differentiate these from drift. The results showed a strong selection for specific haplotypes identified by single site frequency and haplotype homozygosity methods. A set of common selection signatures was identified in both AGP and CGP that supports the hypothesis that separate breeding programs within similar environments coalesce on the fixation of the same key haplotypes. Signatures unique to each breeding program were observed. These results raise the possibility that selection analysis can allow the identification of favorable alleles to enhance directed breeding approaches.

Trisiloxane Surfactants Negatively Affect Reproductive Behaviors and Enhance Viral Replication in Honey Bees.

Trisiloxane surfactants are often applied in formulated adjuvant products to blooming crops, including almonds, exposing the managed honey bees (Apis mellifera) used for pollination of these crops and persisting in colony matrices, such as bee bread. Despite this, little is known regarding the effects of trisiloxane surfactants on important aspects of colony health, such as reproduction. In the present study, we use laboratory assays to examine how exposure to field-relevant concentrations of three trisiloxane surfactants found in commonly used adjuvant formulations affect queen oviposition rates, worker interactions with the queen, and worker susceptibility to endogenous viral pathogens. Trisiloxane surfactants were administered at 5 mg/kg in pollen supplement diet for 14 days. No effects on worker behavior or physiology could be detected, but our results demonstrate that hydroxy-capped trisiloxane surfactants can negatively affect queen oviposition and methyl-capped trisiloxane surfactants cause increased replication of Deformed Wing Virus in workers, suggesting that trisiloxane surfactant use while honey bees are foraging may negatively impact colony longevity and growth. Environ Toxicol Chem 2024;43:222-233. © 2023 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Single-cell dissection of aggression in honeybee colonies.

Understanding how genotypic variation results in phenotypic variation is especially difficult for collective behaviour because group phenotypes arise from complex interactions among group members. A genome-wide association study identified hundreds of genes associated with colony-level variation in honeybee aggression, many of which also showed strong signals of positive selection, but the influence of these 'colony aggression genes' on brain function was unknown. Here we use single-cell (sc) transcriptomics and gene regulatory network (GRN) analyses to test the hypothesis that genetic variation for colony aggression influences individual differences in brain gene expression and/or gene regulation. We compared soldiers, which respond to territorial intrusion with stinging attacks, and foragers, which do not. Colony environment showed stronger influences on soldier-forager differences in brain gene regulation compared with brain gene expression. GRN plasticity was strongly associated with colony aggression, with larger differences in GRN dynamics detected between soldiers and foragers from more aggressive relative to less aggressive colonies. The regulatory dynamics of subnetworks composed of genes associated with colony aggression genes were more strongly correlated with each other across different cell types and brain regions relative to other genes, especially in brain regions involved with olfaction and vision and multimodal sensory integration, which are known to mediate bee aggression. These results show how group genetics can shape a collective phenotype by modulating individual brain gene regulatory network architecture.

Impact of multiple selective breeding programs on genetic diversity in soybean germplasm.

KEY MESSAGE: Independent soybean breeding programs shape genetic diversity from unimproved germplasm to modern cultivars in similar ways, but distinct breeding populations retain unique genetic variation, preserving additional diversity. From the domestication of wild soybean (Glycine soja Sieb. & Zucc.), over 3,000 years ago, to the modern soybean (Glycine max L. Merr) cultivars that provide much of the world's oil and protein, soybean populations have undergone fundamental changes. We evaluated the molecular impact of breeding and selection using 391 soybean accessions including US cultivars and their progenitors from the USDA Soybean Germplasm Collection (CGP), plus two new populations specifically developed to increase genetic diversity and high yield in two alternative gene pools: one derived from exotic G. max germplasm (AGP) and one derived from G. soja (SGP). Reduction in nucleotide genetic diversity (π) was observed with selection within gene pools, but artificial selection in the AGP maintained more diversity than in the CGP. The highest FST levels were seen between ancestral and elite lines in all gene pools, but specific nucleotide-level patterns varied between gene pools. Population structure analyses support that independent selection resulted in high-yielding elite lines with similar allelic compositions in the AGP and CGP. SGP, however, produced elite progeny that were well differentiated from, but lower yielding than, CGP elites. Both the AGP and SGP retained a significant number of private alleles that are absent in CGP. We conclude that the genomic diversity shaped by multiple selective breeding programs can result in gene pools of highly productive elite lines with similar allelic compositions in a genome-wide perspective. Breeding programs with different ancestral lines, however, can retain private alleles representing unique genetic diversity.

Russian honey bee genotype identification through enhanced marker panel set.

Russian honey bees (RHB) are a breeding population developed by USDA-ARS as an effort to provide Varroa-resistant honey bees to beekeepers. The selection strategy for this breeding population was the first in honey bees to incorporate genetic stock identification (GSI). The original GSI approach has been in use for over a decade, and though effective, novel technologies and analytical approaches recently developed provide an opportunity for improvement. Here we outline a novel genotyping assay that capitalizes on the markers used in the GSI as well as new loci recently identified in a whole genome pooled study of commercial honey bee stocks. Our approach utilizes a microfluidic platform and machine learning analyses to arrive at an accurate, high throughput assay. This novel approach provides an improved tool that can be readily incorporated into breeding decisions towards healthier more productive bees.

Parallel mechanisms of visual memory formation across distinct regions of the honey bee brain.

Visual learning is vital to the behavioral ecology of the Western honey bee (Apis mellifera). Honey bee workers forage for floral resources, a behavior that requires the learning and long-term memory of visual landmarks, but how these memories are mapped to the brain remains poorly understood. To address this gap in our understanding, we collected bees that successfully learned visual associations in a conditioned aversion paradigm and compared gene expression correlates of memory formation in the mushroom bodies, a higher-order sensory integration center classically thought to contribute to learning, as well as the optic lobes, the primary visual neuropil responsible for sensory transduction of visual information. We quantified expression of CREB and CaMKII, two classical genetic markers of learning, and fen-1, a gene specifically associated with punishment learning in vertebrates. As expected, we found substantial involvement of the mushroom bodies for all three markers but additionally report the involvement of the optic lobes across a similar time course. Our findings imply the molecular involvement of a sensory neuropil during visual associative learning parallel to a higher-order brain region, furthering our understanding of how a tiny brain processes environmental signals.

Effects of aversive conditioning on expression of physiological stress in honey bees (Apis mellifera).

Stress is defined as any deviation from an organism's baseline physiological levels. Therefore, introduction of new stimuli and information, such as in learning, can be defined as a stressor. A large body of research exists examining the role that stress plays in learning, but virtually none addresses whether or not learning itself is a measurable cause of stress. The current study used a wide variety of learning centric stress responses. Researchers examined changes in expression of ten stress and learning related genes in various physiological systems in domesticated honey bees (Apis mellifera) as a result of exposure to an aversive conditioning task. Gene expression was examined using quantitative real-time polymerase chain reaction following the learning task. Results indicate that learning affects expression of some stress related genes.

Genome-wide patterns of differentiation within and among U.S. commercial honey bee stocks.

BACKGROUND: The population genetics of U.S. honey bee stocks remain poorly characterized despite the agricultural importance of Apis mellifera as the major crop pollinator. Commercial and research-based breeding programs have made significant improvements of favorable genetic traits (e.g. production and disease resistance). The variety of bees produced by artificial selection provides an opportunity to characterize the genetic diversity and regions of the genome undergoing selection in commonly managed stocks. RESULTS: Pooled sequencing of eight honey bee stocks found strong genetic similarity among six of the stocks. Two stocks, Pol-line and Hilo, showed significant differentiation likely due to their intense and largely closed breeding for resistance to the parasitic Varroa mite. Few variants were identified as being specific to any one stock, indicating potential admixture among the sequenced stocks. Juxtaposing the underlying genetic variation of stocks selected for disease- and parasite-resistance behavior, we identified genes and candidate regions putatively associated with resistance regulated by hygienic behavior. CONCLUSION: This study provides important insights into the distinct genetic characteristics and population diversity of honey bee stocks used in the United States, and provides further evidence of high levels of admixture in commercially managed honey bee stocks. Furthermore, breeding efforts to enhance parasite resistance in honey bees may have created unique genetic profiles. Genomic regions of interest have been highlighted for potential future work related to developing genetic markers for selection of disease and parasite resistance traits. Due to the vast genomic similarities found among stocks in general, our findings suggest that additional data regarding gene expression, epigenetic and regulatory information are needed to more fully determine how stock phenotypic diversity is regulated.

Genomic regions influencing aggressive behavior in honey bees are defined by colony allele frequencies.

For social animals, the genotypes of group members affect the social environment, and thus individual behavior, often indirectly. We used genome-wide association studies (GWAS) to determine the influence of individual vs. group genotypes on aggression in honey bees. Aggression in honey bees arises from the coordinated actions of colony members, primarily nonreproductive "soldier" bees, and thus, experiences evolutionary selection at the colony level. Here, we show that individual behavior is influenced by colony environment, which in turn, is shaped by allele frequency within colonies. Using a population with a range of aggression, we sequenced individual whole genomes and looked for genotype-behavior associations within colonies in a common environment. There were no significant correlations between individual aggression and specific alleles. By contrast, we found strong correlations between colony aggression and the frequencies of specific alleles within colonies, despite a small number of colonies. Associations at the colony level were highly significant and were very similar among both soldiers and foragers, but they covaried with one another. One strongly significant association peak, containing an ortholog of the Drosophila sensory gene dpr4 on linkage group (chromosome) 7, showed strong signals of both selection and admixture during the evolution of gentleness in a honey bee population. We thus found links between colony genetics and group behavior and also, molecular evidence for group-level selection, acting at the colony level. We conclude that group genetics dominates individual genetics in determining the fatal decision of honey bees to sting.

Developmental plasticity shapes social traits and selection in a facultatively eusocial bee.

Developmental plasticity generates phenotypic variation, but how it contributes to evolutionary change is unclear. Phenotypes of individuals in caste-based (eusocial) societies are particularly sensitive to developmental processes, and the evolutionary origins of eusociality may be rooted in developmental plasticity of ancestral forms. We used an integrative genomics approach to evaluate the relationships among developmental plasticity, molecular evolution, and social behavior in a bee species (Megalopta genalis) that expresses flexible sociality, and thus provides a window into the factors that may have been important at the evolutionary origins of eusociality. We find that differences in social behavior are derived from genes that also regulate sex differentiation and metamorphosis. Positive selection on social traits is influenced by the function of these genes in development. We further identify evidence that social polyphenisms may become encoded in the genome via genetic changes in regulatory regions, specifically in transcription factor binding sites. Taken together, our results provide evidence that developmental plasticity provides the substrate for evolutionary novelty and shapes the selective landscape for molecular evolution in a major evolutionary innovation: Eusociality.

Colonization history and population differentiation of the Honey Bees (Apis mellifera L.) in Puerto Rico.

Honey bees (Apis mellifera L.) are the primary commercial pollinators across the world. The subspecies A. m. scutellata originated in Africa and was introduced to the Americas in 1956. For the last 60 years, it hybridized successfully with European subspecies, previous residents in the area. The result of this hybridization was called Africanized honey bee (AHB). AHB has spread since then, arriving to Puerto Rico (PR) in 1994. The honey bee population on the island acquired a mosaic of features from AHB or the European honey bee (EHB). AHB in Puerto Rico shows a major distinctive characteristic, docile behavior, and is called gentle Africanized honey bees (gAHB). We used 917 SNPs to examine the population structure, genetic differentiation, origin, and history of range expansion and colonization of gAHB in PR. We compared gAHB to populations that span the current distribution of A. mellifera worldwide. The gAHB population is shown to be a single population that differs genetically from the examined populations of AHB. Texas and PR groups are the closest genetically. Our results support the hypothesis that the Texas AHB population is the source of gAHB in Puerto Rico.

A soft selective sweep during rapid evolution of gentle behaviour in an Africanized honeybee.

Highly aggressive Africanized honeybees (AHB) invaded Puerto Rico (PR) in 1994, displacing gentle European honeybees (EHB) in many locations. Gentle AHB (gAHB), unknown anywhere else in the world, subsequently evolved on the island within a few generations. Here we sequence whole genomes from gAHB and EHB populations, as well as a North American AHB population, a likely source of the founder AHB on PR. We show that gAHB retains high levels of genetic diversity after evolution of gentle behaviour, despite selection on standing variation. We observe multiple genomic loci with significant signatures of selection. Rapid evolution during colonization of novel habitats can generate major changes to characteristics such as morphological or colouration traits, usually controlled by one or more major genetic loci. Here we describe a soft selective sweep, acting at multiple loci across the genome, that occurred during, and may have mediated, the rapid evolution of a behavioural trait.

Social signals and aversive learning in honey bee drones and workers.

The dissemination of information is a basic element of group cohesion. In honey bees (Apis mellifera Linnaeus 1758), like in other social insects, the principal method for colony-wide information exchange is communication via pheromones. This medium of communication allows multiple individuals to conduct tasks critical to colony survival. Social signaling also establishes conflict at the level of the individual who must trade-off between attending to the immediate environment or the social demand. In this study we examined this conflict by challenging highly social worker honey bees, and less social male drone honey bees undergoing aversive training by presenting them with a social stress signal (isopentyl acetate, IPA). We utilized IPA exposure methods that caused lower learning performance in appetitive learning in workers. Exposure to isopentyl acetate (IPA) did not affect performance of drones and had a dose-specific effect on worker response, with positive effects diminishing at higher IPA doses. The IPA effects are specific because non-social cues, such as the odor cineole, improve learning performance in drones, and social homing signals (geraniol) did not have a discernible effect on drone or worker performance. We conclude that social signals do generate conflict and that response to them is dependent on signal relevance to the individual as well as the context. We discuss the effect of social signal on learning both related to its social role and potential evolutionary history.

Measuring individual locomotor rhythms in honey bees, paper wasps and other similar-sized insects.

Circadian rhythms in social insects are highly plastic and are modulated by multiple factors. In addition, complex behaviors such as sun-compass orientation and time learning are clearly regulated by the circadian system in these organisms. Despite these unique features of social insect clocks, the mechanisms as well as the functional and evolutionary relevance of these traits remain largely unknown. Here we show a modification of the Drosophila activity monitoring (DAM) system that allowed us to measure locomotor rhythms of the honey bee, Apis mellifera (three variants; gAHB, carnica and caucasica), and two paper wasps (Polistes crinitus and Mischocyttarus phthisicus). A side-by-side comparison of the endogenous period under constant darkness (free-running period) led us to the realization that these social insects exhibit significant deviations from the Earth's 24 h rotational period as well as a large degree of inter-individual variation compared with Drosophila. Experiments at different temperatures, using honey bees as a model, revealed that testing the endogenous rhythm at 35°C, which is the hive's core temperature, results in average periods closer to 24 h compared with 25°C (23.8 h at 35°C versus 22.7 h at 25°C). This finding suggests that the degree of tuning of circadian temperature compensation varies among different organisms. We expect that the commercial availability, cost-effectiveness and integrated nature of this monitoring system will facilitate the growth of the circadian field in these social insects and catalyze our understanding of the mechanisms as well as the functional and evolutionary relevance of circadian rhythms.

Individual responsiveness to shock and colony-level aggression in honey bees: evidence for a genetic component.

The phenotype of the social group is related to phenotypes of individuals that form that society. We examined how honey bee colony aggressiveness relates to individual response of male drones and foraging workers. Although the natural focus in colony aggression has been on the worker caste, the sterile females engaged in colony maintenance and defense, males carry the same genes. We measured aggressiveness scores of colonies and examined components of individual aggressive behavior in workers and haploid sons of workers from the same colony. We describe for the first time, that males, although they have no stinger, do bend their abdomen (abdominal flexion) in a posture similar to stinging behavior of workers in response to electric shock. Individual worker sting response and movement rates in response to shock were significantly correlated with colony scores. In the case of drones, sons of workers from the same colonies, abdominal flexion significantly correlated but their movement rates did not correlate with colony aggressiveness. Furthermore, the number of workers responding at increasing levels of voltage exhibits a threshold-like response, whereas the drones respond in increasing proportion to shock. We conclude that there are common and caste-specific components to aggressive behavior in honey bees. We discuss implications of these results on social and behavioral regulation and genetics of aggressive response.

Aversive conditioning in honey bees (Apis mellifera anatolica): a comparison of drones and workers.

Honey bees provide a model system to elucidate the relationship between sociality and complex behaviors within the same species, as females (workers) are highly social and males (drones) are more solitary. We report on aversive learning studies in drone and worker honey bees (Apis mellifera anatolica) in escape, punishment and discriminative punishment situations. In all three experiments, a newly developed electric shock avoidance assay was used. The comparisons of expected and observed responses were performed with conventional statistical methods and a systematic randomization modeling approach called object oriented modeling. The escape experiment consisted of two measurements recorded in a master-yoked paradigm: frequency of response and latency to respond following administration of shock. Master individuals could terminate an unavoidable shock triggered by a decrementing 30 s timer by crossing the shuttlebox centerline following shock activation. Across all groups, there was large individual response variation. When assessing group response frequency and latency, master subjects performed better than yoked subjects for both workers and drones. In the punishment experiment, individuals were shocked upon entering the shock portion of a bilaterally wired shuttlebox. The shock portion was spatially static and unsignalled. Only workers effectively avoided the shock. The discriminative punishment experiment repeated the punishment experiment but included a counterbalanced blue and yellow background signal and the side of shock was manipulated. Drones correctly responded less than workers when shock was paired with blue. However, when shock was paired with yellow there was no observable difference between drones and workers.

Mutational bias for body size in rhabditid nematodes.

Mutational bias is a potentially important agent of evolution, but it is difficult to disentangle the effects of mutation from those of natural selection. Mutation-accumulation experiments, in which mutations are allowed to accumulate at very small population size, thus minimizing the efficiency of natural selection, are the best way to separate the effects of mutation from those of selection. Body size varies greatly among species of nematode in the family rhabditidae; mutational biases are both a potential cause and a consequence of that variation. We report data on the cumulative effects of mutations that affect body size in three species of rhabditid nematode that vary fivefold in adult size. Results are very consistent with previous studies of mutations underlying fitness in the same strains: two strains of Caenorhabditis briggsae decline in body size about twice as fast as two strains of C. elegans, with a concomitant higher point estimate of the genomic mutation rate; the confamilial Oscheius myriophila is intermediate. There is an overall mutational bias, such that mutations reduce size on average, but the bias appears consistent between species. The genetic correlation between mutations that affect size and those underlying fitness is large and positive, on average.

Cumulative effects of spontaneous mutations for fitness in Caenorhabditis: role of genotype, environment and stress.

It is often assumed that the mutation rate is an evolutionarily optimized property of a taxon. The relevant mutation rate is for mutations that affect fitness, U, but the strength of selection on the mutation rate depends on the average effect of a mutation. Determination of U is complicated by the possibility that mutational effects depend on the particular environmental context in which the organism exists. It has been suggested that the effects of deleterious mutations are typically magnified in stressful environments, but most studies confound genotype with environment, so it is unclear to what extent environmental specificity of mutations is specific to a particular starting genotype. We report a study designed to separate effects of species, genotype, and environment on the degradation of fitness resulting from new mutations. Mutations accumulated for >200 generations at 20 degrees in two strains of two species of nematodes that differ in thermal sensitivity. Caenorhabditis briggsae and C. elegans have similar demography at 20 degrees, but C. elegans suffers markedly reduced fitness at 25 degrees. We find little evidence that mutational properties differ depending on environmental conditions and mutational correlations between environments are close to those expected if effects were identical in both environments.