Differential expression of mitochondrial genes between queens and workers during caste determination in the honeybee Apis mellifera

Acetylcholine and choline in honey bee (Apis mellifera) worker brood food are seasonal and age-dependent

Nursing honeybees produce brood food with millimolar concentrations of acetylcholine (ACh), which is synthesized through head gland secretions mixed with honey stomach contents. While we previously demonstrated the necessity of ACh for proper larval development, the dynamics of ACh levels throughout ontogenesis and their seasonal variations have remained unclear until now. Our HPLC analysis reveals dependencies of choline and ACh levels on larval development days (LDDs), influenced by seasonal (April-September) variations. Median ACh concentrations peak on LDD 2, declining significantly toward cell capping, while choline levels are lowest during the initial LDDs, rising markedly toward cell capping. Seasonal patterns show peak ACh levels from April to June and a low in August, paralleling choline's peak in July and low in August. This seasonality holds consistently across multiple years (2020-2022) and colonies, despite potential variations in colony performance and environmental conditions. Our analysis found no correlation between temperature, sunshine, precipitation, or favourable foraging days and ACh/choline levels, suggesting the involvement of additional factors. These findings underscore the seasonal fluctuation of ACh levels and its potential implications for the genetic programs governing winter bee development.

Genetic Drift and Purifying Selection Shaped Mitochondrial Genome Variation in the High Royal Jelly-Producing Honeybee Strain (Apis mellifera ligustica)

Mitochondrial genomes (mitogenomes) are involved in cellular energy metabolism and have been shown to undergo adaptive evolution in organisms with increased energy-consuming activities. The genetically selected high royal jelly-producing bees (RJBs, Apis mellifera ligustica) in China can produce 10 times more royal jelly, a highly nutritional and functional food, relative to unselected Italian bees (ITBs). To test for potential adaptive evolution of RJB mitochondrial genes, we sequenced mitogenomes from 100 RJBs and 30 ITBs. Haplotype network and phylogenetic analysis indicate that RJBs and ITBs are not reciprocally monophyletic but mainly divided into the RJB- and ITB-dominant sublineages. The RJB-dominant sublineage proportion is 6-fold higher in RJBs (84/100) than in ITBs (4/30), which is mainly attributable to genetic drift rather than positive selection. The RJB-dominant sublineage exhibits a low genetic diversity due to purifying selection. Moreover, mitogenome abundance is not significantly different between RJBs and ITBs, thereby rejecting the association between mitogenome copy number and royal jelly-producing performance. Our findings demonstrate low genetic diversity levels of RJB mitogenomes and reveal genetic drift and purifying selection as potential forces driving RJB mitogenome evolution.

Transcriptomic, Morphological, and Developmental Comparison of Adult Honey Bee Queens (Apis mellifera) Reared From Eggs or Worker Larvae of Differing Ages

Queens and workers are very distinct phenotypes that develop from the same genome. Larvae from worker cells up to 3.5 d old can be transferred to larger queen cells and will subsequently be reared as queens and develop into functional queens. This has become a very popular queen rearing practice in contemporary apiculture. Here we used RNA-Seq to study the consequences of rearing queens from transplanted worker larvae on the transcriptome of the adult queens. We found that queens reared from transferred older larvae developed slower, weighted less, and had fewer ovarioles than queens reared from transferred eggs, indicating queens were cryptically intercaste. RNA-Seq analysis revealed differentially expressed genes between queens reared from transferred larvae compared with queens reared from transferred eggs: the older the larvae transferred, the greater the number of differentially expressed genes. Many of the differentially expressed genes had functions related to reproduction, longevity, immunity, or metabolism, suggesting that the health and long-term viability of queens was compromised. Our finds verify the previous studies that adult queens reared from older transferred larvae were of lower quality than queens reared from transferred eggs or younger larvae.

Effects of larval Age at Grafting and Juvenile Hormone on Morphometry and Reproductive Quality Parameters of in Vitro Reared Honey Bees (Hymenoptera: Apidae)

The honey bee queen plays a central role in the Apis mellifera L. (Hymenoptera: Apidae) colony, and her high reproductive capacity is fundamental for building up the workforce of a colony. Caste development in honey bee females involves elaborate physiological pathways unleashed at the beginning of the first larval instars, with juvenile hormone (JH) playing a crucial role. Here we took advantage of established in vitro rearing techniques to conduct a 2 × 2 experimental design and test initial rearing age (young vs old) and JH treatment (JH III vs solvent control) to enlighten the role of nutrient quality and JH in shaping honey bee female fertility, morphological features related to queenliness, and key physiological parameters (hemolymph vitellogenin/Vg, sugar levels, and Vg transcript levels). Our results show that while the age at initial larval rearing had major impacts on external morphology development, where younger larvae exhibited a higher probability to develop into queen-like adults morphotypes, the JH application during the larval stage improved physiological pathways related to ovary development and metabolism during the ontogenic development. We detected that the supplementation of queen larvae with JH promoted important benefits regarding queen fertility as the increase of ovariole number and vg levels at hemolymph, both crucial factors at eggs production. The data presented here provide guidance in efforts to improve honey bee queen quality, especially in light of frequent episodes of queen failures in the beekeeping industry.

The Role of Brood in Eusocial Hymenoptera

Study of social traits in offspring traditionally reflects on interactions in simple family groups, with famous examples including parent-offspring conflict and sibling rivalry in birds and mammals. In contrast, studies of complex social groups such as the societies of ants, bees, and wasps focus mainly on adults and, in particular, on traits and interests of queens and workers. The social role of developing individuals in complex societies remains poorly understood. We attempt to fill this gap by illustrating that development in social Hymenoptera constitutes a crucial life stage with important consequences for the individual as well as the colony. We begin by describing the complex social regulatory network that modulates development in Hymenoptera societies. By highlighting the inclusive fitness interests of developing individuals, we show that they may differ from those of other colony members. We then demonstrate that offspring have evolved specialized traits that allow them to play a functional, cooperative role within colonies and give them the potential power to act toward increasing their inclusive fitness. We conclude by providing testable predictions for investigating the role of brood in colony interactions and giving a general outlook on what can be learned from studying offspring traits in hymenopteran societies.

Abscisic acid enhances cold tolerance in honeybee larvae

The natural composition of nutrients present in food is a key factor determining the immune function and stress responses in the honeybee (Apis mellifera). We previously demonstrated that a supplement of abscisic acid (ABA), a natural component of nectar, pollen, and honey, increases honeybee colony survival overwinter. Here we further explored the role of ABA in in vitro-reared larvae exposed to low temperatures. Four-day-old larvae (L4) exposed to 25°C for 3 days showed lower survival rates and delayed development compared to individuals growing at a standard temperature (34°C). Cold-stressed larvae maintained higher levels of ABA for longer than do larvae reared at 34°C, suggesting a biological significance for ABA. Larvae fed with an ABA-supplemented diet completely prevent the low survival rate due to cold stress and accelerate adult emergence. ABA modulates the expression of genes involved in metabolic adjustments and stress responses: Hexamerin 70b, Insulin Receptor Substrate, Vitellogenin, and Heat Shock Proteins 70. AmLANCL2, the honeybee ABA receptor, is also regulated by cold stress and ABA. These results support a role for ABA increasing the tolerance of honeybee larvae to low temperatures through priming effects.

Identification of Differentially Expressed Genes In Deltamethrin-Resistant Culex pipiens quinquefasciatus

Culex quinquefasciatus is one of China's major house-dwelling mosquito species and an important vector of filariasis and encephalitis. Chemical treatments represent one of the most successful approaches for comprehensive mosquito prevention and control. However, the widespread use of chemical pesticides has led to the occurrence and development of insecticide resistance. Therefore, in-depth studies of resistance to insecticides are of vital importance. In this study, we performed a gene expression analysis to investigate genes from Cx. quinquefasciatus that may confer pyrethroid resistance. We aimed to understand the mechanisms of Cx. quinquefasciatus resistance to pyrethroid insecticides and provide insights into insect resistance management. Using a resistance bioassay, we determined the deltamethrin LC₅₀ values (lethal concentration required to kill 50% of the population) for Cx. quinquefasciatus larvae in the F₂₁, F₂₃, F₂₄, F₂₆, F₂₇, and F₃₀ generations. The 7 tested strains exhibited pesticide resistance that was 25.25 to 87.83 times higher than that of the SanYa strain. Moreover, the expression of the OBPjj7a (odorant-binding protein OBPjj7a), OBP28 (odorant-binding protein OBP28), and E2 (ubiquitin-conjugating enzyme) genes was positively correlated with deltamethrin resistance ( R² = 0.836, P = 0.011; R² = 0.788, P = 0.018; and R² = 0.850, P = 0.009, respectively) in Cx. quinquefasciatus. The expression of 4 additional genes, H/ACA, S19, SAR2, and PGRP, was not correlated with deltamethrin resistance. In summary, this study identified 3 Cx. quinquefasciatus genes with potential involvement in deltamethrin resistance, and these results may provide a theoretical basis for the control of mosquito resistance and insights into resistance detection.

Ties between ageing plasticity and reproductive physiology in honey bees (Apis mellifera) reveal a positive relation between fecundity and longevity as consequence of advanced social evolution

Honey bees (Apis mellifera) are the best studied model of ageing among the social insects. As in other social insects, the reproductive queen far outlives her non-reproductive workers despite developing from the same genome in the same colony environment. Thus, the different social roles of the two female castes are critical for the profound phenotypic plasticity. In several special cases, such as the reproductive workers of Apis mellifera capensis, within-caste plasticity enables further studies of the fecundity-longevity syndrome in honey bees. At present, molecular evidence suggests that a reorganization of physiological control pathways may facilitate longevity of reproductive individuals. However, the social role and social environment of the different colony members are also very important and one of the key future questions is how much social facilitation versus internal regulation is responsible for the positive association between fecundity and longevity in honey bees.

Mitochondrial structure and dynamics as critical factors in honey bee (Apis mellifera L.) caste development

The relationship between nutrition and phenotype is an especially challenging question in cases of facultative polyphenism, like the castes of social insects. In the honey bee, Apis mellifera, unexpected modifications in conserved signaling pathways revealed the hypoxia response as a possible mechanism underlying the regulation of body size and organ growth. Hence, the current study was designed to investigate possible causes of why the three hypoxia core genes are overexpressed in worker larvae. Parting from the hypothesis that this has an endogenous cause and is not due to differences in external oxygen levels we investigated mitochondrial numbers and distribution, as well as mitochondrial oxygen consumption rates in fat body cells of queen and worker larvae during the caste fate-critical larval stages. By immunofluorescence and electron microscopy we found higher densities of mitochondria in queen larval fat body, a finding further confirmed by a citrate synthase assay quantifying mitochondrial functional units. Oxygen consumption measurements by high-resolution respirometry revealed that queen larvae have higher maximum capacities of ATP production at lower physiological demand. Finally, the expression analysis of mitogenesis-related factors showed that the honey bee TFB1 and TFB2 homologs, and a nutritional regulator, ERR, are overexpressed in queen larvae. These results are strong evidence that the differential nutrition of queen and worker larvae by nurse bees affects mitochondrial dynamics and functionality in the fat body of these larvae, hence explaining their differential hypoxia response.

Molecular mechanisms of phenotypic plasticity in social insects

Polyphenism in insects, whereby a single genome expresses different phenotypes in response to environmental cues, is a fascinating biological phenomenon. Social insects are especially intriguing examples of phenotypic plasticity because division of labor results in the development of extreme morphological phenotypes, such as the queen and worker castes. Although sociality evolved independently in ants, bees, wasps and termites, similar genetic pathways regulate phenotypic plasticity in these different groups of social insects. The insulin/insulin-like growth signaling (IIS) plays a key role in this process. Recent research reveals that IIS interacts with other pathways including target of rapamycin (TOR), epidermal growth factor receptor (Egfr), juvenile hormone (JH) and vitellogenin (Vg) to regulate caste differentiation.

MicroRNAs in Honey Bee Caste Determination

The cellular mechanisms employed by some organisms to produce contrasting morphological and reproductive phenotypes from the same genome remains one of the key unresolved issues in biology. Honeybees (Apis mellifera) use differential feeding and a haplodiploid sex determination system to generate three distinct organismal outcomes from the same genome. Here we investigate the honeybee female and male caste-specific microRNA and transcriptomic molecular signatures during a critical time of larval development. Both previously undetected and novel miRNAs have been discovered, expanding the inventory of these genomic regulators in invertebrates. We show significant differences in the microRNA and transcriptional profiles of diploid females relative to haploid drone males as well as between reproductively distinct females (queens and workers). Queens and drones show gene enrichment in physio-metabolic pathways, whereas workers show enrichment in processes associated with neuronal development, cell signalling and caste biased structural differences. Interestingly, predicted miRNA targets are primarily associated with non-physio-metabolic genes, especially neuronal targets, suggesting a mechanistic disjunction from DNA methylation that regulates physio-metabolic processes. Accordingly, miRNA targets are under-represented in methylated genes. Our data show how a common set of genetic elements are differentially harnessed by an organism, which may provide the remarkable level of developmental flexibility required.

Sphingolipids, Transcription Factors, and Conserved Toolkit Genes: Developmental Plasticity in the Ant Cardiocondyla obscurior

Developmental plasticity allows for the remarkable morphological specialization of individuals into castes in eusocial species of Hymenoptera. Developmental trajectories that lead to alternative caste fates are typically determined by specific environmental stimuli that induce larvae to express and maintain distinct gene expression patterns. Although most eusocial species express two castes, queens and workers, the ant Cardiocondyla obscurior expresses diphenic females and males; this provides a unique system with four discrete phenotypes to study the genomic basis of developmental plasticity in ants. We sequenced and analyzed the transcriptomes of 28 individual C. obscurior larvae of known developmental trajectory, providing the first in-depth analysis of gene expression in eusocial insect larvae. Clustering and transcription factor binding site analyses revealed that different transcription factors and functionally distinct sets of genes are recruited during larval development to induce the four alternative trajectories. In particular, we found complex patterns of gene regulation pertaining to sphingolipid metabolism, a conserved molecular pathway involved in development, obesity, and aging.

Age at which larvae are orphaned determines their development into typical or rebel workers in the honeybee (Apis mellifera L.)

In the honeybee, diploid larvae fed with royal jelly develop into reproductive queens, whereas larvae fed with royal jelly for three days only and subsequently with honey and pollen develop into facultatively sterile workers. A recent study showed that worker larvae fed in a queenless colony develop into another female polyphenic form: rebel workers. These rebel workers are more queenlike and have greater reproductive potential than normal workers. However, it was unclear whether larvae orphaned at any time during their feeding period can develop into rebels. To answer this question, the anatomical features of newly emerged workers reared in queenless conditions at different ages during the larval period were evaluated. Our results showed that larvae orphaned during the final four or more days of their feeding life develop into rebel workers with more ovarioles in their ovaries, smaller hypopharyngeal glands, and larger mandibular and Dufour’s glands compared with typical workers with low reproductive potential that were reared with a queen or orphaned at the third to last or a later day of feeding life.

Quantitative differences in nourishment affect caste-related physiology and development in the paper wasp Polistes metricus

The distinction between worker and reproductive castes of social insects is receiving increased attention from a developmental rather than adaptive perspective. In the wasp genus Polistes, colonies are founded by one or more females, and the female offspring that emerge in that colony are either non-reproducing workers or future reproductives of the following generation (gynes). A growing number of studies now indicate that workers emerge with activated reproductive physiology, whereas the future reproductive gynes do not. Low nourishment levels for larvae during the worker-rearing phase of the colony cycle and higher nourishment levels for larvae when gynes are reared are now strongly suspected of playing a major role in this difference.

Here, we present the results of a laboratory rearing experiment in which Polistes metricus single foundresses were held in environmental conditions with a higher level of control than in any previously published study, and the amount of protein nourishment made available to feed larvae was the only input variable. Three experimental feeding treatments were tested: restricted, unrestricted, and hand-supplemented. Analysis of multiple response variables shows that wasps reared on restricted protein nourishment, which would be the case for wasps reared in field conditions that subsequently become workers, tend toward trait values that characterize active reproductive physiology. Wasps reared on unrestricted and hand-supplemented protein, which replicates higher feeding levels for larvae in field conditions that subsequently become gynes, tend toward trait values that characterize inactive reproductive physiology. Although the experiment was not designed to test for worker behavior per se, our results further implicate activated reproductive physiology as a developmental response to low larval nourishment as a fundamental aspect of worker behavior in Polistes.

Dimorphic ovary differentiation in honeybee (Apis mellifera) larvae involves caste-specific expression of homologs of ark and buffy cell death genes

The establishment of the number of repeated structural units, the ovarioles, in the ovaries is one of the critical events that shape caste polyphenism in social insects. In early postembryonic development, honeybee (Apis mellifera) larvae have a pair of ovaries, each one consisting of almost two hundred ovariole primordia. While practically all these ovarioles continue developing in queen-destined larvae, they undergo massive programmed cell death (PCD) in worker-destined larvae. So as to gain insight into the molecular basis of this fundamental process in caste differentiation we used quantitative PCR (qPCR) and fluorescent in situ hybridization (FISH) to investigate the expression of the Amark and Ambuffy genes in the ovaries of the two honeybee castes throughout the fifth larval instar. These are the homologs of ark and buffy Drosophila melanogaster genes, respectively, involved in activating and inhibiting PCD. Caste-specific expression patterns were found during this time-window defining ovariole number. Amark transcript levels were increased when ovariole resorption was intensified in workers, but remained at low levels in queen ovaries. The transcripts were mainly localized at the apical end of all the worker ovarioles, but appeared in only a few queen ovarioles, thus strongly suggesting a function in mediating massive ovariolar cell death in worker larvae. Ambuffy was mainly expressed in the peritoneal sheath cells covering each ovariole. The levels of Ambuffy transcripts increased earlier in the developing ovaries of queens than in workers. Consistent with a protective role against cell death, Ambuffy transcripts were localized in practically all queen ovarioles, but only in few worker ovarioles. The results are indicative of a functional relationship between the expression of evolutionary conserved cell death genes and the morphological events leading to caste-specific ovary differentiation in a social insect.

Biased gene expression in early honeybee larval development

Background

Female larvae of the honeybee (Apis mellifera) develop into either queens or workers depending on nutrition. This nutritional stimulus triggers different developmental trajectories, resulting in adults that differ from each other in physiology, behaviour and life span.

Results

To understand how these trajectories are established we have generated a comprehensive atlas of gene expression throughout larval development. We found substantial differences in gene expression between worker and queen-destined larvae at 6 hours after hatching. Some of these early changes in gene expression are maintained throughout larval development, indicating that caste-specific developmental trajectories are established much earlier than previously thought. Within our gene expression data we identified processes that potentially underlie caste differentiation. Queen-destined larvae have higher expression of genes involved in transcription, translation and protein folding early in development with a later switch to genes involved in energy generation. Using RNA interference, we were able to demonstrate that one of these genes, hexamerin 70b, has a role in caste differentiation. Both queen and worker developmental trajectories are associated with the expression of genes that have alternative splice variants, although only a single variant of a gene tends to be differentially expressed in a given caste.

Conclusions

Our data, based on the biases in gene expression early in development together with published data, supports the idea that caste development in the honeybee consists of two phases; an initial biased phase of development, where larvae can still switch to the other caste by differential feeding, followed by commitment to a particular developmental trajectory.

A cell line resource derived from honey bee (Apis mellifera) embryonic tissues

A major hindrance to the study of honey bee pathogens or the effects of pesticides and nutritional deficiencies is the lack of controlled in vitro culture systems comprised of honey bee cells. Such systems are important to determine the impact of these stress factors on the developmental and cell biology of honey bees. We have developed a method incorporating established insect cell culture techniques that supports sustained growth of honey bee cells in vitro. We used honey bee eggs mid to late in their embryogenesis to establish primary cultures, as these eggs contain cells that are progressively dividing. Primary cultures were initiated in modified Leibovitz’s L15 medium and incubated at 32^°C. Serial transfer of material from several primary cultures was maintained and has led to the isolation of young cell lines. A cell line (AmE-711) has been established that is composed mainly of fibroblast-type cells that form an adherent monolayer. Most cells in the line are diploid (2n = 32) and have the Apis mellifera karyotype as revealed by Giemsa stain. The partial sequence for the mitochondrial-encoded cytochrome c oxidase subunit I (Cox 1) gene in the cell line is identical to those from honey bee tissues and a consensus sequence for A. mellifera. The population doubling time is approximately 4 days. Importantly, the cell line is continuously subcultured every 10–14 days when split at a 1:3 ratio and is cryopreserved in liquid nitrogen. The cell culture system we have developed has potential application for studies aimed at honey bee development, genetics, pathogenesis, transgenesis, and toxicology.

Transcriptome comparison between honey bee queen- and worker-destined larvae

Caste differentiation in the female honey bee is one of the most intriguing polyphenism phenomena. This developmental switch depends on the differential expression of entire suites of the genes involved in the larval fate between the queens and workers. In this study, we compared the transcriptome differences between full-sister queen- (QL) and worker-destined larvae (WL) using high-throughput RNA-Seq. QL and WL at fourth (L4) and fifth instar (L5) were used to prepare four libraries and to generate 50,191,699 (QL4), 57,628,541 (WL4), 56,613,619 (QL5), and 58,626,829 (WL5) usable reads, which were assembled into groups of 7,952, 7,993, 7,971, and 8,023 genes, respectively. The transcriptome changes were investigated using the DEGs Package (DEGseq), which resulted in more than 4,500 differentially expressed genes (DEGs) between the castes. Eight of the DEGs were verified by quantitative real-time RT-PCR (qRT-PCR), and the results supported our sequencing data. All of the DEGs were analysed using Web Gene Ontology Annotation Plot (WEGO) and then mapped using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. These results suggest that over 70% of the DEGs in each instar were more highly expressed in QL than in WL, possibly suggesting that the QL genes had higher transcriptional activity than the WL genes during differentiation. The same gene set is active (but differentially expressed) in both castes, which in turn result in dimorphic females. The L4 stage is a very active gene expression period for both QL and WL before their pupal stage. The activity of the mTOR (a target of rapamycin) encoding gene in the mTOR signalling pathway is higher in QL4 than in WL4, and this difference was no longer present by the L5 feeding stage. The genes down-stream of mTOR maintained this change at the L5 stage. These results could contribute to an in-depth study of the candidate genes during honey bee caste differentiation and improve our current understanding of the polyphenism phenomenon in insects.

Physiological variation as a mechanism for developmental caste-biasing in a facultatively eusocial sweat bee

Social castes of eusocial insects may have arisen through an evolutionary modification of an ancestral reproductive ground plan, such that some adults emerge from development physiologically primed to specialize on reproduction (queens) and others on maternal care expressed as allo-maternal behaviour (workers). This hypothesis predicts that variation in reproductive physiology should emerge from ontogeny and underlie division of labour. To test these predictions, we identified physiological links to division of labour in a facultatively eusocial sweat bee, Megalopta genalis. Queens are larger, have larger ovaries and have higher vitellogenin titres than workers. We then compared queens and workers with their solitary counterparts-solitary reproductive females and dispersing nest foundresses-to investigate physiological variation as a factor in caste evolution. Within dyads, body size and ovary development were the best predictors of behavioural class. Queens and dispersers are larger, with larger ovaries than their solitary counterparts. Finally, we raised bees in social isolation to investigate the influence of ontogeny on physiological variation. Body size and ovary development among isolated females were highly variable, and linked to differences in vitellogenin titres. As these are key physiological predictors of social caste, our results provide evidence for developmental caste-biasing in a facultatively eusocial bee.

Epigenetics in social insects: a new direction for understanding the evolution of castes

Epigenetic modifications to DNA, such as DNA methylation, can expand a genome's regulatory flexibility, and thus may contribute to the evolution of phenotypic plasticity. Recent work has demonstrated the importance of DNA methylation in alternative queen and worker “castes” in social insects, particularly honeybees. Social insects are an excellent system for addressing questions about epigenetics and evolution because: (1) they have dramatic caste polyphenisms that appear to be tied to differential methylation, (2) DNA methylation is widespread in various groups of social insects, and (3) there are intriguing connections between the social environment and DNA methylation in many species, from insects to mammals. In this article, we review research on honeybees, and, when available, other social insects, on DNA methylation and queen and worker caste differences. We outline a conceptual framework for the effects of methylation on caste determination in honeybees that may help guide studies of epigenetic regulation in other polyphenic taxa. Finally, we suggest future paths of study for social insect epigenetic research, including the importance of comparative studies of DNA methylation on a broader range of species, and highlight some key unanswered mechanistic questions about how DNA methylation affects gene regulation.

Representational Difference Analysis (RDA) reveals differential expression of conserved as well as novel genes during caste-specific development of the honey bee (Apis mellifera L.) ovary

In highly eusocial insects, such as the honey bee, Apis mellifera, the reproductive bias has become embedded in morphological caste differences. These are most expressively denoted in ovary size, with adult queens having large ovaries consisting of 150-200 ovarioles each, while workers typically have only 1-20 ovarioles per ovary. This morphological differentiation is a result of hormonal signals triggered by the diet change in the third larval instar, which eventually generate caste-specific gene expression patterns. To reveal these we produced differential gene expression libraries by Representational Difference Analysis (RDA) for queen and worker ovaries in a developmental stage when cell death is a prominent feature in the ovarioles of workers, whereas all ovarioles are maintained and extend in length in queens. In the queen library, 48% of the gene set represented homologs of known Drosophila genes, whereas in the worker ovary, the largest set (59%) were ESTs evidencing novel genes, not even computationally predicted in the honey bee genome. Differential expression was confirmed by quantitative RT-PCR for a selected gene set, denoting major differences for two queen and two worker library genes. These included two unpredicted genes located in chromosome 11 (Group11.35 and Group11.31, respectively) possibly representing long non-coding RNAs. Being candidates as modulators of ovary development, their expression and functional analysis should be a focal point for future studies.

Global identification of transcription start sites in the genome of Apis mellifera using 5'LongSAGE

The precise identification of the transcription start sites (TSSs) of genes in the honeybee genome will be helpful for inferring start codons and for determining promoter elements. The 5'SAGE approach provides a powerful tool for identifying TSSs in the sequenced genome. The main purpose of this study is to identify the actual TSSs of expressed genes as well as the usage of different TSSs in the Apis mellifera genome. We performed a 5'LongSAGE (5'LS) analysis for the adult drone head, and the TSSs of the expressed genes were determined by mapping the 5'LS tag sequences to the honeybee genome. A total of 8,280 unique 19 bp 5'LS tag sequences were identified that corresponded to 3,655 predicted genes. Out of these tags, 4,998 tags (60.4%) were mapped to a region from -1,000 bp to +100 bp of the start codon of 2,301 reference coding sequences. Notably, we observed that 28-47% of the 3,655 honeybee genes initiated transcription from alternative TSSs. The TSS consensus pattern of the honeybee genes, DT(rich) PyPu(G(rich))(T/A)(T(rich))(3), was obtained by aligning the sequences flanking the 5'LS-TSSs. We also identified three new genes in the regions downstream of 5'LS tags and validated 21 TSSs using RT-PCR amplification. Additionally, 17 genes identified by the 5'LS tags were associated with the Gene Ontology term "behavior." Mapping of the 5'LS tags on the genome not only provided direct evidence of expression for in silico predicted genes but also allowed for the identification of previously unrecognized, novel exons and alternative TSSs.

Differential expression of hypoxia pathway genes in honey bee (Apis mellifera L.) caste development

Diphenism in social bees is essentially contingent on nutrient-induced cellular and systemic physiological responses resulting in divergent gene expression patterns. Analyses of juvenile hormone (JH) titers and functional genomics assays of the insulin-insulin-like signaling (IIS) pathway and its associated branch, target-of-rapamycin (TOR), revealed systemic responses underlying honey bee (Apis mellifera) caste development. Nevertheless, little attention has been paid to cellular metabolic responses. Following up earlier investigations showing major caste differences in oxidative metabolism and mitochondrial physiology, we herein identified honey bee homologs of hypoxia signaling factors, HIFα/Sima, HIFβ/Tango and PHD/Fatiga and we investigated their transcript levels throughout critical stages of larval development. Amsima, Amtango and Amfatiga showed correlated transcriptional activity, with two peaks of occurring in both queens and workers, the first one shortly after the last larval molt and the second during the cocoon-spinning phase. Transcript levels for the three genes were consistently higher in workers. As there is no evidence for major microenvironmental differences in oxygen levels within the brood nest area, this appears to be an inherent caste character. Quantitative PCR analyses on worker brain, ovary, and leg imaginal discs showed that these tissues differ in transcript levels. Being a highly conserved pathway and linked to IIS/TOR, the hypoxia gene expression pattern seen in honey bee larvae denotes that the hypoxia pathway has undergone a transformation, at least during larval development, from a response to environmental oxygen concentrations to an endogenous regulatory factor in the diphenic development of honey bee larvae.

Next-generation small RNA sequencing for microRNAs profiling in the honey bee Apis mellifera

MicroRNAs (miRNAs) are key regulators in various physiological and pathological processes via post-transcriptional regulation of gene expression. The honey bee (Apis mellifera) is a key model for highly social species, and its complex social behaviour can be interpreted theoretically as changes in gene regulation, in which miRNAs are thought to be involved. We used the SOLiD sequencing system to identify the repertoire of miRNAs in the honey bee by sequencing a mixed small RNA library from different developmental stages. We obtained a total of 36,796,459 raw sequences; of which 5,491,100 short sequences were fragments of mRNA and other noncoding RNAs (ncRNA), and 1,759,346 reads mapped to the known miRNAs. We predicted 267 novel honey bee miRNAs representing 380,182 short reads, including eight miRNAs of other insects in 14,107,583 genome-mapped sequences. We verified 50 of them using stem-loop reverse-transcription PCR (RT-PCR), in which 35 yielded PCR products. Cross-species analyses showed 81 novel miRNAs with homologues in other insects, suggesting that they were authentic miRNAs and have similar functions. The results of this study provide a basis for studies of the miRNA-modulating networks in development and some intriguing phenomena such as caste differentiation in A. mellifera.

Differential protein expression in honeybee (Apis mellifera L.) larvae: underlying caste differentiation

Honeybee (Apis mellifera) exhibits divisions in both morphology and reproduction. The queen is larger in size and fully developed sexually, while the worker bees are smaller in size and nearly infertile. To better understand the specific time and underlying molecular mechanisms of caste differentiation, the proteomic profiles of larvae intended to grow into queen and worker castes were compared at 72 and 120 hours using two dimensional electrophoresis (2-DE), network, enrichment and quantitative PCR analysis. There were significant differences in protein expression between the two larvae castes at 72 and 120 hours, suggesting the queen and the worker larvae have already decided their fate before 72 hours. Specifically, at 72 hours, queen intended larvae over-expressed transketolase, aldehyde reductase, and enolase proteins which are involved in carbohydrate metabolism and energy production, imaginal disc growth factor 4 which is a developmental related protein, long-chain-fatty-acid CoA ligase and proteasome subunit alpha type 5 which metabolize fatty and amino acids, while worker intended larvae over-expressed ATP synthase beta subunit, aldehyde dehydrogenase, thioredoxin peroxidase 1 and peroxiredoxin 2540, lethal (2) 37 and 14-3-3 protein epsilon, fatty acid binding protein, and translational controlled tumor protein. This differential protein expression between the two caste intended larvae was more pronounced at 120 hours, with particular significant differences in proteins associated with carbohydrate metabolism and energy production. Functional enrichment analysis suggests that carbohydrate metabolism and energy production and anti-oxidation proteins play major roles in the formation of caste divergence. The constructed network and validated gene expression identified target proteins for further functional study. This new finding is in contrast to the existing notion that 72 hour old larvae has bipotential and can develop into either queen or worker based on epigenetics and can help us to gain new insight into the time of departure as well as caste trajectory influencing elements at the molecular level.

Differential gene expression and protein abundance evince ontogenetic bias toward castes in a primitively eusocial wasp

Polistes paper wasps are models for understanding conditions that may have characterized the origin of worker and queen castes and, therefore, the origin of paper wasp sociality. Polistes is “primitively eusocial” by virtue of having context-dependent caste determination and no morphological differences between castes. Even so, Polistes colonies have a temporal pattern in which most female larvae reared by the foundress become workers, and most reared by workers become future-reproductive gynes. This pattern is hypothesized to reflect development onto two pathways, which may utilize mechanisms that regulate diapause in other insects. Using expressed sequence tags (ESTs) for Polistes metricus we selected candidate genes differentially expressed in other insects in three categories: 1) diapause vs. non-diapause phenotypes and/or worker vs. queen differentiation, 2) behavioral subcastes of worker honey bees, and 3) no a priori expectation of a role in worker/gyne development. We also used a non-targeted proteomics screen to test for peptide/protein abundance differences that could reflect larval developmental divergence. We found that foundress-reared larvae (putative worker-destined) and worker-reared larvae (putative gyne-destined) differed in quantitative expression of sixteen genes, twelve of which were associated with caste and/or diapause in other insects, and they also differed in abundance of nine peptides/proteins. Some differentially-expressed genes are involved in diapause regulation in other insects, and other differentially-expressed genes and proteins are involved in the insulin signaling pathway, nutrient metabolism, and caste determination in highly social bees. Differential expression of a gene and a peptide encoding hexameric storage proteins is especially noteworthy. Although not conclusive, our results support hypotheses of 1) larval developmental pathway divergence that can lead to caste bias in adults and 2) nutritional differences as the foundation of the pathway divergence. Finally, the differential expression in Polistes larvae of genes and proteins also differentially expressed during queen vs. worker caste development in honey bees may indicate that regulatory mechanisms of caste outcomes share similarities between primitively eusocial and advanced eusocial Hymenoptera.

DNA methylation is widespread and associated with differential gene expression in castes of the honeybee, Apis mellifera

The recent, unexpected discovery of a functional DNA methylation system in the genome of the social bee Apis mellifera underscores the potential importance of DNA methylation in invertebrates. The extent of genomic DNA methylation and its role in A. mellifera remain unknown, however. Here we show that genes in A. mellifera can be divided into 2 distinct classes, one with low-CpG dinucleotide content and the other with high-CpG dinucleotide content. This dichotomy is explained by the gradual depletion of CpG dinucleotides, a well-known consequence of DNA methylation. The loss of CpG dinucleotides associated with DNA methylation also may explain the unusual mutational patterns seen in A. mellifera that lead to AT-rich regions of the genome. A detailed investigation of this dichotomy implicates DNA methylation in A. mellifera development. High-CpG genes, which are predicted to be hypomethylated in germlines, are enriched with functions associated with developmental processes, whereas low-CpG genes, predicted to be hypermethylated in germlines, are enriched with functions associated with basic biological processes. Furthermore, genes more highly expressed in one caste than another are overrepresented among high-CpG genes. Our results highlight the potential significance of epigenetic modifications, such as DNA methylation, in developmental processes in social insects. In particular, the pervasiveness of DNA methylation in the genome of A. mellifera provides fertile ground for future studies of phenotypic plasticity and genomic imprinting.

Genetic and genomic analyses of the division of labour in insect societies

Division of labour--individuals specializing in different activities--features prominently in the spectacular success of the social insects. Until recently, genetic and genomic analyses of division of labour were limited to just a few species. However, research on an ever-increasing number of species has provided new insight, from which we highlight two results. First, heritable influences on division of labour are more pervasive than previously imagined. Second, different forms of division of labour, in lineages in which eusociality has arisen independently, have evolved through changes in the regulation of highly conserved molecular pathways associated with several basic life-history traits, including nutrition, metabolism and reproduction.

Genome-wide analysis reveals differences in brain gene expression patterns associated with caste and reproductive status in honey bees (Apis mellifera)

A key characteristic of eusocial species is reproductive division of labour. Honey bee colonies typically have a single reproductive queen and thousands of sterile workers. Adult queens differ dramatically from workers in anatomy, physiology, behaviour and lifespan. Young female workers can activate their ovaries and initiate egg laying; these 'reproductive' workers differ from sterile workers in anatomy, physiology, and behaviour. These differences, however, are on a much smaller scale than those observed between the queen and worker castes. Here, we use microarrays to monitor expression patterns of several thousand genes in the brains of same-aged virgin queens, sterile workers, and reproductive workers. We found large differences in expression between queens and both worker groups (~2000 genes), and much smaller differences between sterile and reproductive workers (221 genes). The expression patterns of these 221 genes in reproductive workers are more queen-like, and may represent a core group of genes associated with reproductive physiology. Furthermore, queens and reproductive workers preferentially up-regulate genes associated with the nurse bee behavioural state, which supports the hypothesis of an evolutionary link between worker division of labour and molecular pathways related to reproduction. Finally, several functional groups of genes associated with longevity in other species are significantly up-regulated in queens. Identifying the genes that underlie the differences between queens, sterile workers, and reproductive workers will allow us to begin to characterize the molecular mechanisms underlying the evolution of social behaviour and large-scale remodelling of gene networks associated with polyphenisms.

Insights into social insects from the genome of the honeybee Apis mellifera

Here we report the genome sequence of the honeybee Apis mellifera, a key model for social behaviour and essential to global ecology through pollination. Compared with other sequenced insect genomes, the A. mellifera genome has high A+T and CpG contents, lacks major transposon families, evolves more slowly, and is more similar to vertebrates for circadian rhythm, RNA interference and DNA methylation genes, among others. Furthermore, A. mellifera has fewer genes for innate immunity, detoxification enzymes, cuticle-forming proteins and gustatory receptors, more genes for odorant receptors, and novel genes for nectar and pollen utilization, consistent with its ecology and social organization. Compared to Drosophila, genes in early developmental pathways differ in Apis, whereas similarities exist for functions that differ markedly, such as sex determination, brain function and behaviour. Population genetics suggests a novel African origin for the species A. mellifera and insights into whether Africanized bees spread throughout the New World via hybridization or displacement.

Caste development and reproduction: a genome-wide analysis of hallmarks of insect eusociality

The honey bee queen and worker castes are a model system for developmental plasticity. We used established expressed sequence tag information for a Gene Ontology based annotation of genes that are differentially expressed during caste development. Metabolic regulation emerged as a major theme, with a caste-specific difference in the expression of oxidoreductases vs. hydrolases. Motif searches in upstream regions revealed group-specific motifs, providing an entry point to cis-regulatory network studies on caste genes. For genes putatively involved in reproduction, meiosis-associated factors came out as highly conserved, whereas some determinants of embryonic axes either do not have clear orthologs (bag of marbles, gurken, torso), or appear to be lacking (trunk) in the bee genome. Our results are the outcome of a first genome-based initiative to provide an annotated framework for trends in gene regulation during female caste differentiation (representing developmental plasticity) and reproduction.

Genes of the antioxidant system of the honey bee: annotation and phylogeny

Antioxidant enzymes perform a variety of vital functions including the reduction of life-shortening oxidative damage. We used the honey bee genome sequence to identify the major components of the honey bee antioxidant system. A comparative analysis of honey bee with Drosophila melanogaster and Anopheles gambiae shows that although the basic components of the antioxidant system are conserved, there are important species differences in the number of paralogs. These include the duplication of thioredoxin reductase and the expansion of the thioredoxin family in fly; lack of expansion of the Theta, Delta and Omega GST classes in bee and no expansion of the Sigma class in dipteran species. The differential expansion of antioxidant gene families among honey bees and dipteran species might reflect the marked differences in life history and ecological niches between social and solitary species.

Differential expression of cytochrome c oxidase subunit III gene in castes of the termite Reticulitermes santonensis

Social insects such as termites live in colonies in which cooperation is assumed by all individuals developing into castes to which specific tasks are allocated. Little has been reported about molecular aspects underlying termite caste-specific gene expression. Genetic regulation has recently been hypothesized to govern caste-specific traits and physiology in social insects. Cytochrome c oxidase (COX) has been shown to be an interesting candidate for expression study in insects. We used the cytochrome c oxidase subunit III gene (COXIII) that was cloned from mRNA in a lower termite, Reticulitermes santonensis De Feytaud (Isoptera; Rhinotermitidae). The full-length cDNA encodes a protein of 262 amino acids that shows high degree of homology with other insects COXIIIs. Reverse transcriptase-PCR and real-time PCR were performed to compare gene expression between larvae, workers, nymphs and soldiers. Analyses performed on head cDNAs revealed that COXIII is differentially expressed between castes. The level of COXIII is caste-regulated with an increase in workers (approximately 1.9-fold) and nymphs (approximately 2.8-fold) and a decrease in soldiers (0.8-fold) compared to the expression level in larvae (1.0-fold). These results may emphasize the physiological importance of COX in the termite brain at different developmental stages.

Differential gene expression and phenotypic plasticity in behavioural castes of the primitively eusocial wasp, Polistes canadensis

Understanding how a single genome can produce a variety of different phenotypes is of fundamental importance in evolutionary and developmental biology. One of the most striking examples of phenotypic plasticity is the female caste system found in eusocial insects, where variation in reproductive (queens) and non-reproductive (workers) phenotypes results in a broad spectrum of caste types, ranging from behavioural through to morphological castes. Recent advances in genomic techniques allow novel comparisons on the nature of caste phenotypes to be made at the level of the genes in organisms for which there is little genome information, facilitating new approaches in studying social evolution and behaviour. Using the paper wasp Polistes canadensis as a model system, we investigated for the first time how behavioural castes in primitively eusocial insect societies are associated with differential expression of shared genes. We found that queens and newly emerged females express gene expression patterns that are distinct from each other whilst workers generally expressed intermediate patterns, as predicted by Polistes biology. We compared caste-associated genes in P. canadensis with those expressed in adult queens and workers of more advanced eusocial societies, which represent four independent origins of eusociality. Nine genes were conserved across the four taxa, although their patterns of expression and putative functions varied. Thus, we identify several genes that are putatively of evolutionary importance in the molecular biology that underlies a number of caste systems of independent evolutionary origin.

The role of mitochondrial respiration in physiological and evolutionary adaptation

Aerobic mitochondria serve as the power sources of eukaryotes by producing ATP through oxidative phosphorylation (OXPHOS). The enzymes involved in OXPHOS are multisubunit complexes encoded by both nuclear and mitochondrial DNA. Thus, regulation of respiration is necessarily a highly coordinated process that must organize production, assembly and function of mitochondria to meet an organism's energetic needs. Here I review the role of OXPHOS in metabolic adaptation and diversification of higher animals. On a physiological timescale, endocrine-initiated signaling pathways allow organisms to modulate respiratory enzyme concentration and function under changing environmental conditions. On an evolutionary timescale, mitochondrial enzymes are targets of natural selection, balancing cytonuclear coevolutionary constraints against physiological innovation. By synthesizing our knowledge of biochemistry, physiology and evolution of respiratory regulation, I propose that we can now explore questions at the interface of these fields, from molecular translation of environmental cues to selection on mitochondrial haplotype variation.

The evolutionary origin and elaboration of sociality in the aculeate Hymenoptera: maternal effects, sib-social effects, and heterochrony

We discuss the evolutionary origin and elaboration of sociality using an indirect genetic effects perspective. Indirect genetic effects models simultaneously consider zygotic genes, genes expressed in social partners (especially mothers and siblings), and the interactions between them. Incorporation of these diverse genetic effects should lead to more realistic models of social evolution. We first review haplodiploidy as a factor that promotes the evolution of eusociality. Social insect biologists have doubted the importance of relatedness asymmetry caused by haplodiploidy and focused on other predisposing factors such as maternal care. However; indirect effects theory shows that maternal care evolves more readily in haplodiploids, especially with inbreeding and despite multiple mating. Because extended maternal care is believed to be a precondition for the evolution of eusociality, the evolutionary bias towards maternal care in haplodiploids may result in a further bias towards eusociality in these groups. Next, we compare kin selection and parental manipulation and then briefly review additional hypotheses for the evolutionary origin of eusociality. We present a verbal model for the evolutionary origin and elaboration of sib-social care from maternal care based on the modification of the timing of expression of maternal care behaviors. Specifically, heterochrony genes cause maternal care behaviors to be expressed prereproductively towards siblings instead of postreproductively towards offspring. Our review demonstrates that both maternal effect genes (expressed in a parental manipulation manner) and direct effect zygotic genes (expressed in an offspring control manner) are likely involved in the evolution of eusociality. We conclude by describing theoretical and empirical advances with indirect genetic effects and sociogenomics, and we provide specific quantitative genetic and genomic predictions from our heterochrony model for the evolutionary origin and elaboration of eusociality.

Gene expression patterns associated with queen honey bee longevity

The oxidative stress theory of aging proposes that accumulation of oxidative damage is the main proximate cause of aging and that lifespan is determined by the rate at which this damage occurs. Two predictions from this theory are that long-lived organisms produce fewer ROS or have increased antioxidant production. Based in these predictions, molecular mechanisms to promote longevity could include either changes in the regulation of mitochondrial genes that affect ROS production or elevated expression of antioxidant genes. We explored these possibilities in the honey bee, a good model for the study of aging because it has a caste system in which the same genome produces both a long-lived queen and a short-lived worker. We measured mRNA levels for genes encoding eight of the most prominent antioxidant enzymes and five mitochondrial proteins involved in respiration. The expression of antioxidant genes generally decreased with age in queens, but not in workers. Expression of most mitochondrial genes, in particular CytC, was higher in young queens, but these genes showed a faster age-related decline relative to workers. One exception to this trend was COX-I in thorax. This resulted in higher COX-I/CytC ratios in old queens compared to old workers, which suggests caste-specific differences in mitochondrial function that might be related to the caste-specific differences in longevity. Queen honey bee longevity appears to have evolved via mechanisms other than increased antioxidant gene expression.

Differential gene expression in queen-worker caste determination in bumble-bees

Investigating how differential gene expression underlies caste determination in the social Hymenoptera is central to understanding how variation in gene expression underlies adaptive phenotypic diversity. We investigated for the first time the association between differential gene expression and queen-worker caste determination in the bumble-bee Bombus terrestris. Using suppression subtractive hybridization we isolated 12 genes that were differentially expressed in queen- and worker-destined larvae. We found that the sets of genes underlying caste differences in larvae and adults failed to overlap greatly. We also found that B. terrestris shares some of the genes whose differential expression is associated with caste determination in the honeybee, Apis mellifera, but their expression patterns were not identical. Instead, we found B. terrestris to exhibit a novel pattern, whereby most genes upregulated (i.e. showing relatively higher levels of expression) in queen-destined larvae early in development were upregulated in worker-destined larvae late in development. Overall, our results suggest that caste determination in B. terrestris involves a difference not so much in the identity of genes expressed by queen- and worker-destined larvae, but primarily in the relative timing of their expression. This conclusion is of potential importance in the further study of phenotypic diversification via differential gene expression.

Molecular characterization of a cDNA encoding prophenoloxidase and its expression in Apis mellifera

Phenoloxidase (PO), a melanin-synthesizing enzyme known to play an important role in insect defense, is found as a zymogen (ProPO) in hemolymph and cuticle, where it is activated by proteolysis. We characterized the first proPO cDNA in an eusocial insect, the Apis mellifera honey bee. The AmproPO cDNA contains an ORF of 2079 bp encoding 693 amino acids, and is composed of 9 exons and 8 introns. Southern blot of digested genomic DNA suggested that only one copy of the proPO gene is present in A. mellifera. The molecular mass of the deduced ProPO and the active enzyme was predicted to be 80.1 and 74.4 kDa, respectively. The calculated pI was 6.28. BLASTp search of the deduced amino acid sequence, and neighbor-joining analysis, showed similarity with ProPOs from other insects, ranging from 47% to 63%. Protein signature analyses revealed four conserved regions, including the two copper binding sites characteristic of arthropod ProPOs. RT-PCR and Southern blot showed the highest amount of AmproPO transcripts in workers whole body, followed by queens and drones. Expression was also detected in hemocytes and integument. Real time RT-PCR showed higher amounts of AmproPO transcripts in adults and older pupae than in younger pupae and larvae, suggesting a function of AmproPO in adult exoskeleton melanization and differentiation.

Sociogenomics: social life in molecular terms

Spectacular progress in molecular biology, genome-sequencing projects and genomics makes this an appropriate time to attempt a comprehensive understanding of the molecular basis of social life. Promising results have already been obtained in identifying genes that influence animal social behaviour and genes that are implicated in social evolution. These findings - derived from an eclectic mix of species that show varying levels of sociality - provide the foundation for the integration of molecular biology, genomics, neuroscience, behavioural biology and evolutionary biology that is necessary for this endeavour.

Gene expression profiles among immature and adult reproductive castes of the termite Reticulitermes flavipes

Array-based genomic studies were conducted with the goal of identifying immature (i.e. nymph) and adult reproductive caste-biased gene expression in the termite Reticulitermes flavipes. Using cDNA macro-arrays, we identified thirty-four nymph-biased genes falling into eight ontogenic categories. Based on gene expression profiles among diverse castes and developmental stages (determined by quantitative PCR), several important trends emerged. These findings highlight the importance of several developmental and survival-based factors among immature and adult termite reproductives, including: vitellogenesis, nutrient storage, juvenile hormone sequestration, ribosomal translational and filtering mechanisms, fatty acid biosynthesis, apoptosis inhibition, and both endogenous and symbiont cellulase-assisted nutrition. These findings are highly significant as they are the first to elucidate the molecular biology underlying termite reproductive caste differentiation and reproductive caste-specific biology. Other gene expression results are in agreement with previous findings that suggest roles for vitellogenin-like haemolymph proteins in soldier caste differentiation.

Differential gene expression between alate and dealate queens in the red imported fire ant, Solenopsis invicta Buren (Hymenoptera: Formicidae)

The transition of fire ant queens from alates to dealates, following a mating flight, is associated with numerous important physiological changes. A molecular analysis of gene expression differences that occur between alates and dealates was performed using the suppression subtractive hybridization (SSH) method. 983 SSH clones were arrayed and screened by dot blot hybridization, followed by Northern blot analysis for selected clones. Gene expression profiles throughout fire ant development were determined using semi-quantitative reverse transcriptase polymerase chain reactions (RT-PCR). The cytochrome c oxidase subunit II and STARS (striated muscle activator of Rho signaling) transcripts were expressed at higher levels in dealates compared to alates and may be involved in the programmed cell death of the flight muscles. Three different vitellogenin genes and two unique yellow g-like genes were identified that may be closely associated with the reproductive system and/or nutrient transport. Two putative antibacterial peptides, abaecin and hymenoptaecin precursors, were highly expressed in dealate queens, suggesting that they are present as an immune system component during this important stage of fire ant development. The genes identified in this study may be utilized as novel targets for fire ant control and will also provide molecular markers for studies of other social insects.

Tissue-specific expression of male-transmitted mitochondrial DNA and its implications for rates of molecular evolution in Mytilus mussels (Bivalvia: Mytilidae)

Mytilus and other bivalves exhibit an unusual system of mitochondrial DNA (mtDNA) transmission termed doubly uniparental inheritance (DUI). Specifically, males transmit the mtDNA they have received from their fathers to their sons. Females transmit their mother's mtDNA to both sons and daughters. Males are normally heteroplasmic and females are normally homoplasmic, but not exclusively. This system is associated with an unusual pattern of molecular evolution. Male-transmitted mtDNA (M type) evolves faster than female-transmitted (F type) mtDNA. Relatively relaxed selection on the M type has been proposed as an explanation for this phenomenon. To further evaluate the selective forces acting upon the M-type genome, we used RT-PCR to determine where it is expressed. M-type mtDNA expression was detected in all gonad samples and in 50% of somatic tissues of males, and in a single female tissue. F-type mtDNA expression was detected in all female tissues, all male somatic tissues, and all but one male gonad sample. We argue that the expression of M-type mtDNA in male somatic and male gonad tissues has implications for the strength of selection acting upon it.