Long Non-Coding RNAs in Insects

Review on Long Non-Coding RNAs as Biomarkers and Potentially Therapeutic Targets for Bacterial Infections

The confrontation between humans and bacteria is ongoing, with strategies for combating bacterial infections continually evolving. With the advancement of RNA sequencing technology, non-coding RNAs (ncRNAs) associated with bacterial infections have garnered significant attention. Recently, long ncRNAs (lncRNAs) have been identified as regulators of sterile inflammatory responses and cellular defense against live bacterial pathogens. They are involved in regulating host antimicrobial immunity in both the nucleus and cytoplasm. Increasing evidence indicates that lncRNAs are critical for the intricate interactions between host and pathogen during bacterial infections. This paper emphatically elaborates on the potential applications of lncRNAs in clinical hallmarks, cellular damage, immunity, virulence, and drug resistance in bacterial infections in greater detail. Additionally, we discuss the challenges and limitations of studying lncRNAs in the context of bacterial infections and highlight clear directions for this promising field.

Insights into the role of non-coding RNAs in the development of insecticide resistance in insects

Pest control heavily relies on chemical pesticides has been going on for decades. However, the indiscriminate use of chemical pesticides often results in the development of resistance in pests. Almost all pests have developed some degree of resistance to pesticides. Research showed that the mechanisms of insecticide resistance in insects encompass metabolic resistance, behavioral resistance, penetration resistance and target-site resistance. Research on the these mechanisms has been mainly focused on the cis-regulatory or trans-regulatory for the insecticide resistance-related genes, with less attention paid to non-coding RNAs (ncRNAs), such as microRNA (miRNA), long non-coding RNA (lncRNA), and circular RNA (circRNA). There has been increased studies focus on understanding how these ncRNAs are involved in post-transcriptional regulation of insecticide resistance-related genes. Besides, the formatted endogenous RNA (ceRNA) regulatory networks (lncRNA/circRNA-miRNA-mRNA) has been identified as a key player in governing insect resistance formation. This review delves into the functions and underlying mechanisms of miRNA, lncRNA, and circRNA in regulating insect resistance. ncRNAs orchestrate insect resistance by modulating the expression of detoxification enzyme genes, insecticide target genes, as well as receptor genes, effectively regulating both target-site, metabolic and penetration resistance in insects. It also explores the regulatory mechanisms of ceRNA networks in the development of resistance. By enhancing our understanding of the mechanisms of ncRNAs in insecticide resistance, it will not only provide valuable insights into the new mechanisms of insecticide resistance but also help to enrich new directions in ncRNAs gene regulation research.

A Combinatorial Single-Molecule Real-Time and Illumina Sequencing Analysis of Postembryonic Gene Expression in the Asian Citrus Psyllid Diaphorina citri

Huanglongbing (HLB) is a systemic plant disease caused by 'Candidatus Liberibacter asiaticus (CLas)' and transmitted by Diaphorina citri. D. citri acquires the CLas bacteria in the nymph stage and transmits it in the adult stage, indicating that molting from the nymph to adult stages is crucial for HLB transmission. However, the available D. citri reference genomes are incomplete, and gene function studies have been limited to date. In the current research, PacBio single-molecule real-time (SMRT) and Illumina sequencing were performed to investigate the transcriptome of D. citri nymphs and adults. In total, 10,641 full-length, non-redundant transcripts (FLNRTs), 594 alternative splicing (AS) events, 4522 simple sequence repeats (SSRs), 1086 long-coding RNAs (lncRNAs), 281 transcription factors (TFs), and 4459 APA sites were identified. Furthermore, 3746 differentially expressed genes (DEGs) between nymphs and adults were identified, among which 30 DEGs involved in the Hippo signaling pathway were found. Reverse transcription-quantitative PCR (RT-qPCR) further validated the expression levels of 12 DEGs and showed a positive correlation with transcriptome data. Finally, the spatiotemporal expression pattern of genes involved in the Hippo signaling pathway exhibited high expression in the D. citri testis, ovary, and egg. Silencing of the D. citri transcriptional co-activator (DcYki) gene significantly increased D. citri mortality and decreased the cumulative molting. Our results provide useful information and a reliable data resource for gene function research of D. citri.

Role of CYP9E2 and a long non-coding RNA gene in resistance to a spinosad insecticide in the Colorado potato beetle, Leptinotarsa decemlineata

Spinosads are insecticides used to control insect pests, especially in organic farming where limited tools for pest management exist. However, resistance has developed to spinosads in economically important pests, including Colorado potato beetle (CPB), Leptinotarsa decemlineata. In this study, we used bioassays to determine spinosad sensitivity of two field populations of CPB, one from an organic farm exposed exclusively to spinosad and one from a conventional farm exposed to a variety of insecticides, and a reference insecticide naïve population. We found the field populations exhibited significant levels of resistance compared with the sensitive population. Then, we compared transcriptome profiles between the two field populations to identify genes associated primarily with spinosad resistance and found a cytochrome P450, CYP9E2, and a long non-coding RNA gene, lncRNA-2, were upregulated in the exclusively spinosad-exposed population. Knock-down of these two genes simultaneously in beetles of the spinosad-exposed population using RNA interference (RNAi) resulted in a significant increase in mortality when gene knock-down was followed by spinosad exposure, whereas single knock-downs of each gene produced smaller effects. In addition, knock-down of the lncRNA-2 gene individually resulted in significant reduction in CYP9E2 transcripts. Finally, in silico analysis using an RNA-RNA interaction tool revealed that CYP9E2 mRNA contains multiple binding sites for the lncRNA-2 transcript. Our results imply that CYP9E2 and lncRNA-2 jointly contribute to spinosad resistance in CPB, and lncRNA-2 is involved in regulation of CYP9E2 expression. These results provide evidence that metabolic resistance, driven by overexpression of CYP and lncRNA genes, contributes to spinosad resistance in CPB.

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

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

LncRNA and Protein Expression Profiles Reveal Heart Adaptation to High-Altitude Hypoxia in Tibetan Sheep

The Tibetan sheep has an intricate mechanism of adaptation to low oxygen levels, which is influenced by both genetic and environmental factors. The heart plays a crucial role in the adaptation of Tibetan sheep to hypoxia. In the present study, we utilized transcriptomic and proteomic technologies to comprehensively analyze and identify the long non-coding RNAs (lncRNAs), genes, proteins, pathways, and gene ontology (GO) terms associated with hypoxic adaptation in Tibetan sheep at three different altitudes (2500 m, 3500 m, and 4500 m). By integrating the differentially expressed (DE) lncRNA target genes, differentially expressed proteins (DEPs), and differentially expressed genes (DEGs), we were able to identify and characterize the mechanisms underlying hypoxic adaptation in Tibetan sheep. Through this integration, we identified 41 shared genes/proteins, and functional enrichment analyses revealed their close association with lipid metabolism, glycolysis/gluconeogenesis, and angiogenesis. Additionally, significant enrichment was observed in important pathways such as the PPAR signaling pathway, glycolysis/gluconeogenesis, the oxoacid metabolic process, and angiogenesis. Furthermore, the co-expression network of lncRNAs and mRNAs demonstrated that lncRNAs (MSTRG.4748.1, ENSOART00020025894, and ENSOART00020036371) may play a pivotal role in the adaptation of Tibetan sheep to the hypoxic conditions of the plateau. In conclusion, this study expands the existing database of lncRNAs and proteins in Tibetan sheep, and these findings may serve as a reference for the prevention of altitude sickness in humans.

Core genes of biomineralization and cis-regulatory long non-coding RNA regulate shell growth in bivalves

INTRODUCTION

Bivalve molluscs are abundant in marine and freshwater systems and contribute essential ecosystem services. They are characterized by an exuberant diversity of biomineralized shells and typically have two symmetric valves (a.k.a shells), but oysters (Ostreidae), some clams (Anomiidae and Chamidae) and scallops (Pectinida) have two asymmetrical valves. Predicting and modelling the likely consequences of ocean acidification on bivalve survival, biodiversity and aquaculture makes understanding shell biomineralization and its regulation a priority.

OBJECTIVES

This study aimed to a) exploit the atypical asymmetric shell growth of some bivalves and through comparative analysis of the genome and transcriptome pinpoint candidate biomineralization-related genes and regulatory long non-coding RNAs (LncRNAs) and b) demonstrate their roles in regulating shell biomineralization/growth.

METHODS

Meta-analysis of genomes, de novo generated mantle transcriptomes or transcriptomes and proteomes from public databases for six asymmetric to symmetric bivalve species was used to identify biomineralization-related genes. Bioinformatics filtering uncovered genes and regulatory modules characteristic of bivalves with asymmetric shells and identified candidate biomineralization-related genes and lncRNAs with a biased expression in asymmetric valves. A shell regrowth model in oyster and gene silencing experiments, were used to characterize candidate gene function.

RESULTS

Shell matrix genes with asymmetric expression in the mantle of the two valves were identified and unique cis-regulatory lncRNA modules characterized in Ostreidae. LncRNAs that regulate the expression of the tissue inhibitor of metalloproteinases gene family (TIMPDR) and of the shell matrix protein domain family (SMPDR) were identified. In vitro and in vivo silencing experiments revealed the candidate genes and lncRNA were associated with divergent shell growth rates and modified the microstructure of calcium carbonate (CaCO3) crystals.

CONCLUSION

LncRNAs are putative regulatory factors of the bivalve biomineralization toolbox. In the Ostreidae family of bivalves biomineralization-related genes are cis-regulated by lncRNA and modify the planar growth rate and spatial orientation of crystals in the shell.

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.

Enthralling genetic regulatory mechanisms meddling insecticide resistance development in insects: role of transcriptional and post-transcriptional events

Insecticide resistance in insects severely threatens both human health and agriculture, making insecticides less compelling and valuable, leading to frequent pest management failures, rising input costs, lowering crop yields, and disastrous public health. Insecticide resistance results from multiple factors, mainly indiscriminate insecticide usage and mounted selection pressure on insect populations. Insects respond to insecticide stress at the cellular level by modest yet significant genetic propagations. Transcriptional, co-transcriptional, and post-transcriptional regulatory signals of cells in organisms regulate the intricate processes in gene expressions churning the genetic information in transcriptional units into proteins and non-coding transcripts. Upregulation of detoxification enzymes, notably cytochrome P450s (CYPs), glutathione S-transferases (GSTs), esterases [carboxyl choline esterase (CCE), carboxyl esterase (CarE)] and ATP Binding Cassettes (ABC) at the transcriptional level, modification of target sites, decreased penetration, or higher excretion of insecticides are the noted insect physiological responses. The transcriptional regulatory pathways such as AhR/ARNT, Nuclear receptors, CncC/Keap1, MAPK/CREB, and GPCR/cAMP/PKA were found to regulate the detoxification genes at the transcriptional level. Post-transcriptional changes of non-coding RNAs (ncRNAs) such as microRNAs (miRNA), long non-coding RNAs (lncRNA), and epitranscriptomics, including RNA methylation, are reported in resistant insects. Additionally, genetic modifications such as mutations in the target sites and copy number variations (CNV) are also influencing insecticide resistance. Therefore, these cellular intricacies may decrease insecticide sensitivity, altering the concentrations or activities of proteins involved in insecticide interactions or detoxification. The cellular episodes at the transcriptional and post-transcriptional levels pertinent to insecticide resistance responses in insects are extensively covered in this review. An overview of molecular mechanisms underlying these biological rhythms allows for developing alternative pest control methods to focus on insect vulnerabilities, employing reverse genetics approaches like RNA interference (RNAi) technology to silence particular resistance-related genes for sustained insect management.

Genome-wide analysis of long noncoding RNAs and their association in regulating the metamorphosis of the Sarcophaga peregrina (Diptera: Sarcophagidae)

BACKGROUND

The flesh fly, Sarcophaga peregrina (Diptera: Sarcophagidae), is an important hygiene pest, that causes myiasis in humans and other mammals, typically livestock, and as a vector for various parasitic agents, including bacteria, viruses, and parasites. The role of long non-coding RNAs (lncRNAs) in regulating gene expression during metamorphosis of the flesh fly has not been well established.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we performed genome-wide identification and characterization of lncRNAs from the early pupal stage (1-days pupae), mid-term pupal stage (5-days pupae), and late pupal stage (9-days pupae) of S. peregrina by RNA-seq, and a total of 6921 lncRNAs transcripts were identified. RT-qPCR and enrichment analyses revealed the differentially expressed lncRNAs (DE lncRNAs) that might be associated with insect metamorphosis development. Furthermore, functional analysis revealed that the DE lncRNA (SP_lnc5000) could potentially be involved in regulating the metamorphosis of S. peregrina. RNA interference of SP_lnc5000 caused reduced expression of metamorphosis-related genes in 20-hydroxyecdysone (20E) signaling (Br-c, Ftz-F1), cuticle tanning pathway (TH, DOPA), and chitin related pathway (Cht5). Injection of dsSP_lnc5000 in 3rd instar larvae of S. peregrina resulted in deformed pupae, stagnation of pupal-adult metamorphosis, and a decrease in development time of pupal, pupariation rates and eclosion rates. Hematoxylin-eosin staining (H&E), scanning electron microscope (SEM) observation and cuticle hydrocarbons (CHCs) analysis indicated that SP_lnc5000 had crucial roles in the metamorphosis developmental by modulating pupal cuticular development.

CONCLUSIONS/SIGNIFICANCE

We established that the lncRNA SP_lnc5000 potentially regulates the metamorphosis of S. peregrina by putatively affecting the structure and composition of the pupal cuticle. This study enhances our understanding of lncRNAs as regulators of metamorphosis in S. peregrina, and provide valuable insights into the identification of potential targets for vector control and the development of effective strategies for controlling the spread of myiasis and parasitic diseases.

Examining parent-of-origin effects on transcription and RNA methylation in mediating aggressive behavior in honey bees (Apis mellifera)

Conflict between genes inherited from the mother (matrigenes) and the father (patrigenes) is predicted to arise during social interactions among offspring if these genes are not evenly distributed among offspring genotypes. This intragenomic conflict drives parent-specific transcription patterns in offspring resulting from parent-specific epigenetic modifications. Previous tests of the kinship theory of intragenomic conflict in honey bees (Apis mellifera) provided evidence in support of theoretical predictions for variation in worker reproduction, which is associated with extreme variation in morphology and behavior. However, more subtle behaviors - such as aggression - have not been extensively studied. Additionally, the canonical epigenetic mark (DNA methylation) associated with parent-specific transcription in plant and mammalian model species does not appear to play the same role as in honey bees, and thus the molecular mechanisms underlying intragenomic conflict in this species is an open area of investigation. Here, we examined the role of intragenomic conflict in shaping aggression in honey bee workers through a reciprocal cross design and Oxford Nanopore direct RNA sequencing. We attempted to probe the underlying regulatory basis of this conflict through analyses of parent-specific RNA m6A and alternative splicing patterns. We report evidence that intragenomic conflict occurs in the context of honey bee aggression, with increased paternal and maternal allele-biased transcription in aggressive compared to non-aggressive bees, and higher paternal allele-biased transcription overall. However, we found no evidence to suggest that RNA m6A or alternative splicing mediate intragenomic conflict in this species.

Combined PacBio Iso-Seq and Illumina RNA-Seq Analysis of the Tuta absoluta (Meyrick) Transcriptome and Cytochrome P450 Genes

Tuta absoluta (Meyrick) is a devastating invasive pest worldwide. The abamectin and chlorantraniliprole complex have become an alternative option for chemical control because they can enhance insecticidal activity and delay increased drug resistance. Notably, pests are inevitably resistant to various types of insecticides, and compound insecticides are no exception. To identify potential genes involved in the detoxification of abamectin and chlorantraniliprole complex in T. absoluta, PacBio SMRT-seq transcriptome sequencing and Illumina RNA-seq analysis of abamectin and chlorantraniliprole complex-treated T. absoluta were performed. We obtained 80,492 non-redundant transcripts, 62,762 (77.97%) transcripts that were successfully annotated, and 15,524 differentially expressed transcripts (DETs). GO annotation results showed that most of these DETs were involved in the biological processes of life-sustaining activities, such as cellular, metabolic, and single-organism processes. The KEGG pathway enrichment results showed that the pathways related to glutathione metabolism, fatty acid and amino acid synthesis, and metabolism were related to the response to abamectin and chlorantraniliprole complex in T. absoluta. Among these, 21 P450s were differentially expressed (11 upregulated and 10 downregulated). The qRT-PCR results for the eight upregulated P450 genes after abamectin and chlorantraniliprole complex treatment were consistent with the RNA-Seq data. Our findings provide new full-length transcriptional data and information for further studies on detoxification-related genes in T. absoluta.

Editorial: Role of Non-Coding RNAs in Animals

The importance of non-coding RNAs (ncRNAs), such as microRNAs (miRNA), long non-coding RNAs (lncRNA), and circular RNAs (circRNA), in gene regulation is increasingly being appreciated in many species [...].

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

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

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

BACKGROUND

The Asian tiger mosquito, Aedes albopictus, is one of the most invasive species and a vector of numerous arboviruses. The deleterious effects of long-term and inappropriate use of chemical pesticides have stimulated the exploration of new, environmentally friendly control strategies. Non-coding RNAs (ncRNAs) have been proven to participate in almost all biological processes of insects.

METHODS

In this study, circular RNAs (circRNAs) and microRNAs (miRNAs) covering five developmental stages [egg, early larvae, late larvae, pupae, adult (female and male)] of A. albopictus were obtained using whole-transcriptome sequencing technology. Combined with long non-coding RNAs (lncRNAs) from previous research, circRNA/lncRNA‒miRNA‒mitochondrial RNA (mRNA) networks were constructed.

RESULTS

A total of 1434 circRNAs and 208 miRNAs were identified. More differentially expressed circRNAs (DE circRNAs) and miRNAs (DE miRNAs) were found in the egg versus early larvae comparison group. Functional enrichment analysis demonstrated that most of the circRNA/lncRNA‒miRNA‒mRNA networks were involved in chitin metabolism. Hub genes of each circRNA/lncRNA‒miRNA‒mRNA network were screened out, which can be used as novel targets to disturb the molting process of A. albopictus.

CONCLUSIONS

Regulatory relationships obtained from competing endogenous RNA (ceRNA) networks provide more information to manipulate the metamorphosis process and are helpful for developing effective and sustainable methods to control mosquitoes.

The complex role of transcription factor GAGA in germline death during Drosophila spermatogenesis: transcriptomic and bioinformatic analyses

The GAGA protein (also known as GAF) is a transcription factor encoded by the Trl gene in D. melanogaster. GAGA is involved in the regulation of transcription of many genes at all stages of fly development and life. Recently, we investigated the participation of GAGA in spermatogenesis and discovered that Trl mutants experience massive degradation of germline cells in the testes. Trl underexpression induces autophagic death of spermatocytes, thereby leading to reduced testis size. Here, we aimed to determine the role of the transcription factor GAGA in the regulation of ectopic germline cell death. We investigated how Trl underexpression affects gene expression in the testes. We identified 15,993 genes in three biological replicates of our RNA-seq analysis and compared transcript levels between hypomorphic Trl ^R85/Trl ³⁶² and Oregon testes. A total of 2,437 differentially expressed genes were found, including 1,686 upregulated and 751 downregulated genes. At the transcriptional level, we detected the development of cellular stress in the Trl-mutant testes: downregulation of the genes normally expressed in the testes (indicating slowed or abrogated spermatocyte differentiation) and increased expression of metabolic and proteolysis-related genes, including stress response long noncoding RNAs. Nonetheless, in the Flybase Gene Ontology lists of genes related to cell death, autophagy, or stress, there was no enrichment with GAGA-binding sites. Furthermore, we did not identify any specific GAGA-dependent cell death pathway that could regulate spermatocyte death. Thus, our data suggest that GAGA deficiency in male germline cells leads to an imbalance of metabolic processes, impaired mitochondrial function, and cell death due to cellular stress.

Epigenetics, genomics imprinting and non-coding RNAs

Epigenetic traits are heritable phenotypes caused by alterations in chromosomes rather than DNA sequences. The actual epigenetic expression of the somatic cells of a species is identical, however, they may show distinct subtleties in various cell types in which they may be affected. Several recent studies demonstrated that the epigenetic system plays a very important role in regulating all biological natural processes in the body from birth to death. We outline the essential elements of epigenetics, genomic imprinting, and non-coding RNAs in this mini-review.

Systematic identification and functional analysis of long noncoding RNAs involved in indoxacarb resistance in Spodoptera litura

Long noncoding RNAs (lncRNAs) are noncoding transcripts that are more than 200 nucleotides long. They play essential roles in regulating a variety of biological processes in many species, including insects, and some lncRNAs have been found to be associated with insecticide resistance. However, the characteristics and biological functions of lncRNAs involved in indoxacarb resistance are unknown in Spodoptera litura. We performed RNA sequencing in the SS, InRS, and FInRS of S. litura and identified 11 978 lncRNAs, including 3 136 intergenic lncRNAs, 7 393 intronic lncRNAs, and 1 449 anti-sense lncRNAs. Compared with the SS, 51 lncRNAs were upregulated and 134 lncRNAs were downregulated in the two resistant strains, and 908 differentially expressed mRNAs were predicted as the target genes of the 185 differentially expressed lncRNAs. Further analysis showed that 112 of differentially expressed lncRNAs may be associated with indoxacarb resistance by regulating the expression of 14 P450s, seven CCEs, one GST, six UGTs, five ABC transporters, and 24 cuticle protein genes, and 79 of differentially expressed lncRNAs may regulate the expression of 14 detoxification genes and 19 cuticle protein genes to participate in indoxacarb resistance by sponging 10 microRNAs. Interestingly, 47 of differentially expressed lncRNAs may mediate indoxacarb resistance through both lncRNA-mRNA and lncRNA-miRNA-mRNA regulatory pathways. Furthermore, quantitative PCR, RNA interference, and indoxacarb bioassay analyses indicated that overexpressed LNC_004867 and LNC_006576 were involved in indoxacarb resistance. This study provides comprehensive information for lncRNAs of S. litura, and presents evidence that lncRNAs have key roles in conferring insecticide resistance in S. litura.

Advances in the Immune Regulatory Role of Non-Coding RNAs (miRNAs and lncRNAs) in Insect-Pathogen Interactions

Insects are by far the most abundant and diverse living organisms on earth and are frequently prone to microbial attacks. In other to counteract and overcome microbial invasions, insects have in an evolutionary way conserved and developed immune defense mechanisms such as Toll, immune deficiency (Imd), and JAK/STAT signaling pathways leading to the expression of antimicrobial peptides. These pathways have accessory immune effector mechanisms, such as phagocytosis, encapsulation, melanization, nodulation, RNA interference (RNAi), lysis, autophagy, and apoptosis. However, pathogens evolved strategies that circumvent host immune response following infections, which may have helped insects further sophisticate their immune response mechanisms. The involvement of ncRNAs in insect immunity is undeniable, and several excellent studies or reviews have investigated and described their roles in various insects. However, the functional analyses of ncRNAs in insects upon pathogen attacks are not exhaustive as novel ncRNAs are being increasingly discovered in those organisms. This article gives an overview of the main insect signaling pathways and effector mechanisms activated by pathogen invaders and summarizes the latest findings of the immune modulation role of both insect- and pathogen-encoded ncRNAs, especially miRNAs and lncRNAs during insect–pathogen crosstalk.

Systematic identification and characterization of long noncoding RNAs (lncRNAs) during Aedes albopictus development

Background

Aedes albopictus originated in the tropical forests of Southeast Asia and can currently be found on all continents. As one of the main arboviral vectors, the control of Ae. albopictus requires novel strategies, informed by a deep knowledge of its biology. Little is known regarding mosquito long noncoding RNAs (lncRNAs), which are transcripts longer than 200 nucleotides that lack protein-coding potential and have roles in developmental regulation.

Results

Based on RNA-seq data from five developmental time points, eggs, early larvae, late larvae, pupae, and adults (female and male) of Ae. albopictus, 21,414 lncRNAs were characterized in this study. Differential expression analysis revealed that lncRNAs exhibited developmental stage specificity. The expression of most lncRNAs was upregulated at the onset of metamorphosis developmental stages. More differentially expressed lncRNAs were observed between eggs and early larvae. Weighted gene co-expression network analysis (WGCNA) further confirmed that the expression patterns of lncRNAs were obviously correlated with specific developmental time points. Functional annotation using co-expression analysis revealed that lncRNAs may be involved in the regulation of metamorphic developmental transitions of Ae. albopictus. The hub lncRNAs and hub gene clusters were identified for each module that were highly associated with specific developmental time points.

Conclusions

The results of this study will facilitate future researches to elucidate the regulatory mechanisms of lncRNAs in the development of Ae. albopictus and utilize lncRNAs to assist with mosquito control.

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

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

Phenotypic Plasticity: What Has DNA Methylation Got to Do with It?

How does one genome give rise to multiple, often markedly different, phenotypes in response to an environmental cue? This phenomenon, known as phenotypic plasticity, is common amongst plants and animals, but arguably the most striking examples are seen in insects. Well-known insect examples include seasonal morphs of butterfly wing patterns, sexual and asexual reproduction in aphids, and queen and worker castes of eusocial insects. Ultimately, we need to understand how phenotypic plasticity works at a mechanistic level; how do environmental signals alter gene expression, and how are changes in gene expression translated into novel morphology, physiology and behaviour? Understanding how plasticity works is of major interest in evolutionary-developmental biology and may have implications for understanding how insects respond to global change. It has been proposed that epigenetic mechanisms, specifically DNA methylation, are the key link between environmental cues and changes in gene expression. Here, we review the available evidence on the function of DNA methylation of insects, the possible role(s) for DNA methylation in phenotypic plasticity and also highlight key outstanding questions in this field as well as new experimental approaches to address these questions.

The long non-coding RNA landscape of Candida yeast pathogens

Long non-coding RNAs (lncRNAs) constitute a poorly studied class of transcripts with emerging roles in key cellular processes. Despite efforts to characterize lncRNAs across a wide range of species, these molecules remain largely unexplored in most eukaryotic microbes, including yeast pathogens of the Candida clade. Here, we analyze thousands of publicly available sequencing datasets to infer and characterize the lncRNA repertoires of five major Candida pathogens: Candida albicans, Candida tropicalis, Candida parapsilosis, Candida auris and Candida glabrata. Our results indicate that genomes of these species encode hundreds of lncRNAs that show levels of evolutionary constraint intermediate between those of intergenic genomic regions and protein-coding genes. Despite their low sequence conservation across the studied species, some lncRNAs are syntenic and are enriched in shared sequence motifs. We find co-expression of lncRNAs with certain protein-coding transcripts, hinting at potential functional associations. Finally, we identify lncRNAs that are differentially expressed during infection of human epithelial cells for four of the studied species. Our comprehensive bioinformatic analyses of Candida lncRNAs pave the way for future functional characterization of these transcripts.

Long non-coding RNAs (lncRNAs) play roles in key cellular processes, but remain largely unexplored in fungal pathogens such as Candida. Here, Hovhannisyan and Gabaldón analyze thousands of sequencing datasets to infer and characterize the lncRNA repertoires of five Candida species, paving the way for their future functional characterization.

Long non-coding RNAs as possible therapeutic targets in protozoa, and in Schistosoma and other helminths

Long non-coding RNAs (lncRNAs) emerged in the past 20 years due to massive amounts of scientific data regarding transcriptomic analyses. They have been implicated in a plethora of cellular processes in higher eukaryotes. However, little is known about lncRNA possible involvement in parasitic diseases, with most studies only detecting their presence in parasites of human medical importance. Here, we review the progress on lncRNA studies and their functions in protozoans and helminths. In addition, we show an example of knockdown of one lncRNA in Schistosoma mansoni, SmLINC156349, which led to in vitro parasite adhesion, motility, and pairing impairment, with a 20% decrease in parasite viability and 33% reduction in female oviposition. Other observed phenotypes were a decrease in the proliferation rate of both male and female worms and their gonads, and reduced female lipid and vitelline droplets that are markers for well-developed vitellaria. Impairment of female worms' vitellaria in SmLINC156349-silenced worms led to egg development deficiency. All those results demonstrate the great potential of the tools and methods to characterize lncRNAs as potential new therapeutic targets. Further, we discuss the challenges and limitations of current methods for studying lncRNAs in parasites and possible solutions to overcome them, and we highlight the future directions of this exciting field.

There and Back Again: Hox Clusters Use Both DNA Strands

Bilaterian animals operate the clusters of Hox genes through a rich repertoire of diverse mechanisms. In this review, we will summarize and analyze the accumulated data concerning long non-coding RNAs (lncRNAs) that are transcribed from sense (coding) DNA strands of Hox clusters. It was shown that antisense regulatory RNAs control the work of Hox genes in cis and trans, participate in the establishment and maintenance of the epigenetic code of Hox loci, and can even serve as a source of regulatory peptides that switch cellular energetic metabolism. Moreover, these molecules can be considered as a force that consolidates the cluster into a single whole. We will discuss the examples of antisense transcription of Hox genes in well-studied systems (cell cultures, morphogenesis of vertebrates) and bear upon some interesting examples of antisense Hox RNAs in non-model Protostomia.