Heliconius Melpomene Research Articles

Major Improvements to the Heliconius melpomene Genome Assembly Used to Confirm 10 Chromosome Fusion Events in 6 Million Years of Butterfly Evolution.

The Heliconius butterflies are a widely studied adaptive radiation of 46 species spread across Central and South America, several of which are known to hybridize in the wild. Here, we present a substantially improved assembly of the Heliconius melpomene genome, developed using novel methods that should be applicable to improving other genome assemblies produced using short read sequencing. First, we whole-genome-sequenced a pedigree to produce a linkage map incorporating 99% of the genome. Second, we incorporated haplotype scaffolds extensively to produce a more complete haploid version of the draft genome. Third, we incorporated ∼20x coverage of Pacific Biosciences sequencing, and scaffolded the haploid genome using an assembly of this long-read sequence. These improvements result in a genome of 795 scaffolds, 275 Mb in length, with an N50 length of 2.1 Mb, an N50 number of 34, and with 99% of the genome placed, and 84% anchored on chromosomes. We use the new genome assembly to confirm that the Heliconius genome underwent 10 chromosome fusions since the split with its sister genus Eueides, over a period of about 6 million yr.

Evolutionary Novelty in a Butterfly Wing Pattern through Enhancer Shuffling.

An important goal in evolutionary biology is to understand the genetic changes underlying novel morphological structures. We investigated the origins of a complex wing pattern found among Amazonian Heliconius butterflies. Genome sequence data from 142 individuals across 17 species identified narrow regions associated with two distinct red colour pattern elements, dennis and ray. We hypothesise that these modules in non-coding sequence represent distinct cis-regulatory loci that control expression of the transcription factor optix, which in turn controls red pattern variation across Heliconius. Phylogenetic analysis of the two elements demonstrated that they have distinct evolutionary histories and that novel adaptive morphological variation was created by shuffling these cis-regulatory modules through recombination between divergent lineages. In addition, recombination of modules into different combinations within species further contributes to diversity. Analysis of the timing of diversification in these two regions supports the hypothesis of introgression moving regulatory modules between species, rather than shared ancestral variation. The dennis phenotype introgressed into Heliconius melpomene at about the same time that ray originated in this group, while ray introgressed back into H. elevatus much more recently. We show that shuffling of existing enhancer elements both within and between species provides a mechanism for rapid diversification and generation of novel morphological combinations during adaptive radiation.

Efficient Strategies for Calculating Blockwise Likelihoods Under the Coalescent.

The inference of demographic history from genome data is hindered by a lack of efficient computational approaches. In particular, it has proved difficult to exploit the information contained in the distribution of genealogies across the genome. We have previously shown that the generating function (GF) of genealogies can be used to analytically compute likelihoods of demographic models from configurations of mutations in short sequence blocks (Lohse et al. 2011). Although the GF has a simple, recursive form, the size of such likelihood calculations explodes quickly with the number of individuals and applications of this framework have so far been mainly limited to small samples (pairs and triplets) for which the GF can be written by hand. Here we investigate several strategies for exploiting the inherent symmetries of the coalescent. In particular, we show that the GF of genealogies can be decomposed into a set of equivalence classes that allows likelihood calculations from nontrivial samples. Using this strategy, we automated blockwise likelihood calculations for a general set of demographic scenarios in Mathematica. These histories may involve population size changes, continuous migration, discrete divergence, and admixture between multiple populations. To give a concrete example, we calculate the likelihood for a model of isolation with migration (IM), assuming two diploid samples without phase and outgroup information. We demonstrate the new inference scheme with an analysis of two individual butterfly genomes from the sister species Heliconius melpomene rosina and H. cydno.

Pollen feeding proteomics: Salivary proteins of the passion flower butterfly, Heliconius melpomene

While most adult Lepidoptera use flower nectar as their primary food source, butterflies in the genus Heliconius have evolved the novel ability to acquire amino acids from consuming pollen. Heliconius butterflies collect pollen on their proboscis, moisten the pollen with saliva, and use a combination of mechanical disruption and chemical degradation to release free amino acids that are subsequently re-ingested in the saliva. Little is known about the molecular mechanisms of this complex pollen feeding adaptation. Here we report an initial shotgun proteomic analysis of saliva from Heliconius melpomene. Results from liquid-chromatography tandem mass-spectrometry confidently identified 31 salivary proteins, most of which contained predicted signal peptides, consistent with extracellular secretion. Further bioinformatic annotation of these salivary proteins indicated the presence of four distinct functional classes: proteolysis (10 proteins), carbohydrate hydrolysis (5), immunity (6), and “housekeeping” (4). Additionally, six proteins could not be functionally annotated beyond containing a predicted signal sequence. The presence of several salivary proteases is consistent with previous demonstrations that Heliconius saliva has proteolytic capacity. It is likely that these proteins play a key role in generating free amino acids during pollen digestion. The identification of proteins functioning in carbohydrate hydrolysis is consistent with Heliconius butterflies consuming nectar, like other lepidopterans, as well as pollen. Immune-related proteins in saliva are also expected, given that ingestion of pathogens is a likely route to infection. The few “housekeeping” proteins are likely not true salivary proteins and reflect a modest level of contamination that occurred during saliva collection. Among the unannotated proteins were two sets of paralogs, each seemingly the result of a relatively recent tandem duplication. These results offer a first glimpse into the molecular foundation of Heliconius pollen feeding and provide a substantial advance towards comprehensively understanding this striking evolutionary novelty.

Genome-wide identification and characterization of Fox genes in the silkworm, Bombyx mori.

The forkhead box (Fox) transcription factor family has a characteristic of forkhead domain, a winged DNA-binding domain. The Fox genes have been classified into 23 subfamilies, designated FoxA to FoxS, of which the FoxR and FoxS subfamilies are specific to vertebrates. In this review, using whole-genome scanning, we identified 17 distinct Fox genes distributed on 13 chromosomes of the silkworm, Bombyx mori. A phylogenetic tree showed that the silkworm Fox genes could be classified into 13 subfamilies. The FoxK subfamily is specifically absent from the silkworm, although it is present in other lepidopteran insects, including Danaus plexippus and Heliconius melpomene. Microarray data revealed that the Fox genes have distinct expression patterns in the tissues on day 3 of the 5th instar larva. A Gene Ontology analysis suggested that the Fox genes have roles in cellular components, molecular functions, and biological processes, except in pore complex biogenesis. An analysis of the selective pressure on the proteins indicated that most of the amino acid sites in the Fox proteins are undergoing strong purifying selection. Here, we summarize the general characteristics of the Fox genes in the silkworm, which should support further functional studies of the silkworm Fox proteins.

The taxonomic identity of Heliconius melpomene f. pyritosa var. fumigata Zikán (Lepidoptera: Nymphalidae)

The taxonomic identity of the butterfly Heliconius melpomene f. pyritosa var. fumigata Zikán, 1937 (Nymphalidae: Heliconiinae) is discussed based on the discovery of the specimen based on which this name was assigned. The specimen is not a variation of Eueides tales surdus Stichel, 1903, as previously stated, but is in fact a variation of H. melpomene f. pyritosa, which is a synonym of Heliconius erato amalfreda Riffarth, 1901.

The Lepidoptera Odorant Binding Protein gene family: Gene gain and loss within the GOBP/PBP complex of moths and butterflies

Butterflies and moths differ significantly in their daily activities: butterflies are diurnal while moths are largely nocturnal or crepuscular. This life history difference is presumably reflected in their sensory biology, and especially the balance between the use of chemical versus visual signals. Odorant Binding Proteins (OBP) are a class of insect proteins, at least some of which are thought to orchestrate the transfer of odor molecules within an olfactory sensillum (olfactory organ), between the air and odor receptor proteins (ORs) on the olfactory neurons. A Lepidoptera specific subclass of OBPs are the GOBPs and PBPs; these were the first OBPs studied and have well documented associations with olfactory sensilla. We have used the available genomes of two moths, Manduca sexta and Bombyx mori, and two butterflies, Danaus plexippus and Heliconius melpomene, to characterize the GOBP/PBP genes, attempting to identify gene orthologs and document specific gene gain and loss. First, we identified the full repertoire of OBPs in the M. sexta genome, and compared these with the full repertoire of OBPs from the other three lepidopteran genomes, the OBPs of Drosophila melanogaster and select OBPs from other Lepidoptera. We also evaluated the tissue specific expression of the M. sexta OBPs using an available RNAseq databases. In the four lepidopteran species, GOBP2 and all PBPs reside in single gene clusters; in two species GOBP1 is documented to be nearby, about 100 kb from the cluster; all GOBP/PBP genes share a common gene structure indicating a common origin. As such, the GOBP/PBP genes form a gene complex. Our findings suggest that (1) the lepidopteran GOBP/PBP complex is a monophyletic lineage with origins deep within Lepidoptera phylogeny, (2) within this lineage PBP gene evolution is much more dynamic than GOBP gene evolution, and (3) butterflies may have lost a PBP gene that plays an important role in moth pheromone detection, correlating with a shift from olfactory (moth) to visual (butterfly) communication, at least regarding long distance mate recognition. These findings will be clarified by additional lepidopteran genomic data, but the observation that moths and butterflies share most of the PBP/GOBP genes suggests that they also share common chemosensory-based behavioral pathways.

Complete mitochondrial genome sequence of Heliconius melpomene rosina (Insecta: Lepidoptera: Nymphalidae)

The Postman Butterfly (Heliconius melpomene) is one of the heliconiine butterflies found in the Central and South America. In this study, the complete mitochondrial genome sequence of its subspecies, Heliconius melpomene rosina, is determined for the first time. Results show that this circular genome is 15,327 bp in length, and consists of 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs) and a putative control region (D-loop). The genome organization, nucleotide composition and codon usage are identical to those reported from other butterfly mitochondrial genomes. The whole nucleotide composition is 39.33% of A, 42.33% of T, 10.88% of C and 7.46% of G, with a relatively lower level of G and an extremely higher AT content of 81.66%. Most of the PCGs initiate with the ATN start codons except for ND1 and COX1 genes, which separately start with TTG and CGA. Totally 13 PCGs terminate with the standard canonical stop codons of TAA/TAG/TGA or an incomplete stop codon of T– –. The mitochondrial genome sequence reported here would be potentially useful for the phylogenetic studies in butterflies.

Estimation of the spontaneous mutation rate in Heliconius melpomene.

We estimated the spontaneous mutation rate in Heliconius melpomene by genome sequencing of a pair of parents and 30 of their offspring, based on the ratio of number of de novo heterozygotes to the number of callable site-individuals. We detected nine new mutations, each one affecting a single site in a single offspring. This yields an estimated mutation rate of 2.9 × 10−9 (95% confidence interval, 1.3 × 10−9–5.5 × 10−9), which is similar to recent estimates in Drosophila melanogaster, the only other insect species in which the mutation rate has been directly estimated. We infer that recent effective population size of H. melpomene is about 2 million, a substantially lower value than its census size, suggesting a role for natural selection reducing diversity. We estimate that H. melpomene diverged from its Müllerian comimic H. erato about 6 Ma, a somewhat later date than estimates based on a local molecular clock.

Phylogeography of Heliconius cydno and its closest relatives: disentangling their origin and diversification.

The origins of the extraordinary diversity within the Neotropics have long fascinated biologists and naturalists. Yet, the underlying factors that have given rise to this diversity remain controversial. To test the relative importance of Quaternary climatic change and Neogene tectonic and paleogeographic reorganizations in the generation of biodiversity, we examine intraspecific variation across the Heliconius cydno radiation and compare this variation to that within the closely related Heliconius melpomene and Heliconius timareta radiations. Our data, which consist of both mtDNA and genome-scan data from nearly 2250 amplified fragment length polymorphism (AFLP) loci, reveal a complex history of differentiation and admixture at different geographic scales. Both mtDNA and AFLP phylogenies suggest that H.timareta and H.cydno are probably geographic extremes of the same radiation that probably diverged from H.melpomene prior to the Pliocene-Pleistocene boundary, consistent with hypotheses of diversification that rely on geological events in the Pliocene. The mtDNA suggests that this radiation originated in Central America or the northwestern region of South America, with a subsequent colonization of the eastern and western slopes of the Andes. Our genome-scan data indicate significant admixture among sympatric H.cydno/H.timareta and H.melpomene populations across the extensive geographic ranges of the two radiations. Within H.cydno, both mtDNA and AFLP data indicate significant population structure at local scales, with strong genetic differences even among adjacent H.cydno colour pattern races. These genetic patterns highlight the importance of past geoclimatic events, intraspecific gene flow, and local population differentiation in the origin and establishment of new adaptive forms.

Molecular evolution of the odorant and gustatory receptor genes in lepidopteran insects: implications for their adaptation and speciation.

Lepidoptera (comprised of butterflies and moths) is one of the largest groups of insects, including more than 160,000 described species. Chemoreception plays important roles in the adaptation of these species to a wide range of niches, e.g., plant hosts, egg-laying sites, and mates. This study investigated the molecular evolution of the lepidopteran odorant (Or) and gustatory receptor (Gr) genes using recently identified genes from Bombyx mori, Danaus plexippus, Heliconius melpomene, Plutella xylostella, Heliothis virescens, Manduca sexta, Cydia pomonella, and Spodoptera littoralis. A limited number of cases of large lineage-specific gene expansion are observed (except in the P. xylostella lineage), possibly due to selection against tandem gene duplication. There has been strong purifying selection during the evolution of both lepidopteran odorant and gustatory genes, as shown by the low ω values estimated through CodeML analysis, ranging from 0.0093 to 0.3926. However, purifying selection has been relaxed on some amino acid sites in these receptors, leading to sequence divergence, which is a precursor of positive selection on these sequences. Signatures of positive selection were detected only in a few loci from the lineage-specific analysis. Estimation of gene gains and losses suggests that the common ancestor of the Lepidoptera had fewer Or genes compared to extant species and an even more reduced number of Gr genes, particularly within the bitter receptor clade. Multiple gene gains and a few gene losses occurred during the evolution of Lepidoptera. Gene family expansion may be associated with the adaptation of lepidopteran species to plant hosts, especially after angiosperm radiation. Phylogenetic analysis of the moth sex pheromone receptor genes suggested that chromosomal translocations have occurred several times. New sex pheromone receptors have arisen through tandem gene duplication. Positive selection was detected at some amino acid sites predicted to be in the extracellular and transmembrane regions of the newly duplicated genes, which might be associated with the evolution of the new pheromone receptors.

Neighboring genes shaping a single adaptive mimetic trait

The colorful wing patterns of Heliconius butterflies represent an excellent system in which to study the genetic and developmental control of adaptation and convergence. Using qRT-PCR and in situ hybridization on developing wings of the co-mimic species Heliconius melpomene and Heliconius erato, we have profiled the expression of three candidate genes located in the genomic locus controlling red color pattern variation. We found convergent domains of gene expression in H. melpomene and H. erato associated with red wing elements in the two genes optix and kinesin. During early pupal development of both species, the expression of optix perfectly associated with all red pattern elements whereas that of kinesin was specifically correlated with the presence of the red forewing band. These results provide evidence for the use of these two tightly linked patterning genes, acting together to create convergent wing phenotypes in Heliconius and constituting a hotspot of adaptation.

Multiple recent co-options of Optix associated with novel traits in adaptive butterfly wing radiations

BackgroundWhile the ecological factors that drive phenotypic radiations are often well understood, less is known about the generative mechanisms that cause the emergence and subsequent diversification of novel features. Heliconius butterflies display an extraordinary diversity of wing patterns due in part to mimicry and sexual selection. Identifying the genetic drivers of this crucible of evolution is now within reach, as it was recently shown that cis-regulatory variation of the optix transcription factor explains red pattern differences in the adaptive radiations of the Heliconius melpomene and Heliconius erato species groups.ResultsHere, we compare the developmental expression of the Optix protein across a large phylogenetic sample of butterflies and infer that its color patterning role originated at the base of the neotropical passion-vine butterfly clade (Lepidoptera, Nymphalidae, Tribe: Heliconiini), shortly predating multiple Optix-driven wing pattern radiations in the speciose Heliconius and Eueides genera. We also characterize novel Optix and Doublesex expression in the male-specific pheromone wing scales of the basal heliconiines Dryas and Agraulis, thus illustrating that within the Heliconinii lineage, Optix has been evolutionarily redeployed in multiple contexts in association with diverse wing features.ConclusionsOur findings reveal that the repeated co-option of Optix into various aspects of wing scale specification was associated with multiple evolutionary novelties over a relatively short evolutionary time scale. In particular, the recruitment of Optix expression in colored scale cell precursors was a necessary condition to the explosive diversification of passion-vine butterfly wing patterns. The novel deployment of a gene followed by spatial modulation of its expression in a given cell type could be a common mode of developmental innovation for triggering phenotypic radiations.

Transposable element evolution in Heliconius suggests genome diversity within Lepidoptera.

BackgroundTransposable elements (TEs) have the potential to impact genome structure, function and evolution in profound ways. In order to understand the contribution of transposable elements (TEs) to Heliconius melpomene, we queried the H. melpomene draft sequence to identify repetitive sequences.ResultsWe determined that TEs comprise ~25% of the genome. The predominant class of TEs (~12% of the genome) was the non-long terminal repeat (non-LTR) retrotransposons, including a novel SINE family. However, this was only slightly higher than content derived from DNA transposons, which are diverse, with several families having mobilized in the recent past. Compared to the only other well-studied lepidopteran genome, Bombyx mori, H. melpomene exhibits a higher DNA transposon content and a distinct repertoire of retrotransposons. We also found that H. melpomene exhibits a high rate of TE turnover with few older elements accumulating in the genome.ConclusionsOur analysis represents the first complete, de novo characterization of TE content in a butterfly genome and suggests that, while TEs are able to invade and multiply, TEs have an overall deleterious effect and/or that maintaining a small genome is advantageous. Our results also hint that analysis of additional lepidopteran genomes will reveal substantial TE diversity within the group.

Female Behaviour Drives Expression and Evolution of Gustatory Receptors in Butterflies

Secondary plant compounds are strong deterrents of insect oviposition and feeding, but may also be attractants for specialist herbivores. These insect-plant interactions are mediated by insect gustatory receptors (Grs) and olfactory receptors (Ors). An analysis of the reference genome of the butterfly Heliconius melpomene, which feeds on passion-flower vines (Passiflora spp.), together with whole-genome sequencing within the species and across the Heliconius phylogeny has permitted an unprecedented opportunity to study the patterns of gene duplication and copy-number variation (CNV) among these key sensory genes. We report in silico gene predictions of 73 Gr genes in the H. melpomene reference genome, including putative CO2, sugar, sugar alcohol, fructose, and bitter receptors. The majority of these Grs are the result of gene duplications since Heliconius shared a common ancestor with the monarch butterfly or the silkmoth. Among Grs but not Ors, CNVs are more common within species in those gene lineages that have also duplicated over this evolutionary time-scale, suggesting ongoing rapid gene family evolution. Deep sequencing (∼1 billion reads) of transcriptomes from proboscis and labial palps, antennae, and legs of adult H. melpomene males and females indicates that 67 of the predicted 73 Gr genes and 67 of the 70 predicted Or genes are expressed in these three tissues. Intriguingly, we find that one-third of all Grs show female-biased gene expression (n = 26) and nearly all of these (n = 21) are Heliconius-specific Grs. In fact, a significant excess of Grs that are expressed in female legs but not male legs are the result of recent gene duplication. This difference in Gr gene expression diversity between the sexes is accompanied by a striking sexual dimorphism in the abundance of gustatory sensilla on the forelegs of H. melpomene, suggesting that female oviposition behaviour drives the evolution of new gustatory receptors in butterfly genomes.

Cytochrome P450‐encoding genes from the Heliconius genome as candidates for cyanogenesis

Cytochrome P450s are important both in the metabolism of xenobiotics and the production of compounds such as cyanogenic glucosides, which insects use in their defence. In the present study, we use transcriptomic and genomic information to isolate and name P450-encoding genes from the butterfly Heliconius melpomene. We classify each of the putative genes into its appropriate superfamily and compare the distribution of P450s across sequenced insects. We also identify homologues of two P450s known to be involved in cyanogenesis in the six-spot Burnet moth, Zygaena filipendulae. Classification of Heliconius P450s should be an important step in the dissection of their role in the exploitation of their host plant, the passion vine Passiflora.

Genomic architecture of adaptive color pattern divergence and convergence in Heliconius butterflies

Identifying the genetic changes driving adaptive variation in natural populations is key to understanding the origins of biodiversity. The mosaic of mimetic wing patterns in Heliconius butterflies makes an excellent system for exploring adaptive variation using next-generation sequencing. In this study, we use a combination of techniques to annotate the genomic interval modulating red color pattern variation, identify a narrow region responsible for adaptive divergence and convergence in Heliconius wing color patterns, and explore the evolutionary history of these adaptive alleles. We use whole genome resequencing from four hybrid zones between divergent color pattern races of Heliconius erato and two hybrid zones of the co-mimic Heliconius melpomene to examine genetic variation across 2.2 Mb of a partial reference sequence. In the intergenic region near optix, the gene previously shown to be responsible for the complex red pattern variation in Heliconius, population genetic analyses identify a shared 65-kb region of divergence that includes several sites perfectly associated with phenotype within each species. This region likely contains multiple cis-regulatory elements that control discrete expression domains of optix. The parallel signatures of genetic differentiation in H. erato and H. melpomene support a shared genetic architecture between the two distantly related co-mimics; however, phylogenetic analysis suggests mimetic patterns in each species evolved independently. Using a combination of next-generation sequencing analyses, we have refined our understanding of the genetic architecture of wing pattern variation in Heliconius and gained important insights into the evolution of novel adaptive phenotypes in natural populations.

Introgression of wing pattern alleles and speciation via homoploid hybridization inHeliconiusbutterflies: a review of evidence from the genome

The diverse Müllerian mimetic wing patterns of neotropical Heliconius (Nymphalidae) have been proposed to be not only aposematic signals to potential predators, but also intra- and interspecific recognition signals that allow the butterflies to maintain their specific identities, and which perhaps drive the process of speciation, as well. Adaptive features under differential selection that also serve as cues for assortative mating have been referred to as 'magic traits', which can drive ecological speciation. Such traits are expected to exhibit allelic differentiation between closely related species with ongoing gene flow, whereas unlinked neutral traits are expected to be homogenized to a greater degree by introgression. However, recent evidence suggests that interspecific hybridization among Heliconius butterflies may have resulted in adaptive introgression of these very same traits across species boundaries, and in the evolution of new species by homoploid hybrid speciation. The theory and data supporting various aspects of the apparent paradox of 'magic trait' introgression are reviewed, with emphasis on population genomic comparisons of Heliconius melpomene and its close relatives.

Evolutionary Ecological Genomics

The recognition that evolution can happen on ecological timescales (Hairston et al. 2005; Pelletier et al. 2009; Ellner et al. 2011; Becks et al. 2012) has prompted the integration of ecology and evolution, while easier access to high-throughput sequencing technologies has increased the number of genetic nonmodel species entering the ‘omics’ era (e.g. Turner et al. 2010; Colbourne et al. 2011; Jones et al. 2012). We are now in a position to identify the genetic basis of adaptation and the mechanisms of adaptive responses in the wild. It nonetheless remains a challenge to go beyond descriptive measures of patterns of genetic variation and to identify the evolutionary processes driving species adaptation and evolution. This special issue represents a broad cross-section of research into evolutionary adaptation at the genetic level. The approaches used vary from classic QTL studies to RAD sequencing and RNAseq. Given the rapid advance of sequencing technology, we fully expect that ‘genomic’ as defined here will be merely ‘genetic’ in a few years, but we nonetheless hope that the results and methods described in this special issue will serve as a blueprint for future work in this field.

Multiple ancient horizontal gene transfers and duplications in lepidopteran species

Eukaryotic horizontal gene transfer (HGT) events are increasingly being discovered yet few reports have summarized multiple occurrences in a wide range of species. We systematically investigated HGT events in the order Lepidoptera by employing a series of filters. Bombyx mori, Danaus plexippus and Heliconius melpomene had 13, 12 and 12 HGTs, respectively, from bacteria and fungi. These HGTs contributed a total of 64 predicted genes: 22 to B. mori, 22 to D. plexippus and 20 to H. melpomene. Several new genes were generated by post-transfer duplications. Post-transfer duplication of a suite of functional HGTs has rarely been reported in higher organisms. The distributional patterns of paralogues for certain genes differed in the three species, indicating potential independent duplication or loss events. All of these HGTs had homologues expressed in some other lepidopterans, indicating ancient transfer events. Most HGTs were involved in the metabolism of sugar and amino acids. These HGTs appeared to have experienced amelioration, purifying selection and accelerated evolution to adapt to the background genome of the recipient. The discovery of ancient, massive HGTs and duplications in lepidopterans and their adaptive evolution provides further insights into the evolutionary significance of the events from donors to multicellular host recipients.