Subtle Introgression Footprints at the End of the Speciation Continuum in a Clade of Heliconius Butterflies.

Quantifying gene flow between lineages at different stages of the speciation continuum is central to understanding speciation. Heliconius butterflies have undergone an adaptive radiation in wing color patterns driven partly by natural selection for local mimicry. Color patterns are also known to be used as assortative mating cues. Therefore, wing pattern divergence is considered to play a role in speciation. A corollary is that mimicry between closely related species may be associated with hybridization and interfere with reproductive isolation. Here, we take a multifaceted approach to explore speciation history, species boundaries, and traits involved in species differentiation between the two closely related species, Heliconius hecale and Heliconius ismenius. We focus on geographic regions where the two species mimic each other and contrast this with geographic regions where they do not mimic each other. To examine population history and patterns of gene flow, we tested and compared a four-population model accounting for linked selection. This model suggests that the two species have remained isolated for a large part of their history, yet with a small amount of gene exchange. Accordingly, signatures of genomic introgression were small except at a major wing pattern allele and chemosensing genes and stronger in the mimetic populations compared with nonmimetic populations. Behavioral assays confirm that visual confusion exists but that short-range cues determine strong sexual isolation. Tests for chemical differentiation between species identified major differences in putative pheromones which likely mediate mate choice and the maintenance of species differences.

A high-density linkage map reveals broad- and fine-scale sex differences in recombination in the hihi (stitchbird; Notiomystis cincta).

Recombination, the process of DNA exchange between homologous chromosomes during meiosis, plays a major role in genomic diversity and evolutionary change. Variation in recombination rate is widespread despite recombination often being essential for progression of meiosis. One such variation is heterochiasmy, where recombination rates differ between sexes. Heterochiasmy has been observed across broad taxonomic groups, yet it remains an evolutionary enigma. We used Lep-MAP3, a pedigree-based software that is efficient in handling large datasets, to generate linkage maps for the hihi or stitchbird (Notiomystis cincta), utilising information from >36 K SNPs and 36 families. We constructed 29 linkage maps, including for the previously unscaffolded Z chromosome. The hihi is an endangered passerine endemic to Aotearoa New Zealand that is sexually dimorphic and exhibits high levels of sexual conflict, including sperm competition. Patterns in recombination in the hihi are consistent with those in other birds, including higher recombination rates in micro-chromosomes. Heterochiasmy in the hihi is male-biased, in line with predictions of the Haldane-Huxley rule, with the male linkage map being 15% longer. Micro-chromosomes exhibit heterochiasmy to a greater extent, contrary to that reported in other birds. At the intra-chromosomal level, heterochiasmy is higher nearer to chromosome ends and in gene-rich regions. Regions of extreme heterochiasmy are enriched for genes implicated in cell structure. This study adds an important contribution in assessing evolutionary theories of heterochiasmy and provides a framework for future studies investigating fine-scale heterochiasmy.

A cry for kelp: Evidence for polyphenolic inhibition of Oxford Nanopore sequencing of brown algae.

Genomic resources have yielded unprecedented insights into ecological and evolutionary processes, not to mention their importance in economic and conservation management of specific organisms. However, the field of macroalgal genomics is hampered by difficulties in the isolation of suitable DNA. Even when DNA that appears high quality by standard metrics has been isolated, such samples may not perform well during the sequencing process. We here have compared Oxford Nanopore long-read sequencing results for three species of macroalgae to those of nonmacroalgal species and determined that when using macroalgal samples, sequencing activity declined rapidly, resulting in reduced sequencing yield. Chemical analysis of macroalgal DNA that would be considered suitable for sequencing revealed that DNA derived from dried macroalgae was enriched for polyphenol-DNA adducts (DNA with large polyphenols chemically attached to it), which may have led to sequencing inhibition. Of note, we observed the strongest evidence of sequencing inhibition and reduced sequence output when using samples dried using silica gel-suggesting that such storage approaches may not be appropriate for samples destined for Oxford Nanopore sequencing. Our findings have wide-ranging implications for the generation of genomic resources from macroalgae and suggest a need to develop new storage methods that are more amenable to Oxford Nanopore sequencing or to use fresh flash-frozen tissue wherever possible for genome sequencing.

MtING2 encodes an ING domain PHD finger protein which affects Medicago growth, flowering, global patterns of H3K4me3, and gene expression.

Flowering of the reference legume Medicago truncatula is promoted by winter cold (vernalization) followed by long-day photoperiods (VLD) similar to winter annual Arabidopsis. However, Medicago lacks FLC and CO, key regulators of Arabidopsis VLD flowering. Most plants have two INHIBITOR OF GROWTH (ING) genes (ING1 and ING2), encoding proteins with an ING domain with two anti-parallel alpha-helices and a plant homeodomain (PHD) finger, but their genetic role has not been previously described. In Medicago, Mting1 gene-edited mutants developed and flowered normally, but an Mting2-1 Tnt1 insertion mutant and gene-edited Mting2 mutants had developmental abnormalities including delayed flowering particularly in VLD, compact architecture, abnormal leaves with extra leaflets but no trichomes, and smaller seeds and barrels. Mting2 mutants had reduced expression of activators of flowering, including the FT-like gene MtFTa1, and increased expression of the candidate repressor MtTFL1c, consistent with the delayed flowering of the mutant. MtING2 overexpression complemented Mting2-1, but did not accelerate flowering in wild type. The MtING2 PHD finger bound H3K4me2/3 peptides weakly in vitro, but analysis of gene-edited mutants indicated that it was dispensable to MtING2 function in wild-type plants. RNA sequencing experiments indicated that >7000 genes are mis-expressed in the Mting2-1 mutant, consistent with its strong mutant phenotypes. Interestingly, ChIP-seq analysis identified >5000 novel H3K4me3 locations in the genome of Mting2-1 mutants compared to wild type R108. Overall, our mutant study has uncovered an important physiological role of a plant ING2 gene in development, flowering, and gene expression, which likely involves an epigenetic mechanism.

Genomic signatures of inbreeding depression for a threatened Aotearoa New Zealand passerine.

For small and isolated populations, the increased chance of mating between related individuals can result in a substantial reduction in individual and population fitness. Despite the increasing availability of genomic data to measure inbreeding accurately across the genome, inbreeding depression studies for threatened species are still scarce due to the difficulty of measuring fitness in the wild. Here, we investigate inbreeding and inbreeding depression for the extensively monitored Tiritiri Matangi island population of a threatened Aotearoa New Zealand passerine, the hihi (Notiomystis cincta). First, using a custom 45 k single nucleotide polymorphism (SNP) array, we explore genomic inbreeding patterns by inferring homozygous segments across the genome. Although all individuals have similar levels of ancient inbreeding, highly inbred individuals are affected by recent inbreeding, which can probably be explained by bottleneck effects such as habitat loss after European arrival and their translocation to the island in the 1990s. Second, we investigate genomic inbreeding effects on fitness, measured as lifetime reproductive success, and its three components, juvenile survival, adult annual survival and annual reproductive success, in 363 hihi. We find that global inbreeding significantly affects juvenile survival but none of the remaining fitness traits. Finally, we employ a genome-wide association approach to test the locus-specific effects of inbreeding on fitness, and identify 13 SNPs significantly associated with lifetime reproductive success. Our findings suggest that inbreeding depression does impact hihi, but at different genomic scales for different traits, and that purging has therefore failed to remove all variants with deleterious effects from this population of conservation concern.

Tracing the introduction of the invasive common myna using population genomics.

The common myna (Acridotheres tristis) is one of the most invasive bird species in the world, yet its colonisation history is only partly understood. We identified the introduction history and population structure, and quantified the genetic diversity of myna populations from the native range in India and introduced populations in New Zealand, Australia, Fiji, Hawaii, and South Africa, based on thousands of single nucleotide polymorphism markers in 814 individuals. We were able to identify the source population of mynas in several invasive locations: mynas from Fiji and Melbourne, Australia, were likely founded by individuals from a subpopulation in Maharashtra, India, while mynas in Hawaii and South Africa were likely independently founded by individuals from other localities in India. Our findings suggest that New Zealand mynas were founded by individuals from Melbourne, which, in turn, were founded by individuals from Maharashtra. We identified two genetic clusters among New Zealand mynas, divided by New Zealand's North Island's axial mountain ranges, confirming previous observations that mountains and thick forests may form barriers to myna dispersal. Our study provides a foundation for other population and invasion genomic studies and provides useful information for the management of this invasive species.

The gene cortex controls mimicry and crypsis in butterflies and moths.

The wing patterns of butterflies and moths (Lepidoptera) are diverse and striking examples of evolutionary diversification by natural selection. Lepidopteran wing colour patterns are a key innovation, consisting of arrays of coloured scales. We still lack a general understanding of how these patterns are controlled and whether this control shows any commonality across the 160,000 moth and 17,000 butterfly species. Here, we use fine-scale mapping with population genomics and gene expression analyses to identify a gene, cortex, that regulates pattern switches in multiple species across the mimetic radiation in Heliconius butterflies. cortex belongs to a fast-evolving subfamily of the otherwise highly conserved fizzy family of cell-cycle regulators, suggesting that it probably regulates pigmentation patterning by regulating scale cell development. In parallel with findings in the peppered moth (Biston betularia), our results suggest that this mechanism is common within Lepidoptera and that cortex has become a major target for natural selection acting on colour and pattern variation in this group of insects.

Genomic data of different resolutions reveal consistent inbreeding estimates but contrasting homozygosity landscapes for the threatened Aotearoa New Zealand hihi.

Inbreeding can lead to a loss of heterozygosity in a population and when combined with genetic drift may reduce the adaptive potential of a species. However, there is uncertainty about whether resequencing data can provide accurate and consistent inbreeding estimates. Here, we performed an in-depth inbreeding analysis for hihi (Notiomystis cincta), an endemic and nationally vulnerable passerine bird of Aotearoa New Zealand. We first focused on subsampling variants from a reference genome male, and found that low-density data sets tend to miss runs of homozygosity (ROH) in some places and overestimate ROH length in others, resulting in contrasting homozygosity landscapes. Low-coverage resequencing and 50 K SNP array densities can yield comparable inbreeding results to high-coverage resequencing approaches, but the results for all data sets are highly dependent on the software settings employed. Second, we extended our analysis to 10 hihi where low-coverage whole genome resequencing, RAD-seq and SNP array genotypes are available. We inferred ROH and individual inbreeding to evaluate the relative effects of sequencing depth versus SNP density on estimating inbreeding coefficients and found that high rates of missingness downwardly bias both the number and length of ROH. In summary, when using genomic data to evaluate inbreeding, studies must consider that ROH estimates are heavily dependent on analysis parameters, data set density and individual sequencing depth.

Selection and gene flow shape genomic islands that control floral guides.

Genomes of closely-related species or populations often display localized regions of enhanced relative sequence divergence, termed genomic islands. It has been proposed that these islands arise through selective sweeps and/or barriers to gene flow. Here, we genetically dissect a genomic island that controls flower color pattern differences between two subspecies of Antirrhinum majus, A.m.striatum and A.m.pseudomajus, and relate it to clinal variation across a natural hybrid zone. We show that selective sweeps likely raised relative divergence at two tightly-linked MYB-like transcription factors, leading to distinct flower patterns in the two subspecies. The two patterns provide alternate floral guides and create a strong barrier to gene flow where populations come into contact. This barrier affects the selected flower color genes and tightly-linked loci, but does not extend outside of this domain, allowing gene flow to lower relative divergence for the rest of the chromosome. Thus, both selective sweeps and barriers to gene flow play a role in shaping genomic islands: sweeps cause elevation in relative divergence, while heterogeneous gene flow flattens the surrounding "sea," making the island of divergence stand out. By showing how selective sweeps establish alternative adaptive phenotypes that lead to barriers to gene flow, our study sheds light on possible mechanisms leading to reproductive isolation and speciation.

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

Hybrid zones, whereby divergent lineages come into contact and eventually hybridize, can provide insights on the mechanisms involved in population differentiation and reproductive isolation, and ultimately speciation. Suture zones offer the opportunity to compare these processes across multiple species. In this paper we use reduced-complexity genomic data to compare the genetic and phenotypic structure and hybridization patterns of two mimetic butterfly species, Ithomia salapia and Oleria onega (Nymphalidae: Ithomiini), each consisting of a pair of lineages differentiated for their wing colour pattern and that come into contact in the Andean foothills of Peru. Despite similarities in their life history, we highlight major differences, both at the genomic and phenotypic level, between the two species. These differences include the presence of hybrids, variations in wing phenotype, and genomic patterns of introgression and differentiation. In I. salapia, the two lineages appear to hybridize only rarely, whereas in O. onega the hybrids are not only more common, but also genetically and phenotypically more variable. We also detected loci statistically associated with wing colour pattern variation, but in both species these loci were not over-represented among the candidate barrier loci, suggesting that traits other than wing colour pattern may be important for reproductive isolation. Our results contrast with the genomic patterns observed between hybridizing lineages in the mimetic Heliconius butterflies, and call for a broader investigation into the genomics of speciation in Ithomiini - the largest radiation of mimetic butterflies.

Genetic architecture and balancing selection: the life and death of differentiated variants.

Balancing selection describes any form of natural selection, which results in the persistence of multiple variants of a trait at intermediate frequencies within populations. By offering up a snapshot of multiple co-occurring functional variants and their interactions, systems under balancing selection can reveal the evolutionary mechanisms favouring the emergence and persistence of adaptive variation in natural populations. We here focus on the mechanisms by which several functional variants for a given trait can arise, a process typically requiring multiple epistatic mutations. We highlight how balancing selection can favour specific features in the genetic architecture and review the evolutionary and molecular mechanisms shaping this architecture. First, balancing selection affects the number of loci underlying differentiated traits and their respective effects. Control by one or few loci favours the persistence of differentiated functional variants by limiting intergenic recombination, or its impact, and may sometimes lead to the evolution of supergenes. Chromosomal rearrangements, particularly inversions, preventing adaptive combinations from being dissociated are increasingly being noted as features of such systems. Similarly, due to the frequency of heterozygotes maintained by balancing selection, dominance may be a key property of adaptive variants. High heterozygosity and limited recombination also influence associated genetic load, as linked recessive deleterious mutations may be sheltered. The capture of deleterious elements in a locus under balancing selection may reinforce polymorphism by further promoting heterozygotes. Finally, according to recent genomewide scans, balanced polymorphism might be more pervasive than generally thought. We stress the need for both functional and ecological studies to characterize the evolutionary mechanisms operating in these systems.

Chromosomal rearrangements maintain a polymorphic supergene controlling butterfly mimicry.

Supergenes are tight clusters of loci that facilitate the co-segregation of adaptive variation, providing integrated control of complex adaptive phenotypes. Polymorphic supergenes, in which specific combinations of traits are maintained within a single population, were first described for 'pin' and 'thrum' floral types in Primula and Fagopyrum, but classic examples are also found in insect mimicry and snail morphology. Understanding the evolutionary mechanisms that generate these co-adapted gene sets, as well as the mode of limiting the production of unfit recombinant forms, remains a substantial challenge. Here we show that individual wing-pattern morphs in the polymorphic mimetic butterfly Heliconius numata are associated with different genomic rearrangements at the supergene locus P. These rearrangements tighten the genetic linkage between at least two colour-pattern loci that are known to recombine in closely related species, with complete suppression of recombination being observed in experimental crosses across a 400-kilobase interval containing at least 18 genes. In natural populations, notable patterns of linkage disequilibrium (LD) are observed across the entire P region. The resulting divergent haplotype clades and inversion breakpoints are found in complete association with wing-pattern morphs. Our results indicate that allelic combinations at known wing-patterning loci have become locked together in a polymorphic rearrangement at the P locus, forming a supergene that acts as a simple switch between complex adaptive phenotypes found in sympatry. These findings highlight how genomic rearrangements can have a central role in the coexistence of adaptive phenotypes involving several genes acting in concert, by locally limiting recombination and gene flow.

Selection of Valid Reference Genes for Reverse Transcription Quantitative PCR Analysis in Heliconius numata (Lepidoptera: Nymphalidae).

Identifying the genetic basis of adaptive variation is challenging in non-model organisms and quantitative real time PCR. is a useful tool for validating predictions regarding the expression of candidate genes. However, comparing expression levels in different conditions requires rigorous experimental design and statistical analyses. Here, we focused on the neotropical passion-vine butterflies Heliconius, non-model species studied in evolutionary biology for their adaptive variation in wing color patterns involved in mimicry and in the signaling of their toxicity to predators. We aimed at selecting stable reference genes to be used for normalization of gene expression data in RT-qPCR analyses from developing wing discs according to the minimal guidelines described in Minimum Information for publication of Quantitative Real-Time PCR Experiments (MIQE). To design internal RT-qPCR controls, we studied the stability of expression of nine candidate reference genes (actin, annexin, eF1α, FK506BP, PolyABP, PolyUBQ, RpL3, RPS3A, and tubulin) at two developmental stages (prepupal and pupal) using three widely used programs (GeNorm, NormFinder and BestKeeper). Results showed that, despite differences in statistical methods, genes RpL3, eF1α, polyABP, and annexin were stably expressed in wing discs in late larval and pupal stages of Heliconius numata This combination of genes may be used as a reference for a reliable study of differential expression in wings for instance for genes involved in important phenotypic variation, such as wing color pattern variation. Through this example, we provide general useful technical recommendations as well as relevant statistical strategies for evolutionary biologists aiming to identify candidate-genes involved adaptive variation in non-model organisms.

Origin and maintenance of large ribosomal RNA gene repeat size in mammals.

The genes encoding ribosomal RNA are highly conserved across life and in almost all eukaryotes are present in large tandem repeat arrays called the rDNA. rDNA repeat unit size is conserved across most eukaryotes but has expanded dramatically in mammals, principally through the expansion of the intergenic spacer region that separates adjacent rRNA coding regions. Here, we used long-read sequence data from representatives of the major amniote lineages to determine where in amniote evolution rDNA unit size increased. We find that amniote rDNA unit sizes fall into two narrow size classes: "normal" (∼11-20 kb) in all amniotes except monotreme, marsupial, and eutherian mammals, which have "large" (∼35-45 kb) sizes. We confirm that increases in intergenic spacer length explain much of this mammalian size increase. However, in stark contrast to the uniformity of mammalian rDNA unit size, mammalian intergenic spacers differ greatly in sequence. These results suggest a large increase in intergenic spacer size occurred in a mammalian ancestor and has been maintained despite substantial sequence changes over the course of mammalian evolution. This points to a previously unrecognized constraint on the length of the intergenic spacer, a region that was thought to be largely neutral. We finish by speculating on possible causes of this constraint.

The genome of a globally invasive passerine, the common myna, Acridotheres tristis.

In an era of global climate change, biodiversity conservation is receiving increased attention. Conservation efforts are greatly aided by genetic tools and approaches, which seek to understand patterns of genetic diversity and how they impact species health and their ability to persist under future climate regimes. Invasive species offer vital model systems in which to investigate questions regarding adaptive potential, with a particular focus on how changes in genetic diversity and effective population size interact with novel selection regimes. The common myna (Acridotheres tristis) is a globally invasive passerine and is an excellent model species for research both into the persistence of low-diversity populations and the mechanisms of biological invasion. To underpin research on the invasion genetics of this species, we present the genome assembly of the common myna. We describe the genomic landscape of this species, including genome wide allelic diversity, methylation, repeats, and recombination rate, as well as an examination of gene family evolution. Finally, we use demographic analysis to identify that some native regions underwent a dramatic population increase between the two most recent periods of glaciation, and reveal artefactual impacts of genetic bottlenecks on demographic analysis.

Vocal convergence and social proximity shape the calls of the most basal Passeriformes, New Zealand Wrens.

Despite extensive research on avian vocal learning, we still lack a general understanding of how and when this ability evolved in birds. As the closest living relatives of the earliest Passeriformes, the New Zealand wrens (Acanthisitti) hold a key phylogenetic position for furthering our understanding of the evolution of vocal learning because they share a common ancestor with two vocal learners: oscines and parrots. However, the vocal learning abilities of New Zealand wrens remain unexplored. Here, we test for the presence of prerequisite behaviors for vocal learning in one of the two extant species of New Zealand wrens, the rifleman (Acanthisitta chloris). We detect the presence of unique individual vocal signatures and show how these signatures are shaped by social proximity, as demonstrated by group vocal signatures and strong acoustic similarities among distantly related individuals in close social proximity. Further, we reveal that rifleman calls share similar phenotypic variance ratios to those previously reported in the learned vocalizations of the zebra finch, Taeniopygia guttata. Together these findings provide strong evidence that riflemen vocally converge, and though the mechanism still remains to be determined, they may also suggest that this vocal convergence is the result of rudimentary vocal learning abilities.

Fine-scale survey of X chromosome copy number variants and indels underlying intellectual disability.

Copy number variants and indels in 251 families with evidence of X-linked intellectual disability (XLID) were investigated by array comparative genomic hybridization on a high-density oligonucleotide X chromosome array platform. We identified pathogenic copy number variants in 10% of families, with mutations ranging from 2 kb to 11 Mb in size. The challenge of assessing causality was facilitated by prior knowledge of XLID-associated genes and the ability to test for cosegregation of variants with disease through extended pedigrees. Fine-scale analysis of rare variants in XLID families leads us to propose four additional genes, PTCHD1, WDR13, FAAH2, and GSPT2, as candidates for XLID causation and the identification of further deletions and duplications affecting X chromosome genes but without apparent disease consequences. Breakpoints of pathogenic variants were characterized to provide insight into the underlying mutational mechanisms and indicated a predominance of mitotic rather than meiotic events. By effectively bridging the gap between karyotype-level investigations and X chromosome exon resequencing, this study informs discussion of alternative mutational mechanisms, such as noncoding variants and non-X-linked disease, which might explain the shortfall of mutation yield in the well-characterized International Genetics of Learning Disability (IGOLD) cohort, where currently disease remains unexplained in two-thirds of families.

Assembly of female and male hihi genomes (stitchbird; Notiomystis cincta) enables characterization of the W chromosome and resources for conservation genomics.

A high-quality reference genome can be a valuable resource for threatened species by providing a foundation to assess their evolutionary potential to adapt to future pressures such as environmental change. We assembled the genome of a female hihi (Notiomysits cincta), a threatened passerine bird endemic to Aotearoa New Zealand. The assembled genome is 1.06 Gb, and is of high quality and highly contiguous, with a contig N50 of 7.0 Mb, estimated QV of 44 and a BUSCO completeness of 96.8%. A male assembly of comparable quality was generated in parallel. A population linkage map was used to scaffold the autosomal contigs into chromosomes. Female and male sequence coverage and comparative genomics analyses were used to identify Z-, and W-linked contigs. In total, 94.6% of the assembly length was assigned to putative nuclear chromosome scaffolds. Native DNA methylation was highly correlated between sexes, with the W chromosome contigs more highly methylated than autosomal chromosomes and Z contigs. 43 differentially methylated regions were identified, and these may represent interesting candidates for the establishment or maintenance of sex differences. By generating a high-quality reference assembly of the heterogametic sex, we have created a resource that enables characterization of genome-wide diversity and facilitates the investigation of female-specific evolutionary processes. The reference genomes will form the basis for fine-scale assessment of the impacts of low genetic diversity and inbreeding on the adaptive potential of the species and will therefore enable tailored and informed conservation management of this threatened taonga (treasured) species.

Brassinosteroid coordinates cell layer interactions in plants via cell wall and tissue mechanics.

Growth coordination between cell layers is essential for development of most multicellular organisms. Coordination may be mediated by molecular signaling and/or mechanical connectivity between cells, but how genes modify mechanical interactions between layers is unknown. Here we show that genes driving brassinosteroid synthesis promote growth of internal tissue, at least in part, by reducing mechanical epidermal constraint. We identified a brassinosteroid-deficient dwarf mutant in the aquatic plant Utricularia gibba with twisted internal tissue, likely caused by mechanical constraint from a slow-growing epidermis. We tested this hypothesis by showing that a brassinosteroid mutant in Arabidopsis enhances epidermal crack formation, indicative of increased tissue stress. We propose that by remodeling cell walls, brassinosteroids reduce epidermal constraint, showing how genes can control growth coordination between layers by means of mechanics.

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.