lncRNA profile of Apis mellifera and its possible role in behavioural transition from nurses to foragers

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.

Regulatory roles of long non-coding RNAs in short-term heat stress in adult worker bees

Long non-coding RNAs (lncRNAs) are crucial modulators of post-transcriptional gene expression regulation, cell fate determination, and disease development. However, lncRNA functions during short-term heat stress in adult worker bees are poorly understood. Here, we performed deep sequencing and bioinformatic analyses of honeybee lncRNAs. RNA interference was performed by using siRNA targeting the most highly expressed lncRNA. The silencing effect on lncRNA and the relative expression levels of seven heat shock protein (HSP) genes, were subsequently examined. Overall, 7,842 lncRNAs and 115 differentially expressed lncRNAs (DELs) were identified in adult worker bees following heat stress exposure. Structural analysis revealed that the overall expression abundance, length of transcripts, exon number, and open reading frames of lncRNAs were lower than those of mRNAs. GO analysis revealed that the target genes were mainly involved in "metabolism," "protein folding," "response to stress," and "signal transduction" pathways. KEGG analysis indicated that the "protein processing in endoplasmic reticulum" and "longevity regulating pathway-multiple species" pathways were most enriched. Quantitative real-time polymerase chain reaction (qRT-PCR) detection of the selected DELs confirmed the reliability of the sequencing data. Moreover, the siRNA experiment indicated that feeding siRNA yielded a silencing efficiency of 77.51% for lncRNA MSTRG.9645.5. Upon silencing this lncRNA, the expression levels of three HSP genes were significantly downregulated (p < 0.05), whereas those of three other HSP genes were significantly upregulated (p < 0.05). Our results provide a new perspective for understanding the regulatory mechanisms of lncRNAs in adult worker bees under short-term heat stress.

Unique spatially and temporary-regulated/sex-specific expression of a long ncRNA, Nb-1, suggesting its pleiotropic functions associated with honey bee lifecycle

Honey bees are social insects, and each colony member has unique morphological and physiological traits associated with their social tasks. Previously, we identified a long non-coding RNA from honey bees, termed Nb-1, whose expression in the brain decreases associated with the age-polyethism of workers and is detected in some neurosecretory cells and octopaminergic neurons, suggesting its role in the regulation of worker labor transition. Herein, we investigated its spatially and temporary-regulated/sex-specific expression. Nb-1 was expressed as an abundant maternal RNA during oogenesis and embryogenesis in both sexes. In addition, Nb-1 was expressed preferentially in the proliferating neuroblasts of the mushroom bodies (a higher-order center of the insect brain) in the pupal brains, suggesting its role in embryogenesis and mushroom body development. On the contrary, Nb-1 was expressed in a drone-specific manner in the pupal and adult retina, suggesting its role in the drone visual development and/or sense. Subcellular localization of Nb-1 in the brain during development differed depending on the cell type. Considering that Nb-1 is conserved only in Apidae, our findings suggest that Nb-1 potentially has pleiotropic functions in the expression of multiple developmental, behavioral, and physiological traits, which are closely associated with the honey bee lifecycle.

MicroRNA Targets PAP1 to Mediate Melanization in Plutella xylostella (Linnaeus) Infected by Metarhizium anisopliae

MicroRNAs (miRNAs) play a pivotal role in important biological processes by regulating post-transcriptional gene expression and exhibit differential expression patterns during development, immune responses, and stress challenges. The diamondback moth causes significant economic damage to crops worldwide. Despite substantial advancements in understanding the molecular biology of this pest, our knowledge regarding the role of miRNAs in regulating key immunity-related genes remains limited. In this study, we leveraged whole transcriptome resequencing data from Plutella xylostella infected with Metarhizium anisopliae to identify specific miRNAs targeting the prophenoloxidase-activating protease1 (PAP1) gene and regulate phenoloxidase (PO) cascade during melanization. Seven miRNAs (pxy-miR-375-5p, pxy-miR-4448-3p, pxy-miR-279a-3p, pxy-miR-3286-3p, pxy-miR-965-5p, pxy-miR-8799-3p, and pxy-miR-14b-5p) were screened. Luciferase reporter assays confirmed that pxy-miR-279a-3p binds to the open reading frame (ORF) and pxy-miR-965-5p to the 3' untranslated region (3' UTR) of PAP1. Our experiments demonstrated that a pxy-miR-965-5p mimic significantly reduced PAP1 expression in P. xylostella larvae, suppressed PO activity, and increased larval mortality rate. Conversely, the injection of pxy-miR-965-5p inhibitor could increase PAP1 expression and PO activity while decreasing larval mortality rate. Furthermore, we identified four LncRNAs (MSTRG.32910.1, MSTRG.7100.1, MSTRG.6802.1, and MSTRG.22113.1) that potentially interact with pxy-miR-965-5p. Interference assays using antisense oligonucleotides (ASOs) revealed that silencing MSTRG.7100.1 and MSTRG.22113.1 increased the expression of pxy-miR-965-5p. These findings shed light on the potential role of pxy-miR-965-5p in the immune response of P. xylostella to M. anisopliae infection and provide a theoretical basis for biological control strategies targeting the immune system of this pest.

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.

Non-coding RNA in the gut of the blood-feeding parasitic worm, Haemonchus contortus

The intestine of Haemonchus contortus is an essential tissue that has been indicated to be a major target for the prevention of haemonchosis caused by this parasitic nematode of small ruminants. Biological peculiarities of the intestine warrant in-depth exploitation, which can be leveraged for future disease control efforts. Here, we determined the intestinal ncRNA (lncRNA, circRNA and miRNA) atlas using whole-transcriptome sequencing and bioinformatics approaches. In total, 4846 novel lncRNA, 982 circRNA, 96 miRNA (65 known and 31 novel) and 8821 mRNA were identified from the H. contortus intestine. The features of lncRNA, circRNA and miRNA were fully characterized. Comparison of miRNA from the intestines and extracellular vesicles supported the speculation that the miRNA from the latter were of intestinal origin in H. contortus. Further function analysis suggests that the cis-lncRNA targeted genes were involved in protein binding, intracellular anatomical structure, organelle and cellular process, whereas the circRNA parental genes were mainly enriched in molecular function categories, such as ribonucleotide binding, nucleotide binding, ATP binding and carbohydrate derivative binding. The miRNA target genes were related to the cellular process, cellular response to stimulus, cellular protein modification process and signal transduction. Moreover, competing endogenous RNA network analysis revealed that the majority of lncRNA, circRNA and mRNA only have one or two binding sites with specific miRNA. Lastly, randomly selected circRNA, lncRNA and miRNA were verified successfully using RT-PCR. Collectively, these data provide the most comprehensive compilation of intestinal transcripts and their functions, and it will be helpful to decipher the biological and molecular complexity of the intestine and lay the foundation for further functional research.

Expression profiles of lncRNAs, miRNAs, and mRNAs and interaction analysis indicate their potential involvement during testicular fusion in Spodoptera litura

Testicular fusion of Spodoptera litura occures during metamorphosis, which benefits sperms development. Previous research identified involvement of ECM-integrin interaction pathways, MMPs in testicular fusion, but the regulatory mechanism remains unclear. RNA-seq was performed to analyze long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) in testes, aiming to uncover potential regulatory mechanisms of testicular fusion. 2150 lncRNAs, 2742 targeted mRNAs, and 347 miRNAs were identified in testes at three different developmental stages. Up-regulated DElncRNAs and DEmRNAs, as well as down-regulated DEmiRNAs, were observed during testicular fusion, while the opposite expression pattern was observed after fusion. Enrichment analysis of DEmRNAs revealed that cAMP signal pathway, ECM remodeling enzymes, ECM-integrin interaction pathways, and cell adhesion molecules were potentially associated with testicular fusion. The identified DElncRNA-DEmiRNA-DEmRNA regulatory network related to cAMP signal pathway, ECM remodeling enzymes suggests their roles during testicular fusion. Our research will provide new targets for studying the mechanism of testicular fusion.

Transcriptomic Analysis of Metarhizium anisopliae-Induced Immune-Related Long Non-Coding RNAs in Polymorphic Worker Castes of Solenopsis invicta

Long non-coding RNAs (lncRNAs) represent a class of RNA molecules that do not encode proteins. Generally studied for their regulatory potential in model insects, relatively little is known about their immunoregulatory functions in different castes of eusocial insects, including Solenopsis invicta, a notoriously invasive insect pest. In the current study, we used Metarhizium anisopliae, an entomopathogenic fungus, to infect the polymorphic worker castes (Major and Minor Workers) and subjected them to RNA sequencing at different intervals (6, 24, and 48 h post-infection (hpi)). Comprehensive bioinformatic analysis identified 5719 (1869 known and 3850 novel) lncRNAs in all libraries. Genomic characteristics analysis showed that S. invicta lncRNAs exhibited structural similarities with lncRNAs from other eusocial insects, including lower exon numbers, shorter intron and exon lengths, and a lower expression profile. A comparison of lncRNAs in major and minor worker ants revealed that several lncRNAs were exclusively expressed in one worker caste and remained absent in the other. LncRNAs such as MSTRG.12029.1, XR_005575440.1 (6 h), MSTRG.16728.1, XR_005575440.1 (24 h), MSTRG.20263.41, and MSTRG.11994.5 (48 h) were only present in major worker ants, while lncRNAs such as MSTRG.8896.1, XR_005574239.1 (6 h), MSTRG.20289.8, XR_005575051.1 (24 h), MSTRG.20289.8, and MSTRG.6682.1 (48 h) were only detected in minor workers. Additionally, we performed real-time quantitative PCR and experimentally validated these findings. Functional annotation of cis-acting lncRNAs in major worker ants showed that lncRNAs targeted genes such as serine protease, trypsin, melanization protease-1, spaetzle-3, etc. In contrast, apoptosis and autophagy-related genes were identified as targets of lncRNAs in minor ants. Lastly, we identified several lncRNAs as precursors of microRNAs (miRNAs), such as miR-8, miR-14, miR-210, miR-6038, etc., indicating a regulatory relationship between lncRNAs, miRNAs, and mRNAs in antifungal immunity. These findings will serve as a genetic resource for lncRNAs in polymorphic eusocial ants and provide a theoretical basis for exploring the function of lncRNAs from a unique and novel perspective.

Phenotypic and genetic characterization of Africanized Apis mellifera colonies with natural tolerance to Varroa destructor and contrasting defensive behavior

Africanized Apis mellifera colonies with promising characteristics for beekeeping have been detected in northern Argentina (subtropical climate) and are considered of interest for breeding programs. Integral evaluation of this feral material revealed high colony strength and resistance/tolerance to brood diseases. However, these Africanized honeybees (AHB) also showed variable negative behavioral traits for beekeeping, such as defensiveness, tendency to swarm and avoidance behavior. We developed a protocol for the selection of AHB stocks based on defensive behavior and characterized contrasting colonies for this trait using NGS technologies. For this purpose, population and behavioral parameters were surveyed throughout a beekeeping season in nine daughter colonies obtained from a mother colony (A1 mitochondrial haplotype) with valuable characteristics (tolerance to the mite Varroa destructor, high colony strength and low defensiveness). A Defensive Behavior Index was developed and tested in the colonies under study. Mother and two daughter colonies displaying contrasting defensive behavior were analyzed by ddRADseq. High-quality DNA samples were obtained from 16 workers of each colony. Six pooled samples, including two replicates of each of the three colonies, were processed. A total of 12,971 SNPs were detected against the reference genome of A. mellifera, 142 of which showed significant differences between colonies. We detected SNPs in coding regions, lncRNA, miRNA, rRNA, tRNA, among others. From the original data set, we also identified 647 SNPs located in protein-coding regions, 128 of which are related to 21 genes previously associated with defensive behavior, such as dop3 and dopR2, CaMKII and ADAR, obp9 and obp10, and members of the 5-HT family. We discuss the obtained results by considering the influence of polyandry and paternal lineages on the defensive behavior in AHB and provide baseline information to use this innovative molecular approach, ddRADseq, to assist in the selection and evaluation of honey bee stocks showing low defensive behavior for commercial uses.

Social insect transcriptomics and the molecular basis of caste diversity

Studies of gene expression provide fundamentally important information on the molecular mechanisms underlying variation in phenotype. Recent technological advances have allowed for the robust study of gene expression through analysis of whole transcriptomes. Here, we review current advances in social insect transcriptomics and discuss their implications in understanding phenotypic diversity. Recent transcriptomic studies provide detailed inventories of the genes involved in producing distinct phenotypes in social species. These investigations have identified key genes and networks involved in producing distinct social insect castes. Nevertheless, questions concerning the evolution of gene expression patterns remain. We suggest a path forward for studying gene expression in future studies of biological systems.

Genome-wide identification of lncRNAs associated with viral infection in Spodoptera frugiperda

The role of lncRNAs in immune defence has been demonstrated in many multicellular and unicellular organisms. However, investigation of the identification and characterization of long non-coding RNAs (lncRNAs) involved in the insect immune response is still limited. In this study, we used RNA sequencing (RNA-seq) to investigate the expression profiles of lncRNAs and mRNAs in the fall armyworm Spodoptera frugiperda in response to virus infection. To assess the tissue- and virus-specificity of lncRNAs, we analysed and compared their expression profiles in haemocytes and fat body of larvae infected with two entomopathogenic viruses with different lifestyles, i.e. the polydnavirus HdIV (Hyposoter didymator IchnoVirus) and the densovirus JcDV (Junonia coenia densovirus). We identified 1883 candidate lncRNAs, of which 529 showed differential expression following viral infection. Expression profiles differed considerably between samples, indicating that many differentially expressed (DE) lncRNAs showed virus- and tissue-specific expression patterns. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment and target prediction analyses indicated that DE-LncRNAs were mainly enriched in metabolic process, DNA replication and repair, immune response, metabolism of insect hormone and cell adhesion. In addition, we identified three DE-lncRNAs potentially acting as microRNA host genes, suggesting that they participate in gene regulation by producing miRNAs in response to virus infection. This study provides a catalogue of lncRNAs expressed in two important immune tissues and potential insight into their roles in the antiviral defence in S. frugiperda. The results may help future in-depth functional studies to better understand the biological function of lncRNAs in interaction between viruses and the fall armyworm.

Recent Advances and Future Potential of Long Non-Coding RNAs in Insects

Over the last decade, long non-coding RNAs (lncRNAs) have witnessed a steep rise in interest amongst the scientific community. Because of their functional significance in several biological processes, i.e., alternative splicing, epigenetics, cell cycle, dosage compensation, and gene expression regulation, lncRNAs have transformed our understanding of RNA's regulatory potential. However, most knowledge concerning lncRNAs comes from mammals, and our understanding of the potential role of lncRNAs amongst insects remains unclear. Technological advances such as RNA-seq have enabled entomologists to profile several hundred lncRNAs in insect species, although few are functionally studied. This article will review experimentally validated lncRNAs from different insects and the lncRNAs identified via bioinformatic tools. Lastly, we will discuss the existing research challenges and the future of lncRNAs in insects.

Feeding Asian honeybee queens with European honeybee royal jelly alters body color and expression of related coding and non-coding RNAs

Background and aims: The Asian honeybee (Apis cerana) and the European honeybee (Apis mellifera) are reproductively isolated. Previous studies reported that exchanging the larval food between the two species, known as nutritional crossbreeding, resulted in obvious changes in morphology, physiology and behavior. This study explored the molecular mechanisms underlying the honeybee nutritional crossbreeding. Methods: This study used full nutritional crossbreeding technology to rear A. cerana queens by feeding them with an A. mellifera royal jelly-based diet in an incubator. The body color and the expression of certain genes, microRNA, lncRNA, and circRNA among nutritional crossbred A. cerana queens (NQ), and control A. cerana queens (CQ) were compared. The biological functions of two target genes, TPH1 and KMO, were verified using RNA interference. Results: Our results showed that the NQ's body color turned yellow compared to the black control queens. Whole transcriptome sequencing results showed that a total of 1484, 311, 92, and 169 DEGs, DElncRNAs, DEmiRNAs, and DEcircRNAs, respectively, were identified in NQ and CQ, in which seven DEGs were enriched for three key pathways (tryptophan, tyrosine, and dopamine) involved in melanin synthesis. Interestingly, eight DElncRNAs and three DEmiRNAs were enriched into the key pathways regulating the above key DEGs. No circRNAs were enriched into these key pathways. Knocking down two key genes (KMO and TPH1) resulted in altered body color, suggesting that feeding NQ's an RNAi-based diet significantly downregulated the expression of TPH1 and KMO in 4-day-old larvae, which confirmed the function of key DEGs in the regulation of honeybee body color. Conclusion: These findings reveal that the larval diets from A. mellifera could change the body color of A. cerana, perhaps by altering the expression of non-coding RNAs and related key genes. This study serves as a model of epigenetic regulation in insect body color induced by environmental factors.

Mblk-1 regulates sugar responsiveness in honey bee (Apis mellifera) foragers

Brain transcriptional regulatory network for behavior demonstrates that brain gene expression in the honey bee can be accurately predicted from the expression transcription factors (TFs), but roles for specific TFs are less understood. Mushroom bodies (MBs) are important for learning, memory and sensory integration in the honey bee brain. A TFs, Mblk-1, expressed preferentially in the large-type Kenyon cells of the honeybee MBs is predicted to be involved in brain function by regulating transcription of its target genes in honey bee. However, its function and the mechanism of regulation in behavior of honey bee is still obscure. Here we show that Mblk-1 had significantly higher expression in the brains of forager bees relative to nurse bees. Mblk-1 was significantly inhibited in bees fed small interfering RNA. In addition, inhibition of Mblk-1 decreased sucrose responsiveness in foragers. Finally, we determined that Mblk-1 regulated the messenger RNA of AmGR1. These findings suggest that Mblk-1 may target AmGR1 to regulate the sucrose responsiveness of foragers.

Emerging Roles and Mechanisms of lncRNA FOXD3-AS1 in Human Diseases

Numerous long noncoding RNAs (lncRNAs) have been identified as powerful regulators of human diseases. The lncRNA FOXD3-AS1 is a novel lncRNA that was recently shown to exert imperative roles in the initialization and progression of several diseases. Emerging studies have shown aberrant expression of FOXD3-AS1 and close correlation with pathophysiological traits of numerous diseases, particularly cancers. More importantly, FOXD3-AS1 was also found to ubiquitously impact a range of biological functions. This study aims to summarize the expression, associated clinicopathological features, major functions and molecular mechanisms of FOXD3-AS1 in human diseases and to explore its possible clinical applications.

Integrative Analysis of lncRNA-mRNA Co-expression Provides Novel Insights Into the Regulation of Developmental Transitions in Female Varroa destructor

Varroa destructor is a major pathogenic driver of the Western honeybee colony losses globally. Understanding the developmental regulation of V. destructor is critical to develop effective control measures. Development is a complex biological process regulated by numerous genes and long non-coding RNAs (lncRNAs); however, the underlying regulation of lncRNAs in the development of V. destructor remains unknown. In this study, we analyzed the RNA sequencing (RNA-Seq) data derived from the four stages of female V. destructor in the reproductive phase (i.e., egg, protonymph, deutonymph, and adult). The identified differentially expressed mRNAs and lncRNAs exhibited a stage-specific pattern during developmental transitions. Further functional enrichment established that fat digestion and absorption, ATP-binding cassette (ABC) transporters, mitogen-activated protein kinase (MAPK) signaling pathway, and ubiquitin-proteasome pathway play key roles in the maturation of female V. destructor. Moreover, the lncRNAs and mRNAs of some pivotal genes were significantly upregulated at the deutonymph stage, such as cuticle protein 65/6.4/63/38 and mucin 5AC, suggesting that deutonymph is the key stage of metamorphosis development and pathogen resistance acquisition for female V. destructor. Our study provides novel insights into a foundational understanding of V. destructor biology.

Understanding of Waggle Dance in the Honey Bee (Apis mellifera) from the Perspective of Long Non-Coding RNA

The ethological study of dance behaviour has yielded some findings since Karl Von Frisch discovered and interpreted the ‘dance language’ in the honey bee. However, the function and role of long non-coding RNAs on dance behaviour are hardly known until now. In this study, the differential expression patterns of lncRNAs in the brains of waggling dancers and non-dancing bees were analysed by RNA sequencing. Furthermore, lncRNA-mRNA association analysis was constructed to decipher the waggle dance. The results of RNA sequencing indicated that a total of 2877 lncRNAs and 9647 mRNAs were detected from honey bee brains. Further comparison analysis displayed that two lncRNAs, MSTRG.6803.3 and XR_003305156.1, may be involved in the waggle dance. The lncRNA-mRNA association analysis showed that target genes of differentially expressed lncRNAs in the brains between waggling dancers and non-dancing bees were mainly annotated in biological processes related to metabolic process, signalling and response to stimulus and in molecular function associated with signal transducer activity, molecular transducer activity and binding. Nitrogen metabolism was likely implicated in the modulation of the waggle dance. Our findings contribute to further understanding the occurrence and development of waggle dance.

Genome-wide analysis of long non-coding RNAs and their association with wing development in Aphis citricidus (Hemiptera: Aphididae)

Long non-coding RNAs (lncRNAs) play critical roles in the various physiological processes of insects. The wing is a successful adaptation allowing insects to escape from unfavorable environments, while information on lncRNAs related to wing development is limited. In this study, we constructed 12 libraries from two RNA-seq comparisons: 4th instar winged nymphs versus winged adults and 4th instar wingless nymphs versus wingless adults in the brown citrus aphid Aphis citricidus, to identify the wing development-associated lncRNAs. A total of 2914 lncRNAs were identified and 50 lncRNAs were differentially expressed during the 4th instar winged nymphs to winged adults transition, and 28 lncRNAs changed during the 4th instar wingless nymphs to wingless adults transition. The differentially expressed lncRNAs were grouped into six clusters according to the expression patterns in the combined two-winged morphs. lncRNA Ac_lnc54106.1 was up-regulated during 4th instar winged nymphs to winged adults transition, but a lack of change during the 4th instar wingless nymphs to wingless adults transition implied a critical role in the specific regulation of wing development. RNA interference of Ac_lnc54106.1 resulted in malformed wings. Targets prediction, expression patterns, and RNAi assay results showed that Ac_lnc54106.1 may target the PiggyBac transposable element-derived protein 4 (PGBD4) gene, decrease expression of the canonical wing development-related genes, and finally regulate wing development. The systematic identification of lncRNAs in an aphid increases our understanding of how non-coding RNA mediates the wing plasticity of insects.

Identification and comparative analysis of long non-coding RNAs in the brain of fire ant queens in two different reproductive states

Background

Many long non-coding RNAs (lncRNAs) have been extensively identified in higher eukaryotic species. The function of lncRNAs has been reported to play important roles in diverse biological processes, including developmental regulation and behavioral plasticity. However, there are no reports of systematic characterization of long non-coding RNAs in the fire ant Solenopsis invicta.

Results

In this study, we performed a genome-wide analysis of lncRNAs in the brains of S. invicta from RNA-seq. In total, 1,393 novel lncRNA transcripts were identified in the fire ant. In contrast to the annotated lncRNA transcripts having at least two exons, novel lncRNAs are monoexonic transcripts with a shorter length. Besides, the transcriptome from virgin alate and dealate mated queens were analyzed and compared. The results showed 295 differentially expressed mRNA genes (DEGs) and 65 differentially expressed lncRNA genes (DELs) between virgin and mated queens, of which 17 lncRNAs were highly expressed in the virgin alates and 47 lncRNAs were highly expressed in the mated dealates. By identifying the DEL:DEG pairs with a high association in their expression (Spearman’s |rho|> 0.8 and p-value < 0.01), many DELs were co-regulated with DEGs after mating. Furthermore, several remarkable lncRNAs (MSTRG.6523, MSTRG.588, and nc909) that were found to associate with particular coding genes may play important roles in the regulation of brain gene expression in reproductive transition in fire ants.

Conclusion

This study provides the first genome-wide identification of S. invicta lncRNAs in the brains in different reproductive states. It will contribute to a fuller understanding of the transcriptional regulation underpinning reproductive changes.

Genome-wide identification and characterization of long non-coding RNAs in Tribolium castaneum

Long non-coding RNAs (lncRNAs) are poorly understood in insects. In this study, we performed genome-wide analysis of lncRNAs in Tribolium castaneum by RNA-seq. In total, 4516 lncRNA transcripts corresponding to 3917 genes were identified from late embryos, early larvae, late larvae, early pupae, late pupae and early adults of T. castaneum, including 3152 novel lncRNAs and 1364 known lncRNAs. These lncRNAs have few exons and transcripts, and are short in length. During development, they exhibited nine different expression patterns. Functionally, they can act either by targeting messenger RNAs (1813 lncRNAs) and lncRNAs (45 lncRNAs) or as micro RNA (miRNA) precursors (46 lncRNAs). LncRNAs were observed to target the metabolic enzymes of glycolysis, TCA cycle and amino acids, demonstrating that lncRNAs control metabolism by regulating metabolic enzymes. Moreover, lncRNAs were shown to participate in cell differentiation and development via their targets. As miRNA precursors, lncRNAs could participate in the ecdysone signaling pathway. This study provides comprehensive information for lncRNAs of T. castaneum, and will promote functional analysis and target identification of lncRNAs in the insect.

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.

Identification of long noncoding RNAs reveals the effects of dinotefuran on the brain in Apis mellifera (Hymenopptera: Apidae)

Background

Dinotefuran (CAS No. 165252–70-0), a neonicotinoid insecticide, has been used to protect various crops against invertebrate pests and has been associated with numerous negative sublethal effects on honey bees. Long noncoding RNAs (lncRNAs) play important roles in mediating various biological and pathological processes, involving transcriptional and gene regulation. The effects of dinotefuran on lncRNA expression and lncRNA function in the honey bee brain are still obscure.

Results

Through RNA sequencing, a comprehensive analysis of lncRNAs and mRNAs was performed following exposure to 0.01 mg/L dinotefuran for 1, 5, and 10 d. In total, 312 lncRNAs and 1341 mRNAs, 347 lncRNAs and 1458 mRNAs, and 345 lncRNAs and 1155 mRNAs were found to be differentially expressed (DE) on days 1, 5 and 10, respectively. Gene set enrichment analysis (GSEA) indicated that the dinotefuran-treated group showed enrichment in carbohydrate and protein metabolism and immune-inflammatory responses such as glycine, serine and threonine metabolism, pentose and glucuronate interconversion, and Hippo and transforming growth factor-β (TGF-β) signaling pathways. Moreover, the DE lncRNA TCONS_00086519 was shown by fluorescence in situ hybridization (FISH) to be distributed mainly in the cytoplasm, suggesting that it may serve as a competing endogenous RNA and a regulatory factor in the immune response to dinotefuran.

Conclusion

This study characterized the expression profile of lncRNAs upon exposure to neonicotinoid insecticides in young adult honey bees and provided a framework for further study of the role of lncRNAs in honey bee growth and the immune response.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07811-y.

(Epi)Genetic Mechanisms Underlying the Evolutionary Success of Eusocial Insects

Eusocial insects, such as bees, ants, and wasps of the Hymenoptera and termites of the Blattodea, are able to generate remarkable diversity in morphology and behavior despite being genetically uniform within a colony. Most eusocial insect species display caste structures in which reproductive ability is possessed by a single or a few queens while all other colony members act as workers. However, in some species, caste structure is somewhat plastic, and individuals may switch from one caste or behavioral phenotype to another in response to certain environmental cues. As different castes normally share a common genetic background, it is believed that much of this observed within-colony diversity results from transcriptional differences between individuals. This suggests that epigenetic mechanisms, featured by modified gene expression without changing genes themselves, may play an important role in eusocial insects. Indeed, epigenetic mechanisms such as DNA methylation, histone modifications and non-coding RNAs, have been shown to influence eusocial insects in multiple aspects, along with typical genetic regulation. This review summarizes the most recent findings regarding such mechanisms and their diverse roles in eusocial insects.

Transcriptome analysis reveals potential function of long non-coding RNAs in 20-hydroxyecdysone regulated autophagy in Bombyx mori

Background

20-hydroxyecdysone (20E) plays important roles in insect molting and metamorphosis. 20E-induced autophagy has been detected during the larval–pupal transition in different insects. In Bombyx mori, autophagy is induced by 20E in the larval fat body. Long non-coding RNAs (lncRNAs) function in various biological processes in many organisms, including insects. Many lncRNAs have been reported to be potential for autophagy occurrence in mammals, but it has not been investigated in insects.

Results

RNA libraries from the fat body of B. mori dissected at 2 and 6 h post-injection with 20E were constructed and sequenced, and comprehensive analysis of lncRNAs and mRNAs was performed. A total of 1035 lncRNAs were identified, including 905 lincRNAs and 130 antisense lncRNAs. Compared with mRNAs, lncRNAs had longer transcript length and fewer exons. 132 lncRNAs were found differentially expressed at 2 h post injection, compared with 64 lncRNAs at 6 h post injection. Thirty differentially expressed lncRNAs were common at 2 and 6 h post-injection, and were hypothesized to be associated with the 20E response. Target gene analysis predicted 6493 lncRNA-mRNA cis pairs and 42,797 lncRNA-mRNA trans pairs. The expression profiles of LNC_000560 were highly consistent with its potential target genes, Atg4B, and RNAi of LNC_000560 significantly decreased the expression of LNC_000560 and Atg4B. These results indicated that LNC_000560 was potentially involved in the 20E-induced autophagy of the fat body by regulating Atg4B.

Conclusions

This study provides the genome-wide identification and functional characterization of lncRNAs associated with 20E-induced autophagy in the fat body of B. mori. LNC_000560 and its potential target gene were identified to be related to 20-regulated autophagy in B. mori. These results will be helpful for further studying the regulatory mechanisms of lncRNAs in autophagy and other biological processes in this insect model.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07692-1.

Long Non-Coding RNAs in Insects

Only a small subset of all the transcribed RNAs are used as a template for protein translation, whereas RNA molecules that are not translated play a very important role as regulatory non-coding RNAs (ncRNAs). Besides traditionally known RNAs (ribosomal and transfer RNAs), ncRNAs also include small non-coding RNAs (sncRNAs) and long non-coding RNAs (lncRNAs). The lncRNAs, which were initially thought to be junk, have gained a great deal attention because of their regulatory roles in diverse biological processes in animals and plants. Insects are the most abundant and diverse group of animals on this planet. Recent studies have demonstrated the role of lncRNAs in almost all aspects of insect development, reproduction, and genetic plasticity. In this review, we describe the function and molecular mechanisms of the mode of action of different insect lncRNAs discovered up to date.

Systematic and computational identification of Androctonus crassicauda long non-coding RNAs

The potential function of long non-coding RNAs in regulating neighbor protein-coding genes has attracted scientists' attention. Despite the important role of lncRNAs in biological processes, a limited number of studies focus on non-model animal lncRNAs. In this study, we used a stringent step-by-step filtering pipeline and machine learning-based tools to identify the specific Androctonus crassicauda lncRNAs and analyze the features of predicted scorpion lncRNAs. 13,401 lncRNAs were detected using pipeline in A. crassicauda transcriptome. The blast results indicated that the majority of these lncRNAs sequences (12,642) have no identifiable orthologs even in closely related species and those considered as novel lncRNAs. Compared to lncRNA prediction tools indicated that our pipeline is a helpful approach to distinguish protein-coding and non-coding transcripts from RNA sequencing data of species without reference genomes. Moreover, analyzing lncRNA characteristics in A. crassicauda uncovered that lower protein-coding potential, lower GC content, shorter transcript length, and less number of isoform per gene are outstanding features of A. crassicauda lncRNAs transcripts.

Forager age and foraging state, but not cumulative foraging activity, affect biogenic amine receptor gene expression in the honeybee mushroom bodies

Foraging behavior is crucial for the development of a honeybee colony. Biogenic amines are key mediators of learning and the transition from in-hive tasks to foraging. Foragers vary considerably in their behavior, but whether and how this behavioral diversity depends on biogenic amines is not yet well understood. For example, forager age, cumulative foraging activity or foraging state may all be linked to biogenic amine signaling. Furthermore, expression levels may fluctuate depending on daytime. We tested if these intrinsic and extrinsic factors are linked to biogenic amine signaling by quantifying the expression of octopamine, dopamine and tyramine receptor genes in the mushroom bodies, important tissues for learning and memory. We found that older foragers had a significantly higher expression of Amdop1, Amdop2, AmoctαR1, and AmoctβR1 compared to younger foragers, whereas Amtar1 showed the opposite pattern. Surprisingly, our measures of cumulative foraging activity were not related to the expression of the same receptor genes in the mushroom bodies. Furthermore, we trained foragers to collect sucrose solution at a specific time of day and tested if the foraging state of time-trained foragers affected receptor gene expression. Bees engaged in foraging had a higher expression of Amdop1 and AmoctβR3/4 than inactive foragers. Finally, the expression of Amdop1, Amdop3, AmoctαR1, and Amtar1 also varied with daytime. Our results show that receptor gene expression in forager mushroom bodies is complex and depends on both intrinsic and extrinsic factors.

Insecticides cause transcriptional alterations of endocrine related genes in the brain of honey bee foragers

Bees are exposed to endocrine active insecticides. Here we assessed expressional alteration of marker genes indicative of endocrine effects in the brain of honey bees. We exposed foragers to chlorpyrifos, cypermethrin and thiacloprid and assessed the expression of genes after exposure for 24 h, 48 h and 72 h. Chlorpyrifos caused the strongest expressional changes at 24 h characterized by induction of vitellogenin, major royal jelly protein (mrjp) 2 and 3, insulin-like peptide (ilp1), alpha-glucosidase (hbg3) and sima, and down-regulation of buffy. Cypermethrin caused minor induction of mrjp1, mrjp2, mmp1 and ilp1. The sima transcript showed down-regulation at 48 h and up-regulation at 72 h. Exposure to thiacloprid caused down-regulation of vitellogenin, mrjp1 and sima at 24 h, and hbg3 at 72 h, as well as induction of ilp1 at 48 h. The buffy transcript was down-regulated at 24 h and up-regulated at 48 h. Despite compound-specific expression patterns, each insecticide altered the expression of some of the suggested endocrine system related genes. Our study suggests that expressional changes of genes prominently expressed in nurse or forager bees, including down-regulation of buffy and mrjps and up-regulation of hbg3 and ilp1 may serve as indicators for endocrine activity of insecticides in foragers.

The neonicotinoid thiacloprid causes transcriptional alteration of genes associated with mitochondria at environmental concentrations in honey bees

Thiacloprid is widely used in agriculture and may affect pollinators. However, its molecular effects are poorly known. Here, we report the global gene expression profile in the brain of honey bee foragers assessed by RNA-sequencing. Bees were exposed for 72 h to nominal concentrations of 25 and 250 ng/bee via sucrose solution. Determined residue concentrations by LC-MS/MS were 0.59 and 5.49 ng/bee, respectively. Thiacloprid exposure led to 5 and 71 differentially expressed genes (DEGs), respectively. Nuclear genes encoding mitochondrial ribosomal proteins and enzymes involved in oxidative phosphorylation, as well as metabolism enzymes and transporters were altered at 5.49 ng/bee. Kyoto Encylopedia of Genes and Genomes (KEGG) analysis revealed that mitochondrial ribosome proteins, mitochondrial oxidative phosphorylation, pyrimidine, nicotinate and nicotinamide metabolism and additional metabolic pathways were altered. Among 21 genes assessed by RT-qPCR, the transcript of farnesol dehydrogenase involved in juvenile hormone III synthesis was significantly down-regulated. Transcripts of cyp6a14-like and apolipophorin-II like protein, cytochrome oxidase (cox17) and the non-coding RNA (LOC102654625) were significantly up-regulated at 5.49 ng/bee. Our findings indicate that thiacloprid causes transcriptional changes of genes prominently associated with mitochondria, particularly oxidative phosphorylation. This highlight potential effects of this neonicotinoid on energy metabolism, which may compromise bee foraging and thriving populations at environmentally relevant concentrations.

Genome-Wide Identification of Long Non-coding RNAs in the Gravid Ectoparasite Varroa destructor

Long non-coding RNAs (lncRNAs) emerge as critical regulators across a wide variety of biological functions in living organisms. However, to date, no systematic characterization of lncRNAs has been investigated in the ectoparasitic mite Varroa destructor, the most severe biotic threat to honey bees worldwide. Here, we performed an initial genome-wide identification of lncRNAs in V. destructor via high-throughput sequencing technology and reported, for the first time, the transcriptomic landscape of lncRNAs in the devastating parasite. By means of a lncRNA identification pipeline, 6,645 novel lncRNA transcripts, encoded by 3,897 gene loci, were identified, including 2,066 sense lncRNAs, 2,772 lincRNAs, and 1,807 lncNATs. Compared with protein-coding mRNAs, V. destructor lncRNAs are shorter in terms of full length, as well as of the ORF length, contain less exons, and express at lower level. GO term and KEGG pathway enrichment analyses of the lncRNA target genes demonstrated that these predicted lncRNAs may be potentially responsible for the regulatory functions of cellular and biological progresses in the reproductive phase of V. destructor. To our knowledge, this is the first catalog of lncRNA profile in the parasitiformes species, providing a valuable resource for genetic and genomic studies. Understanding the characteristics and features of lncRNAs in V. destructor would promote sustainable parasite control.

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.

Functional genomics of parental care of insects

Parental care was likely the first step most lineages made towards sociality. However, the molecular mechanisms that generate parental care are not broadly characterized. Insects are important as an evolutionary independent group from classic models of parental care, such as, house mice. They provide an opportunity to test the generality of our understanding. With this review, I survey the functional genomics of parental care of insects, summarize several recent advances in the broader framework for studying and understanding parental care, and finish with suggested priorities for further research. Although there are too few studies to draw definitive conclusions, I argue that natural selection appears to be rewiring existing gene networks to produce parental care, that the epigenetic mechanisms influencing parental care are not well understood, and, as an interesting early consensus, that genes strongly associated with carer/offspring interactions appear biased towards proteins that are secreted. I summarize the studies that have functionally validate candidate genes and highlight the increasing need to perform this work. I finish with arguments for both conceptual and practical changes moving forward. I argue that future work can increase the use of predictive frameworks, broaden its definition of conservation of mechanism to gene networks rather than single genes, and increase the use of more established comparative methods. I further highlight the practical considerations of standardizing analyses and reporting, increasing the sampling of both carers and offspring, better characterizing gene regulatory networks, better characterizing taxonomically restricted genes and any consistent role they have underpinning parental care, and using factorial designs to disentangle the influence of multiple variables on the expression of parental care.