Gene Duplication and Gene Expression Changes Play a Role in the Evolution of Candidate Pollen Feeding Genes in Heliconius Butterflies

Colour polymorphism associated with a gene duplication in male wood tiger moths

Colour is often used as an aposematic warning signal, with predator learning expected to lead to a single colour pattern within a population. However, there are many puzzling cases where aposematic signals are also polymorphic. The wood tiger moth, Arctia plantaginis, displays bright hindwing colours associated with unpalatability, and males have discrete colour morphs which vary in frequency between localities. In Finland, both white and yellow morphs can be found, and these colour morphs also differ in behavioural and life-history traits. Here, we show that male colour is linked to an extra copy of a yellow family gene that is only present in the white morphs. This white-specific duplication, which we name valkea, is highly upregulated during wing development. CRISPR targeting valkea resulted in editing of both valkea and its paralog, yellow-e, and led to the production of yellow wings. We also characterise the pigments responsible for yellow, white, and black colouration, showing that yellow is partly produced by pheomelanins, while black is dopamine-derived eumelanin. Our results add to a growing number of studies on the genetic architecture of complex and seemingly paradoxical polymorphisms, and the role of gene duplications and structural variation in adaptive evolution.

takeout gene expression is associated with temporal kin recognition

A key component of parental care is avoiding killing and eating one's own offspring. Many organisms commit infanticide but switch to parental care when their own offspring are expected, known as temporal kin recognition. It is unclear why such types of indirect kin recognition are so common across taxa. One possibility is that temporal kin recognition may evolve through alteration of simple mechanisms, such as co-opting mechanisms that influence the regulation of timing and feeding in other contexts. Here, we determine whether takeout, a gene implicated in coordinating feeding, influences temporal kin recognition in Nicrophorus orbicollis. We found that takeout expression was not associated with non-parental feeding changes resulting from hunger, or a general transition to the full parental care repertoire. However, beetles that accepted and provided care to their offspring had a higher takeout expression than beetles that committed infanticide. Together, these data support the idea that the evolution of temporal kin recognition may be enabled by co-option of mechanisms that integrate feeding behaviour in other contexts.

In silico structural and functional characterization of Antheraea mylitta cocoonase

Background

Cocoonase is a serine protease present in sericigenous insects and majorly involved in dissolving of sericin protein allowing moth to escape. Cocoon structure is made up of sericin protein which holds fibroin filaments together. Cocoonase enzyme hydrolyzes sericin protein without harming the fibroin. However, until date, no detailed characterization of cocoonase enzyme and its presence in wild silk moth Antheraea mylitta has been carried out. Therefore, current study aimed for detailed characterization of amplified cocoonase enzyme, secondary and tertiary structure prediction, sequence and structural alignment, phylogenetic analysis, and computational validation. Several computational tools such as ProtParam, Iterative Threading Assembly Refinement (I-TASSER), PROCHECK, SAVES v6.0, TM-align, Molecular Evolutionary Genetics Analysis (MEGA) X, and Figtree were employed for characterization of cocoonase protein.

Results

The present study elucidates about the isolation of RNA, cDNA preparation, PCR amplification, and in silico characterization of cocoonase from Antheraea mylitta. Here, total RNA was isolated from head region of A. mylitta, and gene-specific primers were designed using Primer3 followed by PCR-based amplification and sequencing. The newly constructed 377-bp length sequence of cocoonase was subjected to in silico characterization. In silico study of A. mylitta cocoonase showed 26% similarity to A. pernyi strain Qing-6 cocoonase using Blastp and belongs to member of chymotrypsin-like serine protease superfamily. From phylogenetic study, it was found that A. mylitta cocoonase sequence is closely related to A. pernyi cocoonase sequence.

Conclusions

The present study revealed about the detailed in silico characterization of cocoonase gene and encoded protein obtained from A. mylitta head region. The results obtained infer the presence of cocoonase enzyme in the wild silkworm A. mylitta and can be used for cocoon degumming which will be a valuable and cost-effective strategy in silk industry.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43141-022-00367-8.

Heliconius Butterflies Host Characteristic and Phylogenetically Structured Adult-Stage Microbiomes

Lepidoptera (butterflies and moths) are diverse and ecologically important, yet we know little about how they interact with microbes as adults. Due to metamorphosis, the form and function of their adult-stage microbiomes might be very different from those of microbiomes in the larval stage (caterpillars). We studied adult-stage microbiomes of Heliconius and closely related passion-vine butterflies (Heliconiini), which are an important model system in evolutionary biology. To characterize the structure and dynamics of heliconiine microbiomes, we used field collections of wild butterflies, 16S rRNA gene sequencing, quantitative PCR, and shotgun metagenomics. We found that Heliconius butterflies harbor simple and abundant bacterial communities that are moderately consistent among conspecific individuals and over time. Heliconiine microbiomes also exhibited a strong signal of the host phylogeny, with a major distinction between Heliconius and other butterflies. These patterns were largely driven by differing relative abundances of bacterial phylotypes shared among host species and genera, as opposed to the presence or absence of host-specific phylotypes. We suggest that the phylogenetic structure in heliconiine microbiomes arises from conserved host traits that differentially filter microbes from the environment. While the relative importance of different traits remains unclear, our data indicate that pollen feeding (unique to Heliconius) is not a primary driver. Using shotgun metagenomics, we also discovered trypanosomatids and microsporidia to be prevalent in butterfly guts, raising the possibility of antagonistic interactions between eukaryotic parasites and colocalized gut bacteria. Our discovery of characteristic and phylogenetically structured microbiomes provides a foundation for tests of adult-stage microbiome function, a poorly understood aspect of lepidopteran biology.IMPORTANCE Many insects host microbiomes with important ecological functions. However, the prevalence of this phenomenon is unclear because in many insect taxa, microbiomes have been studied in only part of the life cycle, if at all. A prominent example is butterflies and moths, in which the composition and functional role of adult-stage microbiomes are largely unknown. We comprehensively characterized microbiomes in adult passion-vine butterflies. Butterfly-associated bacterial communities are generally abundant in guts, consistent within populations, and composed of taxa widely shared among hosts. More closely related butterflies harbor more similar microbiomes, with the most dramatic shift in microbiome composition occurring in tandem with a suite of ecological and life history traits unique to the genus Heliconius Butterflies are also frequently infected with previously undescribed eukaryotic parasites, which may interact with bacteria in important ways. These findings advance our understanding of butterfly biology and insect-microbe interactions generally.

Pollen feeding in Heliconius butterflies: the singular evolution of an adaptive suite

Major evolutionary transitions can be triggered by behavioural novelty, and are often associated with 'adaptive suites', which involve shifts in multiple co-adapted traits subject to complex interactions. Heliconius butterflies represent one such example, actively feeding on pollen, a behaviour unique among butterflies. Pollen feeding permits a prolonged reproductive lifespan, and co-occurs with a constellation of behavioural, neuroanatomical, life history, morphological and physiological traits that are absent in closely related, non-pollen-feeding genera. As a highly tractable system, supported by considerable ecological and genomic data, Heliconius are an excellent model for investigating how behavioural innovation can trigger a cascade of adaptive shifts in multiple diverse, but interrelated, traits. Here, we synthesize current knowledge of pollen feeding in Heliconius, and explore potential interactions between associated, putatively adaptive, traits. Currently, no physiological, morphological or molecular innovation has been explicitly linked to the origin of pollen feeding, and several hypothesized links between different aspects of Heliconius biology remain poorly tested. However, resolving these uncertainties will contribute to our understanding of how behavioural innovations evolve and subsequently alter the evolutionary trajectories of diverse traits impacting resource acquisition, life history, senescence and cognition.

Cocoonase is indispensable for Lepidoptera insects breaking the sealed cocoon

Many insects spin cocoons to protect the pupae from unfavorable environments and predators. After emerging from the pupa, the moths must escape from the sealed cocoons. Previous works identified cocoonase as the active enzyme loosening the cocoon to form an escape-hatch. Here, using bioinformatics tools, we show that cocoonase is specific to Lepidoptera and that it probably existed before the occurrence of lepidopteran insects spinning cocoons. Despite differences in cocooning behavior, we further show that cocoonase evolved by purification selection in Lepidoptera and that the selection is more intense in lepidopteran insects spinning sealed cocoons. Experimentally, we applied gene editing techniques to the silkworm Bombyx mori, which spins a dense and sealed cocoon, as a model of lepidopteran insects spinning sealed cocoons. We knocked out cocoonase using the CRISPR/Cas9 system. The adults of homozygous knock-out mutants were completely formed and viable but stayed trapped and died naturally in the cocoon. This is the first experimental and phenotypic evidence that cocoonase is the determining factor for breaking the cocoon. This work led to a novel silkworm strain yielding permanently intact cocoons and provides a new strategy for controlling the pests that form cocoons.

Spinning and cocooning are the instincts of many insects, providing a shelter to the residing pupae to resist adverse factors. After the metamorphosis of pupa into adult, the adult must break open the cocoon to emerge, which is called decocooning. We have bioinformatically identified that cocoonase is specific to Lepidoptera, and demonstrated that it is the determining factor for breaking the sealed cocoon experimentally for the first time. This work led to a novel silkworm strain yielding permanently intact cocoons and provides a new strategy for controlling the pests that form cocoons and for breeding.

Prokaryotic Genome Expansion Is Facilitated by Phages and Plasmids but Impaired by CRISPR

Viruses and plasmids can introduce novel DNA into bacterial cells, thereby creating an opportunity for genome expansion; conversely, CRISPR, the prokaryotic adaptive immune system, which targets and eliminates foreign DNAs, may impair genome expansions. Recent studies presented conflicting results over the impact of CRISPR on genome expansion. In this study, we constructed a comprehensive dataset of prokaryotic genomes and identified their associations with viruses and plasmids. We found that genomes associated with viruses and/or plasmids were significantly larger than those without, indicating that both viruses and plasmids contribute to genome expansion. Genomes were increasingly larger with increasing numbers of associated viruses or plasmids. Conversely, genomes with CRISPR systems were significantly smaller than those without, indicating that CRISPR has a negative impact on genome size. These results confirmed that on evolutionary timescales, viruses and plasmids facilitate genome expansion, while CRISPR impairs such a process in prokaryotes. Furthermore, our results also revealed that CRISPR systems show a preference for targeting viruses over plasmids.

Drift and Directional Selection Are the Evolutionary Forces Driving Gene Expression Divergence in Eye and Brain Tissue of Heliconius Butterflies

Investigating gene expression evolution over micro- and macroevolutionary timescales will expand our understanding of the role of gene expression in adaptation and speciation. In this study, we characterized the evolutionary forces acting on gene expression levels in eye and brain tissue of five Heliconius butterflies with divergence times of ∼5-12 MYA. We developed and applied Brownian motion (BM) and Ornstein-Uhlenbeck (OU) models to identify genes whose expression levels are evolving through drift, stabilizing selection, or a lineage-specific shift. We found that 81% of the genes evolve under genetic drift. When testing for branch-specific shifts in gene expression, we detected 368 (16%) shift events. Genes showing a shift toward upregulation have significantly lower gene expression variance than those genes showing a shift leading toward downregulation. We hypothesize that directional selection is acting in shifts causing upregulation, since transcription is costly. We further uncovered through simulations that parameter estimation of OU models is biased when using small phylogenies and only becomes reliable with phylogenies having ≥ 50 taxa. Therefore, we developed a new statistical test based on BM to identify highly conserved genes (i.e., evolving under strong stabilizing selection), which comprised 3% of the orthoclusters. In conclusion, we found that drift is the dominant evolutionary force driving gene expression evolution in eye and brain tissue in Heliconius Nevertheless, the higher proportion of genes evolving under directional than under stabilizing selection might reflect species-specific selective pressures on vision and the brain that are necessary to fulfill species-specific requirements.

Evolutionary and structural analyses uncover a role for solvent interactions in the diversification of cocoonases in butterflies

Multi-omic approaches promise to supply the power to detect genes underlying disease and fitness-related phenotypes. Optimal use of the resulting profusion of data requires detailed investigation of individual candidate genes, a challenging proposition. Here, we combine transcriptomic and genomic data with molecular modelling of candidate enzymes to characterize the evolutionary history and function of the serine protease cocoonase. Heliconius butterflies possess the unique ability to feed on pollen; recent work has identified cocoonase as a candidate gene in pollen digestion. Cocoonase was first described in moths, where it aids in eclosure from the cocoon and is present as a single copy gene. In heliconiine butterflies it is duplicated and highly expressed in the mouthparts of adults. At least six copies of cocoonase are present in Heliconius melpomene and copy number varies across H. melpomene sub-populations. Most cocoonase genes are under purifying selection, however branch-site analyses suggest cocoonase 3 genes may have evolved under episodic diversifying selection. Molecular modelling of cocoonase proteins and examination of their predicted structures revealed that the active site region of each type has a similar structure to trypsin, with the same predicted substrate specificity across types. Variation among heliconiine cocoonases instead lies in the outward-facing residues involved in solvent interaction. Thus, the neofunctionalization of cocoonase duplicates appears to have resulted from the need for these serine proteases to operate in diverse biochemical environments. We suggest that cocoonase may have played a buffering role in feeding during the diversification of Heliconius across the neotropics by enabling these butterflies to digest protein from a range of biochemical milieux.

Genome-wide patterns of copy number variations in Spodoptera litura

Spodoptera litura is a polyphagous pest and can feed on more than 100 species of plants, causing great damage to agricultural production. The SNP results showed that there were gene exchanges between different regions. To explore the variations of larger segments in S. litura genome, we used genome resequencing samples from 14 regions of China, India, and Japan to study the copy number variations (CNVs). We identified 3976 CNV events and 1581 unique copy number variation regions (CNVRs) occupying the 108.5 Mb genome of S. litura. A total of 5527 genes that overlapped with CNVRs were detected. Selection signal analysis identified 19 shared CNVRs and 105 group-specific CNVRs, whose related genes were involved in various biological processes in S. litura. We constructed the first CNVs map in S. litura genome, and our findings will be valuable for understanding the genomic variations and population differences of S. litura.

Copy Number Variation and Expression Analysis Reveals a Nonorthologous Pinta Gene Family Member Involved in Butterfly Vision

Vertebrate (cellular retinaldehyde-binding protein) and Drosophila (prolonged depolarization afterpotential is not apparent [PINTA]) proteins with a CRAL-TRIO domain transport retinal-based chromophores that bind to opsin proteins and are necessary for phototransduction. The CRAL-TRIO domain gene family is composed of genes that encode proteins with a common N-terminal structural domain. Although there is an expansion of this gene family in Lepidoptera, there is no lepidopteran ortholog of pinta. Further, the function of these genes in lepidopterans has not yet been established. Here, we explored the molecular evolution and expression of CRAL-TRIO domain genes in the butterfly Heliconius melpomene in order to identify a member of this gene family as a candidate chromophore transporter. We generated and searched a four tissue transcriptome and searched a reference genome for CRAL-TRIO domain genes. We expanded an insect CRAL-TRIO domain gene phylogeny to include H. melpomene and used 18 genomes from 4 subspecies to assess copy number variation. A transcriptome-wide differential expression analysis comparing four tissue types identified a CRAL-TRIO domain gene, Hme CTD31, upregulated in heads suggesting a potential role in vision for this CRAL-TRIO domain gene. RT-PCR and immunohistochemistry confirmed that Hme CTD31 and its protein product are expressed in the retina, specifically in primary and secondary pigment cells and in tracheal cells. Sequencing of eye protein extracts that fluoresce in the ultraviolet identified Hme CTD31 as a possible chromophore binding protein. Although we found several recent duplications and numerous copy number variants in CRAL-TRIO domain genes, we identified a single copy pinta paralog that likely binds the chromophore in butterflies.