Patterns of conservation and change in honey bee developmental genes

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.

Seeking Sense in the Hox Gene Cluster

The Hox gene cluster, responsible for patterning of the head-tail axis, is an ancestral feature of all bilaterally symmetrical animals (the Bilateria) that remains intact in a wide range of species. We can say that the Hox cluster evolved successfully only once since it is commonly the same in all groups, with labial-like genes at one end of the cluster expressed in the anterior embryo, and Abd-B-like genes at the other end of the cluster expressed posteriorly. This review attempts to make sense of the Hox gene cluster and to address the following questions. How did the Hox cluster form in the protostome-deuterostome last common ancestor, and why was this with a particular head-tail polarity? Why is gene clustering usually maintained? Why is there collinearity between the order of genes along the cluster and the positions of their expressions along the embryo? Why do the Hox gene expression domains overlap along the embryo? Why have vertebrates duplicated the Hox cluster? Why do Hox gene knockouts typically result in anterior homeotic transformations? How do animals adapt their Hox clusters to evolve new structural patterns along the head-tail axis?

Genome-Wide Identification and Expression Profiling of the Wnt Gene Family in Three Rice Planthoppers: Sogatella furcifera, Laodelphax striatellus, and Nilaparvata lugens

The Wnt gene family plays essential roles in regulating many developmental processes, including the maintenance of stem cells, cell division, and cell migration. The number of Wnt genes varies among species. Due to the diversity and importance of their functions, the Wnt gene family has gained extensive research interest in various animal species from invertebrates to vertebrates. However, knowledge of the Wnt gene family is limited in rice planthoppers. Three planthopper species, the white-backed planthopper (Sogatella furcifera Horvath), the small brown planthopper (Laodelphax striatellus Fallén) and the brown planthopper (Nilaparvata lugens Stål) (Hemiptera: Delphacidae), are devastating specialist pests of rice and cause serious damage to rice plants. To better study the evolution and function of the Wnt gene family in rice planthoppers, we identified 8 Wnt family genes in three rice planthoppers with both genomic and extensive transcriptomic resources available. We conducted a systematic analysis of the three kinds of rice planthoppers and analyzed the dynamic patterns of gene conservation, as well as Wnt gene loss and duplication. The expression profiles in different developmental stages of S. furcifera and different adult organs and tissues of L. striatellus provide preliminary functional implications for the Wnt genes in rice planthopper. This study presents the first genome-wide study of the Wnt gene family in rice planthoppers, and our findings provide insights into Wnt function and evolution in rice planthoppers.

The organizing role of Wnt signaling pathway during arthropod posterior growth

Wnt signaling pathways are recognized for having major roles in tissue patterning and cell proliferation. In the last years, remarkable progress has been made in elucidating the molecular and cellular mechanisms that underlie sequential segmentation and axial elongation in various arthropods, and the canonical Wnt pathway has emerged as an essential factor in these processes. Here we review, with a comparative perspective, the current evidence concerning the participation of this pathway during posterior growth, its degree of conservation among the different subphyla within Arthropoda and its relationship with the rest of the gene regulatory network involved. Furthermore, we discuss how this signaling pathway could regulate segmentation to establish this repetitive pattern and, at the same time, probably modulate different cellular processes precisely coupled to axial elongation. Based on the information collected, we suggest that this pathway plays an organizing role in the formation of the body segments through the regulation of the dynamic expression of segmentation genes, via controlling the caudal gene, at the posterior region of the embryo/larva, that is necessary for the correct sequential formation of body segments in most arthropods and possibly in their common segmented ancestor. On the other hand, there is insufficient evidence to link this pathway to axial elongation by controlling its main cellular processes, such as convergent extension and cell proliferation. However, conclusions are premature until more studies incorporating diverse arthropods are carried out.

Deformed Mediated Larval Incisor Lobe Development Causes Differing Feeding Behavior between Oriental Armyworm and Fall Armyworm

Mandibular incisor lobes are important for insect feeding behavior, living habits and niche. However, the molecular regulation of insect incisor lobe development remains unknown. In this study, we found that two maize pests, oriental armyworm Mythimna separata and fall armyworm Spodoptera frugiperda, have different feeding patterns in maize, which are closely associated with the different development patterns of their incisor lobes. Different from first to sixth instar S. frugiperda, which feed on leaf tissues and whorls with sharp incisor lobes, older instars of M. separata feed from leaf margins with no incisor lobes. Hox gene Deformed (Dfd) is important for head appendages, but its function in incisor lobe development is not clear. Here, Dfds were identified from two armyworm species, and both were expressed highly in heads and eggs. Interestingly, the expression levels of MsDfd were relatively high in larval mandibles and decreased dramatically from fourth-instar mandibles in M. separata. Knockdown of MsDfd resulted in malformed mandibles with no incisor lobe in M. separata, making the larvae unable to perform window-feeding. However, RNAi of SfDfd did not affect the mandibles and window-feeding pattern of S. frugiperda, indicating the different roles of Dfd in these two species. Moreover, the mortality of new first instar M. separata increased after feeding dsMsDfd but did not for S. frugiperda feeding dsSfDfd. These findings revealed that Dfd mediated the larval mandibular incisor lobe morphology, affecting its feeding pattern in M. separata, broadening the knowledge of Dfd functions in insect mandibles and feeding behavior.

Annotation of Hox cluster and Hox cofactor genes in the Asian citrus psyllid, Diaphorina citri, reveals novel features

Hox genes and their cofactors are essential developmental genes specifying regional identity in animals. Hox genes have a conserved arrangement in clusters in the same order in which they specify identity along the anterior-posterior axis. A few insect species have breaks in the cluster, but these are exceptions. We annotated the 10 Hox genes of the Asian citrus psyllid Diaphorina citri, and found a split in its Hox cluster between the Deformed and Sex combs reduced genes - the first time a break at this position has been observed in an insect Hox cluster. We also annotated D. citri orthologs of the Hox cofactor genes homothorax, PKNOX and extradenticle and found an additional copy of extradenticle in D. citri that appears to be a retrogene. Expression data and sequence conservation suggest that the extradenticle retrogene may have retained the original extradenticle function and allowed divergence of the parental extradenticle gene.

The Nasonia pair-rule gene regulatory network retains its function over 300 million years of evolution

Insect segmentation is a well-studied and tractable system with which to investigate the genetic regulation of development. Though insects segment their germband using a variety of methods, modelling work implies that a single gene regulatory network can underpin the two main types of insect segmentation. This means limited genetic changes are required to explain significant differences in segmentation mode between different insects. This idea needs to be tested in a wider variety of species, and the nature of the gene regulatory network (GRN) underlying this model has not been tested. Some insects, e.g. Nasonia vitripennis and Apis mellifera segment progressively, a pattern not examined in previous studies of this segmentation model, producing stripes at different times progressively through the embryo, but not from a segment addition zone. Here, we aim to understand the GRNs patterning Nasonia using a simulation-based approach. We found that an existing model of Drosophila segmentation ( Clark, 2017) can be used to recapitulate the progressive segmentation of Nasonia, if provided with altered inputs in the form of expression of the timer genes Nv-caudal and Nv-odd paired. We predict limited topological changes to the pair-rule network and show, by RNAi knockdown, that Nv-odd paired is required for morphological segmentation. Together this implies that very limited changes to the Drosophila network are required to simulate Nasonia segmentation, despite significant differences in segmentation modes, implying that Nasonia use a very similar version of an ancestral GRN used by Drosophila, which must therefore have been conserved for at least 300 million years.

Faster-haplodiploid evolution under divergence-with-gene-flow: simulations and empirical data from pine-feeding hymenopterans

Although haplodiploidy is widespread in nature, the evolutionary consequences of this mode of reproduction are not well characterized. Here, we examine how genome-wide hemizygosity and a lack of recombination in haploid males affects genomic differentiation in populations that diverge via natural selection while experiencing gene flow. First, we simulated diploid and haplodiploid "genomes" (500-kb loci) evolving under an isolation-with-migration model with mutation, drift, selection, migration, and recombination; and examined differentiation at neutral sites both tightly and loosely linked to a divergently selected site. So long as there is divergent selection and migration, sex-limited hemizygosity and recombination cause elevated differentiation (i.e., produce a "faster-haplodiploid effect") in haplodiploid populations relative to otherwise equivalent diploid populations, for both recessive and codominant mutations. Second, we used genome-wide SNP data to model divergence history and describe patterns of genomic differentiation between sympatric populations of Neodiprion lecontei and N. pinetum, a pair of pine sawfly species (order: Hymenoptera; family: Diprionidae) that are specialized on different pine hosts. These analyses support a history of continuous gene exchange throughout divergence and reveal a pattern of heterogeneous genomic differentiation that is consistent with divergent selection on many unlinked loci. Third, using simulations of haplodiploid and diploid populations evolving according to the estimated divergence history of N. lecontei and N. pinetum, we found that divergent selection would lead to higher differentiation in haplodiploids. Based on these results, we hypothesize that haplodiploids undergo divergence-with-gene-flow and sympatric speciation more readily than diploids.

Two quantitative trait loci are associated with recapping of Varroa destructor-infested brood cells in Apis mellifera mellifera

Recapping of Varroa destructor-infested brood cells is a trait that has recently attracted interest in honey bee breeding to select mite-resistant Apis mellifera colonies. To investigate the genetic architecture of this trait, we evaluated a sample of A. mellifera mellifera colonies (N = 155) from Switzerland and France and performed a genome-wide association study, using a pool of 500 workers per colony for next-generation sequencing. The results revealed that two QTL were significantly (P < 0.05) associated with recapping of V. destructor-infested brood cells. The best-associated QTL is located on chromosome 5 in a region previously found to be associated with grooming behaviour, a resistance trait against V. destructor, in A. mellifera and Apis cerana. The second best-associated QTL is located on chromosome 4 in an intron of the Dscam gene, which is involved in neuronal wiring. Previous research demonstrated that genes involved in neuronal wiring are associated with recapping and varroa sensitive hygiene. Therefore, our study confirms the role of a gene region on chromosome 5 in social immunity and simultaneously provides novel insights into genetic interactions between common mite resistance traits in honey bees.

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.

Evolution of a cytoplasmic determinant: evidence for the biochemical basis of functional evolution of the novel germ line regulator oskar

Germ line specification is essential in sexually reproducing organisms. Despite their critical role, the evolutionary history of the genes that specify animal germ cells is heterogeneous and dynamic. In many insects, the gene oskar is required for the specification of the germ line. However, the germ line role of oskar is thought to be a derived role resulting from co-option from an ancestral somatic role. To address how evolutionary changes in protein sequence could have led to changes in the function of Oskar protein that enabled it to regulate germ line specification, we searched for oskar orthologs in 1,565 publicly available insect genomic and transcriptomic data sets. The earliest-diverging lineage in which we identified an oskar ortholog was the order Zygentoma (silverfish and firebrats), suggesting that oskar originated before the origin of winged insects. We noted some order-specific trends in oskar sequence evolution, including whole gene duplications, clade-specific losses, and rapid divergence. An alignment of all known 379 Oskar sequences revealed new highly conserved residues as candidates that promote dimerization of the LOTUS domain. Moreover, we identified regions of the OSK domain with conserved predicted RNA binding potential. Furthermore, we show that despite a low overall amino acid conservation, the LOTUS domain shows higher conservation of predicted secondary structure than the OSK domain. Finally, we suggest new key amino acids in the LOTUS domain that may be involved in the previously reported Oskar−Vasa physical interaction that is required for its germ line role.

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.

Annotation of segmentation pathway genes in the Asian citrus psyllid, Diaphorina citri

Insects have a segmented body plan that is established during embryogenesis when the anterior-posterior (A-P) axis is divided into repeated units by a cascade of gene expression. The cascade is initiated by protein gradients created by translation of maternally provided mRNAs, localized at the anterior and posterior poles of the embryo. Combinations of these proteins activate specific gap genes to divide the embryo into distinct regions along the anterior-posterior axis. Gap genes then activate pair-rule genes, which are usually expressed in parts of every other segment. The pair-rule genes, in turn, activate expression of segment polarity genes in a portion of each segment. The segmentation genes are generally conserved among insects, although there is considerable variation in how they are deployed. We annotated 25 segmentation gene homologs in the Asian citrus psyllid, Diaphorina citri. Most of the genes expected to be present in D. citri based on their phylogenetic distribution in other insects were identified and annotated. Two exceptions were eagle and invected, which are present in at least some hemipterans, but were not found in D. citri. Many of the segmentation pathway genes are likely to be essential for D. citri development, and thus they may be useful targets for gene-based pest control methods.

The Power of Drosophila melanogaster for Modeling Neonicotinoid Effects on Pollinators and Identifying Novel Mechanisms

Neonicotinoids are the most widely used insecticides in the world and are implicated in the widespread population declines of insects including pollinators. Neonicotinoids target nicotinic acetylcholine receptors which are expressed throughout the insect central nervous system, causing a wide range of sub-lethal effects on non-target insects. Here, we review the potential of the fruit fly Drosophila melanogaster to model the sub-lethal effects of neonicotinoids on pollinators, by utilizing its well-established assays that allow rapid identification and mechanistic characterization of these effects. We compare studies on the effects of neonicotinoids on lethality, reproduction, locomotion, immunity, learning, circadian rhythms and sleep in D. melanogaster and a range of pollinators. We also highlight how the genetic tools available in D. melanogaster, such as GAL4/UAS targeted transgene expression system combined with RNAi lines to any gene in the genome including the different nicotinic acetylcholine receptor subunit genes, are set to elucidate the mechanisms that underlie the sub-lethal effects of these common pesticides. We argue that studying pollinators and D. melanogaster in tandem allows rapid elucidation of mechanisms of action, which translate well from D. melanogaster to pollinators. We focus on the recent identification of novel and important sublethal effects of neonicotinoids on circadian rhythms and sleep. The comparison of effects between D. melanogaster and pollinators and the use of genetic tools to identify mechanisms make a powerful partnership for the future discovery and testing of more specific insecticides.

High-Quality Genome Assemblies Reveal Long Non-coding RNAs Expressed in Ant Brains

Ants are an emerging model system for neuroepigenetics, as embryos with virtually identical genomes develop into different adult castes that display diverse physiology, morphology, and behavior. Although a number of ant genomes have been sequenced to date, their draft quality is an obstacle to sophisticated analyses of epigenetic gene regulation. We reassembled de novo high-quality genomes for two ant species, Camponotus floridanus and Harpegnathos saltator. Using long reads enabled us to span large repetitive regions and improve genome contiguity, leading to comprehensive and accurate protein-coding annotations that facilitated the identification of a Gp-9-like gene as differentially expressed in Harpegnathos castes. The new assemblies also enabled us to annotate long non-coding RNAs in ants, revealing caste-, brain-, and developmental-stage-specific long non-coding RNAs (lncRNAs) in Harpegnathos. These upgraded genomes, along with the new gene annotations, will aid future efforts to identify epigenetic mechanisms of phenotypic and behavioral plasticity in ants.

The torso-like gene functions to maintain the structure of the vitelline membrane in Nasonia vitripennis, implying its co-option into Drosophila axis formation

Axis specification is a fundamental developmental process. Despite this, the mechanisms by which it is controlled across insect taxa are strikingly different. An excellent example of this is terminal patterning, which in Diptera such as Drosophila melanogaster occurs via the localized activation of the receptor tyrosine kinase Torso. In Hymenoptera, however, the same process appears to be achieved via localized mRNA. How these mechanisms evolved and what they evolved from remains largely unexplored. Here, we show that torso-like, known for its role in Drosophila terminal patterning, is instead required for the integrity of the vitelline membrane in the hymenopteran wasp Nasonia vitripennis. We find that other genes known to be involved in Drosophila terminal patterning, such as torso and Ptth, also do not function in Nasonia embryonic development. These findings extended to orthologues of Drosophila vitelline membrane proteins known to play a role in localizing Torso-like in Drosophila; in Nasonia these are instead required for dorso–ventral patterning, gastrulation and potentially terminal patterning. Our data underscore the importance of the vitelline membrane in insect development, and implies phenotypes caused by knockdown of torso-like must be interpreted in light of its function in the vitelline membrane. In addition, our data imply that the signalling components of the Drosophila terminal patterning systems were co-opted from roles in regulating moulting, and co-option into terminal patterning involved the evolution of a novel interaction with the vitelline membrane protein Torso-like.

This article has an associated First Person interview with the first author of the paper.

Summary: In the parasitic wasp Nasonia, Tsl, a key component of the process that defines the termini of the embryo of Drosophila, has a function in the structure of the vitelline membrane.

A homeotic shift late in development drives mimetic color variation in a bumble bee

Natural phenotypic radiations, with their high diversity and convergence, are well-suited for informing how genomic changes translate to natural phenotypic variation. New genomic tools enable discovery in such traditionally nonmodel systems. Here, we characterize the genomic basis of color pattern variation in bumble bees (Hymenoptera, Apidae, Bombus), a group that has undergone extensive convergence of setal color patterns as a result of Müllerian mimicry. In western North America, multiple species converge on local mimicry patterns through parallel shifts of midabdominal segments from red to black. Using genome-wide association, we establish that a cis-regulatory locus between the abdominal fate-determining Hox genes, abd-A and Abd-B, controls the red-black color switch in a western species, Bombus melanopygus Gene expression analysis reveals distinct shifts in Abd-B aligned with the duration of setal pigmentation at the pupal-adult transition. This results in atypical anterior Abd-B expression, a late developmental homeotic shift. Changing expression of Hox genes can have widespread effects, given their important role across segmental phenotypes; however, the late timing reduces this pleiotropy, making Hox genes suitable targets. Analysis of this locus across mimics and relatives reveals that other species follow independent genetic routes to obtain the same phenotypes.

Evolution of the Torso activation cassette, a pathway required for terminal patterning and moulting

Embryonic terminal patterning and moulting are critical developmental processes in insects. In Drosophila and Tribolium both of these processes are regulated by the Torso-activation cassette (TAC). The TAC consists of a common receptor, Torso, ligands Trunk and prothoracicotropic hormone (PTTH), and the spatially restricted protein Torso-like, with combinations of these elements acting mechanistically to activate the receptor in different developmental contexts. In order to trace the evolutionary history of the TAC we determined the presence or absence of TAC components in the genomes of arthropods. Our analyses reveal that Torso, Trunk and PTTH are evolutionarily labile components of the TAC with multiple individual or combined losses occurring in the arthropod lineages leading to and within the insects. These losses are often correlated, with both ligands and receptor missing from the genome of the same species. We determine that the PTTH gene evolved in the common ancestor of Hemiptera and Holometabola, and is missing from the genomes of a number of species with experimentally demonstrated PTTH activity, implying another molecule may be involved in ecdysis in these species. In contrast, the torso-like gene is a common component of pancrustacean genomes.

A comparison of honeybee (Apis mellifera) queen, worker and drone larvae by RNA-Seq

Honeybees (Apis mellifera) have haplodiploid sex determination: males develop from unfertilized eggs and females develop from fertilized ones. The differences in larval food also determine the development of females. Here we compared the total somatic gene expression profiles of 2-day and 4-day-old drone, queen and worker larvae by RNA-Seq. The results from a co-expression network analysis on all expressed genes showed that 2-day-old drone and worker larvae were closer in gene expression profiles than 2-day-old queen larvae. This indicated that for young larvae (2-day-old) environmental factors such as larval diet have a greater effect on gene expression profiles than ploidy or sex determination. Drones had the most distinct gene expression profiles at the 4-day larval stage, suggesting that haploidy, or sex dramatically affects the gene expression of honeybee larvae. Drone larvae showed fewer differences in gene expression profiles at the 2-day and 4-day time points than the worker and queen larval comparisons (598 against 1190 and 1181), suggesting a different pattern of gene expression regulation during the larval development of haploid males compared to diploid females. This study indicates that early in development the queen caste has the most distinct gene expression profile, perhaps reflecting the very rapid growth and morphological specialization of this caste compared to workers and drones. Later in development the haploid male drones have the most distinct gene expression profile, perhaps reflecting the influence of ploidy or sex determination on gene expression.

lncRNA profile of Apis mellifera and its possible role in behavioural transition from nurses to foragers

BACKGROUND

The behavioural transition from nurses to foragers in honey bees is known to be affected by intrinsic and extrinsic factors, including colony demography, hormone levels, brain chemistry and structure, and gene expression in the brain. However, the molecular mechanism underlying this behavioural transition of honey bees is still obscure.

RESULTS

Through RNA sequencing, we performed a comprehensive analysis of lncRNAs and mRNAs in honey bee nurses and foragers. Nurses and foragers from both typical colonies and single-cohort colonies were used to prepare six libraries to generate 49 to 100 million clear reads per sample. We obtained 6863 novel lncRNAs, 1480 differentially expressed lncRNAs between nurses and foragers, and 9308 mRNAs. Consistent with previous studies, lncRNAs showed features distinct from mRNAs, such as shorter lengths, lower exon numbers, and lower expression levels compared to mRNAs. Bioinformatic analysis showed that differentially expressed genes were mostly involved in the regulation of sensory-related events, such as olfactory receptor activity and odorant binding, and enriched Wnt and FoxO signaling pathways. Moreover, we found that lncRNAs TCONS_00356023, TCONS_00357367, TCONS_00159909 and mRNAs dop1, Kr-h1 and HR38 may play important roles in behavioural transition in honey bees.

CONCLUSION

This study characterized the expression profile of lncRNAs in nurses and foragers and provided a framework for further study of the role of lncRNAs in honey bee behavioural transition.

A genetic switch for worker nutrition-mediated traits in honeybees

Highly social insects are characterized by caste dimorphism, with distinct size differences of reproductive organs between fertile queens and the more or less sterile workers. An abundance of nutrition or instruction via diet-specific compounds has been proposed as explanations for the nutrition-driven queen and worker polyphenism. Here, we further explored these models in the honeybee (Apis mellifera) using worker nutrition rearing and a novel mutational screening approach using the clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) method. The worker nutrition-driven size reduction of reproductive organs was restricted to the female sex, suggesting input from the sex determination pathway. Genetic screens on the sex determination genes in genetic females for size polyphenism revealed that doublesex (dsx) mutants display size-reduced reproductive organs irrespective of the sexual morphology of the organ tissue. In contrast, feminizer (fem) mutants lost the response to worker nutrition-driven size control. The first morphological worker mutants in honeybees demonstrate that the response to nutrition relies on a genetic program that is switched “ON” by the fem gene. Thus, the genetic instruction provided by the fem gene provides an entry point to genetically dissect the underlying processes that implement the size polyphenism.

In honeybees, nutrition drives dimorphic size development of reproductive organs in fertile queens and sterile workers. A study using the first induced morphological mutants in honeybees demonstrates that this developmental plasticity requires a genetic program that is switched on by the “feminizer” gene.

In honeybees, nutrition drives dimorphic size development of reproductive organs in fertile queens and sterile workers. The first induced morphological mutants in honeybees demonstrate that this developmental plasticity requires a genetic program that is switched “ON” by the feminizer (fem) gene.

Genome-Wide Identification and Expression Profiling of Wnt Family Genes in the Silkworm, Bombyx mori

Wnt is a family of conserved glycoproteins that participate in a variety of important biological processes including embryo development, cell proliferation and differentiation, and tissue regeneration. The Wnt family is a metazoan novelty found in all animal phyla. Studies have revealed that the number of Wnt genes varies among species, presumably due to reproduction and loss of genes during evolution. However, a comprehensive inventory of Wnt genes in Lepidoptera is lacking. In this study, we identified the repertoire of Wnt genes in the silkworm and seven other species of Lepidoptera and obtained eight Wnt genes (Wnt1, Wnt5–Wnt7, Wnt9–Wnt11, and WntA) in each species. Four of these Wnt genes are clustered in two orientations (5′-Wnt9-Wnt1-Wnt6-Wnt10-3′ and 5′-Wnt10-Wnt6-Wnt1-Wnt9-3′) in both moths and butterflies. Transcript analysis of Wnt in silkworm embryonic stages showed that each BmWnt gene had a unique expression pattern during embryological development. Analysis of a larval stage revealed differential expression of Wnt family members in diverse tissues. Our study provides an overview of the Wnt family in Lepidoptera and will inspire further functional study of the Wnt genes in the silkworm.

The trigger (and the restriction) of Torso RTK activation

The Torso pathway is an ideal model of receptor tyrosine kinase systems, in particular when addressing questions such as how receptor activity is turned on and, equally important, how it is restricted, how different outcomes can be generated from a single signal, and the extent to which gene regulation by signalling pathways relies on the relief of transcriptional repression. In this regard, we considered it pertinent to single out the fundamental notions learned from the Torso pathway beyond the specificities of this system (Furriols and Casanova 2003 EMBO J. 22, 1947-1952. ( doi:10.1093/emboj/cdg224 )). Since then, the Torso system has gained relevance and its implications beyond its original involvement in morphogenesis and into many disciplines such as oncogenesis, hormone control and neurobiology are now acknowledged. Thus, we believe that it is timely to highlight additional notions supported by new findings and to draw attention to future prospects. Given the late development of research in the field, we wish to devote this review to the events leading to the activation of the Torso receptor, the main focus of our most recent work.

Dynamic transcriptome landscape of Asian domestic honeybee (Apis cerana) embryonic development revealed by high-quality RNA sequencing

Background

Honeybee development consists of four stages: embryo, larva, pupa and adult. Embryogenesis, a key process of cell division and differentiation, takes 3 days in honeybees. However, the embryonic transcriptome and the dynamic regulation of embryonic transcription are still largely uncharacterized in honeybees, especially in the Asian honeybee (Apis cerana). Here, we employed high-quality RNA-seq to explore the transcriptome of Asian honeybee embryos at three ages, approximately 24, 48 and 72 h (referred to as Day1, Day2 and Day3, respectively).

Results

Nine embryo samples, three from each age, were collected for RNA-seq. According to the staging scheme of honeybee embryos and the morphological features we observed, our Day1, Day2 and Day3 embryos likely corresponded to the late stage four, stage eight and stage ten development stages, respectively. Hierarchical clustering and principal component analysis showed that same-age samples were grouped together, and the Day2 samples had a closer relationship with the Day3 samples than the Day1 samples. Finally, a total of 18,284 genes harboring 55,646 transcripts were detected in the A. cerana embryos, of which 44.5% consisted of the core transcriptome shared by all three ages of embryos. A total of 4088 upregulated and 3046 downregulated genes were identified among the three embryo ages, of which 2010, 3177 and 1528 genes were upregulated and 2088, 2294 and 303 genes were downregulated from Day1 to Day2, from Day1 to Day3 and from Day2 to Day3, respectively. The downregulated genes were mostly involved in cellular, biosynthetic and metabolic processes, gene expression and protein localization, and macromolecule modification; the upregulated genes mainly participated in cell development and differentiation, tissue, organ and system development, and morphogenesis. Interestingly, several biological processes related to the response to and detection of light stimuli were enriched in the first-day A. cerana embryogenesis but not in the Apis mellifera embryogenesis, which was valuable for further investigations.

Conclusions

Our transcriptomic data substantially expand the number of known transcribed elements in the A. cerana genome and provide a high-quality view of the transcriptome dynamics of A. cerana embryonic development.

Electronic supplementary material

The online version of this article (10.1186/s12861-018-0169-1) contains supplementary material, which is available to authorized users.

Analysis of the genome of the New Zealand giant collembolan (Holacanthella duospinosa) sheds light on hexapod evolution

BACKGROUND

The New Zealand collembolan genus Holacanthella contains the largest species of springtails (Collembola) in the world. Using Illumina technology we have sequenced and assembled a draft genome and transcriptome from Holacanthella duospinosa (Salmon). We have used this annotated assembly to investigate the genetic basis of a range of traits critical to the evolution of the Hexapoda, the phylogenetic position of H. duospinosa and potential horizontal gene transfer events.

RESULTS

Our genome assembly was ~375 Mbp in size with a scaffold N50 of ~230 Kbp and sequencing coverage of ~180×. DNA elements, LTRs and simple repeats and LINEs formed the largest components and SINEs were very rare. Phylogenomics (370,877 amino acids) placed H. duospinosa within the Neanuridae. We recovered orthologs of the conserved sex determination genes thought to play a role in sex determination. Analysis of CpG content suggested the absence of DNA methylation, and consistent with this we were unable to detect orthologs of the DNA methyltransferase enzymes. The small subunit rRNA gene contained a possible retrotransposon. The Hox gene complex was broken over two scaffolds. For chemosensory ability, at least 15 and 18 ionotropic glutamate and gustatory receptors were identified, respectively. However, we were unable to identify any odorant receptors or their obligate co-receptor Orco. Twenty-three chitinase-like genes were identified from the assembly. Members of this multigene family may play roles in the digestion of fungal cell walls, a common food source for these saproxylic organisms. We also detected 59 and 96 genes that blasted to bacteria and fungi, respectively, but were located on scaffolds that otherwise contained arthropod genes.

CONCLUSIONS

The genome of H. duospinosa contains some unusual features including a Hox complex broken over two scaffolds, in a different manner to other arthropod species, a lack of odorant receptor genes and an apparent lack of environmentally responsive DNA methylation, unlike many other arthropods. Our detection of candidate horizontal gene transfer candidates confirms that this phenomenon is occurring across Collembola. These findings allow us to narrow down the regions of the arthropod phylogeny where key innovations have occurred that have facilitated the evolutionary success of Hexapoda.

The honeybee as a model insect for developmental genetics

Honeybees are an important component of modern agricultural systems, and a fascinating and scientifically engrossing insect. Honeybees are not commonly used as model systems for understanding development in insects despite their importance in agriculture. Honeybee embryogenesis, while being superficially similar to Drosophila, is molecularly very different, especially in axis formation and sex determination. In later development, much of honeybee biology is modified by caste development, an as yet poorly understood, but excellent, system to study developmental plasticity. In adult stages, developmental plasticity of the ovaries, related to reproductive constraint exhibits another aspect of plasticity. Here they review the tools, current knowledge and opportunities in honeybee developmental biology, and provide an updated embryonic staging scheme to support future studies.

Molecular control of gut formation in the spider Parasteatoda tepidariorum

The development of a digestive system is an essential feature of bilaterians. Studies of the molecular control of gut formation in arthropods have been studied in detail in the fruit fly Drosophila melanogaster. However, little is known in other arthropods, especially in noninsect arthropods. To better understand the evolution of arthropod alimentary system, we investigate the molecular control of gut development in the spider Parasteatoda tepidariorum (Pt), the primary chelicerate model species for developmental studies. Orthologs of the ectodermal genes Pt-wingless (Pt-wg) and Pt-hedgehog (Pt-hh), of the endodermal genes, Pt-serpent (Pt-srp) and Pt-hepatocyte-nuclear factor-4 (Pt-hnf4) and of the mesodermal gene Pt-twist (Pt-twi) are expressed in the same germ layers during spider gut development as in D. melanogaster. Thus, our expression data suggest that the downstream molecular components involved in gut development in arthropods are conserved. However, Pt-forkhead (Pt-fkh) expression and function in spiders is considerably different from its D. melanogaster ortholog. Pt-fkh is expressed before gastrulation in a cell population that gives rise to endodermal and mesodermal precursors, suggesting a possible role for this factor in specification of both germ layers. To test this hypothesis, we knocked down Pt-fkh via RNA interference. Pt-fkh RNAi embryos not only fail to develop a proper gut, but also lack the mesodermal Pt-twi expressing cells. Thus, in spiders Pt-fkh specifies endodermal and mesodermal germ layers. We discuss the implications of these findings for the evolution and development of gut formation in Ecdysozoans.

MicroRNA-34 directly targets pair-rule genes and cytoskeleton component in the honey bee

MicroRNAs (miRNAs) are key regulators of developmental processes, such as cell fate determination and differentiation. Previous studies showed Dicer knockdown in honeybee embryos disrupt the processing of functional mature miRNAs and impairs embryo patterning. Here we investigated the expression profiles of miRNAs in honeybee embryogenesis and the role of the highly conserved miR-34-5p in the regulation of genes involved in insect segmentation. A total of 221 miRNAs were expressed in honey bee embryogenesis among which 97 mature miRNA sequences have not been observed before. Interestingly, we observed a switch in dominance between the 5-prime and 3-prime arm of some miRNAs in different embryonic stages; however, most miRNAs present one dominant arm across all stages of embryogenesis. Our genome-wide analysis of putative miRNA-target networks and functional pathways indicates miR-34-5p is one of the most conserved and connected miRNAs associated with the regulation of genes involved in embryonic patterning and development. In addition, we experimentally validated that miR-34-5p directly interacts to regulatory elements in the 3'-untranslated regions of pair-rule (even-skipped, hairy, fushi-tarazu transcription factor 1) and cytoskeleton (actin5C) genes. Our study suggests that miR-34-5p may regulate the expression of pair-rule and cytoskeleton genes during early development and control insect segmentation.

The evolution of insect germline specification strategies

The establishment of the germline is essential for sexually reproducing organisms. In animals, there are two major strategies to specify the germline: maternal provision and zygotic induction. The molecular basis of the maternal provision mode has been well characterized in several model organisms (fly, frog, fish, and nematode), while that of the zygotic induction mode has mainly been studied in mammalian models such as the mouse. Shifts in germline determination modes occur unexpectedly frequently and many such shifts have occurred several times among insects. Given their general tractability and rapidly increasing genomic and genetic tools applicable to many species, the insects present a uniquely powerful model system for understanding major transitions in reproductive strategies, and developmental processes in general.

Transcriptome Analysis of Honeybee (Apis Mellifera) Haploid and Diploid Embryos Reveals Early Zygotic Transcription during Cleavage

In honeybees, the haplodiploid sex determination system promotes a unique embryogenesis process wherein females develop from fertilized eggs and males develop from unfertilized eggs. However, the developmental strategies of honeybees during early embryogenesis are virtually unknown. Similar to most animals, the honeybee oocytes are supplied with proteins and regulatory elements that support early embryogenesis. As the embryo develops, the zygotic genome is activated and zygotic products gradually replace the preloaded maternal material. The analysis of small RNA and mRNA libraries of mature oocytes and embryos originated from fertilized and unfertilized eggs has allowed us to explore the gene expression dynamics in the first steps of development and during the maternal-to-zygotic transition (MZT). We localized a short sequence motif identified as TAGteam motif and hypothesized to play a similar role in honeybees as in fruit flies, which includes the timing of early zygotic expression (MZT), a function sustained by the presence of the zelda ortholog, which is the main regulator of genome activation. Predicted microRNA (miRNA)-target interactions indicated that there were specific regulators of haploid and diploid embryonic development and an overlap of maternal and zygotic gene expression during the early steps of embryogenesis. Although a number of functions are highly conserved during the early steps of honeybee embryogenesis, the results showed that zygotic genome activation occurs earlier in honeybees than in Drosophila based on the presence of three primary miRNAs (pri-miRNAs) (ame-mir-375, ame-mir-34 and ame-mir-263b) during the cleavage stage in haploid and diploid embryonic development.

Comparative analysis of zygotic developmental genes in Rhodnius prolixus genome shows conserved features on the tracheal developmental pathway

Most of the in-depth studies on insect developmental genetic have been carried out in the fruit fly Drosophila melanogaster, an holometabolous insect, so much more still remains to be studied in hemimetabolous insects. Having Rhodnius prolixus sequenced genome available, we search for orthologue genes of zygotic signaling pathways, segmentation, and tracheogenesis in the R. prolixus genome and in three species of Triatoma genus transcriptomes, concluding that there is a high level of gene conservation. We also study the function of two genes required for tracheal system development in D. melanogaster - R. prolixus orthologues: trachealess (Rp-trh) and empty spiracles (Rp-ems). From that we see that Rp-trh is required for early tracheal development since Rp-trh RNAi shows that the primary tracheal branches fail to form. On the other hand, Rp-ems is implied in the proper formation of the posterior tracheal branches, in a similar way to D. melanogaster. These results represent the initial characterization of the genes involved in the tracheal development of an hemimetabolous insect building a bridge between the current genomic era and V. Wigglesworth's classical studies on insects' respiratory system physiology.

The genomes of two key bumblebee species with primitive eusocial organization

Background

The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats.

Results

We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits.

Conclusions

These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0623-3) contains supplementary material, which is available to authorized users.

Transcriptomic study of the red palm weevil Rhynchophorus ferrugineus embryogenesis

The red palm weevil (RPW), Rhynchophorus ferrugineus (Coleoptera: Curculionidae), is an invasive, concealed and destructive tissue borer, and it becomes a lethal pest of the palm family of plants and has been reported to attack 20 palm species around the globe. Here we report a systematic transcriptomic study on embryogenesis of RPW, where we analyze the transcriptomes across five developmental stages of RPW embryogenesis, involving four embryonic stages (E1, E2, E3 and E4) and one larval stage (L1). Using the RNA-seq and next-generation platforms, we generated 80 to 91 million reads for each library and assemble 22 532 genes that are expressed at different embryonic stages. Among the total transcripts from the five embryonic development stages, we found that 30.45 % are differentially expressed, 10.10 % show stage-specificity and even a larger fraction, 62.88 %, exhibit constitutive expression in all the stages. We also analyzes the expression dynamics of several conserved signaling pathways (such as Hedgehog, JAK-STAT, Notch, TGF-β, Ras/MAPK and Wnt), as well as key developmental genes, including those related to apoptosis, axis formation, Hox complex, neurogenesis and segmentation. The datasets provide an essential resource for gene annotation and RPW functional genomics, including studies by using tools and concepts from multiple disciplines, such as development, physiology, biochemistry, molecular biology and genetics.

Evidence for the plasticity of arthropod signal transduction pathways

Metazoans are known to contain a limited, yet highly conserved, set of signal transduction pathways that instruct early developmental patterning mechanisms. Genomic surveys that have compared gene conservation in signal transduction pathways between various insects and Drosophila support the conclusion that these pathways are conserved in evolution. However, the degree to which individual components of signal transduction pathways vary among more divergent arthropods is not known. Here, we report our results of a survey of the genome of the two-spotted spider mite Tetranychus urticae, using a set of 294 Drosophila orthologs of genes that function in signal transduction. We find a third of all genes surveyed absent from the spider mite genome. We also identify several novel duplications that have not been previously reported for a chelicerate. In comparison with previous insect surveys, Tetranychus contains a decrease in overall gene conservation, as well as an unusual ratio of ligands to receptors and other modifiers. These findings suggest that gene loss and duplication among components of signal transduction pathways are common among arthropods and suggest that signal transduction pathways in arthropods are more evolutionarily labile than previously hypothesized.

Capturing embryonic development from metamorphosis: how did the terminal patterning signalling pathway of Drosophila evolve?

The Torso receptor tyrosine kinase has two crucial roles in Drosophila melanogaster development. One is in the control of insect moulting, which is regulated by the neuropeptide hormone PTTH (prothoracicotropic hormone). PTTH activates ERK signalling via Torso in the prothoracic gland to stimulate ecdysone secretion. Torso also has a role in control of one of the earliest events in embryogenesis in Drosophila; patterning of the embryonic termini. Here Torso is activated by a different, but related, peptide called Trunk. During terminal patterning another protein, Torso-like, has a key role in mediating activation of Torso by Trunk. Torso-like is also expressed in the prothoracic gland and null-mutants have defective developmental timing in Drosophila. This function, however, has been recently shown to be independent of Torso and PTTH. We refer to these proteins, Trunk, PTTH, Torso and Torso-like, as the Torso-activation module. Outside Drosophila we see that the genes encoding the Torso-activation module have a complex phylogenetic history, with different origins and multiple losses of components of this signalling pathway during arthropod evolution. This, together with expression and functional data in a range of insects, leads us to propose that the terminal patterning pathway in Drosophila and Tribolium arose through co-option of PTTH/Trunk and Torso, which has a role in developmental timing, into a new context, and that Torso-like was recruited specifically in the ovary to modulate the specificity of this pathway.

In-depth proteomics characterization of embryogenesis of the honey bee worker (Apis mellifera ligustica)

Identifying proteome changes of honey bee embryogenesis is of prime importance for unraveling the molecular mechanisms that they underlie. However, many proteomic changes during the embryonic period are not well characterized. We analyzed the proteomic alterations over the complete time course of honey bee worker embryogenesis at 24, 48, and 72 h of age, using mass spectrometry-based proteomics, label-free quantitation, and bioinformatics. Of the 1460 proteins identified the embryo of all three ages, the core proteome (proteins shared by the embryos of all three ages, accounting for 40%) was mainly involved in protein synthesis, metabolic energy, development, and molecular transporter, which indicates their centrality in driving embryogenesis. However, embryos at different developmental stages have their own specific proteome and pathway signatures to coordinate and modulate developmental events. The young embryos (<24 h) stronger expression of proteins related to nutrition storage and nucleic acid metabolism may correlate with the cell proliferation occurring at this stage. The middle aged embryos (24-48 h) enhanced expression of proteins associated with cell cycle control, transporters, antioxidant activity, and the cytoskeleton suggest their roles to support rudimentary organogenesis. Among these proteins, the biological pathways of aminoacyl-tRNA biosynthesis, β-alanine metabolism, and protein export are intensively activated in the embryos of middle age. The old embryos (48-72 h) elevated expression of proteins implicated in fatty acid metabolism and morphogenesis indicate their functionality for the formation and development of organs and dorsal closure, in which the biological pathways of fatty acid metabolism and RNA transport are highly activated. These findings add novel understanding to the molecular details of honey bee embryogenesis, in which the programmed activation of the proteome matches with the physiological transition observed during embryogenesis. The identified biological pathways and key node proteins allow for further functional analysis and genetic manipulation for both the honey bee embryos and other eusocial insects.

Components of the dorsal-ventral pathway also contribute to anterior-posterior patterning in honeybee embryos (Apis mellifera)

Background

A key early step in embryogenesis is the establishment of the major body axes; the dorsal-ventral (DV) and anterior-posterior (AP) axes. Determination of these axes in some insects requires the function of different sets of signalling pathways for each axis. Patterning across the DV axis requires interaction between the Toll and Dpp/TGF-β pathways, whereas patterning across the AP axis requires gradients of bicoid/orthodenticle proteins and the actions of a hierarchy of gene transcription factors. We examined the expression and function of Toll and Dpp signalling during honeybee embryogenesis to assess to the role of these genes in DV patterning.

Results

Pathway components that are required for dorsal specification in Drosophila are expressed in an AP-restricted pattern in the honeybee embryo, including Dpp and its receptor Tkv. Components of the Toll pathway are expressed in a more conserved pattern along the ventral axis of the embryo. Late-stage embryos from RNA interference (RNAi) knockdown of Toll and Dpp pathways had both DV and AP patterning defects, confirmed by staining with Am-sna, Am-zen, Am-eve, and Am-twi at earlier stages. We also identified two orthologues of dorsal in the honeybee genome, with one being expressed during embryogenesis and having a minor role in axis patterning, as determined by RNAi and the other expressed during oogenesis.

Conclusions

We found that early acting pathways (Toll and Dpp) are involved not only in DV patterning but also AP patterning in honeybee embryogenesis. Changes to the expression patterns and function of these genes may reflect evolutionary changes in the placement of the extra-embryonic membranes during embryogenesis with respect to the AP and DV axes.

Dual mode of embryonic development is highlighted by expression and function of Nasonia pair-rule genes

Embryonic anterior-posterior patterning is well understood in Drosophila, which uses 'long germ' embryogenesis, in which all segments are patterned before cellularization. In contrast, most insects use 'short germ' embryogenesis, wherein only head and thorax are patterned in a syncytial environment while the remainder of the embryo is generated after cellularization. We use the wasp Nasonia (Nv) to address how the transition from short to long germ embryogenesis occurred. Maternal and gap gene expression in Nasonia suggest long germ embryogenesis. However, the Nasonia pair-rule genes even-skipped, odd-skipped, runt and hairy are all expressed as early blastoderm pair-rule stripes and late-forming posterior stripes. Knockdown of Nv eve, odd or h causes loss of alternate segments at the anterior and complete loss of abdominal segments. We propose that Nasonia uses a mixed mode of segmentation wherein pair-rule genes pattern the embryo in a manner resembling Drosophila at the anterior and ancestral Tribolium at the posterior. DOI: http://dx.doi.org/10.7554/eLife.01440.001.

The evolving role of the orphan nuclear receptor ftz-f1, a pair-rule segmentation gene

Segmentation is a critical developmental process that occurs by different mechanisms in diverse taxa. In insects, there are three common modes of embryogenesis-short-, intermediate-, and long-germ development-which differ in the number of segments specified at the blastoderm stage. While genes involved in segmentation have been extensively studied in the long-germ insect Drosophila melanogaster (Dm), it has been found that their expression and function in segmentation in short- and intermediate-germ insects often differ. Drosophila ftz-f1 encodes an orphan nuclear receptor that functions as a maternally expressed pair-rule segmentation gene, responsible for the formation of alternate body segments during Drosophila embryogenesis. Here we investigated the expression and function of ftz-f1 in the short-germ beetle, Tribolium castaneum (Tc). We found that Tc-ftz-f1 is expressed in stripes in Tribolium embryos. These stripes overlap alternate Tc-Engrailed (Tc-En) stripes, indicative of a pair-rule expression pattern. To test whether Tc-ftz-f1 has pair-rule function, we utilized embryonic RNAi, injecting double-stranded RNA corresponding to Tc-ftz-f1 coding or non-coding regions into early Tribolium embryos. Knockdown of Tc-ftz-f1 produced pair-rule segmentation defects, evidenced by loss of expression of alternate En stripes. In addition, a later role for Tc-ftz-f1 in cuticle formation was revealed. These results identify a new pair-rule gene in Tribolium and suggest that its role in segmentation may be shared among holometabolous insects. Interestingly, while Tc-ftz-f1 is expressed in pair-rule stripes, the gene is ubiquitously expressed in Drosophila embryos. Thus, the pair-rule function of ftz-f1 is conserved despite differences in expression patterns of ftz-f1 genes in different lineages. This suggests that ftz-f1 expression changed after the divergence of lineages leading to extant beetles and flies, likely due to differences in cis-regulatory sequences. We propose that the dependence of Dm-Ftz-F1 on interaction with the homeodomain protein Ftz which is expressed in stripes in Drosophila, loosened constraints on Dm-ftz-f1 expression, allowing for ubiquitous expression of this pair-rule gene in Drosophila.

To "bee or not to bee" male or female? An educational primer for use with "the Am-tra2 gene is an essential regulator of female splice regulation at two levels of the sex determination hierarchy of the honeybee"

An article by Nissen et al. in the November 2012 issue of GENETICS emphasizes the importance of alternative splicing in the sex determination cascade of the honeybee Apis mellifera. This study demonstrates the application of reverse transcriptase PCR and RNA interference screens as genetic tools to better understand the regulation of transcription and splicing. It also provides the opportunity to explore the evolutionary origins of genes by considering the functions of orthologs and paralogs in different species. This Primer article provides background information and explanations of the concepts and findings of Nissen et al. and offers discussion questions for use in the classroom.

The pea aphid (Acyrthosiphon pisum) genome encodes two divergent early developmental programs

The pea aphid (Acyrthosiphon pisum) can reproduce either sexually or asexually (parthenogenetically), giving rise, in each case, to almost identical adults. These two modes of reproduction are accompanied by differences in ovarian morphology and the developmental environment of the offspring, with sexual forms producing eggs that are laid, whereas asexual development occurs within the mother. Here we examine the effect each mode of reproduction has on the expression of key maternal and axis patterning genes; orthodenticle (otd), hunchback (hb), caudal (cad) and nanos (nos). We show that three of these genes (Ap-hb, Ap-otd and Ap-cad) are expressed differently between the sexually and asexually produced oocytes and embryos of the pea aphid. We also show, using immunohistochemistry and cytoskeletal inhibitors, that Ap-hb RNA is localized differently between sexually and asexually produced oocytes, and that this is likely due to differences in the 3' untranslated regions of the RNA. Furthermore, Ap-hb and Ap-otd have extensive expression domains in early sexually produced embryos, but are not expressed at equivalent stages in asexually produced embryos. These differences in expression likely correspond with substantial changes in the gene regulatory networks controlling early development in the pea aphid. These data imply that in the evolution of parthenogenesis a new program has evolved to control the development of asexually produced embryos, whilst retaining the existing, sexual, developmental program. The patterns of modification of these developmental processes mirror the changes that we see in developmental processes between species, in that early acting pathways in development are less constrained, and evolve faster, than later ones. We suggest that the evolution of the novel asexual development pathway in aphids is not a simple modification of an ancestral system, but the evolution of two very different developmental mechanisms occurring within a single species.

Unscrambling butterfly oogenesis

Background

Butterflies are popular model organisms to study physiological mechanisms underlying variability in oogenesis and egg provisioning in response to environmental conditions. Nothing is known, however, about; the developmental mechanisms governing butterfly oogenesis, how polarity in the oocyte is established, or which particular maternal effect genes regulate early embryogenesis. To gain insights into these developmental mechanisms and to identify the conserved and divergent aspects of butterfly oogenesis, we analysed a de novo ovarian transcriptome of the Speckled Wood butterfly Pararge aegeria (L.), and compared the results with known model organisms such as Drosophila melanogaster and Bombyx mori.

Results

A total of 17306 contigs were annotated, with 30% possibly novel or highly divergent sequences observed. Pararge aegeria females expressed 74.5% of the genes that are known to be essential for D. melanogaster oogenesis. We discuss the genes involved in all aspects of oogenesis, including vitellogenesis and choriogenesis, plus those implicated in hormonal control of oogenesis and transgenerational hormonal effects in great detail. Compared to other insects, a number of significant differences were observed in; the genes involved in stem cell maintenance and differentiation in the germarium, establishment of oocyte polarity, and in several aspects of maternal regulation of zygotic development.

Conclusions

This study provides valuable resources to investigate a number of divergent aspects of butterfly oogenesis requiring further research. In order to fully unscramble butterfly oogenesis, we also now also have the resources to investigate expression patterns of oogenesis genes under a range of environmental conditions, and to establish their function.

The pea aphid uses a version of the terminal system during oviparous, but not viviparous, development

Background

In most species of aphid, female nymphs develop into either sexual or asexual adults depending on the length of the photoperiod to which their mothers were exposed. The progeny of these sexual and asexual females, in turn, develop in dramatically different ways. The fertilized oocytes of sexual females begin embryogenesis after being deposited on leaves (oviparous development) while the oocytes of asexual females complete embryogenesis within the mother (viviparous development). Compared with oviparous development, viviparous development involves a smaller transient oocyte surrounded by fewer somatic epithelial cells and a smaller early embryo that comprises fewer cells. To investigate whether patterning mechanisms differ between the earliest stages of the oviparous and viviparous modes of pea aphid development, we examined the expression of pea aphid orthologs of genes known to specify embryonic termini in other insects.

Results

Here we show that pea aphid oviparous ovaries express torso-like in somatic posterior follicle cells and activate ERK MAP kinase at the posterior of the oocyte. In addition to suggesting that some posterior features of the terminal system are evolutionarily conserved, our detection of activated ERK in the oocyte, rather than in the embryo, suggests that pea aphids may transduce the terminal signal using a mechanism distinct from the one used in Drosophila. In contrast with oviparous development, the pea aphid version of the terminal system does not appear to be used during viviparous development, since we did not detect expression of torso-like in the somatic epithelial cells that surround either the oocyte or the blastoderm embryo and we did not observe restricted activated ERK in the oocyte.

Conclusions

We suggest that while oviparous oocytes and embryos may specify posterior fate through an aphid terminal system, viviparous oocytes and embryos employ a different mechanism, perhaps one that does not rely on an interaction between the oocyte and surrounding somatic cells. Together, these observations provide a striking example of a difference in the fundamental events of early development that is both environmentally induced and encoded by the same genome.

Canonical terminal patterning is an evolutionary novelty

Patterning of the terminal regions of the Drosophila embryo is achieved by an exquisitely regulated signal that passes between the follicle cells of the ovary, and the developing embryo. This pathway, however, is missing or modified in other insects. Here we trace the evolution of this pathway by examining the origins and expression of its components. The three core components of this pathway: trunk, torso and torso-like have different evolutionary histories and have been assembled step-wise to form the canonical terminal patterning pathway of Drosophila and Tribolium. Trunk, torso and a gene unrelated to terminal patterning, prothoraciotrophic hormone (PTTH), show an intimately linked evolutionary history, with every holometabolous insect, except the honeybee, possessing both PTTH and torso genes. Trunk is more restricted in its phylogenetic distribution, present only in the Diptera and Tribolium and, surprisingly, in the chelicerate Ixodes scapularis, raising the possibility that trunk and torso evolved earlier than previously thought. In Drosophila torso-like restricts the activation of the terminal patterning pathway to the poles of the embryo. Torso-like evolved in the pan-crustacean lineage, but based on expression of components of the canonical terminal patterning system in the hemimetabolous insect Acyrthosiphon pisum and the holometabolous insect Apis mellifera, we find that the canonical terminal-patterning system is not active in these insects. We therefore propose that the ancestral function of torso-like is unrelated to terminal patterning and that torso-like has become co-opted into terminal patterning in the lineage leading to Coleoptera and Diptera. We also show that this co-option has not resulted in changes to the molecular function of this protein. Torso-like from the pea aphid, honeybee and Drosophila, despite being expressed in different patterns, are functionally equivalent. We propose that co-option of torso-like into restricting the activity of trunk and torso facilitated the final step in the evolution of this pathway; the capture of transcriptional control of target genes such as tailless and huckebein by this complex and novel patterning pathway.

Variation and constraint in Hox gene evolution

Despite enormous body plan variation, genes regulating embryonic development are highly conserved. Here, we probe the mechanisms that predispose ancient regulatory genes to reutilization and diversification rather than evolutionary loss. The Hox gene fushi tarazu (ftz) arose as a homeotic gene but functions as a pair-rule segmentation gene in Drosophila. ftz shows extensive variation in expression and protein coding regions but has managed to elude loss from arthropod genomes. We asked what properties prevent this loss by testing the importance of different protein motifs and partners in the developing CNS, where ftz expression is conserved. Drosophila Ftz proteins with mutated protein motifs were expressed under the control of a neurogenic-specific ftz cis-regulatory element (CRE) in a ftz mutant background rescued for segmentation defects. Ftz CNS function did not require the variable motifs that mediate differential cofactor interactions involved in homeosis or segmentation, which vary in arthropods. Rather, CNS function did require the shared DNA-binding homeodomain, which plays less of a role in Ftz segmentation activity. The Antennapedia homeodomain substituted for Ftz homeodomain function in the Drosophila CNS, but full-length Antennapedia did not rescue CNS defects. These results suggest that a core CNS function retains ftz in arthropod genomes. Acquisition of a neurogenic CRE led to ftz expression in unique CNS cells, differentiating its role from neighboring Hox genes, rendering it nonredundant. The inherent flexibility of modular CREs and protein domains allows for stepwise acquisition of new functions, explaining broad retention of regulatory genes during animal evolution.

RNA localization in the honeybee (Apis mellifera) oocyte reveals insights about the evolution of RNA localization mechanisms

Subcellular localization of RNAs is a critical biological process for generation of cellular asymmetries for many cell types and a critical step in axis determination during the early development of animals. We have identified transcripts localized to the anterior and posterior of honeybee oocyte using laser capture microscopy and microarray analysis. Analysis of orthologous transcripts in Drosophila indicates that many do not show a conserved pattern of localization. By microinjecting fluorescently labeled honeybee transcripts into Drosophila egg chambers we show that these RNAs become localized in a similar manner to their localization in honeybee oocytes, indicating conservation of the localization machinery. Thus while the mechanisms for localizing RNA are conserved, the complement of localized RNAs are not. We propose that this complement of localized RNAs may change relatively rapidly through the loss or evolution of signal sequences detected by the conserved localization machinery, and show this has occurred in one transcript that is localized in a novel way in the honeybee. Our proposal, that the acquisition of novel RNA localization is relatively easy to evolve, has implications for the evolution of symmetry breaking mechanisms that trigger axis formation and development in animal embryos.