A micro-RNA is the effector gene of a classic evolutionary hotspot locus

Abstract

In Lepidoptera (butterflies and moths), the genomic region around the gene cortex is a 'hotspot' locus, repeatedly used to generate intraspecific melanic wing color polymorphisms across 100-million-years of evolution. However, the identity of the effector gene regulating melanic wing color within this locus remains unknown. Here, we show that none of the four candidate protein-coding genes within this locus, including cortex, serve as major effectors. Instead, a micro-RNA (miRNA), mir-193, serves as the major effector across three deeply diverged lineages of butterflies, and its function is conserved in Drosophila. In Lepidoptera, mir-193 is derived from a gigantic long non-coding RNA, ivory, and it functions by directly repressing multiple pigmentation genes. We show that a miRNA can drive repeated instances of adaptive evolution in animals.

Fig. 1.. mir-193 is the major melanic color regulator in the cortex locus.

(A) Biogenesis and silencing machinery of miRNAs. (B) Disrupting Drosha or Dicer processing sites using CRISPR-Cas9 inhibits the biogenesis of mature miRNAs. (C) The core cortex locus in B. anynana showing two highly conserved miRNAs and four candidate protein-coding genes. (D) Homozygous mutant lines of mir-193 and mir-2788 in B. anynana and their corresponding genotypes. (E) Expression levels of the guide strand of mir-193, miR-193–3p, and the two mature strands of mir-2788, across the corresponding mutant lines and WT. n=3–4; ns: not significant; *: p<0.05; **: p<0.01; ***: p<0.001; Expression levels with the same letter are not significantly different from each other; Error bar: SEM. (F) Phylogenetic placement of the three model butterfly species in the tree of lepidopterans previously mapped to the cortex locus. Phylogeny is from (28). Mosaic knock-outs (mKOs) of mir-193 in (G) Pieris canidia and (H) Papilio polytes. Images were horizontally flipped when necessary.

Fig. 2.. The lncRNA ivory functions as primary mir-193.

(A) A gigantic lncRNA, ivory, was found in the RNA-seq data of three mir-193 mutants, with a deeply conserved transcription start site (TSS) whose chromatin accessibility increased during the larval-pupal transition (shaded), and an undefined 3’ terminus. No open chromatin was found around the miRNA region (shaded). (B) Time series expression of miRNA mature strands, miR-193–3p, miR-2788–5p, and ivory. (C) Expression levels of ivory across the miRNA mutants and WT. n=3–4; ns: not significant; *: p<0.05; **: p<0.01; ***: p<0.001; Error bar: SEM. (D-D””) Spatial expression (HCR) of pri-mir-193, pri-mir-2788, and ivory in the B. anynana ‘eyespot’ wing color pattern and (E-E”’) their expression signals within individual nuclei. (F) mir-193 and ivory TSS mosaic knock-out (mKO) phenotypes in B. anynana, P. canidia, and P. polytes. Images were horizontally flipped when necessary.

Fig. 3.. mir-193 directly targets multiple pigmentation genes.

(A) Differentially expressed genes (padj<0.01) across sib-paired female mir-193 m4 mutant and WT wing tissues across wing development. Genes highlighted were previously associated with butterfly pigmentation, or potential color regulators. Candidates for in-vitro validation are in bold (B) Four putative binding sites of miR-193–3p, the guide strand of mir-193, were found in CDS and/or 3’UTR regions of three candidate genes, ebony (e), Esp1, and yellow-e3 (yel-e3) that are upregulated in Day6 mutant wings. Full sequence complementarity (allowing G:U wobble base-pairing) between the seed region (nucleotides 2–8 from the 5’ terminus of the miRNA guide strand) and the miRNA targets are highlighted. (C) Dual luciferase reporter assay was used to validate the miRNA-target silencing across the four predicted binding sites in-vitro, with a concentration gradient of miR-193–3p mimic. n=6; *: p<0.05; **: p<0.01; ***: p<0.001; Error bar: SEM.

Fig. 4.. mir-193 is an ancestral melanic color regulator.

(A) Genomic context of D. melanogaster mir-193. (B) Phenotypes of pannier Gal4 driver (pnr-Gal4) transgenic lines expressing either a mir-193 sponge with a sequence of seed-complementary binding sites for the guide strand miR-193–3p (reduced mir-193), or extra mir-193 precursors (increased mir-193), or a mir-193 sponge with scrambled miR-193–3p binding sites (control). (C) Functional conservation of mir-193 across a broader animal phylogeny where mir-193 is deeply conserved. Phylogeny is from (28).