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

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.

Long non-coding RNAs as promising targets for controlling disease vector mosquitoes

Hematophagous female mosquitoes are important vectors of numerous devastating human diseases, posing a major public health threat. Effective prevention and control of mosquito-borne diseases rely considerably on progress in understanding the molecular mechanisms of various life activities, and accordingly, the molecules that regulate the various life activities of mosquitoes are potential targets for implementing future vector control strategies. Many long non-coding RNAs (lncRNAs) have been identified in mosquitoes and significant progress has been made in determining their functions. Here, we present a comprehensive overview of the research advances on mosquito lncRNAs, including their molecular identification, function, and interaction with other non-coding RNAs, as well as their synergistic regulatory roles in mosquito life activities. We also highlight the potential roles of competitive endogenous RNAs in mosquito growth and development, as well as in insecticide resistance and virus-host interactions. Insights into the biological functions and mechanisms of lncRNAs in mosquito life activities, viral replication, pathogenesis, and transmission will contribute to the development of novel drugs and safe vaccines.

Non-Coding RNAs Potentially Involved in Pyrethroid Resistance of Anopheles funestus Population in Western Kenya

Backgrounds The resurgence of Anopheles funestus , a dominant vector of human malaria in western Kenya was partly attributed to insecticide resistance. However, evidence on the molecular basis of pyrethroid resistance in western Kenya is limited. Noncoding RNAs (ncRNAs) form a vast class of RNAs that do not code for proteins and are ubiquitous in the insect genome. Here, we demonstrated that multiple ncRNAs could play a potential role in An. funestus resistance to pyrethroid in western Kenya. Materials and Methods Anopheles funestus mosquitoes were sampled by aspiration methods in Bungoma, Teso, Siaya, Port Victoria and Kombewa in western Kenya. The F1 progenies were exposed to deltamethrin (0.05%), permethrin (0.75%), DDT (4%) and pirimiphos-methyl (0.25%) following WHO test guidelines. A synergist assay using piperonyl butoxide (PBO) (4%) was conducted to determine cytochrome P450s' role in pyrethroid resistance. RNA-seq was conducted on a combined pool of specimens that were resistant and unexposed, and the results were compared with those of the FANG susceptible strain. This approach aimed to uncover the molecular mechanisms underlying pyrethroid resistance. Results Pyrethroid resistance was observed in all the sites with an average mortality rate of 57.6%. Port Victoria had the highest level of resistance to permethrin (MR=53%) and deltamethrin (MR=11%) pyrethroids. Teso had the lowest level of resistance to permethrin (MR=70%) and deltamethrin (MR=87%). Resistance to DDT was observed only in Kombewa (MR=89%) and Port Victoria (MR=85%). A full susceptibility to P-methyl (0.25%) was observed in all the sites. PBO synergist assay revealed high susceptibility (>98%) to the pyrethroids in all the sites except for Port Victoria (MR=96%, n=100). Whole transcriptomic analysis showed that most of the gene families associated with pyrethroid resistance comprised non-coding RNAs (67%), followed by imipenemase (10%),cytochrome P450s (6%), cuticular proteins (5%), olfactory proteins (4%), glutathione S-transferases (3%), UDP-glycosyltransferases (2%), ATP-binding cassettes (2%) and carboxylesterases(1%). Conclusions This study unveils the molecular basis of insecticide resistance in An. funestus in western Kenya, highlighting for the first time the potential role of non-coding RNAs in pyrethroid resistance. Targeting non-coding RNAs for intervention development could help in insecticide resistance management.

Regulatory Roles of Long Non-Coding RNAs Relevant to Antioxidant Enzymes and Immune Responses of Apis cerana Larvae Following Ascosphaera apis Invasion

Long non-coding RNAs (lncRNAs) play an essential part in controlling gene expression and a variety of biological processes such as immune defense and stress-response. However, whether and how lncRNAs regulate responses of Apis cerana larvae to Ascosphaera apis invasion has remained unclear until now. Here, the identification and structural analysis of lncRNAs in the guts of A. cerana worker larvae were conducted, and the expression profile of larval lncRNAs during the A. apis infection process was then analyzed, followed by an investigation of the regulatory roles of differentially expressed lncRNAs (DElncRNAs) in the host response. In total, 76 sense lncRNAs, 836 antisense lncRNAs, 184 intron lncRNAs, 362 bidirectional lncRNAs, and 2181 intron lncRNAs were discovered in the larval guts. Additionally, 30 known and 9 novel lncRNAs were potential precursors for 36 and 11 miRNAs, respectively. In the three comparison groups, 386, 351, and 272 DElncRNAs were respectively identified, indicating the change in the overall expression pattern of host lncRNAs following the A. apis invasion. Analysis of cis-acting effect showed that DElncRNAs in the 4-, 5-, and 6-day-old comparison groups putatively regulated 55, 30, and 20 up- and down-stream genes, respectively, which were involved in a series of crucial functional terms and pathways, such as MAPK signaling pathway, and cell process. Analysis showed that 31, 8, and 11 DElncRNAs as potential antisense lncRNAs may interact with 26, 8, and 9 sense-strand mRNAs. Moreover, investigation of the competing endogenous RNA (ceRNA) network indicated that 148, 283, and 257 DElncRNAs were putatively regulated. The expression of target genes by targeting corresponding DEmiRNAs included those associated with antioxidant enzymes and immune responses. These results suggested that DElncRNAs played a potential part in the larval guts responding to the A. apis infection through a cis-acting manner and ceRNA mechanisms. Our findings deepen our understanding of interactions between A. cerana larvae and A. apis and offer a basis for clarifying the DElncRNA-mediated mechanisms underlying the host response to fungal invasion.

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.

Cuticular competing endogenous RNAs regulate insecticide penetration and resistance in a major agricultural pest

BACKGROUND

The continuously developing pesticide resistance is a great threat to agriculture and human health. Understanding the mechanisms of insecticide resistance is a key step in dealing with the phenomenon. Insect cuticle is recently documented to delay xenobiotic penetration which breaks the previous stereotype that cuticle is useless in insecticide resistance, while the underlying mechanism remains scarce.

RESULTS

Here, we find the integument contributes over 40.0% to insecticide resistance via different insecticide delivery strategies in oriental fruit fly. A negative relationship exists between cuticle thickening and insecticide penetration in resistant/susceptible, also in field strains of oriental fruit fly which is a reason for integument-mediated resistance. Our investigations uncover a regulator of insecticide penetration that miR-994 mimic treatment causes cuticle thinning and increases susceptibility to malathion, whereas miR-994 inhibitor results in opposite phenotypes. The target of miR-994 is a most abundant cuticle protein (CPCFC) in resistant/susceptible integument expression profile, which possesses capability of chitin-binding and influences the cuticle thickness-mediated insecticide penetration. Our analyses find an upstream transcriptional regulatory signal of miR-994 cascade, long noncoding RNA (lnc19419), that indirectly upregulates CPCFC in cuticle of the resistant strain by sponging miR-994. Thus, we elucidate the mechanism of cuticular competing endogenous RNAs for regulating insecticide penetration and demonstrate it also exists in field strain of oriental fruit fly.

CONCLUSIONS

We unveil a regulatory axis of lnc19419 ~ miR-994 ~ CPCFC on the cuticle thickness that leads to insecticide penetration resistance. These findings indicate that competing endogenous RNAs regulate insecticide resistance by modulating the cuticle thickness and provide insight into the resistance mechanism in insects.

Detoxification enzyme is involved in the temperature effect on the toxicity of tetrachlorantraniliprole to Plutella xylostella

The efficacy of insecticides is usually influenced by temperature. Insecticides can be divided into "positive", "negative" and "non-effect" temperature coefficient insecticides (TCI). To assess the temperature-dependent effect of tetrachlorantraniliprole (TET) on Plutella xylostella Linnaeus and to elucidate the mechanism of temperature affects TET toxicity, we determined the toxicity of TET against P. xylostella from 15 °C to 35 °C by leaf dipping method. Moreover, we compared the transcriptome data of the third-instar larvae treated by TET, chlorfenapyr (CHL, non-effect TCI), and the control group at 15, 25, 35 °C, respectively. The results showed that the toxicity of TET against P. xylostella increased with increasing temperature from 15 °C to 35 °C. A total of 21 differential expressed genes (DEGs) of detoxification enzymes were screened by RNA-seq, in which 10 up-regulated genes (3 UGTs, 2 GSTs, 5 P450s) may involve the positive temperature effect of TET, and their expression patterns were consistent with qPCR results. Furthermore, the enzyme activities of GSTs and UGTs significantly increased after TET was treated at 15 °C. Especially, the temperature coefficient (TC) of TET was significantly reduced mixed with UGTs enzyme inhibitor 5-NI. Overall, TET showed higher insecticidal activity with increasing temperature, in which detoxifying enzymes associated with regulation of the positive temperature effect of TET on P. xylostella, such as UGTs, GSTs and P450s, are strongly involved. The transcriptome data provide in-depth information to understand the TET mechanism against diamondback moth. Most importantly, we identified detoxification enzymes that might be involved in regulating TET's positive temperature effect process, and contributed to efficient pest management.

Genome-wide identification of lncRNAs associated with viral infection in Spodoptera frugiperda

The role of lncRNAs in immune defence has been demonstrated in many multicellular and unicellular organisms. However, investigation of the identification and characterization of long non-coding RNAs (lncRNAs) involved in the insect immune response is still limited. In this study, we used RNA sequencing (RNA-seq) to investigate the expression profiles of lncRNAs and mRNAs in the fall armyworm Spodoptera frugiperda in response to virus infection. To assess the tissue- and virus-specificity of lncRNAs, we analysed and compared their expression profiles in haemocytes and fat body of larvae infected with two entomopathogenic viruses with different lifestyles, i.e. the polydnavirus HdIV (Hyposoter didymator IchnoVirus) and the densovirus JcDV (Junonia coenia densovirus). We identified 1883 candidate lncRNAs, of which 529 showed differential expression following viral infection. Expression profiles differed considerably between samples, indicating that many differentially expressed (DE) lncRNAs showed virus- and tissue-specific expression patterns. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment and target prediction analyses indicated that DE-LncRNAs were mainly enriched in metabolic process, DNA replication and repair, immune response, metabolism of insect hormone and cell adhesion. In addition, we identified three DE-lncRNAs potentially acting as microRNA host genes, suggesting that they participate in gene regulation by producing miRNAs in response to virus infection. This study provides a catalogue of lncRNAs expressed in two important immune tissues and potential insight into their roles in the antiviral defence in S. frugiperda. The results may help future in-depth functional studies to better understand the biological function of lncRNAs in interaction between viruses and the fall armyworm.

Genome-wide identification and functional analysis of long non-coding RNAs in Chilo suppressalis reveal their potential roles in chlorantraniliprole resistance

Long non-coding RNAs, referred to as lncRNAs, perform essential functions in some biological processes, including reproduction, metamorphosis, and other critical life functions. Yet, lncRNAs are poorly understood in pesticide resistance, and no reports to date have characterized which lncRNAs are associated with chlorantraniliprole resistance in Chilo suppressalis. Here, RNA-seq was performed on two strains of C. suppressalis exposed to chlorantraniliprole: one is a susceptible strain (S), and the other is a resistant strain (R). In total, 3,470 lncRNAs were identified from 40,573 merged transcripts in six libraries, including 1,879 lincRNAs, 245 intronic lncRNAs, 853 sense lncRNAs, and 493 antisense lncRNAs. Moreover, differential expression analysis revealed 297 and 335 lncRNAs upregulated in S and R strains, respectively. Differentially expressed (DE) lncRNAs are usually assumed to be involved in the chlorantraniliprole resistance in C. suppressalis. As potential targets, adjacent protein-coding genes (within <1000 kb range upstream or downstream of DE lncRNAs), especially detoxification enzyme genes (cytochrome P450s, carboxyl/cholinesterases/esterases, and ATP-binding cassette transporter), were analyzed. Furthermore, the strand-specific RT-PCR was conducted to confirm the transcript orientation of randomly selected 20 DE lincRNAs, and qRT-PCR was carried out to verify the expression status of 8 out of them. MSTRG.25315.3, MSTRG.25315.6, and MSTRG.7482.1 were upregulated in the R strain. Lastly, RNA interference and bioassay analyses indicated overexpressed lincRNA MSTRG.7482.1 was involved in chlorantraniliprole resistance. In conclusion, we represent, for the first time, the genome-wide identification of chlorantraniliprole-resistance-related lncRNAs in C. suppressalis. It elaborates the views underlying the mechanism conferring chlorantraniliprole resistance in lncRNAs.

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.

Analysis of Long Non-Coding RNA-Mediated Regulatory Networks of Plutella xylostella in Response to Metarhizium anisopliae Infection

Long non-coding RNAs (lncRNAs) represent a diverse class of RNAs that are structurally similar to messenger RNAs (mRNAs) but do not encode proteins. Growing evidence suggests that in response to biotic and abiotic stresses, the lncRNAs play crucial regulatory roles in plants and animals. However, the potential role of lncRNAs during fungal infection has yet to be characterized in Plutella xylostella, a devastating pest of cruciferous crops. In the current study, we performed a strand-specific RNA sequencing of Metarhizium anisopliae-infected (Px36hT, Px72hT) and uninfected (Px36hCK, Px72hCK) P. xylostella fat body tissues. Comprehensive bioinformatic analysis revealed a total of 5665 and 4941 lncRNAs at 36 and 72-h post-infection (hpi), including 563 (Px36hT), 532 (Px72hT) known and 5102 (Px36hT), 4409 (Px72hT) novel lncRNA transcripts. These lncRNAs shared structural similarities with their counterparts in other species, including shorter exon and intron length, fewer exon numbers, and a lower expression profile than mRNAs. LncRNAs regulate the expression of neighboring protein-coding genes by acting in a cis and trans manner. Functional annotation and pathway analysis of cis-acting lncRNAs revealed their role in several immune-related genes, including Toll, serpin, transferrin, βGRP etc. Furthermore, we identified multiple lncRNAs acting as microRNA (miRNA) precursors. These miRNAs can potentially regulate the expression of mRNAs involved in immunity and development, suggesting a crucial lncRNA-miRNA-mRNA complex. Our findings will provide a genetic resource for future functional studies of lncRNAs involved in P. xylostella immune responses to M. anisopliae infection and shed light on understanding insect host-pathogen interactions.

Physiological and Molecular Response Modifications by Ultraviolet-C Radiation in Plutella xylostella and Its Compatibility with Cordyceps fumosorosea

Ultraviolet-C (UV-C) radiation significantly impacts living organisms. UV-C radiation can also be used as a pest management tool. Therefore, this study was designed to investigate the effect of UV-C radiation on the physiology and gene expression level of Plutella xylostella, a destructive vegetable pest. Results showed that, after exposure to UV-C radiation for 3, 6, 12, and 24 h, the activity of SOD (superoxide dismutase) and CAT (catalase) of P. xylostella increased, while the activity of PPO (polyphenol oxidase), POD (peroxidase), AChE (acetylcholinesterase), CarE (carboxylesterase), and ACP (acid phosphatase) decreased with increased exposure time. Correlation coefficient analyses indicated that the activity of CAT correlated positively, while PPO and CarE correlated negatively, with exposure time. Gene regulation analysis via qRT-PCR confirmed a significant increase in regulation in CAT, CarE, and PPO-related genes. We also investigated the effect of UV-C exposure on the virulence of Cordyceps fumosorosea against P. xylostella. Here, results indicated that when the fungal treatment was applied to larvae before UV-C radiation, the virulence of C. fumosorosea was significantly reduced. However, this decline in virulence of C. fumosorosea due to UV-C exposure remained only for one generation, and no effect was observed on secondary infection. On the other hand, when larvae were exposed to UV-C radiation before fungal application, the mortality rate significantly increased as the exposure time to UV-C radiation increased. From the current study, it could be concluded that UV-C exposure suppressed the immunity to P. xylostella, which later enhanced the virulence of entomopathogenic fungi. Moreover, the study also suggested that UV irradiation is an effective pest management tool that could be incorporated into pest management strategies, which could help reduce pesticide application, be economically beneficial for the farmer, and be environmentally safe.

Genome-wide identification of long non-coding (lncRNA) in Nilaparvata lugens's adaptability to resistant rice

Background

The brown planthopper (BPH), Nilaparvata lugens (Stål), is a very destructive pest that poses a major threat to rice plants worldwide. BPH and rice have developed complex feeding and defense strategies in the long-term co-evolution.

Methods

To explore the molecular mechanism of BPH's adaptation to resistant rice varieties, the lncRNA expression profiles of two virulent BPH populations were analyzed. The RNA-seq method was used to obtain the lncRNA expression data in TN1 and YHY15.

Results

In total, 3,112 highly reliable lncRNAs in TN1 and YHY15 were identified. Compared to the expression profiles between TN1 and YHY15, 157 differentially expressed lncRNAs, and 675 differentially expressed mRNAs were identified. Further analysis of the possible regulation relationships between differentially expressed lncRNAs and differentially expressed mRNAs, identified three pair antisense targets, nine pair cis-regulation targets, and 3,972 pair co-expressed targets. Function enriched found arginine and proline metabolism, glutathione metabolism, and carbon metabolism categories may significantly affect the adaptability in BPH when it is exposed to susceptible and resistant rice varieties. Altogether, it provided scientific data for the study of lncRNA regulation of brown planthopper resistance to rice. These results are helpful in the development of new control strategies for host defense against BPH and breeding rice for high yield.

Genome-Wide Identification of the Long Noncoding RNAs of Tribolium castaneum in Response to Terpinen-4-ol Fumigation

Simple Summary

Long noncoding RNAs (lncRNAs) are important regulatory factors in multiple biological processes, including genomic imprinting, cancer, RNA interference, and protein translation. Several lncRNAs can respond to insecticides. However, lncRNA functions associated with terpinen-4-ol resistance in the red flour beetle (Tribolium castaneum) have not yet been identified. In previous work, we found terpinen-4-ol to have strong fumigation activity against store-product pests. As a pesticide from plants, terpinen-4-ol shows nearly no residual danger to the environment; however, resistance is inevitable if people use terpinen-4-ol immoderately. To avoid resistance to terpinen-4-ol occurring in the red flour beetle, we deeply sequenced and tried to find some lncRNAs that can regulate target mRNA expression to reduce terpinen-4-ol.

Abstract

Long noncoding RNAs (lncRNAs) are important regulatory factors in multiple biological processes, and several lncRNAs are known to respond to insecticides. However, the lncRNA functions that are associated with terpinen-4-ol resistance in the red flour beetle (Tribolium castaneum) have not yet been identified. In this study, we determined the differentially transcribed lncRNAs between fumigated and control experimental groups. In the six libraries that underwent RNA sequencing, 34,546 transcripts were identified, including 8267 novel lncRNAs, 4155 novel mRNAs, 1151 known lncRNAs, and 20,973 known mRNAs. Among these, we found that the expression of 1858 mRNAs and 1663 lncRNAs was significantly different in the fumigated group compared with the control group. Among the differentially transcribed lncRNAs, 453 were up-regulated and 1210 were down-regulated lncRNAs. In addition, we identified the regulatory function targets of the lncRNAs. Functionally, all lncRNAs and target genes associated with terpinen-4-ol metabolism were enriched in several metabolic pathways, like the ATP-binding cassette transporter, pentose interconversion, and glucuronate interconversion. To the best of our knowledge, this study represents the first global identification of lncRNAs and their potential association with terpinen-4-ol metabolism in the red flour beetle. These results will provide reference information for studies on the resistance to terpinen-4-ol and other essential oil compounds and chemical pesticides, as well as an understanding of other biological processes in T. castaneum.

Transcriptome analysis reveals the potential roles of long non-coding RNAs in feed efficiency of chicken

Feed efficiency is an important economic trait and reduces the production costs per unit of animal product. Up to now, few studies have conducted transcriptome profiling of liver tissue in feed efficiency-divergent chickens (Ross vs native breeds). Also, molecular mechanisms contributing to differences in feed efficiency are not fully understood, especially in terms of long non-coding RNAs (lncRNAs). Hence, transcriptome profiles of liver tissue in commercial and native chicken breeds were analyzed. RNA-Seq data along with bioinformatics approaches were applied and a series of lncRNAs and target genes were identified. Furthermore, protein-protein interaction network construction, co-expression analysis, co-localization analysis of QTLs and functional enrichment analysis were used to functionally annotate the identified lncRNAs. In total, 2,290 lncRNAs were found (including 1,110 annotated, 593 known and 587 novel), of which 53 (including 39 known and 14 novel), were identified as differentially expressed genes between two breeds. The expression profile of lncRNAs was validated by RT-qPCR. The identified novel lncRNAs showed a number of characteristics similar to those of known lncRNAs. Target prediction analysis showed that these lncRNAs have the potential to act in cis or trans mode. Functional enrichment analysis of the predicted target genes revealed that they might affect the differences in feed efficiency of chicken by modulating genes associated with lipid metabolism, carbohydrate metabolism, growth, energy homeostasis and glucose metabolism. Some gene members of significant modules in the constructed co-expression networks were reported as important genes related to feed efficiency. Co-localization analysis of QTLs related to feed efficiency and the identified lncRNAs suggested several candidates to be involved in residual feed intake. The findings of this study provided valuable resources to further clarify the genetic basis of regulation of feed efficiency in chicken from the perspective of lncRNAs.

Common Features in lncRNA Annotation and Classification: A Survey

Long non-coding RNAs (lncRNAs) are widely recognized as important regulators of gene expression. Their molecular functions range from miRNA sponging to chromatin-associated mechanisms, leading to effects in disease progression and establishing them as diagnostic and therapeutic targets. Still, only a few representatives of this diverse class of RNAs are well studied, while the vast majority is poorly described beyond the existence of their transcripts. In this review we survey common in silico approaches for lncRNA annotation. We focus on the well-established sets of features used for classification and discuss their specific advantages and weaknesses. While the available tools perform very well for the task of distinguishing coding sequence from other RNAs, we find that current methods are not well suited to distinguish lncRNAs or parts thereof from other non-protein-coding input sequences. We conclude that the distinction of lncRNAs from intronic sequences and untranslated regions of coding mRNAs remains a pressing research gap.

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.

Identification and the potential roles of long non-coding RNAs in regulating acetyl-CoA carboxylase ACC transcription in spirotetramat-resistant Aphis gossypii

Long non-coding RNAs (lncRNAs) represent the largest class of non-coding transcripts. They act a pivotal part in various insect developmental processes and stress responses. However, the investigation of lncRNA functions in insecticide resistant remains at an early phase. Herein, we conducted whole-transcriptome RNA sequencing for two cotton aphid (Aphis gossypii Glover) strains, i.e., insecticide-susceptible (SS) and spirotetramat-resistant (SR). We discovered 6059 lncRNAs in the RNA-Seq data, and 874 lncRNAs showed differential expression. In addition, 5 lncRNAs among 874 lncRNAs were predicted as targets of acetyl-CoA carboxylase (ACC). Reverse transcription real-time quantitative PCR (RT-qPCR) combined with RNA interference (RNAi) confirmed that selected ACC lncRNA was related to the expression of ACC. Moreover, we also identified two transcription factors, i.e., C/EBP and C/EBPzeta, that regulate the transcription level of ACC lncRNA. These results provide a good basis for the study of cotton aphid lncRNA functions in insecticide resistance development.

Regulation of GSTu1-mediated insecticide resistance in Plutella xylostella by miRNA and lncRNA

The evolution of resistance to insecticides is well known to be closely associated with the overexpression of detoxifying enzymes. Although the role of glutathione S-transferase (GST) genes in insecticide resistance has been widely reported, the underlying regulatory mechanisms are poorly understood. Here, one GST gene (GSTu1) and its antisense transcript (lnc-GSTu1-AS) were identified and cloned, and both of them were upregulated in several chlorantraniliprole-resistant Plutella xylostella populations. GSTu1 was confirmed to be involved in chlorantraniliprole resistance by direct degradation of this insecticide. Furthermore, we demonstrated that lnc-GSTu1-AS interacted with GSTu1 by forming an RNA duplex, which masked the binding site of miR-8525-5p at the GSTu1-3′UTR. In summary, we revealed that lnc-GSTu1-AS maintained the mRNA stability of GSTu1 by preventing its degradation that could have been induced by miR-8525-5p and thus increased the resistance of P. xylostella to chlorantraniliprole. Our findings reveal a new noncoding RNA-mediated pathway that regulates the expression of detoxifying enzymes in insecticide-resistant insects and offer opportunities for the further understanding of the mechanisms of insecticide and drug resistance.

The development of insecticide resistance in insect pests is a worldwide concern and a major problem in agriculture. Understanding the genetics of insecticide resistance is critical for effective crop protection. Plutella. xylostella (L.), a major pest of cruciferous crops, has developed resistance to almost all kinds of insecticide, and has become one of the most resistant pests in the world. Overexpression of detoxification enzymes is closely associated with insecticide resistance, but researches on their regulatory mechanism are still very limited. Here, GSTu1 was identified to be upregulated in several chlorantraniliprole-resistant P. xylostella populations and was confirmed to be involved in chlorantraniliprole resistance by direct degradation of this insecticide. Further, lnc-GSTu1-AS transcribed from the opposite DNA strand to GSTu1 was identified to be able to enhance the mRNA stability of GSTu1 by blocking miRNA activity, and thus increased the resistance of P. xylostella to chlorantraniliprole. The results provide further insights into the mechanisms underlying metabolic resistance.

Comprehensive Transcriptome Analyses in Sea Louse Reveal Novel Delousing Drug Responses Through MicroRNA regulation

The role of miRNAs in pharmacological responses through gene regulation related to drug metabolism and the detoxification system has recently been determined for terrestrial species. However, studies on marine ectoparasites have scarcely been conducted to investigate the molecular mechanisms of pesticide resistance. Herein, we explored the sea louse Caligus rogercresseyi miRNome responses exposed to delousing drugs and the interplaying with coding/non-coding RNAs. Drug sensitivity in sea lice was tested by in vitro bioassays for the pesticides azamethiphos, deltamethrin, and cypermethrin. Ectoparasites strains with contrasting susceptibility to these compounds were used. Small-RNA sequencing was conducted, identifying 2776 novel annotated miRNAs, where 163 mature miRNAs were differentially expressed in response to the drug testing. Notably, putative binding sites for miRNAs were found in the ADME genes associated with the drugs' absorption, distribution, metabolism, and excretion. Interactions between the miRNAs and long non-coding RNAs (lncRNAs) were also found, suggesting putative molecular gene regulation mechanisms. This study reports putative miRNAs correlated to the coding/non-coding RNAs modulation, revealing novel pharmacological mechanisms associated with drug resistance in sea lice species.

Identification of differentially expressed microRNAs under imidacloprid exposure in Sitobion miscanthi

Imidacloprid is a neonicotinoid that targets sucking pests, such as aphids and the green leaf bug and has been widely applied in wheat fields to control wheat aphids in China. To investigate the involvement of miRNAs in imidacloprid resistance, we sequenced small RNA libraries of Sitobion miscanthi Fabricius, across two different treatments using Illumina short-read sequencing technology. As a result, 265 microRNAs (miRNAs), of which 242 were known and 23 were novel, were identified. Quantitative analysis of miRNA levels showed that 23 miRNAs were significantly up-regulated, and 54 miRNAs were significantly down-regulated in the nymphs of S. miscanthi treated with imidacloprid in comparison with those of the control. Modulation of the abundances of differentially expressed miRNAs, smi-miR-316, smi-miR-1000, and smi-miR-iab-4 by the addition of the corresponding antagomir/inhibitor to the artificial diet significantly changed the susceptibility of S. miscanthi to imidacloprid. Subsequently, the post-transcriptional regulatory mechanism was conducted, smi-miR-278 and smi-miR-316 were confirmed to be participated in the post-transcriptional regulation of nAChRα1A and CYP4CJ6, respectively. The results suggested that miRNAs differentially expressed in response to imidacloprid could play a critical regulatory role in the metabolism of S. miscanthi to imidacloprid.

Comparison of long non-coding RNAs in adipose and muscle tissues between seven indigenous Chinese and the Yorkshire pig breeds

lncRNAs play crucial roles in fat metabolism in animals. Previously, we have compared the mRNA transcriptome profiles between seven fat-type Chinese pig breeds and one lean-type Western breed (Yorkshire, YY). The associations between differentially expressed (DE) genes and phenotypical traits were investigated. In the present study, to further explore the underlying regulatory mechanisms, lncRNAs were sequenced and compared between YY and Chinese indigenous breeds. The results showed 9114 and 7538 DE lncRNAs between at least one Chinese breed and the YY breed in the adipose and muscle tissue respectively. KEGG enrichment analysis revealed that the target genes of these DE lncRNAs mainly influenced the glucolipid metabolism, which is an important process affecting meat quality. Correlation analyses between the DE lncRNA and DE mRNA genes related to meat quality and growth traits were performed. The results showed that LTCONS_00073280 was associated with intramuscular fat content. Four lncRNAs (LTCONS_00101781, LTCONS_00037879, LTCONS_00088260 and LTCONS-00128343) might mediate backfat thickness. Overall, this study provides candidate lncRNAs that potentially affect meat quality, which might be useful for molecular breeding of pig breeds in future.

Transcriptome analysis reveals potential function of long non-coding RNAs in 20-hydroxyecdysone regulated autophagy in Bombyx mori

Background

20-hydroxyecdysone (20E) plays important roles in insect molting and metamorphosis. 20E-induced autophagy has been detected during the larval–pupal transition in different insects. In Bombyx mori, autophagy is induced by 20E in the larval fat body. Long non-coding RNAs (lncRNAs) function in various biological processes in many organisms, including insects. Many lncRNAs have been reported to be potential for autophagy occurrence in mammals, but it has not been investigated in insects.

Results

RNA libraries from the fat body of B. mori dissected at 2 and 6 h post-injection with 20E were constructed and sequenced, and comprehensive analysis of lncRNAs and mRNAs was performed. A total of 1035 lncRNAs were identified, including 905 lincRNAs and 130 antisense lncRNAs. Compared with mRNAs, lncRNAs had longer transcript length and fewer exons. 132 lncRNAs were found differentially expressed at 2 h post injection, compared with 64 lncRNAs at 6 h post injection. Thirty differentially expressed lncRNAs were common at 2 and 6 h post-injection, and were hypothesized to be associated with the 20E response. Target gene analysis predicted 6493 lncRNA-mRNA cis pairs and 42,797 lncRNA-mRNA trans pairs. The expression profiles of LNC_000560 were highly consistent with its potential target genes, Atg4B, and RNAi of LNC_000560 significantly decreased the expression of LNC_000560 and Atg4B. These results indicated that LNC_000560 was potentially involved in the 20E-induced autophagy of the fat body by regulating Atg4B.

Conclusions

This study provides the genome-wide identification and functional characterization of lncRNAs associated with 20E-induced autophagy in the fat body of B. mori. LNC_000560 and its potential target gene were identified to be related to 20-regulated autophagy in B. mori. These results will be helpful for further studying the regulatory mechanisms of lncRNAs in autophagy and other biological processes in this insect model.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07692-1.

Genome-wide identification of long non-coding RNAs (lncRNAs) associated with malathion resistance in Bactrocera dorsalis

BACKGROUND

Long non-coding RNAs (lncRNAs) play important roles in the regulation of biological processes and have been identified in many species including insects. However, the association between lncRNAs and pesticide resistance in insect species such as Bactrocera dorsalis is unknown.

RESULTS

RNA-seq was performed on malathion resistant (MR1) and susceptible (MS) strains of B. dorsalis and a total of 6171 lncRNAs transcripts were identified. These included 3728 lincRNAs, 653 antisense lncRNAs, 1402 intronic lncRNAs, and 388 sense lncRNAs. A total of 40 and 52 upregulated lncRNAs were found in females and males of the MR1 strain compared to 54 and 49 in the same sexes of the MS strain, respectively. Twenty-seven of these lncRNAs showed the same trend of expression in both females and males in the MR1 strain, in which 15 lncRNAs were upregulated and 12 were downregulated. RT-qPCR results indicated that the differentially expressed lncRNAs were associated with malathion resistance. The lnc15010.10 and lnc3774.2 were highly expressed in the cuticle of the MR1 strain, indicating that these two lncRNAs may be related to malathion resistance. RNAi of lnc3774.2 and a bioassay showed that malathion resistance was possibly influenced by changes in the B. dorsalis cuticle.

CONCLUSION

LncRNAs of B. dorsalis potentially related to the malathion resistance were identified. Two lncRNAs appear to influence malathion resistance via modulating the structure, or components, of the cuticle. © 2021 Society of Chemical Industry.

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.

Full length sequencing reveals novel transcripts of detoxification genes along with related alternative splicing events and lncRNAs in Phyllotreta striolata

The striped flea beetle, Phyllotreta striolata (Fabricius), damages crops in the Brassicaceae. The genetic data for this pest are insufficient to reveal its insecticide resistance mechanisms or to develop molecular markers for resistance monitoring. We used PacBio Iso-Seq technology to sequence the full-length transcriptome of P. striolata. After isoform sequence clustering and removal of redundant transcripts, a total of 41,293 transcripts were obtained, and 35,640 of these were annotated in the database of gene products. Structure analysis uncovered 4,307 alternative splicing events, and 3,836 sequences were recognized as lncRNAs. Transcripts with the complete coding region of important detoxification enzymes were further classified. There were 57 transcripts of P450s distributed in CYP2, CYP3, CYP4, and Mito CYP clades, 29 transcripts of ESTs from 4 functional groups, 17 transcripts of GSTs classified into 5 families, 51 transcripts of ABCs distributed in 6 families, and 19 transcripts of UGTs. Twenty-five lncRNAs were predicted to be regulators of these detoxification genes. Full-length transcriptome sequencing is an efficient method for molecular study of P. striolata and it is also useful for gene function analysis.

FastD: Fast detection of insecticide target-site mutations and overexpressed detoxification genes in insect populations from RNA-Seq data

Target-site mutations and detoxification gene overexpression are two major mechanisms conferring insecticide resistance. Molecular assays applied to detect these resistance genetic markers are time-consuming and with high false-positive rates. RNA-Seq data contains information on the variations within expressed genomic regions and expression of detoxification genes. However, there is no corresponding method to detect resistance markers at present. Here, we collected 66 reported resistance mutations of four insecticide targets (AChE, VGSC, RyR, and nAChR) from 82 insect species. Next, we obtained 403 sequences of the four target genes and 12,665 sequences of three kinds of detoxification genes including P450s, GSTs, and CCEs. Then, we developed a Perl program, FastD, to detect target-site mutations and overexpressed detoxification genes from RNA-Seq data and constructed a web server for FastD (http://www.insect-genome.com/fastd). The estimation of FastD on simulated RNA-Seq data showed high sensitivity and specificity. We applied FastD to detect resistant markers in 15 populations of six insects, Plutella xylostella, Aphis gossypii, Anopheles arabiensis, Musca domestica, Leptinotarsa decemlineata and Apis mellifera. Results showed that 11 RyR mutations in P. xylostella, one nAChR mutation in A. gossypii, one VGSC mutation in A. arabiensis and five VGSC mutations in M. domestica were found to be with frequency difference >40% between resistant and susceptible populations including previously confirmed mutations G4946E in RyR, R81T in nAChR and L1014F in VGSC. And 49 detoxification genes were found to be overexpressed in resistant populations compared with susceptible populations including previously confirmed detoxification genes CYP6BG1, CYP6CY22, CYP6CY13, CYP6P3, CYP6M2, CYP6P4 and CYP4G16. The candidate target-site mutations and detoxification genes were worth further validation. Resistance estimates according to confirmed markers were consistent with population phenotypes, confirming the reliability of this program in predicting population resistance at omics-level.

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.

Analysis of dynamic and widespread lncRNA and miRNA expression in fetal sheep skeletal muscle

The sheep is an economically important animal, and there is currently a major focus on improving its meat quality through breeding. There are variations in the growth regulation mechanisms of different sheep breeds, making fundamental research on skeletal muscle growth essential in understanding the regulation of (thus far) unknown genes. Skeletal muscle development is a complex biological process regulated by numerous genes and non-coding RNAs, including microRNAs (miRNAs) and long non-coding RNAs (lncRNAs). In this study, we used deep sequencing data from sheep longissimus dorsi (LD) muscles sampled at day 60, 90, and 120 of gestation, as well as at day 0 and 360 following birth, to identify and examine the lncRNA and miRNA temporal expression profiles that regulate sheep skeletal myogenesis. We stained LD muscles using histological sections to analyse the area and circumference of muscle fibers from the embryonic to postnatal development stages. Our results showed that embryonic skeletal muscle growth can be characterized by time. We obtained a total of 694 different lncRNAs and compared the differential expression between the E60 vs. E90, E90 vs. E120, E120 vs. D0, and D0 vs. D360 lncRNA and gene samples. Of the total 701 known sheep miRNAs we detected, the following showed a wide range of expression during the embryonic stage: miR-2387, miR-105, miR-767, miR-432, and miR-433. We propose that the detected lncRNA expression was time-specific during the gestational and postnatal stages. GO and KEGG analyses of the genes targeted by different miRNAs and lncRNAs revealed that these significantly enriched processes and pathways were consistent with skeletal muscle development over time across all sampled stages. We found four visual lncRNA–gene regulatory networks that can be used to explore the function of lncRNAs in sheep and may be valuable in helping improve muscle growth. This study also describes the function of several lncRNAs that interact with miRNAs to regulate myogenic differentiation.

Multiple Known Mechanisms and a Possible Role of an Enhanced Immune System in Bt-Resistance in a Field Population of the Bollworm, Helicoverpa zea: Differences in Gene Expression with RNAseq

Several different agricultural insect pests have developed field resistance to Bt (Bacillus thuringiensis) proteins (ex. Cry1Ac, Cry1F, etc.) expressed in crops, including corn and cotton. In the bollworm, Helicoverpa zea, resistance levels are increasing; recent reports in 2019 show up to 1000-fold levels of resistance to Cry1Ac, a major insecticidal protein in Bt-crops. A common method to analyze global differences in gene expression is RNA-seq. This technique was used to measure differences in global gene expression between a Bt-susceptible and Bt-resistant strain of the bollworm, where the differences in susceptibility to Cry1Ac insecticidal proteins were 100-fold. We found expected gene expression differences based on our current understanding of the Bt mode of action, including increased expression of proteases (trypsins and serine proteases) and reduced expression of Bt-interacting receptors (aminopeptidases and cadherins) in resistant bollworms. We also found additional expression differences for transcripts that were not previously investigated, i.e., transcripts from three immune pathways-Jak/STAT, Toll, and IMD. Immune pathway receptors (ex. PGRPs) and the IMD pathway demonstrated the highest differences in expression. Our analysis suggested that multiple mechanisms are involved in the development of Bt-resistance, including potentially unrecognized pathways.

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.

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.

lincRNA_Tc13743.2-miR-133-5p-TcGSTm02 regulation pathway mediates cyflumetofen resistance in Tetranychus cinnabarinus

Differential expression of metabolic detoxification enzymes is an important mechanism involved in pesticide/acaricide resistance of mite pests. The competing endogenous RNA hypothesis offers a new opportunity to investigate post-transcriptional regulation of those genes. In this study, 4454 long non-coding RNAs were identified in the carmine spider mite Tetranychus cinnabarinus by transcriptome sequencing. Software-based predictions indicated that a long intergenic non-coding RNA, (lincRNA)_Tc13743.2 and a detoxification enzyme gene, TcGSTm02, both contained a microRNA (miR-133-5p) response element. Over-expression of lincRNA_Tc13743.2 and TcGSTm02 were detected in a cyflumetofen-resistant T. cinnabarinus strain (CyR), whereas down-regulation of miR-133-5p was observed in this strain. Conversely, up-regulation of miR-133-5p could inhibit TcGSTm02 expression levels, and both lincRNA_Tc13743.2 and TcGSTm02 were significantly enriched in miR-133-5p biotin-avidin pull-down assays. RNA-binding protein immunoprecipitation assay showed that lincRNA_Tc13743.2 and TcGSTm02 bound to a silencing complex containing miR-133-5p. Moreover, a luciferase reporter assay based on a human cell line revealed that over-expression of lincRNA_Tc13743.2 could significantly reduce the inhibition exerted by miR-133-5p through the TcGSTm02 3'UTR. In addition, co-localization of lincRNA_Tc13743.2 and miR-133-5p was detected using fluorescence in situ hybridization, suggesting that lincRNA_Tc13743.2 interacts directly with miR-133-5p in spider mites. More importantly, silencing the expression of lincRNA_Tc13743.2 significantly reduced the expression levels of TcGSTm02 and increased the sensitivity of spider mites to cyflumetofen. Our data show that lincRNA_Tc13743.2 up-regulates TcGSTm02 expression by competing for miR-133-5p binding, demonstrating that a lincRNA_Tc13743.2-miR-133-5p-TcGSTm02 pathway mediates cyflumetofen resistance in mites.

Suppression of Transferrin Expression Enhances the Susceptibility of Plutella xylostella to Isaria cicadae

Transferrins (Trfs) are multifunctional proteins with key functions in iron transport. In the present study, a Trf (PxTrf) from Plutella xylostella was identified and characterized. The PxTrf consisted of a 2046-bp open reading frame, which encoded a 681 amino acid protein with a molecular weight of 73.43 kDa and had an isoelectric point of 7.18. Only a single iron domain was predicted in the N-lobe of PxTrf. Although PxTrf was expressed ubiquitously, the highest levels of expression were observed in the fourth instar larvae. PxTrf transcript level was highest in fat bodies among various tissues. The PxTrf transcript levels increased significantly after the stimulation of pathogens. A decrease in PxTrf expression via RNA interference enhanced the susceptibility of P. xylostella to the Isaria cicadae fungus and inhibited hemocyte nodulation in response to the fungal challenge. In addition, a considerable increase in the pupation rate was observed in larvae treated with double-stranded PxTrf (dsPxTrf). Overall, according to the results, PxTrf may participate in P. xylostella immunity against fungal infection and insect development.

Post-transcriptional regulation of insect metamorphosis and oogenesis

Metamorphic transformation from larvae to adults along with the high fecundity is key to insect success. Insect metamorphosis and reproduction are governed by two critical endocrines, juvenile hormone (JH), and 20-hydroxyecdysone (20E). Recent studies have established a crucial role of microRNA (miRNA) in insect metamorphosis and oogenesis. While miRNAs target genes involved in JH and 20E-signaling pathways, these two hormones reciprocally regulate miRNA expression, forming regulatory loops of miRNA with JH and 20E-signaling cascades. Insect metamorphosis and oogenesis rely on the coordination of hormones, cognate genes, and miRNAs for precise regulation. In addition, the alternative splicing of genes in JH and 20E-signaling pathways has distinct functions in insect metamorphosis and oogenesis. We, therefore, focus in this review on recent advances in post-transcriptional regulation, with the emphasis on the regulatory role of miRNA and alternative splicing, in insect metamorphosis and oogenesis. We will highlight important new findings of miRNA interactions with hormonal signaling and alternative splicing of JH receptor heterodimer gene Taiman.

Identification of Long Noncoding RNA Associated ceRNA Networks in Rosacea

Rosacea is a chronic and relapsing inflammatory cutaneous disorder with highly variable prevalence worldwide that adversely affects the health of patients and their quality of life. However, the molecular characterization of each rosacea subtype is still unclear. Furthermore, little is known about the role of long noncoding RNAs (lncRNAs) in the pathogenesis or regulatory processes of this disorder. In the current study, we established lncRNA-mRNA coexpression networks for three rosacea subtypes (erythematotelangiectatic, papulopustular, and phymatous) and performed their functional enrichment analyses using Gene Onotology, KEGG, GSEA, and WGCNA. Compared to the control group, 13 differentially expressed lncRNAs and 525 differentially expressed mRNAs were identified in the three rosacea subtypes. The differentially expressed genes identified were enriched in four signaling pathways and the GO terms found were associated with leukocyte migration. In addition, we found nine differentially expressed lncRNAs in all three rosacea subtype-related networks, including NEAT1 and HOTAIR, which may play important roles in the pathology of rosacea. Our study provided novel insights into lncRNA-mRNA coexpression networks to discover the molecular mechanisms involved in rosacea development that can be used as future targets of rosacea diagnosis, prevention, and treatment.

Long-noncoding RNAs (lncRNAs) in drug metabolism and disposition, implications in cancer chemo-resistance

Drug metabolism is an orchestrated process in which drugs are metabolized and disposed through a series of specialized enzymes and transporters. Alterations in the expression and/or activity of these enzymes and transporters can affect the bioavailability (pharmacokinetics, or PK) and therapeutic efficacy (pharmacodynamics, or PD) of drugs. Recent studies have suggested that the long non-coding RNAs (lncRNAs) are highly relevant to drug metabolism and drug resistance, including chemo-resistance in cancers, through the regulation of drug metabolism and disposition related genes. This review summarizes the regulation of enzymes, transporters, or regulatory proteins involved in drug metabolism by lncRNAs, with a particular emphasis on drug metabolism and chemo-resistance in cancer patients. The perspective strategies to integrate multi-dimensional pharmacogenomics data for future in-depth analysis of drug metabolism related lncRNAs are also proposed. Understanding the role of lncRNAs in drug metabolism will not only facilitate the identification of novel regulatory mechanisms, but also enable the discovery of lncRNA-based biomarkers and drug targets to personalize and improve the therapeutic outcome of patients, including cancer patients.

Comparative analysis of dsRNA-induced lncRNAs in three kinds of insect species

Long noncoding RNAs (lncRNAs) that have immune responses to various stimuli have been identified in some insects. One type of pathogen-associated molecular pattern, double-stranded RNA (dsRNA), can trigger the RNA interference (RNAi) pathway and immune response. Interestingly, there has been no research into characterizing the relationship between lncRNA and dsRNA-induced RNAi pathways. In this study, dsRNA-induced lncRNAs were investigated in two species of lepidopteran insects, Helicoverpa armigera and Plutella xylostella, and one species of coleopteran insects, Tribolium castaneum. Between untreated group and dsRNA-induced group; 3,463 H. armigera, 6,245 P. xylostella, and 3,067 T. castaneum differentially expressed lncRNAs were identified while 156 H. armigera, 247 P. xylostella, 415 T. castaneum lncRNAs and their putative target genes showed consistent changes in gene expression. In T. castaneum, most target genes of the differentially expressed lncRNAs are enriched in the cyclic adenosine monophosphate signaling pathway, ABC transporters, and Janus kinase-signal transducers and activators of the transcription signaling pathway. Conversely, in H. armigera and P. xylostella, the differentially expressed lncRNAs were mainly enriched in the metabolic, digestive, and synthetic signaling pathways. This result indicates that dsRNA-induced lncRNA is species-dependent. We also found that both Dicer-2 and the lncRNA that targets Dicer-2 were significantly upregulated after dsRNA treatment in P. xylostella, indicating that some lncRNAs may be involved in the regulation of the core RNAi pathway in insects. Our results are the first to identify a relationship between lncRNAs and dsRNA in various insect species with different RNAi efficiencies. These results provide a reference for future study of the dsRNA-induced RNAi pathway and different RNAi efficiencies among insect species.

Deregulated lncRNA expression profile in the mouse lung adenocarcinomas with KRAS-G12D mutation and P53 knockout

Recent studies have demonstrated that aberrant long non‐coding RNAs (lncRNAs) expression are suggested to be closely associated with multiple human diseases, lung cancer included. However, the roles of lncRNAs in lung cancer are not well understood. In this study, we used microarrays to investigate the aberrantly expressed lncRNAs in the mouse lung adenocarcinoma with P53 knockout and the Kras^G12D mutation. Results revealed that 6424 lncRNAs were differentially expressed (≥ 2‐fold change, P < .05). Two hundred and ten lncRNAs showed more than 8‐fold change and conserved across human and were further analysed in the primary mouse lung adenocarcinoma KP cells, which were isolated from the p53 knockout and the Kras^G12D mutation mice. Among all the 210 lncRNAs, 11 lncRNAs' expression was regulated by P53, 33 lncRNAs by KRAS and 13 lncRNAs by hypoxia in the primary KP cells, respectively. NONMMUT015812, which was remarkably up‐regulated in the mouse lung adenocarcinoma and negatively regulated by the P53 re‐expression, was detected to analyse its cellular function. Results showed that knockdown of NONMMUT015812 by shRNAs decreased proliferation and migration abilities of KP cells. Among those aberrantly expressed lncRNAs in the mouse lung adenocarcinoma, NONMMUT015812 was a potential oncogene.

A long non-coding RNA regulates cadherin transcription and susceptibility to Bt toxin Cry1Ac in pink bollworm, Pectinophora gossypiella

Extensive planting of transgenic crops producing insecticidal proteins from the bacterium Bacillus thuringiensis (Bt) has spurred increasingly rapid evolution of resistance in pests. In the pink bollworm, Pectinophora gossypiella, a devastating global pest, resistance to Bt toxin Cry1Ac produced by transgenic cotton is linked with mutations in a gene (PgCad1) encoding a cadherin protein that binds Cry1Ac in the larval midgut. We previously reported a long non-coding RNA (lncRNA) in intron 20 of cadherin alleles associated with both resistance and susceptibility to Cry1Ac. Here we tested the hypothesis that reducing expression of this lncRNA decreases transcription of PgCad1 and susceptibility to Cry1Ac. Quantitative RT-PCR showed that feeding susceptible neonates small interfering RNAs (siRNAs) targeting this lncRNA but not PgCad1 decreased the abundance of transcripts of both the lncRNA and PgCad1. Moreover, neonates fed the siRNAs had lower susceptibility to Cry1Ac. The results imply that the lncRNA increases transcription of PgCad1 and susceptibility of pink bollworm to Cry1Ac. The results suggest that disruption of lncRNA expression could be a novel mechanism of pest resistance to Bt toxins.

Identification and RNAi-based function analysis of chitinase family genes in diamondback moth, Plutella xylostella

BACKGROUND

Insect chitinases play a vital part in chitin degradation in exoskeletons and gut linings during the molting process, and therefore are considered potential targets for new insecticide designs or RNA interference (RNAi)-based pest management. Systematic functional analysis of chitinase genes has already been conducted in several insect pests, but not Plutella xylostella.

RESULTS

In this study, 13 full-length chitinase transcripts were obtained in P. xylostella. Developmental and tissue-specific expression pattern analysis revealed that seven chitinase transcripts were periodically expressed during molting stage and mainly expressed in the integument or midgut, including PxCht3, PxCht5, PxCht6-2, PxCht7, PxCht8, PxCht10 and PxCht-h. RNAi-mediated knockdown of these specific expressed genes revealed that PxCht5 and PxCht10 were essential in larval molting, pupation and eclosion, and PxCht7 was indispensable only in eclosion. No significant effects were observed on insect survival or normal development when the rest chitinase transcripts were suppressed by RNAi.

CONCLUSION

Our results indicated the function of P. xylostella chitinase family genes during the molting process, and may provide potential targets for RNAi-based management of P. xylostella. © 2018 Society of Chemical Industry.

Long non‑coding RNA CASC2 enhances berberine‑induced cytotoxicity in colorectal cancer cells by silencing BCL2

Berberine, a natural isoquinoline alkaloid derived from Berberis species, has been reported to have anticancer effects. However, the mechanisms of action in human colorectal cancer (CRC) are not well established to date. In the present study, the cell cytotoxicity effect of berberine on human CRC cells, as well as the possible mechanisms involved, was investigated. The results of the cell viability and apoptosis assay revealed that treatment of CRC cells with berberine resulted in inhibition of cell viability and activation of cell apoptosis in a concentration-dependent manner. To reveal the underlying mechanism of berberine-induced anti-tumor activity and cell apoptosis, RNA-sequencing followed by reverse-transcription quantitative PCR were performed. In addition, RNA immunoprecipitation, chromatin immunoprecipitation and western blot analysis were used to identify the functional regulation of CASC2/EZH2/BCL2 axis in berberine-induced CRC cell apoptosis. The data revealed that lncRNA CASC2 was upregulated by berberine treatment. Gain- or loss-of-function assays suggested that lncRNA CASC2 was required for the berberine-induced inhibition of cell viability and activation of cell apoptosis. Subsequently, the downstream antiapoptotic gene BCL2 was identified as a functional target of the berberine/CASC2 mechanism, as BCL2 reversed the berberine/CASC2-induced cell cytotoxicity. lncRNA CASC2 silenced BCL2 expression by binding to the promoter region of BCL2 in an EZH2-dependent manner. In summary, berberine may be a novel therapeutic agent for CRC and lncRNA CASC2 may serve as an important therapeutic target to improve the anticancer effect of berberine.

The emerging importance of noncoding RNAs in the insecticide tolerance, with special emphasis on Plutella xylostella (Lepidoptera: Plutellidae)

Recently generated high-throughput sequencing data sets have shed light on the important regulatory roles of noncoding RNA (ncRNA) molecules in the development of higher organisms. Nowadays it is well-known that regulatory ncRNAs can bind complementary RNA or DNA sequences and recruit chromatin remodelers to selectively modulate gene expression. Consequently, genome sequencing and transcriptomics technologies are now being used to reveal hidden associations among ncRNAs and distinct biological mechanisms. This is the case for the diamondback moth Plutella xylostella, a worldwide pest known to infest cruciferous crops and to display resistance to most insecticides, including Bacillus thuringiensis (Bt) based biopesticides. In P. xylostella, it is thought that ncRNAs could play important roles in both development and insecticide resistance. This review will highlight recent insights into the roles of ncRNAs in P. xylostella and related lepidopterans, and will outline genetic engineering technologies which might be used to design efficient ncRNA-based pest control strategies. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs.

Differential lncRNA expression profiles reveal the potential roles of lncRNAs in antiviral immune response of Crassostrea gigas

Long noncoding RNAs (lncRNAs) may play widespread roles in various biological processes. However, systematic profiles of lncRNAs in the biological responses of Pacific Oyster (Crassostrea gigas) to pathogen infection have not yet been demonstrated. Here, we have conducted an exhaustive comparative transcriptome analysis using a bioinformatics approach to exam the functions of lncRNAs response to Ostreid herpesvirus 1μVar (OsHV-1μVar) challenge. In total, 101 differentially expressed lncRNAs (DE-lncRNA) during OsHV-1μVar infections were identified. Compared with differentially expressed mRNAs (DE-mRNA), DE-lncRNAs are shorter in terms of overall length but longer in terms of exon length. These lncRNAs shared similar characteristics with previously reported invertebrate lncRNAs, such as relatively low GC content, low exon number and low sequence conservation, but low expression level were not observed. 20 DE-lncRNAs are typically co-expressed with their neighboring genes annotated as GO terms (GO: 0044237), indicating that these lncRNAs are involved in binding and cellular process functions in cis mode. The weighted gene co-expression network (WGCNA) analysis resulted in 15 modules. The highlighted blue module was specifically demonstrated a co-expression relationship between 14 DE-lncRNAs and 17 immune-related DE-mRNAs (IR-DE-mRNA). Three hub lncRNAs within this module were co-expressed with one hub IR-DE-mRNA involved in fibrinogen-related protein. It was speculated that lncRNAs is extensively involved in oyster antiviral innate immune system. The present study will facilitate subsequently experimental studies to unravel the function of lncRNAs in marine invertebrate response to pathogen infection.

Altered Expression of Differential Genes in Thoracic Spinal Cord Involved in Experimental Cholestatic Itch Mouse Model

The spinal origin of cholestatic itch in experimental obstructive jaundice mouse model remains poorly understood. In this study, the jaundice model was established by bile duct ligation (BDL) in mice, and differential gene expression patterns were analyzed in the lower thoracic spinal cord involved in cholestatic pruritus after BDL operation using high-throughput RNA sequencing. At 21st day after BDL, the expression levels of ENSRNOG00000060523, ENSRNOG00000058405 and ENSRNOG00000055193 mRNA were significantly up-regulated, and those of ENSRNOG00000042197, ENSRNOG00000008478, ENSRNOGOOOOOO19607, ENSRNOG00000020647, ENSRNOG00000046289, Gemin8, Serpina3n and Trim63 mRNA were significantly down-regulated in BDL group. The RNAseq data of selected mRNAs were validated by RT-qPCR. The expression levels of ENSRNOG00000042197, ENSRNOG00000008478, ENSRNOGOOOOOO 19607, ENSRNOG00000020647, ENSRNOG00000046289 and Serpina3n mRNA were significantly down-regulated in BDL group. This study suggested that cholestatic pruritus in experimental obstructive jaundice mouse model is related with in the changes of gene expression profiles in spinal cord.

Genome-Wide Identification and Functional Prediction of Novel Drought-Responsive lncRNAs in Pyrus betulifolia

Increasing evidence shows that long noncoding RNAs (lncRNAs) play important roles in developmental regulation and many other biological processes in plants. However, identification of lncRNAs in Pyrus betulifolia is limited compared with studies of functional gene expression. Using high-throughput sequencing technology, the transcriptome of P. betulifolia under drought stress was analyzed to identify lncRNAs. A total of 14,478 lncRNAs were identified, of which 251 were found to be drought-responsive. The putative target genes of these differentially expressed lncRNAs were significantly enriched in metabolic processes, organic substance metabolic processes, macromolecule metabolic processes, and heterocyclic compound binding. Real-time quantitative polymerase chain reaction validation suggested that the results of the RNA sequencing data analysis were reliable. This study will provide genetic resources for pear breeding and provide reference to other pomological studies.

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.

Altered expression of differential gene and lncRNA in the lower thoracic spinal cord on different time courses of experimental obstructive jaundice model accompanied with altered peripheral nociception in rats

The spinal origin of jaundice-induced altered peripheral nociceptive response poorly understood. In the current study, we aimed to first validate rats with bile duct ligation (BDL) as a jaundice model accompanied by altered peripheral nociceptive response, and then to analyze differential gene and lncRNA expression patterns in the lower thoracic spinal cord on different time courses after BDL operation by using high-throughput RNA sequencing. The differentially expressed genes (DEGs) identified using reverse transcription-quantitative polymerase chain reaction (RT-qPCR) analysis, followed by clustering analysis, Gene Ontology analysis and pathway analysis. As a result, a total of 2033 lncRNAs were differentially expressed 28d after BDL, in which 1545 probe sets were up-regulated and 488 probe sets were down-regulated, whereas a total of 2800 mRNAs were differentially expressed, in which 1548 probe sets were up-regulated and 1252 probe sets were down-regulated. The RNAseq data of select mRNAs and lncRNAs was validated by RT-qPCR. 28d after BDL, the expressions of lncRNA NONRATT002335 and NONRATT018085 were significantly up-regulated whereas the expression of lncRNA NONRATT025415, NONRATT025388 and NONRATT025409 was significantly down-regulated. 14d after BDL, the expressions of lncRNA NONRATT002335 and NONRATT018085 were significantly up-regulated; the expression of lncRNA NONRATT025415, NONRATT025388 and NONRATT025409 was significantly down-regulated. In conclusion, the present study showed that jaundice accompanied with decreased peripheral nociception involved in the changes of gene and lncRNA expression profiles in spinal cord. These findings extend current understanding of spinal mechanism for obstructive jaundice accompanied by decreased peripheral nociception.