Understanding the fate of the hungry caterpillar: Transcriptional and DNA methylation programs during metamorphosis of the butterfly Bicyclus anynana (RNA-Seq).
Ontology highlight
ABSTRACT:
ORGANISM(S): Bicyclus anynana
PROVIDER: GSE276477 | GEO | 2025/10/13
REPOSITORIES: GEO
Similar Datasets
Project description:An immobile pupal phase in butterfly metamorphosis is central for the production of an airborne adult from the earlier caterpillar feeding stage. Extensive body plan re-organization occurs within the pupal case that is entirely dependent upon nutrients accumulated during caterpillar growth, and sensitive to external environmental conditions. Here, we use DNA and RNA sequencing of a laboratory butterfly population to produce a time course dataset describing changes in DNA methylation and in the transcriptome during progression from the 5th instar caterpillar larval stage through to the latest stage of pupation when the mature butterfly is about to eclose from the pupal case.
2025-10-13 | GSE276476 | GEO
Project description:Disrupting butterfly caterpillar microbiomes does not impact their survival and development
| PRJEB35305 | ENA
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
2019-05-07 | GSE121542 | GEO
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
2019-04-30 | GSE121541 | GEO
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
2019-04-30 | GSE121734 | GEO
Project description:Butterfly wings display a diversity of cell types, including large polyploid scale cells, yet the molecular basis of such diversity is poorly understood. To explore scale cell diversity at a transcriptomic level, we employ single-cell RNA sequencing of ∼5,200 large cells (>6 μm) from 22.5- to 25-h male pupal forewings of the butterfly Bicyclus anynana. Using unsupervised clustering, followed by in situ hybridization, immunofluorescence, and CRISPR-Cas9 editing of candidate genes, we annotate various cell types on the wing. We identify genes marking non-innervated scale cells, pheromone-producing glandular cells, and innervated sensory cell types. We show that senseless, a zinc-finger transcription factor, and HR38, a hormone receptor, determine the identity, size, and color of different scale cell types and are important regulators of scale cell differentiation. This dataset and the identification of various wing cell-type markers provide a foundation to compare and explore scale cell-type diversification across arthropod species.
2024-04-25 | GSE241536 | GEO
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.
2012-06-30 | E-GEOD-38084 | biostudies-arrayexpress
Project description:We investigated gene expression levels in Heliconius erato butterflies with divergent wing patterns across a 656KB genomic interval linked to the red color pattern wing polymorphism. This included comparison of expression between two H. erato color pattern populations (H. e. petiverana and a H.e. etylus x H. himera hybrid) across three sections of the forewing that differed in pigmentation (the basal, mid, and distal wing sections) and five different stages of pupal development (Day 1, 3, 5 pupae and ommochrome and melanin pigmentation stages). These results allowed us to determine whether certain genes in this interval were differentially expressed between the wing pattern elements, and, therefore, potentially responsible for adaptive color pattern variation in these butterflies. Forewings from a total of 29 individuals, covering three biological replicates of five developmental time points for each of the two H. erato distinct phenotypes were dissected, with the only exception being that there were only two replicates of the day 1 hybrid phenotype. Individuals were reared at 25˚C and dissected at the following stages: a) day 1 = 12 hr after pupation; b) day 3 = 60 hr after pupation; c) day 5 = 108 hr after pupation; d) early ommochrome = ~ 156 hr after pupation when red scales in forewing partially mature, showing a pale orange color; and e) early melanin = ~ 180 hr after pupation melanic scales begin to turn black and are present primarily at the center of the wing. Using wing veins as landmarks, each forewing was cut into three sections corresponding to the color pattern boundaries: basal (F1), middle (F2), and apical (F3). A eight custom-designed Roche NimbleGen 12x135K format microarrays with probes spanning a 656,307bp genomic region (Roche NimbleGen Inc., Madison, Wisconsin, United States) were used to hybridize double stranded cDNA from 87 tissue samples. Repetitive sequence elements found more than five times across all currently available H. erato genomic sequences, including the probed region as well as additional genomic BAC sequences, were masked from the tiling region. The remaining unmasked non-repetitive genomic sequences were tiled using 60 bp probes staggered every 13 bp on average, with slight modifications to ensure probe quality, for a total of 48,547 probes. Microarry design and printing was performed by Roche NimbleGen. cDNA labeling, hybridization, and array scanning was performed by the City of Hope Microarray Facility (Duarte, California, United States). In addition to the probes from the red color pattern intervals, the arrays also include 40,763 probes across two other genomic intervals not addressed in this study, 45,046 probes representing 12450 transcripts from a recent transcriptome assembly at 1-6X coverage, and 3248 random probes. Results from the transcriptome and other color pattern intervals will be published separately, however, we analyzed all probes together for array normalization and quality control.
2011-12-02 | E-GEOD-30221 | biostudies-arrayexpress
Project description:Limacodidae is a family of lepidopteran insects comprising more than 1,500 species, over half of which produce pain-inducing defensive venom components in the larval stage, but little is known about their venom. Recently, we characterised proteinaceous toxins from the Australian limacodid caterpillar Doratifera vulnerans, but it is unknown if the venom of this species is typical of other Limacodidae. Here, we use single animal transcriptomics and venom proteomics to investigate the venom of an iconic limacodid, the North American saddleback caterpillar Acharia stimulea. We identified 65 confirmed venom polypeptides, grouped into 31 different families. Neurohormones, knottins, and homologues of the immune signaller Diedel make up the majority of A. stimulea venom, indicating strong similarities to D. vulnerans venom, despite their geographic separation. One notable difference is the presence of RF-amide peptide toxins in A. stimulea venom. Synthetic versions of these RF-amide toxins displayed potent activation of the human neuropeptide FF1 receptor, insecticidal activity when injected into Drosophila melanogaster, and moderate ability to inhibit the larval development of the parasitic nematode Haemonchus contortus. This study provides insights into venom evolution and toxin activity in Limacodidae, and will allow the synthetic production of A. stimulea peptide toxins towards functional and structural characterisation.
2023-06-01 | PXD039980 | Pride
Project description:Expression of activated Ras along with alteration in levels of non-coding RNA hsromega in the eye disc results in lowered ecdysone levels in post pupation stages, resulting in early pupal death. This microarray analysis of total pupal RNAs was undertaken to know how could expression of activated Ras with altered levels of hsromega RNA in eye discs could reduce ecdysone levels
2016-04-28 | GSE80703 | GEO