Seasonal phenotype-specific transcriptional reprogramming during metamorphosis in the European map butterfly Araschnia levana

Cereal weevils' antimicrobial peptides: at the crosstalk between development, endosymbiosis and immune response

Interactions between animals and microbes are ubiquitous in nature and strongly impact animal physiology. These interactions are shaped by the host immune system, which responds to infections and contributes to tailor the associations with beneficial microorganisms. In many insects, beneficial symbiotic associations not only include gut commensals, but also intracellular bacteria, or endosymbionts. Endosymbionts are housed within specialized host cells, the bacteriocytes, and are transmitted vertically across host generations. Host-endosymbiont co-evolution shapes the endosymbiont genome and host immune system, which not only fights against microbial intruders, but also ensures the preservation of endosymbionts and the control of their load and location. The cereal weevil Sitophilus spp. is a remarkable model in which to study the evolutionary adaptation of the immune system to endosymbiosis owing to its binary association with a unique, relatively recently acquired nutritional endosymbiont, Sodalis pierantonius. This Gram-negative bacterium has not experienced the genome size shrinkage observed in long-term endosymbioses and has retained immunogenicity. We focus here on the sixteen antimicrobial peptides (AMPs) identified in the Sitophilus oryzae genome and their expression patterns in different tissues, along host development or upon immune challenges, to address their potential functions in the defensive response and endosymbiosis homeostasis along the insect life cycle. This article is part of the theme issue 'Sculpting the microbiome: how host factors determine and respond to microbial colonization'.

A transcriptomic atlas underlying developmental plasticity of seasonal forms of Bicyclus anynana butterflies

Organisms residing in regions with alternating seasons often develop different phenotypes, or forms, in each season. These forms are often adaptations to each season and result from an altered developmental response to specific environmental cues such as temperature. While multiple studies have examined form-specific gene expression profiles in a diversity of species, little is known about how environments and developmental transitions, cued by hormone pulses, alter post-transcriptional patterns. In this study, we examine how gene expression, alternative splicing, and miRNA-mediated gene silencing in Bicyclus anynana butterfly hindwing tissue, varies across two rearing temperatures at four developmental timepoints. These timepoints flank two temperature-sensitive periods that coincide with two pulses of the insect hormone 20E. Our results suggest that developmental transitions, coincident with 20E pulses, elicit a greater impact on all these transcriptomic patterns than rearing temperatures per se. More similar transcriptomic patterns are observed pre-20E pulses than those observed post-20E pulses. We also found functionally distinct sets of differentially expressed and differentially spliced genes in the seasonal forms. Furthermore, around 10% of differentially expressed genes are predicted to be direct targets of, and regulated by, differentially expressed miRNAs between the seasonal forms. Many differentially expressed genes, miRNAs, or differentially spliced genes potentially regulate eyespot size plasticity, and we validated the differential splicing pattern of one such gene, daughterless. We present a comprehensive and interactive transcriptomic atlas of the hindwing tissue of both seasonal forms of B. anynana throughout development, a model organism of seasonal plasticity.

The European Map Butterfly Araschnia levana as a Model to Study the Molecular Basis and Evolutionary Ecology of Seasonal Polyphenism

The European map butterfly Araschnia levana is a well-known example of seasonal polyphenism. Spring and summer imagoes exhibit distinct morphological phenotypes. Key environmental factors responsible for the expression of different morphs are day length and temperature. Larval exposure to light for more than 16 h per day entails direct development and results in the adult f. prorsa summer phenotype. Less than 15.5 h per day increasingly promotes diapause and the adult f. levana spring phenotype. The phenotype depends on the timing of the release of 20-hydroxyecdysone in pupae. Release within the first days after pupation potentially inhibits the default "levana-gene-expression-profile" because pre-pupae destined for diapause or subitaneous development have unique transcriptomic programs. Moreover, multiple microRNAs and their targets are differentially regulated during the larval and pupal stages, and candidates for diapause maintenance, duration, and phenotype determination have been identified. However, the complete pathway from photoreception to timekeeping and diapause or subitaneous development remains unclear. Beside the wing polyphenism, the hormonal and epigenetic modifications of the two phenotypes also include differences in biomechanical design and immunocompetence. Here, we discuss research on the physiological and molecular basis of polyphenism in A. levana, including hormonal control, epigenetic regulation, and the effect of ecological parameters on developmental fate.

Ontogenetic Basis of Among-Generation Differences in Size-Related Traits in a Polyphenic Butterfly

Seasonal polyphenisms are cases in which individuals representing generations occurring in different times of the year systematically differ in their morphological, physiological, and/or behavioral traits. Such differences are often assumed to constitute adaptive responses to seasonally varying environments, but the evidence for this is still scarce. The adaptive character of the response would be corroborated by the pattern in which the decision about choosing a particular seasonal phenotype is made before the onset of respective environmental conditions (anticipatory plasticity). Alternatively, the between-generation differences can be caused by immediate effects of seasonally varying environments (responsive plasticity). Here we reared the larvae of the seasonally polymorphic map butterfly Araschnia levana under two different photoperiodic regimes, which provided different seasonal cues. These two treatments induced direct development and diapause pathways, respectively. Replicating the experiment at different temperatures and levels of host plant quality allowed us to evaluate both the anticipatory and the responsive components of the associated plastic changes in life-history traits. Larvae representing the direct development pathway invariably had higher growth rates and shorter development periods, although the difference between the developmental pathways was smaller at inferior host quality. Body size differences between the developmental pathways turned out to be less consistent, as the natural pattern of higher pupal mass of the directly developing individuals could only be reproduced at lower rearing temperature. Though being considerably modified by immediate environmental effects, the between-generation differences in size, growth rates, and larval are largely based on anticipatory plasticity (= responses to photoperiodic cues) and should be treated as seasonal adaptations in A. levana. In a more general context, we show how investigating the proximate basis of size differences can serve the purpose of identifying the limits of phenotypic plasticity in juvenile growth schedules.

The complete mitochondrial genome and phylogenetic analysis of the European map butterfly Araschnia levana (Insecta: Lepidoptera: Nymphalidae)

The European map butterfly Araschnia levana (Linnaeus, 1758) is a species showing extreme seasonal polyphenism. The complete 15,207 bp circular A. levana mitogenome consisting of 81.6% AT nucleotides, was assembled by Illumina genome skimming. It includes 22 tRNAs, 13 protein-coding genes, 2 rRNAs, and a control region in the typical butterfly gene order. Araschnia levana COX1 features an atypical CGA start codon and ATP6, COX1, COX2, ND1, ND3, and ND4 have incomplete stop codons completed by 3′A residues added to the mRNA. Phylogenetic reconstruction places A. levana as a basal lineage within tribe Nymphalini, consistent with previous phylogenetic hypotheses.

Functional insights from the GC-poor genomes of two aphid parasitoids, Aphidius ervi and Lysiphlebus fabarum

BACKGROUND

Parasitoid wasps have fascinating life cycles and play an important role in trophic networks, yet little is known about their genome content and function. Parasitoids that infect aphids are an important group with the potential for biological control. Their success depends on adapting to develop inside aphids and overcoming both host aphid defenses and their protective endosymbionts.

RESULTS

We present the de novo genome assemblies, detailed annotation, and comparative analysis of two closely related parasitoid wasps that target pest aphids: Aphidius ervi and Lysiphlebus fabarum (Hymenoptera: Braconidae: Aphidiinae). The genomes are small (139 and 141 Mbp) and the most AT-rich reported thus far for any arthropod (GC content: 25.8 and 23.8%). This nucleotide bias is accompanied by skewed codon usage and is stronger in genes with adult-biased expression. AT-richness may be the consequence of reduced genome size, a near absence of DNA methylation, and energy efficiency. We identify missing desaturase genes, whose absence may underlie mimicry in the cuticular hydrocarbon profile of L. fabarum. We highlight key gene groups including those underlying venom composition, chemosensory perception, and sex determination, as well as potential losses in immune pathway genes.

CONCLUSIONS

These findings are of fundamental interest for insect evolution and biological control applications. They provide a strong foundation for further functional studies into coevolution between parasitoids and their hosts. Both genomes are available at https://bipaa.genouest.org.

Seasonal phenotype-specific expression of microRNAs during metamorphosis in the European map butterfly Araschnia levana

The European map butterfly (Araschnia levana) is a well-known example of seasonal polyphenism because the spring and summer imagoes exhibit distinct morphological phenotypes. The day length and temperature during larval and prepupal development determine whether spring or summer imagoes emerge after metamorphosis. Inspired by the fundamentally different transcriptomic profiles in prepupae developing from larvae exposed to long days or short days, we postulate that posttranscriptional epigenetic regulators such as microRNAs (miRNAs) may contribute to the epigenetic control of seasonal polyphenism in A. levana. To test this hypothesis, we used microarrays containing over 2,000 insect miRNAs to identify candidate regulators that are differentially expressed in last-instar larvae or pupae developing under long-day or short-day conditions. We used our transcriptomic database to identify potential 3'-untranslated regions of messenger RNAs to predict miRNA targets by considering both base pair complementarity and minimum free energy hybridization. This approach resulted in the identification of multiple targets of miRNAs that were differentially regulated in polyphenic morphs of A. levana including a candidate (miR-2856-3p) regulating the previously identified diapause bioclock protein gene. In conclusion, the expression profiling of miRNAs provided insights into their possible involvement in seasonal polyphenism of A. levana and offer an important resource for further studies.

Silencing of the DNA methyltransferase 1 associated protein 1 (DMAP1) gene in the invasive ladybird Harmonia axyridis implies a role of the DNA methyltransferase 1-DMAP1 complex in female fecundity

The invasive harlequin ladybird Harmonia axyridis is a textbook example of polymorphism and polyphenism as the temperature during egg development determines the frequency of melanic morphs and the number and size of black spots in nonmelanic morphs. Recent concepts in evolutionary biology suggest that epigenetic mechanisms can translate environmental stimuli into heritable phenotypic changes. To investigate whether epigenetic mechanisms influence the penetrance and expressivity of colour morphs in H. axyridis, we used RNA interference to silence key enzymes required for DNA methylation and histone modification. We found that neither of these epigenetic mechanisms affected the frequency of different morphs, but there was a significant impact on life-history traits such as longevity and fecundity. Strikingly, we found that silencing the gene encoding for DNA methyltransferase 1 associated protein 1 (DMAP1) severely reduced female fecundity, which correlated with an abundance of degenerated ovaries in DMAP1-knockdown female beetles. Finally, we observed significant differences in DMAP1 expression when we compared native and invasive H. axyridis populations with a biocontrol strain differing in egg-laying capacity, suggesting that the DNA methyltransferase 1-DMAP1 complex may influence the invasive performance of this ladybird.

Immunological larval polyphenism in the map butterfly Araschnia levana reveals the photoperiodic modulation of immunity

The bivoltine European map butterfly (Araschnia levana) displays seasonal polyphenism characterized by the formation of two remarkably distinct dorsal wing phenotypes: The spring generation (A. levana levana) is predominantly orange with black spots and develops from diapause pupae, whereas the summer generation (A. levana prorsa) has black, white, and orange bands and develops from subitaneous pupae. The choice between spring or summer imagoes is regulated by the photoperiod during larval and prepupal development, but polyphenism in the larvae has not been investigated before. Recently, it has been found that the prepupae of A. levana display differences in immunity-related gene expression, so we tested whether larvae destined to become spring (short-day) or summer (long-day) morphs also display differences in innate immunity. We measured larval survival following the injection of a bacterial entomopathogen (Pseudomonas entomophila), the antimicrobial activity in their hemolymph and the induced expression of selected genes encoding antimicrobial peptides (AMPs). Larvae of the short-day generation died significantly later, exhibited higher antibacterial activity in the hemolymph, and displayed higher induced expression levels of AMPs than those of the long-day generation. Our study expands the seasonal polyphenism of A. levana beyond the morphologically distinct spring and summer imagoes to include immunological larval polyphenism that reveals the photoperiodic modulation of immunity. This may reflect life-history traits that manifest as trade-offs between immunity and fecundity.

Conserved roles of Osiris genes in insect development, polymorphism and protection

Much of the variation among insects is derived from the different ways that chitin has been moulded to form rigid structures, both internal and external. In this study, we identify a highly conserved expression pattern in an insect-only gene family, the Osiris genes, that is essential for development, but also plays a significant role in phenotypic plasticity and in immunity/toxicity responses. The majority of Osiris genes exist in a highly syntenic cluster, and the cluster itself appears to have arisen very early in the evolution of insects. We used developmental gene expression in the fruit fly, Drosophila melanogaster, the bumble bee, Bombus terrestris, the harvester ant, Pogonomyrmex barbatus, and the wood ant, Formica exsecta, to compare patterns of Osiris gene expression both during development and between alternate caste phenotypes in the polymorphic social insects. Developmental gene expression of Osiris genes is highly conserved across species and correlated with gene location and evolutionary history. The social insect castes are highly divergent in pupal Osiris gene expression. Sets of co-expressed genes that include Osiris genes are enriched in gene ontology terms related to chitin/cuticle and peptidase activity. Osiris genes are essential for cuticle formation in both embryos and pupae, and genes co-expressed with Osiris genes affect wing development. Additionally, Osiris genes and those co-expressed seem to play a conserved role in insect toxicology defences and digestion. Given their role in development, plasticity, and protection, we propose that the Osiris genes play a central role in insect adaptive evolution.

Oral immune priming with Bacillus thuringiensis induces a shift in the gene expression of Tribolium castaneum larvae

Background

The phenomenon of immune priming, i.e. enhanced protection following a secondary exposure to a pathogen, has now been demonstrated in a wide range of invertebrate species. Despite accumulating phenotypic evidence, knowledge of its mechanistic underpinnings is currently very limited. Here we used the system of the red flour beetle, Tribolium castaneum and the insect pathogen Bacillus thuringiensis (Bt) to further our molecular understanding of the oral immune priming phenomenon. We addressed how ingestion of bacterial cues (derived from spore supernatants) of an orally pathogenic and non-pathogenic Bt strain affects gene expression upon later challenge exposure, using a whole-transcriptome sequencing approach.

Results

Whereas gene expression of individuals primed with the orally non-pathogenic strain showed minor changes to controls, we found that priming with the pathogenic strain induced regulation of a large set of distinct genes, many of which are known immune candidates. Intriguingly, the immune repertoire activated upon priming and subsequent challenge qualitatively differed from the one mounted upon infection with Bt without previous priming. Moreover, a large subset of priming-specific genes showed an inverse regulation compared to their regulation upon challenge only.

Conclusions

Our data demonstrate that gene expression upon infection is strongly affected by previous immune priming. We hypothesise that this shift in gene expression indicates activation of a more targeted and efficient response towards a previously encountered pathogen, in anticipation of potential secondary encounter.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3705-7) contains supplementary material, which is available to authorized users.