Ecdysteroid-dependent expression of the tweedle and peroxidase genes during adult cuticle formation in the honey bee, Apis mellifera

Genome-wide annotation and comparative analysis revealed conserved cuticular protein evolution among non-biting midges with varied environmental adaptability

Chironomidae, non-biting midges, a diverse and abundant insect group in global aquatic ecosystems, represent an exceptional model for investigating genetic adaptability mechanisms in aquatic insects due to their extensive species diversity and resilience to various environmental conditions. The cuticle in insects acts as the primary defense against ecological pressures. Cuticular Proteins (CPs) determine cuticle characteristics, facilitating adaptation to diverse challenges. However, systematic annotation of CP genes has only been conducted for one Chironomidae species, Propsilocerus akamusi, by our team. In this study, we expanded this annotation by identifying CP genes in eight additional Chironomidae species, covering all Chironomidae species with available genome data. We identified a total of 889 CP genes, neatly categorized into nine CP families: 215 CPR RR1 genes, 272 CPR RR2 genes, 23 CPR RR3 genes, 21 CPF genes, 16 CPLCA genes, 19 CPLCG genes, 28 CPLCP genes, 77 CPAP genes, and 37 Tweedle genes. Subsequently, we conducted a comprehensive phylogenetic analysis of CPs within the Chironomidae family. This expanded annotation of CP genes across diverse Chironomidae species significantly contributes to our understanding of their remarkable adaptability.

The cuticle proteome of a planktonic crustacean

The cuticles of arthropods provide an interface between the organism and its environment. Thus, the cuticle's structure influences how the organism responds to and interacts with its surroundings. Here, we used label-free quantification proteomics to provide a proteome of the moulted cuticle of the aquatic crustacean Daphnia magna, which has long been a prominent subject of studies on ecology, evolution, and developmental biology. We detected a total of 278 high-confidence proteins. Using protein sequence domain and functional enrichment analyses, we identified chitin-binding structural proteins and chitin-modifying enzymes as the most abundant protein groups in the cuticle proteome. Structural cuticular protein families showed a similar distribution to those found in other arthropods and indicated proteins responsible for the soft and flexible structure of the Daphnia cuticle. Finally, cuticle protein genes were also clustered as tandem gene arrays in the D. magna genome. The cuticle proteome presented here will be a valuable resource to the Daphnia research community, informing genome annotations and investigations on diverse topics such as the genetic basis of interactions with predators and parasites.

Enzymes in "Green" Synthetic Chemistry: Laccase and Lipase

Enzymes play an important role in numerous natural processes and are increasingly being utilized as environmentally friendly substitutes and alternatives to many common catalysts. Their essential advantages are high catalytic efficiency, substrate specificity, minimal formation of byproducts, and low energy demand. All of these benefits make enzymes highly desirable targets of academic research and industrial development. This review has the modest aim of briefly overviewing the classification, mechanism of action, basic kinetics and reaction condition effects that are common across all six enzyme classes. Special attention is devoted to immobilization strategies as the main tools to improve the resistance to environmental stress factors (temperature, pH and solvents) and prolong the catalytic lifecycle of these biocatalysts. The advantages and drawbacks of methods such as macromolecular crosslinking, solid scaffold carriers, entrapment, and surface modification (covalent and physical) are discussed and illustrated using numerous examples. Among the hundreds and possibly thousands of known and recently discovered enzymes, hydrolases and oxidoreductases are distinguished by their relative availability, stability, and wide use in synthetic applications, which include pharmaceutics, food and beverage treatments, environmental clean-up, and polymerizations. Two representatives of those groups-laccase (an oxidoreductase) and lipase (a hydrolase)-are discussed at length, including their structure, catalytic mechanism, and diverse usage. Objective representation of the current status and emerging trends are provided in the main conclusions.

Tweedle gene family of Drosophila suzukii (Matsumura) larva enhances the basal tolerance to cold and hypoxia

BACKGROUND

Drosophila suzukii (Matsumura) is considered a quarantine pest in the A2 list because it causes serious infection and huge economic losses. Cold and controlled atmosphere treatments have been used to control immature stage pests in fresh fruits. Herein, the basal tolerance response of D. suzukii egg, larva and pupa to cold and hypoxia stress were studied, and underlying transcriptome mechanisms in the larva were pinpointed.

RESULTS

The third instar was more tolerant than 12-h-old egg and 8-day-old pupa when treated at 3 °C + 1% O2 for 7 days, with 34.00% ± 5.22% larval survival. Hypoxia influenced the effect of cold treatment on D. suzukii. Larval survival decreased at 3 °C + 1% O2 , but increased at 0 °C + 1% O2 . Survival increased with temperature between 0 and 5 °C + 1% O2 , but decreased significantly at 25 °C + 1% O2 . RNA-sequencing results showed that the Tweedle (Twdl) family was upregulated and uniquely enriched in larvae treated at 3 °C + 1% O2 . In addition, RNA interference-mediated silencing of a key Twdl gene reduced the survival rate after cold and hypoxia treatment.

CONCLUSION

Hypoxia was able to influence the effect of cold treatment on the survival of D. suzukii positively or negatively. Structural constituents of the chitin-based cuticle, in particular Twdl genes, body morphogenesis, and ATP synthesis-coupled proton transport were involved in the tolerance to cold and hypoxia. In future, the Twdl gene could be used as a nanocarrier delivering RNA pesticides to control D. suzukii in the field and so prevent its worldwide spread. © 2023 Society of Chemical Industry.

The inhibition of ecdysone signal pathway was the key of pyriproxyfen poisoning for silkworm, Bombyx mori

Pyriproxyfen is a juvenile hormone-like pesticide. Once intake occurs, it leads to a series of poisoning characters consequences in silkworm, Bombyx mori (ID: 7091, Lepidoptera), such as non- cocooning, non-pupation, production of low-active eggs, and extended stages. However, the poisoning mechanism is still unclear. Here, silkworms were fed mulberry leaves soaked with different pyriproxyfen concentrations, and the heads were dissected for transcriptome analysis, while the hemolymph was used for determinations of ecdysone and juvenile hormone titers. As a result, after conjoint analysis of 3 feeding groups and a control group, 555 differentially expressed genes (DEGs) were obtained, which were mainly involved in hormone metabolism, glycometabolism and protein metabolism. Meanwhile, 119 genes were significantly correlated with the pyriproxyfen concentrations, and they were mainly involved in drug metabolism and glycometabolism. The ecdysone titers in several feeding groups were significantly lower than those of the control group, while juvenile hormone was not detected in all groups, including the control and feeding groups. Correspondingly, due to activation of the juvenile hormone signaling pathway by pyriproxyfen, key genes in the ecdysone synthesis pathway were downregulated, and a large number of downstream genes were up- or downregulated. In addition, nearly all genes in the detoxification pathway were upregulated. These results suggested that, affected by the juvenile hormone signaling pathway, ecdysone titers decreased and further affected a series of downstream processes, and this was the key reason for pyriproxyfen poisoning in silkworm, B. mori, which could lay a foundation for the study of pyriproxyfen resistance in silkworm.

Genome-wide annotation of cuticular protein genes in non-biting midge Propsilocerus akamusi and transcriptome analysis of their response to heavy metal pollution

The insect cuticle is a sophisticated chitin-protein extracellular structure for mutable functions. The cuticles varied their structures and properties in different species, and the same species but in different regions or at different stages, to fill the requirements of different functions. The alteration of cuticle structures may also be induced due to challenges by some environmental crises, such as pollution exposures. The physical properties of the cuticle were determined by the cuticle proteins (CPs) they contain. The cuticle proteins are large protein groups in all insects, which are commonly divided into different families according to their conserved protein sequence motifs. Although Chironomidae is an abundant and universal insect in global aquatic ecosystems and a popular model for aquatic toxicology, no systematic annotation of CPs was done for any species in Chironomidae before. In this work, we annotated the CP genes of Propsilocerus akamusi, the most abundant Chironomidae species in Asia. A total of 160 CP genes were identified, and 97 of them could be well classified into eight CP families: 76 CPR genes can be subdivided into three groups (further divided into three subgroups: 36 RR1 genes, 37 RR2 genes, and 3 RR3 genes), 2 CPF genes, 3 CPLCA genes, 1 CPLCG gene, 8 CPAP genes, and 3 Tweedle genes. Additionally, we analyzed the response of P. akamusi CP genes at expression level to Cu exposure, which is related to the high heavy metal tolerance and the earlier onset of pupariation in heavy metal polluted water.

Cuticular protein genes involve heat acclimation of insect larvae under global warming

Cuticular proteins (CPs) play important roles in insect growth and development. However, it is unknown whether CPs are related to heat tolerance. Cnaphalocrocis medinalis, a serious pest of rice, occurs in summer and exhibits strong adaptability to high temperature, but the underlying mechanism is unclear. Here, the role of CP genes in heat acclimation was studied. Heat tolerance of the heat-acclimated larvae was significantly stronger than the unacclimated larvae. The cuticular protein content in the heat-acclimated larvae was higher than that of the unacclimated larvae. 191 presumed CP genes of C. medinalis (CmCPs) were identified. Expression patterns of 14 CmCPs were different between the heat acclimated (S39) and unacclimated (S27) larvae under heat stress. CmCPs were specifically expressed in epidermis and the head except CmCPR20 mainly expressed in Malpighian tubules. CmCPR20 was upregulated in S39 while downregulated in S27, but CmTweedle1 and CmCPG1 were upregulated in S27 and downregulated in S39. RNAi CmTweedle1 or CmCPG1 remarkably decreased heat tolerance and cuticular protein content of the heat-acclimated larvae but not the unacclimated larvae. RNAi CmCPR20 decreased heat tolerance and cuticular protein content of the unacclimated larvae but not the heat-acclimated larvae. CmTweedle1 and CmCPG1 genes involve heat acclimation of C. medinalis.

Tweedle proteins form extracellular two-dimensional structures defining body and cell shape in Drosophila melanogaster

Tissue function and shape rely on the organization of the extracellular matrix (ECM) produced by the respective cells. Our understanding of the underlying molecular mechanisms is limited. Here, we show that extracellular Tweedle (Twdl) proteins in the fruit fly Drosophila melanogaster form two adjacent two-dimensional sheets underneath the cuticle surface and above a distinct layer of dityrosinylated and probably elastic proteins enwrapping the whole body. Dominant mutations in twdl genes cause ectopic spherical aggregation of Twdl proteins that recruit dityrosinylated proteins at their periphery within lower cuticle regions. These aggregates perturb parallel ridges at the surface of epidermal cells that have been demonstrated to be crucial for body shaping. In one scenario, hence, this disorientation of epidermal ridges may explain the squatty phenotype of Twdl mutant larvae. In an alternative scenario, this phenotype may be due to the depletion of the dityrosinylated and elastic layer, and the consequent weakening of cuticle resistance against the internal hydrostatic pressure. According to Barlow's formula describing the distribution of internal pressure forces in pipes in dependence of pipe wall material properties, it follows that this reduction in turn causes lateral expansion at the expense of the antero-posterior elongation of the body.

Identification and temporal expression profiles of cuticular proteins in the endoparasitoid wasp, Microplitis mediator

Recently, parasitoid wasp species Microplitis mediator has evoked increasing research attention due to its possible use in the control of Lepidoptera insects. Because insect development involves changes in cuticle composition, identification and expression analysis of M. mediator cuticular proteins may clarify the mechanisms involved in parasite development processes. We found 70 cuticular proteins from the M. mediator transcriptome and divided them into seven distinct families. Expression profiling indicated that most of these cuticular protein genes have expression peaks specific for one particular developmental stage of M. mediator. Eggs and pupae have the highest number of transcriptionally active cuticular protein genes (47 and 52 respectively). Only 12 of these genes maintained high expression activity during late larval development. Functional analysis of two larval proteins, MmCPR3 and MmCPR14, suggested their important role in the proper organization of the cuticle layers of larvae. During M. mediator larval development, normal cuticle formation can be supported by a limited number of cuticular proteins.

A Novel miRNA-hlo-miR-2-Serves as a Regulatory Factor That Controls Molting Events by Targeting CPR1 in Haemaphysalis longicornis Nymphs

Successful completion of the molting process requires new epidermal growth and ecdysis of the old cuticle in Haemaphysalis longicornis (H. longicornis). MicroRNAs (miRNAs) participate in the development of organisms by inhibiting the expression of their target mRNAs. In this study, a novel tick-specific miRNA was identified and denoted hlo-miR-2 that serves as a novel regulator of molting events in H. longicornis nymphs by targeting a cuticular protein. The full length of this cuticular protein was first obtained and named it CPR1. A qRT-PCR analysis showed that hlo-miR-2 and CPR1 exhibit significant tissue and temporal specificity and that their transcription levels are negatively correlated during the molting process. CPR1, as a direct target of hlo-miR-2, was identified by a luciferase reporter assay in vitro. Agomir treatment indicated that the overexpression of hlo-miR-2 significantly reduced the protein expression level of CPR1, decreased the molting rate and delayed the molting time point in H. longicornis nymphs. RNA interference (RNAi) experiments demonstrated that CPR1 was significantly associated with the molting process in H. longicornis nymphs. Phenotypic rescue experiments convincingly showed that hlo-miR-2 participated in molting events by targeting CPR1 in H. longicornis nymphs. In summary, we present evidence demonstrating that miRNAs constitute a novel important regulator of molting events in addition to hormones. The described functional evidence implicating CPR1 in molting events contributes to an improved understanding of the distinct functions of the CPR family in ticks and will aid the development of a promising application of cuticular protein RNAi in tick control.

Prepupae and pupae transcriptomic characterization of Trichogramma chilonis

The egg parasitoid, Trichogramma chilonis, has significant control effects on agriculture and forestry pests and is widely employed in southern China for the biological control of lepidopteran pests. In this study, transcriptomic analysis was used to gain a clear understanding of the molecular changes in prepupae and pupae of T. chilonis. A total of 16.88 Gb of clean data were obtained and finally assembled into 43,136 unigenes, 18,880 of which were annotated. After FPKM standardization, 117 and 838 specific expression genes were found in prepupae and pupae, respectively. There were 3129 differentially expressed genes between prepupae and pupae. Compared to pupae, 806 genes were up-regulated and 2323 were down-regulated in prepupae. Background on the T. chilonis transcriptome, the enriched GO function and KEGG pathway analysis of DEGs were considered. As indicated by GO classification, up-regulated genes were mainly involved in chitin metabolism, cell adhesion and endocytic, while most down-regulated genes were involved in synthesis of cell components, ion transport and biological regulation. KEGG enrichment analysis showed that 458 DEGs were enriched in 94 metabolic pathways. DEGs involved in nucleotide replication and transcription, substance metabolism, insect hormone biosynthesis, cell growth and death, reproductive metabolism, circadian rhythms and signal transduction pathways were up-regulated or down-regulated to different degrees, indicating that these genes played important roles during the process of metamorphosis in T. chilonis. This study provides a rich data source for the future study of T. chilonis molecular and biological mechanisms. A large number of genes related to metamorphosis were found based on comparison analysis between prepupae and pupae transcriptomes. This study lays a good foundation for in-depth study of gene transcription and regulation mechanisms during T. chilonis metamorphosis.

Inhibition of fatty acid desaturation impairs cuticle differentiation in Drosophila melanogaster

Previously, we showed that inhibition of the activity of fatty acid desaturases (Desat) perturbs signalling of the developmental timing hormone ecdysone in the fruit fly Drosophila melanogaster. To understand the impact of this effect on cuticle differentiation, a process regulated by ecdysone, we analysed the cuticle of D. melanogaster larvae fed with the Desat inhibitor CA10556. In these larvae, the expression of most of the key cuticle genes is normal or slightly elevated at day one of CA10556 feeding. As an exception, expression of twdlM coding for a yet uncharacterised cuticle protein is completely suppressed. The cuticle of these larvae appears to be normal at the morphological level. However, these animals are sensitive to desiccation, a trait that according to our data, among others, may be associated with reduced TwdlM amounts. At day two of CA10556 feeding, expression of most of the cuticle genes tested including twdlM is suppressed. Expression of cpr47Eb coding for a chitin-binding protein is, by contrast, highly elevated suggesting that Cpr47Eb participates at a specific compensation program. Overall, the cuticle of these larvae is thinner than the cuticle of control larvae. Taken together, lipid desaturation is necessary for a coordinated deployment of a normal cuticle differentiation program.

Exploring integument transcriptomes, cuticle ultrastructure, and cuticular hydrocarbons profiles in eusocial and solitary bee species displaying heterochronic adult cuticle maturation

Differences in the timing of exoskeleton melanization and sclerotization are evident when comparing eusocial and solitary bees. This cuticular maturation heterochrony may be associated with life style, considering that eusocial bees remain protected inside the nest for many days after emergence, while the solitary bees immediately start outside activities. To address this issue, we characterized gene expression using large-scale RNA sequencing (RNA-seq), and quantified cuticular hydrocarbon (CHC) through gas chromatography-mass spectrometry in comparative studies of the integument (cuticle plus its underlying epidermis) of two eusocial and a solitary bee species. In addition, we used transmission electron microscopy (TEM) for studying the developing cuticle of these and other three bee species also differing in life style. We found 13,200, 55,209 and 30,161 transcript types in the integument of the eusocial Apis mellifera and Frieseomelitta varia, and the solitary Centris analis, respectively. In general, structural cuticle proteins and chitin-related genes were upregulated in pharate-adults and newly-emerged bees whereas transcripts for odorant binding proteins, cytochrome P450 and antioxidant proteins were overrepresented in foragers. Consistent with our hypothesis, a distance correlation analysis based on the differentially expressed genes suggested delayed cuticle maturation in A. mellifera in comparison to the solitary bee. However, this was not confirmed in the comparison with F. varia. The expression profiles of 27 of 119 genes displaying functional attributes related to cuticle formation/differentiation were positively correlated between A. mellifera and F. varia, and negatively or non-correlated with C. analis, suggesting roles in cuticular maturation heterochrony. However, we also found transcript profiles positively correlated between each one of the eusocial species and C. analis. Gene co-expression networks greatly differed between the bee species, but we identified common gene interactions exclusively between the eusocial species. Except for F. varia, the TEM analysis is consistent with cuticle development timing adapted to the social or solitary life style. In support to our hypothesis, the absolute quantities of n-alkanes and unsaturated CHCs were significantly higher in foragers than in the earlier developmental phases of the eusocial bees, but did not discriminate newly-emerged from foragers in C. analis. By highlighting differences in integument gene expression, cuticle ultrastructure, and CHC profiles between eusocial and solitary bees, our data provided insights into the process of heterochronic cuticle maturation associated to the way of life.

Hormonal control and target genes of ftz-f1 expression in the honeybee Apis mellifera: a positive loop linking juvenile hormone, ftz-f1, and vitellogenin

Ftz-f1 is an orphan member of the nuclear hormone receptor superfamily. A 20-hydroxyecdysone pulse allows ftz-f1 gene expression, which then regulates the activity of downstream genes involved in major developmental progression events. In honeybees, the expression of genes like vitellogenin (vg), prophenoloxidase and juvenile hormone-esterase during late pharate-adult development is known to be hormonally controlled in both queens and workers by increasing juvenile hormone (JH) titres in the presence of declining levels of ecdysteroids. Since Ftz-f1 is known for mediating intracellular JH signalling, we hypothesized that ftz-f1 could mediate JH action during the pharate-adult development of honeybees, thus controlling the expression of these genes. Here, we show that ftz-f1 has caste-specific transcription profiles during this developmental period, with a peak coinciding with the increase in JH titre, and that its expression is upregulated by JH and downregulated by ecdysteroids. RNAi-mediated knock down of ftz-f1 showed that the expression of genes essential for adult development (e.g. vg and cuticular genes) depends on ftz-f1 expression. Finally, a double-repressor hypothesis-inspired vg gene knock-down experiment suggests the existence of a positive molecular loop between JH, ftz-f1 and vg.

Developmental changes in gene expression and enzyme activities of anabolic and catabolic enzymes for storage carbohydrates in the honeybee, Apis mellifera

Glycogen and trehalose are important sources of energy in insects. The expression of genes encoding the key metabolic enzymes-glycogen synthase (GS), glycogen phosphorylase (GP), trehalose-6-phosphate synthase (TPS-1), soluble trehalase (Tre-1) and membrane-bound trehalase (Tre-2)-was analyzed in 12 developmental stages of Apis mellifera worker brood. The content of GS and GP proteins, TPS activity, total trehalase activity, and the activity of Tre-1 and Tre-2 were determined. Transcript quantity was not always correlated with the content of the encoded GS or GP protein. The correlation was higher for GS (r = 0.797) than GP (r = 0.651). The expression of the glycogen synthase gene (gs) and the glycogen phosphorylase gene (gp) was high in 4- and 7-day-old larvae and in pupae, excluding the last pupal stage. The expression of the tps-1 gene was highest in the mid-pupal stage and contributed to higher enzyme activity in that stage. The expression of the tre-1 gene was higher than the expression of the tre-2 gene throughout development. In newly hatched workers, the expression of genes encoding catabolic enzymes of both carbohydrates, gp and tre-1, was higher than the expression of genes encoding anabolic enzymes. The results of this study suggest that sugar metabolism genes have somewhat different control mechanisms during larval development and metamorphosis.

Comparative antennal transcriptome of Apis cerana cerana from four developmental stages

Apis cerana cerana, an important endemic honey bee species in China, possesses valuable characteristics such as a sensitive olfactory system, good foraging ability, and strong resistance to parasitic mites. Here, we performed transcriptome sequencing of the antenna, the major chemosensory organ of the bee, using an Illumina sequencer, to identify typical differentially expressed genes (DEGs) in adult worker bees of different ages, namely, T1 (1 day); T2 (10 days); T3 (15 days); and T4 (25 days). Surprisingly, the expression levels of DEGs changed significantly between the T1 period and the other three periods. All the DEGs were classified into 26 expression profiles by trend analysis. Selected trend clusters were analyzed, and valuable information on gene expression patterns was obtained. We found that the expression levels of genes encoding cuticle proteins declined after eclosion, while those of immunity-related genes increased. In addition, genes encoding venom proteins and major royal jelly proteins were enriched at the T2 stage; small heat shock proteins showed significantly higher expression at the T3 stage; and some metabolism-related genes were more highly expressed at the T4 stage. The DEGs identified in this study may serve as a valuable resource for the characterization of expression patterns of antennal genes in A. cerana cerana. Furthermore, this study provides insights into the relationship between labor division in social bees and gene function.

Differentially expressed genes in the head of the 2nd instar pre-molting larvae of the nm2 mutant of the silkworm, Bombyx mori

Molting is an important physiological process in the larval stage of Bombyx mori and is controlled by various hormones and peptides. The silkworm mutant that exhibits the phenotype of non-molting in the 2^nd instar (nm2) is incapable of molting in the 2^nd instar and dies after seven or more days. The ecdysone titer in the nm2 mutant is lower than that in the wildtype, and the mutant can be rescued by feeding with 20E and cholesterol. The results of positional cloning indicated that structural alteration of BmCPG10 is responsible for the phenotype of the nm2 mutant. To explore the possible relationship between BmCPG10 and the ecdysone titer as well as the genes affected by BmCPG10, digital gene expression (DGE) profile analysis was conducted in the nm2 mutant, with the wildtype strain C603 serving as the control. The results revealed 1727 differentially expressed genes, among which 651 genes were upregulated and 1076 were downregulated in nm2. BLASTGO analysis showed that these differentially expressed genes were involved in various biological processes, cellular components and molecular functions. KEGG analysis indicated an enrichment of these differentially expressed genes in 240 pathways, including metabolic pathways, pancreatic secretion, protein digestion and absorption, fat digestion and absorption and glycerolipid metabolism. To verify the accuracy of the DGE results, quantitative reverse transcription PCR (qRT-PCR) was performed, focusing on key genes in several related pathways, and the results were highly consistent with the DGE results. Our findings indicated significant differences in cuticular protein genes, ecdysone biosynthesis genes and ecdysone-related nuclear receptors genes, but no significant difference in juvenile hormone and chitin biosynthesis genes was detected. Our research findings lay the foundation for further research on the formation mechanism of the nm2 mutant.

RNAi-Mediated Functional Analysis of Bursicon Genes Related to Adult Cuticle Formation and Tanning in the Honeybee, Apis mellifera

Bursicon is a heterodimeric neurohormone that acts through a G protein-coupled receptor named rickets (rk), thus inducing an increase in cAMP and the activation of tyrosine hydroxylase, the rate-limiting enzyme in the cuticular tanning pathway. In insects, the role of bursicon in the post-ecdysial tanning of the adult cuticle and wing expansion is well characterized. Here we investigated the roles of the genes encoding the bursicon subunits during the adult cuticle development in the honeybee, Apis mellifera. RNAi-mediated knockdown of AmBurs α and AmBurs β bursicon genes prevented the complete formation and tanning (melanization/sclerotization) of the adult cuticle. A thinner, much less tanned cuticle was produced, and ecdysis toward adult stage was impaired. Consistent with these results, the knockdown of bursicon transcripts also interfered in the expression of genes encoding its receptor, AmRk, structural cuticular proteins, and enzymes in the melanization/sclerotization pathway, thus evidencing roles for bursicon in adult cuticle formation and tanning. Moreover, the expression of AmBurs α, AmBurs β and AmRk is contingent on the declining ecdysteroid titer that triggers the onset of adult cuticle synthesis and deposition. The search for transcripts of AmBurs α, AmBurs β and candidate targets in RNA-seq libraries prepared with brains and integuments strengthened our data on transcript quantification through RT-qPCR. Together, our results support our premise that bursicon has roles in adult cuticle formation and tanning, and are in agreement with other recent studies pointing for roles during the pharate-adult stage, in addition to the classical post-ecdysial ones.

PROTEINS OF THE INTEGUMENTARY SYSTEM OF THE HONEYBEE, Apis mellifera

The integument of insects and other arthropods is composed of an inner basal lamina coated by the epidermis, which secretes the bulk of the outer integument layer, the cuticle. The genome sequencing of several insect species has allowed predicting classes of proteins integrating the cuticle. However, only a small proportion of them, as well as other proteins in the integumentary system, have been validated. Using two-dimensional gel electrophoresis coupled with mass spectrometry, we identified 45 different proteins in a total of 112 selected gel spots derived from thoracic integument samples of developing honeybee workers, including 14 cuticular proteins (AmelCPR 3, AmelCPR 12, AmelCPR 16, AmelCPR 27, apidermin 2, apidermin 3, endocuticle structural glycoprotein SgAbd-8-like, LOC100577363, LOC408365, LOC413679, LOC725454, LOC100576916, LOC725838, and peritrophin 3-C analogous). Gene ontology functional analysis revealed that the higher proportions of the identified proteins have molecular functions related to catalytic and structural molecule activities, are involved in metabolic biological processes, and pertain to the protein class of structural or cytoskeletal proteins and hydrolases. It is noteworthy that 26.7% of the identified proteins, including five cuticular proteins, were revealed as protein species resulting from allelic isoforms or derived from posttranslational modifications. Also, 66.7% of the identified cuticular proteins were expressed in more than one developmental phase, thus indicating that they are part of the larval, pupal, and adult cuticle. Our data provide experimental support for predicted honeybee gene products and new information on proteins expressed in the developing integument.

miR-71 and miR-263 Jointly Regulate Target Genes Chitin synthase and Chitinase to Control Locust Molting

Chitin synthase and chitinase play crucial roles in chitin biosynthesis and degradation during insect molting. Silencing of Dicer-1 results in reduced levels of mature miRNAs and severely blocks molting in the migratory locust. However, the regulatory mechanism of miRNAs in the molting process of locusts has remained elusive. In this study, we found that in chitin metabolism, two crucial enzymes, chitin synthase (CHS) and chitinase (CHT) were regulated by miR-71 and miR-263 during nymph molting. The coding sequence of CHS1 and the 3'-untranslated region of CHT10 contain functional binding sites for miR-71 and miR-263, respectively. miR-71/miR-263 displayed cellular co-localization with their target genes in epidermal cells and directly interacted with CHS1 and CHT10 in the locust integument, respectively. Injections of miR-71 and miR-263 agomirs suppressed the expression of CHS1 and CHT10, which consequently altered chitin production of new and old cuticles and resulted in a molting-defective phenotype in locusts. Unexpectedly, reduced expression of miR-71 and miR-263 increased CHS1 and CHT10 mRNA expression and led to molting defects similar to those induced by miRNA delivery. This study reveals a novel function and balancing modulation pattern of two miRNAs in chitin biosynthesis and degradation, and it provides insight into the underlying molecular mechanisms of the molting process in locusts.

Cuticular protein LmTwdl1 is involved in molt development of the migratory locust

The cuticle, an essential structure for insects, is produced from cuticular proteins and chitin via a series of biochemical reactions. Tweedle genes are important members of the cuticular protein family and have four conserved motifs binding to chitin. Tweedle family genes have been found to play a profound effect on cuticle development. Here, we report that the cuticular protein gene LmTwdl1 of Locusta migratoria belongs to the Tweedle family. In situ hybridization showed that LmTwdl1 is localized to epidermal cells of the cuticle. The expression patterns of LmTwdl1 showed low expression in the cuticle during the early and middle stages of the fifth-instar nymphs; in contrast, its expression rapidly increased in the late stages of fifth-instar nymphs. We performed RNA interference to examine the function of LmTwdl1 in locusts. Silencing of LmTwdl1 resulted in high mortality during the molting process before the next stage. Also, the epicuticle of nymphs failed to molt, tended to be thinner and the arrangement of chitin in the procuticle appeared to be disordered compare to the control group. These results demonstrate that LmTwdl1 plays a critical role in molting, which contributes to a better understanding of the distinct functions of the Tweedle family in locusts.

Insights into the dynamics of hind leg development in honey bee (Apis mellifera L.) queen and worker larvae - A morphology/differential gene expression analysis

Phenotypic plasticity is a hallmark of the caste systems of social insects, expressed in their life history and morphological traits. These are best studied in bees. In their co-evolution with angiosperm plants, the females of corbiculate bees have acquired a specialized structure on their hind legs for collecting pollen. In the highly eusocial bees (Apini and Meliponini), this structure is however only present in workers and absent in queens. By means of histological sections and cell proliferation analysis we followed the developmental dynamics of the hind legs of queens and workers in the fourth and fifth larval instars. In parallel, we generated subtractive cDNA libraries for hind leg discs of queen and worker larvae by means of a Representational Difference Analysis (RDA). From the total of 135 unique sequences we selected 19 for RT-qPCR analysis, where six of these were confirmed as differing significantly in their expression between the two castes in the larval spinning stage. The development of complex structures such as the bees' hind legs, requires diverse patterning mechanisms and signaling modules, as indicated by the set of differentially expressed genes related with cell adhesion and signaling pathways.

Tweedle cuticular protein BmCPT1 is involved in innate immunity by participating in recognition of Escherichia coli

Bombyx mori, a lepidopteran insect, is one of the earliest models for pattern recognition of Gram-negative bacteria, which may induce the IMD pathway for production of antibacterial peptides. So far, several recognition proteins have been reported in B. mori. However, the connection between pattern recognition of Gram negative bacteria and activation of BmRelish1, a transcription factor controlled by the IMD pathway remains largely unknown. In the present study, we identify BmCPT1, a cuticle protein bearing a Tweedle domain. Its gene expression is co-regulated by NF-kappaB and juvenile hormone signals. BmCPT1 is induced by Escherichia coli in fat bodies and hemocytes, but is constitutively expressed in the epidermis. In vitro binding assays indicate that BmCPT1 protein recognizes and binds to E. coli peptidoglycan. Post-transcriptionally modified BmCPT1 in the hemolymph binds to E. coli cells through interactions with peptidoglycan recognition protein-5 (BmPGRP5) and lipopolysaccharide binding protein (BmLBP). Transgenic overexpression of BmCPT1 causes the upregulated expression of BmRelish1 and clear induction of two gloverin genes. Therefore, BmCPT1 may work along with BmPGRP-S5 and BmLBP to recognize E. coli in the hemolymph and indirectly activate BmRelish1 to induce antimicrobial peptide synthesis.

Developmental regulation of ecdysone receptor (EcR) and EcR-controlled gene expression during pharate-adult development of honeybees (Apis mellifera)

Major developmental transitions in multicellular organisms are driven by steroid hormones. In insects, these, together with juvenile hormone (JH), control development, metamorphosis, reproduction and aging, and are also suggested to play an important role in caste differentiation of social insects. Here, we aimed to determine how EcR transcription and ecdysteroid titers are related during honeybee postembryonic development and what may actually be the role of EcR in caste development of this social insect. In addition, we expected that knocking-down EcR gene expression would give us information on the participation of the respective protein in regulating downstream targets of EcR. We found that in Apis mellifera females, EcR-A is the predominantly expressed variant in postembryonic development, while EcR-B transcript levels are higher in embryos, indicating an early developmental switch in EcR function. During larval and pupal stages, EcR-B expression levels are very low, while EcR-A transcripts are more variable and abundant in workers compared to queens. Strikingly, these transcript levels are opposite to the ecdysteroid titer profile. 20-hydroxyecdysone (20E) application experiments revealed that low 20E levels induce EcR expression during development, whereas high ecdysteroid titers seem to be repressive. By means of RNAi-mediated knockdown (KD) of both EcR transcript variants we detected the differential expression of 234 poly-A⁺ transcripts encoding genes such as CYPs, MRJPs and certain hormone response genes (Kr-h1 and ftz-f1). EcR-KD also promoted the differential expression of 70 miRNAs, including highly conserved ones (e.g., miR-133 and miR-375), as well honeybee-specific ones (e.g., miR-3745 and miR-3761). Our results put in evidence a broad spectrum of EcR-controlled gene expression during postembryonic development of honeybees, revealing new facets of EcR biology in this social insect.

Gene expression in five salmon louse (Lepeophtheirus salmonis, Krøyer 1837) tissues

The Atlantic salmon, Salmo salar L, is an important species both for traditional fishery and fish farming. Many Atlantic salmon stocks have been declining and a suspected main contributor to this decline is the salmon louse (Lepeophtheirus salmonis); a parasitic copepod living off the salmonid hosts epidermal tissues and blood. Contributing to the growing body of knowledge on the molecular biology of the salmon louse we have utilized a microarray containing 11,100 salmon louse genes to study the gene expression patterns in selected tissues. This approach has yielded information about potential functions of the transcripts and tissues. Microarray analyses were preformed on subcuticular and frontal (neuronal and gland enriched tissue) tissues, as well as gut, ovary and testes of adult lice. Tissue specific transcriptomes were evident, allowing us to address main traits of functional partitioning between tissues and providing valuable insight into the biology of the louse. The results furthermore represent an important tool and resource for further experiments.

Exoskeleton formation in Apis mellifera: cuticular hydrocarbons profiles and expression of desaturase and elongase genes during pupal and adult development

Cuticular hydrocarbons (CHCs) are abundant in the superficial cuticular layer (envelope) of insects where they play roles as structural, anti-desiccation and semiochemical compounds. Many studies have investigated the CHC composition in the adult insects. However, studies on the profiles of these compounds during cuticle formation and differentiation are scarce and restrict to specific stages of a few insect species. We characterized the CHCs developmental profiles in the honeybee workers during an entire molting cycle (from pupal-to-adult ecdyses) and in mature adults (forager bees). Gas chromatography/mass spectrometry (GC/MS) analysis revealed remarkable differences in the relative quantities of CHCs, thus discriminating pupae, developing and newly-ecdysed adults, and foragers from each other. In parallel, the honeybee genome database was searched for predicted gene models using known amino acid sequences of insect enzymes catalyzing lipid desaturation (desaturases) or elongation (elongases) as queries in BLASTP analysis. The expression levels of six desaturase genes and ten elongase genes potentially involved in CHC biosynthesis were determined by reverse transcription and real time polymerase chain reaction (RT-qPCR) in the developing integument (cuticle and subjacent epidermis). Aiming to predict roles for these genes in CHC biosynthesis, the developmental profiles of CHCs and desaturase/elongase transcript levels were evaluated using Spearman correlation coefficient. This analysis pointed to differential roles for these gene products in the biosynthesis of certain CHC classes. Based on the assumption that homologous proteins may share a similar function, phylogenetic trees were reconstructed as an additional strategy to predict functions and evolutionary relationships of the honeybee desaturases and elongases. Together, these approaches highlighted the molecular complexity underlying the formation of the lesser known layer of the cuticular exoskeleton, the envelope.

Genes involved in thoracic exoskeleton formation during the pupal-to-adult molt in a social insect model, Apis mellifera

Background

The insect exoskeleton provides shape, waterproofing, and locomotion via attached somatic muscles. The exoskeleton is renewed during molting, a process regulated by ecdysteroid hormones. The holometabolous pupa transforms into an adult during the imaginal molt, when the epidermis synthe3sizes the definitive exoskeleton that then differentiates progressively. An important issue in insect development concerns how the exoskeletal regions are constructed to provide their morphological, physiological and mechanical functions. We used whole-genome oligonucleotide microarrays to screen for genes involved in exoskeletal formation in the honeybee thoracic dorsum. Our analysis included three sampling times during the pupal-to-adult molt, i.e., before, during and after the ecdysteroid-induced apolysis that triggers synthesis of the adult exoskeleton.

Results

Gene ontology annotation based on orthologous relationships with Drosophila melanogaster genes placed the honeybee differentially expressed genes (DEGs) into distinct categories of Biological Process and Molecular Function, depending on developmental time, revealing the functional elements required for adult exoskeleton formation. Of the 1,253 unique DEGs, 547 were upregulated in the thoracic dorsum after apolysis, suggesting induction by the ecdysteroid pulse. The upregulated gene set included 20 of the 47 cuticular protein (CP) genes that were previously identified in the honeybee genome, and three novel putative CP genes that do not belong to a known CP family. In situ hybridization showed that two of the novel genes were abundantly expressed in the epidermis during adult exoskeleton formation, strongly implicating them as genuine CP genes. Conserved sequence motifs identified the CP genes as members of the CPR, Tweedle, Apidermin, CPF, CPLCP1 and Analogous-to-Peritrophins families. Furthermore, 28 of the 36 muscle-related DEGs were upregulated during the de novo formation of striated fibers attached to the exoskeleton. A search for cis-regulatory motifs in the 5′-untranslated region of the DEGs revealed potential binding sites for known transcription factors. Construction of a regulatory network showed that various upregulated CP- and muscle-related genes (15 and 21 genes, respectively) share common elements, suggesting co-regulation during thoracic exoskeleton formation.

Conclusions

These findings help reveal molecular aspects of rigid thoracic exoskeleton formation during the ecdysteroid-coordinated pupal-to-adult molt in the honeybee.