GSE23292
High-throughput sequencing of small RNAs in Heliconius melpomene
Background: Heliconius butterflies are an excellent model system for studies of adaptive convergent and divergent phenotypic traits. Wing colour patterns are used as signals to both predators and potential mates and are inherited in a Mendelian manner. The underlying genetic mechanisms of pattern formation have been studied for many years and shed light on broad issues, such as the repeatability of evolution. In Heliconius melpomene, the yellow hindwing bar is controlled by the HmYb locus and several genes in this region show expression pattern differences across races. MicroRNAs (miRNAs) are important post-transcriptional regulators of gene expression that have key roles in many biological processes, including development. It seems likely that miRNAs could act as downstream regulators of genes involved in wing development, patterning and pigmentation. For this reason we characterised miRNAs in developing butterfly wings and examined differences in their expression between colour pattern races. Results: We sequenced small RNA libraries from two colour pattern races and detected 142 Heliconius miRNAs with homology to others found in miRBase. Several highly abundant miRNAs appeared to be differentially expressed between colour pattern races and this was tested further in different developing pupal wing stages using Northern blots. These revealed that differences in expression were due to developmental stage rather than colour pattern. Assembly of sequenced reads to the HmYb region identified miR-193 and miR-2788; located 2380bp apart in an intergenic region. A search for miRNAs in all available H. melpomene BAC sequences (~2.5Mb) did not reveal any other miRNA genes and no novel miRNAs were predicted. There were several regions where other small RNA sequences assembled to the HmYb region and appeared to be differentially expressed.These might represent other regulatory RNAs. Conclusions: Here we describe the first butterfly miRNAs and characterise their expression in developing wings. Some show differences in expression across developing pupal stages. Two miRNAs were located in the HmYb region. Future work will examine the expression of these miRNAs in different colour pattern races and identify miRNA targets among wing patterning genes. Overall design: High-throughput sequencing of Heliconius melpomene endogenous small RNAs. Size fractionated small RNA from total RNA extracts of two different Heliconius melpomene races (Heliconius melpomene melpomene and Heliconius melpomene rosina) were isolated from wing tissue using miRVana kit. 100µg RNA from 11  individuals of different developmental stages was pooled for each race as follows: 4.1% larval stage <1; 2% larval stage 1-1.75; 2.9% larval stage  2-2.5; 22% larval stage 2.75-3; 19% larval stage > 3; 25% early pupae; 25%  mid-melanin pupae. Sequences were ligated to adapters, purified again and reverse transcribed. After PCR amplification the sample was subjected to Solexa/Illumina high throughput pyrosequencing. Please see www.illumina.com for details of the sequencing technology.