Sequence and expression characteristics of long noncoding RNAs in honey bee caste development--potential novel regulators for transgressive ovary size

In silico identification and functional study of long non-coding RNA involved in acute hepatopancreatic necrosis disease caused by Vibrio parahaemolyticus infection in white shrimp, Litopenaeus vannamei

Acute hepatopancreatic necrosis disease (AHPND) caused by toxin-producing Vibrio parahaemolyticus (VpAHPND) has severely affected shrimp production. Long non-coding RNA (lncRNA), a regulatory non-coding RNA, which can play important function in shrimp disease responses. This study aimed to identify and investigate the role of lncRNA involved in VpAHPND infection in Pacific white shrimp, Litopenaeus vannamei. From a total of 368,736 de novo assembled transcripts, 67,559 were identified as putative lncRNAs, and only 72 putative lncRNAs showed differential expression between VpAHPND-infected and normal shrimp. The six candidate lncRNAs were validated for their expression profiles during VpAHPND infection and tissue distribution using RT-qPCR. The role of lnc2088 in response to VpAHPND infection was investigated through RNA interference. The result indicated that the suppression of lnc2088 expression led to an increase in shrimp mortality after VpAHPND infection. To explore the set of genes involved in lnc2088 knockdown, RNA sequencing was performed. A total of 275 differentially expressed transcripts were identified in the hepatopancreas of lnc2088 knockdown shrimp. The expression profiles of five candidate metabolic and immune-related genes were validated in lnc2088 knockdown and VpAHPND-infected shrimp. The result showed that the expression of ChiNAG was significantly increased, while that of NCBP1, WIPF2, and NFKB1 was significantly downregulated in ds2088-injected shrimp. Additionally, the expression of NFKB1, NCBP1 and WIPF2 was significantly increased, whereas that of ChiNAG and CUL5 were significantly decreased after infection with VpAHPND. Our work identified putative lncRNA profiles in L. vannamei in response to VpAHPND infection and investigated the role of lncRNA in shrimp immunity.

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.

Untangling the gordian knot: The intertwining interactions between developmental hormone signaling and epigenetic mechanisms in insects

Hormones control developmental and physiological processes, often by regulating the expression of multiple genes simultaneously or sequentially. Crosstalk between hormones and epigenetics is pivotal to dynamically coordinate this process. Hormonal signals can guide the addition and removal of epigenetic marks, steering gene expression. Conversely, DNA methylation, histone modifications and non-coding RNAs can modulate regional chromatin structure and accessibility and regulate the expression of numerous (hormone-related) genes. Here, we provide a review of the interplay between the classical insect hormones, ecdysteroids and juvenile hormones, and epigenetics. We summarize the mode-of-action and roles of these hormones in post-embryonic development, and provide a general overview of epigenetic mechanisms. We then highlight recent advances on the interactions between these hormonal pathways and epigenetics, and their involvement in development. Furthermore, we give an overview of several 'omics techniques employed in the field. Finally, we discuss which questions remain unanswered and possible avenues for future research.

Transcriptome-wide analysis uncovers regulatory elements of the antennal transcriptome repertoire of bumblebee at different life stages

Bumblebees are crucial pollinators, providing essential ecosystem services and global food production. The success of pollination services relies on the interaction between sensory organs and the environment. The antenna functions as a versatile multi-sensory organ, pivotal in mediating chemosensory/olfactory information, and governs adaptive responses to environmental changes. Despite an increasing number of RNA-sequencing studies on insect antenna, comprehensive antennal transcriptome studies at the different life stages were not elucidated systematically. Here, we quantified the expression profile and dynamics of coding/microRNA genes of larval head and antennal tissues from early- and late-stage pupa to the adult of Bombus terrestris as suitable model organism among pollinators. We further performed Pearson correlation analyses on the gene expression profiles of the antennal transcriptome from larval head tissue to adult stages, exploring both positive and negative expression trends. The positively correlated coding genes were primarily enriched in sensory perception of chemical stimuli, ion transport, transmembrane transport processes and olfactory receptor activity. Negatively correlated genes were mainly enriched in organic substance biosynthesis and regulatory mechanisms underlying larval body patterning and the formation of juvenile antennal structures. As post-transcriptional regulators, miR-1000-5p, miR-13b-3p, miR-263-5p and miR-252-5p showed positive correlations, whereas miR-315-5p, miR-92b-3p, miR-137-3p, miR-11-3p and miR-10-3p exhibited negative correlations in antennal tissue. Notably, based on the inverse expression relationship, positively and negatively correlated microRNA (miRNA)-mRNA target pairs revealed that differentially expressed miRNAs predictively targeted genes involved in antennal development, shaping antennal structures and regulating antenna-specific functions. Our data serve as a foundation for understanding stage-specific antennal transcriptomes and large-scale comparative analysis of transcriptomes in different insects.

Epigenetic studies in insects and the valproic acid perspective

Abstract Valproic acid in association with sodium valproate (VPA) is an important anticonvulsant drug used for decades to treat neurological disorders. VPA also acts as an epigenetic modulator by inhibiting histone deacetylases, permitting histone acetylation, affecting the DNA and histone methylation status and gene expression, and inducing chromatin remodeling. Insects represent an important animal model for studies in several areas of science. Their high phenotypic plasticity makes them alternative models for epigenetic studies. This brief review emphasizes recent reports on insect epigenetics and the contribution of studies on the VPA action in insects, including effects on epigenetic markers, extending the pharmacological understanding of the potential of this drug, and demonstrating the usefulness of insects as an alternative animal model to drug studies.

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

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.

Interaction Analysis of lncRNA and mRNA Based on the Full-Length Transcriptome of the Nymph-to-Adult Developmental Transition of Sogatella furcifera

Little is known on how long noncoding RNAs (lncRNAs) and mRNAs cooperatively participate in regulating the nymph-to-adult development transition of Sogatella furcifera. Herein, lncRNA and mRNA libraries were constructed in three different developmental stages of S. furcifera, namely, prior to (PE), during (DE), and after (AE) ecdysis. Overall, 4649 lncRNAs were identified and divided into intergenic (53.90%), intronic (1.33%), sense (8.99%), antisense (21.75%), and bidirectional (3.94%) lncRNAs. Moreover, 795 differentially expressed lncRNAs were identified. Specifically, upon comparing PE and DE, 2719 target mRNAs were predicted for 574 lncRNAs. Upon comparing PE and AE, 2816 target mRNAs were predicted for 627 lncRNAs. Finally, upon comparing DE and AE, 51 target mRNAs were predicted for 35 lncRNAs. Kyoto Encyclopedia of Genes and Genome functional enrichment analysis revealed that the target genes of 795 lncRNAs were enriched in metabolic pathways, amino sugar and nucleotide sugar metabolism, and fatty acid metabolism. Subsequently, interaction analysis indicated that MSTRG.16086.1, MSTRG.16087.1, and MSTRG.2447.1 were functionally associated with cuticle protein and chitin biosynthesis. Finally, 11 differentially expressed lncRNAs were significantly enriched in 3rd and 4th instar nymphs. Our findings suggest that lncRNAs play a critical regulatory role during the molting of S. furcifera.

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.

Comprehensive investigation and regulatory function of lncRNAs engaged in western honey bee larval immune response to Ascosphaera apis invasion

Ascosphaera apis is a fungal pathogen that exclusively infects bee larvae, causing chalkbrood disease, which results in severe damage for beekeeping industry. Long non-coding RNAs (lncRNAs) are versatile regulators in various biological processes such as immune defense and host-pathogen interaction. However, expression pattern and regulatory role of lncRNAs involved in immune response of bee host to A. apis invasion is still very limited. Here, the gut tissues of Apis mellifera ligustica 4-, 5-, and 6-day-old larvae inoculated by A. apis spores (AmT1, AmT2, and AmT3 groups) and corresponding un-inoculated larval guts (AmCK1, AmCK2, and AmCK3 groups) were prepared and subjected to deep sequencing, followed by identification of lncRNAs, analysis of differentially expressed lncRNAs (DElncRNAs), and investigation of competing endogenous RNA (ceRNA) network. In total, 3,746 A. m. ligustica lncRNAs were identified, including 78 sense lncRNAs, 891 antisense lncRNAs, 1,893 intergenic lncRNAs, 346 bidirectional lncRNAs, and 210 intronic lncRNAs. In the 4-, 5-, and 6- comparison groups, 357, 236, and 505 DElncRNAs were discovered. Additionally, 217, 129, and 272 DElncRNAs were respectively predicted to regulate neighboring genes via cis-acting manner, and these targets were associated with a series of GO terms and KEGG pathways of great importance, such as response to stimulus and Jak-STAT signaling pathway. Moreover, 197, 95, and 356 DElncRNAs were observed to target 10, eight, and 21 DEmiRNAs and further target 147, 79, and 315 DEmRNAs, forming complex regulatory networks. Further investigation suggested that these targets were engaged in several key cellular and humoral immune pathways, such as phagosome and MAPK signaling pathway. Ultimately, the expression trends of nine randomly selected DElncRNAs were verified by RT-qPCR, confirming the authenticity and reliability of our transcriptome data. Findings in this current work not only provide candidate DElncRNAs for functional study, but also lay a foundation for unclosing the mechanism underlying DElncRNA-regulated larval immune responses to A. apis invasion.

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.

Long noncoding RNAs: Emerging players regulating innate immune memory in the red flour beetle

A variety of strategies have been evolved to eradicate invading microbes. Phagocytes have developed in vertebrates and invertebrates to confer a non-specific immune response to pathogens. Besides, vertebrates have evolved lymphocytes to develop memory cells that can quickly respond upon the next exposure to the same pathogen. Although lymphocytes are absent in invertebrates, historical evidence, dating back to the 1920s, indicated the presence of immune memory in invertebrates. However, the concept of long-lasting non-specific defense predominated until recent evidence has been introduced in the first decade of the 21st century. Although more evidence has been introduced later, the molecular mechanism underlying the innate immune memory is largely undefined in invertebrates. Long noncoding RNAs (lncRNAs) have demonstrated a role in regulating various biological processes, including immune response. In this review, we will explore the potential role of lncRNAs in developing innate immune memory in the red flour beetle (Tribolium castaneum).

Characterization of Long Non-Coding RNAs in the Bollworm, Helicoverpa zea, and Their Possible Role in Cry1Ac-Resistance

Multiple insect pest species have developed field resistance to Bt-transgenic crops. There has been a significant amount of research on protein-coding genes that contribute to resistance, such as the up-regulation of protease activity or altered receptors. However, our understanding of the role of non-protein-coding mechanisms in Bt-resistance is minimal, as is also the case for resistance to chemical pesticides. To address this problem relative to Bt, RNA-seq was used to examine statistically significant, differential gene expression between a Cry1Ac-resistant (~100-fold resistant) and Cry1Ac-susceptible strain of Helicoverpa zea, a prevalent caterpillar pest in the USA. Significant differential expression of putative long non-coding RNAs (lncRNAs) was found in the Cry1Ac-resistant strain (58 up- and 24 down-regulated gene transcripts with an additional 10 found only in resistant and four only in susceptible caterpillars). These lncRNAs were examined as potential pseudogenes and for their genomic proximity to coding genes, both of which can be indicative of regulatory relationships between a lncRNA and coding gene expression. A possible pseudogenic lncRNA was found with similarities to a cadherin. In addition, putative lncRNAs were found significantly proximal to a serine protease, ABC transporter, and CYP coding genes, potentially involved in the mechanism of Bt and/or chemical insecticide resistance. Characterization of non-coding genetic mechanisms in Helicoverpa zea will improve the understanding of the genomic evolution of insect resistance, improve the identification of specific regulators of coding genes in general (some of which could be important in resistance), and is the first step for potentially targeting these regulators for pest control and resistance management (using molecular approaches, such as RNAi and others).

lnc94638 is a testis-specific long non-coding RNA involved in spermatozoa formation in Zeugodacus cucurbitae (Coquillett)

Long non-coding RNAs (lncRNAs) generally display tissue-specific distributions, and testis-specific lncRNAs form the highest proportion of lncRNAs in many species. Here, we presented a detailed analysis of testis-specific lncRNAs in the melon fly, Zeugodacus cucurbitae, a highly destructive insect pest of cucurbitaceous and other related crops. Most testis-specific lncRNAs were found to be long intergenic non-coding RNAs (lincRNA). The size distribution of these lncRNAs ranged between 600 and 1000 nucleotides. Testis-specific lncRNAs that harboured one isoform number and two exons were the most abundant. Compared to other male tissues, the testis had more highly expressed lncRNAs. The quantitative real-time polymerase chain reaction results of 10 randomly selected testis-specific lncRNAs showed expression patterns consistent with RNA-seq data. Further analysis of the most highly expressed testis-specific lncRNA, lnc94638, was undertaken. Fluorescent in situ hybridization assays localized lnc94638 to the apical region of the testis that contains mature spermatozoa. RNA interference-mediated knockdown of lnc94638 expression reduced spermatozoa numbers and impaired the fertility of Z. cucurbitae male. This study provides a catalogue of testis-specific lncRNAs, shows that the testis-specific lnc94638 is involved in spermatogenesis and has the potential to be used for treating male sterility.

Contribution of Epigenetic Mechanisms in the Regulation of Environmentally-Induced Polyphenism in Insects

Simple Summary

Polyphenism is a widespread phenomenon in insects that allows organisms to produce alternative and discrete phenotypes in response to environmental conditions. Epigenetic mechanisms, including histone post-translational modifications, DNA methylation and non-coding RNAs, are essential mechanisms that can promote rapid and flexible changes in the expression of transcriptional programs associated with the production of alternative phenotypes. This review summarizes knowledge regarding the contribution of those mechanisms in the regulation of the most-studied examples of polyphenism in insects.

Abstract

Many insect species display a remarkable ability to produce discrete phenotypes in response to changes in environmental conditions. Such phenotypic plasticity is referred to as polyphenism. Seasonal, dispersal and caste polyphenisms correspond to the most-studied examples that are environmentally-induced in insects. Cues that induce such dramatic phenotypic changes are very diverse, ranging from seasonal cues, habitat quality changes or differential larval nutrition. Once these signals are perceived, they are transduced by the neuroendocrine system towards their target tissues where gene expression reprogramming underlying phenotypic changes occur. Epigenetic mechanisms are key regulators that allow for genome expression plasticity associated with such developmental switches. These mechanisms include DNA methylation, chromatin remodelling and histone post-transcriptional modifications (PTMs) as well as non-coding RNAs and have been studied to various extents in insect polyphenism. Differential patterns of DNA methylation between phenotypes are usually correlated with changes in gene expression and alternative splicing events, especially in the cases of dispersal and caste polyphenism. Combinatorial patterns of histone PTMs provide phenotype-specific epigenomic landscape associated with the expression of specific transcriptional programs, as revealed during caste determination in honeybees and ants. Alternative phenotypes are also usually associated with specific non-coding RNA profiles. This review will provide a summary of the current knowledge of the epigenetic changes associated with polyphenism in insects and highlights the potential for these mechanisms to be key regulators of developmental transitions triggered by environmental cues.

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.

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.

Caste-specific gene expression underlying the differential adult brain development in the honeybee Apis mellifera

Apis mellifera adult workers feature more developed key brain regions than queens, which allows them to cope with the broad range of duties they need to perform in a colony. However, at the end of larval development, the brain of queens is largely more developed than that of workers. Major morphogenetic changes take place after metamorphosis that shift caste-specific brain development. Here, we tested the hypothesis that this phenomenon is hormonally governed and involves differential gene expression. Our molecular screening approach revealed a set of differentially expressed genes in Pp (first pharate-adult phase) brains between castes mainly coding for tissue remodelling and energy-converting proteins (e.g. hex 70a and ATPsynβ). An in-depth qPCR analysis of the transcriptional behaviour during pupal and pharate-adult developmental stage in both castes and in response to artificially augmented hormone titres of 18 genes/variants revealed that: i. subtle differences in hormone titres between castes might be responsible for the differential expression of the EcR and insulin/insulin-like signalling (IIS) pathway genes; ii. the morphogenetic activity of the IIS in brain development must be mediated by ILP-2, iii. which together with the tum, mnb and caspase system, can constitute the molecular effectors of the caste-specific opposing brain developmental trajectories.

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.

Genome-wide characterization and developmental expression profiling of long non-coding RNAs in Sogatella furcifera

The genome-wide characterization of long non-coding RNA (lncRNA) in insects demonstrates their importance in fundamental biological processes. Essentially, an in-depth understanding of the functional repertoire of lncRNA in insects is pivotal to insect resources utilization and sustainable pest control. Using a custom bioinformatics pipeline, we identified 1861 lncRNAs encoded by 1852 loci in the Sogatella furcifera genome. We profiled lncRNA expression in different developmental stages and observed that the expression of lncRNAs is more highly temporally restricted compared to protein-coding genes. More up-regulated Sogatella furcifera lncRNA expressed in the embryo, 4th and 5th instars, suggesting that increased lncRNA levels may play a role in these developmental stages. We compared the relationship between the expression of Sogatella furcifera lncRNA and its nearest protein gene and found that lncRNAs were more correlated to their downstream coding neighbors on the opposite strand. Our genome-wide profiling of lncRNAs in Sogatella furcifera identifies exciting candidates for characterization of lncRNAs, and also provides information on lncRNA regulation during insect development.

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.

Genome-wide analysis of long non-coding RNAs in adult tissues of the melon fly, Zeugodacus cucurbitae (Coquillett)

Background

Long non-coding RNAs (lncRNAs) are involved in many fundamental biological processes, such as transcription regulation, protein degradation, and cell differentiation. Information on lncRNA in the melon fly, Zeugodacus cucurbitae (Coquillett) is currently limited.

Results

We constructed 24 RNA-seq libraries from eight tissues (midgut, Malpighian tubules, fat body, ovary, and testis) of Z. cucurbitae adults. A total of 3124 lncRNA transcripts were identified. Among those, 1464 were lincRNAs, 1037 were intronic lncRNAs, 301 were anti-sense lncRNAs, and 322 were sense lncRNAs. The majority of lncRNAs contained two exons and one isoform. Differentially expressed lncRNAs were analyzed between tissues, and Malpighian tubules versus testis had the largest number. Some lncRNAs exhibited strong tissue specificity. Specifically expressed lncRNAs were identified and filtered in tissues of female and male Z. cucurbitae based on their expression levels. Four midgut-specific lncRNAs were validated by quantitative real-time polymerase chain reaction (RT-qPCR), and the data were consistent with RNA-seq data. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of targets of midgut-specific lncRNAs indicated an enrichment of the metabolic process.

Conclusions

This was the first systematic identification of lncRNA in the melon fly. Expressions of lncRNAs in multiple adult tissues were evaluated by quantitative transcriptomic analysis. These qualitative and quantitative analyses of lncRNAs, especially the tissue-specific lncRNAs in Z. cucurbitae, provide useful data for further functional studies.

Long Non-coding RNA Derived from lncRNA-mRNA Co-expression Networks Modulates the Locust Phase Change

Long non-coding RNAs (lncRNAs) regulate various biological processes ranging from gene expression to animal behavior. Although protein-coding genes, microRNAs, and neuropeptides play important roles in the regulation of phenotypic plasticity in migratory locust, empirical studies on the function of lncRNAs in this process remain limited. Here, we applied high-throughput RNA-seq to compare the expression patterns of lncRNAs and mRNAs in the time course of locust phase change. We found that lncRNAs responded more rapidly at the early stages of phase transition. Functional annotations demonstrated that early changed lncRNAs employed different pathways in isolation and crowding phases to cope with changes in the population density. Two overlapping hub lncRNA loci in the crowding and isolation networks were screened for functional verification. One of them, LNC1010057, was validated as a potential regulator of locust phase change. This work offers insights into the molecular mechanism underlying locust phase change and expands the scope of lncRNA functions in animal behavior.

Comprehensive Identification and Characterization of Long Non-coding RNAs Associated With Rice Black-Streaked Dwarf Virus Infection in Laodelphax striatellus (Fallén) Midgut

Long non-coding RNAs (lncRNAs) are involved in a variety of biological functions through transcriptional and post-transcriptional regulation. However, little is known about their functions in the process of insect mediated virus transmission. In the present study, we analyzed using RNA-Seq, the lncRNAs that were differentially expressed in response to Rice black-streaked dwarf virus (RBSDV) infection in Laodelphax striatellus (Fallén) midgut. A total of 13,927 lncRNAs were identified and over 69% were assigned to intergenic regions. Among them, 176 lncRNAs were differentially expressed and predicted to target 168 trans-regulatory genes. Ten differentially expressed lncRNAs were selected and their expression changes were validated by RT-qPCR. KEGG analysis showed that these target genes were enriched in the essential biological process, such as Purine metabolism, Valine, leucine and isoleucine degradation, and Fatty acid elongation. The expression levels of the differentially expressed lncRNAs and the predicted target genes that were significantly enriched in the Human papillomavirus infection pathway were analyzed by RT-qPCR. The results showed that several lncRNAs were co-expressed with their target genes. One of the lncRNAs called MSTRG15394 and its target gene, encoding a secreted protease inhibitor (PI), were up-regulated at the transcriptional level after RBSDV infection. Knockdown of MSTRG15394 could down-regulate the PI expression at mRNA level. Inhibition of either MSTRG15394 or PI expression by RNA interference promoted RBSDV accumulation in L. striatellus midgut. Our finding provides new insights into the function of lncRNAs in regulating virus infection in an important insect vector.

The nuclear and mitochondrial genomes of Frieseomelitta varia - a highly eusocial stingless bee (Meliponini) with a permanently sterile worker caste

BACKGROUND

Most of our understanding on the social behavior and genomics of bees and other social insects is centered on the Western honey bee, Apis mellifera. The genus Apis, however, is a highly derived branch comprising less than a dozen species, four of which genomically characterized. In contrast, for the equally highly eusocial, yet taxonomically and biologically more diverse Meliponini, a full genome sequence was so far available for a single Melipona species only. We present here the genome sequence of Frieseomelitta varia, a stingless bee that has, as a peculiarity, a completely sterile worker caste.

RESULTS

The assembly of 243,974,526 high quality Illumina reads resulted in a predicted assembled genome size of 275 Mb composed of 2173 scaffolds. A BUSCO analysis for the 10,526 predicted genes showed that these represent 96.6% of the expected hymenopteran orthologs. We also predicted 169,371 repetitive genomic components, 2083 putative transposable elements, and 1946 genes for non-coding RNAs, largely long non-coding RNAs. The mitochondrial genome comprises 15,144 bp, encoding 13 proteins, 22 tRNAs and 2 rRNAs. We observed considerable rearrangement in the mitochondrial gene order compared to other bees. For an in-depth analysis of genes related to social biology, we manually checked the annotations for 533 automatically predicted gene models, including 127 genes related to reproductive processes, 104 to development, and 174 immunity-related genes. We also performed specific searches for genes containing transcription factor domains and genes related to neurogenesis and chemosensory communication.

CONCLUSIONS

The total genome size for F. varia is similar to the sequenced genomes of other bees. Using specific prediction methods, we identified a large number of repetitive genome components and long non-coding RNAs, which could provide the molecular basis for gene regulatory plasticity, including worker reproduction. The remarkable reshuffling in gene order in the mitochondrial genome suggests that stingless bees may be a hotspot for mtDNA evolution. Hence, while being just the second stingless bee genome sequenced, we expect that subsequent targeting of a selected set of species from this diverse clade of highly eusocial bees will reveal relevant evolutionary signals and trends related to eusociality in these important pollinators.

Effects of larval Age at Grafting and Juvenile Hormone on Morphometry and Reproductive Quality Parameters of in Vitro Reared Honey Bees (Hymenoptera: Apidae)

The honey bee queen plays a central role in the Apis mellifera L. (Hymenoptera: Apidae) colony, and her high reproductive capacity is fundamental for building up the workforce of a colony. Caste development in honey bee females involves elaborate physiological pathways unleashed at the beginning of the first larval instars, with juvenile hormone (JH) playing a crucial role. Here we took advantage of established in vitro rearing techniques to conduct a 2 × 2 experimental design and test initial rearing age (young vs old) and JH treatment (JH III vs solvent control) to enlighten the role of nutrient quality and JH in shaping honey bee female fertility, morphological features related to queenliness, and key physiological parameters (hemolymph vitellogenin/Vg, sugar levels, and Vg transcript levels). Our results show that while the age at initial larval rearing had major impacts on external morphology development, where younger larvae exhibited a higher probability to develop into queen-like adults morphotypes, the JH application during the larval stage improved physiological pathways related to ovary development and metabolism during the ontogenic development. We detected that the supplementation of queen larvae with JH promoted important benefits regarding queen fertility as the increase of ovariole number and vg levels at hemolymph, both crucial factors at eggs production. The data presented here provide guidance in efforts to improve honey bee queen quality, especially in light of frequent episodes of queen failures in the beekeeping industry.

Genome-Wide Identification of Long Non-Coding RNAs and Their Regulatory Networks Involved in Apis mellifera ligustica Response to Nosema ceranae Infection

Long non-coding RNAs (lncRNAs) are a diverse class of transcripts that structurally resemble mRNAs but do not encode proteins, and lncRNAs have been proven to play pivotal roles in a wide range of biological processes in animals and plants. However, knowledge of expression patterns and potential roles of honeybee lncRNA response to Nosema ceranae infection is completely unknown. Here, we performed whole transcriptome strand-specific RNA sequencing of normal midguts of Apis mellifera ligustica workers (Am7CK, Am10CK) and N. ceranae-inoculated midguts (Am7T, Am10T), followed by comprehensive analyses using bioinformatic and molecular approaches. A total of 6353 A. m. ligustica lncRNAs were identified, including 4749 conserved lncRNAs and 1604 novel lncRNAs. These lncRNAs had minimal sequence similarities with other known lncRNAs in other species; however, their structural features were similar to counterparts in mammals and plants, including shorter exon and intron length, lower exon number, and lower expression level, compared with protein-coding transcripts. Further, 111 and 146 N. ceranae-responsive lncRNAs were identified from midguts at 7-days post-inoculation (dpi) and 10 dpi compared with control midguts. Twelve differentially expressed lncRNAs (DElncRNAs) were shared by Am7CK vs. Am7T and Am10CK vs. Am10T comparison groups, while the numbers of unique DElncRNAs were 99 and 134, respectively. Functional annotation and pathway analysis showed that the DElncRNAs may regulate the expression of neighboring genes by acting in cis and trans fashion. Moreover, we discovered 27 lncRNAs harboring eight known miRNA precursors and 513 lncRNAs harboring 2257 novel miRNA precursors. Additionally, hundreds of DElncRNAs and their target miRNAs were found to form complex competitive endogenous RNA (ceRNA) networks, suggesting that these DElncRNAs may act as miRNA sponges. Furthermore, DElncRNA-miRNA-mRNA networks were constructed and investigated, the results demonstrated that a portion of the DElncRNAs were likely to participate in regulating the host material and energy metabolism as well as cellular and humoral immune host responses to N. ceranae invasion. Our findings revealed here offer not only a rich genetic resource for further investigation of the functional roles of lncRNAs involved in the A. m. ligustica response to N. ceranae infection, but also a novel insight into understanding the host-pathogen interaction during honeybee microsporidiosis.

Filtered reproductive long non-coding RNAs by genome-wide analyses of goat ovary at different estrus periods

Background

The goat is an important farm animal. Reproduction is an important process of goat farming. The ovary is the most important reproductive organ for goats. In recent years, an increasing number of long non-coding RNAs (lncRNAs) have been implicated in the regulation of mammal reproduction. However, there are few studies on the function of lncRNAs in reproduction, particularly lncRNAs in the ovary.

Results

The sequencing of goat ovaries generated 1,122,014,112 clean reads, and 4926 lncRNAs and 1454 TUCPs (transcripts of uncertain coding potential) were identified for further analysis by using the coding potential analysis software, CNCI, CPC and Pfam-sca. There were 115 /22 differential lncRNAs /TUCPs transcripts between the ovaries of the luteal phase and the follicular phase. We predicted the related genes of lncRNA /TUCP based on co-expression and co-localization methods. In total, 2584 /904 genes were predicted by co-expression, and 326/73 genes were predicted by co-localization. The functions of these genes were further analyzed with GO and KEGG analysis. The results showed that lncRNAs /TUCPs, which are highly expressed in goat ovaries in the luteal phase, are mainly associated with the synthesis of progesterone, and we filtered the lncRNAs /TUCPs, such as XR_001918177.1 and TUCP_001362, which may regulate the synthesis of progesterone; lncRNAs /TUCPs, which are highly expressed in goat ovaries in the follicular phase, are mainly associated with oogenesis and the maturation of oocytes, and we filtered the lncRNAs /TUCPs that may regulate the oogenesis and maturation of oocyte, such as XR_001917388.1 and TUCP_000849.

Conclusion

The present study provided the genome expression profile of lncRNAs /TUCPs in goat ovaries at different estrus periods and filtered the potential lncRNAs /TUCPs associated with goat reproduction. These results are helpful to further study the molecular mechanisms of goat reproduction.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5268-7) contains supplementary material, which is available to authorized users.

Autosomal and Mitochondrial Adaptation Following Admixture: A Case Study on the Honeybees of Reunion Island

The honeybee population of the tropical Reunion Island is a genetic admixture of the Apis mellifera unicolor subspecies, originally described in Madagascar, and of European subspecies, mainly A. m. carnica and A. m. ligustica, regularly imported to the island since the late 19th century. We took advantage of this population to study genetic admixing of the tropical-adapted indigenous and temperate-adapted European genetic backgrounds. Whole genome sequencing of 30 workers and 6 males from Reunion, compared with samples from Europe, Madagascar, Mauritius, Rodrigues, and the Seychelles, revealed the Reunion honeybee population to be composed on an average of 53.2 ± 5.9% A. m. unicolor nuclear genomic background, the rest being mainly composed of A. m. carnica and to a lesser extent A. m. ligustica. In striking contrast to this, only 1 out of the 36 honeybees from Reunion had a mitochondrial genome of European origin, suggesting selection has favored the A. m. unicolor mitotype, which is possibly better adapted to the island’s bioclimate. Local ancestry was determined along the chromosomes for all Reunion samples, and a test for preferential selection for the A. m. unicolor or European background revealed 15 regions significantly associated with the A. m. unicolor lineage and 9 regions with the European lineage. Our results provide insights into the long-term consequences of introducing exotic specimen on the nuclear and mitochondrial genomes of locally adapted populations.

The transcriptomic changes associated with the development of social parasitism in the honeybee Apis mellifera capensis

Social insects are characterized by the division of labor. Queens usually dominate reproduction, whereas workers fulfill non-reproductive age-dependent tasks to maintain the colony. Although workers are typically sterile, they can activate their ovaries to produce their own offspring. In the extreme, worker reproduction can turn into social parasitism as in Apis mellifera capensis. These intraspecific parasites occupy a host colony, kill the resident queen, and take over the reproductive monopoly. Because they exhibit a queenlike behavior and are also treated like queens by the fellow workers, they are so-called pseudoqueens. Here, we compare the development of parasitic pseudoqueens and social workers at different time points using fat body transcriptome data. Two complementary analysis methods-a principal component analysis and a time course analysis-led to the identification of a core set of genes involved in the transition from a social worker into a highly fecund parasitic pseudoqueen. Comparing our results on pseudoqueens with gene expression data of honeybee queens revealed many similarities. In addition, there was a set of specific transcriptomic changes in the parasitic pseudoqueens that differed from both, queens and social workers, which may be typical for the development of the social parasitism in A. m. capensis.

Genome-wide identification of long non-coding RNA genes and their association with insecticide resistance and metamorphosis in diamondback moth, Plutella xylostella

Long non-coding RNA (lncRNA) is a class of noncoding RNA >200 bp in length that has essential roles in regulating a variety of biological processes. Here, we constructed a computational pipeline to identify lncRNA genes in the diamondback moth (Plutella xylostella), a major insect pest of cruciferous vegetables. In total, 3,324 lncRNAs corresponding to 2,475 loci were identified from 13 RNA-Seq datasets, including samples from parasitized, insecticide-resistant strains and different developmental stages. The identified P. xylostella lncRNAs had shorter transcripts and fewer exons than protein-coding genes. Seven out of nine randomly selected lncRNAs were validated by strand-specific RT-PCR. In total, 54-172 lncRNAs were specifically expressed in the insecticide resistant strains, among which one lncRNA was located adjacent to the sodium channel gene. In addition, 63-135 lncRNAs were specifically expressed in different developmental stages, among which three lncRNAs overlapped or were located adjacent to the metamorphosis-associated genes. These lncRNAs were either strongly or weakly co-expressed with their overlapping or neighboring mRNA genes. In summary, we identified thousands of lncRNAs and presented evidence that lncRNAs might have key roles in conferring insecticide resistance and regulating the metamorphosis development in P. xylostella.

Genome-wide identification of lncRNAs associated with chlorantraniliprole resistance in diamondback moth Plutella xylostella (L.)

BACKGROUND

Long noncoding RNAs (lncRNAs) are now considered important regulatory factors, with a variety of biological functions in many species including insects. Some lncRNAs have the ability to show rapid responses to diverse stimuli or stress factors and are involved in responses to insecticide. However, there are no reports to date on the characterization of lncRNAs associated with chlorantraniliprole resistance in Plutella xylostella.

RESULTS

Nine RNA libraries constructed from one susceptible (CHS) and two chlorantraniliprole-resistant P. xylostella strains (CHR, ZZ) were sequenced, and 1309 lncRNAs were identified, including 877 intergenic lncRNAs, 190 intronic lncRNAs, 76 anti-sense lncRNAs and 166 sense-overlapping lncRNAs. Of the identified lncRNAs, 1059 were novel. Furthermore, we found that 64 lncRNAs were differentially expressed between CHR and CHS and 83 were differentially expressed between ZZ and CHS, of which 22 were differentially expressed in both CHR and ZZ. Most of the differentially expressed lncRNAs were hypothesized to be associated with chlorantraniliprole resistance in P. xylostella. The targets of lncRNAs via cis- (<10 kb upstream and downstream) or trans- (Pearson's correlation, r > 0.9 or < -0.9, P < 0.05) regulatory effects were also identified; many of the differently expressed lncRNAs were correlated with various important protein-coding genes involved in insecticide resistance, such as the ryanodine receptor, uridine diphosphate glucuronosyltransferase (UGTs), cytochrome P450, esterase and the ATP-binding cassette transporter.

CONCLUSIONS

Our results represent the first global identification of lncRNAs associated with chlorantraniliprole resistance in P. xylostella. These results will facilitate future studies of the regulatory mechanisms of lncRNAs in chlorantraniliprole and other insecticide resistance and in other biological processes in P. xylostella.

Endogenous RNAi Pathways Are Required in Neurons for Dauer Formation in Caenorhabditis elegans

Animals can adapt to unfavorable environments through changes in physiology or behavior. In the nematode, Caenorhabditis elegans, environmental conditions perceived early in development determine whether the animal enters either the reproductive cycle, or enters into an alternative diapause stage named dauer. Here, we show that endogenous RNAi pathways play a role in dauer formation in crowding (high pheromone), starvation, and high temperature conditions. Disruption of the Mutator proteins or the nuclear Argonaute CSR-1 result in differential dauer-deficient phenotypes that are dependent upon the experienced environmental stress. We provide evidence that the RNAi pathways function in chemosensory neurons for dauer formation, upstream of the TGF-β and insulin signaling pathways. In addition, we show that Mutator MUT-16 expression in a subset of individual pheromone-sensing neurons is sufficient for dauer formation in high pheromone conditions, but not in starvation or high temperature conditions. Furthermore, we also show that MUT-16 and CSR-1 are required for expression of a subset of G proteins with functions in the detection of pheromone components. Together, our data suggest a model where Mutator-amplified siRNAs that associate with the CSR-1 pathway promote expression of genes required for the detection and signaling of environmental conditions to regulate development and behavior in C. elegans This study highlights a mechanism whereby RNAi pathways mediate the link between environmental stress and adaptive phenotypic plasticity in animals.

Long Non-Coding RNAs Regulating Immunity in Insects

Recent advances in modern technology have led to the understanding that not all genetic information is coded into protein and that the genomes of each and every organism including insects produce non-coding RNAs that can control different biological processes. Among RNAs identified in the last decade, long non-coding RNAs (lncRNAs) represent a repertoire of a hidden layer of internal signals that can regulate gene expression in physiological, pathological, and immunological processes. Evidence shows the importance of lncRNAs in the regulation of host–pathogen interactions. In this review, an attempt has been made to view the role of lncRNAs regulating immune responses in insects.

Differential gene expression underlying ovarian phenotype determination in honey bee, Apis mellifera L., caste development

Adult honey bee queens and workers drastically differ in ovary size. This adult ovary phenotype difference becomes established during the final larval instar, when massive programmed cell death leads to the degeneration of 95-99% of the ovariole anlagen in workers. The higher juvenile hormone (JH) levels in queen larvae protect the ovaries against such degeneration. To gain insights into the molecular architecture underlying this divergence critical for adult caste fate and worker sterility, we performed a microarray analysis on fourth and early fifth instar queen and worker ovaries. For the fourth instar we found nine differentially expressed genes (DEGs) with log₂FC > 1.0, but this number increased to 56 in early fifth-instar ovaries. We selected 15 DEGs for quantitative PCR (RT-qPCR) analysis. Nine differed significantly by the variables caste and/or development. Interestingly, genes with enzyme functions were higher expressed in workers, while those related to transcription and signaling had higher transcript levels in queens. For the RT-qPCR confirmed genes we analyzed their response to JH. This revealed a significant up-regulation for two genes, a short chain dehydrogenase reductase (sdr) and a heat shock protein 90 (hsp90). Five other genes, including hsp60 and hexamerin 70b (hex70b), were significantly down-regulated by JH. The sdr gene had previously come up as differentially expressed in other transcriptome analyses on honey bee larvae and heat shock proteins are frequently involved in insect hormone responses, this making them interesting candidates for further functional assays.

Chapter 6. Comparative genomic approaches to investigate molecular traits specific to social insects

Ageing is a feature of nearly all known organisms and, by its connection to survival, appears to trade off with fecundity. However, in some organisms such as in queens of social insects, this negative relation appears reversed and individuals live long and reproduce much. Since new experimental techniques, transcriptomes and genomes of many social insects have recently become available, a comparison of these data in a phylogenetic framework becomes feasible. This allows the study of general trends, species specific oddities and evolutionary dynamics of the molecular properties and changes which underlie ageing, fecundity and the reversal of this negative association. In the framework of social insect evolution, we review the most important recent insights, computational methods, their applications and data resources which are available.

Structure Prediction: New Insights into Decrypting Long Noncoding RNAs

Long noncoding RNAs (lncRNAs), which form a diverse class of RNAs, remain the least understood type of noncoding RNAs in terms of their nature and identification. Emerging evidence has revealed that a small number of newly discovered lncRNAs perform important and complex biological functions such as dosage compensation, chromatin regulation, genomic imprinting, and nuclear organization. However, understanding the wide range of functions of lncRNAs related to various processes of cellular networks remains a great experimental challenge. Structural versatility is critical for RNAs to perform various functions and provides new insights into probing the functions of lncRNAs. In recent years, the computational method of RNA structure prediction has been developed to analyze the structure of lncRNAs. This novel methodology has provided basic but indispensable information for the rapid, large-scale and in-depth research of lncRNAs. This review focuses on mainstream RNA structure prediction methods at the secondary and tertiary levels to offer an additional approach to investigating the functions of lncRNAs.

Systematic Identification and Characterization of Long Non-Coding RNAs in the Silkworm, Bombyx mori

Long noncoding RNAs (lncRNAs) are emerging as important regulators in various biological processes. However, to date, no systematic characterization of lncRNAs has been reported in the silkworm Bombyx mori. In the present study, we generated eighteen RNA-seq datasets with relatively high depth. Using an in-house designed lncRNA identification pipeline, 11,810 lncRNAs were identified for 5,556 loci. Among these lncRNAs, 474 transcripts were intronic lncRNAs (ilncRNAs), 6,250 transcripts were intergenic lncRNAs (lincRNAs), and 5,086 were natural antisense lncRNAs (lncNATs). Compared with protein-coding mRNAs, silkworm lncRNAs are shorter in terms of full length but longer in terms of exon and intron length. In addition, lncRNAs exhibit a lower level of sequence conservation, more repeat sequences overlapped and higher tissue specificity than protein-coding mRNAs in the silkworm. We found that 69 lncRNA transcripts from 33 gene loci may function as miRNA precursors, and 104 lncRNA transcripts from 72 gene loci may act as competing endogenous RNAs (ceRNAs). In total, 49.47% of all gene loci (2,749/5,556) for which lncRNAs were identified showed sex-biased expression. Co-expression network analysis resulted in 19 modules, 12 of which revealed relatively high tissue specificity. The highlighted darkgoldenrod module was specifically associated with middle and posterior silk glands, and the hub lncRNAs within this module were co-expressed with proteins involved in translation, translocation, and secretory processes, suggesting that these hub lncRNAs may function as regulators of the biosynthesis, translocation, and secretion of silk proteins. This study presents the first comprehensive genome-wide analysis of silkworm lncRNAs and provides an invaluable resource for genetic, evolutionary, and genomic studies of B. mori.

Genome wide discovery of long intergenic non-coding RNAs in Diamondback moth (Plutella xylostella) and their expression in insecticide resistant strains

Long non-coding RNAs (lncRNAs) play important roles in genomic imprinting, cancer, differentiation and regulation of gene expression. Here, we identified 3844 long intergenic ncRNAs (lincRNA) in Plutella xylostella, which is a notorious pest of cruciferous plants that has developed field resistance to all classes of insecticides, including Bacillus thuringiensis (Bt) endotoxins. Further, we found that some of those lincRNAs may potentially serve as precursors for the production of small ncRNAs. We found 280 and 350 lincRNAs that are differentially expressed in Chlorpyrifos and Fipronil resistant larvae. A survey on P. xylostella midgut transcriptome data from Bt-resistant populations revealed 59 altered lincRNA in two resistant strains compared with the susceptible population. We validated the transcript levels of a number of putative lincRNAs in deltamethrin-resistant larvae that were exposed to deltamethrin, which indicated that this group of lincRNAs might be involved in the response to xenobiotics in this insect. To functionally characterize DBM lincRNAs, gene ontology (GO) enrichment of their associated protein-coding genes was extracted and showed over representation of protein, DNA and RNA binding GO terms. The data presented here will facilitate future studies to unravel the function of lincRNAs in insecticide resistance or the response to xenobiotics of eukaryotic cells.

Genome-wide characterization of long intergenic non-coding RNAs (lincRNAs) provides new insight into viral diseases in honey bees Apis cerana and Apis mellifera

BACKGROUND

Long non-coding RNAs (lncRNAs) are a class of RNAs that do not encode proteins. Recently, lncRNAs have gained special attention for their roles in various biological process and diseases.

RESULTS

In an attempt to identify long intergenic non-coding RNAs (lincRNAs) and their possible involvement in honey bee development and diseases, we analyzed RNA-seq datasets generated from Asian honey bee (Apis cerana) and western honey bee (Apis mellifera). We identified 2470 lincRNAs with an average length of 1011 bp from A. cerana and 1514 lincRNAs with an average length of 790 bp in A. mellifera. Comparative analysis revealed that 5 % of the total lincRNAs derived from both species are unique in each species. Our comparative digital gene expression analysis revealed a high degree of tissue-specific expression among the seven major tissues of honey bee, different from mRNA expression patterns. A total of 863 (57 %) and 464 (18 %) lincRNAs showed tissue-dependent expression in A. mellifera and A. cerana, respectively, most preferentially in ovary and fat body tissues. Importantly, we identified 11 lincRNAs that are specifically regulated upon viral infection in honey bees, and 10 of them appear to play roles during infection with various viruses.

CONCLUSIONS

This study provides the first comprehensive set of lincRNAs for honey bees and opens the door to discover lincRNAs associated with biological and hormone signaling pathways as well as various diseases of honey bee.

LncRNA-ID: Long non-coding RNA IDentification using balanced random forests

MOTIVATION

Long non-coding RNAs (lncRNAs), which are non-coding RNAs of length above 200 nucleotides, play important biological functions such as gene expression regulation. To fully reveal the functions of lncRNAs, a fundamental step is to annotate them in various species. However, as lncRNAs tend to encode one or multiple open reading frames, it is not trivial to distinguish these long non-coding transcripts from protein-coding genes in transcriptomic data.

RESULTS

In this work, we design a new tool that calculates the coding potential of a transcript using a machine learning model (random forest) based on multiple features including sequence characteristics of putative open reading frames, translation scores based on ribosomal coverage, and conservation against characterized protein families. The experimental results show that our tool competes favorably with existing coding potential computation tools in lncRNA identification.

AVAILABILITY AND IMPLEMENTATION

The scripts and data can be downloaded at https://github.com/zhangy72/LncRNA-ID.

Insights into the dynamics of hind leg development in honey bee (Apis mellifera L.) queen and worker larvae - A morphology/differential gene expression analysis

Phenotypic plasticity is a hallmark of the caste systems of social insects, expressed in their life history and morphological traits. These are best studied in bees. In their co-evolution with angiosperm plants, the females of corbiculate bees have acquired a specialized structure on their hind legs for collecting pollen. In the highly eusocial bees (Apini and Meliponini), this structure is however only present in workers and absent in queens. By means of histological sections and cell proliferation analysis we followed the developmental dynamics of the hind legs of queens and workers in the fourth and fifth larval instars. In parallel, we generated subtractive cDNA libraries for hind leg discs of queen and worker larvae by means of a Representational Difference Analysis (RDA). From the total of 135 unique sequences we selected 19 for RT-qPCR analysis, where six of these were confirmed as differing significantly in their expression between the two castes in the larval spinning stage. The development of complex structures such as the bees' hind legs, requires diverse patterning mechanisms and signaling modules, as indicated by the set of differentially expressed genes related with cell adhesion and signaling pathways.

What do studies of insect polyphenisms tell us about nutritionally-triggered epigenomic changes and their consequences?

Many insects are capable of remarkable changes in biology and form in response to their environment or diet. The most extreme example of these are polyphenisms, which are when two or more different phenotypes are produced from a single genotype in response to the environment. Polyphenisms provide a fascinating opportunity to study how the environment affects an animal’s genome, and how this produces changes in form. Here we review the current state of knowledge of the molecular basis of polyphenisms and what can be learnt from them to understand how nutrition may influence our own genomes.

DNA methylation in social insects: how epigenetics can control behavior and longevity

In eusocial insects, genetically identical individuals can exhibit striking differences in behavior and longevity. The molecular basis of such phenotypic plasticity is of great interest to the scientific community. DNA methylation, as well as other epigenetic signals, plays an important role in modulating gene expression and can therefore establish, sustain, and alter organism-level phenotypes, including behavior and life span. Unlike DNA methylation in mammals, DNA methylation in insects, including eusocial insects, is enriched in gene bodies of actively expressed genes. Recent investigations have revealed the important role of gene body methylation in regulating gene expression in response to intrinsic and environmental factors. In this review, we summarize recent advances in DNA methylation research and discuss its significance in our understanding of the epigenetic underpinnings of behavior and longevity.

Developmental regulation of ecdysone receptor (EcR) and EcR-controlled gene expression during pharate-adult development of honeybees (Apis mellifera)

Major developmental transitions in multicellular organisms are driven by steroid hormones. In insects, these, together with juvenile hormone (JH), control development, metamorphosis, reproduction and aging, and are also suggested to play an important role in caste differentiation of social insects. Here, we aimed to determine how EcR transcription and ecdysteroid titers are related during honeybee postembryonic development and what may actually be the role of EcR in caste development of this social insect. In addition, we expected that knocking-down EcR gene expression would give us information on the participation of the respective protein in regulating downstream targets of EcR. We found that in Apis mellifera females, EcR-A is the predominantly expressed variant in postembryonic development, while EcR-B transcript levels are higher in embryos, indicating an early developmental switch in EcR function. During larval and pupal stages, EcR-B expression levels are very low, while EcR-A transcripts are more variable and abundant in workers compared to queens. Strikingly, these transcript levels are opposite to the ecdysteroid titer profile. 20-hydroxyecdysone (20E) application experiments revealed that low 20E levels induce EcR expression during development, whereas high ecdysteroid titers seem to be repressive. By means of RNAi-mediated knockdown (KD) of both EcR transcript variants we detected the differential expression of 234 poly-A⁺ transcripts encoding genes such as CYPs, MRJPs and certain hormone response genes (Kr-h1 and ftz-f1). EcR-KD also promoted the differential expression of 70 miRNAs, including highly conserved ones (e.g., miR-133 and miR-375), as well honeybee-specific ones (e.g., miR-3745 and miR-3761). Our results put in evidence a broad spectrum of EcR-controlled gene expression during postembryonic development of honeybees, revealing new facets of EcR biology in this social insect.

Epigenetics as an answer to Darwin's "special difficulty"

Epigenetic modifications produce distinct phenotypes from the same genome through genome-wide transcriptional control. Recently, DNA methylation in honeybees and histone modifications in ants were found to assist the formation of caste phenotypes during development and adulthood. This insight allows us to revisit one of Darwin’s greatest challenges to his natural selection theory; the derivation of multiple forms of sterile workers within eusocial species. Differential feeding of larvae creates two distinct developmental paths between queens and workers, with workers further refined by pheromone cues. Flexible epigenetic control provides a mechanism to interpret the milieu of social cues that create distinct worker sub-caste phenotypes. Recent findings suggest a distinct use for DNA methylation before and after adult emergence. Further, a comparison of genes that are differentially methylated and transcriptionally altered upon pheromone signaling suggests that epigenetics can play a key role in mediating pheromone signals to derive sub-caste phenotypes. Epigenetic modifications may provide a molecular mechanism to Darwin’s ”special difficulty” and explain the emergence of multiple sub-phenotypes among sterile individuals.

Eusocial insects as emerging models for behavioural epigenetics

Understanding the molecular basis of how behavioural states are established, maintained and altered by environmental cues is an area of considerable and growing interest. Epigenetic processes, including methylation of DNA and post-translational modification of histones, dynamically modulate activity-dependent gene expression in neurons and can therefore have important regulatory roles in shaping behavioural responses to environmental cues. Several eusocial insect species - with their unique displays of behavioural plasticity due to age, morphology and social context - have emerged as models to investigate the genetic and epigenetic underpinnings of animal social behaviour. This Review summarizes recent studies in the epigenetics of social behaviour and offers perspectives on emerging trends and prospects for establishing genetic tools in eusocial insects.

Dimorphic ovary differentiation in honeybee (Apis mellifera) larvae involves caste-specific expression of homologs of ark and buffy cell death genes

The establishment of the number of repeated structural units, the ovarioles, in the ovaries is one of the critical events that shape caste polyphenism in social insects. In early postembryonic development, honeybee (Apis mellifera) larvae have a pair of ovaries, each one consisting of almost two hundred ovariole primordia. While practically all these ovarioles continue developing in queen-destined larvae, they undergo massive programmed cell death (PCD) in worker-destined larvae. So as to gain insight into the molecular basis of this fundamental process in caste differentiation we used quantitative PCR (qPCR) and fluorescent in situ hybridization (FISH) to investigate the expression of the Amark and Ambuffy genes in the ovaries of the two honeybee castes throughout the fifth larval instar. These are the homologs of ark and buffy Drosophila melanogaster genes, respectively, involved in activating and inhibiting PCD. Caste-specific expression patterns were found during this time-window defining ovariole number. Amark transcript levels were increased when ovariole resorption was intensified in workers, but remained at low levels in queen ovaries. The transcripts were mainly localized at the apical end of all the worker ovarioles, but appeared in only a few queen ovarioles, thus strongly suggesting a function in mediating massive ovariolar cell death in worker larvae. Ambuffy was mainly expressed in the peritoneal sheath cells covering each ovariole. The levels of Ambuffy transcripts increased earlier in the developing ovaries of queens than in workers. Consistent with a protective role against cell death, Ambuffy transcripts were localized in practically all queen ovarioles, but only in few worker ovarioles. The results are indicative of a functional relationship between the expression of evolutionary conserved cell death genes and the morphological events leading to caste-specific ovary differentiation in a social insect.