A honey bee hexamerin, HEX 70a, is likely to play an intranuclear role in developing and mature ovarioles and testioles

Revealing Changes in Ovarian and Hemolymphatic Metabolites Using Widely Targeted Metabolomics between Newly Emerged and Laying Queens of Honeybee (Apis mellifera)

The queen bee is a central and pivotal figure within the colony, serving as the sole fertile female responsible for its reproduction. The queen possesses an open circulatory system, with her ovaries immersed in hemolymph. A continuous and intricate transportation and interchange of substances exist between the ovaries and hemolymph of queen bees. To determine the characteristic metabolites in the hemolymph and ovary, as well as understand how their rapid metabolism contributes to the process of egg-laying by queens, we reared Apis mellifera queens from three different age groups: newly emerged queen (NEQ), newly laying queen (NLQ), and old laying queen (OLQ). Using widely targeted metabolomics, our study revealed that the laying queen (NLQ and OLQ) exhibited faster fatty acid metabolism, up-regulated expression of antioxidants, and significant depletion of amino acids compared to the NEQ. This study revealed that the levels of carnitine and antioxidants (GSH, 2-O-α-D-glucopyranosyl-L-ascorbic acid, L-ascorbic acid 2-phosphate, etc.) in the NLQ and OLQ were significantly higher compared to NEQ. However, most of the differentially expressed amino acids, such as L-tryptophan, L-tyrosine, L-aspartic acid, etc., detected in NLQ and OLQ were down-regulated compared to the NEQ. Following egg-laying, pathways in the queens change significantly, e.g., Tryptophan metabolism, Tyrosine metabolism, cAMP signaling pathway, etc. Our results suggest that carnitine and antioxidants work together to maintain the redox balance of the queen. Additionally, various amino acids are responsible for maintaining the queen's egg production.

Genome-wide analysis of the chromatin sites targeted by HEX 70a storage protein in the honeybee brain and fat body

Hexamerins, the proteins massively stored in the larval haemolymph of insects, are gradually used throughout metamorphosis as a source of raw material and energy for the development of adult tissues. Such behaviour defined hexamerins as storage proteins. Immunofluorescence experiments coupled with confocal microscopy show a hexamerin, HEX 70a, in the nucleus of the brain and fat body cells from honeybee workers, an unexpected localization for a storage protein. HEX 70a colocalizes with fibrillarin, a nucleolar-specific protein and H3 histone, thus suggesting a potential role as a chromatin-binding protein. This was investigated through chromatin immunoprecipitation and high-throughput DNA sequencing (ChIP-seq). The significant HEX 70a-DNA binding sites were mainly localized at the intergenic, promoter and intronic regions. HEX 70a targeted DNA stretches mapped to the genomic regions encompassing genes with relevant functional attributes. Several HEX 70a targeted genes were associated with H3K27ac or/and H3K27me3, known as active and repressive histone marks. Brain and fat body tissues shared a fraction of the HEX 70 targeted genes, and tissue-specific targets were also detected. The presence of overrepresented DNA motifs in the binding sites is consistent with specific HEX 70a-chromatin association. In addition, a search for HEX 70a targets in RNA-seq public libraries of fat bodies from nurses and foragers revealed differentially expressed targets displaying hex 70a-correlated developmental expression, thus supporting a regulatory activity for HEX 70a. Our results support the premise that HEX 70a is a moonlighting protein that binds chromatin and has roles in the brain and fat body cell nuclei, apart from its canonical role as a storage protein.

H3K4me1 Modification Functions in Caste Differentiation in Honey Bees

Honey bees are important species for the study of epigenetics. Female honey bee larvae with the same genotype can develop into phenotypically distinct organisms (sterile workers and fertile queens) depending on conditions such as diet. Previous studies have shown that DNA methylation and histone modification can establish distinct gene expression patterns, leading to caste differentiation. It is unclear whether the histone methylation modification H3K4me1 can also impact caste differentiation. In this study, we analyzed genome-wide H3K4me1 modifications in both queen and worker larvae and found that H3K4me1 marks are more abundant in worker larvae than in queen larvae at both the second and fourth instars, and many genes associated with caste differentiation are differentially methylated. Notably, caste-specific H3K4me1 in promoter regions can direct worker development. Thus, our results suggest that H3K4me1 modification may act as an important regulatory factor in the establishment and maintenance of caste-specific transcriptional programs in honey bees; however, the potential influence of other epigenetic modifications cannot be excluded.

Identification of runs of homozygosity in Western honey bees (Apis mellifera) using whole-genome sequencing data

Runs of homozygosity (ROH) are continuous homozygous segments that arise through the transmission of haplotypes that are identical by descent. The length and distribution of ROH segments provide insights into the genetic diversity of populations and can be associated with selection signatures. Here, we analyzed reconstructed whole-genome queen genotypes, from a pool-seq data experiment including 265 Western honeybee colonies from Apis mellifera mellifera and Apis mellifera carnica. Integrating individual ROH patterns and admixture levels in a dynamic population network visualization allowed us to ascertain major differences between the two subspecies. Within A. m. mellifera, we identified well-defined substructures according to the genetic origin of the queens. Despite the current applied conservation efforts, we pinpointed 79 admixed queens. Genomic inbreeding (F ROH) strongly varied within and between the identified subpopulations. Conserved A. m. mellifera from Switzerland had the highest mean F ROH (3.39%), while queens originating from a conservation area in France, which were also highly admixed, showed significantly lower F ROH (0.45%). The majority of A. m. carnica queens were also highly admixed, except 12 purebred queens with a mean F ROH of 2.33%. Within the breed-specific ROH islands, we identified 14 coding genes for A. m. mellifera and five for A. m. carnica, respectively. Local adaption of A. m. mellifera could be suggested by the identification of genes involved in the response to ultraviolet light (Crh-BP, Uvop) and body size (Hex70a, Hex70b), while the A. m. carnica specific genes Cpr3 and Cpr4 are most likely associated with the lighter striping pattern, a morphological phenotype expected in this subspecies. We demonstrated that queen genotypes derived from pooled workers are useful tool to unravel the population dynamics in A. mellifera and provide fundamental information to conserve native honey bees.

Silencing of ATP Synthase β Impairs Egg Development in the Leafhopper Scaphoideus titanus, Vector of the Phytoplasma Associated with Grapevine Flavescence Dorée

Scaphoideus titanus (Hemiptera: Cicadellidae) is the natural vector of Flavescence dorée phytoplasma, a quarantine pest of grapevine with severe impact on European viticulture. RNA interference (RNAi) machinery components are present in S. titanus transcriptome and injection of ATP synthase β dsRNAs into adults caused gene silencing, starting three days post injection (dpi) up to 20 dpi, leading to decrease cognate protein. Silencing of this gene in the closely related leafhopper Euscelidiusvariegatus previously showed female sterility and lack of mature eggs in ovaries. Here, alteration of developing egg morphology in S. titanus ovaries as well as overexpression of hexamerin transcript (amino acid storage protein) and cathepsin L protein (lysosome proteinase) were observed in dsATP-injected females. To evaluate RNAi-specificity, E.variegatus was used as dsRNA-receiving model-species. Different doses of two sets of dsRNA-constructs targeting distinct portions of ATP synthase β gene of both species induced silencing, lack of egg development, and female sterility in E. variegatus, indicating that off-target effects must be evaluated case by case. The effectiveness of RNAi in S. titanus provides a powerful tool for functional genomics of this non-model species and paves the way toward RNAi-based strategies to limit vector population, despite several technical and regulatory constraints that still need to be overcome to allow open field application.

miRNA-34 and miRNA-210 target hexamerin genes enhancing their differential expression during early brain development of honeybee (Apis mellifera) castes

During the honeybee larval stage, queens develop larger brains than workers, with morphological differentiation appearing at the fourth larval phase (L4), just after a boost in nutritional difference both prospective females experience. The molecular promoters of this caste-specific brain development are already ongoing in previous larval phases. Transcriptomic analyses revealed a set of differentially expressed genes in the L3 brains of queens and workers, which represents the early molecular response to differential feeding females receive during larval development. Three genes of this set, hex70b, hex70c and hex110, are more highly transcribed in the brain of workers than in queens. The microRNAs miR-34, miR-210 and miR-317 are in higher levels in the queens' brain at the same phase of larval development. Here, we tested the hypothesis that the brain of workers expresses higher levels of hexamerins than that of queens during key phases of larval development and that this differential hexamerin genes expression is further enhanced by the repressing activity of miR-34, miR-210 and miR-317. Our transcriptional analyses showed that hex70b, hex70c and hex110 genes are differentially expressed in the brain of L3 and L4 larval phases of honeybee queens and workers. In silico reconstructed miRNA-mRNA interaction networks were validated using luciferase assays, which showed miR-34 and miR-210 negatively regulate hex70b and hex110 genes by directly and redundantly binding their 3'UTR (untranslated region) sequences. Taken together, our results suggest that miR-34 and miR-210 act together promoting differential brain development in honeybee castes by downregulating the expression of the putative antineurogenic hexamerin genes hex70b and hex110.

Prospects in Connecting Genetic Variation to Variation in Fertility in Male Bees

Bees are economically and ecologically important pollinating species. Managed and native bee species face increasing pressures from human-created stressors such as habitat loss, pesticide use, and introduced pathogens. There has been increasing attention towards how each of these factors impacts fertility, especially sperm production and maintenance in males. Here, we turn our attention towards another important factor impacting phenotypic variation: genetics. Using honey bees as a model, we explore the current understanding of how genetic variation within and between populations contributes to variation in sperm production, sperm maintenance, and insemination success among males. We conclude with perspectives and future directions in the study of male fertility in honey bees and non-Apis pollinators more broadly, which still remain largely understudied.

Caste-specific gene expression underlying the differential adult brain development in the honeybee Apis mellifera

Apis mellifera adult workers feature more developed key brain regions than queens, which allows them to cope with the broad range of duties they need to perform in a colony. However, at the end of larval development, the brain of queens is largely more developed than that of workers. Major morphogenetic changes take place after metamorphosis that shift caste-specific brain development. Here, we tested the hypothesis that this phenomenon is hormonally governed and involves differential gene expression. Our molecular screening approach revealed a set of differentially expressed genes in Pp (first pharate-adult phase) brains between castes mainly coding for tissue remodelling and energy-converting proteins (e.g. hex 70a and ATPsynβ). An in-depth qPCR analysis of the transcriptional behaviour during pupal and pharate-adult developmental stage in both castes and in response to artificially augmented hormone titres of 18 genes/variants revealed that: i. subtle differences in hormone titres between castes might be responsible for the differential expression of the EcR and insulin/insulin-like signalling (IIS) pathway genes; ii. the morphogenetic activity of the IIS in brain development must be mediated by ILP-2, iii. which together with the tum, mnb and caspase system, can constitute the molecular effectors of the caste-specific opposing brain developmental trajectories.

Gene expression during larval caste determination and differentiation in intermediately eusocial bumblebees, and a comparative analysis with advanced eusocial honeybees

The queen‐worker caste system of eusocial insects represents a prime example of developmental polyphenism (environmentally‐induced phenotypic polymorphism) and is intrinsic to the evolution of advanced eusociality. However, the comparative molecular basis of larval caste determination and subsequent differentiation in the eusocial Hymenoptera remains poorly known. To address this issue within bees, we profiled caste‐associated gene expression in female larvae of the intermediately eusocial bumblebee Bombus terrestris. In B. terrestris, female larvae experience a queen‐dependent period during which their caste fate as adults is determined followed by a nutrition‐sensitive period also potentially affecting caste fate but for which the evidence is weaker. We used mRNA‐seq and qRT‐PCR validation to isolate genes differentially expressed between each caste pathway in larvae at developmental stages before and after each of these periods. We show that differences in gene expression between caste pathways are small in totipotent larvae, then peak after the queen‐dependent period. Relatively few novel (i.e., taxonomically‐restricted) genes were differentially expressed between castes, though novel genes were significantly enriched in late‐instar larvae in the worker pathway. We compared sets of caste‐associated genes in B. terrestris with those reported from the advanced eusocial honeybee, Apis mellifera, and found significant but relatively low levels of overlap of gene lists between the two species. These results suggest both the existence of low numbers of shared toolkit genes and substantial divergence in caste‐associated genes between Bombus and the advanced eusocial Apis since their last common eusocial ancestor.

Protein Stores Regulate When Reproductive Displays Begin in the Male Caribbean Fruit Fly

Many animals exhibit reproductive behavior that requires expenditure of valuable nutrients. In males of many species, competitive energetically demanding displays and the development of sexual ornaments require prior accumulation of nutrient stores. Males must coordinate nutrient stores with ornament development and reproductive displays or they risk depleting their resources mid-development or mid-display, reducing their chance of mating. Males may use nutrient stores to regulate their reproductive behavior. Amino acid reserves may be important for reproduction, but the roles of amino acid stores in initiating maturation and reproductive behavior are less studied than fat stores. Insects store amino acids as hexamerin storage proteins. Many fly species use a specific hexamerin, larval serum protein 2 (LSP-2), as both a juvenile storage medium and to store protein consumed after adult eclosion. Protein stored as LSP-2 has previously been suggested to regulate reproduction in females, but no role has been proposed for LSP-2 in regulating male maturation. We use males of the Caribbean fruit fly, Anastrepha suspensa, a species with nutrient-intensive male sexual displays to test whether LSP-2 stores regulate male reproductive displays. We fed adult A. suspensa males a diet with or without protein, then assayed these males for lsp-2 transcript abundance via qRT-PCR, LSP-2 protein abundance via Western blot, and reproductive display behavior via observation. We found that adult males with ad libitum dietary protein had greater lsp-2 transcript and protein abundance, earlier sexual display behavior, and were more likely to exhibit sexual display behavior than protein-deprived adult males. We show that lsp-2 knockdown via RNAi decreases the proportion of males exhibiting reproductive displays, particularly early in the onset of reproductive behavior. Our results suggest circulating LSP-2 protein stores regulate reproductive behavior in A. suspensa males, consistent with protein stores modulating reproduction in males with expensive reproductive strategies. Our results are consistent with hexamerin storage proteins performing dual roles of protein storage and protein signaling. Our work also has substantial practical applications because tephritid flies are a pest group and the timing and expression of male reproductive displays in this group are important for control efforts using the sterile insect technique.

A novel hexamerin with an unexpected contribution to the prophenoloxidase activation system of the Chinese oak silkworm, Antheraea pernyi

Hexamerin was originally identified as a storage protein but later confirmed to be involved in many physiological processes. In the present study, we cloned and characterized a novel hexamerin complementary DNA sequence from the Chinese oak silkworm, Antheraea pernyi (Ap-hexamerin), which shows high homology with reported insect methionine-rich hexamerins. The tissue distribution and time course of expression demonstrated that Ap-hexamerin was predominantly synthesized in the fat body and the expression level was significantly increased in response to the microbial challenge, suggesting the relevance of Ap-hexamerin to immune responses. In further immune functional studies, Ap-hexamerin was confirmed to take part in the upregulation of prophenoloxidase (PPO) activation in A. pernyi haemolymph triggered by pathogen-associated molecular patterns (PAMPs). Additional molecular interaction analysis revealed that Ap-hexamerin is capable of binding the PAMPs used in the phenoloxidase assay, suggesting hexamerin in A. pernyi may positively regulate haemolymph PPO activation, acting as a pattern recognition protein.

Developmental expression and evolution of hexamerin and haemocyanin from Folsomia candida (Collembola)

Haemocyanins constitute a group of copper-containing respiratory proteins, and hexamerins were derived from hexapod haemocyanin but lost the ability to transport oxygen and serve as storage proteins. Although hexamerins have been reported in most insect species, none of them has been identified in Collembola, one of the most primitive hexapod lineages, thereby preventing us from exploring relevant evolutionary scenarios regarding the origin and evolution of hexamerins in hexapods. Here we report on collembolan hexamerins for the first time, and investigated the temporal expression profiles of hexamerin and haemocyanin in the collembolan Folsomia candida. Haemocyanin was expressed over the entire life cycle, with higher expression at the embryonic stage than at other stages, whereas hexamerin expression was restricted to embryos, unlike insect hexamerins, which are generally expressed from larval to adult stages. A phylogenetic analysis and molecular clock estimation suggested that all investigated hexapod hexamerins have a single and ancient origin (~423 Ma), coincident with the rise of atmospheric oxygen levels in the Silurian-Devonian period, indicating a physiological link between molecular evolution and Palaeozoic oxygen changes.

Caste-based differential transcriptional expression of hexamerins in response to a juvenile hormone analog in the red imported fire ant (Solenopsis invicta)

The reproductive ground plan hypothesis proposes that gene networks regulating foraging behavior and reproductive female physiology in social insects emerged from ancestral gene and endocrine factor networks. Expression of storage proteins such as vitellogenins and hexamerins is an example of this co-option. Hexamerins, through their role modulating juvenile hormone availability, are involved in caste determination in termites. The genome of the fire ant (Solenopsis invicta) encodes four hexamerin genes, hexamerin-like (LOC105192919, hereafter called hexamerin 1), hexamerin (LOC105204474, hereafter called hexamerin 2), arylphorin subunit alpha-like, and arylphorin subunit beta. In this study, a phylogenetic analysis of the S. invicta hexamerins determined that each predicted protein clustered with one of the orthologous Apis mellifera hexamerins. Gene expression analyses by RT-qPCR revealed differential expression of the hexamerins between queens and workers, and between specific task-allocated workers (nurses and foragers). Queens and nurses had significantly higher expression of all genes when compared to foragers. Hexamerin 1 was expressed at higher levels in queens, while hexamerin 2 and arylphorin subunit beta were expressed at significantly higher levels in nurses. Arylphorin subunit alpha-like showed no significant difference in expression between virgin queens and nurses. Additionally, we analyzed the relationship between the expression of hexamerin genes and S-hydroprene, a juvenile hormone analog. Significant changes in hexamerin expression were recorded in nurses, virgin queens, and foragers 12 h after application of the analog. Hexamerin 1 and arylphorin subunit alpha-like expression were significantly lower after analog application in virgin queens. In foragers, hexamerin 2 and arylphorin subunit beta were significantly lower after analog application, while in nurses expression of all genes were significantly lower after analog application. Our results suggest that in S. invicta hexamerin genes could be associated with reproductive division of labor and task-allocation of workers.

Comparative transcriptomics of social insect queen pheromones

Queen pheromones are chemical signals that mediate reproductive division of labor in eusocial animals. Remarkably, queen pheromones are composed of identical or chemically similar compounds in some ants, wasps and bees, even though these taxa diverged >150MYA and evolved queens and workers independently. Here, we measure the transcriptomic consequences of experimental exposure to queen pheromones in workers from two ant and two bee species (genera: Lasius, Apis, Bombus), and test whether they are similar across species. Queen pheromone exposure affected transcription and splicing at many loci. Many genes responded consistently in multiple species, and the set of pheromone-sensitive genes was enriched for functions relating to lipid biosynthesis and transport, olfaction, production of cuticle, oogenesis, and histone (de)acetylation. Pheromone-sensitive genes tend to be evolutionarily ancient, positively selected, peripheral in the gene coexpression network, hypomethylated, and caste-specific in their expression. Our results reveal how queen pheromones achieve their effects, and suggest that ants and bees use similar genetic modules to achieve reproductive division of labor.

Brain gene expression analyses in virgin and mated queens of fire ants reveal mating-independent and socially regulated changes

Transcriptomes of dissected brains from virgin alate and dealate mated queens from polygyne fire ants (Solenopsis invicta) were analyzed and compared. Thirteen genes were upregulated in mated queen brain, and nine were downregulated. While many of the regulated genes were either uncharacterized or noncoding RNAs, those annotated genes included two hexamerin proteins, astakine neuropeptide, serine proteases, and serine protease inhibitors. We found that for select differentially expressed genes in the brain, changes in gene expression were most likely driven by the changes in physiological state (i.e., age, nutritional status, or dominance rank) or in social environment (released from influence of primer pheromone). This was concluded because virgins that dealated after being separated from mated queens showed similar patterns of gene expression in the brain as those of mated queens for hexamerin 1, astakine, and XR_850909. Abaecin (XR_850725), however, appears upregulated only after mating. Therefore, our findings contribute to distinguish how specific gene networks, especially those influenced by queen primer pheromone, are regulated in queen ants. Additionally, to identify brain signaling pathways, we mined the fire ant genome and compiled a list of G-protein-coupled receptors (GPCRs). The expression level of GPCRs and other genes in the "genetic toolkit" in the brains of virgin alates and mated dealate queens is reported.

The Neuropeptide Corazonin Controls Social Behavior and Caste Identity in Ants

Social insects are emerging models to study how gene regulation affects behavior because their colonies comprise individuals with the same genomes but greatly different behavioral repertoires. To investigate the molecular mechanisms that activate distinct behaviors in different castes, we exploit a natural behavioral plasticity in Harpegnathos saltator, where adult workers can transition to a reproductive, queen-like state called gamergate. Analysis of brain transcriptomes during the transition reveals that corazonin, a neuropeptide homologous to the vertebrate gonadotropin-releasing hormone, is downregulated as workers become gamergates. Corazonin is also preferentially expressed in workers and/or foragers from other social insect species. Injection of corazonin in transitioning Harpegnathos individuals suppresses expression of vitellogenin in the brain and stimulates worker-like hunting behaviors, while inhibiting gamergate behaviors, such as dueling and egg deposition. We propose that corazonin is a central regulator of caste identity and behavior in social insects.

Gut symbiotic bacteria stimulate insect growth and egg production by modulating hexamerin and vitellogenin gene expression

Recent studies have suggested that gut symbionts modulate insect development and reproduction. However, the mechanisms by which gut symbionts modulate host physiologies and the molecules involved in these changes are unclear. To address these questions, we prepared three different groups of the insect Riptortus pedestris: Burkholderia gut symbiont-colonized (Sym) insects, Burkholderia-non-colonized (Apo) insects, and Burkholderia-depleted (Sym^Burk-) insects, which were fed tetracycline. When the hemolymph proteins of three insects were analyzed by SDS-PAGE, the hexamerin-α, hexamerin-β and vitellogenin-1 proteins of Sym-adults were highly expressed compared to those of Apo- and Sym^Burk--insects. To investigate the expression patterns of these three genes during insect development, we measured the transcriptional levels of these genes. The hexamerin-β gene was specifically expressed at all nymphal stages, and its expression was detected 4-5 days earlier in Sym-insect nymphs than that in Apo- and Sym^Burk--insects. However, the hexamerin-α and vitellogenin-1 genes were only expressed in adult females, and they were also detected 6-7 days earlier and were 2-fold higher in Sym-adult females than those in the other insects. Depletion of hexamerin-β by RNA interference in 2nd instar Sym-nymphs delayed adult emergence, whereas hexamerin-α and vitellogenin-1 RNA interference in 5th instar nymphs caused loss of color of the eggs of Sym-insects. These results demonstrate that the Burkholderia gut symbiont modulates host development and egg production by regulating production of these three hemolymph storage proteins.

Differential gene expression underlying ovarian phenotype determination in honey bee, Apis mellifera L., caste development

Adult honey bee queens and workers drastically differ in ovary size. This adult ovary phenotype difference becomes established during the final larval instar, when massive programmed cell death leads to the degeneration of 95-99% of the ovariole anlagen in workers. The higher juvenile hormone (JH) levels in queen larvae protect the ovaries against such degeneration. To gain insights into the molecular architecture underlying this divergence critical for adult caste fate and worker sterility, we performed a microarray analysis on fourth and early fifth instar queen and worker ovaries. For the fourth instar we found nine differentially expressed genes (DEGs) with log₂FC > 1.0, but this number increased to 56 in early fifth-instar ovaries. We selected 15 DEGs for quantitative PCR (RT-qPCR) analysis. Nine differed significantly by the variables caste and/or development. Interestingly, genes with enzyme functions were higher expressed in workers, while those related to transcription and signaling had higher transcript levels in queens. For the RT-qPCR confirmed genes we analyzed their response to JH. This revealed a significant up-regulation for two genes, a short chain dehydrogenase reductase (sdr) and a heat shock protein 90 (hsp90). Five other genes, including hsp60 and hexamerin 70b (hex70b), were significantly down-regulated by JH. The sdr gene had previously come up as differentially expressed in other transcriptome analyses on honey bee larvae and heat shock proteins are frequently involved in insect hormone responses, this making them interesting candidates for further functional assays.

20-hydroxyecdysone mediates fat body arylphorin regulation during development of rice moth, Corcyra cephalonica

Arylphorin hexamerins are one of the major insect storage proteins involved in diverse functions during metamorphosis. However, their regulation during development is not elucidated so far. In the present study, we documented 20-hydroxyecdysone (20E)-mediated regulation of arylphorin expression in the fat body of the stored grain pest, Corcyra cephalonica. Based on the differential developmental expression and 20E-induced transcriptional as well as translational level alterations of arylphorin, we isolated the 5' upstream region of the gene to analyze regulatory motifs. Promoter motif analysis revealed the presence of ecdysone response element (ERE). Transient transfection studies showed the functionality of the ERE. Enzyme mobility shift experiments with radiolabelled, cold and mutated probes indicate ERE-nuclear factor binding. This study is the first to report transcriptional regulation of arylphorins by 20E in lepdopteran insect species.

Candidate biomarkers for mosquito age-grading identified by label-free quantitative analysis of protein expression in Aedes albopictus females

We applied a "shotgun" approach based on nanoliquid chromatography-high resolution mass spectrometry associated to label free quantification (LFQ) to identify proteins varying with age, independently from the physiological state, in Aedes albopictus, a mosquito species which in the last decades invaded temperate regions in North America and Europe, creating concerns for associated high nuisance and risk of arbovirus transmission. The combined "shotgun" and LFQ approach was shown to be highly suitable to simultaneously compare several biological samples, as needed in a study aimed to analyze different age-groups and physiological states of adult mosquito females. The results obtained represent the first wide-scale analysis of protein expression in Ae. albopictus females: >1000 and 665 proteins were identified from few micrograms of crude protein extracts of mosquito heads and thoraxes, respectively. Six of these proteins were shown to significantly vary from 2- to 16-day-old females, independently from their physiological state (i.e. virgin, mated, host-seeking, blood-fed, and gravid).

BIOLOGICAL SIGNIFICANCE

Mosquito-borne diseases, such as malaria, dengue and other arboviroses, are a persistent cause of global mortality and morbidity, affecting hundreds of thousands of people. Billions of people living in tropical areas are at risk of being bitten every day by an infective mosquito female and the spread of tropical species such as Aedes albopictus to temperate areas is creating alarm in the northern hemisphere. Mosquito longevity is a critical factor affecting mosquito-borne pathogen transmission cycles and the mosquito capacity to transmit pathogens. However, large scale analyses of the age structure of mosquito field populations is hampered by the lack of optimal age-grading approaches. Our findings open new perspectives for the development of reliable, simple and cheap protein-based assays to age-grade Ae. albopictus females and, most likely, other mosquito species of higher medical relevance, such as the main dengue vector, Aedes aegypti, and the major Afrotropical malaria vectors. These assays would greatly contribute to epidemiological studies aimed at defining the actual vectorial capacity of a given mosquito species. Moreover, they would be very valuable in assessing the effectiveness of mosquito control interventions based on the relative ratio between young and old individuals before and after the intervention.

Expression and evolution of hexamerins from the tobacco hornworm, Manduca sexta, and other Lepidoptera

Hexamerins are large hemolymph-proteins that accumulate during the late larval stages of insects. Hexamerins have emerged from hemocyanin, but have lost the ability to bind oxygen. Hexamerins are mainly considered as storage proteins for non-feeding stages, but may also have other functions, e.g. in cuticle formation, transport and immune response. The genome of the hornworm Manduca sexta harbors six hexamerin genes. Two of them code for arylphorins (Msex2.01690, Msex2.15504) and two genes correspond to a methionine-rich hexamerin (Msex2.10735) and a moderately methionine-rich hexamerin (Msex2.01694), respectively. Two other genes do not correspond to any known hexamerin and distantly resemble the arylphorins (Msex2.01691, Msex2.01693). Five of the six hexamerin genes are clustered within ∼45 kb on scaffold 00023, which shows conserved synteny in various lepidopteran genomes. The methionine-rich hexamerin gene is located at a distinct site. M. sexta and other Lepidoptera have lost the riboflavin-binding hexamerin. With the exception of Msex2.01691, which displays low mRNA levels throughout the life cycle, all hexamerins are most highly expressed during pre-wandering phase of the 5th larval instar of M. sexta, supporting their role as storage proteins. Notably, Msex2.01691 is most highly expressed in the brain, suggesting a divergent function. Phylogenetic analyses showed that hexamerin evolution basically follows insect systematics. Lepidoptera display an unparalleled diversity of hexamerins, which exceeds that of other hexapod orders. In contrast to previous analyses, the lepidopteran hexamerins were found monophyletic. Five distinct types of hexamerins have been identified in this order, which differ in terms of amino acid composition and evolutionary history: i. the arylphorins, which are rich in aromatic amino acids (∼20% phenylalanine and tyrosine), ii. the distantly related arylphorin-like hexamerins, iii. the methionine-rich hexamerins, iv. the moderately methionine rich hexamerins, and v. the riboflavin-binding hexamerins.

Identification of novel allergen in edible insect, Gryllus bimaculatus and its cross-reactivity with Macrobrachium spp. allergens

Edible insects have recently been promoted as a source of protein and have a high nutrition value. Identification of allergens and cross-reactivity between Macrobrachium spp. and the field cricket (Gryllus bimaculatus) is necessary for food safety control and to assist in the diagnosis and therapy of allergy symptoms. Denaturing polyacrylamide gel electrophoresis (SDS-PAGE) was used to separate proteins. Allergens were determined and identified by IgE-immunoblotting with pooled sera from prawn-allergic patients (n=16) and LC-MS/MS. Arginine kinase (AK) and glyceraldehyde 3-phosphate dehydrogenase (GAPDH) were determined as the important allergens in muscle of Macrobrachium rosenbergii whereas, hemocyanin (HC) was identified as an allergen in Macrobrachium spp. The allergens in Macrobrachium lanchesteri were identified as AK and HC. In addition, hexamerin1B (HEX1B) was identified as a novel and specific allergen in G. bimaculatus. The important allergen in G. bimaculatus and Macrobrachium spp. is AK and was found to cross-react between both species.

Transcriptome analysis reveals signature of adaptation to landscape fragmentation

We characterize allelic and gene expression variation between populations of the Glanville fritillary butterfly (Melitaea cinxia) from two fragmented and two continuous landscapes in northern Europe. The populations exhibit significant differences in their life history traits, e.g. butterflies from fragmented landscapes have higher flight metabolic rate and dispersal rate in the field, and higher larval growth rate, than butterflies from continuous landscapes. In fragmented landscapes, local populations are small and have a high risk of local extinction, and hence the long-term persistence at the landscape level is based on frequent re-colonization of vacant habitat patches, which is predicted to select for increased dispersal rate. Using RNA-seq data and a common garden experiment, we found that a large number of genes (1,841) were differentially expressed between the landscape types. Hexamerin genes, the expression of which has previously been shown to have high heritability and which correlate strongly with larval development time in the Glanville fritillary, had higher expression in fragmented than continuous landscapes. Genes that were more highly expressed in butterflies from newly-established than old local populations within a fragmented landscape were also more highly expressed, at the landscape level, in fragmented than continuous landscapes. This result suggests that recurrent extinctions and re-colonizations in fragmented landscapes select a for specific expression profile. Genes that were significantly up-regulated following an experimental flight treatment had higher basal expression in fragmented landscapes, indicating that these butterflies are genetically primed for frequent flight. Active flight causes oxidative stress, but butterflies from fragmented landscapes were more tolerant of hypoxia. We conclude that differences in gene expression between the landscape types reflect genomic adaptations to landscape fragmentation.

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.

Culex quinquefasciatus storage proteins

Insect storage proteins accumulate at high levels during larval development of holometabolous insects. During metamorphosis they are degraded, supplying energy and amino acids for the completion of adult development. The genome of Culex quinquefasciatus contains eleven storage protein-coding genes. Their transcripts are more abundant in larvae than in pupae and in adults. In fact, only four of these genes are transcribed in adults, two of which in blood-fed adult females but not in adult males. Transcripts corresponding to all Cx. quinquefasciatus storage proteins were detected by RT-PCR, while mass spectrometric analysis of larval and pupal proteins identified all storage proteins with the exception of one encoded by Cq LSP1.8. Our results indicate that the identified Cx. quinquefasciatus storage protein-coding genes are candidates for identifying regulatory sequences for the development of molecular tools for vector control.

KONAGAbase: a genomic and transcriptomic database for the diamondback moth, Plutella xylostella

Background

The diamondback moth (DBM), Plutella xylostella, is one of the most harmful insect pests for crucifer crops worldwide. DBM has rapidly evolved high resistance to most conventional insecticides such as pyrethroids, organophosphates, fipronil, spinosad, Bacillus thuringiensis, and diamides. Therefore, it is important to develop genomic and transcriptomic DBM resources for analysis of genes related to insecticide resistance, both to clarify the mechanism of resistance of DBM and to facilitate the development of insecticides with a novel mode of action for more effective and environmentally less harmful insecticide rotation. To contribute to this goal, we developed KONAGAbase, a genomic and transcriptomic database for DBM (KONAGA is the Japanese word for DBM).

Description

KONAGAbase provides (1) transcriptomic sequences of 37,340 ESTs/mRNAs and 147,370 RNA-seq contigs which were clustered and assembled into 84,570 unigenes (30,695 contigs, 50,548 pseudo singletons, and 3,327 singletons); and (2) genomic sequences of 88,530 WGS contigs with 246,244 degenerate contigs and 106,455 singletons from which 6,310 de novo identified repeat sequences and 34,890 predicted gene-coding sequences were extracted. The unigenes and predicted gene-coding sequences were clustered and 32,800 representative sequences were extracted as a comprehensive putative gene set. These sequences were annotated with BLAST descriptions, Gene Ontology (GO) terms, and Pfam descriptions, respectively. KONAGAbase contains rich graphical user interface (GUI)-based web interfaces for easy and efficient searching, browsing, and downloading sequences and annotation data. Five useful search interfaces consisting of BLAST search, keyword search, BLAST result-based search, GO tree-based search, and genome browser are provided. KONAGAbase is publicly available from our website (http://dbm.dna.affrc.go.jp/px/) through standard web browsers.

Conclusions

KONAGAbase provides DBM comprehensive transcriptomic and draft genomic sequences with useful annotation information with easy-to-use web interfaces, which helps researchers to efficiently search for target sequences such as insect resistance-related genes. KONAGAbase will be continuously updated and additional genomic/transcriptomic resources and analysis tools will be provided for further efficient analysis of the mechanism of insecticide resistance and the development of effective insecticides with a novel mode of action for DBM.

Transcriptome analysis of the Asian honey bee Apis cerana cerana

Background

The Eastern hive honey bee, Apis cerana cerana is a native and widely bred honey bee species in China. Molecular biology research about this honey bee species is scarce, and genomic information for A. c. cerana is not currently available. Transcriptome and expression profiling data for this species are therefore important resources needed to better understand the biological mechanisms of A. c. cerana. In this study, we obtained the transcriptome information of A. c. cerana by RNA-sequencing and compared gene expression differences between queens and workers of A. c. cerana by digital gene expression (DGE) analysis.

Results

Using high-throughput Illumina RNA sequencing we obtained 51,581,510 clean reads corresponding to 4.64 Gb total nucleotides from a single run. These reads were assembled into 46,999 unigenes with a mean length of 676 bp. Based on a sequence similarity search against the five public databases (NR, Swissport, GO, COG, KEGG) with a cut-off E-value of 10⁻⁵ using BLASTX, a total of 24,630 unigenes were annotated with gene descriptions, gene ontology terms, or metabolic pathways. Using these transcriptome data as references we analyzed the gene expression differences between the queens and workers of A. c. cerana using a tag-based digital gene expression method. We obtained 5.96 and 5.66 million clean tags from the queen and worker samples, respectively. A total of 414 genes were differentially expressed between them, with 189 up-regulated and 225 down-regulated in queens.

Conclusions

Our transcriptome data provide a comprehensive sequence resource for future A. c. cerana study, establishing an important public information platform for functional genomic studies in A. c. cerana. Furthermore, the DGE data provide comprehensive gene expression information for the queens and workers, which will facilitate our understanding of the molecular mechanisms of the different physiological aspects of the two castes.

Nuclear Immunolocalization of Hexamerins in the Fat Body of Metamorphosing Honey Bees

Hexamerins are storage proteins with primordial functions in insect metamorphosis. They are actively secreted by the larval fat body and stored in the hemolymph. During metamorphosis, they return to the fat body to be processed. For decades, these proteins were thought to exclusively function as an amino acid source for tissue reconstruction during the non-feeding pupal and pharate adult stages and, in some species, for egg production. Recently, new findings have linked the hexamerins to caste polyphenism and gonad development in social insects. To explore the roles of hexamerins during the honey bee metamorphosis, we used specific antibodies in expression analysis by western blot, in situ immunolocalization by confocal laser-scanning microscopy and in vivo injections to lower their endogenous levels. Our expression analysis highlighted the changing expression patterns in the fat body and hemolymph during development, which is consistent with the temporal dynamics of hexamerin secretion, storage and depletion. Confocal microscopy showed hexamerin expression in the cytoplasm of both types of fat body cells, trophocytes and oenocytes. Notably, hexamerin foci were also found in the nuclei of these cells, thus confirming our western blot analysis of fat body nuclear-enriched fractions. We also observed that the decrease in soluble hexamerins in antibody-treated pharate adults led to a precocious adult ecdysis, perhaps in response to the lack (or decrease) in hexamerin-derived amino acids. Taken together, these findings indicate that hexamerins have other functions in addition to their well-established role as amino acid sources for development.

Potential costs of bacterial infection on storage protein gene expression and reproduction in queenless Apis mellifera worker bees on distinct dietary regimes

Insects are able to combat infection by initiating an efficient immune response that involves synthesizing antimicrobial peptides and a range of other defense molecules. These responses may be costly to the organism, resulting in it exploiting endogenous resources to maintain homeostasis or support defense to the detriment of other physiological needs. We used queenless worker bees on distinct dietary regimes that may alter hemolymph protein storage and ovary activation to investigate the physiological costs of infection with Serratia marcescens. The expression of the genes encoding the storage proteins vitellogenin and hexamerin 70a, the vitellogenin receptor, and vasa (which has a putative role in reproduction), was impaired in the infected bees. This impairment was mainly evident in the bees fed beebread, which caused significantly higher expression of these genes than did royal jelly or syrup, and this was confirmed at the vitellogenin and hexamerin 70a protein levels. Beebread was also the only diet that promoted ovary activation in the queenless bees, but this activation was significantly impaired by the infection. The expression of the genes encoding the storage proteins apolipophorins-I and -III and the lipophorin receptor was not altered by infection regardless the diet provided to the bees. Similarly, the storage of apolipophorin-I in the hemolymph was only slightly impaired by the infection, independently of the supplied diet. Taken together these results indicate that, infection demands a physiological cost from the transcription of specific protein storage-related genes and from the reproductive capacity.