Contrasting genomic and phenotypic outcomes of hybridization between pairs of mimetic butterfly taxa across a suture zone

Spatiotemporal dynamics in butterfly hybrid zones

Evaluating whether hybrid zones are stable or mobile can provide novel insights for evolution and conservation biology. Butterflies exhibit high sensitivity to environmental changes and represent an important model system for the study of hybrid zone origins and maintenance. Here, we review the literature exploring butterfly hybrid zones, with a special focus on their spatiotemporal dynamics and the potential mechanisms that could lead to their movement or stability. We then compare different lines of evidence used to investigate hybrid zone dynamics and discuss the strengths and weaknesses of each approach. Our goal with this review is to reveal general conditions associated with the stability or mobility of butterfly hybrid zones by synthesizing evidence obtained using different types of data sampled across multiple regions and spatial scales. Finally, we discuss spatiotemporal dynamics in the context of a speciation/divergence continuum, the relevance of hybrid zones for conservation biology, and recommend key topics for future investigation.

Quantitative Analyses of Coupling in Hybrid Zones

In hybrid zones, whether barrier loci experience selection mostly independently or as a unit depends on the ratio of selection to recombination as captured by the coupling coefficient. Theory predicts a sharper transition between an uncoupled and coupled system when more loci affect hybrid fitness. However, the extent of coupling in hybrid zones has rarely been quantified. Here, we use simulations to characterize the relationship between the coupling coefficient and variance in clines across genetic loci. We then reanalyze 25 hybrid zone data sets and find that cline variances and estimated coupling coefficients form a smooth continuum from high variance and weak coupling to low variance and strong coupling. Our results are consistent with low rates of hybridization and a strong genome-wide barrier to gene flow when the coupling coefficient is much greater than 1, but also suggest that this boundary might be approached gradually and at a near constant rate over time.

First chromosome scale genomes of ithomiine butterflies (Nymphalidae: Ithomiini): Comparative models for mimicry genetic studies

The ithomiine butterflies (Nymphalidae: Danainae) represent the largest known radiation of Müllerian mimetic butterflies. They dominate by number the mimetic butterfly communities, which include species such as the iconic neotropical Heliconius genus. Recent studies on the ecology and genetics of speciation in Ithomiini have suggested that sexual pheromones, colour pattern and perhaps hostplant could drive reproductive isolation. However, no reference genome was available for Ithomiini, which has hindered further exploration on the genetic architecture of these candidate traits, and more generally on the genomic patterns of divergence. Here, we generated high-quality, chromosome-scale genome assemblies for two Melinaea species, M. marsaeus and M. menophilus, and a draft genome of the species Ithomia salapia. We obtained genomes with a size ranging from 396 to 503 Mb across the three species and scaffold N50 of 40.5 and 23.2 Mb for the two chromosome-scale assemblies. Using collinearity analyses we identified massive rearrangements between the two closely related Melinaea species. An annotation of transposable elements and gene content was performed, as well as a specialist annotation to target chemosensory genes, which is crucial for host plant detection and mate recognition in mimetic species. A comparative genomic approach revealed independent gene expansions in ithomiines and particularly in gustatory receptor genes. These first three genomes of ithomiine mimetic butterflies constitute a valuable addition and a welcome comparison to existing biological models such as Heliconius, and will enable further understanding of the mechanisms of adaptation in butterflies.

Redundancy analysis, genome-wide association studies and the pigmentation of brown trout (Salmo trutta L.)

The association of molecular variants with phenotypic variation is a main issue in biology, often tackled with genome-wide association studies (GWAS). GWAS are challenging, with increasing, but still limited, use in evolutionary biology. We used redundancy analysis (RDA) as a complimentary ordination approach to single- and multitrait GWAS to explore the molecular basis of pigmentation variation in brown trout (Salmo trutta) belonging to wild populations impacted by hatchery fish. Based on 75,684 single nucleotide polymorphic (SNP) markers, RDA, single- and multitrait GWAS allowed the extraction of 337 independent colour patterning loci (CPLs) associated with trout pigmentation traits, such as the number of red and black spots on flanks. Collectively, these CPLs (i) mapped onto 35 out of 40 brown trout linkage groups indicating a polygenic genomic architecture of pigmentation, (ii) were found to be associated with 218 candidate genes, including 197 genes formerly mentioned in the literature associated to skin pigmentation, skin patterning, differentiation or structure notably in a close relative, the rainbow trout (Onchorhynchus mykiss), and (iii) related to functions relevant to pigmentation variation (e.g., calcium- and ion-binding, cell adhesion). Annotated CPLs include genes with well-known pigmentation effects (e.g., PMEL, SLC45A2, SOX10), but also markers associated with genes formerly found expressed in rainbow or brown trout skins. RDA was also shown to be useful to investigate management issues, especially the dynamics of trout pigmentation submitted to several generations of hatchery introgression.

ClineHelpR: an R package for genomic cline outlier detection and visualization

Background

Patterns of multi-locus differentiation (i.e., genomic clines) often extend broadly across hybrid zones and their quantification can help diagnose how species boundaries are shaped by adaptive processes, both intrinsic and extrinsic. In this sense, the transitioning of loci across admixed individuals can be contrasted as a function of the genome-wide trend, in turn allowing an expansion of clinal theory across a much wider array of biodiversity. However, computational tools that serve to interpret and consequently visualize ‘genomic clines’ are limited, and users must often write custom, relatively complex code to do so.

Results

Here, we introduce the ClineHelpR R-package for visualizing genomic clines and detecting outlier loci using output generated by two popular software packages, bgc and Introgress. ClineHelpR bundles both input generation (i.e., filtering datasets and creating specialized file formats) and output processing (e.g., MCMC thinning and burn-in) with functions that directly facilitate interpretation and hypothesis testing. Tools are also provided for post-hoc analyses that interface with external packages such as ENMeval and RIdeogram.

Conclusions

Our package increases the reproducibility and accessibility of genomic cline methods, thus allowing an expanded user base and promoting these methods as mechanisms to address diverse evolutionary questions in both model and non-model organisms. Furthermore, the ClineHelpR extended functionality can evaluate genomic clines in the context of spatial and environmental features, allowing users to explore underlying processes potentially contributing to the observed patterns and helping facilitate effective conservation management strategies.

Comparative transcriptome analysis at the onset of speciation in a mimetic butterfly-The Ithomiini Melinaea marsaeus

Ecological speciation entails divergent selection on specific traits and ultimately on the developmental pathways responsible for these traits. Selection can act on gene sequences but also on regulatory regions responsible for gene expression. Mimetic butterflies are a relevant system for speciation studies because wing colour pattern (WCP) often diverges between closely related taxa and is thought to drive speciation through assortative mating and increased predation on hybrids. Here, we generate the first transcriptomic resources for a mimetic butterfly of the tribe Ithomiini, Melinaea marsaeus, to examine patterns of differential expression between two subspecies and between tissues that express traits that likely drive reproductive isolation; WCP and chemosensory genes. We sequenced whole transcriptomes of three life stages to cover a large catalogue of transcripts, and we investigated differential expression between subspecies in pupal wing discs and antennae. Eighteen known WCP genes were expressed in wing discs and 115 chemosensory genes were expressed in antennae, with a remarkable diversity of chemosensory protein genes. Many transcripts were differentially expressed between subspecies, including two WCP genes and one odorant receptor. Our results suggest that in M. marsaeus the same genes as in other mimetic butterflies are involved in traits causing reproductive isolation, and point at possible candidates for the differences in those traits between subspecies. Differential expression analyses of other developmental stages and body organs and functional studies are needed to confirm and expand these results. Our work provides key resources for comparative genomics in mimetic butterflies, and more generally in Lepidoptera.

3-Acetoxy-fatty acid isoprenyl esters from androconia of the ithomiine butterfly Ithomia salapia

Male ithomiine butterflies (Nymphalidae: Danainae) have hairpencils on the forewings (i.e., androconia) that disseminate semiochemicals during courtship. While most ithomiines are known to contain derivatives of pyrrolizidine alkaloids, dihydropyrrolizines, or γ-lactones in these androconia, here we report on a new class of fatty acid esters identified in two subspecies, Ithomia salapia aquinia and I. s. derasa. The major components were identified as isoprenyl (3-methyl-3-butenyl) (Z)-3-acetoxy-11-octadecenoate, isoprenyl (Z)-3-acetoxy-13-octadecenoate (12) and isoprenyl 3-acetoxyoctadecanoate (11) by GC/MS and GC/IR analyses, microderivatizations, and synthesis of representative compounds. The absolute configuration of 12 was determined to be R. The two subspecies differed not only in the composition of the ester bouquet, but also in the composition of more volatile androconial constituents. While some individuals of I. s. aquinia contained ithomiolide A (3), a pyrrolizidine alkaloid derived γ-lactone, I. s. derasa carried the sesquiterpene α-elemol (8) in the androconia. These differences might be important for the reproductive isolation of the two subspecies, in line with previously reported low gene exchange between the two species in regions where they co-occur. Furthermore, the occurrence of positional isomers of unsaturated fatty acid derivatives indicates activity of two different desaturases within these butterflies, Δ9 and Δ11, which has not been reported before in male Lepidoptera.