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

Transcriptional Identification of Differentially Expressed Genes Associated with Division of Labor in Apis Cerana Cerana

Differential Transcriptome Profiles of Heads from Foragers: Comparison Between Apis Mellifera Ligustica and Apis Cerana Cerana

Molecular mechanisms of task allocation in workers of the red imported fire ant, Solenopsis invicta

Gene expression and epigenetics reveal species-specific mechanisms acting upon common molecular pathways in the evolution of task division in bees

Differential protein expression in honeybee (Apis mellifera L.) larvae: underlying caste differentiation

Transcriptome comparison between honey bee queen- and worker-destined larvae

The Diverging Epigenomic Landscapes of Honeybee Queens and Workers Revealed by Multiomic Sequencing.

Insulin signaling is involved in the regulation of worker division of labor in honey bee colonies

Identification of Genes Related to High Royal Jelly Production in the Honey Bee (apis Mellifera) Using Microarray Analysis

Full-Length Transcriptome Profile of Apis cerana Revealed by Nanopore Sequencing

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

Genome-Wide Differential DNA Methylation in Reproductive, Morphological, and Visual System Differences Between Queen Bee and Worker Bee (Apis mellifera)

Changes in gut microbiota and metabolism associated with phenotypic plasticity in the honey bee Apis mellifera

Transcriptome Comparison Between Inactivated And Activated Ovaries Of The Honey Bee Apis Mellifera L.

Analysis of Highly Expressed Genes of Larval Gut of Apis Cerana Cerana

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

Extent and complexity of RNA processing in the development of honey bee queen and worker castes revealed by Nanopore direct RNA sequencing

Transcriptome comparison between inactivated and activated ovaries of the honey bee<i>Apis mellifera</i> L

Sex-specific molecular specialization and activity rhythm dependent gene expression changes in honey bee antennae

Sex-specific molecular specialization and activity rhythm-dependent gene expression in honey bee antennae

A comparison of honeybee (Apis mellifera) queen, worker and drone larvae by RNA-Seq