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.

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

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

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

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

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

Differential Expression of Mirnas Related to Caste Differentiation in the Honey Bee, Apis Mellifera

A Global Survey of the Full-Length Transcriptome of Apis mellifera by Single-Molecule Long-Read Sequencing

Genome‐wide Analysis Reveals Differences in Brain Gene Expression Patterns Associated with Caste and Reproductive Status in Honey Bees (apis Mellifera)

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

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

RNA m6A Modification Functions in Larval Development and Caste Differentiation in Honeybee (Apis mellifera)

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

Intronic Non-CG DNA Hydroxymethylation and Alternative Mrna Splicing in Honey Bees

Recipe for A Busy Bee: Micrornas in Honey Bee Caste Determination

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

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

Next-generation small RNA sequencing for microRNAs profiling in the honey bee Apis mellifera.

Intragenomic conflict underlies extreme phenotypic plasticity in queen-worker caste determination in honey bees ( )

Next-Generation Small Rna Sequencing For Micrornas Profiling In Apis Mellifera: Comparison Between Nurses And Foragers

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

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