MicroRNA signatures characterizing caste-independent ovarian activity in queen and worker honeybees (Apis mellifera L.)

Deformed wing virus genotypes A and B do not elicit immunologically different responses in naïve honey bee hosts

Iflavirus aladeformis (Picornavirales: Iflaviridae), commonly known as deformed wing virus(DWV), in association with Varroa destructor Anderson and Trueman (Mesostigmata: Varroidae), is a leading factor associated with honey bee (Apis mellifera L. [Hymenoptera: Apidae]) deaths. The virus and mite have a near global distribution, making it difficult to separate the effect of one from the other. The prevalence of two main DWV genotypes (DWV-A and DWV-B) has changed over time, leading to the possibility that the two strains elicit a different immune response by the host. Here, we use a honey bee population naïve to both the mite and the virus to investigate if honey bees show a different immunological response to DWV genotypes. We examined the expression of 19 immune genes by reverse transcription quantitative PCR (RT-qPCR) and analysed small RNA after experimental injection with DWV-A and DWV-B. We found no evidence that DWV-A and DWV-B elicit different immune responses in honey bees. RNA interference genes were up-regulated during DWV infection, and small interfering RNA (siRNA) responses were proportional to viral loads yet did not inhibit DWV accumulation. The siRNA response towards DWV was weaker than the response to another honey bee pathogen, Triatovirus nigereginacellulae (Picornavirales: Dicistroviridae; black queen cell virus), suggesting that DWV is comparatively better at evading host antiviral defences. There was no evidence for the production of virus-derived Piwi-interacting RNAs (piRNAs) in response to DWV. In contrast to previous studies, and in the absence of V. destructor, we found no evidence that DWV has an immunosuppressive effect. Overall, our results advance our understanding of the immunological effect that DWV in isolation elicits in honey bees.

Regulation of insect behavior by non-coding RNAs

The adaptation of insects to environments relies on a sophisticated set of behaviors controlled by molecular and physiological processes. Over the past several decades, accumulating studies have unveiled the roles of non-coding RNAs (ncRNAs) in regulating insect behaviors. ncRNAs assume particularly pivotal roles in the behavioral plasticity of insects by rapidly responding to environmental stimuli. ncRNAs also contribute to the maintenance of homeostasis of insects by fine-tuning the expression of target genes. However, a comprehensive review of ncRNAs' roles in regulating insect behaviors has yet to be conducted. Here, we present the recent progress in our understanding of how ncRNAs regulate various insect behaviors, including flight and movement, social behavior, reproduction, learning and memory, and feeding. We refine the intricate mechanisms by which ncRNAs modulate the function of neural, motor, reproductive, and other physiological systems, as well as gene expression in insects like fruit flies, social insects, locusts, and mosquitos. Furthermore, we discuss potential avenues for future studies in ncRNA-mediated insect behaviors.

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.

Ame-miR-980-3p participates in autophagy-mediated midgut remodelling in Apis mellifera via targeting Atg2B

Autophagy is a process that serves to degrade damaged proteins and organelles, thereby promoting cell homeostasis, differentiation, development and survival. Many miRNAs have been found to have regulatory roles in autophagy. In insects, it has been shown that autophagy is involved in hormone-regulated programmed cell death during metamorphic midgut remodelling. However, whether this is also true during the remodelling of the honey bee midgut is unclear. In the present study, we explored the relationship between autophagy and midgut remodelling and sought to identify miRNAs involved in this physiological process. We found that autophagy occurred during midgut remodelling and that the inhibition of autophagy resulted in midgut dysplasia in prepupae. Differentially expressed miRNAs enriched in the autophagy signalling pathway during midgut remodelling were identified by small RNA-seq. Ame-miR-980-3p, which targets the autophagy-related gene Atg2B, was screened out. Furthermore, abnormal expression of ame-miR-980-3p in the pupal stage led to the thinning of the midgut wall of newly emerged bees (NE). When ame-miR-980-3p expression was inhibited, the intestinal villi of NE bees became significantly shorter and sparse, and the lipid signal in the peritrophic matrix of Pb almost disappeared, indicating that the adult midgut was underdeveloped and the lipid absorption ability was weakened. Taken together, ame-miR-980-3p targeted Atg2B to participate in the regulation of midgut autophagy in the pupae, and the abnormal expression of ame-miR-980-3p would interfere with cell proliferation and death in the process of midgut remodelling, hinder the formation of adult midgut and eventually lead to adult midgut dysplasia and affect the lipid absorption function of the midgut in Apis mellifera.

Age- and caste-independent piRNAs in the germline and miRNA profiles linked to caste and fecundity in the ant Temnothorax rugatulus

Social insects are models for studies of phenotypic plasticity. Ant queens and workers vary in fecundity and lifespan, which are enhanced and extended in queens. Yet, the regulatory mechanisms underlying this variation are not well understood. Ant queens live and reproduce for years, so that they need to protect their germline from transposable element (TE) activity, which may be redundant in short-lived, often sterile workers. We analysed the expression of two protective classes of small RNAs, microRNAs (miRNAs) and Piwi-interacting RNAs (piRNAs), in various tissues, castes and age classes of the ant Temnothorax rugatulus. In queens, piRNAs were highly abundant in ovaries with TEs being their clear targets, with reduced but still detectable piRNA-specific ping-pong signatures in thorax and brains. piRNA pathway activity varied little with age in queens. Moreover, the reduced ovaries of workers also exhibited similar piRNA activity and this not only in young, fertile workers, but also in older foragers with regressed ovaries. Therefore, these ants protect their germline through piRNA activity, regardless of ovarian development, age or caste, even in sterile workers often considered the soma of the superorganism. Our tissue-specific miRNA analysis detected the expression of 304 miRNAs, of which 105 were expressed in all tissues, 10 enriched in the brain, three in the thorax, whereas 83 were ovarian-specific. We identified ovarian miRNAs whose expression was related to caste, fecundity and age, and which likely regulate group-specific gene expression. sRNA shifts in young- to middle-aged queens were minor, suggesting delayed senescence in this reproductive caste.

Exposure to sublethal concentrations of thiacloprid insecticide modulated the expression of microRNAs in honeybees (Apis mellifera L.)

This study aimed to investigate the sublethal effects of thiacloprid on microRNA (miRNA) expression in honeybees (Apis mellifera L.) and the role of a specific miRNA, ame-miR-283-5p, in thiacloprid resistance. The high-throughput small RNA-sequencing was used to analyze global miRNA expression profile changes in honeybees orally exposed to sublethal concentrations of thiacloprid (20 mg/L and 4 mg/L) for 72 h. Thiacloprid at 20 mg/L had no observed adverse effects. In addition, bees were fed with miRNA mimics or inhibitors to increase or decrease ame-miR-283-5p expression, and its effects on P450 gene expression (CYP9Q2 and CYP9Q3) were examined. Thiacloprid susceptibility was also detected. The results showed that treatment with thiacloprid at 20 mg/L and 4 mg/L induced 11 and five differentially expressed miRNAs (DEMs), respectively. Bioinformatic analysis suggested that the DEMs are mainly involved in the regulation of growth and development, metabolism, nerve conduction, and behavior. ame-miR-283-5p was downregulated by both concentrations, which was validated using quantitative real-time reverse transcription PCR analysis. Enhancing ame-miR-283-5p expression significantly inhibited CYP9Q2 mRNA and protein expression, and increased thiacloprid toxicity, while reducing ame-miR-283-5p expression significantly promoted CYP9Q2 expression, and decreased thiacloprid susceptibility. Our study revealed a novel role of miRNAs in insecticide resistance in honeybees.

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.

The Expression of miRNAs Involved in Long-Term Memory Formation in the CNS of the Mollusk Helix lucorum

Mollusks are unique animals with a relatively simple central nervous system (CNS) containing giant neurons with identified functions. With such simple CNS, mollusks yet display sufficiently complex behavior, thus ideal for various studies of behavioral processes, including long-term memory (LTM) formation. For our research, we use the formation of the fear avoidance reflex in the terrestrial mollusk Helix lucorum as a learning model. We have shown previously that LTM formation in Helix requires epigenetic modifications of histones leading to both activation and inactivation of the specific genes. It is known that microRNAs (miRNAs) negatively regulate the expression of genes; however, the role of miRNAs in behavioral regulation has been poorly investigated. Currently, there is no miRNAs sequencing data being published on Helix lucorum, which makes it impossible to investigate the role of miRNAs in the memory formation of this mollusk. In this study, we have performed sequencing and comparative bioinformatics analysis of the miRNAs from the CNS of Helix lucorum. We have identified 95 different microRNAs, including microRNAs belonging to the MIR-9, MIR-10, MIR-22, MIR-124, MIR-137, and MIR-153 families, known to be involved in various CNS processes of vertebrates and other species, particularly, in the fear behavior and LTM. We have shown that in the CNS of Helix lucorum MIR-10 family (26 miRNAs) is the most representative one, including Hlu-Mir-10-S5-5p and Hlu-Mir-10-S9-5p as top hits. Moreover, we have shown the involvement of the MIR-10 family in LTM formation in Helix. The expression of 17 representatives of MIR-10 differentially changes during different periods of LTM consolidation in the CNS of Helix. In addition, using comparative analysis of microRNA expression upon learning in normal snails and snails with deficient learning abilities with dysfunction of the serotonergic system, we identified a number of microRNAs from several families, including MIR-10, which expression changes only in normal animals. The obtained data can be used for further fundamental and applied behavioral research.

miR-12 Derived from Bone Marrow Mesenchymal Stem Cells Accelerates the Development of Human Papillomavirus by Up-Regulating AN1

The miRNA derived from Bone marrow mesenchymal stem cells (BMSCs) have crucial effects on tumors. The tumor could be affected by the abnormal expression of miRNA in human papillomavirus (HPV). Our study aimed to identify the potential brand-new biomarker in order to reveal the pathogenesis of HPV. miRNA derived from BMSCs was detected and identified. The action of miR-12 on biological behavior of HPV was detected. The level of AN1 protein was detected by Western-blot and IHC method. The relationship between miR-12 and AN1 was assessed by bioinformatics analysis and luciferase assay. The tumor cell biological behaviors were evaluated by manipulating miR12 and AN1 level. The tumor volume derived from BMSCs was diminished significantly compared with normal tissues. The tumor volume was bigger after combined injection with Hela cell and miR-12 compared with single injection. The cell proliferative and invasive ability was strengthened after transfection with miR-12mimics. The cell invasive ability was reduced significantly after transfection of si-miR-12. AN1 was a target gene of miR-12 as confirmed by the analysis on bioinformatics and luciferase activity. The phenotype was reversed after the silent presentation of AN1 was disturbed. In conclusion, miR-12 expression is elevated in HPV cells and affects HPV cells through targeting the AN1 signaling pathway.

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.

Temporal analysis of microRNAs associated with wing development in the English grain aphid, Sitobion avenae (F.) (Homoptera: Aphidiae)

Molecular mechanisms underlying wing evolution and development have been a point of scientific inquiry for decades. Phloem-feeding aphids are one of the most devastating global insect pests, where dispersal of winged morphs lead to annual movements, migrations, and range expansions. Aphids show a polyphenic wing dimorphism trait, and offer a model to study the role of environment in determining morphological plasticity of a single genotype. Despite recent progresses in the genetic understanding of wing polyphenism, the influence of environmental cues remains unclear. To investigate the involvement of miRNAs in wing development, we sequenced small RNA libraries of the English grain aphid, Sitobion avenae (F.) across six different developmental stages. As a result, we identified 113 conserved and 193 S. avenae-specific miRNAs. Gene Ontology and KEGG pathway analyses of putative target mRNAs for the six differentially expressed miRNAs are enriched for wing development processes. Dietary uptake of miR-263a, miR-316, and miR-184a agomirs and antagomirs led to significantly higher mortality (>70%) and a lower proportion of winged morphs (<5%). On the other hand, wing malformation was observed in miR-2 and miR-306 agomirs and miR-2 and miR-14 antagomirs, respectively, suggesting their involvement in S. avenae wing morphogenesis. These combined results not only shed light on the regulatory role of miRNAs in wing dimorphism, but also provide potential novel targets for the long-term sustainable management of S. avenae, a devastating global grain pest.

An InR/mir-9a/NlUbx regulatory cascade regulates wing diphenism in brown planthoppers

Wing polymorphism significantly contributes to the ecological success of some insect species. For example, the brown planthopper (BPH) Nilaparvata lugens, which is one of the most destructive rice pests in Asia, can develop into either highly mobile long-winged or highly fecund short-winged adult morphs. A recent study reported a highly provocative result that the Hox gene Ultrabithorax (Ubx) is expressed in BPH forewings and showed that this wing development gene is differentially expressed in nymphs that develop into long-winged versus short-winged morphs. Here, we found that Ubx may be a mir-9a target, and used dual luciferase reporter assays and injected micro RNA (miRNA) mimics and inhibitors to confirm the interactions between mir-9a and NlUbx. We measured the mir-9a and NlUbx expression profiles in nymphs and found that the expression of these two biomolecules was negatively correlated. By rearing BPH nymphs on host rice plants with different nutritional status, we were able to characterize a regulatory cascade between insulin receptor genes, mir-9a, and NlUbx that regulate wing length in BPHs. When host quality was low, NlInR1 expression in the nymph terga increased and NlInR2 expression decreased; this led to a higher mir-9a level, which in turn reduced the NlUbx transcript level and ultimately resulted in longer wing lengths. Beyond extending our understanding of the interplay between host plant status and genetic events that modulate polymorphism, we demonstrated both the upstream signal and miRNA-based regulatory mechanism that control Ubx expression in BPH forewings.

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.

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?'

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.

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.

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.

Single nucleotide polymorphism rs737028527 (G>A) affect miR-1b-3p biogenesis and effects on chicken egg-laying traits

The abundance of miR-1b-3p in the chicken ovary is greater after sexual maturation. In the present study there was assessment of whether a single nucleotine polymorphism (SNP) led to an alteration in expression of reproductive traits. The miR-1b-3p abundance was greatest in ovarian follicles The SNP site of rs737028527 (G > A), located in the 734 bp upstream region of pre-miR-1b-3p, was identified in three different chicken breeds. Results from an association analysis of chicken egg-laying traits indicated the SNP was associated with age at first egg production (AFE) and egg number at 32 and 48 weeks (E32, E48; P < 0.01). Hens with genotype AA had an earlier AFE and greater E32 and E48 than hens with other genotypes. The abundance of mature miR-1b-3p in the hens with the GG genotype was larger than those with the AA genotype (P < 0.01), and the luciferase activity of GG genotype promoter was also greater in birds with the AA than GG genotype (P < 0.05). There was inhibition of the production of the transcription factor bound by the specificity protein 1 (Sp1) as a result of the G-to-A mutation, and the luciferase activity of the GG, but not AA, genotype was markedly increased by Sp1. In conclusion, the SNP, rs737028527 (G> A), affected the abundance of mature miR-1b-3p by Sp1 and was associated with chicken egg-laying traits. Data from the present study allow for an increased understanding of the functions and regulation of miR-1b-3p in ovarian follicle development of hens.

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.

Systematic identification of circular RNAs and corresponding regulatory networks unveil their potential roles in the midguts of eastern honeybee workers

Currently, knowledge of circular RNAs (circRNAs) in insects including honeybee is extremely limited. Here, differential expression profiles and regulatory networks of circRNAs in the midguts of Apis cerana cerana workers were comprehensively investigated using transcriptome sequencing and bioinformatics. In total, 9589 circRNAs (201-800 nt in length) were identified from 8-day-old and 11-day-old workers' midguts (Ac1 and Ac2); among them, 5916 (61.70%) A. cerana cerana circRNAs showed conservation with our previously indentified circRNAs in Apis mellifera ligucstica workers' midguts (Xiong et al., Acta Entomologica Sinica 61:1363-1375, 2018). Five circRNAs were confirmed by RT-PCR and Sanger sequencing. Interestingly, novel_circ_003723, novel_circ_002714, novel_circ_002451, and novel_circ_001980 were highly expressed in both Ac1 and Ac2. In addition, the source genes of circRNAs were involved in 34 GO terms including organelle and cellular process and 141 pathways such as endocytosis and Wnt signaling pathway. Moreover, 55 DEcircRNAs including 34 upregulated and 21 downregulated circRNAs were identified in Ac2 compared with Ac1. circRNA-miRNA regulatory networks indicated that 1060 circRNAs can target 74 miRNAs; additionally, the DEcircRNA-miRNA-mRNA networks suggested that 13 downregulated circRNAs can bind to eight miRNAs and 29 miRNA-targeted mRNAs, while 16 upregulated circRNAs can link to 9 miRNAs and 29 miRNA-targeted mRNAs. These results indicated that DEcircRNAs as ceRNAs may play a comprehensive role in the growth, development, and metabolism of the worker's midgut via regulating source genes and interacting with miRNAs. Notably, eight DEcircRNAs targeting miR-6001-y were likely to be key participants in the midgut development. Our findings not only offer a valuable resource for further studies on A. cerana cerana circRNA and novel insights into understanding the molecular mechanisms underlying the midgut development of eastern honeybee but also provide putative circRNA candidates for functional research in the near future and novel biomarkers for identification of eastern honeybee species including A. cerana cerana and honeybee diseases such as chalkbrood and microsporidiosis.

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.

The microRNA-306/abrupt regulatory axis controls wing and haltere growth in Drosophila

Growth control relies on extrinsic and intrinsic mechanisms that regulate and coordinate the size and pattern of organisms. This control is crucial for a homeostatic development and healthy physiology. The gene networks acting in this process are large and complex: factors involved in growth control are also important in diverse biological processes and these networks include multiple regulators that interact and respond to intra- and extra-cellular inputs that may ultimately converge in the control of the cell cycle. In this work we have studied the function of the Drosophila abrupt gene, coding for a BTB-ZF protein and previously reported to be required for wing vein pattern, in the control of haltere and wing growth. We have found that inactivation of abrupt reduces the size of the wing and haltere. We also found that the microRNA miR-306 controls abrupt expression and that miR-306 and abrupt genetically interact to control wing size. Moreover, the reduced appendage size due to abrupt inactivation is rescued by overexpression of Cyclin-E and by inactivation of dacapo. These findings define a miR-306-abrupt regulatory axis that controls wing and haltere size, whereby miR-306 maintains appropriate levels of abrupt expression which, in turn, regulates the cell cycle. Thus, our results uncover a novel function of abrupt in the regulation of the size of Drosophila appendages during development and contribute to the understanding of the coordination between growth and pattern as well as to the understanding of abrupt oncogenic function in flies.

Reproductive capacity and castes in eusocial stingless bees (Hymenoptera: Apidae)

Eusocial lifestyle is one of the most important transitions in the evolutionary history of some groups of organisms. In bees, there are only two eusocial groups: the honey bees (Apini) and the stingless bees (Meliponini). Despite similarities on the eusocial lifestyles of these taxa, they present profound differences related to caste determination, development, behavior, and reproductive capacity of their members. In most of them the queen has a monopoly on reproduction. However, even though workers are tipically sterile, they can contribute to producing haploid eggs that generate males, or trophic eggs, used as an additional nutrition by the queen.

Integration of lncRNA-miRNA-mRNA reveals novel insights into oviposition regulation in honey bees

BACKGROUND

The honey bee (Apis mellifera) is a highly diverse species commonly used for honey production and pollination services. The oviposition of the honey bee queen affects the development and overall performance of the colony. To investigate the ovary activation and oviposition processes on a molecular level, a genome-wide analysis of lncRNAs, miRNAs and mRNA expression in the ovaries of the queens was performed to screen for differentially expressed coding and noncoding RNAs. Further analysis identified relevant candidate genes or RNAs.

RESULTS

The analysis of the RNA profiles in different oviposition phase of the queens revealed that 740 lncRNAs, 81 miRNAs and 5,481 mRNAs were differently expressed during the ovary activation; 88 lncRNAs, 13 miRNAs and 338 mRNAs were differently expressed during the oviposition inhibition process; and finally, 100 lncRNAs, four miRNAs and 497 mRNAs were differently expressed during the oviposition recovery process. In addition, functional annotation of differentially expressed RNAs revealed several pathways that are closely related to oviposition, including hippo, MAPK, notch, Wnt, mTOR, TGF-beta and FoxO signaling pathways. Furthermore, in the QTL region for ovary size, 73 differentially expressed genes and 14 differentially expressed lncRNAs were located, which are considered as candidate genes affecting ovary size and oviposition. Moreover, a core set of genes served as bridges among different miRNAs were identified through the integrated analysis of lncRNA-miRNA-mRNA network.

CONCLUSION

The observed dramatic expression changes of coding and noncoding RNAs suggest that they may play a critical role in honey bee queens' oviposition. The identified candidate genes for oviposition activation and regulation could serve as a resource for further studies of genetic markers of oviposition in honey bees.

De novo characterization of microRNAs in oriental fruit moth Grapholita molesta and selection of reference genes for normalization of microRNA expression

MicroRNAs (miRNAs) are a group of endogenous non-coding small RNAs that have critical regulatory functions in almost all known biological processes at the post-transcriptional level in a variety of organisms. The oriental fruit moth Grapholita molesta is one of the most serious pests in orchards worldwide and threatens the production of Rosacea fruits. In this study, a de novo small RNA library constructed from mixed stages of G. molesta was sequenced through Illumina sequencing platform and a total of 536 mature miRNAs consisting of 291 conserved and 245 novel miRNAs were identified. Most of the conserved and novel miRNAs were detected with moderate abundance. The miRNAs in the same cluster normally showed correlated expressional profiles. A comparative analysis of the 79 conserved miRNA families within 31 arthropod species indicated that these miRNA families were more conserved among insects and within orders of closer phylogenetic relationships. The KEGG pathway analysis and network prediction of target genes indicated that the complex composed of miRNAs, clock genes and developmental regulation genes may play vital roles to regulate the developmental circadian rhythm of G. molesta. Furthermore, based on the sRNA library of G. molesta, suitable reference genes were selected and validated for study of miRNA transcriptional profile in G. molesta under two biotic and six abiotic experimental conditions. This study systematically documented the miRNA profile in G. molesta, which could lay a foundation for further understanding of the regulatory roles of miRNAs in the development and metabolism in this pest and might also suggest clues to the development of genetic-based techniques for agricultural pest control.

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

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

Evolutionary Developmental Biology (Evo-Devo) Research in Latin America

Famous for its blind cavefish and Darwin's finches, Latin America is home to some of the richest biodiversity hotspots of our planet. The Latin American fauna and flora inspired and captivated naturalists from the nineteenth and twentieth centuries, including such notable pioneers such as Fritz Müller, Florentino Ameghino, and Léon Croizat who made a significant contribution to the study of embryology and evolutionary thinking. But, what are the historical and present contributions of the Latin American scientific community to Evo-Devo? Here, we provide the first comprehensive overview of the Evo-Devo laboratories based in Latin America and describe current lines of research based on endemic species, focusing on body plans and patterning, systematics, physiology, computational modeling approaches, ecology, and domestication. Literature searches reveal that Evo-Devo in Latin America is still in its early days; while showing encouraging indicators of productivity, it has not stabilized yet, because it relies on few and sparsely distributed laboratories. Coping with the rapid changes in national scientific policies and contributing to solve social and health issues specific to each region are among the main challenges faced by Latin American researchers. The 2015 inaugural meeting of the Pan-American Society for Evolutionary Developmental Biology played a pivotal role in bringing together Latin American researchers eager to initiate and consolidate regional and worldwide collaborative networks. Such networks will undoubtedly advance research on the extremely high genetic and phenotypic biodiversity of Latin America, bound to be an almost infinite source of amazement and fascinating findings for the Evo-Devo community.