MicroRNAs of two medically important mosquito species: Aedes aegypti and Anopheles stephensi

Analysis of Amblyomma americanum microRNAs in response to Ehrlichia chaffeensis infection and their potential role in vectorial capacity

Background

MicroRNAs (miRNAs) represent a subset of small noncoding RNAs and carry tremendous potential for regulating gene expression at the post-transcriptional level. They play pivotal roles in distinct cellular mechanisms including inhibition of bacterial, parasitic, and viral infections via immune response pathways. Intriguingly, pathogens have developed strategies to manipulate the host's miRNA profile, fostering environments conducive to successful infection. Therefore, changes in an arthropod host's miRNA profile in response to pathogen invasion could be critical in understanding host-pathogen dynamics. Additionally, this area of study could provide insights into discovering new targets for disease control and prevention. The main objective of the present study is to investigate the functional role of differentially expressed miRNAs upon Ehrlichia chaffeensis, a tick-borne pathogen, infection in tick vector, Amblyomma americanum.

Methods

Small RNA libraries from uninfected and E. chaffeensis-infected Am. americanum midgut and salivary gland tissues were prepared using the Illumina Truseq kit. Small RNA sequencing data was analyzed using miRDeep2 and sRNAtoolbox to identify novel and known miRNAs. The differentially expressed miRNAs were validated using a quantitative PCR assay. Furthermore, a miRNA inhibitor approach was used to determine the functional role of selected miRNA candidates.

Results

The sequencing of small RNA libraries generated >147 million raw reads in all four libraries and identified a total of >250 miRNAs across the four libraries. We identified 23 and 14 differentially expressed miRNAs in salivary glands, and midgut tissues infected with E. chaffeensis, respectively. Three differentially expressed miRNAs (miR-87, miR-750, and miR-275) were further characterized to determine their roles in pathogen infection. Inhibition of target miRNAs significantly decreased the E. chaffeensis load in tick tissues, which warrants more in-depth mechanistic studies.

Conclusions

The current study identified known and novel miRNAs and suggests that interfering with these miRNAs may impact the vectorial capacity of ticks to harbor Ehrlichia. This study identified several new miRNAs for future analysis of their functions in tick biology and tick-pathogen interaction studies.

Analysis of Amblyomma americanum microRNAs in response to Ehrlichia chaffeensis infection and their potential role in vectorial capacity

Background

MicroRNAs (miRNAs) represent a subset of small noncoding RNAs and carry tremendous potential for regulating gene expression at the post-transcriptional level. They play pivotal roles in distinct cellular mechanisms including inhibition of bacterial, parasitic, and viral infections via immune response pathways. Intriguingly, pathogens have developed strategies to manipulate the host's miRNA profile, fostering environments conducive to successful infection. Therefore, changes in an arthropod host's miRNA profile in response to pathogen invasion could be critical in understanding host-pathogen dynamics. Additionally, this area of study could provide insights into discovering new targets for disease control and prevention. The main objective of the present study is to investigate the functional role of differentially expressed miRNAs upon Ehrlichia chaffeensis, a tick-borne pathogen, infection in tick vector, Amblyomma americanum.

Methods

Small RNA libraries from uninfected and E. chaffeensis-infected Am. americanum midgut and salivary gland tissues were prepared using the Illumina Truseq kit. Small RNA sequencing data was analyzed using miRDeep2 and sRNAtoolbox to identify novel and known miRNAs. The differentially expressed miRNAs were validated using a quantitative PCR assay. Furthermore, a miRNA inhibitor approach was used to determine the functional role of selected miRNA candidates.

Results

The sequencing of small RNA libraries generated >147 million raw reads in all four libraries and identified a total of >250 miRNAs across the four libraries. We identified 23 and 14 differentially expressed miRNAs in salivary glands, and midgut tissues infected with E. chaffeensis, respectively. Three differentially expressed miRNAs (miR-87, miR-750, and miR-275) were further characterized to determine their roles in pathogen infection. Inhibition of target miRNAs significantly decreased the E. chaffeensis load in tick tissues, which warrants more in-depth mechanistic studies.

Conclusions

The current study identified known and novel miRNAs and suggests that interfering with these miRNAs may impact the vectorial capacity of ticks to harbor Ehrlichia. This study identified several new miRNAs for future analysis of their functions in tick biology and tick-pathogen interaction studies.

The Genome-Wide Identification, Characterization, and Expression Analysis of the Strictosidine Synthase-like Family in Maize (Zea mays L.)

Maize is often subjected to various environmental stresses. The strictosidine synthase-like (SSL) family is thought to catalyze the key step in the monoterpene alkaloids synthesis pathway in response to environmental stresses. However, the role of ZmSSL genes in maize growth and development and its response to stresses is unknown. Herein, we undertook the systematic identification and analysis of maize SSL genes. Twenty SSL genes were identified in the maize genome. Except for chromosomes 3, 5, 6, and 10, they were unevenly distributed on the remaining 6 chromosomes. A total of 105 SSL genes from maize, sorghum, rice, Aegilops tauschii, and Arabidopsis were divided into five evolutionary groups, and ZmSSL gene structures and conserved protein motifs in the same group were similar. A collinearity analysis showed that tandem duplication plays an important role in the evolution of the SSL family in maize, and ZmSSL genes share more collinear genes in crops (maize, sorghum, rice, and Ae. tauschii) than in Arabidopsis. Cis-element analysis in the ZmSSL gene promoter region revealed that most genes contained many development and stress response elements. We evaluated the expression levels of ZmSSL genes under normal conditions and stress treatments. ZmSSL4-9 were widely expressed in different tissues and were positively or negatively regulated by heat, cold, and infection stress from Colletotrichum graminicola and Cercospora zeina. Moreover, ZmSSL4 and ZmSSL5 were localized in the chloroplast. Taken together, we provide insight into the evolutionary relationships of the ZmSSL genes, which would be useful to further identify the potential functions of ZmSSLs in maize.

microRNA Expression Dynamics in Culicoides sonorensis Biting Midges Following Blood-Feeding

Culicoides sonorensis midges vector multiple livestock arboviruses, resulting in significant economic losses worldwide. Due to the tight association between virus transmission, blood feeding, and egg development, understanding midge physiology is paramount to limiting pathogen transmission. Previous studies have demonstrated the importance of small non-coding RNAs (ncRNAs), specifically microRNAs (miRNAs), in multiple aspects of vector physiology. These small ncRNAs regulate gene expression at the post-transcriptional level and display differential expression during pathogen infection. Due to the lack of annotated miRNAs in the biting midge and associated expression profiles, we used small RNA-Seq and miRDeep2 analyses to determine the Culicoides miRNAs in whole females and midgut tissues in response to blood feeding. Our analyses revealed 76 miRNAs within C. sonorensis composed of 73 orthologous and three candidate novel miRNAs, as well as conserved miRNA clusters. miRNA conservation suggests an interesting evolutionary relationship between miRNA expression and hematophagy in the infraorder Culicomorpha. We also identified multiple blood meal-regulated and tissue-enriched miRNAs. Lastly, we further identified miRNAs with expression patterns potentially associated with virus infection by probing publicly available datasets. Together, our data provide a foundation for future ncRNA work to untangle the dynamics of gene regulation associated with midge physiology.

Genome-Wide Analysis of Strictosidine Synthase-like Gene Family Revealed Their Response to Biotic/Abiotic Stress in Poplar

The strictosidine synthase-like (SSL) gene family is a small plant immune-regulated gene family that plays a critical role in plant resistance to biotic/abiotic stresses. To date, very little has been reported on the SSL gene in plants. In this study, a total of thirteen SSLs genes were identified from poplar, and these were classified into four subgroups based on multiple sequence alignment and phylogenetic tree analysis, and members of the same subgroup were found to have similar gene structures and motifs. The results of the collinearity analysis showed that poplar SSLs had more collinear genes in the woody plants Salix purpurea and Eucalyptus grandis. The promoter analysis revealed that the promoter region of PtrSSLs contains a large number of biotic/abiotic stress response elements. Subsequently, we examined the expression patterns of PtrSSLs following drought, salt, and leaf blight stress, using RT-qPCR to validate the response of PtrSSLs to biotic/abiotic stresses. In addition, the prediction of transcription factor (TF) regulatory networks identified several TFs, such as ATMYB46, ATMYB15, AGL20, STOP1, ATWRKY65, and so on, that may be induced in the expression of PtrSSLs in response to adversity stress. In conclusion, this study provides a solid basis for a functional analysis of the SSL gene family in response to biotic/abiotic stresses in poplar.

Antiviral RNAi Mechanisms to Arboviruses in Mosquitoes: microRNA Profile of Aedes aegypti and Culex quinquefasciatus from Grenada, West Indies

Mosquito-borne arboviruses, such as dengue virus, West Nile virus, Zika virus and yellow fever virus, impose a tremendous cost on the health of populations around the world. As a result, much effort has gone into the study of the impact of these viruses on human infections. Comparatively less effort, however, has been made to study the way these viruses interact with mosquitoes themselves. As ingested arboviruses infect their midgut and subsequently other tissue, the mosquito mounts a multifaceted innate immune response. RNA interference, the central intracellular antiviral defense mechanism in mosquitoes and other invertebrates can be induced and modulated through outside triggers (small RNAs) and treatments (transgenesis or viral-vector delivery). Accordingly, modulation of this facet of the mosquito’s immune system would thereby suggest a practical strategy for vector control. However, this requires a detailed understanding of mosquitoes’ endogenous small RNAs and their effects on the mosquito and viral proliferation. This paper provides an up-to-date overview of the mosquito’s immune system along with novel data describing miRNA profiles for Aedes aegypti and Culex quinquefasiatus in Grenada, West Indies.

MicroRNAs and other small RNAs in Aedes aegypti saliva and salivary glands following chikungunya virus infection

Mosquito saliva facilitates blood feeding through the anti-haemostatic, anti-inflammatory and immunomodulatory properties of its proteins. However, the potential contribution of non-coding RNAs to host manipulation is still poorly understood. We analysed small RNAs from Aedes aegypti saliva and salivary glands and show here that chikungunya virus-infection triggers both the siRNA and piRNA antiviral pathways with limited effects on miRNA expression profiles. Saliva appears enriched in specific miRNA subsets and its miRNA content is well conserved among mosquitoes and ticks, clearly pointing to a non-random sorting and occurrence. Finally, we provide evidence that miRNAs from Ae. aegypti saliva may target human immune and inflammatory pathways, as indicated by prediction analysis and searching for experimentally validated targets of identical human miRNAs. Overall, we believe these observations convincingly support a scenario where both proteins and miRNAs from mosquito saliva are injected into vertebrates during blood feeding and contribute to the complex vector-host-pathogen interactions.

An Insight Into the microRNA Profile of the Ectoparasitic Mite Varroa destructor (Acari: Varroidae), the Primary Vector of Honey Bee Deformed Wing Virus

The remarkably adaptive mite Varroa destructor is the most important honey bee ectoparasite. Varroa mites are competent vectors of deformed wing virus (DWV), and the Varroa-virus complex is a major determinant of annual honey bee colony mortality and collapse. MicroRNAs (miRNAs) are 22-24 nucleotide non-coding RNAs produced by all plants and animals and some viruses that influence biological processes through post-transcriptional regulation of gene expression. Knowledge of miRNAs and their function in mite biology remains limited. Here we constructed small RNA libraries from male and female V. destructor using Illumina's small RNA-Seq platform. A total of 101,913,208 and 91,904,732 small RNA reads (>18 nucleotides) from male and female mites were analyzed using the miRDeep2 algorithm. A conservative approach predicted 306 miRNAs, 18 of which were upregulated and 13 downregulated in female V. destructor compared with males. Quantitative real-time PCR validated the expression of selected differentially-expressed female Varroa miRNAs. This dataset provides a list of potential miRNA targets involved in regulating vital Varroa biological processes and paves the way for developing strategies to target Varroa and their viruses.

Uncovering the Worldwide Diversity and Evolution of the Virome of the Mosquitoes Aedes aegypti and Aedes albopictus

Aedes aegypti, the yellow fever mosquito, and Aedes albopictus, the Asian tiger mosquito, are the most significant vectors of dengue, Zika, and Chikungunya viruses globally. Studies examining host factors that control arbovirus transmission demonstrate that insect-specific viruses (ISVs) can modulate mosquitoes' susceptibility to arbovirus infection in both in vivo and in vitro co-infection models. While research is ongoing to implicate individual ISVs as proviral or antiviral factors, we have a limited understanding of the composition and diversity of the Aedes virome. To address this gap, we used a meta-analysis approach to uncover virome diversity by analysing ~3000 available RNA sequencing libraries representing a worldwide geographic range for both mosquitoes. We identified ten novel viruses and previously characterised viruses, including mononegaviruses, orthomyxoviruses, negeviruses, and a novel bi-segmented negev-like group. Phylogenetic analysis suggests close relatedness to mosquito viruses implying likely insect host range except for one arbovirus, the multi-segmented Jingmen tick virus (Flaviviridae) in an Italian colony of Ae. albopictus. Individual mosquito transcriptomes revealed remarkable inter-host variation of ISVs within individuals from the same colony and heterogeneity between different laboratory strains. Additionally, we identified striking virus diversity in Wolbachia infected Aedes cell lines. This study expands our understanding of the virome of these important vectors. It provides a resource for further assessing the ecology, evolution, and interaction of ISVs with their mosquito hosts and the arboviruses they transmit.

Development of miRNA-Based Approaches to Explore the Interruption of Mosquito-Borne Disease Transmission

MicroRNA (miRNA or miR)-based approaches to interrupt the transmission of mosquito-borne diseases have been explored since 2005. A review of these studies and areas in which to proceed is needed. In this review, significant progress is reviewed at the level of individual miRNAs, and miRNA diversification and relevant confounders are described in detail. Current miRNA studies in mosquitoes include four steps, namely, identifying miRNAs, validating miRNA-pathogen interactions, exploring action mechanisms, and performing preapplication investigations. Notably, regarding the Plasmodium parasite, mosquito miRNAs generally bind to mosquito immunity- or development-related mRNAs, indirectly regulating Plasmodium infection; However, regarding arboviruses, mosquito miRNAs can bind to the viral genome, directly modifying viral replication. Thus, during explorations of miRNA-based approaches, researchers need select an ideal miRNA for investigation based on the mosquito species, tissue, and mosquito-borne pathogen of interest. Additionally, strategies for miRNA-based approaches differ for arboviruses and protozoan parasites.

Hormonal regulation of microRNA expression dynamics in the gut of the yellow fever mosquito Aedes aegypti

The yellow fever mosquito Aedes aegypti is an obligatory blood feeder and a major arboviral disease vector, evoking severe public health concerns worldwide. In adult female mosquitoes, the gut is critical for blood digestion and pathogen entry. We aimed for a systematic exploration of microRNA expression dynamics in the gut during the gonadotrophic cycle. Small RNA libraries were constructed from female mosquito gut tissues at five time points. Unsupervised hierarchical clustering revealed three expression clusters (early, mid and late) peaking at sequential time points - 24, 48 and 72 h posteclosion. Differentially expressed miRNAs were identified at 24 h post-blood meal (PBM). Depletions of Methoprene-tolerant [Met; the juvenile hormone (JH) receptor] and Ecdysone receptor [EcR; the receptor to 20-hydroxyecdysone (20E)] were performed using dsRNA to these genes to investigate impacts on microRNA expressions. Our results suggest that Met-mediated signalling downregulates miRNA expression from the early cluster and upregulates that from the late cluster. EcR signalling either up- or downregulated miRNA levels at 24 h PBM, indicating a differential effect of this receptor in miRNA gene expression. Furthermore, miR-281, which is the most abundant miRNA in the gut tissue, is induced and repressed by Met- and EcR-mediated signalling, respectively. Systematic depletion using synthetic antagomir and phenotype examinations indicate that miR-281 is obligatory for the normal progression of blood digestion, ovarian development and reproduction. Collectively, this study unveils expression dynamics of microRNAs in the female gut tissue during the gonadotrophic cycle and demonstrates that they are affected by JH and 20E signalling.

Improved reference genome of the arboviral vector Aedes albopictus

BACKGROUND

The Asian tiger mosquito Aedes albopictus is globally expanding and has become the main vector for human arboviruses in Europe. With limited antiviral drugs and vaccines available, vector control is the primary approach to prevent mosquito-borne diseases. A reliable and accurate DNA sequence of the Ae. albopictus genome is essential to develop new approaches that involve genetic manipulation of mosquitoes.

RESULTS

We use long-read sequencing methods and modern scaffolding techniques (PacBio, 10X, and Hi-C) to produce AalbF2, a dramatically improved assembly of the Ae. albopictus genome. AalbF2 reveals widespread viral insertions, novel microRNAs and piRNA clusters, the sex-determining locus, and new immunity genes, and enables genome-wide studies of geographically diverse Ae. albopictus populations and analyses of the developmental and stage-dependent network of expression data. Additionally, we build the first physical map for this species with 75% of the assembled genome anchored to the chromosomes.

CONCLUSION

The AalbF2 genome assembly represents the most up-to-date collective knowledge of the Ae. albopictus genome. These resources represent a foundation to improve understanding of the adaptation potential and the epidemiological relevance of this species and foster the development of innovative control measures.

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.

Assessing the Potential Interactions between Cellular miRNA and Arboviral Genomic RNA in the Yellow Fever Mosquito, Aedes aegypti

Although the role of exogenous small interfering RNA (siRNA) and P-element induced wimpy testis (PIWI)-interacting RNA (piRNA) pathways in mosquito antiviral immunity is increasingly better understood, there is still little knowledge regarding the role of mosquito cellular microRNA (miRNA). Identifying direct interactions between the mosquito miRNAs and the RNA genome of arboviruses and choosing the relevant miRNA candidates to explore resulting antiviral mechanisms are critical. Here, we carried out genomic analyses to identify Aedes aegypti miRNAs that potentially interact with various lineages and genotypes of chikungunya, dengue, and Zika viruses. By using prediction tools with distinct algorithms, several miRNA binding sites were commonly found within different genotypes/and or lineages of each arbovirus. We further analyzed those miRNAs that could target more than one arbovirus, required a low energy threshold to form miRNA-viralRNA (vRNA) complexes, and predicted potential RNA structures using RNAhybrid software. We predicted miRNA candidates that might participate in regulating arboviral replication in Ae. aegypti. Even without any experimental validation, which should be done as a next step, this study can shed further light on the role of miRNA in mosquito innate immunity and targets for future studies.

Small RNA-Seq Analysis Reveals miRNA Expression Dynamics Across Tissues in the Malaria Vector, Anopheles gambiae

Malaria continues to be a major global health problem, where disease transmission is deeply linked to the repeated blood feeding nature of the anautogenous mosquito. Given the tight link between blood feeding and disease transmission, understanding basic biology behind mosquito physiology is a requirement for developing effective vector-borne disease control strategies. In the mosquito, numerous loss of function studies with notable phenotypes demonstrate microRNAs (miRNAs) play significant roles in mosquito physiology. While the field appreciates the importance of a handful of miRNAs, we still need global mosquito tissue miRNA transcriptome studies. To address this need, our goal was to determine the miRNA transcriptome for multiple tissues of the pre-vitellogenic mosquito. To this end, by using small RNA-Seq analysis, we determined miRNA transcriptomes in tissues critical for mosquito reproduction and immunity including (i) fat body-abdominal wall enriched tissues, (ii) midguts, (iii) ovaries, and (iv) remaining tissues comprised of the head and thorax. We found numerous examples of miRNAs exhibiting pan-tissue high- or low- expression, tissue exclusion, and tissue enrichment. We also updated and consolidated the miRNA catalog and provided a detailed genome architecture map for the malaria vector, Anopheles gambiae. This study aims to build a foundation for future research on how miRNAs and potentially other small RNAs regulate mosquito physiology as it relates to vector-borne disease transmission.

MicroRNAs from saliva of anopheline mosquitoes mimic human endogenous miRNAs and may contribute to vector-host-pathogen interactions

During blood feeding haematophagous arthropods inject into their hosts a cocktail of salivary proteins whose main role is to counteract host haemostasis, inflammation and immunity. However, animal body fluids are known to also carry miRNAs. To get insights into saliva and salivary gland miRNA repertoires of the African malaria vector Anopheles coluzzii we used small RNA-Seq and identified 214 miRNAs, including tissue-enriched, sex-biased and putative novel anopheline miRNAs. Noteworthy, miRNAs were asymmetrically distributed between saliva and salivary glands, suggesting that selected miRNAs may be preferentially directed toward mosquito saliva. The evolutionary conservation of a subset of saliva miRNAs in Anopheles and Aedes mosquitoes, and in the tick Ixodes ricinus, supports the idea of a non-random occurrence pointing to their possible physiological role in blood feeding by arthropods. Strikingly, eleven of the most abundant An. coluzzi saliva miRNAs mimicked human miRNAs. Prediction analysis and search for experimentally validated targets indicated that miRNAs from An. coluzzii saliva may act on host mRNAs involved in immune and inflammatory responses. Overall, this study raises the intriguing hypothesis that miRNAs injected into vertebrates with vector saliva may contribute to host manipulation with possible implication for vector-host interaction and pathogen transmission.

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.

Evidence that microRNAs are part of the molecular toolkit regulating adult reproductive diapause in the mosquito, Culex pipiens

For many insects, diapause is the primary mechanism for surviving unfavorable seasons. Some aspects of diapause regulation are well known, but we still lack a mechanistic understanding of molecular mechanisms that control the diapause pathway. Accumulating evidence suggests microRNAs regulate diapause in evolutionarily diverse insect species including flesh flies and moths, and, it is likely that microRNAs regulate multiple characteristics of diapause, including arrested egg follicle development and fat hypertrophy, in females of the Northern house mosquito, Culex pipiens. To investigate microRNA regulation of diapause in this species, we measured changes in egg follicle development and total lipid content over 22 days following adult emergence. We also evaluated changes in the abundance of candidate microRNAs associated with these physical changes during the same time frame. We found egg follicle size and lipid content were nearly the same in diapausing and nondiapausing females on the day of adult emergence, and then diverged over time such that by day 22 diapausing females had significantly smaller egg follicles and higher total lipids than their nondiapausing counterparts. Several microRNAs associated with lipid metabolism in insects, including miR-14-3p, miR-277-3p, and miR-305-5p, were underexpressed in diapausing females compared to nondiapausing females on the day of adult emergence, which suggests microRNA regulation occurs ahead of observed changes in these two features of the diapause phenotype. We also found miR-309-3p, miR-375-3p which stimulate ovarian development in other mosquito species, were underexpressed in diapausing females of Cx. pipiens at times after diapause is fully established and may be responsible for the arrest in ovarian development in this species. Taken together, our results demonstrate that changes in the abundance of some microRNAs is associated with phenotypic changes in diapause Cx. pipiens and suggests this epigenetic mechanism is part of the molecular toolkit regulating diapause.

microRNA profiles and functions in mosquitoes

Mosquitoes are incriminated as vectors for many crippling diseases, including malaria, West Nile fever, Dengue fever, and other neglected tropical diseases (NTDs). microRNAs (miRNAs) can interact with multiple target genes to elicit biological functions in the mosquitoes. However, characterization and function of individual miRNAs and their potential targets have not been fully determined to date. We conducted a systematic review of published literature following PRISMA guidelines. We summarize the information about miRNAs in mosquitoes to better understand their metabolism, development, and responses to microorganisms. Depending on the study, we found that miRNAs were dysregulated in a species-, sex-, stage-, and tissue/organ-specific manner. Aberrant miRNA expressions were observed in development, metabolism, host-pathogen interactions, and insecticide resistance. Of note, many miRNAs were down-regulated upon pathogen infection. The experimental studies have expanded the identification of miRNA target from the 3' untranslated regions (UTRs) of mRNAs of mosquitoes to the 5' UTRs of mRNAs of the virus. In addition, we discuss current trends in mosquito miRNA research and offer suggestions for future studies.

A comparative analysis of corpora allata-corpora cardiaca microRNA repertoires revealed significant changes during mosquito metamorphosis

The corpora allata (CA) are a pair of endocrine glands with neural connections to the brain and close association with another neuroendocrine organ, the corpora cardiaca (CC). The CA from adult female Aedes aegypti mosquitoes synthesize fluctuating levels of juvenile hormone (JH), which have been linked to the ovarian development and are influenced by nutritional signals. In this study, we investigated the potential involvement of microRNAs (miRNAs), a type of small non-coding RNAs, in the regulation of gene expression in CA-CC complexes during mosquito reproductive development, at stages with distinct JH biosynthesis patterns. We analyzed the miRNA repertoires expressed in the CA-CC of pupae, sugar-fed and blood-fed female Ae. aegypti. In total, 156 mature miRNAs were detected in the CA-CC, with 84 displaying significant differences in expression among the three CA-CC developmental stages. There were more miRNAs that were expressed in pupae, and decreased or were absent after adult emergence, when compared with changes between CA-CC of sugar and blood-fed females. Analysis of the genes identified as potential targets for the CA-CC miRNA repertoires classified them into the broad categories of metabolism, information storage and processing, and cellular processes and signaling; with genes involved in cellular processes and signaling representing the largest portion. Among them, the signal-transduction mechanisms and intracellular trafficking, secretion and vesicular transport contained almost 55% of the genes' targets. A substantial number of miRNAs were differentially abundant in the libraries of the three developmental stages, and those changes were much more notable when pupae and adult stages were compared. We detected putative binding sites for some of the most abundant miRNAs on genes encoding JH biosynthetic enzymes and CC neuropeptides. These studies should help us to gain a better understanding of the regulation of CA-CC activity mediated by miRNAs during major developmental stages in mosquitoes.

Characterization and potential role of microRNA in the Chinese dominant malaria mosquito Anopheles sinensis (Diptera: Culicidae) throughout four different life stages

BACKGROUND

microRNAs (miRNAs) are one kind of small non-coding RNAs widely distributed in insects. Many studies have shown that miRNAs play critical roles in development, differentiation, apoptosis, and innate immunity. However, there are a few reports describing miRNAs in Anopheles sinensis, the most common, and one of the dominant malaria mosquito in China. Here, we investigated the global miRNA expression profile across four different developmental stages including embryo, larval, pupal, and adult stages using Illumina Hiseq 2500 sequencing.

RESULTS

In total, 164 miRNAs were obtained out of 107.46 million raw sequencing reads. 99 of them identified as known miRNAs, and the remaining 65 miRNAs were considered as novel. By analyzing the read counts of miRNAs in all developmental stages, 95 miRNAs showed stage-specific expression (q < 0.01 and |log2 (fold change)| > 1) in consecutive stages, indicating that these miRNAs may be involved in critical physiological activity during development. Sixteen miRNAs were identified to be commonly dysregulated throughout four developmental stages. Many miRNAs showed stage-specific expression, such as asi-miR-2943 was exclusively expressed in the embryo stage, and asi-miR-1891 could not be detected in larval stage. The expression of six selected differentially expressed miRNAs identified by qRT-PCR were consistent with our sequencing results. Furthermore, 5296 and 1902 target genes were identified for the dysregulated 68 known and 27 novel miRNAs respectively by combining miRanda and RNAhybrid prediction. GO annotation and KEGG pathway analysis for the predicted genes of dysregulated miRNAs revealed that they might be involved in a broad range of biological processes related with the development, such as membrane, organic substance transport and several key pathways including protein processing in endoplasmic reticulum, propanoate metabolism and folate biosynthesis. Thirty-two key miRNAs were identified by microRNA-gene network analysis.

CONCLUSION

The present study represents the first global characterization of An. sinensis miRNAs in its four developmental stages. The presence and differential expression of An. sinensis miRNAs imply that such miRNAs may play critical roles in An. sinensis life cycle. A better understanding of the functions of these miRNAs will have great implication for the effective control of vector population and therefore interrupting malaria transmission.

Circulating microRNAs in malaria infection: bench to bedside

Severe malaria has a poor prognosis with a morbidity rate of 80% in tropical areas. The early parasite detection is one of the effective means to prevent severe malaria of which specific treatment strategies are limited. Many clinical characteristics and laboratory testings have been used for the early diagnosis and prediction of severe disease. However, a few of these factors could be applied to clinical practice. MicroRNAs (miRNAs) were demonstrated as useful biomarkers in many diseases such as malignant diseases and cardiovascular diseases. Recently it was found that plasma miR-451 and miR-16 were downregulated in malaria infection at parasitic stages or with multi-organ failure involvement. MiR-125b, -27a, -23a, -150, 17-92 and -24 are deregulated in malaria patients with multiple organ failures. Here, the current findings of miRNAs were reviewed in relation to clinical severity of malaria infection and emphasized that miRNAs are potential biomarkers for severe malaria infection.

Comparative genomics shows that viral integrations are abundant and express piRNAs in the arboviral vectors Aedes aegypti and Aedes albopictus

Background

Arthropod-borne viruses (arboviruses) transmitted by mosquito vectors cause many important emerging or resurging infectious diseases in humans including dengue, chikungunya and Zika. Understanding the co-evolutionary processes among viruses and vectors is essential for the development of novel transmission-blocking strategies. Episomal viral DNA fragments are produced from arboviral RNA upon infection of mosquito cells and adults. Additionally, sequences from insect-specific viruses and arboviruses have been found integrated into mosquito genomes.

Results

We used a bioinformatic approach to analyse the presence, abundance, distribution, and transcriptional activity of integrations from 425 non-retroviral viruses, including 133 arboviruses, across the presently available 22 mosquito genome sequences. Large differences in abundance and types of viral integrations were observed in mosquito species from the same region. Viral integrations are unexpectedly abundant in the arboviral vector species Aedes aegypti and Ae. albopictus, in which they are approximately ~10-fold more abundant than in other mosquito species analysed. Additionally, viral integrations are enriched in piRNA clusters of both the Ae. aegypti and Ae. albopictus genomes and, accordingly, they express piRNAs, but not siRNAs.

Conclusions

Differences in the number of viral integrations in the genomes of mosquito species from the same geographic area support the conclusion that integrations of viral sequences is not dependent on viral exposure, but that lineage-specific interactions exist. Viral integrations are abundant in Ae. aegypti and Ae. albopictus, and represent a thus far underappreciated component of their genomes. Additionally, the genome locations of viral integrations and their production of piRNAs indicate a functional link between viral integrations and the piRNA pathway. These results greatly expand the breadth and complexity of small RNA-mediated regulation and suggest a role for viral integrations in antiviral defense in these two mosquito species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3903-3) contains supplementary material, which is available to authorized users.

Transcriptome-wide microRNA and target dynamics in the fat body during the gonadotrophic cycle of Aedes aegypti

The mosquito Aedes aegypti is a major vector of numerous viral diseases, because it requires a blood meal to facilitate egg development. The fat body, a counterpart of mammalian liver and adipose tissues, is the metabolic center, playing a key role in reproduction. Therefore, understanding of regulatory networks controlling its functions is critical, and the role of microRNAs (miRNAs) in the process is largely unknown. We aimed to explore miRNA expression and potential targets in the female fat body of Ae. aegypti, as well as their changes postblood meal (PBM). Small RNA library analysis revealed five unique miRNA patterns sequentially expressed at five sampled time points, likely responding to, and affecting, waves of upstream hormonal signals and gene expression in the same period. To link miRNA identities with downstream targets, transcriptome-wide mRNA 3' UTR interaction sites were experimentally determined at 72 h posteclosion and 24 h PBM through Argonaute 1 cross-linking and immunoprecipitation followed by high-throughput sequencing. Several target sites were validated by means of in vitro luciferase assays with wild-type and mutated 3' UTRs for six miRNA families. With established transgenic lines, consistent results were observed with spatiotemporal knockdown of miR-8 and luciferase assays. We further investigated miRNAs potentially regulating various physiological processes based on Clusters of Orthologous Groups functional categories. Hence, the present work comprehensively elucidated miRNA expression and target dynamics in the female mosquito fat body, providing a solid foundation for future functional studies of miRNA regulation during the gonadotrophic cycle.

Identification and characterization of microRNAs in the screwworm flies Cochliomyia hominivorax and Cochliomyia macellaria (Diptera: Calliphoridae)

MicroRNAs (miRNAs) are small noncoding RNAs that modulate gene expression through post-transcriptional regulation. Here, we report the identification and characterization of miRNAs in two closely related screwworm flies with different feeding habits: Cochliomyia hominivorax and Cochliomyia macellaria. The New World screwworm, C. hominivorax, is an obligatory parasite of warm-blooded vertebrates, whereas the secondary screwworm, C. macellaria, is a free-living organism that feeds on decaying organic matter. Here, the small RNA transcriptomes of adults and third-instar larvae of both species were sequenced. A total of 110 evolutionarily conserved miRNAs were identified, and 10 putative precursor miRNAs (pre-miRNAs) were predicted. The relative expression of six selected miRNAs was further investigated, including miRNAs that are related to reproduction and neural processes in other insects. Mature miRNAs were also characterized across an evolutionary time scale, suggesting that the majority of them have been conserved since the emergence of the Arthropoda [540 million years ago (Ma)], Hexapoda (488 Ma) and Brachycera (195 Ma) lineages. This study is the first report of miRNAs for screwworm flies. We also performed a comparative analysis with the hereby predicted miRNAs from the sheep blowfly, Lucilia cuprina. The results presented may advance our understanding of parasitic habits within Calliphoridae and assist further functional studies in blowflies.

Identification and characterization of microRNAs expressed in antennae of Holotrichia parallela Motschulsky and their possible roles in olfactory regulation

MicroRNAs (miRNAs) are small noncoding RNAs that play posttranscriptional, regulatory roles in various biological processes. However, there has been limited investigation into the potential function of miRNAs in olfaction. The coleopteran Holotrichia parallela is an economically important pest, and miRNAs have been identified in only one coleopteran (Tribolium castaneum). Therefore, this study was conducted to identify miRNAs expressed in the antennae of H. parallela and obtain insights into their possible roles in olfaction. By combining deep sequencing and miRDeep2 software, a total of 99 miRNAs, including 76 conserved miRNAs and 23 novel miRNAs, were identified from H. parallela antennae. The 76 conserved miRNAs belong to 63 families and the other 23 may be species specific or tissue specific. The identified miRNAs have many conserved features of miRNAs. Evaluation of the conservation of the identified miRNA families across different species revealed that most of the families are insect specific. The prediction and annotation of targets suggested that 13 of the identified miRNAs participate in olfactory regulation. Gender differences in antennal expression of nine of the olfactory-related miRNAs were confirmed by quantitative real-time PCR.

Mosquito-specific microRNA-1890 targets the juvenile hormone-regulated serine protease JHA15 in the female mosquito gut

Females of the hematophagous mosquito species require a vertebrate blood meal to supply amino acids and other nutrients necessary for egg development, serving as the driving force for the spread of many vector-borne diseases in humans. Blood digestion utilizes both early and late phase serine proteases (SPs) that are differentially regulated at the transcriptional and post-transcriptional level. To uncover the regulatory complexity of SPs in the female mosquito midgut, we investigated involvement of miRNAs in regulating the juvenile hormone (JH)-controlled chymotrypsin-like SP, JHA15. We identified regulatory regions complementary to the mosquito-specific miRNA, miR-1890, within the 3' UTR of JHA15 mRNA. The level of the JHA15 transcript is highest post eclosion and drastically declines post blood meal (PBM), exhibiting an opposite trend to miR-1890 that peaks at 24 h PBM. Depletion of miR-1890 results in defects in blood digestion, ovary development and egg deposition. JHA15 mRNA and protein levels are elevated in female mosquitoes with miR-1890 inhibition. JHA15 RNA interference in the miR-1890 depletion background alleviates miR-1890 depletion phenotypes. The miR-1890 gene is activated by the 20-hydroxyecdysone pathway that involves the ecdysone receptor and the early genes, E74B and Broad Z2. Our study suggests that miR-1890 controls JHA15 mRNA stability in a stage- and tissue- specific manner.

Small RNA Profiling in Dengue Virus 2-Infected Aedes Mosquito Cells Reveals Viral piRNAs and Novel Host miRNAs

In Aedes mosquitoes, infections with arthropod-borne viruses (arboviruses) trigger or modulate the expression of various classes of viral and host-derived small RNAs, including small interfering RNAs (siRNAs), PIWI interacting RNAs (piRNAs), and microRNAs (miRNAs). Viral siRNAs are at the core of the antiviral RNA interference machinery, one of the key pathways that limit virus replication in invertebrates. Besides siRNAs, Aedes mosquitoes and cells derived from these insects produce arbovirus-derived piRNAs, the best studied examples being viruses from the Togaviridae or Bunyaviridae families. Host miRNAs modulate the expression of a large number of genes and their levels may change in response to viral infections. In addition, some viruses, mostly with a DNA genome, express their own miRNAs to regulate host and viral gene expression. Here, we perform a comprehensive analysis of both viral and host-derived small RNAs in Aedes aegypti Aag2 cells infected with dengue virus 2 (DENV), a member of the Flaviviridae family. Aag2 cells are competent in producing all three types of small RNAs and provide a powerful tool to explore the crosstalk between arboviral infection and the distinct RNA silencing pathways. Interestingly, besides the well-characterized DENV-derived siRNAs, a specific population of viral piRNAs was identified in infected Aag2 cells. Knockdown of Piwi5, Ago3 and, to a lesser extent, Piwi6 results in reduction of vpiRNA levels, providing the first genetic evidence that Aedes PIWI proteins produce DENV-derived small RNAs. In contrast, we do not find convincing evidence for the production of virus-derived miRNAs. Neither do we find that host miRNA expression is strongly changed upon DENV2 infection. Finally, our deep-sequencing analyses detect 30 novel Aedes miRNAs, complementing the repertoire of regulatory small RNAs in this important vector species.

Mosquitoes of the Aedes family transmit many important viruses, including dengue virus, between their vertebrate hosts. In the mosquito, the growth of these viruses is limited by the antiviral RNA interference pathway. Key to this pathway is a class of small non-coding RNAs known as small interfering RNAs (siRNAs). In addition, two related but distinct small RNA pathways known as the microRNA (miRNA) and the PIWI-interacting RNA (piRNA) pathway are implicated in regulating virus replication in mosquitoes. Thus, since small RNAs may critically influence the transmission of dengue virus, we set out to analyze the populations of viral and mosquito small RNAs that are produced in infected Aedes mosquito cells. We found that besides the well-known viral siRNAs, dengue virus-derived piRNAs were produced in these cells and we identified the PIWI proteins that these small RNAs rely on. In addition, we found that viral miRNAs were not expressed from the dengue virus genome and that the levels of mosquito miRNAs were barely changed upon infection. Finally, our data allowed for the identification of novel Aedes miRNAs, complementing the repertoire of these important regulatory RNAs in vector mosquitoes.

MicroRNA evolution, expression, and function during short germband development in Tribolium castaneum

MicroRNAs are well-established players in the development of multicellular animals. Most of our understanding of microRNA function in arthropod development comes from studies in Drosophila. Despite their advantages as model systems, the long germband embryogenesis of fruit flies is an evolutionary derived state restricted to several holometabolous insect lineages. MicroRNA evolution and expression across development in animals exhibiting the ancestral and more widespread short germband mode of embryogenesis has not been characterized. We sequenced small RNA libraries of oocytes and successive intervals covering the embryonic development of the short germband model organism, Tribolium castaneum. We analyzed the evolution and temporal expression of the microRNA complement and sequenced libraries of total RNA to investigate the relationships with microRNA target expression. We show microRNA maternal loading and sequence-specific 3′ end nontemplate oligoadenylation of maternally deposited microRNAs that is conserved between Tribolium and Drosophila. We further uncover large clusters encoding multiple paralogs from several Tribolium-specific microRNA families expressed during a narrow interval of time immediately after the activation of zygotic transcription. These novel microRNAs, together with several early expressed conserved microRNAs, target a significant number of maternally deposited transcripts. Comparison with Drosophila shows that microRNA-mediated maternal transcript targeting is a conserved process in insects, but the number and sequences of microRNAs involved have diverged. The expression of fast-evolving and species-specific microRNAs in the early blastoderm of T. castaneum is consistent with previous findings in Drosophila and shows that the unique permissiveness for microRNA innovation at this stage is a conserved phenomenon.

Dynamic expression of miRNAs across immature and adult stages of the malaria mosquito Anopheles stephensi

Background

MicroRNAs are small non-coding RNAs that are involved in various biological processes including insect development. Anopheles stephensi serves as primary vector of malaria parasite in Asia and exhibits holometabolous life cycle that involves four different stages of development. Regulation and role of mosquito miRNAs during various stages of mosquito development remain largely unknown.

Methods

High throughput small RNA sequencing was employed for identification and profiling of miRNAs across immature and adult stages of malaria vector, which were further validated using Northern hybridization and real time PCR. Target prediction and pathway analysis was carried out to understand the role of regulated miRNAs in insect development. Degradome sequencing was employed to identify cleaved targets of some regulated miRNAs. Loss of function strategy was employed for miR-989 to understand its probable role in female reproductive process.

Results

Small RNA sequencing and data analysis revealed 111 and 14 known and novel miRNAs respectively across all stages of Anopheles stephensi. Nine miRNAs showed gender specific regulation across different stages of mosquito development. Analysis of miRNAs revealed regulation of 24 and 26 miRNAs across different stages of male and female mosquito development respectively. mRNA targets and significant pathways targeted by regulated miRNAs were identified for each stage of mosquito development. Degradome sequencing revealed twenty nine cleaved targets of insect miRNAs. MicroRNA-989 showed significant up-regulation in the adult female as compared to adult male mosquito. Knockdown of miR-989 expression in adult female using miRNA specific antagomir affected targets playing roles in protein binding, proteolysis and nucleic acid binding in ovary tissue of female mosquito post blood feeding.

Conclusions

This is the first comprehensive effort to understand regulation of Anopheles stephensi miRNAs across developmental stages of male and female mosquito. Preliminary role of regulated miRNAs in mosquito development was revealed by target prediction and pathway analysis. MicroRNA-989 emerged to have important roles in adult female mosquitoes showing significant up-regulation which was further studied using miR-989 specific antagomir. This study provides insights into mosquito development and reproductive process and has implications for effective control of mosquito population required for reducing spread of mosquito-borne infectious diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-0772-y) contains supplementary material, which is available to authorized users.