Transcriptome sequencing provides novel insights into larval development and sexual dimorphism in the firefly...

Fireflies are regarded as one of the most charismatic beetles due to their bioluminescence and ecological importance as bioindicators of freshwater quality. However, molecular mechanisms of larval development and sexual dimorphism in aquatic species remain poorly understood. Here, we performed multi-stage transcriptomic analysis of the aquatic firefly Aquatica leii across larval instars from L2 to L6, together with adult females and males, with three biological replicates per stage. Using time-series expression clustering, differential expression analysis, and weighted gene co-expression network analysis (WGCNA), we characterized the transcriptional dynamics of continuous larval development and the onset of sex-biased gene expression. We identified a critical transcriptional transition occurred at L5-L6, marked by downregulation of early morphogenetic genes and upregulation of juvenile hormone metabolism, oxidoreductase activity, and muscle contraction genes, indicating a shift from growth to metamorphic preparation. WGCNA identified a module strongly correlated with L6 (R = 0.97) enriched for the same functions, confirming a coordinated late-larval program. Notably, genes exhibiting sex-biased expression in adults were already expressed during late larval stages (L5 and L6), and 123 genes progressively upregulated from L2 to L6 showed enrichment in chitin biosynthesis, heart contraction, and ion transport; among these, six genes maintained high expression in adults with clear male-biased (Alei052192, Alei006658, and Alei087054) or female-biased (Alei003725, Alei096818, and Alei074026) patterns. These findings establish that transcriptional foundations for sexual dimorphism and adult tissue formation are laid during late larval stages, providing the first multi-stage transcriptomic resource for aquatic firefly conservation and breeding.