Identification, characterization and target gene analysis of testicular microRNAs in the oriental fruit fly Bactrocera dorsalis

RNAinsecta: A tool for prediction of precursor microRNA in insects and search for their target in the model organism Drosophila melanogaster

INTRODUCTION AND BACKGROUND

Pre-MicroRNAs are the hairpin loops from which microRNAs are produced that have been found to negatively regulate gene expression in several organisms. In insects, microRNAs participate in several biological processes including metamorphosis, reproduction, immune response, etc. Numerous tools have been designed in recent years to predict novel pre-microRNA using binary machine learning classifiers where prediction models are trained with true and pseudo pre-microRNA hairpin loops. Currently, there are no existing tool that is exclusively designed for insect pre-microRNA detection.

AIM

Application of machine learning algorithms to develop an open source tool for prediction of novel precursor microRNA in insects and search for their miRNA targets in the model insect organism, Drosophila melanogaster.

METHODS

Machine learning algorithms such as Random Forest, Support Vector Machine, Logistic Regression and K-Nearest Neighbours were used to train insect true and false pre-microRNA features with 10-fold Cross Validation on SMOTE and Near-Miss datasets. miRNA targets IDs were collected from miRTarbase and their corresponding transcripts were collected from FlyBase. We used miRanda algorithm for the target searching.

RESULTS

In our experiment, SMOTE performed significantly better than Near-Miss for which it was used for modelling. We kept the best performing parameters after obtaining initial mean accuracy scores >90% of Cross Validation. The trained models on Support Vector Machine achieved accuracy of 92.19% while the Random Forest attained an accuracy of 80.28% on our validation dataset. These models are hosted online as web application called RNAinsecta. Further, searching target for the predicted pre-microRNA in Drosophila melanogaster has been provided in RNAinsecta.

Integrated analysis of miRNA and mRNA expression profiles in response to gut microbiota depletion in the abdomens of female Bactrocera dorsalis

Insect gut microbiota has been reported to participate in regulating host multiple biological processes including metabolism and reproduction. However, the corresponding molecular mechanisms remain largely unknown. Recent studies suggest that microRNAs (miRNAs) are involved in complex interactions between the gut microbiota and the host. Here, we used next-generation sequencing technology to characterize miRNA and mRNA expression profiles and construct the miRNA-gene regulatory network in response to gut microbiota depletion in the abdomens of female Bactrocera dorsalis. A total of 3016 differentially expressed genes (DEGs) and 18 differentially expressed miRNAs (DEMs) were identified. Based on the integrated analysis of miRNA and mRNA sequencing data, 229 negatively correlated miRNA-gene pairs were identified from the miRNA-mRNA network. Gene ontology enrichment analysis indicated that DEMs could target several genes involved in the metabolic process, oxidation-reduction process, oogenesis, and insulin signaling pathway. Finally, real-time quantitative polymerase chain reaction further verified the accuracy of RNA sequencing results. In conclusion, our study provides the profiles of miRNA and mRNA expressions under antibiotics treatment and provides an insight into the roles of miRNAs and their target genes in the interaction between the gut microbiota and its host.

Exploring miRNA-mRNA regulatory modules responding to tannic acid stress in Micromelalopha troglodyta (Graeser) (Lepidoptera: Notodontidae) via small RNA sequencing

MicroRNAs (miRNAs) are small noncoding RNAs (sRNAs) that regulate gene expression by inhibiting translation or degrading mRNA. Although the functions of miRNAs in many biological processes have been reported, there is currently no research on the possible roles of miRNAs in Micromelalopha troglodyta (Graeser) involved in the response of plant allelochemicals. In this article, six sRNA libraries (three treated with tanic acid and three control) from M. troglodyta were constructed using Illumina sequencing. From the results, 312 known and 43 novel miRNAs were differentially expressed. Notably, some of the most abundant miRNAs, such as miR-432, miR-541-3p, and miR-4448, involved in important physiological processes were also identified. To better understand the function of the targeted genes, we performed Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. The results indicated that differentially expressed miRNA targets were involved in metabolism, development, hormone biosynthesis, and immunity. Finally, we visualized a miRNA-mRNA regulatory module that supports the role of miRNAs in host-allelochemical interactions. To our knowledge, this is the first report on miRNAs responding to tannic acid in M. troglodyta. This study provides indispensable information for understanding the potential roles of miRNAs in M. troglodyta and the applications of these miRNAs in M. troglodyta management.

miR-125-3p and miR-276b-3p Regulate the Spermatogenesis of Bactrocera dorsalis by Targeting the orb2 Gene

Bactrocera dorsalis is considered a major threat to horticultural crops. It has evolved resistance against insecticides. It is believed that development of new methods is highly desirable to control this destructive agricultural pest. Sterile insect technique is emerging as a potential tool to control this insect pest by reducing their reproductive ability. Here we report that orb2 has high expression in the testis of B. dorsalis which is the target of miR-125-3p and miR-276b-3p and plays a critical role in the spermatogenesis. Dual luciferase reporter assay using HEKT293 cells demonstrates that orb2 gene is downregulated by miR-125-3p and miR-276b-3p and is a common target of these miRNAs. Dietary treatment of adult male flies separately and in combination of agomir-125-3p (Ago-125-3p) and agomir-276b-3p (Ago-276b-3p) significantly downregulated the mRNA of orb2. The combined treatments of agomirs suppressed the level of mRNA of orb2 significantly more than any single treatment. Altered expression of miR-125-3p and miR-276b-3p significantly decreased the total and live spermatozoa in the testis which ultimately caused reduction in male fertility. Furthermore, we demonstrate that miR-125-3p, miR-276b-3p, and orb2 dsRNA are the novel agents that could be used in a genetic-based sterile insect technique (SIT) to control the B. dorsalis.

RNA interference of Argonaute-1 delays ovarian development in the oriental fruit fly, Bactrocera dorsalis (Hendel)

BACKGROUND

With the development of rapid resistance, new modes of action for pesticides are needed for insect control, such as RNAi-based biopesticides targeting essential genes. To explore the function of Argonaute-1 (Ago-1) and potential miRNAs in ovarian development of Bactrocera dorsalis, an important agricultural pest, and to develop a novel control strategy for the pest, BdAgo-1 was first identified in B. dorsalis.

RESULTS

Spatiotemporal expression analysis indicated that BdAgo-1 had a relatively high transcriptional level in the ovarian tissues of adult female B. dorsalis during the sexual maturation period. RNA interference (RNAi) experiment showed that BdAgo-1 knockdown significantly decreased the expression levels of ovarian development-related genes and delayed ovarian development. Although RNAi-mediated silencing of Ago-1 led to a reduced ovary surface area, a subsequent oviposition assay revealed that the influence was minimal over a longer time period. Small RNA libraries were constructed and sequenced from different ovarian developmental stages of B. dorsalis adults. Among 161 identified miRNAs, 84 miRNAs were differentially expressed during the three developmental stages of the B. dorsalis ovary. BdAgo-1 silencing caused significant down-regulation of seven differentially expressed miRNAs (DEMs) showing relatively high expression levels (>1000 TPM (Transcripts per kilobase of exon model per million mapped reads)). The expression patterns of these seven core DEMs and their putative target genes were analyzed in the ovaries of B. dorsalis.

CONCLUSION

The results indicate that Ago-1 and Ago-1-dependent miRNAs are indispensable for normal ovarian development in B. dorsalis and help identify miRNA targets useful for control of this pest.

Identification and profiling of Bactrocera dorsalis microRNAs and their potential roles in regulating the developmental transitions of egg hatching, molting, pupation and adult eclosion

MicroRNAs (miRNAs) are endogenous small noncoding RNAs (18-25 nt) that are involved in many physiological processes including development, cancer, immunity, apoptosis and host-microbe interactions through post-transcriptional regulation of gene expression. In this study, we measured the profile of small RNAs over the developmental transitions of the oriental fruit fly Bactrocera dorsalis from egg hatching, molting, and pupation to adult eclosion. We identified 250 miRNAs, including 83 known and 167 novel miRNAs, and 47 isomiRNAs. In addition, we identified the miRNAs differentially expressed over the developmental transitions. Interestingly, the miR-309 cluster, the miR-2 cluster/family and the let-7 cluster were among these differentially expressed miRNAs, suggesting a role in the regulation of egg hatching, molting and pupation/adult eclosion, respectively. Moreover, a detailed analysis of the temporal expression patterns of 14 highly expressed miRNAs in the pupal stage revealed three types of expression profiles. Furthermore, injection of a miR-100 mimic in the 3rd instar larvae resulted in a significant decrease in pupation and adult eclosion rates, whereas injection of a miR-317 antagomir resulted in a significant decrease in the pupation rate and a decrease in the pupation time, indicating that miR-100 and miR-317 are involved in the process of pupation. Finally, injection of a miR-100/miR-285 mimic or antagomir in pupae resulted in a significant decrease in the eclosion rate and a significant increase in the prevalence of a partial eclosion phenotype, implying the involvement of miR-100 and miR-285 in the process of adult eclosion. This study identified critical miRNAs involved in the transitions of this important holometabolic model and pest insect B. dorsalis from egg hatching to adult eclosion, thus providing a useful resource for exploring the regulatory role of miRNAs during insect post-embryonic development.

Transcriptome-Wide Identification of MicroRNAs and Analysis of Their Potential Roles in Development of Indian Meal Moth (Lepidoptera:Pyralidae)

MicroRNAs (miRNAs) have been reported to play indispensable roles in regulating various developmental processes via the posttranscriptional repression of target genes in insect species. In the present paper, we studied the miRNAs in Indian meal moth (Plodia interpunctella (Hübener)), one of the most economically important stored grains pests around the world. In total, 12 small RNA libraries from four developmental stages of P. interpunctella were constructed, and 178 known and 23 novel miRNAs were identified. In addition, the expression profiles of these miRNAs were assessed across different developmental stages and miRNAs that were highly expressed in eggs, larvae, pupae, and adults were identified. Specifically, 100, 61, and 52 miRNAs were differentially expressed between eggs and larvae, larvae and pupae, and pupae and adults, respectively. The KEGG and GO analysis of the predicted target genes suggested the essential roles of miRNAs in the regulation of complex development of P. interpunctella. Importantly, we also found a set of miRNAs might be involved in the larval metamorphic molting process, with their expressions increasing and then decreasing during the larva-pupa-adult stages of P. interpunctella. In conclusion, the current paper has discovered numerous miRNAs, and some key miRNAs that might be responsible for regulating development in P. interpunctella. To our knowledge, this is the first study to document miRNAs and their expression patterns in interpunctella, and those findings would lay an important molecular foundation for future functional analysis of these miRNAs in P. interpunctella.

MicroRNA Let-7 targets the ecdysone signaling pathway E75 gene to control larval-pupal development in Bactrocera dorsalis

MicroRNAs (miRNAs) regulate various biological processes during insect development; however, their role in larval-pupal development in oriental fruit fly, Bactrocera dorsalis (Hendel) remains unknown. In the current study, we address the biological function of a conserved miRNA, Bdo-Let-7 in the regulation of BdE75 gene, which belongs to the ecdysone signaling pathway and participates in the larval-pupal development in B. dorsalis. Using dual luciferase reporter assay in HEK293T cells we show that Bdo-Let-7 miRNA interacts with the 3' untranslated region of BdE75 gene and suppresses its expression. The Bdo-Let-7 and BdE75 are also co-expressed in the larval-pupal stages and in different tissues of B. dorsalis. In in vivo experiments, the injection of Bdo-Let-7 agomir and antagomir in third instar larvae down- and up-regulated the expression of BdE75, respectively. The 20-hydroxyecdysone (20E) injection assay shows that 20E up-regulated the expression of Bdo-Let-7 on the 5th day of the larvae. Moreover, abnormal pupation and eclosion were observed after larval Bdo-Let-7 antagomir injection. Based on these results, we show that Bdo-Let-7 regulates the ecdysone signaling pathway through the exact dose of BdE75 gene, and is indispensable for normal larval-pupal development in B. dorsalis.

Integrated Analysis of miRNA and mRNA Expression Profiles Reveals Functional miRNA-Targets in Development Testes of Small Tail Han Sheep

Small Tail Han Sheep is a highly valued local breed in China because of their precocity, perennial estrus, and high fecundity. The average annual lambing rate of ewes is as high as 180–270%, the semen of ram has characteristics of high yield, high density, and good motility. To reveal the key miRNAs and miRNA-targets underlying testis development and spermatogenesis in male Small Tail Han Sheep, integrated analysis of miRNA and mRNA expression profiles in 2-, 6-, and 12-month-old testes was performed by RNA-seq technology and bioinformatics methods. The results showed that total of 153 known sheep miRNAs and 2712 novel miRNAs were obtained in 2-,6 - and 12-month-old Small Tail Han Sheep testes; 5, 1, and 4 differentially expressed (DE) known sheep miRNAs, and 132, 105, and 24 DE novel miRNAs were identified in 2- vs. 6-, 6- vs. 12-, and 2- vs. 12-month-old testes, respectively. We combined miRNA results of this study and the mRNA results obtained in our previous study to predict the target mRNAs of DE known sheep miRNAs; 131, 10, and 15 target mRNAs of DE known sheep miRNAs and 76, 1, and 11 DE miRNA–targets were identified in the three groups, respectively. GO and KEGG analyses showed that: in 2- vs. 6-month-olds, the target genes of DE known sheep miRNAs were involved in 100 biological processes and 11 signaling pathways; in 6- vs. 12-month-olds, the target genes of DE known sheep miRNAs were involved in 4 biological processes; and in 2- vs. 12-month-olds, the target genes of DE known sheep miRNAs were involved in 17 biological processes and 4 signaling pathways. Three miR–target regulatory networks were constructed based on these DE miRNA–targets. The key miRNA-Targets involved in testis development and spermatogenesis were screened. 6 known sheep miRNAs and 6 novel miRNAs were selected to validate the accuracy of miRNA sequencing data by qRT-PCR. The binding sites of oar-miR-379-5p with WNT8A was validated by a dual luciferase reporter gene detection system.

Genome-Wide Screening and Functional Analysis Reveal That the Specific microRNA nlu-miR-173 Regulates Molting by Targeting Ftz-F1 in Nilaparvata lugens

Background: Molting is a crucial physiological behavior during arthropod growth. In the past few years, molting as well as chitin biosynthesis triggered by molting, is subject to regulation by miRNAs. However, how many miRNAs are involved in insect molting at the genome-wide level remains unknown. Results: We deeply sequenced four samples obtained from nymphs at the 2nd-3rd and 4th-5th instars, and then identified 61 miRNAs conserved in the Arthropoda and 326 putative novel miRNAs in the brown planthopper Nilaparvata lugens, a fearful pest of rice. A total of 36 mature miRNAs with significant different expression levels at the genome scale during molting, including 19 conserved and 17 putative novel miRNAs were identified. After comparing the expression profiles, we found that most of the targets of 36 miRNAs showing significantly differential expression were involved in energy and hormone pathways. One of the 17 putative novel miRNAs, nlu-miR-173 was chosen for functional study. nlu-miR-173 acts in 20-hydroxyecdysone signaling through its direct target, N. lugens Ftz-F1(NlFtz-F1), a transcription factor. Furthermore, we found that the transcription of nlu-miR-173 was promoted by Broad-Complex (BR-C), suggesting that its involvement in the 20-hydroxyecdysone pathway contributes to proper molting function. Conclusion: We provided a comprehensive resource of miRNAs associated with insect molting and identified a novel miRNA as a potential target for pest control.

Assembly and Annotation of Transcriptome Provided Evidence of miRNA Mobility between Wheat and Wheat Stem Sawfly

Wheat Stem Sawfly (WSS), Cephus Cinctus Norton (Hymenoptera: Cephidae), is one of the most important pests, causing yield and economic losses in wheat and barley. The lack of information about molecular mechanisms of WSS for defeating plant's resistance prevents application of effective pest control strategies therefore, it is essential to identify the genes and their regulators behind WSS infestations. Long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) are recognized with their regulatory functions on gene expression, tuning protein production by controlling transcriptional and post-transcriptional activities. A transcriptome-guided approach was followed in order to identify miRNAs, lncRNAs, and mRNA of WSS, and their interaction networks. A total of 1,893 were presented here as differentially expressed between larva and adult WSS insects. There were 11 miRNA families detected in WSS transcriptome. Together with the annotation of 1,251 novel mRNAs, the amount of genetic information available for WSS was expanded. The network between WSS miRNAs, lncRNAs, and mRNAs suggested miRNA-mediated regulatory roles of lncRNAs as competing endogenous RNAs. In the light of the previous evidence that small RNA molecules of a pathogen could suppress the immune response of host plant, we analyzed the putative interactions between larvae and wheat at the miRNA level. Overall, this study provides a profile of larva and adult WSS life stages in terms of coding and non-coding elements. These findings also emphasize the potential roles of wheat and larval miRNAs in wheat resistance to infestation and in the suppression of resistance which is critical for the development of effective pest control strategies.

Estradiol regulates expression of miRNAs associated with myogenesis in rainbow trout

17β-Estradiol (E2) is a steroid hormone that negatively affects muscle growth in rainbow trout, but the mechanism associated with this response is not fully understood. To better characterize the effects of E2 on muscle, we identified differentially regulated microRNAs (miRNAs) and muscle atrophy-related transcripts in juvenile rainbow trout exposed to E2. Small RNA-Seq analysis of E2-treated vs. control muscle identified 36 differentially expressed miRNAs including those known to be involved in myogenesis, cell cycle, apoptosis, and cell death. Some important myogenic miRNAs, such as miR-133 and miR-206, are upregulated while others like miR-145 and miR-499, are downregulated. Gene Ontology analysis of the target genes regulated by the miRNAs involved in atrophy and cell cycle indicates that E2 influence leads to expansion of quiescent myogenic precursor cell population to address atrophying mature muscle in rainbow trout during sexual development.

Identification and profiling of miRNAs in overwintering Lissorhoptrus oryzophilus via next-generation sequencing

MicroRNAs (miRNAs) are important regulators of various biological processes in organisms. Insects subjected to abiotic stress can regulate gene expression post-transcriptionally through the use of microRNAs. However, the role of miRNAs in response to cold stress in Lissorhoptrus oryzophilus Kuschel remains unknown. The rice water weevil, L. oryzophilus, is an invasive insect that is able to survive cold winters. To characterize changes in miRNAs in response to overwintering in L. oryzophilus, a comprehensive comparative analysis of microRNAs was performed involving an overwintering and a normal adult. High-throughput Illumina sequencing and bioinformatics analyses revealed 121 conserved and 14 potential novel microRNAs in two small libraries. The novel miRNAs exhibit low expression levels in both libraries. After the expression profiles of the miRNAs in the two libraries were normalized, 36 miRNAs in L. oryzophilus were found to be differentially expressed in response to overwintering. In particular, 14 conserved miRNAs and 6 novel miRNAs were up-regulated, while 15 conserved miRNAs and 1 novel miRNA were down-regulated. In addition, the expression patterns of 11 conserved and potentially novel miRNAs were confirmed by quantitative RT-PCR analysis. Most importantly, this work provides a unique resource of characterized miRNAs for overwintering L. oryzophilus and contributes to studies of the functions of cold-related and other L. oryzophilus miRNAs.

Identification and Characterization of Sex-Biased MicroRNAs in Bactrocera dorsalis (Hendel)

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that regulate various biological processes including sexual dimorphism. The oriental fruit fly Bactrocera dorsalis is one of the most destructive agricultural insect pests in many Asian countries. However, no miRNAs have been identified from the separate sex and gonads to elucidate sex gonad differentiation in B. dorsalis. In this study, we constructed four small RNA libraries from whole body of females, males (except ovaries and testes) and ovaries, testes of B. dorsalis for deep sequencing. The data analysis revealed 183 known and 120 novel miRNAs from these libraries. 18 female-biased and 16 male-biased miRNAs that may be involved in sexual differentiation were found by comparing the miRNA expression profiles in the four libraries. Using a bioinformatic approach, we predicted doublesex (dsx) as a target gene of the female-biased miR-989-3p which is considered as the key switch gene in the sex determination of tephritid insects. This study reveals the first miRNA profile related to the sex differentiation and gives a first insight into sex differences in miRNA expression of B. dorsalis which could facilitate studies of the reproductive organ specific roles of miRNAs.

Testis-specific RNA-Seq of Calligrapha (Chrysomelidae) as a transcriptomic resource for male-biased gene inquiry in Coleoptera

We report the architecture of testis transcriptomes of four closely related species of Calligrapha (Chrysomelidae) beetles, which diverged during the last 3 million years. Five cDNA libraries were sequenced using Illumina HiSeq technology, retrieving 102 884-176 514 assembled contigs, of which ~33-45% of these longer than 499 nt were functionally annotated. Annotation and sequence similarity comparisons of these libraries revealed high homogeneity in gene composition and the presence of several functional candidates related to reproduction or reproductive processes (0.72-1.08% of annotated sequences). Stringent sequence similarity analyses of these transcriptomes against empirically demonstrated male-biased genes in Drosophila melanogaster and Tribolium castaneum allowed the identification of 77 homologues in Calligrapha, possible candidates of male-biased expression. Some of these genes - including CG9313, Tektin-A or tomboy40 - were confirmed as orthologs of these male-biased genes using phylogenetic inference and available model insect data, increasing our confidence that they represent functional homologues too. Our transcriptomes are a valuable transcriptomic resource for the analysis of male-biased genes in Calligrapha, which has the added interest of including several female-only species. But it simultaneously represents a landmark for similar studies in Coleoptera, broadening the taxonomic diversity currently represented by the model species T. castaneum, and incipient genomic data in other herbivorous lineages, including weevils, longhorn beetles and leaf beetles.

miR-8-3p regulates mitoferrin in the testes of Bactrocera dorsalis to ensure normal spermatogenesis

Genetics-enhanced sterile insect techniques (SIT) are promising novel approaches to control Bactrocera dorsalis, the most destructive horticultural pest in East Asia and the Pacific region. To identify novel genetic agents to alter male fertility of B. dorsalis, previous studies investigated miRNA expression in testes of B. dorsalis. One miRNA, miR-8-3p was predicted to bind the 3'UTR of putative B. dorsalis mitoferrin (bmfrn). The ortholog of bmfrn in D. melanogaster is essential for male fertility. Here we show that bmfrn has all conserved amino acid residues of known mitoferrins and is most abundantly expressed in B. dorsalis testes, making miR-8-3p and mitoferrin candidates for genetics-enhanced SIT. Furthermore, using a dual-luciferase reporter system, we show in HeLa cells that miR-8-3p interacts with the 3'UTR of bmfrn. Dietary treatments of adult male flies with miR-8-3p mimic, antagomiR, or bmfrn dsRNA, altered mitoferrin expression in the testes and resulted in reduced male reproductive capacity due to reduced numbers and viability of spermatozoa. We show for the first time that a mitoferrin is regulated by a miRNA and we demonstrate miR-8-3p as well as bmfrn dsRNA to be promising novel agents that could be used for genetics-enhanced SIT.