Next Article in Journal
Identification of Terpenoid Compounds and Toxicity Assays of Essential Oil Microcapsules from Artemisia stechmanniana
Previous Article in Journal
Secondary Rhinaria Contribute to Major Sexual Dimorphism of Antennae in the Aphid Semiaphis heraclei (Takahashi)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Insects
Volume 14
Issue 5
10.3390/insects14050469
insects-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Xiaoxue Fan
1,
Wende Zhang
1,
Sijia Guo
1,
Leran Zhu
1,
Yiqiong Zhang
1,
Haodong Zhao
1,
Xuze Gao
1,
Haibin Jiang
2,
Tianze Zhang
1,
Dafu Chen
1,3,
Rui Guo
1,3,* and
Qingsheng Niu
2,*
1
College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, Fuzhou 350002, China
2
Apiculture Science Institute of Jilin Province, Jilin 132000, China
3
Apitherapy Research Institute of Fujian Province, Fuzhou 350002, China
*
Authors to whom correspondence should be addressed.
Insects 2023, 14(5), 469; https://doi.org/10.3390/insects14050469
Submission received: 17 April 2023 / Revised: 9 May 2023 / Accepted: 12 May 2023 / Published: 16 May 2023
(This article belongs to the Section Insect Molecular Biology and Genomics)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

Currently, little is known about the regulatory role of miRNAs in the development of the Asian honey bee (Apis cerana) gut. On the basis of our previously gained small RNA-seq data, the expression pattern, regulatory network, and putative role of miRNAs in the developmental process of A. cerana worker larval guts were for the first time analyzed. DEmiRNAs were further analyzed. Here, 330, 351, and 321 miRNAs were identified in the 4-, 5-, and 6-day-old larval guts, respectively. Sequences of six miRNAs were verified using Sanger sequencing. Additionally, 15 and 10 differentially expressed miRNAs (DEmiRNAs) were, respectively, detected in the “Ac4 vs. Ac5” and“Ac5 vs. Ac6” comparison groups, and the target genes of these DEmiRNAs were involved in a series of GO terms and KEGG pathways relevant to growth and development, such as Notch and Wnt. miR-6001-y was induced to be activated in the development of larval guts. Based on RT-qPCR, we validated the expression levels of five randomly chosen DEmiRNAs. These findings lay a foundation for clarifying the developmental mechanism of the A. cerana worker larval guts.

Abstract

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

Graphical Abstract

1. Introduction

The honey bee is a kind of representative eusocial insect that plays an essential role in global ecological environment through pollination for a substantial quantity of crops and flowers. Additionally, honey bees can be domesticated for production of abundant and excellent api-products such as royal jelly, honey, propolis, beeswax, and bee venom [1]. Moreover, the honey bee is used as a research model to investigate growth and development [2], neurobiology [3], social behavior [4], host–pathogen interaction [5,6], and epigenetic regulation [7]. Apis cerana is a predominant bee species distributed across eastern, southern, and southeastern Asia and widely used in the beekeeping industry, with tremendous ecological and economic value [8]. In comparison with Apis mellifera, A. cerana has a subseries of advantages such as higher resistance to low temperature, greater ability to collect sporadic nectar sources, longer honey collection period, and stronger resilience to the ectoparasitic mite [9,10]. Apis cerana cerana is an endemic pollinator and commercial farming bee species in China and some other countries in Asia [11].
MicroRNAs (miRNAs) are a kind of small, non-coding RNA (ncRNA) molecule with a distribution range of 18~25 nucleotides (nt) in length, which is evolutionarily conserved and exerts pivotal functions in regulating considerable biological processes at the post-transcriptional level, such as development, metamorphosis, metabolism, immunity, and reproduction [12,13,14]. Increasing studies are suggestive of the participation of miRNAs in modulating the growth and development of various insects [15]. For instance, Shen et al. [15] found that injection of miR-2766 and miR-14 agomirs suppressed the levels of HaTH mRNA and protein and further resulted in defective pupation of Helicoverpa armigera. Wang et al. [16] verified that miR-282-5p in silkworm (Bombyx mori) regulated the expression of chitinase 5 (chi5), which was essential for the normal molting process, and observed that injection of miR-282-5p agomir led to the death of 40% of silkworms due to molting failure. By performing dual-luciferase assay and functional investigation, Zhang et al. [17] detected that miR-309a interacted with the transcription factor pannier (GATA-binding factor A/pnr) and then activated vitellogenin 2 (Vg2) and vitellogenin receptor (VgR) genes, which further facilitated the ovarian development. In honey bees, miRNAs have been shown to be regulators in developmental process [18], oviposition regulation [19], caste determination [20], immune defense [21], as well as organ development [22]. Previously, by performing small RNA-seq (sRNA-seq) and bioinformatics, we discovered 560 miRNAs in the Apis mellifera ligustica larval guts [23] and deciphered the expression profile as well as regulatory network of miRNAs across the developmental process of larval guts [24]. Currently, limited studies have shown that A. cerana miRNAs are potentially involved in response to fungi or viruses [25,26]. For example, Deng et al. [26] identified a total of 23 differently expressed miRNAs (DEmiRNAs) in A. cerana larvae upon CSBV infection and inferred that ame-miR-3759 may play a key part in melanin synthesis during larval response. Recently, we identified 435 miRNAs in the larval guts of A. c. cerana, the nominate subspecies of A. cerana, based on sRNA-seq and bioinformatics and further confirmed the sequences of six miRNAs and the expression of five miRNAs using molecular approaches [27].
The insect gut is a vital tissue responsible for digestion of food, absorption of nutrient, and immune defense [28]. Earlier studies regarding bee guts were mainly concentrated on the intestinal microorganisms [29,30,31]. However, few advancements in the developmental mechanism of the bee guts have been gained until now. Recently, we resolved the expression patterns, regulatory networks, and potential functions of piRNAs, another kind of small ncRNAs, in the larval guts of both A. mellifera and A. cerana [32,33]. However, little is known about the miRNA-regulated developmental mechanism underlying the bee gut at present.
Here, on the basis of our formerly obtained high-quality sRNA-seq datasets, miRNAs in A. c. cerana 4-, 5-, and 6-day-old larval guts were characterized using bioinformatics followed by the molecular validation; additionally, the differential expression pattern of miRNAs and corresponding regulatory networks during the developmental process of larval guts were investigated, with a focus on DEmiRNA-involved networks associated with development-associated signaling pathways as well as cellular and humoral immune. Our data can not only enrich our understanding of the development of the A. cerana larval gut but also lay a foundation for illustration of the mechanism underlying the DEmiRNA-modulated larval gut development.

2. Materials and Methods

2.1. Source of sRNA-seq Datasets

In our previous work, A. c. cerana 4-, 5-, and 6-day-old larvae were reared under laboratory conditions, and the larval gut tissues (Ac4 group: Ac4-1, Ac4-2, Ac4-3; Ac5 group: Ac5-1, Ac5-2, Ac5-3; Ac6 group: Ac6-1, Ac6-2, Ac6-3) were prepared following the established method [34], followed by cDNA library construction, deep sequencing utilizing sRNA-seq technology, and data quality control [27]. The high-quality raw data were deposited in the NCBI SRA database and linked to the BioProject number PRJNA395108.

2.2. Identification and Analysis of miRNAs

A. c. cerana miRNAs were identified according to our previously described protocol [27]: The clean reads were mapped to the reference genome of A. cerana (Assembly ACSNU-2.0), and the mapped clean reads were then aligned to GenBank (https://www.ncbi.nlm.nih.gov/genbank/, accessed on 2 March 2021) and Rfam (11.0) (http://rfam.xfam.org/, accessed on 2 March 2021) databases to remove small ncRNA including rRNA, scRNA, snoRNA, snRNA, and tRNA; further, clean tags of the reference genome were compared with known miRNA precursor sequences in the miRbase database (http://www.miRBase.org, accessed on 2 March 2021) to identify known miRNAs, while novel miRNAs were identified according to whether miRNAs had classical stem-loop structures. For each miRNA, the expression level was normalized to tags per million (TPM) based on the following formula: TPM = T × 106/N (T refers to miRNA clean reads; N refers to total sRNA clean reads).

2.3. Validation of miRNAs by Stem-Loop RT-PCR and Sanger Sequencing

Six miRNAs (miR-14-y, bantam-y, miR-184-y, miR-9-z, miR-8-y, and miR-275-y) were randomly selected for RT-PCR. Specific stem-loop primers and forward primers (F) as well as universal reverse primers (R) for each of these DEmiRNAs were designed using Primer Premier 6 [35] (Table S1) and then synthesized by Sangon Biotech Co., Ltd. (Shanghai, China). Using an RNA extraction kit (Promega, China), total RNA was extracted from 4-, 5-, and 6-day-old A. c. ceranae larvae guts (N = 3) and then divided into two portions. Next, reverse transcription was performed with stem-loop primers, and the resulting cDNA was used as template for PCR. On basis of the method described by Wu et al. [36], PCR amplification was conducted and subjected to 1.8% agarose gel electrophoresis, followed by extraction of the expected fragment, ligation to the pESI-T vector (Yeasen, China), and transformation into Escherichia coli DH5α-competent cells. After microbial PCR, 1 mL of bacterial solution with a positive signal was sent to Sangon Biotech Co., Ltd. (Shanghai, China) for single-end Sanger sequencing.

2.4. Investigation of DEmiRNAs

For each miRNA, the fold change in expression level between “Ac4 vs. Ac5” and “Ac5 vs. Ac6” comparison groups was determined according to the following formula: (TPM in Ac5)/(TPM in Ac4) or (TPM in Ac6)/(TPM in Ac5). Based on the standard of p-value ≤ 0.05 and |log2(Fold change)| ≥ 1, DEmiRNAs in “Ac4 vs. Ac5” and “Ac5 vs. Ac6” comparison groups were filtered out. Upregulated miRNAs refers to expression level of miRNAs between two groups with log2(Fold change) > 1 and downregulated miRNAs with log2(Fold change) < −1. By using the OmicShare platform (https://www.omicshare.com/, accessed on 5 March 2021), Venn analysis as well as expression clustering of the DEmiRNAs in these comparison groups were performed.

2.5. RNA-seq Data Source

Our group previously prepared the 4-, 5-, and 6-day-old larval guts of A. c. cerana followed by construction of strand-specific cDNA libraries and RNA-seq, and the produced raw reads were subjected to strict quality control to generate high-quality clean reads [34,37], which were used to conduct target prediction of DEmiRNAs in this current work. The raw data derived from RNA-seq are available in the NCBI SRA database under the following accession number: SRA456721.

2.6. Prediction and Annotation of DEmiRNA-Targeted mRNAs

Following the procedure outlined by Chen et al. [25], target mRNAs of DEmiRNAs were predicted by utilizing a collaboration of three software, including MiRanda (V3.3a) [38], RNAhybrid (V2.1.2) [39], and TargetFind [40], and the intersection was deemed as reliable targets. Next, these targets were mapped to the Gene Ontology (GO) (http://www.geneontology.org/, accessed on 5 March 2021) and Kyoto Encyclopedia of Genes and Genomes (KEGG) (https://www.kegg.jp, accessed on 5 March 2021) databases to gain corresponding function and pathway annotations by using Blast tool with default parameters, followed by drawing chord diagrams and bubble diagrams with related tools in the OmicShare platform.

2.7. Construction and Investigation of DEmiRNA-DEmRNA Regulatory Networks

Based on the predicted targeting relationships with a threshold of binding free energy less than −15 kcal/mol, regulatory networks between A. c. cerana DEmiRNAs in “Ac4 vs. Ac5” and “Ac5 vs. Ac6” comparison groups and target mRNAs were constructed and then visualized by Cytoscape v.3.2.1 software [41] with default parameters. In order to explore the miRNA-mediated process of growth and development as well as the immune process during the stage of 4-day-old to 6-day-old larvae, following the KEGG database annotations, targets relevant to development-associated signaling pathways as well as cellular and humoral immune pathways were filtered out, followed by construction and visualization of corresponding regulatory networks utilizing the method mentioned above.

2.8. RT-qPCR Verification of DEmiRNAs

The cDNA was synthesized as described in 2.3 for qPCR reaction, which was performed on QuanStudio 3 Fluorescence quantitative PCR instrument (ABI, Los Angeles, CA, USA). The reaction system (20 μL) contained 10 μL of SYBR green dye, 1 μL of forward and reverse primers (2.5 μmol/L), 1 μL of cDNA template, and 7 μL of DEPC water. The reaction conditions were set as follows: 95 °C pre-denaturation for 5 min; 95 °C denaturation for 30 s; 60 °C annealing and extension for 30 s; and a total of 40 cycles each group of qPCR reaction, and the experiment was set for repeating three times. The U6 gene (GenBank accession number: LOC108002035) was used as internal reference. The relative expression level of each DEmiRNA was calculated using the 2−∆∆Ct method [42]. Each experiment was carried out with at least three samples in parallel and was repeated three times. The primers used in this work are detailed in the Table S1.

2.9. Statistical Analysis

Fisher’s exact test was performed to filter out the significant KEGG pathways with the R package 3.3.1 (R Development Core Team, https://www.r-project.org/, accessed on 5 March 2021). The qPCR data are shown as mean ± standard deviation (SD) and subjected to Student’s t-test using SPSS 16.0 (IBM, Armonk, NY, USA) and GraphPad Prism v8.0 software (GraphPad, San Diego, CA, USA, https://www.graphpad.com/features, accessed on 5 May 2022). p < 0.05 was considered statistically significant.

3. Results

3.1. Analysis of miRNAs in the Larval Guts of A. c. cerana

Here, 330, 351, and 321 miRNAs were identified in the 4-, 5-, and 6-day-old larval guts, respectively. Venn analysis suggested that as many as 257 miRNAs were shared by the 4-, 5-, and 6-day-old larval guts, whereas 38, 51, and 36 ones were specific for each group (Figure 1A). In addition, expression clustering analysis showed that the identified A. c. cerana miRNAs presented various trends during the process of development of the larval guts (Figure 1B). Further investigation indicated that one DEmiRNA, namely miR-6001-y, displayed a continuous up-regulation trend in the larval guts, with TPM values of 508.86, 3961.65, and 9390.71, respectively. PCR amplification and Sanger sequencing of six randomly selected miRNAs demonstrated that their sequences were consistent with the predicted sequences based on sRNA-seq data, which verified the authenticity of the A. c. cerana miRNAs identified in this work (Figure 1C).

3.2. Structural Characteristics of miRNAs in 4-, 5-, and 6-Day-Old Larval Guts

The length distribution of miRNAs in Ac4, Ac5, and Ac6 groups ranged from 18 nt to 25 nt; the most enriched length was 22 nt, followed by 23 nt and 18 nt (Figure 2A). In addition, for miRNAs with length ranging from 18 to 25 nt, the first base is always biased toward U (Figure 2B), and miRNAs in the aforementioned three groups had similar bias at each base (Figure 2C).

3.3. Expression Pattern of miRNAs during the Development Process of the A. c. cerana Larval Guts

In the Ac4 vs. Ac5 comparison group, seven up-regulated and eight down-regulated miRNAs were detected; of these, the three most up-regulated miRNAs were miR-1672-x (log2FC = 8.41, p < 0.01), novel-m0032-5p (log2FC = 6.79, p < 0.01), and miR-6001-y (log2FC = 2.96, p < 0.05), while the three most down-regulated ones were novel-m0016-5p (log2FC = −9.25, p < 0.01), novel-m0016-5p (log2FC = −9.25, p < 0.01), and miR-6941-x (log2FC = −3.74, p < 0.05), respectively (Figure 3A). Comparatively, there were four up-regulated and six down-regulated miRNAs in the Ac5 vs. Ac6 comparison group (Figure 3B), and of these, the most up-regulated miRNA was miR-342-y (log2FC = 3.38, p < 0.05), followed by miR-6001-y (log2FC = 1.25, p < 0.05) and novel-m0005-3p (log2FC = 1.16, p < 0.05), whereas the most down-regulated miRNA was miR-26-x (log2FC = −9.41, p < 0.01), followed by miR-21-x (log2FC = −8.73, p < 0.01) and miR-1-z (log2FC = −2.32, p < 0.05) (Figure 3B).

3.4. Investigation and Annotation of DEmiRNA-Targeted mRNAs

DEmiRNAs in the Ac4 vs. Ac5 comparison group were predicted to target 5041 mRNAs, which were annotated to 18 biological-process-associated GO terms such as cellular process (any process that is carried out at the cellular level, but not necessarily restricted to a single), metabolic process (the chemical reactions and pathways, including anabolism and catabolism, by which living organisms transform chemical substances), and biological regulation (any process that modulates a measurable attribute of any biological process, quality or function); 15 cellular component-associated terms such as cell, cell part (refers to various cellular components), and membrane; and 11 molecular-function-associated terms such as binding (refers to binding relationship between different biological molecules), catalytic activity (catalysis of a biochemical reaction at physiological temperatures), and transporter activity (enables the directed movement of substances into, out of, or within a cell or between cells). In contrast, DEmiRNAs within the Ac5 vs. Ac6 comparison group could target 3727 mRNAs, which were annotated to 18 biological-process-related functional terms such as cellular process and metabolic process, 17 cellular-component-associated terms such as cell part and membrane, and 11 molecular-function-associated terms such as catalytic activity and transporter activity (Figure 4).
Additionally, the targets in the Ac4 vs. Ac5 comparison group were shown to be engaged in 130 KEGG pathways associated with metabolism, organismal systems, cellular processes, environmental information processing, and genetic information processing, including the biosynthesis of amino acids, dorso-ventral axis formation, endocytosis, phosphatidylinositol signaling system, and RNA degradation (Figure 5A). In the Ac5 vs. Ac6 comparison group, the targets were implicated in 107 pathways such as the biosynthesis of amino acids, dorso-ventral axis formation, endocytosis, Jak-STAT signaling pathway, and RNA degradation (Figure 5B).

3.5. DEmiRNA-mRNA Regulatory Network Engaged in the Development of A. c. cerana Larval Guts

In the Ac4 vs. Ac5 comparison group, novel-m0032-5p had the highest number of targets (2743), followed by miR-6001-y (2042 targets) and novel-m0049-3p (630 targets); among these targets, each of the 12 mRNA (XM_017052480.1, XM_017052481.1, XM_017052482.1, XM_017052624.1, XM_017052627.1, XM_017052629.1, XM_017054924.1, XM_017058051.1, XM_017062539.1, XM_017062547.1, XM_017062549.1, and XM_017062558.1) was targeted by the five DEmiRNAs. Additionally, miR-6001-y in the Ac5 vs. Ac6 comparison group potentially targeted the most mRNAs (2042), followed by miR-26-x (599 mRNAs) and miR-1-z (466 mRNAs); among these targets, each of the seven mRNA (XM_017048298.1, XM_017048300.1, XM_017057380.1, XM_017057382.1, XM_017064695.1, XM_017064696.1, and XM_017064697.1) was detected to be targeted by six DEmiRNAs.
Further analysis demonstrated that nine DEmiRNAs in the Ac4 vs. Ac5 comparison group could target 43 mRNAs annotated in development-associated signaling pathways including Wnt, Hippo, and Notch (Figure 6A). In contrast, eight DEmiRNAs in the Ac5 vs. Ac6 comparison group could target 31 mRNAs involved in the Wnt, Hippo, and Notch signaling pathways (Figure 6B).

3.6. DEmiRNA-mRNA Regulatory Network Involved in the Cellular and Humoral Immune of A. c. cerana Larval Guts

In the Ac4 vs. Ac5 comparison group, six DEmiRNAs could target four mRNAs annotated in the humoral immune-related signaling pathway Jak-STAT, and nine DEmiRNAs could target 61 mRNAs annotated in cellular immune-related signaling pathways including endocytosis, phagosome, and lysosome (Figure 7A). In the Ac5 vs. Ac6 comparison group, four DEmiRNAs could target seven mRNAs annotated in Jak-STAT, while seven DEmiRNAs could target 34 mRNAs annotated in endocytosis, phagosome, and lysosome (Figure 7B).

3.7. RT-qPCR Detection of DEmiRNAs

RT-qPCR detection showed that the expression trends of five randomly selected DEmiRNAs were in accordance with those in the sequencing results (Figure 8), confirming the trustworthiness of the transcriptome datasets employed in the present study.

4. Discussion

The larval stage of a honey bee lasts until it is six days old. It is hard to dissect intact gut tissues of 1- and 2-day-old bee larvae, which are very small and vulnerable. Hence, 3-day-old larvae were previously removed from the combs and placed in the 24-well culture plates to adapt to the lab conditions for 24 h, and 4-, 5-, and 6-day-old larval guts were selected for sequencing and investigation [27]. Moreover, we previously conducted strand-specific cDNA library-based RNA-seq of 4-, 5-, and 6-day-old larval guts, and the generated datasets included mRNAs and ncRNAs [34,37]. To maintain uniform larval ages in different sequencing projects enables researchers to analyze the interaction relationship between mRNAs and various ncRNAs such as miRNAs, piRNAs, lncRNAs, and circRNAs.
Here, as many as 257 miRNAs were shared in the A. c. cerana Ac4, Ac5, and Ac6 groups, indicating that the majority of miRNAs were co-expressed in 4-, 5-, and 6-day-old larval guts. It is inferred that these shared miRNAs play indispensable roles in the development of larval guts. Additionally, there were 38, 51, and 36 miRNAs specifically expressed in guts of 4-, 5-, and 6-day-old larvae, which is suggestive of their putative functions in the larval guts at different days of the larval stage. Previously, we identified 560 miRNAs in the larval guts of Apis mellifera ligustica, a subspecies of A. mellifera, among which 331 (59%) miRNAs were shared by the 4-, 5-, and 6-day-old larval guts, whereas the quantities of specific miRNAs were 44, 38, and 88, respectively [23]. Together, these results indicated that expression of the majority of miRNAs was essential for the development of the larval guts of both A. c. cerana and A. m. ligustica.
In this current work, following investigation of the structural characteristics, we observed that miRNAs in the 4-, 5-, and 6-day-old larval guts ranged from 18 nt to 25 nt in length, and the first base is always biased toward “U”, which is consistent with miRNAs identified in Equus caballus [43], Bombyx mori [44], and Crassostrea gigas [45]. This was indicative of the high conservation of structural characteristics of miRNAs in various animals and plants. A putative explanation for the “U” bias of first base was that it is more helpful in forming a protein complex with AGO protein, binding to the target gene, and regulating the gene expression [46].
MiRNAs were capable of regulating gene expression by targeting 3′ UTR or the coding sequence of corresponding mRNAs, further modulating numerous life activities [47,48]. In the present study, 15 and 10 DEmiRNAs were predicted to target 5041 and 3727 mRNAs in Ac4 vs. Ac5 and Ac5 vs. Ac6 comparison groups, respectively, involving 44 and 46 GO terms such as binding, cellular process, and catalytic activity as well as 130 and 107 KEGG pathways such as biosynthesis of amino acids, dorso-ventral axis formation, and endocytosis. The results demonstrated that corresponding DEmiRNAs may play regulatory roles in extensive aspects during the developmental process of the A. c. cerana larval guts.
Growth and development of organisms depend on the concerted action of multiple conserved signaling pathways [49,50]. Here, we observed that 10 DEmiRNAs (miR-6001-y, miR-3759-y, miR-3793-x, etc.) in the Ac4 vs. Ac5 comparison group and eight DEmiRNAs (miR-6001-y, miR-1-z, miR-8117-y, etc.) in the Ac5 vs. Ac6 comparison group could target 21 mRNAs and 16 mRNAs relevant to Wnt signaling pathway. In Drosophila, two branches of the Wnt signaling pathway are required for maintenance and regeneration of intestinal stem cells [51,52]. The DEmiRNA-mediated regulation of Wnt in A. c. cerana larvae gut development is expected to be further explored. The Hippo signaling pathway is conserved among various organisms and participates in the regulation of the size of organs [53]. In this current work, seven miRNAs (miR-6001-y, miR-3793-x, novel-m0016-5p, etc.) in the Ac4 vs. Ac5 comparison group and seven miRNAs (miR-6001-y, miR-1-z, miR-8117-y, etc.) in the Ac5 vs. Ac6 comparison group could target 20 and 16 mRNAs associated with the Hippo signaling pathway. The results indicated that corresponding DEmiRNAs may participate in the modulation of the development of larval guts via regulation of Hippo signaling pathways.
Gsk3β, an evolutionarily conserved serine/threonine kinase in all eukaryotes, regulates diverse biological processes from growth and metabolism to cell cycle and apoptosis [54,55]. The Drosophila Shaggy (Sgg) protein is the homologue of glycogen synthase kinase 3β (Gsk3β) [56]. Wu et al. [57] discovered that loss of sgg in Drosophila gave rise to Yorkie inhibition and downregulation of target genes in the Hippo signaling pathway, and sgg regulated the growth of eyes and wings via the Hippo signaling pathway. Here, kinase shaggy-like (XM_017057538.1) was found to be targeted by three DEmiRNAs (miR-6001-y, miR-1638-x, and novel-m0032-5p) in the Ac4 vs. Ac5 comparison group and two DEmiRNAs (miR-6001-y and miR-1-z) in the Ac5 vs. Ac6 comparison group. It demonstrated that Sgg may be continuously regulated by A. c. cerana larvae DEmiRNAs, and furthermore, miR-6001-y might participate in regulating Sgg to mediate growth.
The Son of sevenless (Sos) protein has been suggested to function as a guanine nucleotide transfer factor for Ras and to interact with the receptor tyrosine kinase Sevenless via the Drk protein, a homolog of mammalian Grb2 [58]. Previous studies showed that SOS1-knockout mutants were embryonic lethal [59]. Here, we noticed that the sos (XM_017058556.1) was a common target for miR-6001-x, miR-1638-x, and novel-m0049-3p in the Ac4 vs. Ac5 comparison group as well as miR-1-z in the Ac5 vs. Ac6 comparison group. In summary, these findings suggested that the aforementioned five DEmiRNAs may well regulate the growth and development of A. c. cerana larval guts via regulating sos genes’ expression. However, more efforts are needed to verify their interactions and underlying regulatory mechanisms.
In Bombus terrestris, the expression levels of both Bte-miR-6001-5p and Bte-miR-6001-3p were higher in queen- than in worker-destined late-instar larvae [60]. In A. mellifera, ame-miR-6001-5p has been suggested to play a regulatory part in ecdysone secretion and caste differentiation [61,62]. In this work, miR-6001-y was found to upregulate in both Ac4 vs. Ac5 and Ac5 vs. Ac6 comparison groups, which is suggestive of the activation of miR-6001-y during the development process of larval guts. In addition, we observed that miR-6001-y putatively targeted as many as 2042 mRNAs, involving Wnt, Hippo, and Jak-STAT signaling pathways, suggesting the involvement of miR-6001-y in adjusting these three signaling pathways of great importance. Additional work such as overexpression and knockdown of miR-6001-y in the larval guts via feeding method, following our previously established protocol [36], is required to decipher the function and action mechanism of miR-6001-y.
Different ncRNAs such as lncRNAs and circRNAs, which contain miRNA response elements, are able to interact with miRNAs via competing endogenous networks to further regulate expression of downstream genes [63]. In insects, ncRNA-involved ceRNA networks were shown to affect and development [16], metabolism [64], and immunity [65]. For example, Wang et al. [66] discovered that the nucleus-enriched lncR26319 in silkworm was able to regulate Endophilin A (EndoA), a member of the endophilin family of endocytic proteins, through competitive binding to miR-2834; the lncR26319-miR-2834-EndoA axis was required for Vg endocytosis. Our group previously discovered that 29 up-regulated lncRNAs and 14 up-regulated circRNAs in A. m. ligustica workers’ midguts potentially targeted ame-miR-6001-3p [67,68]. Findings from this work and our previous studies demonstrate that miR-6001 as a hub within ceRNA regulatory networks likely exerts a critical function in modulating the larval guts’ growth and development and workers’ midguts in both A. c. cerana and A. m. ligustica. Recently, we established the platforms for dissecting the functions of miRNA [36] and lncRNA [69] in bee larval guts. Hence, to identify candidate lncRNAs and circRNAs targeted by ace-miR-6001 and investigate the regulatory functions of the lncRNA-ace-miR-6001 axis and circRNA-ace-miR-6001 axis in the process of the development of larval guts is an interesting and promising direction in the near future.

5. Conclusions

In conclusion, the dynamic expression and structural alteration of miRNAs were accompanied by the development of A. c. cerana larval guts; DEmiRNAs were likely to participate in the modulation of the growth and development of larval guts by regulating the expression of target genes via affecting an array of crucial pathways including Wnt and Hippo; miR-6001-y is a potential key regulator in the development of larval guts by ceRNA mechanism or silencing genes. Our data offer a basis for elucidating the developmental mechanism underlying Asian honey bee larval guts.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/insects14050469/s1, Table S1: Primer information.

Author Contributions

R.G. and D.C. designed this research; X.F. and W.Z. contributed to the writing of the article; X.F., W.Z., S.G., L.Z., Y.Z., H.Z., X.G., H.J. and T.Z. conducted experiments and data analyses; R.G., Q.N. and D.C. supervised the study and the preparation of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the National Natural Science Foundation of China (32172792, 31702190), the Earmarked fund for China Agriculture Research System (CARS-44-KXJ7), the Natural Science Foundation of Fujian Province (2022J01131334), the Master Supervisor Team Fund of Fujian Agriculture and Forestry University (Rui Guo), the Special Fund for Science and Technology Innovation of Fujian Agriculture and Forestry University (Rui Guo), the Scientific Research Project of College of Animal Sciences (College of Bee Science) of Fujian Agriculture and Forestry University (Rui Guo), and the Undergraduate Innovation and Entrepreneurship Training Program of Fujian province (202210389128).

Data Availability Statement

All the data are contained within the article.

Acknowledgments

We thank all editors and reviewers for their constructive comments and recommendations.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ruttner, F. Biogeography and Taxonomy of Honeybees; Springer: Berlin, Germany; New York, NY, USA, 2013. [Google Scholar]
  2. Raymann, K.; Moran, N.A. The role of the gut microbiome in health and disease of adult honey bee workers. Curr. Opin. Insect Sci. 2018, 26, 97–104. [Google Scholar] [CrossRef]
  3. Giurfa, M.; Marcout, C.; Hilpert, P.; Thevenot, C.; Rugani, R. An insect brain organizes numbers on a left-to-right mental number line. Proc. Natl. Acad. Sci. USA 2022, 119, e2203584119. [Google Scholar] [CrossRef]
  4. Giurfa, M.; Zhang, S.; Jenett, A.; Menzel, R.; Srinivasan, M.V. The concepts of ‘sameness’ and ‘difference’ in an insect. Nature 2001, 410, 930–933. [Google Scholar] [CrossRef]
  5. Fan, X.; Zhang, W.; Zhang, K.; Zhang, J.; Long, Q.; Wu, Y.; Zhang, K.; Zhu, L.; Chen, D.; Guo, R. In-depth investigation of microRNA-mediated cross-kingdom regulation between Asian honey bee and microsporidian. Front. Microbiol. 2022, 13, 1003294. [Google Scholar] [CrossRef]
  6. Keller, A.; Brandel, A.; Becker, M.C.; Balles, R.; Abdelmohsen, U.R.; Ankenbrand, M.J.; Sickel, W. Wild bees and their nests host Paenibacillus bacteria with functional potential of avail. Microbiome 2018, 6, 229. [Google Scholar] [CrossRef] [PubMed]
  7. Watson, O.T.; Buchmann, G.; Young, P.; Lo, K.; Remnant, E.J.; Yagound, B.; Shambrook, M.; Hill, A.F.; Oldroyd, B.P.; Ashe, A. Abundant small RNAs in the reproductive tissues and eggs of the honey bee, Apis mellifera. BMC Genom. 2022, 23, 257. [Google Scholar] [CrossRef] [PubMed]
  8. Diao, Q.; Sun, L.; Zheng, H.; Zeng, Z.; Wang, S.; Xu, S.; Zheng, H.; Chen, Y.; Shi, Y.; Wang, Y.; et al. Genomic and transcriptomic analysis of the Asian honeybee Apis cerana provides novel insights into honeybee biology. Sci. Rep. 2018, 8, 822. [Google Scholar] [CrossRef]
  9. Winston, M.L. The Biology of the Honey Bee; Harvard University Press: Cambridge, UK, 1991. [Google Scholar]
  10. Hu, X.; Ke, L.; Wang, Z.; Zeng, Z. Dynamic transcriptome landscape of Asian domestic honeybee (Apis cerana) embryonic development revealed by high-quality RNA sequencing. BMC Dev. Biol. 2018, 18, 11. [Google Scholar] [CrossRef]
  11. Ai, H.; Yan, X.; Han, R. Occurrence and prevalence of seven bee viruses in Apis mellifera and Apis cerana apiaries in china. J. Invertebr. Pathol. 2012, 109, 160–164. [Google Scholar] [CrossRef]
  12. Zhang, Q.; Dou, W.; Taning, C.N.T.; Smagghe, G.; Wang, J.J. Regulatory roles of microRNAs in insect pests: Prospective targets for insect pest control. Curr. Opin. Biotechnol. 2021, 70, 158–166. [Google Scholar] [CrossRef] [PubMed]
  13. Yang, Y.; Wang, A.; Zhang, Y.; Xue, C.; Zhao, M.; Zhang, J. Activating pathway of three metabolic detoxification phases via down-regulated endogenous microRNAs, modulates triflumezopyrim tolerance in the small brown planthopper, Laodelphax striatellus (Fallén). Int. J. Biol. Macromol. 2022, 222, 2439–2451. [Google Scholar] [CrossRef] [PubMed]
  14. Yang, M.; Du, B.; Xu, L.; Wang, H.; Wang, Y.; Lin, K.; He, G.; Kang, L. Glutamate-GABA imbalance mediated by miR-8-5p and its STTM regulates phase-related behavior of locusts. Proc. Natl. Acad. Sci. USA 2023, 120, e2215660120. [Google Scholar] [CrossRef] [PubMed]
  15. Shen, Z.J.; Zhu, F.; Liu, Y.J.; Li, Z.; Moural, T.W.; Liu, X.M.; Liu, X. MicroRNAs miR-14 and miR-2766 regulate tyrosine hydroxylase to control larval-pupal metamorphosis in Helicoverpa armigera. Pest Manag. Sci. 2022, 78, 3540–3550. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, Y.; Zou, H.; Lai, J.; Zhang, Z.; Sun, W. The miR-282-5p regulates larval moulting process by targeting chitinase 5 in Bombyx mori. Insect Mol. Biol. 2022, 31, 190–201. [Google Scholar] [CrossRef]
  17. Zhang, Q.; Dou, W.; Taning, C.N.T.; Yu, S.S.; Yuan, G.R.; Shang, F.; Smagghe, G.; Wang, J.J. miR-309a is a regulator of ovarian development in the oriental fruit fly Bactrocera dorsalis. PLoS Genet. 2022, 18, e1010411. [Google Scholar] [CrossRef]
  18. Freitas, F.C.; Pires, C.V.; Claudianos, C.; Cristino, A.S.; Simões, Z.L. MicroRNA-34 directly targets pair-rule genes and cytoskeleton component in the honey bee. Sci. Rep. 2017, 7, 40884. [Google Scholar] [CrossRef]
  19. Chen, X.; Ma, C.; Chen, C.; Lu, Q.; Shi, W.; Liu, Z.; Wang, H.; Guo, H. Integration of lncRNA-miRNA-mRNA reveals novel insights into oviposition regulation in honey bees. PeerJ 2017, 5, e3881. [Google Scholar] [CrossRef]
  20. Guo, X.; Su, S.; Skogerboe, G.; Dai, S.; Li, W.; Li, Z.; Liu, F.; Ni, R.; Guo, Y.; Chen, S.; et al. Recipe for a busy bee: MicroRNAs in honey bee caste determination. PLoS ONE 2013, 8, e81661. [Google Scholar] [CrossRef]
  21. Ma, L.; Liu, L.; Zhao, Y.; Yang, L.; Chen, C.; Li, Z.; Lu, Z. JNK pathway plays a key role in the immune system of the pea aphid and is regulated by microRNA-184. PLoS Pathog. 2020, 16, e1008627. [Google Scholar] [CrossRef]
  22. Huang, D.Y.; Xia, X.L.; Huang, R.; Li, S.; Yuan, D.W.; Liu, S.N. The steroid-induced microRNA let-7 regulates developmental growth by targeting cdc7 in the Drosophila fat body. Insect Sci. 2021, 28, 1621–1632. [Google Scholar] [CrossRef]
  23. Xiong, C.L.; Du, Y.; Chen, D.F.; Zheng, Y.Z.; Fu, Z.M.; Wang, H.P.; Geng, S.H.; Chen, H.Z.; Zhou, D.D.; Wu, S.Z.; et al. Bioinformatic prediction and analysis of miRNAs in the Apis mellifera ligustica larval gut. Chin. J. Appl. Entomol. 2018, 55, 1023–1033. (In Chinese) [Google Scholar]
  24. Guo, R.; Du, Y.; Xiong, C.L.; Zheng, Y.Z.; Fu, Z.M.; Xu, G.J.; Wang, H.P.; Chen, H.Z.; Geng, S.H.; Zhou, D.D.; et al. Differentially expressed microRNA and their regulation networks during the developmental process of Apis mellifera ligustica larval gut. Sci. Agric. Sin. 2018, 51, 4197–4209. (In Chinese) [Google Scholar]
  25. Chen, D.; Du, Y.; Chen, H.; Fan, Y.; Fan, X.; Zhu, Z.; Wang, J.; Xiong, C.; Zheng, Y.; Hou, C.; et al. Comparative identification of microRNAs in Apis cerana cerana workers’ midguts in responseto Nosema ceranae invasion. Insects 2019, 10, 258. [Google Scholar] [CrossRef] [PubMed]
  26. Deng, Y.; Zhao, H.; Shen, S.; Yang, S.; Yang, D.; Deng, S.; Hou, C. Identification of immune response to sacbrood virus infection in Apis cerana under natural condition. Front. Genet. 2020, 11, 587509. [Google Scholar] [CrossRef]
  27. Feng, R.R.; Fu, Z.M.; Du, Y.; Zhang, W.D.; Fan, X.X.; Wang, H.P.; Wan, J.Q.; Zhou, Z.Y.; Kang, Y.X.; Guo, R.; et al. Identification and analysis of micrornas in the larval gut of Apis cerana cerana. Sci. Agric. Sin. 2022, 55, 208–218. (In Chinese) [Google Scholar]
  28. Sharon, G.; Segal, D.; Ringo, J.M.; Hefetz, A.; Zilber-Rosenberg, I.; Rosenberg, E. Commensal bacteria play a role in mating preference of Drosophila melanogaster. Proc. Natl. Acad. Sci. USA 2010, 107, 20051–20056. [Google Scholar] [CrossRef]
  29. Zheng, H.; Steele, M.I.; Leonard, S.P.; Motta, E.V.S.; Moran, N.A. Honey bees as models for gut microbiota research. Lab Anim. 2018, 47, 317–325. [Google Scholar] [CrossRef]
  30. Dosch, C.; Manigk, A.; Streicher, T.; Tehel, A.; Paxton, R.J.; Tragust, S. The gut microbiota can provide viral tolerance in the honey bee. Microorganisms 2021, 9, 871. [Google Scholar] [CrossRef]
  31. Pino, A.; Benkaddour, B.; Inturri, R.; Amico, P.; Vaccaro, S.C.; Russo, N.; Vaccalluzzo, A.; Agolino, G.; Caggia, C.; Miloud, H.; et al. Characterization of Bifidobacterium asteroides Isolates. Microorganisms 2022, 10, 655. [Google Scholar] [CrossRef]
  32. Long, Q.; Sun, M.H.; Fan, X.X.; Cai, Z.B.; Zhang, K.Y.; Wang, S.Y.; Zhang, J.X.; Gu, X.Y.; Song, Y.X.; Chen, D.F.; et al. First identification and investigation of piRNAs in the larval gut of the Asian honeybee, Apis cerana. Insects 2022, 14, 16. [Google Scholar] [CrossRef] [PubMed]
  33. Xu, Y.J.; Long, Q.; Fan, X.X.; Ye, Y.P.; Zhang, K.Y.; Zhang, J.X.; Zhao, H.D.; Yao, Y.T.; Fu, Z.M.; Chen, D.F.; et al. Transcriptome-wide characterization of piRNAs during the developmental process of european honey-bee larval guts. Genes 2022, 13, 1879. [Google Scholar] [CrossRef]
  34. Xu, X.J.; Guo, R.; Luo, Q.; Xiong, C.L.; Liang, Q.; Zhang, C.L.; Zheng, Y.Z.; Zhang, Z.N.; Huang, Z.J.; Zhang, L.; et al. De novo transcriptome assembly for Apis cerana cerana larval gut and identification of SSR molecular markers. Sci. Agric. Sin. 2017, 50, 1157–1167. [Google Scholar]
  35. Singh, V.K.; Mangalam, A.K.; Dwivedi, S.; Naik, S. Primer premier: Program for design of degenerate primers from a protein sequence. BioTechniques 1998, 24, 318–319. [Google Scholar] [CrossRef] [PubMed]
  36. Wu, Y.; Guo, Y.; Fan, X.; Zhao, H.; Zhang, Y.; Guo, S.; Jing, X.; Liu, Z.; Feng, P.; Liu, X.; et al. ame-miR-34 modulates the larval body weight and immune response of Apis mellifera workers to Ascosphara apis invasion. Int. J. Mol. Sci. 2023, 24, 1214. [Google Scholar] [CrossRef] [PubMed]
  37. Guo, R.; Chen, D.; Diao, Q.; Xiong, C.; Zheng, Y.; Hou, C. Transcriptomic investigation of immune responses of the Apis cerana cerana larval gut infected by Ascosphaera apis. J. Invertebr. Pathol. 2019, 166, 107210. [Google Scholar] [CrossRef] [PubMed]
  38. Enright, A.J.; John, B.; Gaul, U.; Tuschl, T.; Sander, C.; Marks, D.S. MicroRNA targets in Drosophila. Genome Biol. 2003, 5, R1. [Google Scholar] [CrossRef] [PubMed]
  39. Krüger, J.; Rehmsmeier, M. RNAhybrid: MicroRNA target prediction easy, fast and flexible. Nucleic Acids Res. 2006, 34, W451–W454. [Google Scholar] [CrossRef]
  40. Allen, E.; Xie, Z.; Gustafson, A.M.; Carrington, J.C. microRNA-directed phasing during trans-acting siRNA biogenesis in plants. Cell 2005, 121, 207–221. [Google Scholar] [CrossRef]
  41. Smoot, M.E.; Ono, K.; Ruscheinski, J.; Wang, P.L.; Ideker, T. Cytoscape 2.8: New features for data integration and network visualization. Bioinformatics 2011, 27, 431–432. [Google Scholar] [CrossRef]
  42. Livak, K.J.; Schmittgen, T.D. Anaysis of relative gene expression data using realtime quantitative PCR and the 2−ΔΔCt method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  43. Kim, M.C.; Lee, S.W.; Ryu, D.Y.; Cui, F.J.; Bhak, J.; Kim, Y. Identification and characterization of microRNAs in normal equine tissues by next generation sequencing. PLoS ONE 2014, 9, e93662. [Google Scholar] [CrossRef]
  44. Pan, Q.L.; Li, T.; He, Q.; Ma, Z.G.; Fan, X.D.; Zhang, X.Y.; Wang, Y.L.; Zhou, Z.Y.; Xu, J.S. Characterization of small RNAs in microsporidian Nosema bombycis. Acta Entomol. Sin. 2015, 58, 1213–1221. [Google Scholar] [CrossRef]
  45. Wang, X.; Wang, W.J.; Luo, Q.H.; Sun, G.H.; Feng, Y.W.; Ma, J.J.; Yang, J.M. Bioinformatics analysis of regulatory non-coding RNA in gonad of Crassostrea gigas. J. Fish. China 2020, 44, 723–734. [Google Scholar] [CrossRef]
  46. Mi, S.; Cai, T.; Hu, Y.; Chen, Y.; Hodges, E.; Ni, F.; Wu, L.; Li, S.; Zhou, H.; Long, C.; et al. Sorting of small RNAs into Arabidopsis argonaute complexes is directed by the 5′ terminal nucleotide. Cell 2008, 133, 116–127. [Google Scholar] [CrossRef] [PubMed]
  47. Shabalina, S.A.; Koonin, E.V. Origins and evolution of eukaryotic RNA interference. Trends Ecol. Evol. 2008, 23, 578–587. [Google Scholar] [CrossRef]
  48. Eichhorn, S.W.; Guo, H.; McGeary, S.E.; Rodriguez-Mias, R.A.; Shin, C.; Baek, D.; Hsu, S.H.; Ghoshal, K.; Villén, J.; Bartel, D.P. mRNA destabilization is the dominant effect of mammalian microRNAs by the time substantial repression ensues. Mol. Cell 2014, 56, 104–115. [Google Scholar] [CrossRef]
  49. Steinhart, Z.; Angers, S. Wnt signaling in development and tissue homeostasis. Development 2018, 145, dev146589. [Google Scholar] [CrossRef]
  50. Choi, K.W. Upstream paths for Hippo signaling in Drosophila organ development. BMB Rep. 2018, 51, 134–142. [Google Scholar] [CrossRef]
  51. Cordero, J.B.; Sansom, O.J. Wnt signalling and its role in stem cell-driven intestinal regeneration and hyperplasia. Acta Physiol. 2012, 204, 137–143. [Google Scholar] [CrossRef] [PubMed]
  52. Mah, A.T.; Yan, K.S.; Kuo, C.J. Wnt pathway regulation of intestinal stem cells. J. Physiol. 2016, 594, 4837–4847. [Google Scholar] [CrossRef]
  53. Russell, J.O.; Camargo, F.D. Hippo signalling in the liver: Role in development, regeneration and disease. Nat. Rev. Gastroenterol. Hepatol. 2022, 19, 297–312. [Google Scholar] [CrossRef] [PubMed]
  54. Patel, P.; Woodgett, J.R. Glycogen synthase kinase 3: A kinase for all pathways? Curr. Top. Dev. Biol. 2017, 123, 277–302. [Google Scholar] [CrossRef]
  55. He, L.; Feng, H.; Yin, B.; Li, W.; Wang, X.; Umar, T.; Gao, H.; Zhou, N.; Qiu, C. Sodium new houttuyfonate induces apoptosis of breast cancer cells via ROS/PDK1/AKT/GSK3β Axis. Cancers 2023, 15, 1614. [Google Scholar] [CrossRef]
  56. Kanuka, H.; Kuranaga, E.; Takemoto, K.; Hiratou, T.; Okano, H.; Miura, M. Drosophila caspase transduces Shaggy/GSK-3β kinase activity in neural precursor development. EMBO J. 2005, 24, 3793–3806. [Google Scholar] [CrossRef]
  57. Wu, H.; Zhu, N.; Liu, J.; Ma, J.; Jiao, R. Shaggy regulates tissue growth through Hippo pathway in Drosophila. Sci. China Life Sci. 2022, 65, 2131–2144. [Google Scholar] [CrossRef] [PubMed]
  58. Karlovich, C.A.; Bonfini, L.; McCollam, L.; Rogge, R.D.; Daga, A.; Czech, M.P.; Banerjee, U. In vivo functional analysis of the Ras exchange factor son of sevenless. Science 1995, 268, 576–579. [Google Scholar] [CrossRef] [PubMed]
  59. Baltanás, F.C.; García-Navas, R.; Santos, E. SOS2 Comes to the Fore: Differential Functionalities in Physiology and Pathology. Int. J. Mol. Sci. 2021, 22, 6613. [Google Scholar] [CrossRef]
  60. Collins, D.H.; Mohorianu, I.; Beckers, M.; Moulton, V.; Dalmay, T.; Bourke, A.F. MicroRNAs associated with caste determination and differentiation in a primitively eusocial insect. Sci. Rep. 2017, 7, 45674. [Google Scholar] [CrossRef]
  61. Shi, Y.Y.; Zheng, H.J.; Pan, Q.Z.; Wang, Z.L.; Zeng, Z.J. Differentially expressed microRNAs between queen and worker larvae of the honey bee (Apis mellifera). Apidologie 2015, 46, 35–45. [Google Scholar] [CrossRef]
  62. Ashby, R.; Foret, S.; Searle, I.; Maleszka, R. MicroRNAs in honey bee caste determination. Sci. Rep. 2016, 6, 18794. [Google Scholar] [CrossRef] [PubMed]
  63. Salmena, L.; Poliseno, L.; Tay, Y.; Kats, L.; Pandolfi, P.P. A ceRNA hypothesis: The Rosetta Stone of a hidden RNA language? Cell 2011, 146, 353–358. [Google Scholar] [CrossRef] [PubMed]
  64. Liu, W.; An, S.; Cheng, P.; Zhang, K.; Gong, M.; Zhang, Z.; Zhang, R. Whole-transcriptome profiling across different developmental stages of Aedes albopictus (Diptera: Culicidae) provides insights into chitin-related non-coding RNA and competing endogenous RNA networks. Parasites Vectors 2023, 16, 33. [Google Scholar] [CrossRef] [PubMed]
  65. Mao, W.; Zeng, Q.; She, L.; Yuan, H.; Luo, Y.; Wang, R.; She, Y.; Wang, W.; Wang, C.; Pan, X. Wolbachia Utilizes lncRNAs to Activate the Anti-Dengue Toll Pathway and Balance Reactive Oxygen Species Stress in Aedes aegypti Through a Competitive Endogenous RNA Network. Front. Cell. Infect. Microbiol. 2022, 11, 823403. [Google Scholar] [CrossRef]
  66. Wang, Y.; Fu, Y.; Cheng, H.; Zhao, C.; Huang, Q.; Chang, M.; Qiu, W.; Shen, Y.; Li, D. lncR26319/miR-2834/EndophilinA axis regulates oogenesis of the silkworm, Bombyx mori. Insect Sci. 2023, 30, 65–80. [Google Scholar] [CrossRef] [PubMed]
  67. Guo, R.; Chen, H.Z.; Xiong, C.L.; Zheng, Y.Z.; Fu, Z.M.; Xu, G.J.; Du, Y.; Wang, H.P.; Geng, S.H.; Zhou, D.D.; et al. Analysis of differentially expressed circular RNAs and their regulation networks during the developmental process of Apis mellifera ligustica worker’s midgut. Sci. Agric. Sin. 2018, 51, 4575–4590. (In Chinese) [Google Scholar]
  68. Guo, R.; Geng, S.H.; Xiong, C.L.; Zheng, Y.Z.; Fu, Z.M.; Wang, H.P.; Du, Y.; Tong, X.Y.; Zhao, H.X.; Chen, D.F. Differential expression analysis of long non-coding RNAs during the developmental process of Apis mellifera ligustica worker’s midgut. Sci. Agric. Sin. 2018, 51, 3600–3613. (In Chinese) [Google Scholar]
  69. Fu, Z.M.; Gu, X.Y.; Hu, Y.; Zhao, H.D.; Zhu, Z.W.; Zhang, H.Y.; Ji, T.; Niu, Q.S.; Chen, D.F.; Guo, R. Lnc13164 regulates immune response of Apis cerana cerana larvae to Ascosphaera apis infection via ace-miR-4968-y. Acta Microbiol. Sin. 2023, 63, 1047–1059. [Google Scholar]
Insects 14 00469 g001 550
Figure 1. Expression analysis and molecular validation of miRNAs identified in the A. c. cerana 4-, 5-, and 6-day-old larval guts. (A) Venn analysis of miRNAs discovered in the 4-, 5-, and 6-day-old larval guts; (B) expression clustering of miRNAs expressed in the 4-, 5-, and 6-day-old larval guts; (C) Sanger sequencing of six miRNAs shared by the 4-, 5-, and 6-day-old larval guts.
Figure 1. Expression analysis and molecular validation of miRNAs identified in the A. c. cerana 4-, 5-, and 6-day-old larval guts. (A) Venn analysis of miRNAs discovered in the 4-, 5-, and 6-day-old larval guts; (B) expression clustering of miRNAs expressed in the 4-, 5-, and 6-day-old larval guts; (C) Sanger sequencing of six miRNAs shared by the 4-, 5-, and 6-day-old larval guts.
Insects 14 00469 g001
Insects 14 00469 g002 550
Figure 2. Structural characteristics of miRNAs in A. c. cerana 4-, 5-, and 6-day-old larval guts. (A) Length distribution of miRNAs; (B) first nucleotide bias of miRNAs; (C) bias of each nucleotide of miRNAs.
Figure 2. Structural characteristics of miRNAs in A. c. cerana 4-, 5-, and 6-day-old larval guts. (A) Length distribution of miRNAs; (B) first nucleotide bias of miRNAs; (C) bias of each nucleotide of miRNAs.
Insects 14 00469 g002
Insects 14 00469 g003 550
Figure 3. Radar maps of DEmiRNAs in the Ac4 vs. Ac5 (A) and Ac5 vs. Ac6 (B) comparison groups. Blue circles represent downregulation, while orange circles represent upregulation. The larger the circles, the greater the difference.
Figure 3. Radar maps of DEmiRNAs in the Ac4 vs. Ac5 (A) and Ac5 vs. Ac6 (B) comparison groups. Blue circles represent downregulation, while orange circles represent upregulation. The larger the circles, the greater the difference.
Insects 14 00469 g003
Insects 14 00469 g004 550
Figure 4. GO terms annotated by the DEmiRNA-targeted mRNAs in the Ac4 vs. Ac5 and Ac5 vs. Ac6 comparison groups. Different colors indicate different functional terms. The numbers inside the chord diagram indicate target mRNAs annotated to corresponding terms.
Figure 4. GO terms annotated by the DEmiRNA-targeted mRNAs in the Ac4 vs. Ac5 and Ac5 vs. Ac6 comparison groups. Different colors indicate different functional terms. The numbers inside the chord diagram indicate target mRNAs annotated to corresponding terms.
Insects 14 00469 g004
Insects 14 00469 g005 550
Figure 5. KEGG pathways annotated by the targets of DEmiRNAs in the Ac4 vs. Ac5 (A) and Ac5 vs. Ac6 (B) comparison groups. Bubbles indicate target mRNAs enriched in corresponding pathways; the larger the bubbles, the greater the amount of target mRNAs, and the smaller the bubbles, the fewer the target mRNAs. Different colors represent the p-values of different pathways; the red color indicates higher p-values, while the purple color indicates lower p-values.
Figure 5. KEGG pathways annotated by the targets of DEmiRNAs in the Ac4 vs. Ac5 (A) and Ac5 vs. Ac6 (B) comparison groups. Bubbles indicate target mRNAs enriched in corresponding pathways; the larger the bubbles, the greater the amount of target mRNAs, and the smaller the bubbles, the fewer the target mRNAs. Different colors represent the p-values of different pathways; the red color indicates higher p-values, while the purple color indicates lower p-values.
Insects 14 00469 g005
Insects 14 00469 g006 550
Figure 6. DEmiRNA–mRNA networks relevant to development-associated signaling pathways in the Ac4 vs. Ac5 (A) and Ac5 vs. Ac6 (B) comparison groups. Orange diamonds represent target mRNAs, while purple circles represent DEmiRNAs. Grey lines indicate potential targeting relationships between DEmiRNAs and target mRNAs.
Figure 6. DEmiRNA–mRNA networks relevant to development-associated signaling pathways in the Ac4 vs. Ac5 (A) and Ac5 vs. Ac6 (B) comparison groups. Orange diamonds represent target mRNAs, while purple circles represent DEmiRNAs. Grey lines indicate potential targeting relationships between DEmiRNAs and target mRNAs.
Insects 14 00469 g006
Insects 14 00469 g007 550
Figure 7. DEmiRNA-mRNA networks involved in the cellular and humoral immune in the Ac4 vs. Ac5 (A) and Ac4 vs. Ac5 (B) comparison groups. Orange circles represent target mRNAs, red circles represent miRNAs, and green circles represent annotated pathway. Grey lines indicate potential targeting relationships between the miRNAs–mRNAs pathway.
Figure 7. DEmiRNA-mRNA networks involved in the cellular and humoral immune in the Ac4 vs. Ac5 (A) and Ac4 vs. Ac5 (B) comparison groups. Orange circles represent target mRNAs, red circles represent miRNAs, and green circles represent annotated pathway. Grey lines indicate potential targeting relationships between the miRNAs–mRNAs pathway.
Insects 14 00469 g007
Insects 14 00469 g008 550
Figure 8. RT-qPCR validation of five DEmiRNAs. Both qPCR data and RNA-seq data were presented as mean ± standard deviation (SD) and subjected to Student’s t test, * p < 0.05; ** p < 0.01.
Figure 8. RT-qPCR validation of five DEmiRNAs. Both qPCR data and RNA-seq data were presented as mean ± standard deviation (SD) and subjected to Student’s t test, * p < 0.05; ** p < 0.01.
Insects 14 00469 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fan, X.; Zhang, W.; Guo, S.; Zhu, L.; Zhang, Y.; Zhao, H.; Gao, X.; Jiang, H.; Zhang, T.; Chen, D.; et al. Expression Profile, Regulatory Network, and Putative Role of microRNAs in the Developmental Process of Asian Honey Bee Larval Guts. Insects 2023, 14, 469. https://doi.org/10.3390/insects14050469

AMA Style

Fan X, Zhang W, Guo S, Zhu L, Zhang Y, Zhao H, Gao X, Jiang H, Zhang T, Chen D, et al. Expression Profile, Regulatory Network, and Putative Role of microRNAs in the Developmental Process of Asian Honey Bee Larval Guts. Insects. 2023; 14(5):469. https://doi.org/10.3390/insects14050469

Chicago/Turabian Style

Fan, Xiaoxue, Wende Zhang, Sijia Guo, Leran Zhu, Yiqiong Zhang, Haodong Zhao, Xuze Gao, Haibin Jiang, Tianze Zhang, Dafu Chen, and et al. 2023. "Expression Profile, Regulatory Network, and Putative Role of microRNAs in the Developmental Process of Asian Honey Bee Larval Guts" Insects 14, no. 5: 469. https://doi.org/10.3390/insects14050469

APA Style

Fan, X., Zhang, W., Guo, S., Zhu, L., Zhang, Y., Zhao, H., Gao, X., Jiang, H., Zhang, T., Chen, D., Guo, R., & Niu, Q. (2023). Expression Profile, Regulatory Network, and Putative Role of microRNAs in the Developmental Process of Asian Honey Bee Larval Guts. Insects, 14(5), 469. https://doi.org/10.3390/insects14050469

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop