The Hypolimnas misippus genome supports a common origin of the W chromosome in Lepidoptera

Moths and butterflies (Lepidoptera) have a heterogametic sex chromosome system with females carrying ZW chromosomes and males ZZ. The lack of W chromosomes in early diverging lepidopteran lineages has led to the suggestion of an ancestral Z0 system in this clade and a B chromosome origin of the W. This contrasts with the canonical model of W chromosome evolution in which the W would have originated from the same homologous autosomal pair as the Z chromosome. Despite the distinct models proposed, the rapid evolution of the W chromosome has hindered the elucidation of its origin. Here, we present high-quality, chromosome-level genome assemblies of two Hypolimnas species (Hypolimnas bolina and Hypolimnas misippus) and use the H. misippus assembly to explore the evolution of W chromosomes in butterflies and moths. We show that in H. misippus the W chromosome has higher similarity to the Z chromosome than any other chromosome, which could suggest a possible origin from the same homologous autosome pair as the Z chromosome. However, using genome assemblies of closely related species (ditrysian lineages) containing assembled W chromosomes, we present contrasting evidence suggesting that the W chromosome might have evolved from a B chromosome instead. Crucially, by using a synteny analysis to infer homology, we show that W chromosomes are likely to share a common evolutionary origin in Lepidoptera. This study highlights the difficulty of studying the evolution of W chromosomes and contributes to better understanding its evolutionary origins. Significance Butterflies and moths have a sex determination system in which females carry two different sex chromosomes, Z and W, while males carry two copies of the Z. The evolutionary origin of the W chromosome has been elusive, with many possible scenarios being suggested, such as the independent evolution of W chromosomes in many butterfly and moth species. Here, we present genome assemblies of two Hypolimnas butterfly species and use one of them to shed light on the evolution of the W chromosome. We show that W chromosomes across butterflies and moths are very similar which suggests a shared common origin.

[1]  J. Mappes,et al.  Colour polymorphism associated with a gene duplication in male wood tiger moths , 2023, bioRxiv.

[2]  Simon H. Martin,et al.  Genome assembly of Danaus chrysippus and comparison with the Monarch Danaus plexippus , 2021, bioRxiv.

[3]  Simon H. Martin,et al.  The Dryas iulia Genome Supports Multiple Gains of a W Chromosome from a B Chromosome in Butterflies , 2021, Genome biology and evolution.

[4]  Travis J. Wheeler,et al.  The Dfam community resource of transposable element families, sequence models, and genome annotations , 2020, Mobile DNA.

[5]  Haowen Zhang,et al.  Haplotype-resolved de novo assembly with phased assembly graphs , 2020, 2008.01237.

[6]  Wenting Wan,et al.  Chromosome‐level reference genome assembly and gene editing of the dead‐leaf butterfly Kallima inachus , 2020, Molecular ecology resources.

[7]  Shane A. McCarthy,et al.  A haplotype-resolved, de novo genome assembly for the wood tiger moth (Arctia plantaginis) through trio binning , 2020, bioRxiv.

[8]  J. Tyson.  Bead-free long fragment LSK109 library preparation v1 , 2020, protocols.io.

[9]  Pasi Rastas,et al.  Lep-Anchor: automated construction of linkage map anchored haploid genomes , 2020, Bioinform..

[10]  M. Dalíková,et al.  Absence of W Chromosome in Psychidae Moths and Implications for the Theory of Sex Chromosome Evolution in Lepidoptera , 2019, Genes.

[11]  Mark Blaxter,et al.  BlobToolKit – Interactive Quality Assessment of Genome Assemblies , 2019, G3: Genes, Genomes, Genetics.

[12]  Emmanuel F. A. Toussaint,et al.  Phylogenomics reveals the evolutionary timing and pattern of butterflies and moths , 2019, Proceedings of the National Academy of Sciences.

[13]  Heng Li,et al.  Fast and accurate long-read assembly with wtdbg2 , 2019, Nature Methods.

[14]  Christophe Klopp,et al.  D-GENIES: dot plot large genomes in an interactive, efficient and simple way , 2018, PeerJ.

[15]  David J. Lohman,et al.  Evolution of Hypolimnas butterflies (Nymphalidae): Out-of-Africa origin and Wolbachia-mediated introgression. , 2018, Molecular phylogenetics and evolution.

[16]  Gerard Talavera,et al.  A Comprehensive and Dated Phylogenomic Analysis of Butterflies , 2018, Current Biology.

[17]  Sanjit S. Batra,et al.  The Juicebox Assembly Tools module facilitates de novo assembly of mammalian genomes with chromosome-length scaffolds for under $1000 , 2018, bioRxiv.

[18]  J. Quick.  Ultra-long read sequencing protocol for RAD004 , 2018 .

[19]  Pasi Rastas,et al.  Lep-MAP3: robust linkage mapping even for low-coverage whole genome sequencing data , 2017, Bioinform..

[20]  Beatriz Vicoso,et al.  The deep conservation of the Lepidoptera Z chromosome suggests a non-canonical origin of the W , 2017, Nature Communications.

[21]  Neva C. Durand,et al.  De novo assembly of the Aedes aegypti genome using Hi-C yields chromosome-length scaffolds , 2016, Science.

[22]  Son K. Pham,et al.  Chromosome assembly of large and complex genomes using multiple references , 2016, bioRxiv.

[23]  J. Mank,et al.  How to make a sex chromosome , 2016, Nature Communications.

[24]  James T. Robinson,et al.  Juicebox Provides a Visualization System for Hi-C Contact Maps with Unlimited Zoom. , 2016, Cell systems.

[25]  Neva C. Durand,et al.  Juicer Provides a One-Click System for Analyzing Loop-Resolution Hi-C Experiments. , 2016, Cell systems.

[26]  Blake A. Simmons,et al.  MaxBin 2.0: an automated binning algorithm to recover genomes from multiple metagenomic datasets , 2016, Bioinform..

[27]  Evgeny M. Zdobnov,et al.  BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs , 2015, Bioinform..

[28]  Åsa K. Björklund,et al.  Tn5 transposase and tagmentation procedures for massively scaled sequencing projects , 2014, Genome research.

[29]  Christina A. Cuomo,et al.  Pilon: An Integrated Tool for Comprehensive Microbial Variant Detection and Genome Assembly Improvement , 2014, PloS one.

[30]  Roland Eils,et al.  circlize implements and enhances circular visualization in R , 2014, Bioinform..

[31]  M. Kirkpatrick,et al.  Sex Determination: Why So Many Ways of Doing It? , 2014, PLoS biology.

[32]  Brian J. Raney,et al.  Ragout—a reference-assisted assembly tool for bacterial genomes , 2014, Bioinform..

[33]  Robert Gentleman,et al.  Software for Computing and Annotating Genomic Ranges , 2013, PLoS Comput. Biol..

[34]  F. Marec,et al.  Rapid turnover of the W chromosome in geographical populations of wild silkmoths, Samia cynthia ssp. , 2013, Chromosome Research.

[35]  David Haussler,et al.  HAL: a hierarchical format for storing and analyzing multiple genome alignments , 2013, Bioinform..

[36]  Heng Li.  Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM , 2013, 1303.3997.

[37]  M. Kirkpatrick,et al.  Are all sex chromosomes created equal? , 2011, Trends in genetics : TIG.

[38]  Masakatsu Watanabe,et al.  B Chromosomes Have a Functional Effect on Female Sex Determination in Lake Victoria Cichlid Fishes , 2011, PLoS genetics.

[39]  David Haussler,et al.  Cactus: Algorithms for genome multiple sequence alignment. , 2011, Genome research.

[40]  Miriah D. Meyer,et al.  Genome-wide synteny through highly sensitive sequence alignment: Satsuma , 2010, Bioinform..

[41]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[42]  G. Roderick,et al.  The joint evolutionary histories of Wolbachia and mitochondria in Hypolimnas bolina , 2009, BMC Evolutionary Biology.

[43]  F. Marec,et al.  Molecular divergence of the W chromosomes in pyralid moths (Lepidoptera) , 2007, Chromosome Research.

[44]  G. Roderick,et al.  Evolution of Male-Killer Suppression in a Natural Population , 2006, PLoS biology.

[45]  W. J. Kent,et al.  BLAT--the BLAST-like alignment tool. , 2002, Genome research.

[46]  E. Dyson,et al.  Wolbachia infection associated with all-female broods in Hypolimnas bolina (Lepidoptera: Nymphalidae): evidence for horizontal transmission of a butterfly male killer , 2002, Heredity.

[47]  V. Lukhtanov.  Sex chromatin and sex chromosome systems in nonditrysian Lepidoptera (Insecta) , 2000 .

[48]  S. Wachtel,et al.  On the evolution of sex determination , 1987, Proceedings of the Royal Society of London. Series B. Biological Sciences.

[49]  R. Vane-Wright,et al.  The polymorphism, mimicry, and host plant relationships of Hypolimnas butterflies , 1977 .

[50]  M. Rothschild,et al.  Hypolimnas bolina (L.), a mimic of danaid butterflies, and its model Euploea core (Cram.) store cardioactive substances , 1977, Nature.

[51]  D. A. Smith.  Phenotypic diversity, mimicry and natural selection in the African butterfly Hypolimnas misippus L. (Lepidoptera: Nymphalidae) , 1976 .

[52]  P. Sheppard,et al.  Absence of crossing-over in female butterflies (Heliconius) , 1975, Heredity.

[53]  Mosè Manni,et al.  BUSCO: Assessing Genome Assembly and Annotation Completeness. , 2019, Methods in molecular biology.

[54]  J. Keilwagen,et al.  GeMoMa: Homology-Based Gene Prediction Utilizing Intron Position Conservation and RNA-seq Data. , 2019, Methods in molecular biology.

[55]  Jill P. Mesirov,et al.  Juicebox Provides a Visualization System for HiC Contact Maps with Unlimited Zoom Citation , 2016 .

[56]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[57]  K. Sahara,et al.  Sex chromosome evolution in moths and butterflies , 2011, Chromosome Research.

[58]  What Are All Those Funny Symbols in a Blast Printout? Blast = Basic Local Alignment Search Tool , 2022 .