Skip to main content

Advertisement

Identification and characterization of microRNAs expressed in the African malaria vector Anopheles funestus life stages using high throughput sequencing

Article metrics

Abstract

Background

Over the past several years, thousands of microRNAs (miRNAs) have been identified in the genomes of various insects through cloning and sequencing or even by computational prediction. However, the number of miRNAs identified in anopheline species is low and little is known about their role. The mosquito Anopheles funestus is one of the dominant malaria vectors in Africa, which infects and kills millions of people every year. Therefore, small RNA molecules isolated from the four life stages (eggs, larvae, pupae and unfed adult females) of An. funestus were sequenced using next generation sequencing technology.

Results

High throughput sequencing of four replicates in combination with computational analysis identified 107 mature miRNA sequences expressed in the An. funestus mosquito. These include 20 novel miRNAs without sequence identity in any organism and eight miRNAs not previously reported in the Anopheles genus but are known in non-anopheles mosquitoes. Finally, the changes in the expression of miRNAs during the mosquito development were determined and the analysis showed that many miRNAs have stage-specific expression, and are co-transcribed and co-regulated during development.

Conclusions

This study presents the first direct experimental evidence of miRNAs in An. funestus and the first profiling study of miRNA associated with the maturation in this mosquito. Overall, the results indicate that miRNAs play important roles during the growth and development. Silencing such molecules in a specific life stage could decrease the vector population and therefore interrupt malaria transmission.

Background

MicroRNAs (miRNAs) comprise a large family of endogenous, evolutionarily conserved, non-coding RNA that post-transcriptionally regulate gene expression in plants and animals [1, 2]. Many studies have indicated that miRNAs play a pivotal role in most critical biological events [317]. Therefore, identifying and characterization of miRNAs to understand their physiological and pathological roles have become popular research topics.

After the discovery of the first miRNA in the nematode, Caenorhabditis elegans, thousands of miRNAs have been identified [1821]. The major approaches to identifying miRNAs include genetic screening, direct cloning, bioinformatic analysis, and deep sequencing [22, 23]. The majority of known miRNAs have been identified through traditional direct cloning, which is both time-consuming and labour-intensive. Bioinformatic analysis is limited because the majority of computational programs require genomic sequence as a template. Next generation sequencing (NGS) has provided an innovative tool to look into the genome with an unprecedented depth of coverage. It allows for a comprehensive coverage of miRNAs in many species [2433]. As a consequence, this new approach has opened the door to functional genomic analyses of non-model species. It is widely used for profiling miRNAs in populations in various developmental stages, in either normal or diseased states [3450].

Although there are over 30 species of Anopheles which transmit malaria in the world [51], miRNAs have so far only been experimentally identified in three malaria vectors, the African vector (Anopheles gambiae) and in two Asian vectors (Anopheles stephensi and Anopheles anthropophagus) [7, 52, 53]. Among the African vectors, Anopheles funestus is one of the most proficient malaria vectors, mainly because of its remarkable ability to populate a wide range of ecological settings across Africa [5458]. Therefore, experimental identification of miRNAs controlling key genes required for mosquitoes to complete their life-cycle will not only help to better understand the vector biology, but it may uncover novel approaches to control this mosquito. In this context, the main objective of this study was to identify miRNAs expressed in the four main life stages (eggs, larvae, pupae and unfed adult females) of the mosquito An. funestus using high throughput sequencing technology and bioinformatics approaches.

Methods

Mosquito strain and rearing condition

The experimental work was performed using a colony of An. funestus (FUMOZ) that originates from southern Mozambique [59]. The mosquitoes were reared in the insectary of the Vector Control Reference Laborator at the National Institute for Communicable Diseases (NICD), Johannesburg, South Africa since 2000. The insectary is kept at 25 °C, 80% relative humidity with a 12-h day/night lighting regime and 45-min dusk/dawn cycles.

RNA extraction and small RNA sequencing

Total RNA was isolated separately from the four different life cycle stages of An. funestus (eggs, larvae, pupae and unfed adult females) using TRIzol reagent (Invitrogen, USA) according to the manufacturer’s protocol. Four different biological replicates for each life stage for RNA extraction were prepared. In order to obtain a large and broad miRNA transcriptome data set, RNA was extracted from 5 egg patches, 100 fourth instar larvae, 100 pupae, and 100 unfed adult females for each single batch. The quality and quantity of resulting RNA was assessed using a spectrophotometer (Nanodrop Technologies), and quality assessment determined by using Agilent 2100 Bioanalyzer RNA Nano 6000 kit. No rRNA depletion was performed. The RNA extracts from the four life stages were sent to Macrogen Inc (South Korea) for small RNA sequencing. Sequencing libraries were generated from 1 μg each of the total RNA samples by ligation of adapter RNAs at 5′ and 3′ ends, followed by reverse transcription and PCR using Illumina TruSeq Small RNA Sample Preparation protocol (Illumina Inc., USA). The libraries were size-selected for sequencing of RNA fragments of 18–30 nucleotides. Sequencing was performed on an Illumina HiSeq 2000 platform to obtain single-end reads of 50 bases. Each batch sequenced independently.

Sequencing data processing and analysis

Read quality check and filtering

Following sequencing, the quality of the four sequenced libraries were checked using FastQC [60]. Later, all sequencing reads with low quality tags and shorter than 18 nucleotides were removed using FASTX-Toolkit [61]. For all following analysis, all reads from each stage were combined into a single input.

Mapping the reads to the reference genome

All reads were mapped to the An. funestus strain FUMOZ genome [GenBank: KB668221]. The sequence reads were mapped to the genome using miRDeep2 mapper module [62].

Small non-coding RNAs detection

All reads were aligned to the RNA families database (Rfam) version 11 [63] using the BLASTn algorithm allowing for a 2 nucleotide mismatch and e-value lower than 0.01 in order to annotate the small non-coding RNAs present in the libraries.

Identification of miRNA sequences

For identification of miRNA sequences present in the four life stage datasets, the miRDeep2 package was employed [62]. The package scripts detect known or novel miRNAs from deep sequencing data. It looks for the pattern that the miRNA processing machinery leaves in the sequencing data. The most important pattern that miRDeep2 considers are clusters of reads that align along the reference genome that is compatible with the mature miRNA sequence, the loop sequence, and the star sequence structure of the miRNA precursor molecule. If such a pattern is found, miRDeep2 cuts out the potential miRNA precursor sequence from the reference genome and utilizes an RNA folding algorithm (randfold) from the Vienna package [64] to assess if the sequence can be folded into a hairpin structure. Furthermore, the prediction software searches for potential cleavage sites of Drosha and Dicer. The phylogenetic conservation and filtering of other known small non-coding RNA species can be optionally used to improve the predictions. All identified miRNAs are named according to their most similar miRNAs in the miRNA database (miRBase) release 21 (June 2014) [1821] match.

Differential expression of the miRNAs

Differentially expressed miRNAs between two sequential life stages (egg-larva, larva-pupa, pupa-unfed adult female) were determined by log2 fold change >2 for the normalized reads +1. This analysis was computed using gtools package for the R programming language.

Results

Sequence quality of the four libraries

Small RNAs were sequenced in quadruplicate using Illumina sequencing from eggs, larvae, pupae and unfed adult females batches to identify An. funestus miRNAs expressed during development. More than 150 million raw reads were produced from each life stage (Table 1). The length distribution of the reads was between 18 and 30 nucleotides (data not shown). After filtering the impurities and reads of length smaller than 18 nucleotides, 75.92, 71.52, 80.85 and 75.39%, high-quality reads had Phred quality values (PQV) of 20 obtained from the eggs, larvae, pupae, and unfed adult females libraries, respectively (Table 1). The PQV has previously been reported to be an indicator of base call accuracy and therefore sequence quality [65].

Table 1 Summary of small RNA sequencing data analysis for the four life stage libraries of An. funestus

Mapping reads from the four libraries

Mapping reads over the unmasked genome represents an unbiased option, allowing the detection of known and still undiscovered miRNAs. The total number of the reads mapped to An. funestus genome constitutes only 62.61, 63.90, 73.85 and 59.42% of the total high-quality reads from eggs, larvae, pupae, and unfed adult females libraries, respectively (Table 1).

Fig. 1
figure1

Small non-coding RNAs annotated from the four life stage libraries of An. funestus. From aligning the high quality reads to the RNA families database (Rfam) version 11. The total reads can be divided into five categories; miRNAs, tRNAs, rRNAs, others and unmatched. The others and unmatched referred to the other class of small non-coding RNAs and the reads not aligned to RNA families database , respectively

Annotation of small non-coding RNAs in the four libraries

In order to annotate other small non-coding RNAs in the four libraries, all clean reads were aligned against Rfam version 11. As shown in Fig. 1, the most abundant class of small non-coding RNAs in the eggs library were miRNAs. However, rRNAs were most abundant in the larvae and pupae, and tRNAs in the unfed adult females library.

Table 2 Known miRNAs identified in the four life stages libraries of An. funestus

Identification of miRNAs in the four life stages libraries of An. funestus

Of the 65 previously known An. gambiae mature miRNAs in miRBase release 21 (June 2014), 64 mature miRNAs were detected in the sequenced short reads of An. funestus (Table 2). The two miR-276 precursors in An. gambiae generate two different mature miRNAs whereas in An. funestus they produce the same mature miRNA. These include all the 46 miRNA sequences computationally predicted in the An. funestus strain FUMOZ genome hosted by the invertebrate vectors of human pathogens database (VectorBase) [66]. The full precursor structure (mature, loop and star sequence) or parts of it were detected. Only star sequences were found for miR-12, miR-306, miR-965, miR-989, miR-1174 and miR-1889 and the precursor locations for these miRNAs were not determined. Reads that do not match any known An. gambiae miRNA sequences from miRBase and mapped to the An. funestus genome were screened for precursor or hairpin structures. The miRDeep2 algorithm was employed to investigate if non-annotated sequences mapping to the An. funestus genome demonstrated folding properties of precursors or hairpins. A total of 43 mature miRNA candidates expressed from 42 precursor candidates were identified (Table 3). Among these miRNAs, 15 miRNAs have been described in the Anopheles genus (An. gambiae and/or An. anthropophagus) but not registered in miRBase database, eight have not been previously reported in the Anopheles genus but are known in other mosquitoes such as Aedes aegypti and/or Culex quinquefasciatus, and the remaining 20 miRNAs have not been described before in any species. For all the mature miRNA candidates, the full or parts of the precursor were detected in the four life stage libraries except for miR-2943 that was detected only in the eggs library. MiRNAs frequently exhibit sequence differences when compared to a reference mature sequence, generating multiple variations that are known as isoforms. New isoforms were identified for miR-927 (miR-927-3p) and miR-980 (miR-980-5p). The new isoforms of miR-980 were expressed only in the eggs library. Our analysis uncovered a third stem-loop precursor for miR-2.

Table 3 Novel miRNAs identified in the four life stages libraries of An. funestus
Fig. 2
figure2

Heatmaps clustering of miRNAs expressed in the four life stage libraries of An. funestus. The clustering was performed on all known (a) and novel (b) miRNAs based on a raw read copy number (sequencing frequency) and the four life stage samples. Each row represents a stage and each column represents one miRNA. The miRNA clustering is shown on top. The colour scale (shown on the left) illustrates the number of the reads of a miRNA across the life stages

Fig. 3
figure3

Fold-change of the known miRNAs during the development of An. funestus. The fold-changes were calculated for each miRNA (x-axis) using normalized reads +1 (y-axis, bars)

Fig. 4
figure4

Fold-change of the novel miRNAs during the development of An. funestus. The fold-changes were calculated for each miRNA (x-axis) using normalized reads +1 (y-axis, bars)

Expression profiling of miRNAs during An. funestus development

To obtain insight into possible stage dependent roles of miRNAs in An. funestus, the expression patterns of miRNAs in different life stages were examined based on the number of reads obtained. The heatmaps (Fig. 2) summarize the expression of the miRNAs identified in the four life stage libraries. The majority of miRNAs were sequenced between 1–66 times. The expression profiles of miRNAs varied from highly specific to ubiquitous during the four stages. Among the known miRNAs, miR-263a was the most frequently expressed miRNA in the eggs, the larvae, and the pupae libraries (>1 million reads), whereas miR-8 was the dominant miRNA in the unfed adult females library (>300 thousand reads). Nevertheless, miR-87 had the lowest expression level in the egg and the pupa, similarly miR-189 in the larva, and miR-309 and miR-929 in the unfed adult females library. The highly expressed miRNA amongst the new miRNAs in the eggs library was miR-2944a (>1.5 million reads), miR-2765 (>20 thousand reads) in the larvae and pupae libraries and miR-999 (>30 thousand reads) in the unfed adult females library. To analyse the changes in the expression of miRNAs during the mosquito life stages, all read counts within a dataset were normalized, then log2 (fold-changes) were calculated between each two stages. The results of the comparisons (egg-larva, larva-pupa and pupa-unfed adult female) of the expression level of all the miRNAs in the four libraries are shown in Figs. 3 and 4. Between the egg to the larva, 29 miRNAs were up-regulated and 12 other miRNAs were down-regulated. Six miRNAs were up-regulated and another 13 down-regulated between the larva and pupa stage. However, seven miRNAs were up-regulated and four down-regulated between the pupa and the unfed adult female.

Discussion

Although thousands of small RNAs have been identified [1821], the challenge remains to fully identify all small RNAs, especially very low abundance miRNAs and to determine their individual functions. The majority of known miRNAs have been identified through the traditional cloning method [67]. The Illumina sequencing approach is one of high throughput technologies by which miRNAs can be detected in any organism without prior sequence or secondary structure information sequence or secondary structure information. Here, 107 miRNA sequences expressed in the four life stages of the African malaria vector An. funestus were identify and characterize using the deep sequencing approach. This method is more powerful than other conventional technologies previously used in mosquitoes [7, 11, 52], as it is able to identify new miRNAs that are beyond the capabilities of older traditional methods.

In this study, approximately 100 million high-quality reads from each life stage by deep sequencing were obtained. The size distribution of sequenced reads showed peaks between 20–30 nucleotides compare to the average miRNA length of 22 nucleotides in animals [1]. Using similar deep sequencing techniques, other studies on insects show small RNA sequence size distributions with a peak at 22 nucleotides [11, 34, 68, 69]. All the high-quality reads were mapped on to the An. funestus genome with low percentage (average of 65%) identity. Such a mapping bias has been reported in several studies [21, 7072]. This approach has the weakness that it might favour alignment ambiguities due to the limited alignment specificity given by the small length of mature miRNAs (18–25 nucleotides) detected by the short reads, and to the size and high complexity of an unmasked reference genome [70].

As expected, most of the known miRNAs identified in An. funestus were highly conserved across diverse animal species, suggesting that the ancient regulatory pathways mediated by evolutionary conserved miRNAs are present in mosquitoes. The full stem-loop precursor structure was detected in all four life stage libraries for miR-8, miR-9c, miR-100, miR-275, miR-277, miR-317, miR-1000 and miR-1890. In the majority of cases, mature miRNAs are more abundant than loop and star sequences. Additionally, miR-309 and miR-988 were not detected in the pupae and unfed adult females libraries, miR-iab-4 was not found in eggs and pupae libraries and miR-308 was not identified from the larvae library. These results indicate that the expression of some miRNAs are probably stage specific. It can be speculated that a miRNA may be involved in regulation of function and dysfunction, differentiation, growth and development of a specific stage [73].

The identification of novel miRNAs is an eminent and challenging problem for the understanding of post-transcriptional gene regulation. The characteristic hairpin structure of miRNA precursors can be used to predict novel miRNAs. With this feature and using miRDeep2, novel miRNAs were predicted by exploring the secondary structure, the Dicer cleavage site, and the minimum free energy of the unannotated small RNA tags that could be mapped to the An. funestus genome. Based on the analysis, 43 mature sequences were identified as novel miRNA candidates because they were not captured in any RNA database. The full precursor structure for some of these novel miRNAs was identified. Detection of the miRNA star is a strong clue, albeit not infallible, for the formation of precursor hairpin structures. This adds weight to the authenticity of the predicted candidates [74, 75]. However, the evolution and function of the star miRNAs remains unclear. Two studies proposed that these star miRNAs might differ from their sense partners by acting on different mRNA targets [76, 77].

The miRBase databases release 21 (June 2014) was searched for homologs to determine whether these novel miRNAs are conserved among other animal species. Some candidates are conserved in other insect species but not in the Anopheles genus, suggesting that these are insect-specific miRNAs. Among the newly identified Anopheles miRNAs, four mosquito specific miRNAs (miR-2940, miR-2941, miR-2942 and miR-2943) were detected [68]. This is the first description of miR-2941 and miR-2943 in the anopheline species using the high-throughput sequencing technology.

The NGS technology was sensitive enough to detect a new variant for miR-927 and miR-980, miRNAs found only in insects. These results are congruent with the existence of multiple variants for both miRNAs in other insects such as Drosophila [18]. A third stem-loop precursor for miR-2 was also detected which produces different mature miR-2. The miR-2 family is widespread in invertebrates, and it is the largest family of miRNAs in the model species Drosophila melanogaster (8 family members), Capitella teleta (7 family members), Apis mellifera (6 family members), Bombyx mori (5 members) [20].

Counting redundant miRNA reads revealed that expression varies significantly among different miRNAs. In the three pre-adult stages, insect-specific miRNA (miR-263a) was found to be the most abundant miRNA. However, miR-8 was the common miRNA in the unfed adult females library. Furthermore, the four libraries shared five out of the top ten most frequently occurring miRNAs: miR-263a, miR-10, miR-184, bantam and miR-8. In Drosophila, miRNA miR-263a confers robustness during development by protecting nascent sense organs from apoptosis [9]. Moreover, a recent study on Nilaparvata lugens reported that miR-263a was found as high abundant miRNAs in the last instar female nymph females [78]. Both miR-8 and miR-184 were reported in the embryos of Drosophila [79, 80], mosquitoes [68, 81], silk worm [82], Schistosoma japonicum [83], zebrafish [84], Asian seabass[85], mouse [8688] and humans [15, 89]. This suggest a conservative developmental function for these two miRNAs across different animals. High number of reads for miR-2941, miR-2943 miR-2944a, and miR-2944b in the eggs library were also observed, which indicate an embryonic roles for these insect specific miRNAs. High expression for miR-2941 and miR-2943 in the embryonic stage was reported previously [53, 90].

The expression of miRNAs varies across different developmental stages [82, 91]. In this study, bantam, miR-8, miR-31 and miR-278 were among the up-regulated miRNAs between the egg and larva stage. In D. melanogaster, overexpression of bantam induces tissue overgrowth. Related to growth, and also in D. melanogaster, miR-8 and miR-278 have been implicated in insulin receptor signaling , thus contributing to regulation of the energy balance [6]. Embryos depleted for miR-31 can complete development, but were affected by severe segmentation defects [92]. In the gregarious phase of locust, canonical miRNAs were expressed at levels between 1.5- and 2-fold higher than in the solitary phase. The most prominent differences were found in miR-276, miR-125, and miR-315 [34]. Interestingly, the same change in these miRNAs between the egg and the larva stage except for miR-315 were observed. As is known, the genes that encode for miRNA are distributed across the chromosome either individually, or in clusters in which two or more miRNA genes are located within a short distance on the same segment of a chromosome [93, 94]. Therefore, it is assumed that miRNA genes located in a gene cluster are first transcribed as a single primary transcript that is subsequently processed to generate the individual miRNAs [7]. In the larvae library, very high relative expression levels of the let-7-complex locus (miR-125, let-7 and miR-100), the cluster of the two mosquito-specifc miRNAs; miR-1174 and miR-1175, miR-278 and miR-307, miR-305 and miR-275, miR-210 and miR-927 were observed. Down-regulation of the miR-92 cluster (miR-92a and miR-92b) and miR-309 cluster (miR-309 and miR-286) were also observed. Similarly in Drosophila, miR-309 and miR-286 (which was processed from a single 1.5 kb primary transcript) displayed a dynamic pattern of expression in the early embryo [95].

Between the larva and the pupa, miR-193, miR-277 and miR-1890 were significantly up-regulated. miR-34, miR-308, miR-375, miR-1175 and and miR-2942 were down-regulated in the pupa after increasing in the larva stage. Again, down-regulation of the miR-92a cluster was noticed. The expression of the members of this cluster has been related to the embryonic development in Ae. aegypti [68] and B. mori [69].

Among the up-regulated miRNA between the pupa and the unfed adult female is miR-1891 the mosquito-specific miRNA which displayed changes in its expression levels in the unfed adult females library, suggesting a significant role for this miRNA during development beside its function in the host response to infection [14]. In the unfed adult females library, up-regulation of miR-34 was noticed compared to significant down-regulation in the pupae library, and a major change in miR-989 expression for the first time. The expression of these miRNAs has been studied in different fly species, and results have shown that miR-989 expression is restricted to females, and predominantly to the ovary in An. anthropophagus, An. stephensi, Ae. aegypti and Drosophila [52, 53, 68, 96], although it was later detected in the midgut of An. gambiae [7]. Both miRNAs (miR-34 and miR-989) with two other miRNAs displayed changes in the expression levels during Plasmodium parasite invasion [7]. For unknown reasons, some miRNAs such as miR-87 and the arthropod-specific miR-929 showed low read counts in the four life stage libraries, suggesting that miRNAs might be involved in other process but not development. Further studies are needed to reveal the function of these miRNAs.

The predicted novel miRNAs exhibited much lower expression levels, consistent with the evidence that non-conserved miRNAs are often expressed at a lower level than conserved miRNAs [9799].

It is important to note that the expression profile analysis in this study are based only on the sequence read counts and required further experimental validation. However, this result is a key step towards improving our understanding of the complexity and regulation mode of miRNAs in mosquitoes. Changes in the expression profiles were noted in all stages indicating a role for these small RNAs in the mosquito maturation. Knockdown or blocking the biogenesis pathway of one of these miRNAs may limit the mosquito’s development at a crucial stage thereby leading to novel approaches to combat this mosquito in the early development stages.

Conclusion

In this study, 107 mature miRNA sequences in the four developing stages (eggs, larvae, pupae, and unfed adult females) of An. funestus were identified, one of the most prevalent malaria vector on the African continent using the high throughput sequencing and described their expression patterns during development.

References

  1. 1.

    Bartel DP. MicroRNAs: genomics, biogenesis, mechanism and function. Cell. 2004;116:281–97.

  2. 2.

    Krol J, Loedige I, Filipowicz W. The widespread regulation of microRNA biogenesis, function and decay. Nat Rev Genet. 2010;11:597–610.

  3. 3.

    Macdonald PM, Struhl G. A molecular gradient in early Drosophila embryos and its role in specifying the body pattern. Nature. 1986;324:537–45.

  4. 4.

    Hutvágner G, Zamore PD. A microRNA in a multiple-turnover RNAi enzyme complex. Science. 2002;297:2056–60.

  5. 5.

    Brennecke J, Hipfner DR, Stark A, Russell RB, Cohen SM. Bantam encodes a developmentally regulated microRNA that controls cell proliferation and regulates the proapoptotic gene hid in Drosophila. Cell. 2003;113:25–36.

  6. 6.

    Teleman AA, Maitra S, Cohen SM. Drosophila lacking microRNA miR-278 are defective in energy homeostasis. Genes Dev. 2006;20:417–22.

  7. 7.

    Winter F, Edaye S, Hüttenhofer A, Brunel C. Anopheles gambiae miRNAs as actors of defence reaction against Plasmodium invasion. Nucleic Acids Res. 2007;35:6953–62.

  8. 8.

    Zhang Y, Zhou X, Ge X, Jiang J, Li M, Jia S, et al. Insect-specific microRNA involved in the development of the silkworm Bombyx mori. PLoS ONE. 2009;4:e4677.

  9. 9.

    Hilgers V, Bushati N, Cohen SM. Drosophila microRNAs 263a/b confer robustness during development by protecting nascent sense organs from apoptosis. PLoS Biol. 2010;8:e1000396.

  10. 10.

    Hussain M, Asgari S. Functional analysis of a cellular microRNA in insect host-ascovirus interaction. J Virol. 2010;84:612–20.

  11. 11.

    Skalsky RL, Vanlandingham DL, Scholle F, Higgs S, Cullen BR. Identification of microRNAs expressed in two mosquito vectors, Aedes albopictus and Culex quinquefasciatus. BMC Genom. 2010;11:119.

  12. 12.

    Zeiner GM, Norman KL, Thomson JM, Hammond SM, Boothroyd JC. Toxoplasma gondii infection specifically increases the levels of key host microRNAs. PLoS ONE. 2010;5:e8742.

  13. 13.

    Dkhil M, Abdel-Baki AA, Delić D, Wunderlich F, Sies H, Al-Quraishy S. Eimeria papillata: upregulation of specific miRNA-species in the mouse jejunum. Exp Parasitol. 2011;127:581–6.

  14. 14.

    Hussain M, Frentiu FD, Moreira LA, O’Neill SL, Asgari S. Wolbachia uses host microRNAs to manipulate host gene expression and facilitate colonization of the dengue vector Aedes aegypti. Proc Natl Acad Sci USA. 2011;108:9250–5.

  15. 15.

    Vallejo DM, Caparros E, Dominguez M. Targeting notch signalling by the conserved miR-8/200 microRNA family in development and cancer cells. EMBO J. 2011;30:756–69.

  16. 16.

    Marco A, Hooks K, Griffiths-Jones S. Evolution and function of the extended miR-2 microRNA family. RNA Biol. 2012;9:242–8.

  17. 17.

    Choi IK, Hyun S. Conserved microRNA miR-8 in fat body regulates innate immune homeostasis in Drosophila. Dev Comp Immunol. 2012;37:50–4.

  18. 18.

    Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ. miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006;34:D140–4.

  19. 19.

    Griffiths-Jones S. miRBase: the microRNA sequence database. Methods Mol Biol. 2006;342:129–38.

  20. 20.

    Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ. miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008;36:D154–8.

  21. 21.

    Kozomara A, Griffiths-Jones S. miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011;39:D152–7.

  22. 22.

    Wang WC, Lin FM, Chang WC, Lin KY, Huang HD, Lin NS. miRExpress: analyzing high-throughput sequencing data for profiling microRNA expression. BMC Bioinform. 2009;10:328.

  23. 23.

    Metzker ML. Sequencing technologies—the next generation. Nat Rev Genet. 2010;11:31–46.

  24. 24.

    Creighton CJ, Benham AL, Zhu H, Khan MF, Reid JG, Nagaraja AK, et al. Discovery of novel microRNAs in female reproductive tract using next generation sequencing. PLoS ONE. 2010;5:e9637.

  25. 25.

    Jung CH, Hansen MA, Makunin IV, Korbie DJ, Mattick JS. Identification of novel non-coding RNAs using profiles of short sequence reads from next generation sequencing data. BMC Genom. 2010;11:77.

  26. 26.

    Hackl M, Jakobi T, Blom J, Doppmeier D, Brinkrolf K, Szczepanowski R, et al. Next-generation sequencing of the Chinese hamster ovary microRNA transcriptome: identification, annotation and profiling of microRNAs as targets for cellular engineering. J Biotechnol. 2011;153:62–75.

  27. 27.

    Havecker ER. Detection of small RNAs and microRNAs using deep sequencing technology. Methods Mol Biol. 2011;732:55–68.

  28. 28.

    Keller A, Backes C, Leidinger P, Kefer N, Boisguerin V, Barbacioru C, et al. Next-generation sequencing identifies novel microRNAs in peripheral blood of lung cancer patients. Mol Biosyst. 2011;7:3187–99.

  29. 29.

    Kong BW. Identification of virus encoding microRNAs using 454 FLX sequencing platform. Methods Mol Biol. 2011;733:81–91.

  30. 30.

    Persson H, Kvist A, Rego N, Staaf J, Vallon-Christersson J, Luts L, et al. Identification of new microRNAs in paired normal and tumor breast tissue suggests a dual role for the ERBB2/Her2 gene. Cancer Res. 2011;71:78–86.

  31. 31.

    Gébelin V, Argout X, Engchuan W, Pitollat B, Duan C, Montoro P, et al. Identification of novel microRNAs in Hevea brasiliensis and computational prediction of their targets. BMC Plant Biol. 2012;12:18.

  32. 32.

    Gunaratne PH, Coarfa C, Soibam B, Tandon A. miRNA data analysis: next-gen sequencing. Methods Mol Biol. 2012;822:273–88.

  33. 33.

    Pérez-Quintero ÁL, Quintero A, Urrego O, Vanegas P, López C. Bioinformatic identification of cassava miRNAs differentially expressed in response to infection by Xanthomonas axonopodis pv. manihotis. BMC Plant Biol. 2012;12:29.

  34. 34.

    Wei Y, Chen S, Yang P, Ma Z, Kang L. Characterization and comparative profiling of the small RNA transcriptomes in two phases of locust. Genome Biol. 2009;10:R6.

  35. 35.

    Buchold GM, Coarfa C, Kim J, Milosavljevic A, Gunaratne PH, Matzuk MM. Analysis of microRNA expression in the prepubertal testis. PLoS ONE. 2010;5:e15317.

  36. 36.

    Nikopoulos K, Gilissen C, Hoischen A, van Nouhuys CE, Boonstra FN, Blokland EAW, et al. Next-generation sequencing of a 40 Mb linkage interval reveals TSPAN12 mutations in patients with familial exudative vitreoretinopathy. Am J Hum Genet. 2010;86:240–7.

  37. 37.

    Borges F, Pereira PA, Slotkin RK, Martienssen RA, Becker JD. MicroRNA activity in the Arabidopsis male germline. J Exp Bot. 2011;62:1611–20.

  38. 38.

    Krawetz SA, Kruger A, Lalancette C, Tagett R, Anton E, Draghici S, et al. A survey of small RNAs in human sperm. Hum Reprod. 2011;26:3401–12.

  39. 39.

    Mohorianu I, Schwach F, Jing R, Lopez-Gomollon S, Moxon S, Szittya G, et al. Profiling of short RNAs during fleshy fruit development reveals stage-specific sRNAome expression patterns. Plant J. 2011;67:232–46.

  40. 40.

    Wu Q, Lu Z, Li H, Lu J, Guo L, Ge Q. Next-generation sequencing of microRNAs for breast cancer detection. J Biomed Biotechnol. 2011;2011:597145.

  41. 41.

    Boeri M, Pastorino U, Sozzi G. Role of microRNAs in lung cancer: microRNA signatures in cancer prognosis. Cancer J. 2012;18:268–74.

  42. 42.

    Cutting AD, Bannister SC, Doran TJ, Sinclair AH, Tizard MVL, Smith CA. The potential role of microRNAs in regulating gonadal sex differentiation in the chicken embryo. Chromosome Res. 2012;20:201–13.

  43. 43.

    Gilabert-Estelles J, Braza-Boils A, Ramon LA, Zorio E, Medina P, Espana F, et al. Role of microRNAs in gynecological pathology. Curr Med Chem. 2012;19:2406–13.

  44. 44.

    Kang M, Zhao Q, Zhu D, Yu J. Characterization of microRNAs expression during maize seed development. BMC Genom. 2012;13:360.

  45. 45.

    Leidner RS, Ravi L, Leahy P, Chen Y, Bednarchik B, Streppel M, et al. The microRNAs, MiR-31 and MiR-375, as candidate markers in Barrett’s esophageal carcinogenesis. Genes Chromosomes Cancer. 2012;51:473–9.

  46. 46.

    Liu F, Peng W, Li Z, Li W, Li L, Pan J, et al. Next-generation small RNA sequencing for microRNAs profiling in Apis mellifera: comparison between nurses and foragers. Insect Mol Biol. 2012;21:297–303.

  47. 47.

    Nikaki A, Piperi C, Papavassiliou AG. Role of microRNAs in gliomagenesis: targeting miRNAs in glioblastoma multiforme therapy. Expert Opin Investig Drugs. 2012;21:1475–88.

  48. 48.

    Shamimuzzaman M, Vodkin L. Identification of soybean seed developmental stage-specific and tissue-specific miRNA targets by degradome sequencing. BMC Genom. 2012;13:310.

  49. 49.

    Wei C, Salichos L, Wittgrove CM, Rokas A, Patton JG. Transcriptome-wide analysis of small RNA expression in early zebrafish development. RNA. 2012;18:915–29.

  50. 50.

    Yao MJ, Chen G, Zhao PP, Lu MH, Jian J, Liu MF, et al. Transcriptome analysis of microRNAs in developing cerebral cortex of rat. BMC Genom. 2012;13:232.

  51. 51.

    Kiszewski A, Mellinger A, Spielman A, Malaney P, Sachs SE, Sachs J. A global index representing the stability of malaria transmission. Am J Trop Med Hyg. 2004;70:486–98.

  52. 52.

    Mead EA, Tu Z. Cloning, characterization, and expression of microRNAs from the Asian malaria mosquito Anopheles stephensi. BMC Genom. 2008;9:244.

  53. 53.

    Liu W, Huang H, Xing C, Li C, Tan F, Liang S. Identification and characterization of the expression profile of microRNAs in Anopheles anthropophagus. Parasit Vectors. 2014;7:159.

  54. 54.

    Shiff CJ, Minjas JN, Hall T, Hunt RH, Lyimo S, Davis JR. Malaria infection potential of anopheline mosquitoes sampled by light trapping indoors in coastal Tanzanian villages. Med Vet Entomol. 1995;9:256–62.

  55. 55.

    Hargreaves K, Koekemoer LL, Brooke BD, Hunt RH, Mthembu J, Coetzee M. Anopheles funestus resistant to pyrethroid insecticides in South Africa. Med Vet Entomol. 2000;14:181–9.

  56. 56.

    Kamau L, Koekemoer LL, Hunt RH, Coetzee M. Anopheles parensis: the main member of the Anopheles funestus species group found resting inside human dwellings in Mwea area of central Kenya toward the end of the rainy season. J Am Mosq Control Assoc. 2003;19:130–3.

  57. 57.

    Coetzee M, Fontenille D. Advances in the study of Anopheles funestus, a major vector of malaria in Africa. Insect Biochem Mol Biol. 2004;34:599–605.

  58. 58.

    Ayala D, Caro-Riaño H, Dujardin JP, Rahola N, Simard F, Fontenille D. Chromosomal and environmental determinants of morphometric variation in natural populations of the malaria vector Anopheles funestus in Cameroon. Infect Genet Evol. 2011;11:940–7.

  59. 59.

    Hunt RH, Brooke BD, Pillay C, Koekemoer LL, Coetzee M. Laboratory selection for and characteristics of pyrethroid resistance in the malaria vector Anopheles funestus. Med Vet Entomol. 2005;19:271–5.

  60. 60.

    FastQC. http://www.bioinformatics.babraham.ac.uk/projects/fastqc.

  61. 61.

    FASTX-Toolkit. http://hannonlab.cshl.edu/fastx_toolkit/.

  62. 62.

    Friedländer MR, Mackowiak SD, Li N, Chen W, Rajewsky N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012;40:37–52.

  63. 63.

    Gardner PP, Daub J, Tate J, Moore BL, Osuch IH, Griffiths-Jones S. Rfam: Wikipedia, clans and the “decimal” release. Nucleic Acids Res. 2011;39:D141–5.

  64. 64.

    Hofacker IL, Fontana W, Stadler PF, Bonhoeffer S, Tacker M, Schuster P. Fast folding and comparison of RNA secondary structures. Monatshefte f Chemie. 1994;125:167–88.

  65. 65.

    Ewing B, Hillier L, Wendl MC, Green P. Base-calling of automated sequencer traces using phred I. accuracy assessment. Genome Res. 1998;8:175–85.

  66. 66.

    The invertebrate vectors of human pathogens database. https://www.vectorbase.org.

  67. 67.

    Tian G, Yin X, Luo H, Xu X, Bolund L, Zhang X, et al. Sequencing bias: comparison of different protocols of microRNA library construction. BMC Biotechnol. 2010;10:64.

  68. 68.

    Li S, Mead EA, Liang S, Tu Z. Direct sequencing and expression analysis of a large number of miRNAs in Aedes aegypti and a multi-species survey of novel mosquito miRNAs. BMC Genomics. 2009;10:581.

  69. 69.

    Liu S, Li D, Li Q, Zhao P, Xiang Z, Xia Q. MicroRNAs of Bombyx mori identified by Solexa sequencing. BMC Genom. 2010;11:148.

  70. 70.

    Cordero F, Beccuti M, Arigoni M, Donatelli S, Calogero RA. Optimizing a massive parallel sequencing workflow for quantitative miRNA expression analysis. PLoS ONE. 2012;7:e31630.

  71. 71.

    Oshlack A, Wakefield MJ. Transcript length bias in RNA-seq data confounds systems biology. Biol Direct. 2009;4:14.

  72. 72.

    Pelaez P, Trejo MS, Iniguez LP, Estrada-Navarrete G, Covarrubias AA, Reyes JL, et al. Identification and characterization of microRNAs in Phaseolus vulgaris by high-throughput sequencing. BMC Genom. 2012;13:83.

  73. 73.

    Tang Y, Liu D, Zhang L, Ingvarsson S, Chen H. Quantitative analysis of miRNA expression in seven human foetal and adult organs. PLoS ONE. 2011;6:e28730.

  74. 74.

    Fahlgren N, Howell MD, Kasschau KD, Chapman EJ, Sullivan CM, Cumbie JS, et al. High-throughput sequencing of Arabidopsis microRNAs: evidence for frequent birth and death of miRNA genes. PLoS ONE. 2007;2:e219.

  75. 75.

    Sunkar R, Zhou X, Zheng Y, Zhang W, Zhu JK. Identification of novel and candidate miRNAs in rice by high throughput sequencing. BMC Plant Biol. 2008;8:25.

  76. 76.

    Zhang B, Pan X, Stellwag EJ. Identification of soybean microRNAs and their targets. Planta. 2008;229:161–82.

  77. 77.

    Chi X, Yang Q, Chen X, Wang J, Pan L, Chen M, et al. Identification and characterization of microRNAs from peanut (Arachis hypogaea L.) by high-throughput sequencing. PLoS ONE. 2011;6:e27530.

  78. 78.

    Chen Q, Lu L, Hua H, Zhou F, Lu L, Lin Y. Characterization and comparative analysis of small RNAs in three small RNA libraries of the brown planthopper (Nilaparvata lugens). PLoS ONE. 2012;7:e32860.

  79. 79.

    Li P, Peng J, Hu J, Xu Z, Xie W, Yuan L. Localized expression pattern of miR-184 in Drosophila. Mol Biol Rep. 2011;38:355–8.

  80. 80.

    Iovino N, Pane A, Gaul U. miR-184 has multiple roles in Drosophila female germline development. Dev Cell. 2009;17:123–33.

  81. 81.

    Zheng Pm Wu, Jb Jy, Gu, Zj Tu, Xg Chen. Isolation, identification and analysis of the expression profile of miRNAs in Aedes albopictus. Nan Fang Yi Ke Da Xue Xue Bao. 2010;30:677–80.

  82. 82.

    Liu S, Zhang L, Li Q, Zhao P, Duan J, Cheng D, et al. MicroRNA expression profiling during the life cycle of the silkworm (Bombyx mori). BMC Genom. 2009;10:455.

  83. 83.

    Huang J, Hao P, Chen H, Hu W, Yan Q, Liu F, et al. Genome-wide identification of Schistosoma japonicum microRNAs using a deep-sequencing approach. PLoS ONE. 2009;4:e8206.

  84. 84.

    Flynt AS, Thatcher EJ, Burkewitz K, Li N, Liu Y, Patton JG. miR-8 microRNAs regulate the response to osmotic stress in zebrafish embryos. J Cell Biol. 2009;185:115–27.

  85. 85.

    Xia JH, He XP, Bai ZY, Yue GH. Identification and characterization of 63 MicroRNAs in the Asian seabass Lates calcarifer. PLoS ONE. 2011;6:e17537.

  86. 86.

    Wenguang Z, Jianghong W, Jinquan L, Yashizawa M. A subset of skin-expressed microRNAs with possible roles in goat and sheep hair growth based on expression profiling of mammalian microRNAs. OMICS. 2007;11:385–96.

  87. 87.

    Juhila J, Sipilä T, Icay K, Nicorici D, Ellonen P, Kallio A, et al. MicroRNA expression profiling reveals miRNA families regulating specific biological pathways in mouse frontal cortex and hippocampus. PLoS ONE. 2011;6:e21495.

  88. 88.

    Wu J, Bao J, Wang L, Hu Y, Xu C. MicroRNA-184 downregulates nuclear receptor corepressor 2 in mouse spermatogenesis. BMC Dev Biol. 2011;11:64.

  89. 89.

    Hyun S, Lee JH, Jin H, Nam J, Namkoong B, Lee G, et al. Conserved MicroRNA miR-8/miR-200 and its target USH/FOG2 control growth by regulating PI3K. Cell. 2009;139:1096–108.

  90. 90.

    Gu J, Hu W, Wu J, Zheng P, Chen M, James AA, et al. miRNA genes of an invasive vector mosquito Aedes albopictus. PLoS ONE. 2013;8:e67638.

  91. 91.

    Xu H, Wang X, Du Z, Li N. Identification of microRNAs from different tissues of chicken embryo and adult chicken. FEBS Lett. 2006;580:3610–6.

  92. 92.

    Leaman D, Chen PY, Fak J, Yalcin A, Pearce M, Unnerstall U, et al. Antisense-mediated depletion reveals essential and specific functions of microRNAs in Drosophila development. Cell. 2005;121:1097–108.

  93. 93.

    Altuvia Y, Landgraf P, Lithwick G, Elefant N, Pfeffer S, Aravin A, et al. Clustering and conservation patterns of human microRNAs. Nucleic Acids Res. 2005;33:2697–706.

  94. 94.

    Zhang X, Azhar G, Wei JY. The expression of microRNA and microRNA clusters in the aging heart. PLoS ONE. 2012;7:e34688.

  95. 95.

    Biemar F, Zinzen R, Ronshaugen M, Sementchenko V, Manak JR, Levine MS. Spatial regulation of microRNA gene expression in the Drosophila embryo. Proc Natl Acad Sci USA. 2005;102:15907–11.

  96. 96.

    Kugler JM, Verma P, Chen YW, Weng R, Cohen SM. miR-989 is required for border cell migration in the Drosophila ovary. PLoS ONE. 2013;8:e67075.

  97. 97.

    Ling KH, Brautigan PJ, Hahn CN, Daish T, Rayner JR, Cheah PS, et al. Deep sequencing analysis of the developing mouse brain reveals a novel microRNA. BMC Genom. 2011;12:176.

  98. 98.

    Inukai S, de Lencastre A, Turner M, Slack F. Novel MicroRNAs differentially expressed during aging in the mouse brain. PLoS ONE. 2012;7:e40028.

  99. 99.

    Zhang JZ, Ai XY, Guo WW, Peng SA, Deng XX, Hu CG. Identification of miRNAs and their target genes using deep sequencing and degradome analysis in trifoliate orange [Poncirus trifoliate (L.) Raf]. Mol Biotechnol. 2012;51:44–57.

  100. 100.

    Castellano L, Rizzi E, Krell J, Di Cristina M, Galizi R, Mori A, et al. The germline of the malaria mosquito produces abundant miRNAs, endo-siRNAs, piRNAs and 29-nt small RNAs. BMC Genom. 2015;16:100.

Download references

Authors’ contributions

MA conducted the experiments and data analyses, interpreted results, and drafted the first version of the manuscript. BLS and HA reared the mosquitoes and carried out the RNA isolation. DM performed the statistical analysis. LLK provided technical support, supervised the experiments and contributed to the editing of the manuscript. AC conceived the project, coordinated the study, provided funding for the study and helped with revision of the manuscript. All authors read and commented on the final manuscript. All authors read and approved the final manuscript.

Acknowledgements

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Availability of data and material

The sequencing reads are available in ftp://ftp.sanbi.ac.za/anopheles funestus_miRNA_sequencing_reads. All miRNA sequences identified in this study are provided in the manuscript.

Funding

This work was supported by the South African Research Chairs Initiative of the Department of Science and Technology and National Research Foundation of South Africa.

Author information

Correspondence to Alan Christoffels.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Allam, M., Spillings, B.L., Abdalla, H. et al. Identification and characterization of microRNAs expressed in the African malaria vector Anopheles funestus life stages using high throughput sequencing. Malar J 15, 542 (2016) doi:10.1186/s12936-016-1591-0

Download citation

Keywords

  • MicroRNA
  • Non-coding RNA
  • Anopheles funestus
  • Vector-borne diseases
  • Malaria

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.