RNA sequencing reveals differential thermal regulation mechanisms between sexes of Glanville fritillary butterfly in the Tianshan Mountains, China

The Glanville fritillary butterfly (Melitaea cinxia; Nymphalidae) has been extensively studied as a model species in metapopulation ecology. We investigated in the earlier studies that female butterflies exhibit higher thermal tolerance than males in the Tianshan Mountains of China. We aim to understand the molecular mechanism of differences of thermal responses between sexes. We used RNA-seq approach and performed de novo assembly of transcriptome to compare the gene expression patterns between two sexes after heat stress. All the reads were assembled into 84,376 transcripts and 72,701 unigenes. The number of differential expressed genes (DEGs) between control and heat shock samples was 175 and 268 for males and females, respectively. Heat shock proteins genes (hsps) were up-regulated in response to heat stress in both males and females. Most of the up-regulated hsps showed higher fold changes in males than in females. Females expressed more ribosomal subunit protein genes, transcriptional elongation factor genes, and methionine-rich storage protein genes, participating in protein synthesis. It indicated that protein synthesis is needed for females to replace the damaged proteins due to heat shock. In addition, aspartate decarboxylase might contribute to thermal tolerance in females. These differences in gene expression may at least partly explain the response to high temperature stress, and the fact that females exhibit higher thermal tolerance.