A single master regulatory gene optix underlies both color and iridescence in butterflies
Ontology highlight
ABSTRACT: It is debated to what extent individual regulatory genes might underlie non-gradualistic evolution of complex morphological traits. We have shown that a single regulatory gene, optix, is required for development of all color and iridescence in nymphalid butterfly wings by CRISPR/Cas9 genome editing. To investigate how wing features are controlled by optix, we used mRNA-seq to compare transcript abundance in wild type and optix knockout wings of V. cardui and J. coenia. We generated an average of 269 million reads 36 bp paired-end Nexseq 500 reads from 16 libraries and identified numerous candidate genes including known genes involved in the insect melanin and ommochrome biosynthetic pathways and suspected components potential involved in scale structure.
Project description:In the eastern United States the buckeye butterfly, Junonia coenia, shows a seasonal wing color polyphenism where adults emerging in the spring are pale brown, while those emerging in the autumn are dark red. This variation can be artificially induced in laboratory colonies, thus making J. coenia a useful model system to examine the developmental basis of phenotypic plasticity. We used RNA-seq to generate the first set of assembled transcripts for this species while simultaneously quantifying relative gene expression associated with development of alternative seasonal color morphs. The assembled consolidated wing transcriptome was 77.55 Mb. 16,251 contigs of over 1000bp in length were assembled, of which 3,145 were differentially expressed between stages and/or color morphs. Depending on the developmental stage, between 547 and 1420 transcripts were significantly differentially expressed between brown and red wing morphs. These extensive differences in gene expression stand in stark contrast to the much smaller numbers found in previous studies on genetic wing pattern variation, and suggest that environmentally induced phenotypic shifts may arise from very broad systemic processes. Overall gene ontology (GO) analyses revealed that genes associated with structural constituents of ribosomes and oxygen transport were significantly upregulated in the pale brown morph, while genes associated with peptidase activity were very significantly upregulated in the dark red morph. Focused analyses of candidate endocrine and pigmentation pathways revealed a number of notable genes upregulated in the red morph, including several ecdysone-related genes and cinnabar, an ommochrome pigment gene implicated in color pattern variation in other butterflies. Surprisingly, we found numerous melanin-related transcripts, including tan and yellow-family genes, strongly upregulated in the red morph, leading us to speculate that red pigmentation in autumn J. coenia may include red or brown melanins in addition to ommochromes. While we identified several endocrine and pigmentation genes as obvious candidates for color morph differentiation, we speculate that the majority of gene expression differences we observed were due to thermal stress response. The buckeye transcriptome provides a basis for further developmental studies of phenotypic plasticity. mRNA profiling of hind wings from 4 developmental stages of two color morphs (Rosa and Linea) of the buckeye butterfly (J. coenia), generated by deep sequencing, in triplicate, using Illumina GAII or HiSeq 2000.
Project description:Background: We studied the chromatin accessibility landscsape in wings during butterfly metamorphosis, and investigate which transcription factors might be driving changes in accessibility Methods: We sequencing the Junonia coenia genome, and we studied chromatin accessibility using ATAC seq in multiple stages of wing development in both forewings and hindwings. For sites showing a large change in accessibility, we investigate which motifs are enriched, and correlate this with changes in gene expression of associated transcription factors. We confirm promising candidates with ChIP-seq Results: We find a highly dynamic landscape, with multiple peaks showing a double increase in accessibility throughout development. We show that transcription factor spineless, but not ecdysone receptor, is highly predictive of opening sites Conclusions: This work provides a characterization of the chromatin dynamics of insect wing metamorphosis, identifies novel candidate chromatin remodeling factors in insects, and provides the first genome assembly of the model butterfly Junonia coenia, with gene and cis-regulatory element annotations
Project description:Background: We studied the chromatin accessibility landscsape in wings during butterfly metamorphosis, and investigate which transcription factors might be driving changes in accessibility Methods: We sequencing the Junonia coenia genome, and we studied chromatin accessibility using ATAC seq in multiple stages of wing development in both forewings and hindwings. For sites showing a large change in accessibility, we investigate which motifs are enriched, and correlate this with changes in gene expression of associated transcription factors. We confirm promising candidates with ChIP-seq Results: We find a highly dynamic landscape, with multiple peaks showing a double increase in accessibility throughout development. We show that transcription factor spineless, but not ecdysone receptor, is highly predictive of opening sites Conclusions: This work provides a characterization of the chromatin dynamics of insect wing metamorphosis, identifies novel candidate chromatin remodeling factors in insects, and provides the first genome assembly of the model butterfly Junonia coenia, with gene and cis-regulatory element annotations
Project description:Background: We studied the chromatin accessibility landscsape in wings during butterfly metamorphosis, and investigate which transcription factors might be driving changes in accessibility Methods: We sequencing the Junonia coenia genome, and we studied chromatin accessibility using ATAC seq in multiple stages of wing development in both forewings and hindwings. For sites showing a large change in accessibility, we investigate which motifs are enriched, and correlate this with changes in gene expression of associated transcription factors. We confirm promising candidates with ChIP-seq Results: We find a highly dynamic landscape, with multiple peaks showing a double increase in accessibility throughout development. We show that transcription factor spineless, but not ecdysone receptor, is highly predictive of opening sites Conclusions: This work provides a characterization of the chromatin dynamics of insect wing metamorphosis, identifies novel candidate chromatin remodeling factors in insects, and provides the first genome assembly of the model butterfly Junonia coenia, with gene and cis-regulatory element annotations
Project description:Heliconius butterfly wing pattern diversity offers a unique opportunity to investigate how natural genetic variation can drive the evolution of complex adaptive phenotypes. Here we took a large-scale transcriptomic approach to identify the network of genes involved in Heliconius wing pattern development and variation. This included applying 147 microarrays representing the Heliconius transcriptome to assay shifts in gene expression across pupal development among several wing pattern morphs of Heliconius erato. We focused in particular on genes differentially expressed relative to the gene optix, which controls red pattern elements in wings. We combined expression results from three hindwing morphs from Peru and from dissected basal to apical wing elements in two forewing morphs to uncover two main classes of genes. First we looked for candidate upstream regulators of optix by determining transcripts expressed differently across basal to apical sections of the forewing prior to optix expression. Second, we assessed how optix regulates downstream gene expression by targeting transcripts with differential expression similar to optix, where expression differs among red wing pattern elements of both the forewing and hindwing.
Project description:Heliconius butterfly wing pattern diversity offers a unique opportunity to investigate how natural genetic variation can drive the evolution of complex adaptive phenotypes. Here we took a large-scale transcriptomic approach to identify the network of genes involved in Heliconius wing pattern development and variation. This included applying 147 microarrays representing the Heliconius transcriptome to assay shifts in gene expression across pupal development among several wing pattern morphs of Heliconius erato. We focused in particular on genes differentially expressed relative to the gene optix, which controls red pattern elements in wings. We combined expression results from three hindwing morphs from Peru and from dissected basal to apical wing elements in two forewing morphs to uncover two main classes of genes. First we looked for candidate upstream regulators of optix by determining transcripts expressed differently across basal to apical sections of the forewing prior to optix expression. Second, we assessed how optix regulates downstream gene expression by targeting transcripts with differential expression similar to optix, where expression differs among red wing pattern elements of both the forewing and hindwing. This study is an analysis of two distinct datasets generated using the same microarray platform. One dataset involved comparative analysis of forewing sections of different color morphs, while the other compared whole hindwings with different color patterns. For the forewing analysis we compared proximal, medial, and distal wing sections of two color pattern morphs: H. erato petiverana and a hybrid H. himera x H. erato etylus. The proximal section in H. erato petiverana is black and the hybrid form orange/red, the medial section is red in H. erato petiverana and pale yellow in the hybrid form, and the distal section is black in both races. For the hindwing analysis, we compared hindwing color pattern gene expression in three races that meet in a hybrid zone in Peru. H. erato emma has a rayed hindwing, H. erato favorinus has a yellow-barred hindwing, and H. erato amphritrite has a black hindwing. Wings were dissected at five time intervals: 1, 3, and 5 days after pupation, when orange/red ommochrome pigments were beginning to be expressed (~7 days after pupation), and when black melanin pigments were starting to pepper the center of the wings (~8 days after pupation). In the forewings, Days 1, 3, and 5 were at 12, 36, and 60 hours post-pupation. In the hindwings these stages were sampled at 24, 48, and 72 hours. Samples hybridized to microarrays included three replicates each of each race, stage, and wing section for forewings (3 replicates x 2 morphs x 3 wing sections x 5 stages, with one replicate wing missing for Day 1 H. e. petiverana = 87 samples) and four replicates of each stage and race for hindwings (4 replicates x 3 races x 5 stages = 60 samples). Total RNA was extracted and converted to cDNA. Cy3-labeling of samples, hybridization, and array scanning was performed according to NimbleGen protocols (2008): for the forewings this was performed at the City of Hope Functional Genomics Core, while the hindwings were run separately at NCSU.Samples were hybridized to NimbleGen HD2 12-plex arrays. These arrays include 12 identical subarrays with 135,000 60 bp probes each, each hybridizing a separate sample. Samples were distributed across arrays to prevent repeat conditions as much as possible and to space similar conditions in different regions of the slide. The array design involved two classes of probes. First there was a tiling component involving 89,310 probes tiled across three genomic intervals. Results from the tiling data were used for the initial discovery of the optix gene and are not the focus of the present study. The second component involved a representation of a set of 12,450 transcript contigs at 1-6X coverage for a total of 40,046 probes, with a mean coverage of 3-4 probes per contig. The number of probes for each contig depended on the ability to create suitable probes according to NimbleGen probe selection criteria and was limited by the small size of some transcripts and the minimum spacing criterion of 15 bp apart. Sequences of low complexity and high repeats with the rest of the genome (>5X representation), determined by comparison against 1.6 MB of genomic sequence available at the time, were avoided for designing probes. An additional 3,248 random probes were placed on the array for quality control.
Project description:Butterfly wing patterns are an important model for studying the genetic basis of morphological evolution. Here we used RNA-seq expression profiling in the butterfly Vanessa cardui to characterize the transcriptional basis of wing pigmentation. This approach identified numerous candidate genes including known and suspected components of the insect melanin and ommochrome biosynthetic pathways.