Next Article in Journal
Mussel-Inspired Fabrication of Konjac Glucomannan/Poly (Lactic Acid) Cryogels with Enhanced Thermal and Mechanical Properties
Next Article in Special Issue
New Method for Differentiation of Granuloviruses (Betabaculoviruses) Based on Multitemperature Single Stranded Conformational Polymorphism
Previous Article in Journal
IFN-β: A Contentious Player in Host–Pathogen Interaction in Tuberculosis
Previous Article in Special Issue
Identification and Characterization of Hyphantria cunea Aminopeptidase N as a Binding Protein of Bacillus thuringiensis Cry1Ab35 Toxin
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 18
Issue 12
10.3390/ijms18122728
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Differentially Expressed miRNAs in Colorado Potato Beetles (Leptinotarsa decemlineata (Say)) Exposed to Imidacloprid

by
Mathieu D. Morin
1,
Pierre J. Lyons
2,
Nicolas Crapoulet
2,
Sébastien Boquel
3 and
Pier Jr Morin
1,*
1
Department of Chemistry and Biochemistry, University of Moncton, 18 Antonine-Maillet Avenue, Moncton, NB E1A 3E9, Canada
2
Atlantic Cancer Research Institute, Pavillon Hôtel-Dieu 35 Providence Street, Moncton, NB E1C 8X3, Canada
3
Fredericton Research and Development Centre, Agriculture and Agri-Food Canada, 850 Lincoln Road, Fredericton, NB E3B 4Z7, Canada
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2017, 18(12), 2728; https://doi.org/10.3390/ijms18122728
Submission received: 3 November 2017 / Revised: 8 December 2017 / Accepted: 12 December 2017 / Published: 16 December 2017
(This article belongs to the Special Issue Molecular Entomology of Insects of Economic Importance)
Browse Figure
Versions Notes

Abstract

:
The Colorado potato beetle (Leptinotarsa decemlineata (Say)) is a significant pest of potato plants that has been controlled for more than two decades by neonicotinoid imidacloprid. L. decemlineata can develop resistance to this agent even though the molecular mechanisms underlying this resistance are not well characterized. MicroRNAs (miRNAs) are short ribonucleic acids that have been linked to response to various insecticides in several insect models. Unfortunately, the information is lacking regarding differentially expressed miRNAs following imidacloprid treatment in L. decemlineata. In this study, next-generation sequencing and quantitative real-time polymerase chain reaction (qRT-PCR) were used to identify modulated miRNAs in imidacloprid-treated versus untreated L. decemlineata. This approach identified 33 differentially expressed miRNAs between the two experimental conditions. Of interest, miR-282 and miR-989, miRNAs previously shown to be modulated by imidacloprid in other insects, and miR-100, a miRNA associated with regulation of cytochrome P450 expression, were significantly modulated in imidacloprid-treated beetles. Overall, this work presents the first report of a miRNA signature associated with imidacloprid exposure in L. decemlineata using a high-throughput approach. It also reveals interesting miRNA candidates that potentially underly imidacloprid response in this insect pest.

1. Introduction

The Colorado potato beetle (CPB) (Leptinotarsa decemlineata (Say)) is a significant insect pest harming potato crops worldwide [1]. CPBs are often considered to be the primary insect associated with potato plant defoliation [2], leading to up to 75% of foliage consumption and impacting revenues for growers [3,4]. Pest control strategies targeting CPBs typically involve pesticides, even though resistance against a variety of such compounds, including spinosad, thiamethoxam, and deltamethrin, has been observed [5,6,7]. Neonicotinoids, which target the nicotinic acetylcholine receptors (nAChRs), are a class of insecticides that are used extensively against CPBs. The neonicotinoid imidacloprid (1-(6-chloro-3-pyridylmethyl)-N-nitroimidazolidin-2-ylideneamine) was registered for such a purpose over 20 years ago, and was initially successful in managing CPB populations [6]. Unfortunately, early studies reported adult CPBs exhibiting as much as a 155-fold increase in imidacloprid resistance in select populations [8]. Efforts have been deployed in recent years to understand the molecular players associated with imidacloprid resistance in CPBs. Multiple studies have notably highlighted the differential expression of cytochrome P450s in imidacloprid-treated CPBs [9,10], positioning these enzymes as targets that could sensitize CBPs to this agent. Nevertheless, the complete molecular picture associated with imidacloprid resistance in CPBs remains incompletely characterized.
MicroRNAs (miRNAs) are short conserved non-coding RNAs capable of regulating the expression of multiple target mRNA transcripts. These molecules are involved in the response to various stresses including starvation, anoxia, and freezing [11,12]. Previous work has also revealed miRNA modulation in insects following exposure to different chemicals including pyrethroids, chlorantraniliprole and fenpropathrin [13,14,15]. In addition, recent studies have highlighted the likely regulation of pyrethroid resistance in the mosquito Culex pipiens pallens (L.) by miRNAs [16,17]. These results thus support the potential involvement of miRNAs in insecticide resistance. Unfortunately, data is currently lacking regarding a signature of differentially expressed miRNAs in response to imidacloprid treatment in CPBs.
The current study was undertaken to characterize the miRNA footprint observed in CPBs following imidacloprid exposure, and to identify novel miRNA-associated molecular leads that could be utilized in the development of alternative strategies to control CPBs using an approach relying on next-generation sequencing and qRT-PCR. The miRNA target prediction was also performed using predictive tools to better assess the likely consequences of these imidacloprid-associated miRNAs. Transcript targets linked to transcriptional regulation and glucose metabolism were observed and are further discussed.

2. Results

2.1. Small RNA Sequence Analysis

High-throughput sequencing resulted in 59,406,387 reads among all samples. Low quality reads were discarded. Reads ranging from 16 to 60 nucleotides were conserved for analysis. Small RNA libraries displayed comparable read distribution profiles in the 16 and 24 nucleotides region, with maximum peaks observed at 21 nucleotides. This approach resulted in 22,534,030 and 13,763,160 total reads for control and imidacloprid-treated insects, respectively (Table 1). Within these reads, 8,825,610 (control) and 6,211,853 (imidacloprid) unique reads were obtained. MiRBase was utilized to identify the small RNAs of interest in L. decemlineata samples using red flour beetle Tribolium castaneum reference databases [18]. The miRNAs were annotated using sRNAbench in sRNAtoolbox [19]. Small RNAs unique reads were mapped to miRNAs (26,447 for control and 24,221 for imidacloprid-treated). SnRNAs (4735 for control and 3616 for imidacloprid-treated), snRNAs (4705 for control and 2937 for imidacloprid-treated) and tRNAs (40,144 for control and 22,115 for imidacloprid-treated) were also observed.

2.2. miRNA Expression in Control and Imidacloprid-Treated L. decemlineata by High-Throughput Sequencing

The most frequent miRNAs detected with this approach were miR-14-3p (95,802 reads in control and 81,633 reads in imidacloprid-exposed insects) and miR-8-3p (68,314 reads in control and 79,134 reads in imidacloprid-exposed insects) (Table 2). Based on next-generation sequencing results, log2 ratios of miRNA expression in control and imidacloprid-exposed L. decemlineata demonstrated 33 differentially expressed miRNAs. A total of 14 upregulated and 19 downregulated miRNAs displaying absolute log2 fold-changes greater than 0.3 (p < 0.05, n = 3) were identified (Table 3).

2.3. qRT-PCR Quantification of Selected miRNAs

Eight miRNAs were further amplified and quantified by qRT-PCR. These miRNAs displayed reduced, elevated or stable levels in insects treated with imidacloprid for 8 h when compared with control insects using deep sequencing. Fold-changes were compared with expression results measured by next-generation sequencing (Figure 1). MiR-92a-5p, miR-133-3p, miR-305-5p, miR-927a-5p and miR-989-3p displayed substantial modulation in imidacloprid-treated versus untreated insects as measured by high-throughput sequencing, with 0.59-fold, 1.23-fold, 0.77-fold, 0.74-fold and 0.57-fold variations, respectively. These changes were comparable to the ones obtained by qRT-PCR. The miR-12-5p levels were also comparable when both methods were used to quantify its expression. MiR-7-5p and miR-9a-5p, stable miRNAs in imidacloprid-treated insects as measured by next-generation sequencing, displayed reduced levels (0.62-fold and 0.73-fold, respectively) when measured by qRT-PCR. Assessment of miRNA expression in insects treated with imidacloprid for 24 h was also performed by qRT-PCR (Figure 1). Differential expression of miR-7-5p, miR-9a-5p, miR-12-5p, miR-92a-5p, and miR-133-5p was observed between insects submitted to imidacloprid treatments of varying lengths. Levels of miR-305-5p, miR-927a-5p and miR-989-3p remained stable between the two treatment groups.

2.4. Prediction of miRNA Targets and Functional Classification

Target prediction tools were able to identify transcripts potentially regulated by 17 L. decemlineata imidacloprid-associated miRNAs. Targets were assessed for the following miRNAs: miR-282-5p, miR-1000-5p, miR-193-3p, miR-124-3p, miR-970-3p, let-7-5p, miR-1-3p, miR-133-3p, miR-263b-5p, miR-92a-3p, miR-305-5p, miR-927a-5p, miR-9c-5p, miR-316-5p, miR-100-5p, miR-315-5p and miR-989-3p (Table 4). Common targets are underlined. Biological processes with predicted transcripts are presented in Table 5. Several target transcripts were associated with transcriptional regulation, whilst others were linked to glucose metabolism.

3. Discussion

The molecular changes underlying neonicotinoid resistance in insects have been investigated with great interest in recent years. Work conducted on various insect models treated with imidacloprid has revealed targets potentially involved in the response to this agent. These targets include select cytochrome P450s, ATP-binding cassette (ABC) transporters, and nAChRs [20,21,22]. Numerous examples exist that highlight the regulation of these imidacloprid-relevant targets by miRNAs [23,24]. While recent studies have started to associate miRNAs with insecticide resistance in various insects [14,25,26], no work so far has attempted to characterize modulated miRNAs in imidacloprid-treated CPBs. Using a high-throughput sequencing approach, the current study reveals a signature of differentially expressed miRNAs in CPBs exposed to this neonicotinoid.
The information is sparse regarding the molecular changes associated with imidacloprid response in insects. Pioneering work performed on beehives submitted to low levels of imidacloprid revealed a set of differentially expressed miRNAs in larvae gathered from imidacloprid-exposed versus unexposed hives [27]. The present study also identified similar changes in expression levels of select miRNAs, including miR-282 upregulation and miR-989 downregulation, in imidacloprid-treated CPBs. MiR-989 has been linked with key physiological processes, as Drosophila melanogaster (Meigen) miR-989 mutants notably exhibit impaired border cell migration [28], and Bombyx mori (L.) strongly upregulates miR-989 during the pupal phases [29]. It is interesting to note that levels of miR-989-3p, like miR-305-5p and miR-927a-5p, remained unchanged in insects treated with imidacloprid for different lengths of time, strengthening its potential involvement in response to this insecticide. MiR-282, on the other hand, can influence D. melanogaster viability and longevity by targeting the adenylate cyclase rutabaga [30]. While modulation of miR-282 and miR-989 in imidacloprid-treated CPBs presented here supports their probable role in response to this pesticide, additional work is required to fully delineate their functions in insects exposed to imidacloprid.
Several reports have highlighted the upregulation of cytochrome P450s in CPBs treated with imidacloprid. A high-throughput sequencing approach performed in imidacloprid-sensitive and imidacloprid-resistant CPBs revealed that transcripts for 41 cytochrome P450s were more than 2-fold upregulated in pesticide-exposed insects. In addition, cytochrome P450 activity was increased in the resistant beetles [10]. Gene expression analysis performed on CPB populations susceptible or resistant to imidacloprid reported differential cytochrome P450 levels. This further strengthens an underlying role for these enzymes in imidacloprid resistance [9]. Several studies have reported miRNA-mediated regulation of cytochrome P450 expression. Recent work conducted in the tobacco aphid Myzus persicae nicotianae (Blackman) showed that expression of CYP6CY3, a cytochrome P450 previously associated with neonicotinoid resistance [31], was regulated by miR-100 [32]. The present study highlighted downregulation of miR-100-5p levels in imidacloprid-treated CPBs, supporting a potential role for the miR-100-cytochrome P450 axis in the imidacloprid response.
Functional annotation was undertaken using the predicted targets of differentially expressed miRNAs following imidacloprid exposure in order to better characterize the potential impact of this compound in CPBs (Table 5). This approach notably highlighted gene transcription as a likely process affected by the modulated miRNAs. Some miRNAs, such as miR-92a-3p and miR-927a-5p, were predicted to impact key players involved in transcriptional regulation, including the cAMP response element binding (Creb) transcription factor and the ecdysone receptor (EcR). Previous work has reported that Creb levels were differentially expressed in honey bees (Apis mellifera carnica (Pollmann)) exposed to various neonicotinoids [33], supporting a potential role for Creb in the response to imidacloprid. Several studies have investigated the effect of imidacloprid treatment on gene expression in multiple insect models. A microarray-based approach performed in the mosquito Aedes aegypti (L.) revealed 344 and 108 differentially expressed genes in imidacloprid-resistant larvae and adults, respectively, when compared to their imidacloprid-sensitive counterparts [34]. In addition, a high-throughput sequencing approach in D. melanogaster that were sensitive or resistant to the same agent highlighted 357 transcripts that were either upregulated or downregulated in imidacloprid-resistant flies [35]. Multiple transcript targets associated with imidacloprid-modulated miRNAs were also linked to glucose metabolism. Previous reports have linked imidacloprid treatment with modulation of various targets involved in glucose homeostasis. Honey bee larvae exposed to imidacloprid exhibited increased transcript levels of multiple genes involved in the glycolytic and gluconeogenic pathways. Transcript levels of phosphoenolpyruvate carboxykinase (PEPCK), a rate-limiting enzyme for the latter cascade, were notably upregulated in imidacloprid-exposed larvae [27]. In addition, a study performed on Madagascar hissing cockroaches (Gromphadorhina portentosa (Schaum)) showed variations in glucose uptake and carbohydrate metabolism in cockroaches exposed to imidacloprid [36]. These leads generate potential processes via which the identified miRNAs can influence imidacloprid response. However, further characterization of the link that exists between these processes and imidacloprid-associated miRNAs is envisioned. Ultimately, such transcripts could be targeted by diverse means including RNA interference (RNAi)-based approaches that form part of a strategy aimed at influencing imidacloprid response in CPBs. Such approaches are currently under investigation in various pests including CPBs [37].
In conclusion, the present study revealed a set of modulated miRNAs in CPBs following exposure to the neonicotinoid imidacloprid, and is aligned with recent work aimed at better understanding the functions associated with miRNAs in these insects [38]. These results present the first example of a miRNA signature in L. decemlineata treated with this agent. Assessment of miRNA expression following different imidacloprid treatment durations also revealed varying levels of select miRNAs in short- versus long-term insecticide exposures. This warrants further investigation of early and late molecular changes associated with response to this agent. Predicted miRNA transcript targets further suggest that differential miRNA expression could impact molecular players involved in key processes. It is important to note that while the miRNA changes uncovered in this study generated valuable information regarding miRNAs underlying insect response to imidacloprid, future work is needed to investigate miRNA signatures observed in response to additional experimental conditions. These conditions include repeated topical imidacloprid exposures, usage of a diet comprised of imidacloprid-treated potato leaves, or field CPBs naturally exposed to this agent and exhibiting various degrees of tolerance against imidacloprid. Further investigations will be undertaken to examine the expression status of enzymes such as cytochrome P450s in CPB populations that do or do not respond to imidacloprid. Overall, this work provides novel information on a potential role for miRNAs in imidacloprid response and resistance, as well as contributing to the growing knowledge of the molecular levers underlying the response to this agent in Colorado potato beetles.

4. Materials and Methods

4.1. Insect Collection and Treatment

Adult Colorado potato beetles that had overwintered were collected in nursery potato fields at the Fredericton Research and Development Centre in Fredericton (NB, Canada) in June 2016. Potato fields where collection occurred were not treated with any types of insecticides. This nevertheless cannot rule out the potential exposure of beetles to insecticides in years prior to the study, nor exposure attributable to insect migration. Beetles were placed in plastic containers (50 per) containing potato leaves, and closed using a lid with insect screening for ventilation. Insects were subsequently transported to the Université de Moncton (Moncton, NB, Canada). A group of 30 insects was placed in an incubator (Thermo Fisher Scientific, Waltham, MA, USA) set at 25 °C for 5 days under 16L:8D cycles. Insects were provided with potato plants. Following acclimation, 5 µL (0.5 µg) of analytical grade imidacloprid (100 µg/mL in acetonitrile, Sigma-Aldrich, St. Louis, MO, USA) was applied topically on the abdomen of 15 beetles, as previously described [10]. A volume of 5 µL acetonitrile was applied on 15 beetles in parallel, which were used as controls. Imidacloprid-treated and control insects were returned to the incubator for 8 h. A similar experimental protocol was undertaken where insects were treated with imidacloprid for 24 h. Insects were active after the incubation period and prior to storage. A parallel bioassay performed on insects treated with increasing imidacloprid doses (0.05 µg to 5.0 µg) and not destined for high-throughput sequencing highlighted a marked impact on insect activity, albeit not viability, at maximum doses of 2.5 µg and 5.0 µg of imidacloprid (LD50 >5.0 µg/insect). All insects destined for sequencing were ultimately placed rapidly in liquid nitrogen and stored at −80 °C until RNA isolation.

4.2. Small RNA Isolation

Small RNA fractions were isolated from control and imidacloprid-treated L. decemlineata using the miRVana miRNA Isolation Kit (Thermo Fisher Scientific, Waltham, MA, USA) as per the manufacturer’s instructions and as described previously [39]. Each isolate was prepared using two L. decemlineata. RNA isolates were generated in triplicates. The quality and integrity of small RNA fractions were confirmed using the High Sensitivity RNA ScreenTape on an Agilent TapeStation 2200 instrument (Agilent Technologies (Santa Clara, CA, USA)).

4.3. Small RNA Library Construction and Sequencing

Small RNA libraries were prepared and sequenced following Ion Torrent protocols and as reported previously [40]. Samples were loaded on an Ion PI Chip v2 and sequenced with an Ion Proton Sequencer (Thermo Fisher Scientific, Waltham, MA, USA). Three control and three imidacloprid-treated libraries were constructed. FASTQ files were generated for each sample and adaptor sequences were trimmed. Reads with less than 16 and more than 60 nucleotides were discarded. Reads with Q scores of less than 20 were not included. Small RNA annotation was performed using the sRNAbench tool from sRNAtoolbox. Red flour beetle Tribolium castaneum was used as a reference [19]. The number of permitted mismatches when mapping to mature miRNAs was set to 3. Default values were used for the remaining parameters. Annotated sequences having less than 10 normalized read counts in either control or insecticide-treated samples were discarded, and log2 fold-changes were generated.

4.4. cDNA Synthesis

cDNA was synthetized as described earlier [41]. L. decemlineata miRNA sequences obtained via next-generation sequencing were used to design stem-loop primers to amplify miR-1-3p, miR-7-5p, miR-9a-5p, miR-12-5p, miR-92a-5p, miR-133-3p, miR-305-5p, miR-927a-5p and miR-989-3p (Table 6). A volume containing 0.5–1 µg of small RNA was combined with 5 μL of 300 nM of miRNA-specific stem-loop primer and the mixture was placed at 95 °C for 5 min, 60 °C for 5 min, and on ice for 1 min. Subsequently, 4 μL of 5× first strand buffer, 2 μL of 0.1 M DTT, 1 μL of 10 mM dNTPs, 2 μL of diethyl pyrocarbonate (DEPC)-treated water and 1 μL of M-MLV reverse transcriptase (Thermo Fisher Scientific, Waltham, MA, USA) were added. Mixtures were incubated at 16 °C for 30 min, 42 °C for 30 min and 85 °C for 5 min. cDNA was serially diluted and used for PCR.

4.5. PCR and qRT-PCR Amplification of miRNAs

Initial PCR reactions were performed to confirm correct product amplification. Reactions consisted of 5 µL of cDNA (10−1), 5.5 µL of DEPC-treated water, 1 µL of 25 μM miRNA-specific forward primer, 1 μL of 25 μM universal reverse primer and 12.5 μL of 2× Taq FroggaMix (FroggaBio, Toronto, ON, Canada) [42]. Forward and universal primers used are presented in Table 6. An initial denaturing step at 95 °C for 5 min was performed, followed by 35 cycles at 95 °C for 15 s and at a gradient of temperatures for 1 min. Products were separated on a 2% agarose gel and visualized using a ChemiDoc MP system (Bio-Rad, Hercules, CA, USA). Products were then sequenced at the Université Laval (Quebec City, QC, Canada) and identities were confirmed using Basic Local Alignment Search Tool (BLAST).
The qRT-PCR reactions were performed by mixing 2.5 µL of diluted cDNA template with 0.5 µL of DEPC-treated water, 1 µL of 5 µM miRNA-specific forward primer, 1 µL of 5 µM universal reverse primer and 5 µL of 2× iTaq Universal SYBR Green Supermix (Bio-Rad, Hercules, CA, USA) prepared in triplicate in 96-well plates. Amplification protocol consisted of an initial denaturing step at 95 °C for 3 min, followed by 40 cycles at 95 °C for 15 s and optimal annealing temperature for each miRNAs (Table 6) for 30 s. Reactions were conducted on a CFX Connect Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA). Efficiencies were calculated for each primer pair with the same protocol as above in serial cDNA dilutions. NormFinder software (Aarhus University Hospital, Aarhus, Denmark) [43] revealed miR-1-3p as the most stably expressed miRNA amongst a group of transcripts in insects exposed to imidacloprid for 8 h as determined by qRT-PCR, and was used as reference transcript.

4.6. miRNA Transcript Targets Prediction

TargetScanFly 6.2 [44] and the fruit fly miRanda algorithm [45] target prediction tools were used to identify likely mRNA candidates regulated by the imidacloprid-associated miRNAs identified by high-throughput sequencing. The top five transcript targets linked with imidacloprid-modulated miRNAs that were available in both prediction tools were obtained. A list of miRNA targets was subsequently generated and functionally annotated with the Database for Annotation, Visualization and Integrated Discovery (DAVID) bioinformatics resources v6.8 [46]. Biological processes associated with these targets were generated.

4.7. Quantification and Statistics

The Cq values following qRT-PCR runs were obtained using the Bio-Rad CFX Manager software (Bio-Rad, Hercules, CA, USA). miRNA levels were normalized using miR-1-3p levels amplified from the same sample. The miRNA levels were measured using the 2−ΔΔCq method [47]. Ratios of normalized miRNA levels in control CPBs to average transcript expression in imidacloprid-treated CPBs were generated. Statistical differences between control and treated conditions were assessed with the Student’s t-test. For high-throughput sequencing, annotated sequences were normalized using the trimmed mean of M-values (TMM) function in edgeR Bioconductor [48], and log2 ratios depicting miRNA expression in insecticide-treated versus control CPBs were calculated. Bioconductor “limma” package, with linear model fit and empirical Bayes statistics, was used to calculate differential expression [49]. Results for miRNAs with average normalized expression (ANE) greater than 10, log2 fold-change higher than 0.3 and p < 0.05 were investigated.

Acknowledgments

Pier Jr Morin is supported by a Discovery Grant from the Natural Sciences and Engineering Research Council of Canada (RGPIN/402222-2012), a Research Grant (EARI 15-007) via the Enabling Agricultural and Research Innovation (EARI) program under the Canada/New Brunswick Growing Forward 2 initiative and an Emerging Project Grant (RIF 2016-036) by the New Brunswick Innovation Foundation.

Author Contributions

Mathieu D. Morin, Nicolas Crapoulet, Sébastien Boquel and Pier Jr Morin conceived and designed the experiments; Mathieu D. Morin and Pierre J. Lyons performed the experiments; Mathieu D. Morin, Pierre J. Lyons and Nicolas Crapoulet analyzed the data; Mathieu D. Morin, Sébastien Boquel, and Pier Jr Morin wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Weber, D. Colorado beetle: Pest on the move. Pestic. Outlook 2003, 14, 256–259. [Google Scholar] [CrossRef]
  2. Radcliffe, E.B.; Lagnaoui, A. Insect pests in potato: Insects. In Potato Biology and Biotechnology: Advances and Perspectives; Vreughenhil, D., Bradshaw, J., Gebhardt, C., Govers, F., Taylor, M., MacKerron, D., Ross, H., Eds.; Elsevier: Amsterdam, The Netherlands, 2007; pp. 543–567. [Google Scholar]
  3. Hare, D.J. Impact of defoliation by the Colorado potato beetle on potato yields. J. Econ. Entomol. 1980, 73, 369–373. [Google Scholar] [CrossRef]
  4. Shields, E.J.; Wyman, J.A. Effect of defoliation at specific growth stages on potato yields. J. Econ. Entomol. 1984, 77, 1194–1199. [Google Scholar] [CrossRef]
  5. Jiang, W.H.; Wang, Z.T.; Xiong, M.H.; Lu, W.P.; Liu, P.; Guo, W.C.; Li, G.Q. Insecticide resistance status of Colorado Potato Beetle (Coleoptera: Chrysomelidae) adults in northern Xinjiang Uygur autonomous region. J. Econ. Entomol. 2010, 103, 1365–1371. [Google Scholar] [CrossRef] [PubMed]
  6. Alyokhin, A.; Baker, M.; Mota-Sanchez, D.; Dively, G.; Grafius, E. Colorado Potato Beetle Resistance to Insecticides. Am. J. Potato Res. 2008, 85, 395–413. [Google Scholar] [CrossRef]
  7. Szendrei, Z.; Grafius, E.; Byrne, A.; Ziegler, A. Resistance to neonicotinoid insecticides in field populations of the Colorado potato beetle (Coleoptera: Chrysomelidae). Pest Manag. Sci. 2012, 68, 941–946. [Google Scholar] [CrossRef] [PubMed]
  8. Zhao, J.Z.; Bishop, B.A.; Grafius, E.J. Inheritance and synergism of resistance to imidacloprid in the Colorado potato beetle (Coleoptera: Chrysomelidae). J. Econ. Entomol. 2000, 93, 1508–1514. [Google Scholar] [CrossRef] [PubMed]
  9. Clements, J.; Schoville, S.; Peterson, N.; Lan, Q.; Groves, R.L. Characterizing Molecular Mechanisms of Imidacloprid Resistance in Select Populations of Leptinotarsa decemlineata in the Central Sands Region of Wisconsin. PLoS ONE 2016, 11, e0147844. [Google Scholar] [CrossRef] [PubMed]
  10. Zhu, F.; Moural, T.W.; Nelson, D.R.; Palli, S.R. A specialist herbivore pest adaptation to xenobiotics through up-regulation of multiple Cytochrome P450s. Sci. Rep. 2016, 6, 20421. [Google Scholar] [CrossRef] [PubMed]
  11. Barrio, L.; Dekanty, A.; Milán, M. MicroRNA-mediated regulation of Dp53 in the Drosophila fat body contributes to metabolic adaptation to nutrient deprivation. Cell Rep. 2014, 8, 528–541. [Google Scholar] [CrossRef] [PubMed]
  12. Lyons, P.J.; Storey, K.B.; Morin, P.J. Expression of miRNAs in response to freezing and anoxia stresses in the freeze tolerant fly Eurosta solidaginis. Cryobiology 2015, 71, 97–102. [Google Scholar] [CrossRef] [PubMed]
  13. Hong, S.; Guo, Q.; Wang, W.; Hu, S.; Fang, F.; Lv, Y.; Yu, J.; Zou, F.; Lei, Z.; Ma, K.; et al. Identification of differentially expressed microRNAs in Culex pipiens and their potential roles in pyrethroid resistance. Insect Biochem. Mol. Biol. 2014, 55, 39–50. [Google Scholar] [CrossRef] [PubMed]
  14. Li, X.; Guo, L.; Zhou, X.; Gao, X.; Liang, P. miRNAs regulated overexpression of ryanodine receptor is involved in chlorantraniliprole resistance in Plutella xylostella (L.). Sci. Rep. 2015, 5, 14095. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, Y.; Xu, Z.; Wu, Q.; Peng, M.; Liu, Y.; Liu, X.; Shi, L.; Shen, G.; Pan, Y.; He, L. Identification of Differentially Expressed microRNAs between the Fenpropathrin Resistant and Susceptible Strains in Tetranychus cinnabarinus. PLoS ONE 2016, 11, e0152924. [Google Scholar] [CrossRef] [PubMed]
  16. Liu, B.; Tian, M.; Guo, Q.; Ma, L.; Zhou, D.; Shen, B.; Sun, Y.; Zhu, C. MiR-932 Regulates Pyrethroid Resistance in Culex pipiens pallens (Diptera: Culicidae). J. Med. Entomol. 2016, 53, 1205–1210. [Google Scholar] [CrossRef] [PubMed]
  17. Tian, M.; Liu, B.; Hu, H.; Li, X.; Guo, Q.; Zou, F.; Liu, X.; Hu, M.; Guo, J.; Ma, L.; et al. MiR-285 targets P450 (CYP6N23) to regulate pyrethroid resistance in Culex pipiens pallens. Parasitol. Res. 2016, 115, 4511–4517. [Google Scholar] [CrossRef] [PubMed]
  18. Kozomara, A.; Griffiths-Jones, S. miRBase: Annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014, 42, D68–D73. [Google Scholar] [CrossRef] [PubMed]
  19. Rueda, A.; Barturen, G.; Lebrón, R.; Gómez-Martín, C.; Alganza, Á.; Oliver, J.L.; Hackenberg, M. sRNAtoolbox: An integrated collection of small RNA research tools. Nucleic Acids Res. 2015, 43, 467–473. [Google Scholar] [CrossRef] [PubMed]
  20. Wang, X.; Sun, H.; Zhang, Y.; Liu, C.; Liu, Z. Transcriptional Changes in nAChRs, Interactive Proteins and P450s in Locusta migratoria manilensis (Orthoptera: Acrididae) CNS in Response to High and Low Oral Doses of Imidacloprid. J. Insect Sci. 2015, 15, 102. [Google Scholar] [CrossRef] [PubMed]
  21. Kang, X.L.; Zhang, M.; Wang, K.; Qiao, X.F.; Chen, M.H. Molecular cloning, expression pattern of multidrug resistance associated protein 1 (MRP1, ABCC1) gene, and the synergistic effects of verapamil on toxicity of two insecticides in the bird cherry-oat aphid. Arch. Insect Biochem. Physiol. 2016, 92, 65–84. [Google Scholar] [CrossRef] [PubMed]
  22. Sun, H.; Pu, J.; Chen, F.; Wang, J.; Han, Z. Multiple ATP-binding cassette transporters are involved in insecticide resistance in the small brown planthopper, Laodelphax striatellus. Insect Mol. Biol. 2017, 26, 343–355. [Google Scholar] [CrossRef] [PubMed]
  23. Haenisch, S.; Werk, A.N.; Cascorbi, I. MicroRNAs and their relevance to ABC transporters. Br. J. Clin. Pharmacol. 2014, 77, 587–596. [Google Scholar] [CrossRef] [PubMed]
  24. Hogan, E.M.; Casserly, A.P.; Scofield, M.D.; Mou, Z.; Zhao-Shea, R.; Johnson, C.W.; Tapper, A.R.; Gardner, P.D. miRNAome analysis of the mammalian neuronal nicotinic acetylcholine receptor gene family. RNA 2014, 20, 1890–1899. [Google Scholar] [CrossRef] [PubMed]
  25. Ma, K.; Li, X.; Hu, H.; Zhou, D.; Sun, Y.; Ma, L.; Zhu, C.; Shen, B. Pyrethroid-resistance is modulated by miR-92a by targeting CpCPR4 in Culex pipiens pallens. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 2017, 203, 20–24. [Google Scholar] [CrossRef] [PubMed]
  26. Zhu, B.; Li, X.; Liu, Y.; Gao, X.; Liang, P. Global identification of microRNAs associated with chlorantraniliprole resistance in diamondback moth Plutella xylostella (L.). Sci. Rep. 2017, 7, 40713. [Google Scholar] [CrossRef] [PubMed]
  27. Derecka, K.; Blythe, M.J.; Malla, S.; Genereux, D.P.; Guffanti, A.; Pavan, P.; Moles, A.; Snart, C.; Ryder, T.; Ortori, C.A.; et al. Transient exposure to low levels of insecticide affects metabolic networks of honeybee larvae. PLoS ONE 2013, 8, e68191. [Google Scholar] [CrossRef] [PubMed]
  28. Kugler, J.M.; Verma, P.; Chen, Y.W.; Weng, R.; Cohen, S.M. miR-989 is required for border cell migration in the Drosophila ovary. PLoS ONE 2013, 8, e67075. [Google Scholar] [CrossRef]
  29. Jagadeeswaran, G.; Zheng, Y.; Sumathipala, N.; Jiang, H.; Arrese, E.L.; Soulages, J.L.; Zhang, W.; Sunkar, R. Deep sequencing of small RNA libraries reveals dynamic regulation of conserved and novel microRNAs and microRNA-stars during silkworm development. BMC Genom. 2010, 11, 52. [Google Scholar] [CrossRef] [PubMed]
  30. Vilmos, P.; Bujna, A.; Szuperák, M.; Havelda, Z.; Várallyay, É.; Szabad, J.; Kucerova, L.; Somogyi, K.; Kristó, I.; Lukácsovich, T.; et al. Viability, longevity, and egg production of Drosophila melanogaster are regulated by the miR-282 microRNA. Genetics 2013, 195, 469–480. [Google Scholar] [CrossRef] [PubMed]
  31. Puinean, A.M.; Foster, S.P.; Oliphant, L.; Denholm, I.; Field, L.M.; Millar, N.S.; Williamson, M.S.; Bass, C. Amplification of a cytochrome P450 gene is associated with resistance to neonicotinoid insecticides in the aphid Myzus persicae. PLoS Genet. 2010, 6, e1000999. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Peng, T.; Pan, Y.; Gao, X.; Xi, J.; Zhang, L.; Ma, K.; Wu, Y.; Zhang, J.; Shang, Q. Reduced abundance of the CYP6CY3-targeting let-7 and miR-100 miRNAs accounts for host adaptation of Myzus persicae nicotianae. Insect Biochem. Mol. Biol. 2016, 75, 89–97. [Google Scholar] [CrossRef] [PubMed]
  33. Christen, V.; Mittner, F.; Fent, K. Molecular Effects of Neonicotinoids in HoneyBees (Apis mellifera). Environ. Sci. Technol. 2016, 50, 4071–4081. [Google Scholar] [CrossRef] [PubMed]
  34. Riaz, M.A.; Chandor-Proust, A.; Dauphin-Villemant, C.; Poupardin, R.; Jones, C.M.; Strode, C.; Régent-Kloeckner, M.; David, J.P.; Reynaud, S. Molecular mechanisms associated with increased tolerance to the neonicotinoid insecticide imidacloprid in the dengue vector Aedes aegypti. Aquat. Toxicol. 2013, 126, 326–337. [Google Scholar] [CrossRef] [PubMed]
  35. Kalajdzic, P.; Oehler, S.; Reczko, M.; Pavlidi, N.; Vontas, J.; Hatzigeorgiou, A.G.; Savakis, C. Use of mutagenesis, genetic mapping and next generation transcriptomics to investigate insecticide resistance mechanisms. PLoS ONE 2012, 7, e40296. [Google Scholar] [CrossRef] [PubMed]
  36. Sawczyn, T.; Dolezych, B.; Klosok, M.; Augustyniak, M.; Stygar, D.; Buldak, R.J.; Kukla, M.; Michalczyk, K.; Karcz-Socha, I.; Zwirska-Korczala, K. Alteration of carbohydrates metabolism and midgut glucose absorption in Gromphadorhina portentosa after subchronic exposure to imidacloprid and fenitrothion. J. Environ. Sci. Health Part A 2012, 47, 1644–1651. [Google Scholar] [CrossRef] [PubMed]
  37. Wan, P.J.; Fu, K.Y.; Lü, F.G.; Wang, X.X.; Guo, W.C.; Li, G.Q. Knocking down a putative Δ(1)-pyrroline-5-carboxylate dehydrogenase gene by RNA interference inhibits flight and causes adult lethality in the Colorado potato beetle Leptinotarsa decemlineata (Say). Pest Manag. Sci. 2015, 71, 1387–1396. [Google Scholar] [CrossRef] [PubMed]
  38. Morin, M.D.; Frigault, J.J.; Lyons, P.J.; Crapoulet, N.; Boquel, S.; Storey, K.B.; Morin, P.J. Amplification and quantification of cold-associated microRNAs in the Colorado potato beetle (Leptinotarsa decemlineata) agricultural pest. Insect Mol. Biol. 2017, 26, 574–583. [Google Scholar] [CrossRef] [PubMed]
  39. Lyons, P.J.; Crapoulet, N.; Storey, K.B.; Morin, P.J. Identification and profiling of miRNAs in the freeze-avoiding gall moth Epiblema scudderiana via next-generation sequencing. Mol. Cell. Biochem. 2015, 410, 155–163. [Google Scholar] [CrossRef] [PubMed]
  40. Lyons, P.J.; Govaere, L.; Crapoulet, N.; Storey, K.B.; Morin, P.J. Characterization of cold-associated microRNAs in the freeze-tolerant gall fly Eurosta solidaginis using high-throughput sequencing. Comp. Biochem. Phys. D 2016, 20, 95–100. [Google Scholar] [CrossRef] [PubMed]
  41. Biggar, K.K.; Kornfield, S.F.; Storey, K.B. Amplification and sequencing of mature microRNAs in uncharacterized animal models using stem-loop reverse transcription-polymerase chain reaction. Anal. Biochem. 2011, 416, 231–233. [Google Scholar] [CrossRef] [PubMed]
  42. Lang-Ouellette, D.; Morin, P.J. Differential expression of miRNAs with metabolic implications in hibernating thirteen-lined ground squirrels, Ictidomys tridecemlineatus. Mol. Cell. Biochem. 2014, 394, 291–298. [Google Scholar] [CrossRef] [PubMed]
  43. Andersen, C.L.; Ledet-Jensen, J.; Ørntoft, T. Normalization of real-time quantitative RT-PCR data: A model based variance estimation approach to identify genes suited for normalization, applied bladder- and colon-cancer data-sets. Cancer Res. 2004, 64, 5245–5250. [Google Scholar] [CrossRef] [PubMed]
  44. Ruby, J.G.; Stark, A.; Johnston, W.K.; Kellis, M.; Bartel, D.P.; Lai, E.C. Evolution, biogenesis, expression, and target predictions of a substantially expanded set of Drosophila microRNAs. Genome Res. 2007, 17, 1850–1864. [Google Scholar] [CrossRef] [PubMed]
  45. Enright, A.J.; John, B.; Gaul, U.; Tuschl, T.; Sander, C.; Marks, D.S. MicroRNA targets in Drosophila. Genome Biol. 2003, 5. [Google Scholar] [CrossRef] [PubMed]
  46. Huang da, W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 445–447. [Google Scholar] [CrossRef] [PubMed]
  47. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  48. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Biogeosciences 2010, 26, 139–140. [Google Scholar] [CrossRef] [PubMed]
  49. Smyth, G.K. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat. Appl. Genet. Mol. Biol. 2004, 3, 1–25. [Google Scholar] [CrossRef] [PubMed]
Ijms 18 02728 g001 550
Figure 1. Next-generation sequencing (NGS) and quantitative real-time polymerase chain reaction (qRT-PCR) expression data of select miRNAs. Expression of eight miRNAs quantified in control and imidacloprid-treated L. decemlineata. Data for qRT-PCR is normalized to transcript levels (mean ± SEM, n = 3) and NGS levels represent fold-changes in normalized read counts (mean ± SEM, n = 3).
Figure 1. Next-generation sequencing (NGS) and quantitative real-time polymerase chain reaction (qRT-PCR) expression data of select miRNAs. Expression of eight miRNAs quantified in control and imidacloprid-treated L. decemlineata. Data for qRT-PCR is normalized to transcript levels (mean ± SEM, n = 3) and NGS levels represent fold-changes in normalized read counts (mean ± SEM, n = 3).
Ijms 18 02728 g001
Table
Table 1. Mapping statistics of reads in control versus imidacloprid-treated L. decemlineata.
Table 1. Mapping statistics of reads in control versus imidacloprid-treated L. decemlineata.
Types of RNAsUnique Reads Control/TreatedPercent Control/TreatedRead Count Control/TreatedPercent Control/Treated
Total reads8,825,610/6,211,853100%/100%22,534,030/13,763,160100%/100%
miRNAs26,447/24,2210.30%/0.39%1,558,719/1,511,2016.92%/10.98%
tRNAs40,144/22,1150.45%/0.36%326,281/125,0851.45%/0.91%
snRNAs4705/29370.05%/0.05%11,161/63020.05%/0.05%
snoRNAs4735/36160.05%/0.06%10,813/69100.05%/0.05%
rRNAs102,238/64,1891.16%/1.03%619,781/172,9692.75%/1.26%
Table
Table 2. miRNAs strongly expressed in L. decemlineata.
Table 2. miRNAs strongly expressed in L. decemlineata.
miRNAsmiRNA SequencesNormalized Expression
miR-14-3pUCAGUCUUUUUCUCUCUCCUAU88,717.57
miR-8-3pUAAUACUGUCAGGUAAAGAUGUC73,723.79
miR-276-3pUAGGAACUUCAUACCGUGCUCU54,681.67
miR-317-3pUGAACACAGCUGGUGGUAUCUCAGU39,445.82
miR-1-3pUGGAAUGUAAAGAAGUAUGGAG16,976.97
miR-2a-3pUAUCACAGCCAGCUUUGAUGAGC16,697.46
miR-281-5pAAGAGAGCUAUCCGUCGACAGU16,606.36
miR-1175-3pUGAGAUUCAACUCCUCCAUCUC14,289.62
miR-13b-3pUAUCACAGCCAUUUUGACGAGU14,017.42
bantam-3pUGAGAUCAUUGUGAAAGCUGAUU13,166.75
Presented average normalized expressions are from all samples characterized by high-throughput sequencing.
Table
Table 3. Differential expression of 33 miRNAs following imidacloprid treatment in L. decemlineata.
Table 3. Differential expression of 33 miRNAs following imidacloprid treatment in L. decemlineata.
miRNAsmiRNA SequencesNormalized Expression Treated/ControlLog2 Fold-Change Treated/Control
miR-iab-8-3pAGGAUACAUUCAGUAUACG27.7753.710.95
miR-252a-3pCCUGCUGCUCAAGUGCUUAUC12.4021.640.80
miR-282-5pUAGCCUCUCCUAGGCUUUGUCU13.2320.320.62
miR-iab-8-5pUUACGUAUACUGAAGGUAUACCGGAC30.5444.060.53
miR-1000-5pAUAUUGUCCUGUCACAGC63.2986.100.44
miR-3849-5pUGACAUUUUAACCAUAGUGCU59.9781.160.44
miR-193-3pUACUGGCCUGUUAAGUCCCAAGU57.3177.140.43
miR-7-3pCAAGGAAUCACUAAUCAUCCCAC36.7949.250.42
miR-124-3pUAAGGCACGCGGUGAAUGCCAAG183.64239.110.38
miR-970-3pUCAUAAGACACACGCGGCUAU337.13427.790.34
miR-276-5pAGCGAGGUAUAGAGUUCCUACGUG2120.162672.970.33
let-7-5pUGAGGUAGUAGGUUGUAUAG4613.215776.290.32
miR-1-3pUGGAAUGUAAAGAAGUAUGGAG15,177.8818,776.070.31
miR-133-3pUUGGUCCCCUUCAACCAGCUGU1771.682187.690.30
miR-263b-5pCUUGGCACUGGAAGAAUUCAC95.9176.77−0.32
miR-281-5pAAGAGAGCUAUCCGUCGACAGU18,510.9114,701.80−0.33
miR-2944c-3pUAUCACAGCCAGUAGUUACC4468.293543.74−0.33
miR-1175-5pAAGUGGAGCAGUGGUCUCUUCAC286.32224.91−0.35
miR-92a-3pAUUGCACUAGUCCCGGCCUA62.5248.60−0.36
miR-305-5pAUUGUACUUCAUCAGGUGCUC7239.925605.19−0.37
miR-iab-4-3pCGGUAUACCUUCAGUAUACGUAAC24.5618.65−0.40
miR-927a-5pUUUAGAAUUCCUACGCUUUA20.1014.98−0.42
miR-9c-5pUCUUUGGUGAUCUAGCCGUGUG405.18297.72−0.44
miR-34-3pCGACCACUAUCCAUACUCCCUCC29.1121.28−0.45
miR-316-5pUGUCUUUUUCCGCUUUGCUGC11,436.588348.15−0.45
miR-998-3pUAGCACCAUGGGAUUCAGCUCA87.9059.03−0.57
miR-100-5pAACCCGUAGAUCCGAACUUGUGGG5341.403528.46−0.60
miR-750-3pCCAGAUCUAACUCUUCCAUAUGACG6089.673795.84−0.68
miR-995-3pUAGCACCACAUGAUUCAGCUUACG551.70340.13−0.70
miR-2796-5pAGGGGUUUCUUUCGGCCUCCAGCG41.1924.70−0.74
miR-92a-5pAGUCCGUGAUGCGUGACAAUAU197.04117.23−0.75
miR-315-5pUUUUGAUUGUUGCUCAGAAAGC21.8012.80−0.77
miR-989-3pUGUGAUGUGACGUAGUGG4764.902748.54−0.79
Table
Table 4. Target prediction of select L. decemlineata miRNAs with TargetScanFly and miRanda.
Table 4. Target prediction of select L. decemlineata miRNAs with TargetScanFly and miRanda.
miRNAsTargetScanFly TargetsMiRanda Targets
miR-282-5pEro1L, Fus, CG14435, CG9515, CcapRMeso18E, Rogdi, Resilin, Nkd, Kraken
miR-1000-5pCG34355, Nplp1, CG10804, Net, CG13384Nplp1, CG10804, Kon, CG13384, Net
miR-193-3pP38b, Cp7Fb, Ana, CG11041, CG11313CG34394, RYBP, CG6707, Pb, CG32736
miR-124-3pSinu, Gli, CG12977, Pk92B, AxsPk92B, CG14299, Sinu, Gli, Cp110
miR-970-3pBtsz, CG15097, Ru, Ace, Rab6CG15097, Ubc-E2H, CG32372, CG42256, StmA
let-7-5pAb, CG18265, A3-3, Apt, Cpr49AaAb, CG5026, CG34118, CG9548, Slam
miR-1-3pCG4297, Tub, CG30457, Hmu, CG6490Par-6, CG18542, Pen, CG5053, Pdm2
miR-133-3pFili, CG17193, CG2774, CG33324, CG9541CG30409, SkpA, Wts, Pde1c, CG32351
miR-263b-5pQkr54B, Bx, Wls, CG32062, CG34339CG2371, Wls, Qkr54B, Bx, GATAe
miR-92a-3pSha, Khc-73, CrebA, CG4297, TuspCrebA, Cpr50Ca, CG8128, CG3077, CG14274
miR-305-5pCG33174, CG11997, CG3287, Gr98d, Nerfin-1CG3287, CG31855, Mi-2, NfI, Eya
miR-927a-5pGprs, CG32245, Kr-h1, CG8485, GrowlCG8485, EcR, CG9850, Aef1, Sfl
miR-9c-5pNerfin-1, Rbp9, CG11206, CG32333, Bru-2Sinu, CadN, CG5746, CG11206, CG9426
miR-316-5pHr39, Rbp9, CG32121, Numb, Bsg25DNkd, Numb, Aret, Cib, RdgC
miR-100-5pE(Pc), Gogo, CG17985, CG31772, DopR2CG3630, CG10979, Gogo, Pc, CG13326
miR-315-5pEip93F, CG15465, CG32333, CG32137, CG32206CG12424, CG14989, CG32137, Rtnl1, CG34126
miR-989-3pTankyrase, Fal, Lac, PGRP-SD, Qkr54BCG12772, Chrw, CG34449, Nedd4, Kni
Table
Table 5. Biological processes predicted to be modulated by select miRNAs with varying expression levels in imidacloprid-treated L. decemlineata.
Table 5. Biological processes predicted to be modulated by select miRNAs with varying expression levels in imidacloprid-treated L. decemlineata.
GO Biological ProcessGenes TargetedmiRNAsp-Value
Sensory perception of pain13116.0 × 10−3
Regulation of transcription, DNA-templated12107.2 × 10−4
Transcription, DNA-templated1191.6 × 10−3
Regulation of glucose metabolic process971.5 × 10−4
Positive regulation of transcription RNA pol II promoter873.8 × 10−3
Border follicle cell migration665.7 × 10−3
Axon guidance652.3 × 10−2
Imaginal disc-derived wing morphogenesis654.1 × 10−2
Muscle organ development565.9 × 10−3
Negative regulation of transcription, DNA-templated552.1 × 10−2
Negative regulation of transcription RNA pol II promoter554.1 × 10−2
Table
Table 6. Sequences, efficiencies, and optimal melting temperatures of primers used in this study.
Table 6. Sequences, efficiencies, and optimal melting temperatures of primers used in this study.
PrimerSequenceEff.Temp.
miR-1-3p5′-ACACTCCAGCTGGGTGGAATGTAAAGAAGTA-3′92.0%62.5 °C
5′-CTCACAGTACGTTGGTATCCTTGTGATGTTCGATGCCATATTGTACTGTGAGCTCCATAC-3′
miR-7-5p5′-ACACTCCAGCTGGGTGGAAGACTAGTGAT-3′96.0%60.0 °C
5′-CTCACAGTACGTTGGTATCCTTGTGATGTTCGATGCCATATTGTACTGTGAGCACAACAA-3′
miR-9a-5p5′-ACACTCCAGCTGGGTCTTTGGTTATCTAG-3′101.4%58.8 °C
5′-CTCACAGTACGTTGGTATCCTTGTGATGTTCGATGCCATATTGTACTGTGAGTCATACAG-3′
miR-12-5p5′-ACACTCCAGCTGGGTGAGTATTACATCAGGT-3′95.6%64.5 °C
5′-CTCACAGTACGTTGGTATCCTTGTGATGTTCGATGCCATATTGTACTGTGAGCAGTACCT-3′
miR-92a-5p5′-ACACTCCAGCTGGGAGTCCGTGATGCGTGAC-3′88.5%56.5°C
5′-CTCACAGTACGTTGGTATCCTTGTGATGTTCGATGCCATATTGTACTGTGAGATATTGTC-3′
miR-133-3p5′-ACACTCCAGCTGGGTTGGTCCCCTTCAACCA-3′84.3%64.5 °C
5′-CTCACAGTACGTTGGTATCCTTGTGATGTTCGATGCCATATTGTACTGTGAG-3′
miR-305-5p5′-ACACTCCAGCTGGGATTGTACTTCATCAGGT-3′91.3%63.4 °C
5′-CTCACAGTACGTTGGTATCCTTGTGATGTTCGATGCCATATTGTACTGTGAGGAGCACCT-3′
miR-927a-5p5′-ACACTCCAGCTGGGTTTAGAATTCCTACGCT-3′105.3%60.9 °C
5′-CTCACAGTACGTTGGTATCCTTGTGATGTTCGATGCCATATTGTACTGTGAGTAAAGCGT-3′
miR-989-3p5′-ACACTCCAGCTGGGTGTGATGTGACGTAGTG-3′99.5%56.9 °C
5′-CTCACAGTACGTTGGTATCCTTGTGATGTTCGATGCCATATTGTACTGTGAGCCACTACG-3′
Universal5′-CTCACAGTACGTTGGTATCCTTGTG-3′--
Top and bottom sequences represent forward and stem-loop primers, respectively. Abbreviations: Efficiency (Eff.); Temperature (Temp.).

Share and Cite

MDPI and ACS Style

Morin, M.D.; Lyons, P.J.; Crapoulet, N.; Boquel, S.; Morin, P.J. Identification of Differentially Expressed miRNAs in Colorado Potato Beetles (Leptinotarsa decemlineata (Say)) Exposed to Imidacloprid. Int. J. Mol. Sci. 2017, 18, 2728. https://doi.org/10.3390/ijms18122728

AMA Style

Morin MD, Lyons PJ, Crapoulet N, Boquel S, Morin PJ. Identification of Differentially Expressed miRNAs in Colorado Potato Beetles (Leptinotarsa decemlineata (Say)) Exposed to Imidacloprid. International Journal of Molecular Sciences. 2017; 18(12):2728. https://doi.org/10.3390/ijms18122728

Chicago/Turabian Style

Morin, Mathieu D., Pierre J. Lyons, Nicolas Crapoulet, Sébastien Boquel, and Pier Jr Morin. 2017. "Identification of Differentially Expressed miRNAs in Colorado Potato Beetles (Leptinotarsa decemlineata (Say)) Exposed to Imidacloprid" International Journal of Molecular Sciences 18, no. 12: 2728. https://doi.org/10.3390/ijms18122728

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop