The gene csd is the primary signal for sexual development in the honeybee and encodes an SR-type protein

Mechanisms of Pathogen and Pesticide Resistance in Honey Bees

Bees are the most important insect pollinators of the crops humans grow, and Apis mellifera, the Western honey bee, is the most commonly managed species for this purpose. In addition to providing agricultural services, the complex biology of honey bees has been the subject of scientific study since the 18th century, and the intricate behaviors of honey bees and ants, fellow hymenopterans, inspired much sociobiological inquest. Unfortunately, honey bees are constantly exposed to parasites, pathogens, and xenobiotics, all of which pose threats to their health. Despite our curiosity about and dependence on honey bees, defining the molecular mechanisms underlying their interactions with biotic and abiotic stressors has been challenging. The very aspects of their physiology and behavior that make them so important to agriculture also make them challenging to study, relative to canonical model organisms. However, because we rely on A. mellifera so much for pollination, we must continue our efforts to understand what ails them. Here, we review major advancements in our knowledge of honey bee physiology, focusing on immunity and detoxification, and highlight some challenges that remain.

LncRNA gene ANTSR coordinates complementary sex determination in the Argentine ant

Animals have evolved various sex determination systems. Here, we describe a newly found mechanism. A long noncoding RNA (lncRNA) transduces complementary sex determination (CSD) signal in the invasive Argentine ant. In this haplodiploid species, we identified a 5-kilobase hyper-polymorphic region underlying CSD: Heterozygous embryos become females, while homozygous and hemizygous embryos become males. Heterozygosity at the CSD locus correlates with higher expression of ANTSR, a gene that overlaps with the CSD locus and specifies an lncRNA transcript. ANTSR knockdown in CSD heterozygotes leads to male development. Comparative analyses indicated that, in Hymenoptera, ANTSR is an ancient yet rapidly evolving gene. This study reveals an lncRNA involved in genetic sex determination, alongside a previously unknown regulatory mechanism underlying sex determination based on complementarity among noncoding alleles.

Zygosity-based sex determination in a butterfly drives hypervariability of Masculinizer

Nature has devised many ways of producing males and females. Here, we report on a previously undescribed mechanism for Lepidoptera that functions without a female-specific gene. The number of alleles or allele heterozygosity in a single Z-linked gene (BaMasc) is the primary sex-determining switch in Bicyclus anynana butterflies. Embryos carrying a single BaMasc allele develop into WZ (or Z0) females, those carrying two distinct alleles develop into ZZ males, while (ZZ) homozygotes initiate female development, have mismatched dosage compensation, and die as embryos. Consequently, selection against homozygotes has favored the evolution of spectacular allelic diversity: 205 different coding sequences of BaMasc were detected in a sample of 246 females. The structural similarity of a hypervariable region (HVR) in BaMasc to the HVR in Apis mellifera csd suggests molecular convergence between deeply diverged insect lineages. Our discovery of this primary switch highlights the fascinating diversity of sex-determining mechanisms and underlying evolutionary drivers.

Serial founder effects slow range expansion in an invasive social insect

Invasive populations often experience founder effects: a loss of genetic diversity relative to the source population, due to a small number of founders. Even where these founder effects do not impact colonization success, theory predicts they might affect the rate at which invasive populations expand. This is because secondary founder effects are generated at advancing population edges, further reducing local genetic diversity and elevating genetic load. We show that in an expanding invasive population of the Asian honey bee (Apis cerana), genetic diversity is indeed lowest at range edges, including at the complementary sex determiner, csd, a locus that is homozygous-lethal. Consistent with lower local csd diversity, range edge colonies had lower brood viability than colonies in the range centre. Further, simulations of a newly-founded and expanding honey bee population corroborate the spatial patterns in mean colony fitness observed in our empirical data and show that such genetic load at range edges will slow the rate of population expansion.

Multivariate Models of Animal Sex: Breaking Binaries Leads to a Better Understanding of Ecology and Evolution

"Sex" is often used to describe a suite of phenotypic and genotypic traits of an organism related to reproduction. However, these traits-gamete type, chromosomal inheritance, physiology, morphology, behavior, etc.-are not necessarily coupled, and the rhetorical collapse of variation into a single term elides much of the complexity inherent in sexual phenotypes. We argue that consideration of "sex" as a constructed category operating at multiple biological levels opens up new avenues for inquiry in our study of biological variation. We apply this framework to three case studies that illustrate the diversity of sex variation, from decoupling sexual phenotypes to the evolutionary and ecological consequences of intrasexual polymorphisms. We argue that instead of assuming binary sex in these systems, some may be better categorized as multivariate and nonbinary. Finally, we conduct a meta-analysis of terms used to describe diversity in sexual phenotypes in the scientific literature to highlight how a multivariate model of sex can clarify, rather than cloud, studies of sexual diversity within and across species. We argue that such an expanded framework of "sex" better equips us to understand evolutionary processes, and that as biologists, it is incumbent upon us to push back against misunderstandings of the biology of sexual phenotypes that enact harm on marginalized communities.

Recognition of polymorphic Csd proteins determines sex in the honeybee

Sex in honeybees, Apis mellifera, is genetically determined by heterozygous versus homo/hemizygous genotypes involving numerous alleles at the single complementary sex determination locus. The molecular mechanism of sex determination is however unknown because there are more than 4950 known possible allele combinations, but only two sexes in the species. We show how protein variants expressed from complementary sex determiner (csd) gene determine sex. In females, the amino acid differences between Csd variants at the potential-specifying domain (PSD) direct the selection of a conserved coiled-coil domain for binding and protein complexation. This recognition mechanism activates Csd proteins and, thus, the female pathway. In males, the absence of polymorphisms establishes other binding elements at PSD for binding and complexation of identical Csd proteins. This second recognition mechanism inactivates Csd proteins and commits male development via default pathway. Our results demonstrate that the recognition of different versus identical variants of a single protein is a mechanism to determine sex.

Masculinizer gene controls male sex determination in the codling moth, Cydia pomonella

The molecular mechanisms of sex determination in moths and butterflies (Lepidoptera) with female heterogamety (WZ/ZZ) are poorly understood, except in the silkworm Bombyx mori. However, the Masculinizer (Masc) gene that controls male development and dosage compensation in B. mori, appears to be conserved in Lepidoptera, as its masculinizing function was recently confirmed in several moth species. In this work, we investigated the role of the Masc gene in sex determination of the codling moth Cydia pomonella (Tortricidae), a globally important pest of pome fruits and walnuts. The gene structure of the C. pomonella Masc ortholog, CpMasc, is similar to B. mori Masc. However, unlike B. mori, we identified 14 splice variants of CpMasc in the available transcriptomes. Subsequent screening for sex specificity and genetic variation using publicly available data and RT-PCR revealed three male-specific splice variants. Then qPCR analysis of these variants revealed sex-biased expression showing a peak only in early male embryos. Knockdown of CpMasc by RNAi during early embryogenesis resulted in a shift from male-to female-specific splicing of the C. pomonella doublesex (Cpdsx) gene, its downstream effector, in ZZ embryos, leading to a strongly female-biased sex ratio. These data clearly demonstrate that CpMasc functions as a masculinizing gene in the sex-determining cascade of C. pomonella. Our study also showed that CpMasc transcripts are provided maternally, as they were detected in unfertilized eggs after oviposition and in mature eggs dissected from virgin females. This finding is unique, as maternal provision of mRNA has rarely been studied in Lepidoptera.

The identification and expression pattern of the sex determination genes and their sex-specific variants in the egg parasitoid Trichogramma dendrolimi Matsumura (Hymenoptera: Trichogrammatidae)

Introduction: Trichogramma wasps are egg parasitoids of agricultural lepidopteran pests. The sex of Trichogramma is determined by its ploidy as well as certain sex ratio distorters, such as the endosymbiotic bacteria Wolbachia spp. and the paternal sex ratio (PSR) chromosome. The sex determination systems of hymenopterans, such as Trichogramma spp., involve cascades of the genes transformer (tra), transformer-2 (tra2), and doublesex (dsx) and are associated with sex-specific tra and dsx splicing. First, these genes and their sex-specific variants must be identified to elucidate the interactions between the sex ratio disorders and the sex determination mechanism of Trichogramma. Methods: Here, we characterized the sex determination genes tra, tra2, and dsx in Trichogramma dendrolimi. Sex-specific tra and dsx variants were detected in cDNA samples obtained from both male and female Trichogramma wasps. They were observed in the early embryos (1-10 h), late embryos (12-20 h), larvae (32 h and 48 h), pre-pupae (96 h), and pupae (144 h, 168 h, 192 h, and 216 h) of both male and female T. dendrolimi offspring. Results: We detected female-specific tra variants throughout the entire early female offspring stage. The male-specific variant began to express at 9-10 h as the egg was not fertilized. However, we did not find any maternally derived, female-specific tra variant in the early male embryo. This observation suggests that the female-specific tra variant expressed in the female embryo at 1-9 h may not have originated from the maternal female wasp. Discussion: The present study might be the first to identify the sex determination genes and sex-specific gene splicing in Trichogramma wasps. The findings of this study lay the foundation for investigating the sex determination mechanisms of Trichogramma and other wasps. They also facilitate sex identification in immature T. dendrolimi and the application of this important egg parasitoid in biological insect pest control programs.

Genetic methods in honey bee breeding

The honey bee Apis mellifera is a rather difficult object for selection due to the peculiarities of its biology. Breeding activities in beekeeping are aimed at obtaining bee colonies with high rates of economically useful traits, such as productivity, resistance to low temperatures and diseases, hygienic behavior, oviposition of the queen, etc. With two apiaries specializing in the breeding of A. m. mellifera and A. m. carnica as examples, the application of genetic methods in the selection of honey bees is considered. The first stage of the work was subspecies identification based on the analysis of the polymorphism of the intergenic mtDNA locus tRNAleu-COII (or COI-COII) and microsatellite nuclear DNA loci Ap243, 4a110, A24, A8, A43, A113, A88, Ap049, A28. This analysis confirmed that the studied colonies correspond to the declared subspecies. In the apiary with A. m. mellifera, hybrid colonies have been identified. A method based on the analysis of polymorphisms of the tRNAleu-COII locus and microsatellite nuclear DNA loci has been developed to identify the dark forest bee A. m. mellifera and does not allow one to differentiate subspecies from C (A. m. carnica and A. m. ligustica) and O (A. m. caucasica) evolutionary lineages from each other. The second stage was the assessment of the allelic diversity of the csd gene. In the apiary containing colonies of A. m. mellifera (N = 15), 20 csd alleles were identified. In the apiary containing colonies of A. m. carnica (N = 44), 41 alleles were identified. Six alleles are shared by both apiaries. DNA diagnostics of bee diseases showed that the studied colonies are healthy. Based on the data obtained, a scheme was developed for obtaining primary material for honey bee breeding, which can subsequently be subjected to selection according to economically useful traits. In addition, the annual assessment of the allelic diversity of the csd gene will shed light on the frequency of formation of new allelic variants and other issues related to the evolution of this gene.

SIMplyBee: an R package to simulate honeybee populations and breeding programs

BACKGROUND

The Western honeybee is an economically important species globally, but has been experiencing colony losses that lead to economical damage and decreased genetic variability. This situation is spurring additional interest in honeybee breeding and conservation programs. Stochastic simulators are essential tools for rapid and low-cost testing of breeding programs and methods, yet no existing simulator allows for a detailed simulation of honeybee populations. Here we describe SIMplyBee, a holistic simulator of honeybee populations and breeding programs. SIMplyBee is an R package and hence freely available for installation from CRAN http://cran.r-project.org/package=SIMplyBee .

IMPLEMENTATION

SIMplyBee builds upon the stochastic simulator AlphaSimR that simulates individuals with their corresponding genomes and quantitative genetic values. To enable honeybee-specific simulations, we extended AlphaSimR by developing classes for global simulation parameters, SimParamBee, for a honeybee colony, Colony, and multiple colonies, MultiColony. We also developed functions to address major honeybee specificities: honeybee genome, haplodiploid inheritance, social organisation, complementary sex determination, polyandry, colony events, and quantitative genetics at the individual- and colony-levels.

RESULTS

We describe its implementation for simulating a honeybee genome, creating a honeybee colony and its members, addressing haplodiploid inheritance and complementary sex determination, simulating colony events, creating and managing multiple colonies at the same time, and obtaining genomic data and honeybee quantitative genetics. Further documentation, available at http://www.SIMplyBee.info , provides details on these operations and describes additional operations related to genomics, quantitative genetics, and other functionalities.

DISCUSSION

SIMplyBee is a holistic simulator of honeybee populations and breeding programs. It simulates individual honeybees with their genomes, colonies with colony events, and individual- and colony-level genetic and breeding values. Regarding the latter, SIMplyBee takes a user-defined function to combine individual- into colony-level values and hence allows for modeling any type of interaction within a colony. SIMplyBee provides a research platform for testing breeding and conservation strategies and their effect on future genetic gain and genetic variability. Future developments of SIMplyBee will focus on improving the simulation of honeybee genomes, optimizing the simulator's performance, and including spatial awareness in mating functions and phenotype simulation. We invite the honeybee genetics and breeding community to join us in the future development of SIMplyBee.

Transformer gene regulates feminization under two complementary sex determination loci in the ant, Vollenhovia emeryi

Organisms that reproduce sexually have evolved well-organized mechanisms to determine two sexes. Some hymenopterans (such as ants, bees, and wasps) have a complementary sex-determination system in which heterozygosity at one CSD locus induces female development, whereas hemi- or homozygosity at the locus induces male development. This system can generate a high cost of inbreeding, as individuals that are homozygous at the locus become sterile, diploid males. On the other hand, some hymenopterans have evolved a multi-locus, complementary, sex-determination system in which heterozygosity in at least one CSD locus induces female development. This system effectively reduces the proportion of sterile diploid males; however, how these multiple
primary signals based on CSD pass through a molecular cascade to regulate downstream genes has remained unclear. To clarify this matter, we used a backcross to investigate the molecular cascade in the ant, Vollenhovia emeryi, with two CSD loci. Here we show by gene disruption that transformer (tra) is necessary for proper feminization. Expression analysis of tra and doublesex (dsx) showed that heterozygosity in at least one of the two CSD loci is sufficient to promote female sex determination. Analysis of overexpression suggested that female-type Tra protein promotes splicing of tra pre-mRNA to female isoform by a positive-regulatory-feedback loop. Our data also showed that tra affects splicing of dsx. We conclude that two-loci sex determination system in V. emeryi evolved based on tra-dsx splicing cascade that is well conserved in other insect species. Finally, we suggest a cascade model to arrive at a binary determination of sex under multiple primary signals.

Drones Do Not Drift between Nests in a Wild Population of Apis cerana

The modes through which individuals disperse prior to reproduction has important consequences for gene flow in populations. In honey bees (Apis sp.), drones (males) reproduce within a short flight range of their natal nest, leaving and returning each afternoon within a narrow mating window. Drones are assumed to return to their natal nests as they depend on workers to feed them. However, in apiaries, drones are reported to regularly make navigation errors and return to a non-natal nest, where they are accepted and fed by unrelated workers. If such a "drone drift" occurred in wild populations, it could facilitate some further degree of dispersal for males, particularly if drones drift into host nests some distance away from their natal nest. Here, we investigated whether drone drift occurs in an invasive population of the Asian honey bee (Apis cerana). Based on the genotypes of 1462 drones from 19 colonies, we found only a single drone that could be considered a candidate drifter (~0.07%). In three other colonies, drones whose genotypes differed from the inferred queen were best explained by recent queen turnover or worker-laying. We concluded that drone drift in this population is low at best, and A. cerana drones either rarely make navigation errors in wild populations or are not accepted into foreign nests when they do so. We therefore confirm that drone dispersal distance is limited to the distance of daily drone flights from natal nests, a key assumption of both colony density estimates based on sampling of drone congregation areas and population genetic models of gene flow in honey bees.

Observation of Genetic Gain with Instrumental Insemination of Honeybee Queens

Controlling mating in the honeybee (Apis mellifera) is part of one of the greatest challenges for the beekeeping industry's genetic selection programs due to specific characteristics of their reproduction. Several techniques for supervising honeybee mating with relative effective control have been developed over the years to allow honeybee selection. As part of this project, we compared the genetic gains for several colony performance traits, obtained using the BLUP-animal method, according to the selection pressure applied in controlled reproduction (directed fertilization versus instrumental insemination). Our results show similar genetic gains for hygienic behavior and honey production between colonies whether queens were fertilized naturally or via instrumental insemination, as well as similar or lower genetic gains for colonies with queens inseminated for spring development. In addition, we noticed greater fragility in queens following insemination. These findings show that instrumental insemination is an effective tool for reproductive control in genetic selection and for estimating breeding values more precisely. However, this technique does not result in queens of superior genetic quality for commercial purposes.

Comparative phylogeography reveals widespread cryptic diversity driven by ecology in Panamanian birds

Widespread species often harbor unrecognized genetic diversity, and investigating the factors associated with such cryptic variation can help us better understand the forces driving diversification. Here, we identify potential cryptic species based on a comprehensive dataset of COI mitochondrial DNA barcodes from 2,333 individual Panamanian birds across 429 species, representing 391 (59%) of the 659 resident landbird species of the country, as well as opportunistically sampled waterbirds. We complement this dataset with additional publicly available mitochondrial loci, such as ND2 and cytochrome b, obtained from whole mitochondrial genomes from 20 taxa. Using barcode identification numbers (BINs), we find putative cryptic species in 19% of landbird species, highlighting hidden diversity in the relatively well-described avifauna of Panama. Whereas some of these mitochondrial divergence events corresponded with recognized geographic features that likely isolated populations, such as the Cordillera Central highlands, the majority (74%) of lowland splits were between eastern and western populations. The timing of these splits are not temporally coincident across taxa, suggesting that historical events, such as the formation of the Isthmus of Panama and Pleistocene climatic cycles, were not the primary drivers of cryptic diversification. Rather, we observed that forest species, understory species, insectivores, and strongly territorial species-all traits associated with lower dispersal ability-were all more likely to have multiple BINs in Panama, suggesting strong ecological associations with cryptic divergence. Additionally, hand-wing index, a proxy for dispersal capability, was significantly lower in species with multiple BINs, indicating that dispersal ability plays an important role in generating diversity in Neotropical birds. Together, these results underscore the need for evolutionary studies of tropical bird communities to consider ecological factors along with geographic explanations, and that even in areas with well-known avifauna, avian diversity may be substantially underestimated.

LAY SUMMARY

- What factors are common among bird species with cryptic diversity in Panama? What role do geography, ecology, phylogeographic history, and other factors play in generating bird diversity?- 19% of widely-sampled bird species form two or more distinct DNA barcode clades, suggesting widespread unrecognized diversity.- Traits associated with reduced dispersal ability, such as use of forest understory, high territoriality, low hand-wing index, and insectivory, were more common in taxa with cryptic diversity. Filogeografía comparada revela amplia diversidad críptica causada por la ecología en las aves de Panamá.

RESUMEN

Especies extendidas frecuentemente tiene diversidad genética no reconocida, y investigando los factores asociados con esta variación críptica puede ayudarnos a entender las fuerzas que impulsan la diversificación. Aquí, identificamos especies crípticas potenciales basadas en un conjunto de datos de códigos de barras de ADN mitocondrial de 2,333 individuos de aves de Panama en 429 especies, representando 391 (59%) de las 659 especies de aves terrestres residentes del país, además de algunas aves acuáticas muestreada de manera oportunista. Adicionalmente, complementamos estos datos con secuencias mitocondriales disponibles públicamente de otros loci, tal como ND2 o citocroma b, obtenidos de los genomas mitocondriales completos de 20 taxones. Utilizando los números de identificación de código de barras (en ingles: BINs), un sistema taxonómico numérico que proporcina una estimación imparcial de la diversidad potencial a nivel de especie, encontramos especies crípticas putativas en 19% de las especies de aves terrestres, lo que destaca la diversidad oculta en la avifauna bien descrita de Panamá. Aunque algunos de estos eventos de divergencia conciden con características geográficas que probablemente aislaron las poblaciones, la mayoría (74%) de la divergencia en las tierras bajas se encuentra entre las poblaciones orientales y occidentales. El tiempo de esta divergencia no coincidió entre los taxones, sugiriendo que eventos históricos tales como la formación del Istmo de Panamá y los ciclos climáticos del pleistoceno, no fueron los principales impulsores de la especiación. En cambio, observamos asociaciones fuertes entre las características ecológicas y la divergencia mitocondriale: las especies del bosque, sotobosque, con una dieta insectívora, y con territorialidad fuerte mostraton múltiple BINs probables. Adicionalmente, el índice mano-ala, que está asociado a la capacidad de dispersión, fue significativamente menor en las especies con BINs multiples, sugiriendo que la capacidad de dispersión tiene un rol importamente en la generación de la diversidad de las aves neotropicales. Estos resultos demonstran la necesidad de que estudios evolutivos de las comunidades de aves tropicales consideren los factores ecológicos en conjunto con las explicaciones geográficos. Palabras clave: biodiversidad tropical, biogeografía, códigos de barras, dispersión, especies crípticas.

In Vitro Comparison of Sex-Specific Splicing Efficiencies of fem Pre-mRNA under Monoallelic and Heteroallelic Conditions of csd, a Master Sex-Determining Gene in the Honeybee

The sexual fate of honeybees is determined by the complementary sex determination (CSD) model: heterozygosity at a single locus (the CSD locus) determines femaleness, while hemizygosity or homozygosity at the CSD locus determines maleness. The csd gene encodes a splicing factor that regulates sex-specific splicing of the downstream target gene feminizer (fem), which is required for femaleness. The female mode of fem splicing occurs only when csd is present in the heteroallelic condition. To gain insights into how Csd proteins are only activated under the heterozygous allelic composition, we developed an in vitro assay system to evaluate the activity of Csd proteins. Consistent with the CSD model, the co-expression of two csd alleles, both of which lack splicing activity under the single-allele condition, restored the splicing activity that governs the female mode of fem splicing. RNA immunoprecipitation quantitative PCR analyses demonstrated that the CSD protein was specifically enriched in several exonic regions in the fem pre-mRNA, and enrichment in exons 3a and 5 was significantly greater under the heterozygous allelic composition than the single-allelic condition. However, in most cases csd expression under the monoallelic condition was capable of inducing the female mode of fem splicing contrary to the conventional CSD model. In contrast, repression of the male mode of fem splicing was predominant under heteroallelic conditions. These results were reproduced by real-time PCR of endogenous fem expression in female and male pupae. These findings strongly suggest that the heteroallelic composition of csd may be more important for the repression of the male splicing mode than for the induction of the female splicing mode of the fem gene.

Theme and variation in the evolution of insect sex determination

The development of dimorphic adult sexes is a critical process for most animals, one that is subject to intense selection. Work in vertebrate and insect model species has revealed that sex determination mechanisms vary widely among animal groups. However, this variation is not uniform, with a limited number of conserved factors. Therefore, sex determination offers an excellent context to consider themes and variations in gene network evolution. Here we review the literature describing sex determination in diverse insects. We have screened public genomic sequence databases for orthologs and duplicates of 25 genes involved in insect sex determination, identifying patterns of presence and absence. These genes and a 3.5 reference set of 43 others were used to infer phylogenies and compared to accepted organismal relationships to examine patterns of congruence and divergence. The function of candidate genes for roles in sex determination (virilizer, female-lethal-2-d, transformer-2) and sex chromosome dosage compensation (male specific lethal-1, msl-2, msl-3) were tested using RNA interference in the milkweed bug, Oncopeltus fasciatus. None of these candidate genes exhibited conserved roles in these processes. Amidst this variation we wish to highlight the following themes for the evolution of sex determination: (1) Unique features within taxa influence network evolution. (2) Their position in the network influences a component's evolution. Our analyses also suggest an inverse association of protein sequence conservation with functional conservation.

Transcriptome changes of Apis mellifera female embryos with fem gene knockout by CRISPR/Cas9

The sex of honey bees is decided by a regulatory cascade comprising of csd, fem and Amdsx. In order to further identify other genes involved in sex determination and differentiation of honey bees in the early stages of embryo development, the CRISPR/Cas9 method was used to knock out fem gene in the embryonic stage of diploid western honey bees, and RNA-seq was used to analyze gene expression changes in the embryo after fem knockout. Finally, we found that the bees had undergone gender changes due to fem knockout. A total of 155 differentially expressed genes (DEGs) were obtained, with 48 up-regulated and 107 down-regulated DEGs in the mutant group compared to the control group. Of them, many genes are related to sex development or differentiation. In addition, 1502 differentially expressed alternative splicing events (DEASEs) related to 1011 genes, including the main honey bee sex-determining genes csd, tra2, fem, and Amdsx, were identified between the mutant group and control group, indicating that fem regulates alternative splicing of a large number of downstream genes. Our results provide valuable clues for further investigating the molecular mechanism of sex determination and differentiation in honey bees.

Sexual dimorphism as a facilitator of worker caste evolution in ants

Ant societies are primarily composed of females, whereby labor is divided into reproductive, or queen, and non-reproductive, or worker, castes. Workers and reproductive queens can differ greatly in behavior, longevity, physiology, and morphology, but queen-worker differences are usually modest relative to the differences in males. Males are short-lived, typically do not provide the colony with labor, often look like a different species, and only occur seasonally. It is these differences that have historically led to their neglect in social insect research, but also why they may facilitate novel phenotypic variation - by increasing the phenotypic variability that is available for selection. In this study, worker variation in multivariate size-shape space paralleled male-queen variation. As worker variation increased within species, so did sexual variation. Across species in two independent genera, using head width as a proxy for body size, sexual size dimorphism correlated with worker polymorphism regardless of whether the ancestral condition was large or small worker/sexual dimorphism. Mounting molecular data support the hypothesis that queen-worker caste determination has co-opted many genes/pathways from sex determination. The molecular evidence, coupled with the observations from this study, leads to the hypothesis that sexual selection and selection on colony-level traits are non-independent, and that sexual dimorphism may even have facilitated the evolution of the distinct worker caste.

The function and evolution of a genetic switch controlling sexually dimorphic eye differentiation in honeybees

Animals develop sex-specific morphological structures that are diverse between organisms. However, understanding the developmental and evolutionary mechanisms governing these traits is still limited and largely restricted to DM domain genes, which are conserved, sex-specific developmental regulators identified in genetic models. Here, we report a sex-specific developmental regulator gene, glubschauge (glu) that selectively regulates sexually dimorphic eye differentiation in honeybees. We found that the sex determination gene feminizer (fem) controls sex-specific splicing of glu transcripts, establishing a genetic switch in which Glu proteins with a zinc finger (ZnF) domain are only expressed in females. We showed that female coding sequence was essential and sufficient for partial feminization. Comparative sequence and functional studies revealed that the evolutionary origination of the genetic switch was followed by the mutational origin of the essential ZnF domain. Our results demonstrate that glu is a newly evolved sex-specific genetic switch for region-specific regulation of a dimorphic character.

A history of the genetic and molecular identification of genes and their functions controlling insect sex determination

The genetics of the sex determination regulatory cascade in Drosophila melanogaster has a fascinating history, interlinked with the foundation of the Genetics discipline itself. The discovery that alternative splicing rather than differential transcription is the molecular mechanism underlying the upstream control of sex differences in the Drosophila model system was surprising. This notion is now fully integrated into the scientific canon, appearing in many genetics textbooks and online education resources. In the last three decades, it was a key reference point for starting evolutionary studies in other insect species by using homology-based approaches. This review will introduce a very brief history of Drosophila genetics. It will describe the genetic and molecular approaches applied for the identifying and cloning key genes involved in sex determination in Drosophila and in many other insect species. These comparative analyses led to supporting the idea that sex-determining pathways have evolved mainly by recruiting different upstream signals/genes while maintaining widely conserved intermediate and downstream regulatory genes. The review also provides examples of the link between technological advances and research achievements, to stimulate reflections on how science is produced. It aims to hopefully strengthen the related historical and conceptual knowledge of general readers of other disciplines and of younger geneticists, often focused on the latest technical-molecular approaches.

Chromatin accessibility-based characterisation of brain gene regulatory networks in three distinct honey bee polyphenisms

The honey bee genome has the capacity to produce three phenotypically distinct organisms (two diploid female castes: queen and worker, and a haploid male drone). Previous studies have implicated metabolic flux acting via epigenetic regulation in directing nutrition-driven phenotypic plasticity in the honey bee. However, the cis-acting DNA regulatory elements that establish tissue and polyphenism -specific epigenomes and gene expression programmes, remain unclear. Using a high resolution multiomic approach including assay for transposase-accessible chromatin by sequencing (ATAC-seq), RNA-seq and ChIP-seq, we produce the first genome-wide maps of the regulatory landscape across all three adult honey bee phenotypes identifying > 5000 regulatory regions in queen, 7500 in worker and 6500 in drone, with the vast majority of these sites located within intronic regions. These regions are defined by positive enrichment of H3K27ac and depletion of H3K4me3 and show a positive correlation with gene expression. Using ATAC-seq footprinting we determine queen, worker and drone -specific transcription factor occupancy and uncover novel phenotype-specific regulatory networks identifying two key nuclear receptors that have previously been implicated in caste-determination and adult behavioural maturation in honey bees; ecdysone receptor and ultraspiracle. Collectively, this study provides novel insights into key gene regulatory networks that are associated with these distinct polyphenisms in the honey bee.

Identification of potential sex determination genes and functional analyses in Neoseiulus californicus under prey stress

BACKGROUND

Phytoseiid mites are important natural enemies of spider mites. Sex-determination mechanism are important basic scientific issues in the reproduction and evolution of predatory mites. Clarifying sex-determination mechanism may provide reference for exploring genetic approach to have the phytoseiid mites produce more female offspring, which could improve their effectiveness as a biological control agent.

RESULTS

We used transcriptome sequencing to identify and characterize 20 putative sex-determination genes in the phytoseiid mite Neoseiulus californicus, a species with uncommon pseudo-arrhenotoky, including doublesex-like (dsx1-like), transformer-2 (tra-2), intersex (ix), and fruitless-like (BTB2). A significant negative correlation was found between prey stress and offspring sex ratio. But the most genes identified showed no difference in expression between the groups with lowest and highest female offspring ratios. The hatching rate and sex ratio of female offspring were reduced when the ix gene was silenced, and the oviposition days and fecundity were reduced when the BTB2 gene was silenced. The fecundity was reduced when the tra2 gene was silenced and the snf gene is essential for oviposition in female. There was no effect on reproduction and female sex determination when silencing the dsx1-like and dsx2-like gene.

CONCLUSION

The genes BTB2, tra2 and snf are involved in oviposition, and ix may be involved in female sex determination and egg formation in Neoseiulus californicus. The results are conductive to further understanding molecular regulatory mechanism of sex determination in predatory mites and may provide a reference for better use of this predatory by producing more females. © 2022 Society of Chemical Industry.

A Comparison of RNA Interference via Injection and Feeding in Honey Bees

RNA interference (RNAi) has been used successfully to reduce target gene expression and induce specific phenotypes in several species. It has proved useful as a tool to investigate gene function and has the potential to manage pest populations and reduce disease pathogens. However, it is not known whether different administration methods are equally effective at interfering with genes in bees. Therefore, we compared the effects of feeding and injection of small interfering RNA (siRNA) on the messenger RNA (mRNA) levels of alpha-aminoadipic semialdehyde dehydrogenase (ALDH7A1), 4-coumarate-CoA ligase (4CL), and heat shock protein 70 (HSP70). Both feeding and injection of siRNA successfully knocked down the gene but feeding required more siRNA than the injection. Our results suggest that both feeding and injection of siRNA effectively interfere with brain genes in bees. The appropriateness of each method would depend on the situation.

Sex- and caste-specific transcriptomes of larval honey bee (Apis mellifera L.) gonads: DMRT A2 and Hsp83 are differentially expressed and regulated by juvenile hormone

The gonads of honey bee, Apis mellifera, queens and drones are each composed of hundreds of serial units, the ovarioles and testioles, while the ovaries of the adult subfertile workers consist of only few ovarioles. We performed a comparative RNA-seq analysis on early fifth-instar (L5F1) larval gonads, which is a critical stage in gonad development of honey bee larvae. A total of 1834 genes were identified as differentially expressed (Padj < 0.01) among the three sex and caste phenotypes. The Gene Ontology analysis showed significant enrichment for metabolism, protein or ion binding, and oxidoreductase activity, and a KEGG analysis revealed metabolic pathways as enriched. In a principal component analysis for the total transcriptomes and hierarchical clustering of the DEGs, we found higher similarity between the queen and worker ovary transcriptomes compared to the drone testis, despite the onset of programmed cell death in the worker ovaries. Four DEGs were selected for RT-qPCR analyses, including their response to juvenile hormone (JH), which is a critical factor in the caste-specific development of the ovaries. Among these, DMRT A2 and Hsp83 were found upregulated by JH and, thus, emerged as potential molecular markers for sex- and caste-specific gonad development in honey bees.

Buffered fitness components: Antagonism between malnutrition and an insecticide in bumble bees

Global insect biodiversity declines due to reduced fitness are linked to interactions between environmental stressors. In social insects, inclusive fitness depends on successful mating of reproductives, i.e. males and queens, and efficient collaborative brood care by workers. Therefore, interactive effects between malnutrition and environmental pollution on sperm and feeding glands (hypopharyngeal glands (HPGs)) would provide mechanisms for population declines, unless buffered against due to their fitness relevance. However, while negative effects for bumble bee colony fitness are known, the effects of malnutrition and insecticide exposure singly and in combination on individuals are poorly understood. Here we show, in a fully-crossed laboratory experiment, that malnutrition and insecticide exposure result in neutral or antagonistic interactions for spermatozoa and HPGs of bumble bees, Bombus terrestris, suggesting strong selection to buffer key colony fitness components. No significant effects were observed for mortality and consumption, but significant negative effects were revealed for spermatozoa traits and HPGs. The combined effects on these parameters were not higher than the individual stressor effects, which indicates an antagonistic interaction between both. Despite the clear potential for additive effects, due to the individual stressors impairing muscle quality and neurological control, simultaneous malnutrition and insecticide exposure surprisingly did not reveal an increased impact compared to individual stressors, probably due to key fitness traits being resilient. Our data support that stressor interactions require empirical tests on a case-by-case basis and need to be regarded in context to understand underlying mechanisms and so adequately mitigate the ongoing decline of the entomofauna.

New insights into the criteria of functional heterozygosity of the Apis mellifera complementary sex determining gene–Discovery of a functional allele pair differing by a single amino acid

The complementary sex determiner (csd) gene is responsible for controlling the sex-determination molecular switch in western honey bees (Apis mellifera): bees that are heterozygous for csd develop into females, whereas bees that are hemizygous or homozygous develop into males. The homozygous diploid males are destroyed at an early stage of their development. It has been proposed that the minimal number of amino acid differences between two csd alleles needed to fully determine femaleness is five and it has also been shown that smaller differences may result in forming an evolutionary intermediate that is not fully capable of female determination, but has increased fitness compared to the homozygous genotype. In this study, we have implemented a terminal restriction length polymorphism-based method of identifying and distinguishing paternal alleles in a given bee colony and assigning them to a particular maternal allele in order to gather information on large number of functional csd pairs and also to identify, to some extent, genotypes that are underrepresented or absent in bee colonies. The main finding of this study is the identification of a fully functional genotype consisting of csd alleles that differed from each other by a one amino acid position. The individuals carrying this genotype expressed only female-specific transcripts of feminizer and double-sex genes. By comparing the sequences differences between the csd pair identified in our study with those described earlier, we conclude that functional heterozygosity of the csd gene is dependent not only on the number of the amino acid differences but also on the sequence context and position of the change. The discovery of a functional allele pair differing by a single amino acid also implies that the generation of a new csd specificity may also occur during a single mutation step with no need for evolutionary intermediates accumulating further mutations.

Behavioral genetics and genomics: Mendel's peas, mice, and bees

The question of the heritability of behavior has been of long fascination to scientists and the broader public. It is now widely accepted that most behavioral variation has a genetic component, although the degree of genetic influence differs widely across behaviors. Starting with Mendel's remarkable discovery of "inheritance factors," it has become increasingly clear that specific genetic variants that influence behavior can be identified. This goal is not without its challenges: Unlike pea morphology, most natural behavioral variation has a complex genetic architecture. However, we can now apply powerful genome-wide approaches to connect variation in DNA to variation in behavior as well as analyses of behaviorally related variation in brain gene expression, which together have provided insights into both the genetic mechanisms underlying behavior and the dynamic relationship between genes and behavior, respectively, in a wide range of species and for a diversity of behaviors. Here, we focus on two systems to illustrate both of these approaches: the genetic basis of burrowing in deer mice and transcriptomic analyses of division of labor in honey bees. Finally, we discuss the troubled relationship between the field of behavioral genetics and eugenics, which reminds us that we must be cautious about how we discuss and contextualize the connections between genes and behavior, especially in humans.

Gene dynamics of haplodiploidy favor eusociality in the Hymenoptera

The problem of whether haplodiploidy is responsible for the frequent evolution of eusociality in the Hymenoptera remains unresolved. The little-known "protected invasion hypothesis" posits that because a male will transmit a new allele for alloparental care to all his daughters under haplodiploidy, such an allele has a higher probability of spreading to fixation under haplodiploidy than under diploidy. This mechanism is investigated using the mating system and lifecycles ancestral to eusocial lineages. It is shown that although haplodiploidy increases the probability of fixation of a new allele, the effect is cancelled by a higher probability of the allele arising in a diploid population. However, the same effect of male haploidy results in a 30% lower threshold amount of reproductive help by a worker necessary to favor eusociality if the sex ratio of dispersing first-brood offspring remains even. This occurs because when first-brood daughters become workers, the sex ratio of dispersing first-brood offspring becomes male-biased, selecting for an overall female-biased first-brood sex ratio. Through this mechanism, haplodiploidy may favor eusociality in the absence of a female-biased sex ratio in dispersing reproductive offspring. The gene-centric approach used here reveals the critical role of male haploidy in structuring the social group.

Absence of complementary sex determination in two Leptopilina species (Figitidae, Hymenoptera) and a reconsideration of its incompatibility with endosymbiont-induced thelytoky

Complementary sex determination (CSD) is a widespread sex determination mechanism in haplodiploid Hymenoptera. Under CSD, sex is determined by the allelic state of one or multiple CSD loci. Heterozygosity at one or more loci leads to female development, whereas hemizygosity of haploid eggs and homozygosity of diploid eggs results in male development. Sexual (arrhenotokous) reproduction normally yields haploid male and diploid female offspring. Under asexual reproduction (thelytoky), diploidized unfertilized eggs develop into females. Thelytoky is often induced by bacterial endosymbionts that achieve egg diploidization by gamete duplication. As gamete duplication leads to complete homozygosity, endosymbiont-induced thelytokous reproduction is presumed to be incompatible with CSD, which relies on heterozygosity for female development. Previously, we excluded CSD in four Asobara (Braconidae) species and proposed a two-step mechanism for Wolbachia-induced thelytoky in Asobara japonica. Here, we conclusively reject CSD in two cynipid wasp species, Leptopilina heterotoma and Leptopilina clavipes. We further show that thelytoky in L. clavipes depends on Wolbachia titer but that diploidization and feminization steps cannot be separated, unlike in A. japonica. We discuss what these results reveal about the sex determination mechanism of L. clavipes and the presumed incompatibility between CSD and endosymbiont-induced thelytoky in the Hymenoptera.

Analysis of Complementary Sex-Determiner (csd) Allele Diversity in Different Honeybee Subspecies from Italy Based on NGS Data

Sexual regulation in Apis mellifera is controlled by the complementary sex-determiner (csd) gene: females (queens and workers) are heterozygous at this locus and males (drones) are hemizygous. When homozygous diploid drones develop, they are eaten by worker bees. High csd allelic diversity in honeybee populations is a priority for colony survival. The focus of this study is to investigate csd variability in the genomic sequence of the hypervariable region (HVR) of the csd gene in honeybee subspecies sampled in Italy. During the summer of 2017 and 2018, worker bees belonging to 125 colonies were sampled. The honeybees belonged to seven different A. mellifera subspecies: A. m. ligustica, A. m. sicula, A. m cecropia, A. m. carnica, A. m. mellifera, Buckfast and hybrid Carnica. Illumina genomic resequencing of all samples was performed and used for the characterization of global variability among colonies. In this work, a pipeline using existing resequencing data to explore the csd gene allelic variants present in the subspecies collection, based on de novo assembly of sequences falling within the HVR region, is described. On the whole, 138 allelic sequences were successfully reconstructed. Among these, 88 different alleles were identified, 68 of which match with csd alleles present in the NCBI GenBank database.

Highly efficient site-specific integration of DNA fragments into the honeybee genome using CRISPR/Cas9

Functional genetic studies in honeybees have been limited to transposon mediated transformation and site directed mutagenesis tools. However, site- and sequence-specific manipulations that insert DNA fragments or replace sequences at specific target sites are lacking. Such tools would enable the tagging of proteins, the expression of reporters and site-specific amino acid changes, which are all gold standard manipulations for physiological, organismal, and genetic studies. However, such manipulations must be very efficient in honeybees since screening and crossing procedures are laborious due to their social organization. Here, we report an accurate and remarkably efficient site-specific integration of DNA-sequences into the honeybee genome using clustered regularly interspaced short palindromic repeat/clustered regularly interspaced short palindromic repeat-associated protein 9-mediated homology-directed repair. We employed early embryonic injections and selected a highly efficient sgRNA in order to insert 294 and 729 bp long DNA sequences into a specific locus at the dsx gene. These sequences were locus-specifically integrated in 57% and 59% of injected bees. Most importantly, 21% and 25% of the individuals lacked the wildtype sequence demonstrating that we generated homozygous mutants in which all cells are affected (no mosaicism). The highly efficient, locus-specific insertions of nucleotide sequences generating homozygous mutants demonstrate that systematic molecular studies for honeybees are in hand that allow somatic mutation approaches via workers or studies in the next generation using queens with their worker progeny. The employment of early embryonic injections and screenings of highly efficient sgRNAs may offer the prospect of highly successful sequence- and locus-specific mutations also in other organisms.

Abundant small RNAs in the reproductive tissues and eggs of the honey bee, Apis mellifera

Background

Polyandrous social insects such as the honey bee are prime candidates for parental manipulation of gene expression in offspring. Although there is good evidence for parent-of-origin effects in honey bees the epigenetic mechanisms that underlie these effects remain a mystery. Small RNA molecules such as miRNAs, piRNAs and siRNAs play important roles in transgenerational epigenetic inheritance and in the regulation of gene expression during development.

Results

Here we present the first characterisation of small RNAs present in honey bee reproductive tissues: ovaries, spermatheca, semen, fertilised and unfertilised eggs, and testes. We show that semen contains fewer piRNAs relative to eggs and ovaries, and that piRNAs and miRNAs which map antisense to genes involved in DNA regulation and developmental processes are differentially expressed between tissues. tRNA fragments are highly abundant in semen and have a similar profile to those seen in the semen of other animals. Intriguingly we also find abundant piRNAs that target the sex determination locus, suggesting that piRNAs may play a role in honey bee sex determination.

Conclusions

We conclude that small RNAs may play a fundamental role in honey bee gametogenesis and reproduction and provide a plausible mechanism for parent-of-origin effects on gene expression and reproductive physiology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08478-9.

Sexual diversification of splicing regulation during embryonic development in honeybees (Apis mellifera), A haplodiploid system

The honeybee is a haplodiploid organism in which sexual development is determined by the complementary sex determiner (csd) gene and realized by sex-specific splicing processes involving the feminizer (fem) gene. We used high throughput transcriptome sequencing (RNA-Seq) to characterize the transcriptional differences between the sexes caused by the fertilization and sex determination processes in honeybee (Apis mellifera) embryos. We identified 758, 372 and 43 differentially expressed genes (DEGs) and 58, 176 and 233 differentially spliced genes (DSGs) in 10-15-h-old, 25-40-h-old and 55-70-h-old female and male embryos, respectively. The early difference in male and female embryos in response to the fertilization and non-fertilization processes resulted mainly in differential expression of genes (758 DEGs vs. 58 DSGs). In the latest sampled embryonic stage, the transcriptional differences between the sexes were dominated by alternative splicing of transcripts (43 DEGs vs. 233 DSGs). Interestingly, differentially spliced transcripts that encode RNA-binding properties were overrepresented in 55-70-h-old embryos, indicating a more diverse regulation via alternative splicing than previous work on the sex determination pathway suggested. These stage- and sex-specific transcriptome data from honeybee embryos provide a comprehensive resource for examining the roles of fertilization and sex determination in developmental programming in a haplodiploid system.

A feminizing switch in a hemimetabolous insect

The mechanism of sex determination remains poorly understood in hemimetabolous insects. Here, in the brown planthopper (BPH), Nilaparvata lugens, a hemipteran rice pest, we identified a feminizing switch or a female determiner (Nlfmd) that encodes a serine/arginine-rich protein. Knockdown of Nlfmd in female nymphs resulted in masculinization of both the somatic morphology and doublesex splicing. The female-specific isoform of Nlfmd, Nlfmd-F, is maternally deposited and zygotically transcribed. Depletion of Nlfmd by maternal RNAi or CRISPR-Cas9 resulted in female-specific embryonic lethality. Knockdown of an hnRNP40 family gene named female determiner 2 (Nlfmd2) also conferred masculinization. In vitro experiments showed that an Nlfmd2 isoform, NlFMD2³⁴⁰, bound the RAAGAA repeat motif in the Nldsx pre-mRNA and formed a protein complex with NlFMD-F to modulate Nldsx splicing, suggesting that NlFMD2 may function as an RNA binding partner of the feminizing switch NlFMD. Our results provide novel insights into the diverse mechanisms of insect sex determination.

Practical Applications of Genomics in Managing Honey bee Health

The honey bee Apis mellifera is a model organism for sociogenomics and one of the most important managed pollinators. High mortalities experienced by honey bee colonies over the past several decades are expected to have a substantive effect on crop pollination and global food security. These threats and the availability of a growing number of genomic resources for the honey bee have motivated research on how genetics and genomics can be practically applied to manage bee health. The authors review 3 such applications: (1) Certification of bee lineages using single-polymorphism markers; (2) breeding bees using marker-assisted selection; (3) diagnosing honey bee stressors using biomarkers.

Female developmental environment delays development of male honeybee (Apis mellifera)

Background

Nutrition and cell size play an important role in the determination of caste differentiation in queen and worker of honeybees (Apis mellifera), whereas the haploid genome dominates the differentiation of drones. However, the effects of female developmental environment on the development of males remain unclear. In this study, young drone larvae were transferred into worker cells (WCs) or remained in drone cells (DCs) to rear drones. The drone larvae were also grafted into queen cells (QCs) for 48 h and then transplanted into drone cells until emerging. Morphological indexes and reproductive organs of these three types of newly emerged drones were measured. Newly emerged drones and third instar drone larvae from WCs, DCs and QCs were sequenced by RNA sequencing (RNA-Seq).

Results

The amount of food remaining in cells of the QC and WC groups was significantly different to that in the DC group at the early larval stage. Morphological results showed that newly emerged DC drones had bigger body sizes and more well-developed reproductive tissues than WC and QC drones, whereas the reproductive tissues of QC drones were larger than those of WC drones. Additionally, whole body gene expression results showed a clear difference among three groups. At larval stage there were 889, 1761 and 1927 significantly differentially expressed genes (DEGs) in WC/DC, QC/DC and WC/QC comparisons, respectively. The number of DEGs decreased in adult drones of these three comparisons [678 (WC/DC), 338 (QC/DC) and 518 (WC/QC)]. A high number of DEGs were involved in sex differentiation, growth, olfaction, vision, mammalian target of rapamycin (mTOR), Wnt signaling pathways, and other processes.

Conclusions

This study demonstrated that the developmental environment of honeybee females can delay male development, which may serve as a model for understanding the regulation of sex differentiation and male development in social insects by environmental factors.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08014-1.

Application of Next Generation Semiconductor-Based Sequencing for the Identification of Apis mellifera Complementary Sex Determiner (csd) Alleles from Honey DNA

Simple Summary

Honey contains traces of the DNA of the honey bees that produced it. This environmental DNA can therefore be used to investigate the genome of the honey bees. In this study, we used a next generation sequencing technology to analyze the variability of a key gene of Apis mellifera L., the complementary sex determiner (csd) gene, using honey environmental DNA as a source of honey bee DNA. This gene determines the sex of the bees. Two different alleles at this locus are needed to produce females whereas males have only one copy of this gene as they are haploid. In case two identical alleles are present in a diploid individual, the larvae are not vital and are discarded by the workers. Therefore, there is an advantage in maintaining a large csd diversity in honey bee populations. In light of the recent decline in honey bee populations, it is important to monitor the allele variability at this gene. The applied methodology provided a new strategy to disclose the genetic diversity at the csd gene at the population-wide level and identify most, if not all, csd alleles present in the colonies in a single analysis.

Abstract

The complementary sex determiner (csd) gene plays an essential role in the sex determination of Apis mellifera L. Females develop only if fertilized eggs have functional heterozygous genotypes at this gene whereas males, being haploids, are hemizygous. Two identical csd alleles produce non vital males. In light of the recent decline in honey bee populations, it is therefore important to monitor the allele variability at this gene. In this study, we tested the application of next generation semiconductor-based sequencing technology (Ion Torrent) coupled with environmental honey DNA as a source of honey bee genome information to retrieve massive sequencing data for the analysis of variability at the hypervariable region (HVR) of the csd gene. DNA was extracted from 12 honey samples collected from honeycombs directly retrieved from 12 different colonies. A specifically designed bioinformatic pipeline, applied to analyze a total of about 1.5 million reads, identified a total of 160 different csd alleles, 55% of which were novel. The average number of alleles per sample was compatible with the number of expected patrilines per colony, according to the mating behavior of the queens. Allele diversity at the csd could also provide information useful to reconstruct the history of the honey.

A simulation study of a honeybee breeding scheme accounting for polyandry, direct and maternal effects on colony performance

BACKGROUND

Efficient breeding programs are difficult to implement in honeybees due to their biological specificities (polyandry and haplo-diploidy) and complexity of the traits of interest, with performances being measured at the colony scale and resulting from the joint effects of tens of thousands of workers (called direct effects) and of the queen (called maternal effects). We implemented a Monte Carlo simulation program of a breeding plan designed specifically for Apis mellifera's populations to assess the impact of polyandry versus monoandry on colony performance, inbreeding level and genetic gain depending on the individual selection strategy considered, i.e. complete mass selection or within-family (maternal lines) selection. We simulated several scenarios with different parameter setups by varying initial genetic variances and correlations between direct and maternal effects, the selection strategy and the polyandry level. Selection was performed on colony phenotypes.

RESULTS

All scenarios showed strong increases in direct breeding values of queens after 20 years of selection. Monoandry led to significantly higher direct than maternal genetic gains, especially when a negative correlation between direct and maternal effects was simulated. However, the relative increase in these genetic gains depended also on their initial genetic variability and on the selection strategy. When polyandry was simulated, the results were very similar with either 8 or 16 drones mated to each queen. Across scenarios, polyandrous mating resulted in equivalent or higher gains in performance than monoandrous mating, but with considerably lower inbreeding rates. Mass selection conferred a ~ 20% increase in performance compared to within-family selection, but was also accompanied by a strong increase in inbreeding levels (25 to 50% higher).

CONCLUSIONS

Our study is the first to compare the long-term effects of polyandrous versus monoandrous mating in honeybee breeding. The latter is an emergent strategy to improve specific traits, such as resistance to varroa, which can be difficult or expensive to phenotype. However, if used during several generations in a closed population, monoandrous mating increases the inbreeding level of queens much more than polyandrous mating, which is a strong limitation of this strategy.

Eusocial insect declines: Insecticide impairs sperm and feeding glands in bumblebees

Insecticides are contributing to global insect declines, thereby creating demand to understand the mechanisms underlying reduced fitness. In the eusocial Hymenoptera, inclusive fitness depends on successful mating of male sexuals (drones) and efficient collaborative brood care by female workers. Therefore, sublethal insecticide effects on sperm and glands used in larval feeding (hypopharyngeal glands (HPG)) would provide key mechanisms for population declines in eusocial insects. However, while negative impacts for bumblebee colony fitness have been documented, the effects of insecticide exposure on individual physiology are less well understood. Here, we show that field-realistic concentrations (4.5-40 ng ml^-1) of the neonicotinoid insecticide thiamethoxam significantly impair Bombus terrestris sperm and HPGs, thereby providing plausible mechanisms underlying bumblebee population decline. In the laboratory, drones and workers were exposed to five thiamethoxam concentrations (4.5 to 1000 ng ml^-1). Then, survival, food consumption, body mass, HPG development, sperm quantity and viability were assessed. At all concentrations, drones were more exposed than workers due to higher food consumption. Increased body mass was observed in drones starting at 20 ng ml^-1 and in workers at 100 ng ml^-1. Furthermore, environmentally realistic concentrations (4.5-40 ng ml^-1) did not significantly affect survival or consumption for either sex. However, thiamethoxam exposure significantly negatively affected both sperm viability and HPG development at all tested concentrations. Therefore, the results indicate a trade-off between survival and fitness components, possibly due to costly detoxification. Since sperm and HPG are corner stones of colony fitness, the data offer plausible mechanisms for bumblebee population declines. To adequately mitigate ongoing biodiversity declines for the eusocial insects, this study suggests it is essential to evaluate the impact of insecticides on fitness parameters of both sexuals and workers.

The neonicotinoid thiamethoxam impairs male fertility in solitary bees, Osmia cornuta

The ongoing loss of global biodiversity is endangering ecosystem functioning and human food security. While environmental pollutants are well known to reduce fertility, the potential effects of common neonicotinoid insecticides on insect fertility remain poorly understood. Here, we show that field-realistic neonicotinoid exposure can drastically impact male insect fertility. In the laboratory, male and female solitary bees Osmia cornuta were exposed to four concentrations of the neonicotinoid thiamethoxam to measure survival, food consumption, and sperm traits. Despite males being exposed to higher dosages of thiamethoxam, females revealed an overall increased hazard rate for survival; suggesting sex-specific differences in toxicological sensitivity. All tested sublethal concentrations (i.e., 1.5, 4.5 and 10 ng g^-1) reduced sperm quantity by 57% and viability by 42% on average, with the lowest tested concentration leading to a reduction in total living sperm by 90%. As the tested sublethal concentrations match estimates of global neonicotinoid pollution, this reveals a plausible mechanism for population declines, thereby reflecting a realistic concern. An immediate reduction in environmental pollutants is required to decelerate the ongoing loss of biodiversity.

Societies to genes: can we get there from here?

Understanding the organization and evolution of social complexity is a major task because it requires building an understanding of mechanisms operating at different levels of biological organization from genes to social interactions. I discuss here, a unique forward genetic approach spanning more than 30 years beginning with human-assisted colony-level selection for a single social trait, the amount of pollen honey bees (Apis mellifera L.) store. The goal was to understand a complex social trait from the social phenotype to genes responsible for observed trait variation. The approach combined the results of colony-level selection with detailed studies of individual behavior and physiology resulting in a mapped, integrated phenotypic architecture composed of correlative relationships between traits spanning anatomy, physiology, sensory response systems, and individual behavior that affect individual foraging decisions. Colony-level selection reverse engineered the architecture of an integrated phenotype of individuals resulting in changes in the social trait. Quantitative trait locus (QTL) studies combined with an exceptionally high recombination rate (60 kb/cM), and a phenotypic map, provided a genotype-phenotype map of high complexity demonstrating broad QTL pleiotropy, epistasis, and epistatic pleiotropy suggesting that gene pleiotropy or tight linkage of genes within QTL integrated the phenotype. Gene expression and knockdown of identified positional candidates revealed genes affecting foraging behavior and confirmed one pleiotropic gene, a tyramine receptor, as a target for colony-level selection that was under selection in two different tissues in two different life stages. The approach presented here has resulted in a comprehensive understanding of the structure and evolution of honey bee social organization.

Sex-Specific Development in Haplodiploid Honeybee Is Controlled by the Female-Embryo-Specific Activation of Thousands of Intronic LncRNAs

Embryonic development depends on a highly coordinated shift in transcription programs known as the maternal-to-zygotic transition (MZT). It remains unclear how haploid and diploid embryo coordinate their genomic activation and embryonic development during MZT in haplodiploid animals. Here, we applied a single-embryo RNA-seq approach to characterize the embryonic transcriptome dynamics in haploid males vs. diploid females of the haplodiploid insect honeybee (Apis mellifera). We observed typical zygotic genome activation (ZGA) occurred in three major waves specifically in female honeybee embryos; haploid genome activation was much weaker and occurred later. Strikingly, we also observed three waves of transcriptional activation for thousands of long non-coding transcripts (lncRNA), 73% of which are transcribed from intronic regions and 65% were specific to female honeybee embryos. These findings support a model in which introns encode thousands of lncRNAs that are expressed in a diploid-embryo-specific and ZGA-triggered manner that may have potential functions to regulate gene expression during early embryonic development in the haplodiploid insect honeybee.

Transcriptional Profiles of Diploid Mutant Apis mellifera Embryos after Knockout of csd by CRISPR/Cas9

Simple Summary

In honey bees, males are haploid while females are diploid, leading to a fundamental difference in genetic materials between the sexes. In order to better control the comparison of gene expression between males and females, diploid mutant males were generated by knocking out the sex-determining gene, complementary sex determiner (csd), in fertilized embryos. The diploid mutant drones had male external morphological features, as well as male gonads. RNA sequencing was performed on the diploid mutant embryos and one-day-old larvae. The transcriptome analysis showed that several female-biased genes, such as worker-enriched antennal (Wat), vitellogenin (Vg), and some venom-related genes, were down-regulated in the diploid mutant males. In contrast, some male-biased genes, like takeout and apolipophorin-III-like protein (A4), were up-regulated. Moreover, the co-expression gene networks suggested that csd might interact very closely with fruitless (fru), feminizer (fem) might have connections with hexamerin 70c (hex70c), and transformer-2 (tra2) might play roles with troponin T (TpnT). Foundational information about the differences in the gene expression caused by sex differentiation was provided in this study. It is believed that this study will pave the ground for further research on the different mechanisms between males and females in honey bees.

Abstract

In honey bees, complementary sex determiner (csd) is the primary signal of sex determination. Its allelic composition is heterozygous in females, and hemizygous or homozygous in males. To explore the transcriptome differences after sex differentiation between males and females, with genetic differences excluded, csd in fertilized embryos was knocked out by CRISPR/Cas9. The diploid mutant males at 24 h, 48 h, 72 h, and 96 h after egg laying (AEL) and the mock-treated females derived from the same fertilized queen were investigated through RNA-seq. Mutations were detected in the target sequence in diploid mutants. The diploid mutant drones had typical male morphological characteristics and gonads. Transcriptome analysis showed that several female-biased genes, such as worker-enriched antennal (Wat), vitellogenin (Vg), and some venom-related genes, were down-regulated in the diploid mutant males. In contrast, some male-biased genes, such as takeout and apolipophorin-III-like protein (A4), had higher expressions in the diploid mutant males. Weighted gene co-expression network analysis (WGCNA) indicated that there might be interactions between csd and fruitless (fru), feminizer (fem) and hexamerin 70c (hex70c), transformer-2 (tra2) and troponin T (TpnT). The information provided by this study will benefit further research on the sex dimorphism and development of honey bees and other insects in Hymenoptera.

A conserved role of the duplicated Masculinizer gene in sex determination of the Mediterranean flour moth, Ephestia kuehniella

Sex determination in the silkworm, Bombyx mori, is based on Feminizer (Fem), a W-linked Fem piRNA that triggers female development in WZ individuals, and the Z-linked Masculinizer (Masc), which initiates male development and dosage compensation in ZZ individuals. While Fem piRNA is missing in a close relative of B. mori, Masc determines sex in several representatives of distant lepidopteran lineages. We studied the molecular mechanisms of sex determination in the Mediterranean flour moth, Ephestia kuehniella (Pyralidae). We identified an E. kuehniella Masc ortholog, EkMasc, and its paralog resulting from a recent duplication, EkMascB. Both genes are located on the Z chromosome and encode a similar Masc protein that contains two conserved domains but has lost the conserved double zinc finger domain. We developed PCR-based genetic sexing and demonstrated a peak in the expression of EkMasc and EkMascB genes only in early male embryos. Simultaneous knock-down experiments of both EkMasc and EkMascB using RNAi during early embryogenesis led to a shift from male- to female-specific splicing of the E. kuehniella doublesex gene (Ekdsx), their downstream effector, in ZZ embryos and resulted in a strong female-biased sex-ratio. Our results thus confirmed the conserved role of EkMasc and/or EkMascB in masculinization. We suggest that the C-terminal proline-rich domain, we have identified in all functionally confirmed Masc proteins, in conjunction with the masculinizing domain, is important for transcriptional regulation of sex determination in Lepidoptera. The function of the Masc double zinc finger domain is still unknown, but appears to have been lost in E. kuehniella.

The sex-determining cascade in the silkworm, Bombyx mori, differs greatly from those of other insects. In B. mori, female development is initiated by Fem piRNA expressed from the W chromosome during early embryogenesis. Fem piRNA silences Masculinizer (Masc) thereby blocking the male pathway resulting in female development. It is currently unknown whether this cascade is conserved across Lepidoptera. In the Mediterranean flour moth, Ephestia kuehniella, we identified an ortholog of Masc and discovered its functional duplication on the Z chromosome, which has not yet been found in any other lepidopteran species. We provide two lines of evidence that the EkMasc and/or EkMascB genes play an essential role in masculinization: (i) they show a peak of expression during early embryogenesis in ZZ but not in WZ embryos and (ii) their simultaneous silencing by RNAi results in female-specific splicing of the E. kuehniella doublesex gene (Ekdsx) in ZZ embryos and in a female-biased sex ratio. Our results suggest a conserved role of the duplicated Masc gene in sex determination of E. kuehniella.

Identification and Functional Analysis of the Sex-Determiner Transformer-2 Homologue in the Freshwater Pearl Mussel, Hyriopsis cumingii

Transformer-2 (Tra-2) is an upstream regulatory element of the sex regulation mechanism in insects and plays a critical role in sex formation. To understand the role of tra-2 in Hyriopsis cumingii, the full-length Hctra-2 (1867 bp) was obtained from the gonads, and sequence alignment with other species showed that HCTRA-2 protein had a highly conserved RRM domain. Phylogenetic analysis showed that the HCTRA-2 protein was a close relative to of the mollusks TRA-2 protein. The qRT-PCR of tissue-specific expression pattern showed that the Hctra-2 was abundant in gonads, and the expression in testes was higher than that in ovaries (p < 0.01). It suggests that Hctra-2 may play a potential regulatory role in gonadal development of H. cumingii. In the early gonadal development, the Hctra-2 expression was the highest on the third day after fertilization and increased slightly from 4 months to 5 months, which may be related to the embryonic sex determination and early gonadal development. In situ hybridization showed that Hctra-2 mRNA signals were present in both male and female gonads. After silencing Hctra-2 by RNAi, the expression levels of Hcfem-1b and Hcdmrt were changed. It is speculated that there may be a certain relationship between them, which plays an important role in the sex regulation of H. cumingii. Our research will help to deepen our understanding of the shellfish sex determination mechanisms.

Comparative Transcriptome Analysis Reveals Sex-Based Differences during the Development of the Adult Parasitic Wasp Cotesia vestalis (Hymenoptera: Braconidae)

The endoparasitic wasp Cotesia vestalis is an important biological agent for controlling the population of Plutella xylostella, a major pest of cruciferous crops worldwide. Though the genome of C. vestalis has recently been reported, molecular mechanisms associated with sexual development have not been comprehensively studied. Here, we combined PacBio Iso-Seq and Illumina RNA-Seq to perform genome-wide profiling of pharate adult and adult development of male and female C. vestalis. Taking advantage of Iso-Seq full-length reads, we identified 14,466 novel transcripts as well as 8770 lncRNAs, with many lncRNAs showing a sex- and stage-specific expression pattern. The differentially expressed gene (DEG) analyses showed 2125 stage-specific and 326 sex-specific expressed genes. We also found that 4819 genes showed 11,856 alternative splicing events through combining the Iso-Seq and RNA-Seq data. The results of comparative analyses showed that most genes were alternatively spliced across developmental stages, and alternative splicing (AS) events were more prevalent in females than in males. Furthermore, we identified six sex-determining genes in this parasitic wasp and verified their sex-specific alternative splicing profiles. Specifically, the characterization of feminizer and doublesex splicing between male and female implies a conserved regulation mechanism of sexual development in parasitic wasps.

Adaptive, caste-specific changes to recombination rates in a thelytokous honeybee population

The ability to clone oneself has clear benefits-no need for mate hunting or dilution of one's genome in offspring. It is therefore unsurprising that some populations of haplo-diploid social insects have evolved thelytokous parthenogenesis-the virgin birth of a female. But thelytokous parthenogenesis has a downside: the loss of heterozygosity (LoH) as a consequence of genetic recombination. LoH in haplo-diploid insects can be highly deleterious because female sex determination often relies on heterozygosity at sex-determining loci. The two female castes of the Cape honeybee, Apis mellifera capensis, differ in their mode of reproduction. While workers always reproduce thelytokously, queens always mate and reproduce sexually. For workers, it is important to reduce the frequency of recombination so as to not produce offspring that are homozygous. Here, we ask whether recombination rates differ between Cape workers and Cape queens that we experimentally manipulated to reproduce thelytokously. We tested our hypothesis that Cape workers have evolved mechanisms that restrain genetic recombination, whereas queens have no need for such mechanisms because they reproduce sexually. Using a combination of microsatellite genotyping and whole-genome sequencing we find that a reduction in recombination is confined to workers only.

Conservation insights from wild bee genetic studies: Geographic differences, susceptibility to inbreeding, and signs of local adaptation

Conserving bees are critical both ecologically and economically. Genetic tools are valuable for monitoring these vital pollinators since tracking these small, fast-flying insects by traditional means is difficult. By surveying the current state of the literature, this review discusses how recent advances in landscape genetic and genomic research are elucidating how wild bees respond to anthropogenic threats. Current literature suggests that there may be geographic differences in the vulnerability of bee species to landscape changes. Populations of temperate bee species are becoming more isolated and more genetically depauperate as their landscape becomes more fragmented, but tropical bee species appear unaffected. These differences may be an artifact of historical differences in land-use, or it suggests that different management plans are needed for temperate and tropical bee species. Encouragingly, genetic studies on invasive bee species indicate that low levels of genetic diversity may not lead to rapid extinction in bees as once predicted. Additionally, next-generation sequencing has given researchers the power to identify potential genes under selection, which are likely critical to species' survival in their rapidly changing environment. While genetic studies provide insights into wild bee biology, more studies focusing on a greater phylogenetic and life-history breadth of species are needed. Therefore, caution should be taken when making broad conservation decisions based on the currently few species examined.

Evolution of sexual development and sexual dimorphism in insects

Most animal species consist of two distinct sexes. At the morphological, physiological, and behavioral levels the differences between males and females are numerous and dramatic, yet at the genomic level they are often slight or absent. This disconnect is overcome because simple genetic differences or environmental signals are able to direct the sex-specific expression of a shared genome. A canonical picture of how this process works in insects emerged from decades of work on Drosophila. But recent years have seen an explosion of molecular-genetic and developmental work on a broad range of insects. Drawing these studies together, we describe the evolution of sexual dimorphism from a comparative perspective and argue that insect sex determination and differentiation systems are composites of rapidly evolving and highly conserved elements.

The direct regulation of Aalbdsx on AalVgR is indispensable for ovarian development in Aedes albopictus

BACKGROUND

Aedes albopictus is an important vector with an extensive worldwide distribution. Only female mosquitoes play a significant role in the transmission of pathogens. Doublesex (dsx) is a central nexus gene in the insect somatic sex determination hierarchy.

RESULTS

In this study, we characterized the full-length sex-specific splicing forms of the Ae. albopictus dsx (Aalbdsx) gene. Then, we identified 15 direct target genes of DSX in adult females using digital gene expression combined with quantitative real-time polymerase chain reaction (qPCR) by performing a chromatin immunoprecipitation (ChIP) assay with specific DSX antibodies. Knockdown of Aalbdsx suppressed ovarian development and decreased the transcript levels of the Aalbdsx target vitellogenin receptor (VgR) gene, whereas vitellogenin (Vg) expression showed an increase in the fat body. Genes in the major Vg regulatory pathway were also up-regulated.

CONCLUSION

Our results suggest that both Vg and VgR are direct target genes of Aalbdsx and that direct regulation of Aalbdsx on VgR is indispensable for ovarian development in Ae. albopictus, which not only provides a reference for the further elucidation of the evolutionarily conserved role of dsx in Ae. albopictus sexual differentiation but also reveals potential molecular targets for application to the development of sterile male mosquitoes to be released for vector control. © 2020 Society of Chemical Industry.