MicroRNAs in Honey Bee Caste Determination

Juvenile hormone induces reproduction via miR-1175-3p in the red imported fire ant, Solenopsis invicta

Juvenile hormone (JH) acts in the regulation of caste differentiation between queens and workers (i.e., with or without reproductive capacity) during vitellin synthesis and oogenesis in social insects. However, the regulatory mechanisms have not yet been elucidated. Here, we identified a highly expressed microRNA (miRNA), miR-1175-3p, in the red imported fire ant, Solenopsis invicta. We found that miR-1175-3p is prominently present in the fat bodies and ovaries of workers. Furthermore, miR-1175-3p interacts with its target gene, broad-complex core (Br-C), in the fat bodies. By utilizing miR-1175-3p agomir, we successfully suppressed the expression of the Br-C protein in queens, resulting in reduced vitellogenin expression, fewer eggs, and poorly developed ovaries. Conversely, decreasing miR-1175-3p levels led to the increased expression of Br-C and vitellogenin in workers, triggering the "re-development" of the ovaries. Moreover, when queens were fed with JH, the expression of miR-1175-3p decreased, whereas the expression of vitellogenin-2 and vitellogenin-3 increased. Notably, the suppression of fertility in queens caused by treatment with agomir miR-1175-3p was completely rescued by the increased vitellogenin expression induced by being fed with JH. These results suggest the critical role of miR-1175-3p in JH-regulated reproduction, shedding light on the molecular mechanism underlying miRNA-mediated fecundity in social insects and providing a novel strategy for managing S. invicta.

Regulation of soldier caste differentiation by microRNAs in Formosan subterranean termite (Coptotermes formosanus Shiraki)

The soldier caste is one of the most distinguished castes inside the termite colony. The mechanism of soldier caste differentiation has mainly been studied at the transcriptional level, but the function of microRNAs (miRNAs) in soldier caste differentiation is seldom studied. In this study, the workers of Coptotermes formosanus Shiraki were treated with methoprene, a juvenile hormone analog which can induce workers to transform into soldiers. The miRNomes of the methoprene-treated workers and the controls were sequenced. Then, the differentially expressed miRNAs (DEmiRs) were corrected with the differentially expressed genes DEGs to construct the DEmiR-DEG regulatory network. Afterwards, the DEmiR-regulated DEGs were subjected to GO enrichment and KEGG enrichment analysis. A total of 1,324 miRNAs were identified, among which 116 miRNAs were screened as DEmiRs between the methoprene-treated group and the control group. A total of 4,433 DEmiR-DEG pairs were obtained. No GO term was recognized as significant in the cellular component, molecular function, or biological process categories. The KEGG enrichment analysis of the DEmiR-regulated DEGs showed that the ribosome biogenesis in eukaryotes and circadian rhythm-fly pathways were enriched. This study demonstrates that DEmiRs and DEGs form a complex network regulating soldier caste differentiation in termites.

Exploring ethylene-related genes in Cannabis sativa: implications for sexual plasticity

KEY MESSAGE

Presented here are model Yang cycle, ethylene biosynthesis and signaling pathways in Cannabis sativa. C. sativa floral transcriptomes were used to predict putative ethylene-related genes involved in sexual plasticity in the species. Sexual plasticity is a phenomenon, wherein organisms possess the ability to alter their phenotypic sex in response to environmental and physiological stimuli, without modifying their sex chromosomes. Cannabis sativa L., a medically valuable plant species, exhibits sexual plasticity when subjected to specific chemicals that influence ethylene biosynthesis and signaling. Nevertheless, the precise contribution of ethylene-related genes (ERGs) to sexual plasticity in cannabis remains unexplored. The current study employed Arabidopsis thaliana L. as a model organism to conduct gene orthology analysis and reconstruct the Yang Cycle, ethylene biosynthesis, and ethylene signaling pathways in C. sativa. Additionally, two transcriptomic datasets comprising male, female, and chemically induced male flowers were examined to identify expression patterns in ERGs associated with sexual determination and sexual plasticity. These ERGs involved in sexual plasticity were categorized into two distinct expression patterns: floral organ concordant (FOC) and unique (uERG). Furthermore, a third expression pattern, termed karyotype concordant (KC) expression, was proposed, which plays a role in sex determination. The study revealed that CsERGs associated with sexual plasticity are dispersed throughout the genome and are not limited to the sex chromosomes, indicating a widespread regulation of sexual plasticity in C. sativa.

Dynamic Evolution of Repetitive Elements and Chromatin States in Apis mellifera Subspecies

In this study, we elucidate the contribution of repetitive DNA sequences to the establishment of social structures in honeybees (Apis mellifera). Despite recent advancements in understanding the molecular mechanisms underlying the formation of honeybee castes, primarily associated with Notch signaling, the comprehensive identification of specific genomic cis-regulatory sequences remains elusive. Our objective is to characterize the repetitive landscape within the genomes of two honeybee subspecies, namely A. m. mellifera and A. m. ligustica. An observed recent burst of repeats in A. m. mellifera highlights a notable distinction between the two subspecies. After that, we transitioned to identifying differentially expressed DNA elements that may function as cis-regulatory elements. Nevertheless, the expression of these sequences showed minimal disparity in the transcriptome during caste differentiation, a pivotal process in honeybee eusocial organization. Despite this, chromatin segmentation, facilitated by ATAC-seq, ChIP-seq, and RNA-seq data, revealed a distinct chromatin state associated with repeats. Lastly, an analysis of sequence divergence among elements indicates successive changes in repeat states, correlating with their respective time of origin. Collectively, these findings propose a potential role of repeats in acquiring novel regulatory functions.

Labor division of worker ants can be controlled by insulin synthesis targeted through miR-279c-5p in Solenopsis invicta (Hymenoptera: Formicidae)

BACKGROUND

In social insects, the labor division of workers is ubiquitous and controlled by genetic and environmental factors. However, how they modulate this coordinately remains poorly understood.

RESULTS

We report miR-279c-5p participation in insulin synthesis and behavioral transition by negatively regulating Rab8A in Solenopsis invicta. Eusocial specific miR-279c-5p is age-associated and highly expressed in nurse workers, and localized in the cytoplasm of neurons, where it is partly co-localized with its target, Rab8A. We determined that miR-279c-5p agomir suppressed Rab8A expression in forager workers, consequently decreasing insulin content, resulting in the behavioral shift to 'nurse-like' behaviors, while the decrease in miR-279c-5p increased Rab8A expression and increased insulin content in nurse workers, leading to the behavioral shift to 'foraging-like' behaviors. Moreover, insulin could rescue the 'foraging behavior' induced by feeding miR-279c-5p to nurse workers. The overexpression and suppression of miR-279c-5p in vivo caused an obvious behavioral transition between foragers and nurses, and insulin synthesis was affected by miR-279c-5p by regulating the direct target Rab8A.

CONCLUSION

We first report that miR-279c-5p is a novel regulator that promotes labor division by negatively regulating the target gene Rab8A by controlling insulin production in ants. This miRNA-mediated mechanism is significant for understanding the behavioral plasticity of social insects between complex factors and potentially provides new targets for controlling red imported fire ants. © 2023 Society of Chemical Industry.

Diverse Regulatory Manners and Potential Roles of lncRNAs in the Developmental Process of Asian Honey Bee (Apis cerana) Larval Guts

Long non-coding RNAs (lncRNAs) are crucial modulators in a variety of biological processes, such as gene expression, development, and immune defense. However, little is known about the function of lncRNAs in the development of Asian honey bee (Apis cerana) larval guts. Here, on the basis of our previously obtained deep-sequencing data from the 4-, 5-, and 6-day-old larval guts of A. cerana workers (Ac4, Ac5, and Ac6 groups), an in-depth transcriptome-wide investigation was conducted to decipher the expression pattern, regulatory manners, and potential roles of lncRNAs during the developmental process of A. cerana worker larval guts, followed by the verification of the relative expression of differentially expressed lncRNAs (DElncRNAs) and the targeting relationships within a competing endogenous RNA (ceRNA) axis. In the Ac4 vs. Ac5 and Ac5 vs. Ac6 comparison groups, 527 and 498 DElncRNAs were identified, respectively, which is suggestive of the dynamic expression of lncRNAs during the developmental process of larval guts. A cis-acting analysis showed that 330 and 393 neighboring genes of the aforementioned DElncRNAs were respectively involved in 29 and 32 functional terms, such as cellular processes and metabolic processes; these neighboring genes were also respectively engaged in 246 and 246 pathways such as the Hedgehog signaling pathway and the Wnt signaling pathway. Additionally, it was found that 79 and 76 DElncRNAs as potential antisense lncRNAs may, respectively, interact with 72 and 60 sense-strand mRNAs. An investigation of competing endogenous RNA (ceRNA) networks suggested that 75 (155) DElncRNAs in the Ac4 vs. Ac5 (Ac5 vs. Ac6) comparison group could target 7 (5) DEmiRNAs and further bind to 334 (248) DEmRNAs, which can be annotated to 33 (29) functional terms and 186 (210) pathways, including 12 (16) cellular- and humoral-immune pathways (lysosome pathway, necroptosis, MAPK signaling pathway, etc.) and 11 (10) development-associated signaling pathways (Wnt, Hippo, AMPK, etc.). The RT-qPCR detection of five randomly selected DElncRNAs confirmed the reliability of the used sequencing data. Moreover, the results of a dual-luciferase reporter assay were indicative of the binding relationship between MSTRG.11294.1 and miR-6001-y and between miR-6001-y and ncbi_107992440. These results demonstrate that DElncRNAs are likely to modulate the developmental process of larval guts via the regulation of the source genes' transcription, interaction with mRNAs, and ceRNA networks. Our findings not only yield new insights into the developmental mechanism underlying A. cerana larval guts, but also provide a candidate ceRNA axis for further functional dissection.

Expression Profile, Regulatory Network, and Putative Role of microRNAs in the Developmental Process of Asian Honey Bee Larval Guts

MiRNAs, as a kind of key regulators in gene expression, play vital roles in numerous life activities from cellular proliferation and differentiation to development and immunity. However, little is known about the regulatory manner of miRNAs in the development of Asian honey bee (Apis cerana) guts. Here, on basis of our previously gained high-quality transcriptome data, transcriptome-wide identification of miRNAs in the larval guts of Apis cerana cerana was conducted, followed by investigation of the miRNAs' differential expression profile during the gut development. In addition to the regulatory network, the potential function of differentially expressed miRNAs (DEmiRNAs) was further analyzed. In total, 330, 351, and 321 miRNAs were identified in the 4-, 5-, and 6-day-old larval guts, respectively; among these, 257 miRNAs were shared, while 38, 51, and 36 ones were specifically expressed. Sequences of six miRNAs were confirmed by stem-loop RT-PCR and Sanger sequencing. Additionally, in the "Ac4 vs. Ac5" comparison group, there were seven up-regulated and eight down-regulated miRNAs; these DEmiRNAs could target 5041 mRNAs, involving a series of GO terms and KEGG pathways associated with growth and development, such as cellular process, cell part, Wnt, and Hippo. Comparatively, four up-regulated and six down-regulated miRNAs detected in the "Ac5 vs. Ac6" comparison group and the targets were associated with diverse development-related terms and pathways, including cell, organelle, Notch and Wnt. Intriguingly, it was noticed that miR-6001-y presented a continuous up-regulation trend across the developmental process of larval guts, implying that miR-6001-y may be a potential essential modulator in the development process of larval guts. Further investigation indicated that 43 targets in the "Ac4 vs. Ac5" comparison group and 31 targets in the "Ac5 vs. Ac6" comparison group were engaged in several crucial development-associated signaling pathways such as Wnt, Hippo, and Notch. Ultimately, the expression trends of five randomly selected DEmiRNAs were verified using RT-qPCR. These results demonstrated that dynamic expression and structural alteration of miRNAs were accompanied by the development of A. c. cerana larval guts, and DEmiRNAs were likely to participate in the modulation of growth as well as development of larval guts by affecting several critical pathways via regulation of the expression of target genes. Our data offer a basis for elucidating the developmental mechanism underlying Asian honey bee larval guts.

Cross-Kingdom Regulation of Plant-Derived miRNAs in Modulating Insect Development

MicroRNAs (miRNAs), a class of non-coding small RNAs, are crucial regulatory factors in plants and animals at the post-transcriptional level. These tiny molecules suppress gene expression by complementary oligonucleotide binding to sites in the target messenger. Recently, the discovery of plant-derived miRNAs with cross-kingdom abilities to regulate gene expression in insects has promoted exciting discussion, although some controversies exist regarding the modulation of insect development by plant-derived miRNAs. Here, we review current knowledge about the mechanisms of miRNA biogenesis, the roles of miRNAs in coevolution between insects and plants, the regulation of insect development by plant-derived miRNAs, the cross-kingdom transport mechanisms of plant-derived miRNAs, and cross-kingdom regulation. In addition, the controversy regarding the modulation of insect development by plant-derived miRNAs also was discussed. Our review provides new insights for understanding complex plant-insect interactions and discovering new strategies for pest management and even crop genetic improvement.

Phenotypic dimorphism between honeybee queen and worker is regulated by complicated epigenetic modifications

Phenotypic dimorphism between queens and workers is an important biological characteristic of honeybees that has been the subject of intensive research. The enormous differences in morphology, lifespan, physiology, and behavior between queens and workers are caused by a complicated set of factors. Epigenetic modifications are considered to play an important role in this process. In this study, we analyzed the differences in chromosome interactions and H3K27ac and H3K4me1 modifications between the queens and workers using high-throughput chromosome conformation capture (Hi-C) and Chromatin immunoprecipitation followed by sequencing (ChIP-Seq) technologies. We found that the queens contain more chromosome interactions and more unique H3K27ac modifications than workers; in contrast, workers have more H3K4me1 modifications than queens. Moreover, we identified Map3k15 as a potential caste gene in queen-worker differentiation. Our results suggest that chromosomal conformation and H3K27ac and H3K4me1 modifications are involved in regulating queen-worker differentiation, which reveals that the queen-worker phenotypic dimorphism is regulated by multiple epigenetic modifications.

Convergent and complementary selection shaped gains and losses of eusociality in sweat bees

Sweat bees have repeatedly gained and lost eusociality, a transition from individual to group reproduction. Here we generate chromosome-length genome assemblies for 17 species and identify genomic signatures of evolutionary trade-offs associated with transitions between social and solitary living. Both young genes and regulatory regions show enrichment for these molecular patterns. We also identify loci that show evidence of complementary signals of positive and relaxed selection linked specifically to the convergent gains and losses of eusociality in sweat bees. This includes two pleiotropic proteins that bind and transport juvenile hormone (JH)-a key regulator of insect development and reproduction. We find that one of these proteins is primarily expressed in subperineurial glial cells that form the insect blood-brain barrier and that brain levels of JH vary by sociality. Our findings are consistent with a role of JH in modulating social behaviour and suggest that eusocial evolution was facilitated by alteration of the proteins that bind and transport JH, revealing how an ancestral developmental hormone may have been co-opted during one of life's major transitions. More broadly, our results highlight how evolutionary trade-offs have structured the molecular basis of eusociality in these bees and demonstrate how both directional selection and release from constraint can shape trait evolution.

Circular RNA ame_circ_2015 Function as microRNA Sponges in Regulating Egg-Laying of Honeybees (Apis mellifera)

Honeybees (Apis mellifera) are critical to maintaining ecological balance and are important pollinators. The oviposition behavior in honeybees is important and complex. Circular RNAs (circRNAs) are found to form circRNA-miRNA crosstalk and play important roles in reproduction processes. Here, dual luciferase reporter was used to confirm the crosstalk between ame_circ_2015 and ame_miR-14-3p. Functional experiments in vitro and in vivo were performed to investigate the biological functions of ame_circ_2015 in egg-laying of queens. The results showed that ame_circ_2015 directly target ame_miR-14-3p, and the expression of ame_circ_2015 was negatively correlated with ame_miR-14-3p expression. Overexpression results showed that ame_circ_2015 promoted the number of eggs laid and knockdown of ame_circ_2015 suppressed the number of eggs laid. It demonstrates that up-regulated ame_circ_2015 promotes the number of eggs laid by sponging ame_miR-14-3p. The study will provide information towards a better understanding of circRNA-miRNA crosstalk in egg-laying in honeybees.

MicroRNAs reshape the immunity of insects in response to bacterial infection

The interaction between bacteria and insects can significantly impact a wide range of different areas because bacteria and insects are widely distributed around the globe. The bacterial-insect interactions have the potential to directly affect human health since insects are vectors for disease transmission, and their interactions can also have economic consequences. In addition, they have been linked to high mortality rates in economically important insects, resulting in substantial economic losses. MicroRNAs (miRNAs) are types of non-coding RNAs involved in regulating gene expression post-transcriptionally. The length of miRNAs ranges from 19 to 22 nucleotides. MiRNAs, in addition to their ability to exhibit dynamic expression patterns, have a diverse range of targets. This enables them to govern various physiological activities in insects, like innate immune responses. Increasing evidence suggests that miRNAs have a crucial biological role in bacterial infection by influencing immune responses and other mechanisms for resistance. This review focuses on some of the most recent and exciting discoveries made in recent years, including the correlation between the dysregulation of miRNA expression in the context of bacterial infection and the progression of the infection. Furthermore, it describes how they profoundly impact the immune responses of the host by targeting the Toll, IMD, and JNK signaling pathways. It also emphasizes the biological function of miRNAs in regulating immune responses in insects. Finally, it also discusses current knowledge gaps about the function of miRNAs in insect immunity, in addition to areas that require more research in the future.

In the battle of the disease: a transcriptomic analysis of European foulbrood-diseased larvae of the Western honey bee (Apis mellifera)

BACKGROUND

European foulbrood is a significant bacterial brood disease of Apis sp. and can cause severe and devastating damages in beekeeping operations. Nevertheless, the epidemiology of its causative agent Melissococcus plutonius has been begun to uncover but the underlying mechanisms of infection and cause of disease still is not well understood. Here, we sought to provide insight into the infection mechanism of EFB employing RNAseq in in vitro reared Apis mellifera larvae of two developmental stages to trace transcriptional changes in the course of the disease, including Paenibacillus alvei secondary infected individuals.

RESULTS

In consideration of the progressing development of the larva, we show that infected individuals incur a shift in metabolic and structural protein-encoding genes, which are involved in metabolism of crucial compounds including all branches of macronutrient metabolism, transport protein genes and most strikingly chitin and cuticle associated genes. These changes underpin the frequently observed developmental retardation in EFB disease. Further, sets of expressed genes markedly differ in different stages of infection with almost no overlap. In an earlier stage of infection, a group of regulators of the melanization response cascade and complement component-like genes, predominantly C-type lectin genes, are up-regulated while a differential expression of immune effector genes is completely missing. In contrast, late-stage infected larvae up-regulated the expression of antimicrobial peptides, lysozymes and prominent bacteria-binding haemocyte receptor genes compared to controls. While we clearly show a significant effect of infection on expressed genes, these changes may partly result from a shift in expression timing due to developmental alterations of infection. A secondary infection with P. alvei elicits a specific response with most of the M. plutonius associated differential immune effector gene expression missing and several immune pathway genes even down-regulated.

CONCLUSION

We conclude that with progressing infection diseased individuals undergo a systemic response with a change of metabolism and their activated immune defence repertoire. Moreover, larvae are capable of adjusting their response to a secondary invasion in late stage infections.

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

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

Extent and complexity of RNA processing in honey bee queen and worker caste development

The distinct honeybee (Apis mellifera) worker and queen castes have become a model for the study of genomic mechanisms of phenotypic plasticity. Here we performed a nanopore-based direct RNA sequencing with exceptionally long reads to compare the mRNA transcripts between queen and workers at three points during their larval development. We found thousands of significantly differentially expressed transcript isoforms (DEIs) between queen and worker larvae. These DEIs were formatted by a flexible splicing system. We showed that poly(A) tails participated in this caste differentiation by negatively regulating the expression of DEIs. Hundreds of isoforms uniquely expressed in either queens or workers during their larval development, and isoforms were expressed at different points in queen and worker larval development demonstrating a dynamic relationship between isoform expression and developmental mechanisms. These findings show the full complexity of RNA processing and transcript expression in honey bee phenotypic plasticity.

•

Honeybee caste differentiation has a complexity of RNA processing

•

Isoforms differentially express between queens and workers during larval development

•

Isoforms are formatted by a flexible alternative splicing system

•

Poly(A) tails are negatively correlated with isoform expression

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Social environment affects sensory gene expression in ant larvae

Social insects depend on communication to regulate social behaviour. This also applies to their larvae, which are commonly exposed to social interactions and can react to social stimulation. However, how social insect larvae sense their environment is not known. Using RNAseq, we characterized expression of sensory-related genes in larvae of the ant Formica fusca, upon exposure to two social environments: isolation without contact to other individuals, and stimulation via the presence of other developing individuals. Expression of key sensory-related genes was higher following social stimulation, and larvae expressed many of the same sensory-related genes as adult ants and larvae of other insects, including genes belonging to the major insect chemosensory gene families. Our study provides first insights into the molecular changes associated with social information perception in social insect larvae.

Regulation of membrane phospholipids during the adult life of worker honey bee

Two female castes that are genetically identical are found in honey bees: workers and queens. Adult female honey bees differ in their morphology and behaviors, but the most intriguing difference between the castes is the difference in their longevity. Queens live for years while workers live generally for weeks. The mechanisms that mediate this extraordinary difference in lifespan remain mostly unknown. Both castes share similar developmental stages and are fed liquid food (i.e. a jelly) during development. However, after emergence, workers begin to feed on pollen while queens are fed the same larval food for their entire life. Pollen has a high content of polyunsaturated fatty acids (PUFA) while royal jelly has negligible amounts. The difference in food during adult life leads to drastic changes in membrane phospholipids of female honey bees, and those changes have been proposed as mechanisms that could explain the difference in lifespan. To provide further details on those mechanisms, we characterized the membrane phospholipids of adult workers at seven different ages covering all life-history stages. Our results suggest that the majority of changes in worker membranes occur in the first four days of adult life. Shortly after emergence, workers increase their level of total phospholipids by producing phospholipids that contained saturated (SFA) and monounsaturated fatty acids (MUFA). From the second day, workers start replacing fatty acid chains from those pre-synthesized molecules with PUFA acquired from pollen. After four days, worker membranes are set and appear to be maintained for the rest of adult life, suggesting that damaged PUFA are replaced effectively. Plasmalogen phospholipids increase continuously throughout worker adult life, suggesting that plasmalogen might help to reduce lipid peroxidation in worker membranes. We postulate that the diet-induced increase in PUFA in worker membranes makes them far more prone to lipid-based oxidative damage compared to queens.

Location-Dependent DNA Methylation Signatures in a Clonal Invasive Crayfish

DNA methylation is an important epigenetic modification that has been repeatedly implied in organismal adaptation. However, many previous studies that have linked DNA methylation patterns to environmental parameters have been limited by confounding factors, such as cell-type heterogeneity and genetic variation. In this study, we analyzed DNA methylation variation in marbled crayfish, a clonal and invasive freshwater crayfish that is characterized by a largely tissue-invariant methylome and negligible genetic variation. Using a capture-based subgenome bisulfite sequencing approach that covers a small, variably methylated portion of the marbled crayfish genome, we identified specific and highly localized DNA methylation signatures for specimens from geographically and ecologically distinct wild populations. These results were replicated both biologically and technically by re-sampling at different time points and by using independent methodology. Finally, we show specific methylation signatures for laboratory animals and for laboratory animals that were reared at a lower temperature. Our results thus demonstrate the existence of context-dependent DNA methylation signatures in a clonal animal.

Thrice out of Asia and the adaptive radiation of the western honey bee

The origin of the western honey bee Apis mellifera has been intensely debated. Addressing this knowledge gap is essential for understanding the evolution and genetics of one of the world’s most important pollinators. By analyzing 251 genomes from 18 native subspecies, we found support for an Asian origin of honey bees with at least three expansions leading to African and European lineages. The adaptive radiation of honey bees involved selection on a few genomic “hotspots.” We found 145 genes with independent signatures of selection across all bee lineages, and these genes were highly associated with worker traits. Our results indicate that a core set of genes associated with worker and colony traits facilitated the adaptive radiation of honey bees across their vast distribution.

Active genic machinery for epigenetic RNA modifications in bees

Epitranscriptomics is an emerging field of investigation dedicated to the study of post-transcriptional RNA modifications. RNA methylations regulate RNA metabolism and processing, including changes in response to environmental cues. Although RNA modifications are conserved from bacteria to eukaryotes, there is little evidence of an epitranscriptomic pathway in insects. Here we identified genes related to RNA m⁶ A (N6-methyladenine) and m⁵ C (5-methylcytosine) methylation machinery in seven bee genomes (Apis mellifera, Melipona quadrifasciata, Frieseomelitta varia, Eufriesea mexicana, Bombus terrestris, Megachile rotundata and Dufourea novaeangliae). In A. mellifera, we validated the expression of methyltransferase genes and found that the global levels of m⁶ A and m⁵ C measured in the fat body and brain of adult workers differ significantly. Also, m⁶ A levels were differed significantly mainly between the fourth larval instar of queens and workers. Moreover, we found a conserved m⁵ C site in the honeybee 28S rRNA. Taken together, we confirm the existence of epitranscriptomic machinery acting in bees and open avenues for future investigations on RNA epigenetics in a wide spectrum of hymenopteran species.

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.

Transcriptome-based insights into gene networks controlling myopia prevention

Myopia (short-sightedness), usually caused by excessive elongation of the eye during development, has reached epidemic proportions worldwide. In animal systems including the chicken model, several treatments have been shown to inhibit ocular elongation and experimental myopia. Although diverse in their apparent mechanism of action, each one leads to a reduction in the rate of ocular growth. We hypothesize that a defined set of retinal molecular changes may underlie growth inhibition, irrespective of the treatment agent used. Accordingly, across five well-established but diverse methods of inhibiting myopia, significant overlap is seen in the retinal transcriptome profile (transcript levels and alternative splicing events) in chicks when analyzed by RNA-seq. Within the two major pathway networks enriched during growth inhibition, that of cell signaling and circadian entrainment, transcription factors form the largest functional grouping. Importantly, a large percentage of those genes forming the defined retinal response are downstream targets of the transcription factor EGR1 which itself shows a universal response to all five growth-inhibitory treatments. This supports EGR1's previously implicated role in ocular growth regulation. Finally, by contrasting our data with human linkage and GWAS studies on refractive error, we confirm the applicability of our study to the human condition. Together, these findings suggest that a universal set of transcriptome changes, which sit within a well-defined retinal network that cannot be bypassed, is fundamental to growth regulation, thus paving a way for designing novel targets for myopia therapies.

Diet-derived transmission of MicroRNAs from host plant into honey bee Midgut

Background

MicroRNA (miRNA) is a class of small noncoding RNAs, which targets on thousands of mRNA and thus plays important roles in many biological processes. It has been reported that miRNA has cross-species regulation functions between parasitoid-host, or plant-animal, etc. For example, several plant miRNAs enter into the honey bees and regulate gene expression. However, whether cross-species regulation function of miRNAs is a universal mechanism remains a debate question.

Results

We have evaluated transmission of miRNAs from sunflower and sedr plants into the midgut of honey bee using RNA-Seq analyses complemented with confirmation by RT-qPCR. The results showed that at least 11 plant miRNAs were found in the midgut of honey bee feeding by sunflower and sedr pollen. Among which, nine miRNAs, including miR-30d, miR-143, miR-148a, miR-21, let-7 g, miR-26a, miR-126, miR-27a, and miR-203, were shared between the sunflower- and sedr-fed honey bees, suggesting they might have essential roles in plant-insect interactions. Moreover, existence of these co-shared miRNAs presents a strong evidence to support the successful transmission of miRNAs into the midgut of the insect. In total, 121 honeybee mRNAs were predicted to be the target of these 11 plant-derived miRNAs. Interestingly, a sedr-derived miRNA, miR-206, targets on 53 honeybee genes. Kyoto Encyclopedia of Genes and Genome (KEGG) analyses showed that these target genes are significantly involved in hippo signaling pathway-fly, Wnt signaling pathway, and N-Glycan biosynthesis.

Conclusions

In summary, these results provide evidence of cross-species regulation function of miRNA between honeybee and flowering host plants, extending our understanding of the molecular interactions between plants and animals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07916-4.

Characterization and target gene analysis of microRNAs in the antennae of the parasitoid wasp Microplitis mediator

MicroRNAs (miRNAs), a class of non-coding single-strand RNA molecules encoded by endogenous genes, are about 21-24 nucleotides long and are involved in the post-transcriptional regulation of gene expression in plants and animals. Generally, the types and quantities of miRNAs in the different tissues of an organism are diverse, and these divergences may be related to their specific functions. Here we have identified 296 known miRNAs and 46 novel miRNAs in the antennae of the parasitoid wasp Microplitis mediator by high-throughput sequencing. Thirty-three miRNAs were predicted to target olfactory-associated genes, including odorant binding proteins (OBPs), chemosensory proteins, odorant receptors (ORs), ionotropic receptors (IRs) and gustatory receptors. Among these, 17 miRNAs were significantly highly expressed in the antennae, four miRNAs were highly expressed both in the antennae and head or wings, while the remaining 12 miRNAs were mainly expressed in the head, thorax, abdomen, legs and wings. Notably, miR-9a-5p and miR-2525-3p were highly expressed in male antennae, whereas miR-1000-5p and novel-miR-13 were enriched in female antennae. The 17 miRNAs highly expressed in antennae are likely to be associated with olfaction, and were predicted to target one OBP (targeted by miR-3751-3p), one IR (targeted by miR-7-5p) and 14 ORs (targeted by 15 miRNAs including miR-6-3p, miR-9a-5p, miR-9b-5p, miR-29-5p, miR-71-5p, miR-275-3p, miR-1000-5p, miR-1000-3p, miR-2525-3p, miR-6012-3p, miR-9719-3p, novel-miR-10, novel-miR-13, novel-miR-14 and novel-miR-28). These candidate olfactory-associated miRNAs are all likely to be involved in chemoreception through the regulation of chemosensory gene expression in the antennae of M. mediator.

Molecular underpinnings of the early brain developmental response to differential feeding in the honey bee Apis mellifera

Brain differential morphogenesis in females is one of the major phenotypic manifestations of caste development in honey bees. Brain diphenism appears at the fourth larval phase as a result of the differential feeding regime developing females are submitted during early phases of larval development. Here, we used a forward genetics approach to test the early brain molecular response to differential feeding leading to the brain diphenism observed at later developmental phases. Using RNA sequencing analysis, we identified 53 differentially expressed genes (DEGs) between the brains of queens and workers at the third larval phase. Since miRNAs have been suggested to play a role in caste differentiation after horizontal and vertical transmission, we tested their potential participation in regulating the DEGs. The miRNA-mRNA interaction network, including the DEGs and the royal- and worker-jelly enriched miRNA populations, revealed a subset of miRNAs potentially involved in regulating the expression of DEGs. The interaction of miR-34, miR-210, and miR-317 with Takeout, Neurotrophin-1, Forked, and Masquerade genes was experimentally confirmed using a luciferase reporter system. Taken together, our results reconstruct the regulatory network that governs the development of the early brain diphenism in honey bees.

The role of epigenetics, particularly DNA methylation, in the evolution of caste in insect societies

Eusocial insects can be defined as those that live in colonies and have distinct queens and workers. For most species, queens and workers arise from a common genome, and so caste-specific developmental trajectories must arise from epigenetic processes. In this review, we examine the epigenetic mechanisms that may be involved in the regulation of caste dimorphism. Early work on honeybees suggested that DNA methylation plays a causal role in the divergent development of queen and worker castes. This view has now been challenged by studies that did not find consistent associations between methylation and caste in honeybees and other species. Evidence for the involvement of methylation in modulating behaviour of adult workers is also inconsistent. Thus, the functional significance of DNA methylation in social insects remains equivocal. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

Organ-specific transcriptome analysis reveals differential gene expression in different castes under natural conditions in Apis cerana

Honeybees are one of the most environmentally important insects, as their pollination of various plant species contributes to the balance among different ecosystems. It has been studied extensively for their unique attribute of forming a caste society. Unlike other insects, honeybees communicate socially by secreting pheromones or by exhibiting specific patterns of motion. In the honeybee industry, the Asian honeybees (Apis cerana) and the Western honeybees (Apis mellifera) are dominant species. However, molecular research on the transcriptomes of A. cerana has not been studied as extensively as those of A. mellifera. Therefore, in this study, caste-specific transcriptional differences were analyzed, which provides a comprehensive analysis of A. cerana. In our dataset, we analyzed gene expression profiles using organs from worker, drone, and queen bees. This gene-expression profile helped us obtain more detailed information related to organ-specific genes, immune response, detoxification mechanisms, venom-specific genes, and ovary development. From our result, we found 4096 transcripts representing different gene-expression pattern in each organ. Our results suggest that caste-specific transcripts of each organ were expressed differently even under natural conditions. These transcriptome-wide analyses provide new insights into A. cerana and that promote honeybee research and conservation.

Genome-wide identification of miRNAs and target regulatory network in the invasive ectoparasitic mite Varroa destructor

Varroa destructor is an ectoparasite mite that attacks bees leading to colony disorders worldwide. microRNAs (miRNAs) are key molecules used by eukaryotes to post-transcriptional control of gene expression. Nevertheless, still lack information aboutV. destructor miRNAs and its regulatory networks. Here, we used an integrative strategy to characterize the miRNAs in the V. destructor mite. We identified 310 precursors that give rise to 500 mature miRNAs, which 257 are likely mite-specific elements. miRNAs showed canonical length ranging between 18 and 25 nucleotides and 5' uracil preference. Top 10 elements concentrated over 80% of total miRNA expression, with bantam alone representing ~50%. We also detected non-templated bases in precursor-derived small RNAs, indicative of miRNA post-transcriptional regulatory mechanisms. Finally, we note that conserved miRNAs control similar processes in different organisms, suggesting a conservative role. Altogether, our findings contribute to the better understanding of the mite biology that can assist future studies on varroosis control.

10-hydroxy-2E-decenoic acid (10HDA) does not promote caste differentiation in Melipona scutellaris stingless bees

In bees from genus Melipona, differential feeding is not enough to fully explain female polyphenism. In these bees, there is a hypothesis that in addition to the environmental component (food), a genetic component is also involved in caste differentiation. This mechanism has not yet been fully elucidated and may involve epigenetic and metabolic regulation. Here, we verified that the genes encoding histone deacetylases HDAC1 and HDAC4 and histone acetyltransferase KAT2A were expressed at all stages of Melipona scutellaris, with fluctuations between developmental stages and castes. In larvae, the HDAC genes showed the same profile of Juvenile Hormone titers-previous reported-whereas the HAT gene exhibited the opposite profile. We also investigated the larvae and larval food metabolomes, but we did not identify the putative queen-fate inducing compounds, geraniol and 10-hydroxy-2E-decenoic acid (10HDA). Finally, we demonstrated that the histone deacetylase inhibitor 10HDA-the major lipid component of royal jelly and hence a putative regulator of honeybee caste differentiation-was unable to promote differentiation in queens in Melipona scutellaris. Our results suggest that epigenetic and hormonal regulations may act synergistically to drive caste differentiation in Melipona and that 10HDA is not a caste-differentiation factor in Melipona scutellaris.

The microRNA miR-14 Regulates Egg-Laying by Targeting EcR in Honeybees (Apis mellifera)

Honeybees (Apis mellifera) are important pollinators and are commonly used for honey production. The oviposition behavior in honeybees is complex and errors in oviposition could affect the development of the bee colony. Recent studies reported that RNA-RNA cross-talk played a critical role in several biological processes, including reproduction. Ecdysone receptor (EcR) and miR-14 were previously reported to play important roles in egg-laying. Moreover, EcR was predicted to be the target gene of miR-14 and may form miR-14-EcR cross-talk. In this study, knocking down and overexpression of miR-14 and EcR in queen model were implemented. The effect of RNA expression of miR-14 and EcR on the number of eggs laid by honeybee queens were analyzed. Further, luciferase assay was used to confirm the target relation between miR-14 and 3'UTR of EcR. The results showed that the expression of miR-14 and EcR was associated with the number of eggs laid by queens. In specific, inhibition of miR-14 expression enhanced the number of eggs laid, while overexpression of EcR enhanced the number of eggs laid. Lastly, we determined that miR-14 directly targets the mRNA of EcR. These findings suggest that the cross-talk of miR-14-EcR plays an important role in the number of eggs laid by honeybee queens.

Gene Duplication in the Honeybee: Patterns of DNA Methylation, Gene Expression, and Genomic Environment

Gene duplication serves a critical role in evolutionary adaptation by providing genetic raw material to the genome. The evolution of duplicated genes may be influenced by epigenetic processes such as DNA methylation, which affects gene function in some taxa. However, the manner in which DNA methylation affects duplicated genes is not well understood. We studied duplicated genes in the honeybee Apis mellifera, an insect with a highly sophisticated social structure, to investigate whether DNA methylation was associated with gene duplication and genic evolution. We found that levels of gene body methylation were significantly lower in duplicate genes than in single-copy genes, implicating a possible role of DNA methylation in postduplication gene maintenance. Additionally, we discovered associations of gene body methylation with the location, length, and time since divergence of paralogous genes. We also found that divergence in DNA methylation was associated with divergence in gene expression in paralogs, although the relationship was not completely consistent with a direct link between DNA methylation and gene expression. Overall, our results provide further insight into genic methylation and how its association with duplicate genes might facilitate evolutionary processes and adaptation.

Effects of Dinotefuran on Brain miRNA Expression Profiles in Young Adult Honey Bees (Hymenopptera: Apidae)

Honey bees are important pollinators of wild plants and crops. MicroRNAs (miRNAs) are endogenous regulators of gene expression. In this study, we initially determined that the lethal concentration 50 (LC50) of dinotefuran was 0.773 mg/l. Then, the expression profiles and differentially expressed miRNAs (DE miRNAs) in honey bee brains after 1, 5, and 10 d of treatment with the lethal concentration 10 (LC10) of dinotefuran were explored via deep small-RNA sequencing and bioinformatics. In total, 2, 23, and 27 DE miRNAs were identified after persistent exposure to the LC10 of dinotefuran for 1, 5, and 10 d, respectively. Some abundant miRNAs, such as ame-miR-375-3p, ame-miR-281-5p, ame-miR-3786-3p, ame-miR-10-5p, and ame-miR-6037-3p, were extremely significantly differentially expressed. Enrichment analysis suggested that the candidate target genes of the DE miRNAs are involved in the regulation of biological processes, cellular processes, and behaviors. These results expand our understanding of the regulatory roles of miRNAs in honey bee Apis mellifera (Hymenopptera: Apidae) responses to neonicotinoid insecticides and facilitate further studies on the functions of miRNAs in honey bees.

Epigenetic Molecular Mechanisms in Insects

Insects are the largest animal group on Earth both in biomass and diversity. Their outstanding success has inspired genetics and developmental research, allowing the discovery of dynamic process explaining extreme phenotypic plasticity and canalization. Epigenetic molecular mechanisms (EMMs) are vital for several housekeeping functions in multicellular organisms, regulating developmental, ontogenetic trajectories and environmental adaptations. In Insecta, EMMs are involved in the development of extreme phenotypic divergences such as polyphenisms and eusocial castes. Here, we review the history of this research field and how the main EMMs found in insects help to understand their biological processes and diversity. EMMs in insects confer them rapid response capacity allowing insect either to change with plastic divergence or to keep constant when facing different stressors or stimuli. EMMs function both at intra as well as transgenerational scales, playing important roles in insect ecology and evolution. We discuss on how EMMs pervasive influences in Insecta require not only the control of gene expression but also the dynamic interplay of EMMs with further regulatory levels, including genetic, physiological, behavioral, and environmental among others, as was earlier proposed by the Probabilistic Epigenesis model and Developmental System Theory.

Special Significance of Non-Drosophila Insects in Aging

Aging is the leading risk factor of human chronic diseases. Understanding of aging process and mechanisms facilitates drug development and the prevention of aging-related diseases. Although many aging studies focus on fruit fly as a canonical insect system, minimal attention is paid to the potentially significant roles of other insects in aging research. As the most diverse group of animals, insects provide many aging types and important complementary systems for aging studies. Insect polyphenism represents a striking example of the natural variation in longevity and aging rate. The extreme intraspecific variations in the lifespan of social insects offer an opportunity to study how aging is differentially regulated by social factors. Insect flight, as an extremely high-intensity physical activity, is suitable for the investigation of the complex relationship between metabolic rate, oxidative stress, and aging. Moreover, as a "non-aging" state, insect diapause not only slows aging process during diapause phase but also affects adult longevity during/after diapause. In the past two decades, considerable progress has been made in understanding the molecular basis of aging regulation in insects. Herein, the recent research progress in non-Drosophila insect aging was reviewed, and its potential utilization in aging in the future was discussed.

Eusociality influences the strength of negative selection on insect genomes

While much of the focus of sociobiology concerns identifying genomic changes that influence social behaviour, we know little about the consequences of social behaviour on genome evolution. It has been hypothesized that social evolution can influence the strength of negative selection via two mechanisms. First, division of labour can influence the efficiency of negative selection in a caste-specific manner; indirect negative selection on worker traits is theoretically expected to be weaker than direct selection on queen traits. Second, increasing social complexity is expected to lead to relaxed negative selection because of its influence on effective population size. We tested these two hypotheses by estimating the strength of negative selection in honeybees, bumblebees, paper wasps, fire ants and six other insects that span the range of social complexity. We found no consistent evidence that negative selection was significantly stronger on queen-biased genes relative to worker-biased genes. However, we found strong evidence that increased social complexity reduced the efficiency of negative selection. Our study clearly illustrates how changes in behaviour can influence patterns of genome evolution by modulating the strength of natural selection.

Brain microRNAs among social and solitary bees

Evolutionary transitions to a social lifestyle in insects are associated with lineage-specific changes in gene expression, but the key nodes that drive these regulatory changes are unknown. We examined the relationship between social organization and lineage-specific microRNAs (miRNAs). Genome scans across 12 bee species showed that miRNA copy-number is mostly conserved and not associated with sociality. However, deep sequencing of small RNAs in six bee species revealed a substantial proportion (20-35%) of detected miRNAs had lineage-specific expression in the brain, 24-72% of which did not have homologues in other species. Lineage-specific miRNAs disproportionately target lineage-specific genes, and have lower expression levels than shared miRNAs. The predicted targets of lineage-specific miRNAs are not enriched for genes with caste-biased expression or genes under positive selection in social species. Together, these results suggest that novel miRNAs may coevolve with novel genes, and thus contribute to lineage-specific patterns of evolution in bees, but do not appear to have significant influence on social evolution. Our analyses also support the hypothesis that many new miRNAs are purged by selection due to deleterious effects on mRNA targets, and suggest genome structure is not as influential in regulating bee miRNA evolution as has been shown for mammalian miRNAs.

Exploring DNA Methylation Diversity in the Honey Bee Brain by Ultra-Deep Amplicon Sequencing

Understanding methylation dynamics in organs or tissues containing many different cell types is a challenging task that cannot be efficiently addressed by the low-depth bisulphite sequencing of DNA extracted from such sources. Here we explored the feasibility of ultra-deep bisulphite sequencing of long amplicons to reveal the brain methylation patterns in three selected honey bee genes analysed across five distinct conditions on the Illumina MiSeq platform. By combing 15 libraries in one run we achieved a very high sequencing depth of 240,000–340,000 reads per amplicon, suggesting that most of the cell types in the honey bee brain, containing approximately 1 million neurons, are represented in this dataset. We found a small number of gene-specific patterns for each condition in individuals of different ages and performing distinct tasks with 80–90% of those were represented by no more than a dozen patterns. One possibility is that such a small number of frequent patterns is the result of differentially methylated epialleles, whereas the rare and less frequent patterns reflect activity-dependent modifications. The condition-specific methylation differences within each gene appear to be position-dependent with some CpGs showing significant changes and others remaining stable in a methylated or non-methylated state. Interestingly, no significant loss of methylation was detected in very old individuals. Our findings imply that these diverse patterns represent a special challenge in the analyses of DNA methylation in complex tissues and organs that cannot be investigated by low-depth genome-wide bisulphite sequencing. We conclude that ultra-deep sequencing of gene-specific amplicons combined with genotyping of differentially methylated epialleles is an effective way to facilitate more advanced neuro-epigenomic studies in honey bees and other insects.

Identification and characterization of miRNAs in an endoparasitoid wasp, Pteromalus puparum

MicroRNAs (miRNAs) are a form of endogenous small noncoding RNAs that regulate protein-coding gene expression at the posttranscriptional level. So far, knowledge of miRNAs in parasitoids remains rudimentary. We investigated miRNAs in Pteromalus puparum, a pupal endoparasitoid wasp with genome and transcriptome sequences completed. In this study, we constructed eight small RNA libraries from selected developmental stages and genders: male embryos, male larvae, male pupae, male adults, mixed-sex embryos, mixed-sex larvae, mixed-sex pupae, and female adults. We identified 254 mature miRNAs with 5p/3p arm features originated from 75 known and 119 novel miRNA genes in P. puparum, 88 of which reside in 26 clusters. The miRNAs in more than half of the clusters exhibit a consistent expression pattern, indicating they were co-transcribed from a long transcript. Comparing miRNA expression in the eight libraries, we found that 84 mature miRNAs were differentially expressed between embryos and larvae, 20 between larvae and pupae, and 26 between pupae and adults. We found some miRNAs were differentially expressed between sexes in embryos (10), larvae (29), pupae (8), and adults (14). Target predictions resulted in 211,571 miRNA-mRNA interactions for 254 different mature miRNAs. These miRNAs may be involved in sexual and developmental regulation of gene expression.

Small RNAs in eucaryotes: new clues for amplifying microRNA benefits

miRNAs, the smallest nucleotide molecules able to regulate gene expression at post transcriptional level, are found in both animals and plants being involved in fundamental processes for growth and development of living organisms. The number of miRNAs has been hypothesized to increase when some organisms specialized the process of mastication and grinding of food. Further to the vertical transmission, miRNAs can undergo horizontal transmission among different species, in particular between plants and animals. In the last years, an increasing number of studies reported that miRNA passage occurs through feeding, and that in animals, plant miRNAs can survive the gastro intestinal digestion and transferred by blood into host cells, where they can exert their functions modulating gene expression. The present review reports studies on miRNAs during evolution, with particular focus on biogenesis and mechanisms regulating their stability in plants and animals. The different biogenesis and post biogenesis modifications allow to discriminate miRNAs of plant origin from those of animal origin, and make it possible to better clarify the controversial question on whether a possible cross-kingdom miRNA transfer through food does exist. The majority of human medicines and supplements derive from plants and a regular consumption of plant food is suggested for their beneficial effects in the prevention of metabolic diseases, cancers, and dietary related disorders. So far, these beneficial effects have been generally attributed to the content of secondary metabolites, whereas mechanisms regarding other components remain unclear. Therefore, in light of the above reported studies miRNAs could result another component for the medical properties of plants. miRNAs have been mainly studied in mammals characterizing their sequences and molecular targets as available in public databases. The herein presented studies provide evidences that miRNA situation is much more complex than the static situation reported in databases. Indeed, miRNAs may have redundant activities, variable sequences, different methods of biogenesis, and may be differently influenced by external and environmental factors. In-depth knowledge of mechanisms of synthesis, regulation and transfer of plant miRNAs to other species can open new frontiers in the therapy of many human diseases, including cancer.

Plant miRNA Cross-Kingdom Transfer Targeting Parasitic and Mutualistic Organisms as a Tool to Advance Modern Agriculture

MicroRNAs (miRNAs), defined as small non-coding RNA molecules, are fine regulators of gene expression. In plants, miRNAs are well-known for regulating processes spanning from cell development to biotic and abiotic stress responses. Recently, miRNAs have been investigated for their potential transfer to distantly related organisms where they may exert regulatory functions in a cross-kingdom fashion. Cross-kingdom miRNA transfer has been observed in host-pathogen relations as well as symbiotic or mutualistic relations. All these can have important implications as plant miRNAs can be exploited to inhibit pathogen development or aid mutualistic relations. Similarly, miRNAs from eukaryotic organisms can be transferred to plants, thus suppressing host immunity. This two-way lane could have a significant impact on understanding inter-species relations and, more importantly, could leverage miRNA-based technologies for agricultural practices. Additionally, artificial miRNAs (amiRNAs) produced by engineered plants can be transferred to plant-feeding organisms in order to specifically regulate their cross-kingdom target genes. This minireview provides a brief overview of cross-kingdom plant miRNA transfer, focusing on parasitic and mutualistic relations that can have an impact on agricultural practices and discusses some opportunities related to miRNA-based technologies. Although promising, miRNA cross-kingdom transfer remains a debated argument. Several mechanistic aspects, such as the availability, transfer, and uptake of miRNAs, as well as their potential to alter gene expression in a cross-kingdom manner, remain to be addressed.

Analysis of miRNAs in the Heads of Different Castes of the Bumblebee Bombus lantschouensis (Hymenoptera: Apidae)

Bumblebees are important insect pollinators for many wildflowers and crops. MicroRNAs (miRNAs) are endogenous non-coding small RNAs that regulate different biological functions in insects. In this study, the miRNAs in the heads of the three castes of the bumblebee Bombus lantschouensis were identified and characterized by small RNA deep sequencing. The significant differences in the expression of miRNAs and their target genes were analyzed. The results showed that the length of the small RNA reads from males, queens, and workers was distributed between 18 and 30 nt, with a peak at 22 nt. A total of 364 known and 89 novel miRNAs were identified from the heads of the three castes. The eight miRNAs with the highest expressed levels in males, queens, and workers were identical, although the order of these miRNAs based on expression differed. The male vs. queen, male vs. worker, and worker vs. queen comparisons identified nine, fourteen, and four miRNAs with significant differences in expression, respectively. The different castes were clustered based on the differentially expressed miRNAs (DE miRNAs), and the expression levels of the DE miRNAs obtained by RT-qPCR were consistent with the read counts obtained through Solexa sequencing. The putative target genes of these DE miRNAs were enriched in 29 Gene Ontology (GO) terms, and catalytic activity was the most enriched GO term, as demonstrated by its association with 2837 target genes in the male vs. queen comparison, 3535 target genes in the male vs. worker comparison, and 2185 target genes in the worker vs. queen comparison. This study highlights the characteristics of the miRNAs in the three B. lantschouensis castes and will aid further studies on the functions of miRNAs in bumblebees.

Methylation and gene expression differences between reproductive and sterile bumblebee workers

Phenotypic plasticity is the production of multiple phenotypes from a single genome and is notably observed in social insects. Multiple epigenetic mechanisms have been associated with social insect plasticity, with DNA methylation being explored to the greatest extent. DNA methylation is thought to play a role in caste determination in Apis mellifera, and other social insects, but there is limited knowledge on its role in other bee species. In this study, we analyzed whole genome bisulfite sequencing and RNA-seq data sets from head tissue of reproductive and sterile castes of the eusocial bumblebee Bombus terrestris. We found that genome-wide methylation in B. terrestris is similar to other holometabolous insects and does not differ between reproductive castes. We did, however, find differentially methylated genes between castes, which are enriched for multiple biological processes including reproduction. However, we found no relationship between differential methylation and differential gene expression or differential exon usage between castes. Our results also indicate high intercolony variation in methylation. These findings suggest that methylation is associated with caste differences but may serve an alternate function, other than direct caste determination in this species. This study provides the first insights into the nature of a bumblebee caste-specific methylome as well as its interaction with gene expression and caste-specific alternative splicing, providing greater understanding of the role of methylation in phenotypic plasticity within social bee species. Future experimental work is needed to determine the function of methylation and other epigenetic mechanisms in insects.

Unique DNA Methylation Profiles Are Associated with cis-Variation in Honey Bees

DNA methylation is an important epigenetic modification that mediates diverse processes such as cellular differentiation, phenotypic plasticity, and genomic imprinting. Mounting evidence suggests that local DNA sequence variation can be associated with particular DNA methylation states, indicating that the interplay between genetic and epigenetic factors may contribute synergistically to the phenotypic complexity of organisms. Social insects such as ants, bees, and wasps have extensive phenotypic plasticity manifested in their different castes, and this plasticity has been associated with variation in DNA methylation. Yet, the influence of genetic variation on DNA methylation state remains mostly unknown. Here we examine the importance of sequence-specific methylation at the genome-wide level, using whole-genome bisulfite sequencing of the semen of individual honey bee males. We find that individual males harbor unique DNA methylation patterns in their semen, and that genes that are more variable at the epigenetic level are also more likely to be variable at the genetic level. DNA sequence variation can affect DNA methylation by modifying CG sites directly, but can also be associated with local variation in cis that is not CG-site specific. We show that covariation in sequence polymorphism and DNA methylation state contributes to the individual-specificity of epigenetic marks in social insects, which likely promotes their retention across generations, and their capacity to influence evolutionary adaptation.

Comparative Identification of MicroRNAs in Apis cerana cerana Workers' Midguts in Responseto Nosema ceranae Invasion

Here, the expression profiles and differentially expressed miRNAs (DEmiRNAs) in the midguts of Apis cerana cerana workers at 7 d and 10 d post-inoculation (dpi) with N. ceranae were investigated via small RNA sequencing and bioinformatics. Five hundred and twenty nine (529) known miRNAs and 25 novel miRNAs were identified in this study, and the expression of 16 predicted miRNAs was confirmed by Stem-loop RT-PCR. A total of 14 DEmiRNAs were detected in the midgut at 7 dpi, including eight up-regulated and six down-regulated miRNAs, while 12 DEmiRNAs were observed in the midgut at 10 dpi, including nine up-regulated and three down-regulated ones. Additionally, five DEmiRNAs were shared, while nine and seven DEmiRNAs were specifically expressed in midguts at 7 dpi and 10 dpi. Gene ontology analysis suggested some DEmiRNAs and corresponding target mRNAs were involved in various functions including immune system processes and response to stimulus. KEGG pathway analysis shed light on the potential functions of some DEmiRNAs in regulating target mRNAs engaged in material and energy metabolisms, cellular immunity and the humoral immune system. Further investigation demonstrated a complex regulation network between DEmiRNAs and their target mRNAs, with miR-598-y, miR-252-y, miR-92-x and miR-3654-y at the center. Our results can facilitate future exploration of the regulatory roles of miRNAs in host responses to N. ceranae, and provide potential candidates for further investigation of the molecular mechanisms underlying eastern honeybee-microsporidian interactions.

Caste-specific microRNA expression in termites: insights into soldier differentiation

Eusocial insects have polyphenic caste systems in which each caste exhibits characteristic morphology and behaviour. In insects, caste systems arose independently in different lineages, such as Isoptera and Hymenoptera. Although partial molecular mechanisms for the development of eusociality in termites have been clarified by the functional analysis of genes and hormones, the contribution of microRNAs (miRNAs) to caste differentiation is unknown. To understand the role of miRNAs in termite caste polyphenism, we performed small RNA sequencing in a subterranean termite (Reticulitermes speratus) and identified the miRNAs that were specifically expressed in the soldier and worker castes. Of the 550 miRNAs annotated in the R. speratus genome, 74 were conserved in insects and 174 were conserved in other termite species. We found that eight miRNAs (mir-1, mir-125, mir-133, mir-2765, mir-87a and three termite-specific miRNAs) are differentially expressed (DE) in soldiers and workers of R. speratus. This differential expression was experimentally verified for five miRNAs by real-time quantitative PCR. Further, four of the eight DE miRNAs in soldier and worker termite castes were also differentially expressed in hymenopteran castes. The finding that Isoptera and Hymenoptera shared several DE miRNAs amongst castes suggests that these miRNAs evolved independently in these phylogenetically distinct lineages.

You are what you eat: Sequence analysis reveals how plant microRNAs may regulate the human genome

BACKGROUND

Nutrigenomic has revolutionized our understanding of nutrition. As plants make up a noticeable part of our diet, in the present study we chose microRNAs of edible plants and investigated if they can perfectly match human genes, indicating potential regulatory functionalities.

METHODS

miRNAs were obtained using the PNRD database. Edible plants were separated and microRNAs in common in at least four of them entered our analysis. Using vmatchPattern, these 64 miRNAs went through four steps of refinement to improve target prediction: Alignment with the whole genome (2581 results), filtered for those in gene regions (1371 results), filtered for exon regions (66 results) and finally alignment with the human CDS (41 results). The identified genes were further analyzed in-silico to find their functions and relations to human diseases.

RESULTS

Four common plant miRNAs were identified to match perfectly with 22 human transcripts. The identified target genes were involved in a broad range of body functions, from muscle contraction to tumor suppression. We could also indicate some connections between these findings and folk herbology and botanical medicine.

CONCLUSIONS

The food that we regularly eat has a great potential in affecting our genome and altering body functions. Plant miRNAs can provide means of designing drugs for a vast range of health problems including obesity and cancer, since they target genes involved in cell cycle (CCNC), digestion (GIPR) and muscular contractions (MYLK). They can also target regions of CDS for which we still have no sufficient information, to help boost our knowledge of the human genome.

Purification and Functional Analysis of Plant Long Noncoding RNAs (lncRNA)

More than 70% of eukaryotic genomes are transcribed into RNA transcripts, the majority of these transcripts are noncoding protein, and their biological functions are largely unknown. Over the last decade, the application of high-throughput sequencing technologies has led to the description of almost all cellular coding and noncoding RNA transcripts except perhaps for those transcripts that are lowly abundant or those present only in specific cells that are underrepresented in sampled tissue(s). An often underrepresented class of noncoding are long noncoding RNAs (lncRNAs), and these often play key regulatory functions for many biological processes such as cell identity and cell division. However, the purification and functional characterization in vitro are still a challenge in both animal and plant experimental systems. Here, we describe in detail methodology for purification of specific cell types, bioinformatic annotation of lncRNAs, and investigation of biological function using the reference plant Arabidopsis thaliana.