Cloning, characterization, and expression of microRNAs from the Asian malaria mosquito, Anopheles stephensi
© Mead and Tu; licensee BioMed Central Ltd. 2008
Received: 16 October 2007
Accepted: 23 May 2008
Published: 23 May 2008
microRNAs (miRNAs) are non-coding RNAs that are now recognized as a major class of gene-regulating molecules widely distributed in metozoans and plants. miRNAs have been found to play important roles in apoptosis, cancer, development, differentiation, inflammation, longevity, and viral infection. There are a few reports describing miRNAs in the African malaria mosquito, Anopheles gambiae, on the basis of similarity to known miRNAs from other species. An. stephensi is the most important malaria vector in Asia and it is becoming a model Anopheline species for physiological and genetics studies.
We report the cloning and characterization of 27 distinct miRNAs from 17-day old An. stephensi female mosquitoes. Seventeen of the 27 miRNAs matched previously predicted An. gambiae miRNAs, offering the first experimental verification of miRNAs from mosquito species. Ten of the 27 are miRNAs previously unknown to mosquitoes, four of which did not match any known miRNAs in any organism. Twenty-five of the 27 Anopheles miRNAs had conserved sequences in the genome of a divergent relative, the yellow fever mosquito Aedes aegypti. Two clusters of miRNAs were found within introns of orthologous genes in An. gambiae, Ae. aegypti, and Drosophila melanogaster. Mature miRNAs were detected in An. stephensi for all of the nine selected miRNAs, including the four novel miRNAs (miR-x1- miR-x4), either by northern blot or by Ribonuclease Protection Assay. Expression profile analysis of eight of these miRNAs revealed distinct expression patterns from early embryo to adult stages in An. stephensi. In both An. stephensi and Ae. aegypti, the expression of miR-x2 was restricted to adult females and predominantly in the ovaries. A significant reduction of miR-x2 level was observed 72 hrs after a blood meal. Thus miR-x2 is likely involved in female reproduction and its function may be conserved among divergent mosquitoes. A mosquito homolog of miR-14, a regulator of longevity and apoptosis in D. melanogaster, represented 25% of all sequenced miRNA clones from 17-day old An. stephensi female mosquitoes. An. stephensi miR-14 displayed a relatively strong signal from late embryonic to adult stages. miR-14 expression is consistent during the adult lifespan regardless of age, sex, and blood feeding status. Thus miR-14 is likely important across all mosquito life stages.
This study provides experimental evidence for 23 conserved and four new microRNAs in An. stephensi mosquitoes. Comparisons between miRNA gene clusters in Anopheles and Aedes mosquitoes, and in D. melanogaster suggest the loss or significant change of two miRNA genes in Ae. aegypti. Expression profile analysis of eight miRNAs, including the four new miRNAs, revealed distinct patterns from early embryo to adult stages in An. stephensi. Further analysis showed that miR-x2 is likely involved in female reproduction and its function may be conserved among divergent mosquitoes. Consistent expression of miR-14 suggests that it is likely important across all mosquito life stages from embryos to aged adults. Understanding the functions of mosquito miRNAs will undoubtedly contribute to a better understanding of mosquito biology including longevity, reproduction, and mosquito-pathogen interactions, which are important to disease transmission.
microRNAs (miRNAs) are non-coding RNAs that are now recognized as a major class of gene-regulating molecules widely distributed in metozoans and plants [1, 2]. Many miRNA genes are transcribed by RNA polymerase II, yielding primary miRNAs (pri-miRNAs) of hundreds to thousands of bases in length [1, 3–5]. A given pri-miRNA can be either monocistronic, containing a sequence for one mature miRNA, or polycistronic, containing a sequence for multiple mature miRNA products [1, 3, 6]. In Drosophila, the pri-miRNA is processed by a Drosha-Pasha complex to yield pre-miRNA, small stem-loop structures that are approximately 70 nucleotides (nt) in length [7, 8]. These stem-loops, or hairpins, are then exported to the cytoplasm and processed by Dicer-1 to form an miRNA:miRNA* duplex [1, 3, 6]. The duplex molecules are separated by a helicase, and based upon the strength of 5' end pairing, one single strand is chosen as the mature miRNA . The opposing strand is referred to as the miRNA* strand, and believed to rapidly degrade following separation . Mature miRNAs associate with an Argonaute protein and bind their mRNA targets, which are often in the 3' untranslated region (UTR), resulting in inhibition of translation or possibly target mRNA degradation in animals . The "seed region" (bases 2–8 from the 5' end) contributes significantly to miRNA-target recognition [9, 10].
miRNAs have been found to play important roles in apoptosis, cancer, development, differentiation, inflammation, longevity, and viral infection [8, 11–17]. Estimates of the extent of miRNA gene regulation vary from 4% of transcripts in the Drosophila ovary  to a third of human genes . It is estimated that approximately 110 different miRNAs are expressed across the different life stages of D. melanogaster . In flies, the adult stage is characterized by significant miRNA expression . In a study to uncover Drosophila miRNAs, Lai et al. (2003) reported 38 putative miRNAs in the African malaria mosquito, Anopheles gambiae, that are conserved with Drosophila miRNAs . There are two additional reports describing An. gambiae miRNAs on the basis of similarity to known miRNAs [22, 23]. However, there is no direct experimental evidence for any of these miRNAs in mosquitoes. We are interested in identifying conserved as well as mosquito-specific microRNAs and exploring their potential functions in mosquito biology and mosquito-pathogen interactions. We carried out our cloning work on 17-day old adult female mosquitoes, which are highly relevant to disease transmission because it takes approximately two weeks for Plasmodium parasites to mature and become infective within a female mosquito . We used An. stephensi because this species is an important malaria vector in Asia and it is becoming a model Anopheline species for physiological and genetics studies. Here we report direct cloning and characterization of 23 conserved and four new miRNAs from the An. stephensi adult female. Comparative analysis uncovered the loss or significant change of two miRNA genes in a divergent mosquito Ae. aegypti. We also determined the expression profile of several selected miRNAs including the four new miRNAs across all life stages of An. stephensi. We performed further expression analysis on two miRNAs that are implicated in mosquito reproduction and longevity.
Classification of Cloned Small RNAs in An. stephensi.
% of Total Clones
miRNAs shown in Table 2 1
Possible miRNAs listed in Additional file 11
Low Quality or Short Sequences 3
Sequence and Characteristics of Cloned miRNAs in An. stephensi.
3R, 10270763 (-)
3R, 29002032 (+)
3R, 24898098 (+)
X, 21450383 (+)
2L, 18991766 (+)
2R, 28234532 (-)
ast-miR-277* 2 (Ia)
2R, 28234583 (-)
2L, 17362444 (-)
2R, 37890092 (-)
2R, 28252007 (-)
3L, 38943098 (+)
ast-miR-8* 2 (Ia)
3L, 38943058 (+)
2L, 15089338 (-)
2R, 13042084 (-)
2L, 18991809 (+)
2R, 28232720 (-)
X, 261196 (-)
2R, 37888125 (-)
3R, 51640581 (-)
2L, 37757111 (-)
2R, 37888805 (-)
3R, 5888649 (-)
X, 1228349 (+)
2R, 55572846 (-)
3L, 2905484 (+)
3R, 21181098 (-)
3R, 5819094 (-)
Confirmation of previously predicted Anopheline miRNAs by direct cloning and northern blot
Six new miRNAs in Anopheles that are similar to known miRNAs from other organisms
Two of the above miRNAs are worth noting. ast-miR-2a matches dme-miR-2a perfectly thus is named as miR-2a. There is an aga-miR-2 in miRBase (derived from two precursors aga-miR-2-1 and aga-miR-2-2), which is reverse complementary to ast-miR-2a. Sequence comparison showed that ast-miR-2a is not derived from aga-miR-2* because of the existence of multiple indels/mismatches between the alignment of aga-miR-2 and aga-miR-2*. The orientation of our cloned ast-mir-2a is consistent with the miRscan prediction based on An. gambiae and Ae. aegypti hairpins (Table 2, score 15.99). In addition, ast-mir-304 is reverse complementary to its top match (dme-miR-304) with 86% identity, which is described in the next section.
Anopheles miR-304 and miR-306 are in two separate miRNA clusters located in introns of protein-coding genes
Interestingly, miR-304 is closely flanked by miR-12 and miR-283 while miR-306 is in a different cluster with miR-9b and miR-79 (Figure 3). Clustering of An. gambiae miR-9b and miR-79 is noted on miRBase, but not miR-306; clustering of An. gambiae miR-304, -12, and -283 is not predicted in miRBase [33, 34]. Both clusters are within introns of protein coding genes. The miR-12, -304, -283 cluster occurs within a conserved gene of unknown function, while the miR-9b, -79, -306 cluster occurs within an ortholog of a gene coding for a Drosophila serine-threonine kinase group protein. The exons flanking each of the miRNA clusters are conserved between An. gambiae, Ae. aegypti, and D. melanogaster, which indicates that the clusters are orthologous. The order of miRNAs within the introns is conserved between An. gambiae and D. melanogaster, again supporting the miRNA status of the Anopheles miR-304 and miR-306. All miRNAs in the two clusters are in the same orientation as the flanking genes, indicating that these miRNAs may be transcribed from the promoters of their respective flanking genes, which is consistent with previous reports [35, 36]. In this regard, it is suggestive that, albeit reverse complementary to each other, the mature miR-304 in Anopheles and Drosophila are both correctly annotated because the orientation of transcription of miR-304 is consistent with the flanking gene in both species.
Discovery of four novel mosquito miRNAs that have no apparent homolog outside of mosquitoes
Expression profile of mosquito miRNAs across different developmental stages
miR-x2 in both An. stephensi and Ae. aegypti is predominantly found in the ovaries and its level is significantly reduced 72 hrs after blood feeding
We decided to carry out a detailed expression analysis of miR-x2 in both An. stephensi and Ae. aegypti. When different An. stephensi tissues were analyzed (Figure 6), miR-x2 showed strong signals in ovaries while no expression was detected in the midgut samples. The expression of miR-x2 was weak in the heads and the "remainders" (thorax plus abdomen without midgut and ovary). Following blood feeding, miR-x2 expression in the ovaries remained high at 24 hours post bloodmeal, but declined sharply by 72 hours post bloodmeal. miR-x2 was hardly detectable among the other tissues in either of the post bloodmeal time points. The same pattern of expression was observed in the distantly related Ae. aegypti (Figure 6). We also confirmed the lack of miR-x2 expression in adult males in both An. stephensi and Ae. aegypti (Figure 6). Thus miR-x2 expression profile is conserved between the two divergent mosquitoes in all samples tested.
Two criteria are critical for demonstrating a valid miRNA [31, 38]. First, expression of an approximately 22 nt RNA should be detected by small RNA cloning or by RNA hybridization methods such as northern blot [31, 38]. Second, the miRNA should be traceable back to a precursor with a hairpin structure [31, 38]. All 27 miRNAs described in this report, including two miRNA* sequences, met the above criteria. The expression was indicated by direct cloning of small RNAs from 17-day old female An. stephensi samples and the precursor hairpin was identified using the genome assembly of a related Anopheles species of the same subgenus, An. gambiae. Additional evidence to validate miRNA status that was proposed in the literature includes proof of processing to a mature form via Dicer, and conservation of the mature sequence and its precursor [31, 38]. For all 27 Anopheles miRNAs, we were able to identify conserved sequences encompassing the entire hairpin structure either from the genome assembly of a divergent mosquito species, Ae. aegypti (25 out of 27) or from D. melanogaster (see Additional file 1). Furthermore, miRscan analysis based on conserved hairpin alignment provided strong support for 26 of the 27 Anopheline miRNAs. In the case of ast-miR-304, which has no conserved sequence in Ae. aegypti and no miRscan score, a closer examination of the miR-304 hairpin from An. gambiae suggests that it meets all of previously described criteria for miRNA structures [31, 32]. Furthermore, the Anopheles miR-304 and the D. melanogaster miR-304 are both in a conserved miRNA cluster and both are flanked by miR-12 and miR-283, thus lending additional support for the validity of Anopheles miR-304. Finally miRNA expression was detected for all nine selected ast-miRNAs during either northern blot or RPA, complementing our cloning results. Thus we feel that the overall support for the presence of these 27 miRNAs in Anopheles mosquitoes is strong. One miRNA, ast-miR-2a, is worth noting. Cloning results for ast-miR-2a, comparisons to dme-miR-2a, and miRscan predictions are all consistent. However, ast-miR-2a is reverse complementary to the aga-miR-2 reported in miRBase.
As mentioned earlier, miR-304 and miR-306 do not have obvious homologs in Ae. aegypti. In both An. gambiae and D. melanogaster, miR-304 is closely flanked by miR-12 and miR-283 while miR-306 is in a different cluster with miR-9b and miR-79 (Figure 3). Both clusters are in the intron of orthologous genes in the two species. Comparisons of the two miRNAs clusters in these three species (Figure 3) would suggest that miR-304 and miR-306 may have either been lost in Ae. aegypti or evolved to significantly different sequences. As the mature miRNA sequences of the Anopheles and Drosophila miR-304 are reverse complementary to each other, it is possible that the miR-304 precursor sequence may have undergone an inversion after the two dipteran lineages separated. Thus comparisons of clear orthologous miRNA clusters will likely shed new light on the evolution (expansion, loss, rearrangement) of miRNAs, which has not been well studied.
The expression profiles of the eight mosquito miRNAs are informative. When D. melanogaster data are available for comparison as in the cases of let-7 and miR-9a, similar expression profiles were found between An. stephensi and D. melanogaster. This is not surprising as these conserved miRNAs are likely to have similar functions in these Dipteran insects. miR-14 was expressed in the same stages in An. stephensi as in D. melanogaster. We provided extended expression analysis on miR-14, the miRNA that represents 25% of the sequenced miRNAs during the cloning of the 17-day old An. stephensi samples. We showed that the miR-14 level increased slightly during embryonic development and remained relatively high through larvae, pupae and adult stages and we did not observe significant changes in adults regardless of age, sex, and blood feeding status. These results do not necessarily imply that mosquito miR-14 is important to longevity, a function of miR-14 demonstrated in D. melanogaster. Nonetheless, it appears that miR-14 is important across different mosquito life stages from embryos to aged adults. Further research on the targets and function of miR-14 in mosquitoes will help determine whether it is important to mosquito longevity. The expression profile of miR-x1 through miR-x4 confirmed our cloning results for these new miRNAs, and provided useful information for future research into their functions. The expression of miR-x2 in An. stephensi was adult specific as well as female specific. miR-x2 was predominantly expressed in the ovary and its level was reduced 72 hrs after blood feeding. These results indicate that miR-x2 is likely involved in An. stephensi female reproduction. The same pattern of expression was shown for miR-x2 in Ae. aegypti, a mosquito that is highly divergent from An. stephensi. Thus the function of miR-x2 may be conserved among divergent mosquitoes.
We have identified 10 miRNAs previously unknown to mosquitoes, four of which did not match any known miRNAs in any organism. These four miRNAs are conserved between An. stephensi, An. gambiae, and Ae. aegypti, suggesting important functions possibly common to all mosquitoes as the Anopheles and Aedes genera are two of the most divergent among all mosquitoes, separated approximately 145–200 million years ago . Considering the modest number of small RNA clones we sequenced during this study and the specific developmental stage of our total RNA source, we suspect that a number of novel miRNAs still await to be discovered in mosquitoes. This study demonstrates the importance of direct cloning and expression profile analysis to the identification and characterization of conserved as well as mosquito-specific microRNAs, some of which will likely regulate genes that significantly affect mosquito biology and perhaps mosquito-pathogen interactions.
Three research articles became available when this manuscript was being reviewed. Two of these papers describe additional miRNAs from Drosophila [40, 41] and a third paper reports cloning of 10 distinct miRNAs and primer extension analysis of an additional 8 miRNAs from An. gambiae midgut samples infected with Plasmodium bergei . Comparisons between the 27 miRNAs reported in this manuscript and those published in the above three papers showed that An. stephensi miR-x1 and miR-x2 are nearly identical to An. gambiae miR-996 and miR-989, respectively  and homologs of these two miRNAs are found in the expanded list of Drosophila miRNAs [40, 41]. Thus miR-x1 and miR-x2 are not unique to mosquitoes. Ast-miR-76 showed a perfect match to the newly identified miR-981 in Drosophila [40, 41], thus ast-miR-76 should be renamed as ast-miR-981. An additional five of the cloned miRNAs reported in Winter et al.  match five of the 27 miRNAs reported in this manuscript. They are let-7, miR-281, miR-34, miR-12, and miR-306. An. stephensi miR-x3 and miR-x4 showed no similarity to any miRNAs in the miRbase which included updates from the above three papers. Thus miR-x3 and miR-x4 remain novel and potentially specific to mosquitoes. The final names of these miRNAs as assigned by miRBase are shown in Table 2. We also compared miRNA expression profiles obtained in our study using northern analysis with the profiles of overlapping miRNAs reported in Winter et al. , which were obtained using primer extension. Winter and colleagues state that miR-989 (miR-x2) is expressed only in midguts of An. gambiae (page 6958). However, we have shown that miR-x2 (miR-989) is predominantly expressed in ovaries in both An. stephensi and Ae. aegypti and miR-x2 was not detected in midguts under our condition. We note that the Figure 2 of Winter et al.  does show a much stronger expression in "leftovers", which include ovaries, than in midguts.
This study provides experimental evidence for 23 conserved and four new microRNAs in An. stephensi mosquitoes. Comparisons between miRNA gene clusters in An. gambiae, Ae. aegypti, and D. melanogaster suggest the loss or significant change of two miRNA genes in Ae. aegypti. Expression profile analysis of eight miRNAs, including the four new miRNAs, revealed distinct patterns from early embryo to adult stages in An. stephensi. Further analysis showed that miR-x2 is likely involved in female reproduction and its function may be conserved among divergent mosquitoes. Consistent expression of miR-14 suggests that miR-14 is likely important across all mosquito life stages from embryos to aged adults. Understanding the functions of mosquito miRNAs will undoubtedly contribute to a better understanding of mosquito biology including longevity, reproduction, and mosquito-pathogen interactions, which are important to disease transmission.
An. stephensi (Indian wild-type strain) and Ae. aegypti (Liverpool strain) mosquitoes were maintained in humidified incubators at 27°C on a 12 hour light:dark cycle.
Small RNA cloning and sequencing
For small RNA cloning, approximately 1000 female An. stephensi adults were collected at 17 days post-emergence. Mosquitoes were fed blood and allowed to oviposit prior to collection. Aged females are the most relevant mosquito life stage for malaria transmission as it takes approximately two weeks for Plasmodium parasites to mature within the mosquito to an infective stage . RNA was isolated using a mirVana small RNA Isolation Kit (Ambion, Austin, TX). Subsequent steps were performed following published protocols  from the David Bartel laboratory (Whitehead Institute, MIT, Cambridge, MA). RNA oligo markers (18-mer, 24-mer; sequence as detailed from ) were purchased from Dharmacon (Lafayette, CO) and 2' ends were deprotected prior to use as indicated by Dharmacon. These oligos were end-labeled with 32P (32P-γ ATP, Perkin-Elmer, Waltham, MA) by T4 polynucleotide kinase (Promega, Madison, WI) and added to the sample to isolate RNAs between 18 and 24 nt in length. Next, 5' and 3' linkers were sequentially ligated to the isolated small RNA as well as the RNA markers. cDNA were generated by RT-PCR using primers derived from the two linker sequences. cDNA were ligated into a 2.1 TOPO TA vector (Invitrogen, Carlsbad, CA). Ligated plasmids were transformed in One-Shot Mach 1-T1R Competent cells (Invitrogen). We did not concatemerize cDNA prior to cloning. Sequencing was performed by Virginia Bioinformatics Institute core facilities and the University of Washington High Throughput Genomics Unit.
Small RNA Sequence Analyses
Sequences obtained from cloning were analyzed with ClustalW  to identify inserts and orientation. After removing low quality sequences determined by inspection of chromatographs, insert sequences were analyzed using BLAST  against different databases as described below to identify matches to known miRNAs, mRNAs, tRNAs, rRNA, and other small RNAs (Table 1). Sequences that were identical to previously predicted An. gambiae miRNAs were identified by comparing with the miRBase miRNA registry  and An. gambiae miRNA predictions listed in  using BLAST . Comparisons to miRBase also revealed sequences that match miRNAs from other organisms but were not reported in mosquitoes. The remaining An. stephensi small RNA sequences that matched the An. gambiae genome assembly were further analyzed to uncover new miRNAs from mosquitoes (Table 2). An. gambiae genome assembly was used to retrieve miRNA precursors for secondary structural analysis because there was little sequence information available for the An. stephensi genome and the two species are in the same subgenus Celia. The matching An. gambiae sequences plus 100 nt flanking sequences were obtained through Ensembl  and lowest energy conformations were generated by Vienna RNAfold . The above mentioned An. gambiae precursor sequences were used to search for conserved sequences in Ae. aegypti, and in some cases, in D. melanogaster. Conserved sequence pairs were then examined via miRscan . One An. stephensi miRNA (ast-miR-304) was unscorable by miRscan (Table 2) and the precursor of this miRNA was examined using the biogenesis criteria described in the literature: 1) the miRNA appears in the stem of a hairpin structure; 2) there is considerable similarity between the miRNA and miRNA*; and 3) few bulges if any are present within the miRNA:miRNA* pairing in the stem [31, 32].
Northern hybridizations shown in Figure 1 were conducted using miRCURY-LNA probes from Exiqon (Vedbaek, Denmark) for eight different miRNAs. Total RNA was isolated from An. stephensi of the same age or developmental stage and used for each blot. Ten micrograms of total RNA was loaded for each sample in a 15% denaturing polyacrylamide gel. After initial size confirmation using let-7 RNA, 19 nt and 23 nt single-stranded DNA oligos were used as size markers in subsequent experiments. Gels were stained with ethidium bromide for verification of equal loading and RNA quality before transferring to a BrightStar-Plus membrane (Ambion). Subsequent steps were based upon Wienholds et al. . Following UV crosslinking, membranes were prehybridized in a rotating oven for 30 minutes at 42°C using ULTRAhyb-Oligo Hybridization Buffer (Ambion), followed by overnight hybridization in the same buffer at 42°C with a final concentration of 0.1 nM digoxigenin (DIG)-labeled antisense miRCURY LNA probe (Exiqon). All probes used in northern blots were designed to hybridize with An. stephensi miRNAs. The membranes were washed twice for 35 minutes in an Ambion-recommended wash buffer (2 × SSC, 0.5% SDS) at 42°C, and then once in low stringency buffer for 5 minutes at room temperature. The membranes were incubated for 30 minutes in blocking buffer followed by a 1 hr incubation with Anti-DIG-alkaline phosphatase fAb (Roche, Basel, Switzerland) in blocking buffer. The membranes were washed three times for 15 minutes in low stringency wash buffer followed by twice with alkaline phosphatase buffer. The membranes were then immersed in CDP-Star solution (Roche) for 5 minutes, and placed inside saran wrap for exposure to X-ray film for 30 minutes.
For the northern blots shown in Figure 5, a similar procedure was followed. The samples were 3, 5, 10, 17, and 24 day old An. stephensi adult females as well as 3, 5, 10, and 17 day old An. stephensi adult males that were maintained on sugar water. The same cohort of adult females, which were fed on blood on day 5 after emergence and allowed to oviposit two days later, were collected at day 10, 17, and 24. RNA isolation and northerns were carried out as detailed above. For tissue sample northern blots (Figure 6), five day old adult An. stephensi females were split into two groups, one maintained on sugar water, the other provided a blood meal then maintained on sugar water. Approximately twenty-four hrs post bloodfeeding, heads, ovaries, midguts and remainders (everything left behind) were collected from 60 bloodfed and 60 non-bloodfed six day old mosquitoes. This procedure was repeated again at 72 hrs after blood feeding with eight day old mosquitoes. The mosquitoes were not permitted to oviposit. All tissues were stored in RNA later (Ambion) during collection, then vortexed and stored in -80°C until RNA isolation as above using a mir Vana miRNA Isolation kit (Ambion). Northern hybridization was performed as describe above except that 5 ug of total RNA were used. The same procedure was followed for Ae. aegypti mosquito tissue analysis. In separate experiments, 10 ug of total RNA isolated from whole body were used to compare the expression of miR-x2 in males and non-bloodfed females of both species. We have performed replicates for all expression profile analyses and obtained consistent expression patterns between replicates.
Ribonuclease Protection Assay (RPA)
RPA was used to examine ast-miR-76. Double-stranded oligos with a T7 promoter 5' of an antisense sequence of the miRNA were produced by annealing two single-stranded oligos. The annealed oligos were used to synthesize an RNA probe that was 9 nt longer than the miRNA, as seen below:
Our design also considered the possibility for future multiplexing by adding extra 4 adenosines immediately following +1 which were indigestible according to the instruction manual of the "mirVana miRNA Probe Construction Kit" (Ambion). This permitted the undigested probe to run at ~29 nt, and digested/protected probe to run at ~24 nt. We utilized a "MEGAscript RNAi kit" (Ambion) to synthesize radiolabeled (32P-UTP, Perkin-Elmer) antisense transcripts. Templates were synthesized for 4 hours, and purified from acrylamide gels by elution and isopropanol precipitation (1× volume) with 1/10 volume 3 M NaOAc and glycogen. 50,000 cpm of probe was added to 5 ug total RNA and denatured at 96°C for 3 minutes followed by 42°C hybridization overnight. Next, the hybridized samples were digested with RNase A/T1 for 45 minutes at 37°C to remove single-stranded RNA, and to trim unhybridized regions of the probes. Afterwards, the samples were ethanol precipitated and resuspended in 5 ul 2 × loading buffer. Resuspended samples were denatured at 96°C for 5 minutes followed by a quick chill on ice before loading. RNA were separated using a 15% denaturing polyacrylamide gel at 150 V. Locations of bands were examined by X-ray film exposure, using an intensifying screen. Size markers were as described for northern blot analysis.
Locked nucleic acid
We would like to thank Randy Saunders for maintaining mosquito colonies, and Ray Miller for improving the resolution of the hairpin images. We thank Jim Biedler and Zach Adelman for reviewing a draft of the manuscript. This work is supported by NIH grants (AI070854 and AI063252) to Z.T.
- Miska EA: How microRNAs control cell division, differentiation and death. Curr Opin Genet Dev. 2005, 15 (5): 563-568. 10.1016/j.gde.2005.08.005.PubMedView ArticleGoogle Scholar
- Du T, Zamore PD: Beginning to understand microRNA function. Cell Research. 2007, 17 (8): 661-663. 10.1038/cr.2007.67.PubMedView ArticleGoogle Scholar
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.PubMedView ArticleGoogle Scholar
- Pasquinelli AE, Hunter S, Bracht J: MicroRNAs: a developing story. Curr Opin Genet Dev. 2005, 15 (2): 200-205. 10.1016/j.gde.2005.01.002.PubMedView ArticleGoogle Scholar
- Zhou X, Ruan J, Wang G, Zhang W: Characterization and identification of MicroRNA core promoters in four model species. PLoS Comput Biol. 2007, 3 (3): e37-10.1371/journal.pcbi.0030037.PubMedPubMed CentralView ArticleGoogle Scholar
- Lau NC, Lim LP, Weinstein EG, Bartel DP: An Abundant Class of Tiny RNAs with Probable Regulatory Roles in Caenorhabditis elegans. Science. 2001, 294: 858-862. 10.1126/science.1065062.PubMedView ArticleGoogle Scholar
- Denli A, Tops BBJ, Plasterk RHA, Ketting RF, Hannon GJ: Processing of primary microRNAs by the Microprocessor complex. Nature. 2004, 432: 231-235. 10.1038/nature03049.PubMedView ArticleGoogle Scholar
- Bushati N, Cohen SM: microRNA Functions. Annual Review of Cell and Developmental Biology. 2007, 23 (1): 175-205. 10.1146/annurev.cellbio.23.090506.123406.PubMedView ArticleGoogle Scholar
- Brennecke J, Stark A, Russell RB, Cohen SM: Principles of microRNA-target recognition. PLoS Biol. 2005, 3 (3): e85-10.1371/journal.pbio.0030085.PubMedPubMed CentralView ArticleGoogle Scholar
- Grimson A, Farh KK-H, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP: MicroRNA Targeting Specificity in Mammals: Determinants beyond Seed Pairing. Mol Cell. 2007, 27: 91-105. 10.1016/j.molcel.2007.06.017.PubMedPubMed CentralView ArticleGoogle Scholar
- Lai EC: miRNAs: Whys and Wherefores of miRNA-Mediated Regulation. Curr Biol. 2005, 15 (12): R458-R460. 10.1016/j.cub.2005.06.015.PubMedView ArticleGoogle Scholar
- Du T, Zamore PD: microPrimer: the biogenesis and function of microRNA. Development. 2005, 132 (21): 4645-4652. 10.1242/dev.02070.PubMedView ArticleGoogle Scholar
- Kim VN: Small RNAs: Classification, Biogenesis, and Function. Mol Cells. 2005, 19 (1): 1-15. 10.1016/j.molcel.2005.05.026.PubMedView ArticleGoogle Scholar
- Ambros V, Chen X: The regulation of genes and genomes by small RNAs. Development. 2007, 134 (9): 1635-1641. 10.1242/dev.002006.PubMedView ArticleGoogle Scholar
- O'Connell RM, Taganov KD, Boldin MP, Cheng G, Baltimore D: MicroRNA-155 is induced during the macrophage inflammatory response. PNAS. 2007, 104 (5): 1604-1609. 10.1073/pnas.0610731104.PubMedPubMed CentralView ArticleGoogle Scholar
- le Sage C, Nagel R, Egan DA, Schrier M, Mesman E, Mangiola A, Anile C, Maira G, Mercatelli N, Ciafrè SA: Regulation of the p27Kip1 tumor suppressor by miR-221 and miR-222 promotes cancer cell proliferation. EMBO J. 2007, 26: 3699-3708. 10.1038/sj.emboj.7601790.PubMedPubMed CentralView ArticleGoogle Scholar
- Cullen BR: Viruses and microRNAs. Nature Genetics. 2006, 38: S25-S30. 10.1038/ng1793.PubMedView ArticleGoogle Scholar
- Nakahara K, Kim K, Sciulli C, Dowd SR, Minden JS, Carthew RW: Targets of microRNA regulation in the Drosophila oocyte proteome. PNAS. 2005, 102 (34): 12023-12028. 10.1073/pnas.0500053102.PubMedPubMed CentralView ArticleGoogle Scholar
- Lewis BP, Burge CB, Bartel DP: Conserved Seed Pairing, Often Flanked by Adenosines, Indicates that Thousands of Human Genes are MicroRNA Targets. Cell. 2005, 120 (1): 15-20. 10.1016/j.cell.2004.12.035.PubMedView ArticleGoogle Scholar
- Lai E, Tomancak P, Williams R, Rubin G: Computational identification of Drosophila microRNA genes. Genome Biology. 2003, 4 (7): R42-10.1186/gb-2003-4-7-r42.PubMedPubMed CentralView ArticleGoogle Scholar
- Aravin AA, Lagos-Quintana M, Yalcin A, Zavolan M, Marks D, Snyder B, Gaasterland T, Meyer J, Tuschl T: The Small RNA Profile during Drosophila melanogaster Development. Dev Cell. 2003, 5 (2): 337-350. 10.1016/S1534-5807(03)00228-4.PubMedView ArticleGoogle Scholar
- Wang X, Zhang J, Li F, Gu J, He T, Zhang X, Li Y: MicroRNA identification based on sequence and structure alignment. Bioinformatics. 2005, 21 (18): 3610-3614. 10.1093/bioinformatics/bti562.PubMedView ArticleGoogle Scholar
- Chatterjee R, Chaudhuri K: An approach for the identification of microRNA with an application to Anopheles gambiae. Acta Biochim Pol. 2006, 53: 303-309.PubMedGoogle Scholar
- U.S. Centers for Disease Control. [http://www.cdc.gov/malaria/biology]
- The miRBase Registry. [http://microrna.sanger.ac.uk]
- Xu P, Vernooy SY, Guo M, Hay BA: The Drosophila microRNA Mir-14 suppresses cell death and is required for normal fat metabolism. Curr Biol. 2003, 13: 790-795. 10.1016/S0960-9822(03)00250-1.PubMedView ArticleGoogle Scholar
- Valoczi A, Hornyik C, Varga N, Burgyan J, Kauppinen S, Havelda Z: Sensitive and specific detection of microRNAs by northern blot analysis using LNA-modified oligonucleotide probes. Nucl Acids Res. 2004, 32 (22): e175-10.1093/nar/gnh171.PubMedPubMed CentralView ArticleGoogle Scholar
- miRscan Web Server. [http://genes.mit.edu/mirscan]
- Lim LP, Lau NC, Weinstein EG, Abdelhakim A, Yekta S, Rhoades MW, Burge CB, Bartel DP: The microRNAs of Caenorhabditis elegans. Genes Dev. 2003, 17: 991-1008. 10.1101/gad.1074403.PubMedPubMed CentralView ArticleGoogle Scholar
- Lim LP, Glasner ME, Yekta S, Burge CB, Bartel DP: Vertebrate MicroRNA Genes. Science. 2003, 299: 1540-10.1126/science.1080372.PubMedView ArticleGoogle Scholar
- Ambros V, Bartel B, Bartel DP, Burge CB, Carrington JC, Chen X, Dreyfuss G, Eddy SR, Griffiths-Jones SAM, Marshall M, Matzke M, Ruvkun G, Tuschl T: A uniform system for microRNA annotation. RNA. 2003, 9: 277-279. 10.1261/rna.2183803.PubMedPubMed CentralView ArticleGoogle Scholar
- Gu Z, Eleswarapu S, Jiang H: Identification and characterization of microRNAs from the bovine adipose tissue and mammary gland. FEBS Lett. 2007, 581: 981-988. 10.1016/j.febslet.2007.01.081.PubMedView ArticleGoogle Scholar
- Griffiths-Jones S: The microRNA Registry. Nucl Acids Res. 2004, 32 (suppl_1): D109-111. 10.1093/nar/gkh023.PubMedPubMed CentralView ArticleGoogle Scholar
- Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucl Acids Res. 2006, 34 (suppl_1): D140-144. 10.1093/nar/gkj112.PubMedPubMed CentralView ArticleGoogle Scholar
- Smalheiser NA: EST analyses predict the existence of a population of chimeric microRNA precursor-mRNA transcripts expressed in normal human and mouse tissues. Genome Biol. 2003, 4 (7): 403-10.1186/gb-2003-4-7-403.PubMedPubMed CentralView ArticleGoogle Scholar
- Kim YK, Kim VN: Processing of intronic microRNAs. EMBO J. 2007, 26 (3): 775-783. 10.1038/sj.emboj.7601512.PubMedPubMed CentralView ArticleGoogle Scholar
- Sempere LF, Sokol NS, Dubrovsky EB, Berger EM, Ambros V: Temporal regulation of microRNA expression in Drosophila melanogaster mediated by hormonal signals and broad-Complex gene activity. Dev Biol. 2003, 259 (1): 9-18. 10.1016/S0012-1606(03)00208-2.PubMedView ArticleGoogle Scholar
- Berezikov E, Cuppen E, Plasterk RHA: Approaches to microRNA discovery. Nat Genet. 2006, 38: 52-57. 10.1038/ng1914.View ArticleGoogle Scholar
- Krzywinski J, Grushko OG, Besansky NJ: Analysis of the complete mitochondrial DNA from Anopheles funestus: An improved dipteran mitochondrial genome annotation and a temporal dimension of mosquito evolution. Mol Phylogenet Evol. 2006, 39 (2): 417-423. 10.1016/j.ympev.2006.01.006.PubMedView ArticleGoogle Scholar
- Ruby JG, Stark A, Johnston WK, Kellis M, Bartel DP, Lai EC: Evolution, biogenesis, expression, and target predictions of a substantially expanded set of Drosophila microRNAs. Genome Res. 2007, 17: 1850-1864. 10.1101/gr.6597907.PubMedPubMed CentralView ArticleGoogle Scholar
- Stark A, Kheradpour P, Parts L, Brennecke J, Hodges E, Hannon GJ, Kellis M: Systematic discovery and characterization of fly microRNAs using 12 Drosophila genomes. Genome Res. 2007, 17 (12): 1865-1879. 10.1101/gr.6593807.PubMedPubMed CentralView ArticleGoogle Scholar
- Winter F, Edaye S, Huttenhofer A, Brunel C: Anopheles gambiae miRNAs as actors of defence reaction against Plasmodium invasion. Nucl Acids Res. 2007, 35 (20): 6953-6962. 10.1093/nar/gkm686.PubMedPubMed CentralView ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position- specific gap penalties and weight matrix choice. Nucl Acids Res. 1994, 22 (22): 4673-4680. 10.1093/nar/22.22.4673.PubMedPubMed CentralView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucl Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- Ensembl. [http://www.ensembl.org]
- Vienna RNA Secondary Structure Prediction. [http://rna.tbi.univie.ac.at/cgi-bin/RNAfold.cgi]
- Wienholds E, Kloosterman WP, Miska E, Alvarez-Saavedra E, Berezikov E, de Bruijn E, Horvitz HR, Kauppinen S, Plasterk RHA: MicroRNA Expression in Zebrafish Embryonic Development. Science. 2005, 309: 310-311. 10.1126/science.1114519.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.