Single amino acid mutation in an ATP-binding cassette transporter gene causes resistance to Bt toxin Cry1Ab in the silkworm, Bombyx mori.

Bt toxins derived from the arthropod bacterial pathogen Bacillus thuringiensis are widely used for insect control as insecticides or in transgenic crops. Bt resistance has been found in field populations of several lepidopteran pests and in laboratory strains selected with Bt toxin. Widespread planting of crops expressing Bt toxins has raised concerns about the potential increase of resistance mutations in targeted insects. By using Bombyx mori as a model, we identified a candidate gene for a recessive form of resistance to Cry1Ab toxin on chromosome 15 by positional cloning. BGIBMGA007792-93, which encodes an ATP-binding cassette transporter similar to human multidrug resistance protein 4 and orthologous to genes associated with recessive resistance to Cry1Ac in Heliothis virescens and two other lepidopteran species, was expressed in the midgut. Sequences of 10 susceptible and seven resistant silkworm strains revealed a common tyrosine insertion in an outer loop of the predicted transmembrane structure of resistant alleles. We confirmed the role of this ATP-binding cassette transporter gene in Bt resistance by converting a resistant silkworm strain into a susceptible one by using germline transformation. This study represents a direct demonstration of Bt resistance gene function in insects with the use of transgenesis.

Comparison of phenotypic value changes in pure lines of Bombyx mori (Lepidoptera: Bombycidae) during consecutive generations following initial selection on cocoon weight.

The experiments reported here were conducted to investigate the effect of selection on three quantitative traits, namely cocoon weight, cocoon shell weight, and cocoon shell percentage, during four generations by rearing six pure breeds of domesticated silkworm, Bombyx mori L. (Lepidoptera: Bombycidae) of Chinese and Japanese origin compared with random unselected groups as controls. All stages of rearing and data recording were performed over four rearing periods, with generations 1-3 during successive spring seasons and generation 4 during the autumn season in year 3. Each pure line contained two groups of selected and random (control) groups. Comparisons included the effect of selection methods, pure line, and generation on the phenotypic values. We found strong main effects of pure line, generation, sex, and group and support for nearly all interactions between these main effects for all three response traits. The results indicated that cocoon weight and cocoon shell weight in the selected group were higher than in the control or nonselected group. Both selected and nonselected groups had the lowest cocoon weight, cocoon shell weight, and cocoon shell percentage in the fourth generation when environmental conditions during the autumn season were less favorable than spring. The cocoon weight and cocoon shell weight averages were higher for nonselected groups in the second and third generations, and for the selected group in the first generation due to the direct effect of selection.

Repression of tyrosine hydroxylase is responsible for the sex-linked chocolate mutation of the silkworm, Bombyx mori.

Pigmentation patterning has long interested biologists, integrating topics in ecology, development, genetics, and physiology. Wild-type neonatal larvae of the silkworm, Bombyx mori, are completely black. By contrast, the epidermis and head of larvae of the homozygous recessive sex-linked chocolate (sch) mutant are reddish brown. When incubated at 30 degrees C, mutants with the sch allele fail to hatch; moreover, homozygous mutants carrying the allele sch lethal (sch(l)) do not hatch even at room temperature (25 degrees C). By positional cloning, we narrowed a region containing sch to 239,622 bp on chromosome 1 using 4,501 backcross (BC1) individuals. Based on expression analyses, the best sch candidate gene was shown to be tyrosine hydroxylase (BmTh). BmTh coding sequences were identical among sch, sch(l), and wild-type. However, in sch the approximately 70-kb sequence was replaced with approximately 4.6 kb of a Tc1-mariner type transposon located approximately 6 kb upstream of BmTh, and in sch(l), a large fragment of an L1Bm retrotransposon was inserted just in front of the transcription start site of BmTh. In both cases, we observed a drastic reduction of BmTh expression. Use of RNAi with BmTh prevented pigmentation and hatching, and feeding of a tyrosine hydroxylase inhibitor also suppressed larval pigmentation in the wild-type strain, pnd(+) and in a pS (black-striped) heterozygote. Feeding L-dopa to sch neonate larvae rescued the mutant phenotype from chocolate to black. Our results indicate the BmTh gene is responsible for the sch mutation, which plays an important role in melanin synthesis producing neonatal larval color.

Advances in the Arms Race Between Silkworm and Baculovirus.

Insects are the largest group of animals. Nearly all organisms, including insects, have viral pathogens. An important domesticated economic insect is the silkworm moth Bombyx mori. B. mori nucleopolyhedrovirus (BmNPV) is a typical baculovirus and a primary silkworm pathogen. It causes major economic losses in sericulture. Baculoviruses are used in biological pest control and as a bioreactor. Silkworm and baculovirus comprise a well-established model of insect-virus interactions. Several recent studies have focused on this model and provided novel insights into viral infections and host defense. Here, we focus on baculovirus invasion, silkworm immune response, baculovirus evasion of host immunity, and enhancement of antiviral efficacy. We also discuss major issues remaining and future directions of research on silkworm antiviral immunity. Elucidation of the interaction between silkworm and baculovirus furnishes a theoretical basis for targeted pest control, enhanced pathogen resistance in economically important insects, and bioreactor improvement.

Circadian regulation of night feeding and daytime detoxification in a formidable Asian pest Spodoptera litura.

Voracious feeding, trans-continental migration and insecticide resistance make Spodoptera litura among the most difficult Asian agricultural pests to control. Larvae exhibit strong circadian behavior, feeding actively at night and hiding in soil during daytime. The daily pattern of larval metabolism was reversed, with higher transcription levels of genes for digestion (amylase, protease, lipase) and detoxification (CYP450s, GSTs, COEs) in daytime than at night. To investigate the control of these processes, we annotated nine essential clock genes and analyzed their transcription patterns, followed by functional analysis of their coupling using siRNA knockdown of interlocked negative feedback system core and repressor genes (SlituClk, SlituBmal1 and SlituCwo). Based on phase relationships and overexpression in cultured cells the controlling mechanism seems to involve direct coupling of the circadian processes to E-boxes in responding promoters. Additional manipulations involving exposure to the neonicotinoid imidacloprid suggested that insecticide application must be based on chronotoxicological considerations for optimal effectiveness.

Lepidopteran wing scales contain abundant cross-linked film-forming histidine-rich cuticular proteins.

Scales are symbolic characteristic of Lepidoptera; however, nothing is known about the contribution of cuticular proteins (CPs) to the complex patterning of lepidopteran scales. This is because scales are resistant to solubilization, thus hindering molecular studies. Here we succeeded in dissolving developing wing scales from Bombyx mori, allowing analysis of their protein composition. We identified a distinctive class of histidine rich (His-rich) CPs (6%-45%) from developing lepidopteran scales by LC-MS/MS. Functional studies using RNAi revealed CPs with different histidine content play distinct and critical roles in constructing the microstructure of the scale surface. Moreover, we successfully synthesized films in vitro by crosslinking a 45% His-rich CP (BmorCPR152) with laccase2 using N-acetyl- dopamine or N-ß-alanyl-dopamine as the substrate. This molecular study of scales provides fundamental information about how such a fine microstructure is constructed and insights into the potential application of CPs as new biomaterials.

Construction and sequence sampling of deep-coverage, large-insert BAC libraries for three model lepidopteran species.

BACKGROUND: Manduca sexta, Heliothis virescens, and Heliconius erato represent three widely-used insect model species for genomic and fundamental studies in Lepidoptera. Large-insert BAC libraries of these insects are critical resources for many molecular studies, including physical mapping and genome sequencing, but not available to date. RESULTS: We report the construction and characterization of six large-insert BAC libraries for the three species and sampling sequence analysis of the genomes. The six BAC libraries were constructed with two restriction enzymes, two libraries for each species, and each has an average clone insert size ranging from 152-175 kb. We estimated that the genome coverage of each library ranged from 6-9 x, with the two combined libraries of each species being equivalent to 13.0-16.3 x haploid genomes. The genome coverage, quality and utility of the libraries were further confirmed by library screening using 6 approximately 8 putative single-copy probes. To provide a first glimpse into these genomes, we sequenced and analyzed the BAC ends of approximately 200 clones randomly selected from the libraries of each species. The data revealed that the genomes are AT-rich, contain relatively small fractions of repeat elements with a majority belonging to the category of low complexity repeats, and are more abundant in retro-elements than DNA transposons. Among the species, the H. erato genome is somewhat more abundant in repeat elements and simple repeats than those of M. sexta and H. virescens. The BLAST analysis of the BAC end sequences suggested that the evolution of the three genomes is widely varied, with the genome of H. virescens being the most conserved as a typical lepidopteran, whereas both genomes of H. erato and M. sexta appear to have evolved significantly, resulting in a higher level of species- or evolutionary lineage-specific sequences. CONCLUSION: The high-quality and large-insert BAC libraries of the insects, together with the identified BACs containing genes of interest, provide valuable information, resources and tools for comprehensive understanding and studies of the insect genomes and for addressing many fundamental questions in Lepidoptera. The sample of the genomic sequences provides the first insight into the constitution and evolution of the insect genomes.

An integrated genetic linkage map for silkworms with three parental combinations and its application to the mapping of single genes and QTL.

BACKGROUND: Bombyx mori, the domesticated silkworm, is a well-studied model insect with great economic and scientific significance. Although more than 400 mutations have been described in silkworms, most have not been identified, especially those affecting economically-important traits. Simple sequence repeats (SSRs) are effective and economical tools for mapping traits and genetic improvement. The current SSR linkage map is of low density and contains few polymorphisms. The purpose of this work was to develop a dense and informative linkage map that would assist in the preliminary mapping and dissection of quantitative trait loci (QTL) in a variety of silkworm strains. RESULTS: Through an analysis of > 50,000 genotypes across new mapping populations, we constructed two new linkage maps covering 27 assigned chromosomes and merged the data with previously reported data sets. The integrated consensus map contains 692 unique SSR sites, improving the density from 6.3 cM in the previous map to 4.8 cM. We also developed 497 confirmed neighboring markers for corresponding low-polymorphism sites, with 244 having polymorphisms. Large-scale statistics on the SSR type were suggestive of highly efficient markers, based upon which we searched 16,462 available genomic scaffolds for SSR loci. With the newly constructed map, we mapped single-gene traits, the QTL of filaments, and a number of ribosomal protein genes. CONCLUSION: The integrated map produced in this study is a highly efficient genetic tool for the high-throughput mapping of single genes and QTL. Compared to previous maps, the current map offers a greater number of markers and polymorphisms; thus, it may be used as a resource for marker-assisted breeding.

KAIKObase: an integrated silkworm genome database and data mining tool.

BACKGROUND: The silkworm, Bombyx mori, is one of the most economically important insects in many developing countries owing to its large-scale cultivation for silk production. With the development of genomic and biotechnological tools, B. mori has also become an important bioreactor for production of various recombinant proteins of biomedical interest. In 2004, two genome sequencing projects for B. mori were reported independently by Chinese and Japanese teams; however, the datasets were insufficient for building long genomic scaffolds which are essential for unambiguous annotation of the genome. Now, both the datasets have been merged and assembled through a joint collaboration between the two groups. DESCRIPTION: Integration of the two data sets of silkworm whole-genome-shotgun sequencing by the Japanese and Chinese groups together with newly obtained fosmid- and BAC-end sequences produced the best continuity (~3.7 Mb in N50 scaffold size) among the sequenced insect genomes and provided a high degree of nucleotide coverage (88%) of all 28 chromosomes. In addition, a physical map of BAC contigs constructed by fingerprinting BAC clones and a SNP linkage map constructed using BAC-end sequences were available. In parallel, proteomic data from two-dimensional polyacrylamide gel electrophoresis in various tissues and developmental stages were compiled into a silkworm proteome database. Finally, a Bombyx trap database was constructed for documenting insertion positions and expression data of transposon insertion lines. CONCLUSION: For efficient usage of genome information for functional studies, genomic sequences, physical and genetic map information and EST data were compiled into KAIKObase, an integrated silkworm genome database which consists of 4 map viewers, a gene viewer, and sequence, keyword and position search systems to display results and data at the level of nucleotide sequence, gene, scaffold and chromosome. Integration of the silkworm proteome database and the Bombyx trap database with KAIKObase led to a high-grade, user-friendly, and comprehensive silkworm genome database which is now available from URL: http://sgp.dna.affrc.go.jp/KAIKObase/.

Genome-wide annotation and comparative analysis of cuticular protein genes in the noctuid pest Spodoptera litura.

Insect cuticle is considered an adaptable and versatile building material with roles in the construction and function of exoskeleton. Its physical properties are varied, as the biological requirements differ among diverse structures and change during the life cycle of the insect. Although the bulk of cuticle consists basically of cuticular proteins (CPs) associated with chitin, the degree of cuticular sclerotization is an important factor in determining its physical properties. Spodoptera litura, the tobacco cutworm, is an important agricultural pest in Asia. Compared to the domestic silkworm, Bombyx mori, another lepidopteran whose CP genes have been well annotated, S. litura has a shorter life cycle, hides in soil during daytime beginning in the 5th instar and is exposed to soil in the pupal stage without the protection of a cocoon. In order to understand how the CP genes may have been adapted to support the characteristic life style of S. litura, we searched its genome and found 287 putative cuticular proteins that can be classified into 9 CP families (CPR with three groups (RR-1, RR-2, RR-3), CPAP1, CPAP3, CPF, CPFL, CPT, CPG, CPCFC and CPLCA), and a collection of unclassified CPs named CPH. There were also 112 cuticular proteins enriched in Histidine residues with content varying from 6% to 30%, comprising many more His-rich cuticular proteins than B. mori. A phylogenetic analysis between S. litura, M. sexta and B. mori uncovered large expansions of RR-1 and RR-2 CPs, forming large gene clusters in different regions of S. litura chromosome 9. We used RNA-seq analysis to document the expression profiles of CPs in different developmental stages and tissues of S. litura. The comparative genomic analysis of CPs between S. litura and B. mori integrated with the unique behavior and life cycle of the two species offers new insights into their contrasting ecological adaptations.

Patterns of male reproductive success in Crepidula fornicata provide new insight for sex allocation and optimal sex change.

The size-advantage model and sex-allocation theory are frequently invoked to explain the evolution and maintenance of sequential hermaphroditism in many taxa. A test of current theory requires quantitative estimates of reproductive success and knowledge of the relationship between reproduction and size for each gender. Reproductive success can be difficult to measure. In species where polyandry occurs, it can be quantified only by determining paternity of offspring. We employed microsatellite loci to establish paternity for 12 families of Crepidula fornicata, where a family is defined as a single female, her brood, and the males stacked on top of her. Genetic data were analyzed and paternity was assigned to a single potential father for more than 83% of the offspring tested. Estimates of reproductive success revealed that one male within the family fathered the majority of offspring and that he was usually the largest male and the one closest to the brooding female. The dominant male's success also tended to decrease as the number of mature males within the family increased. Our results suggest that sperm competition could be a driving force in determining male reproductive success and the timing of sex change in C. fornicata.

Construction of a single nucleotide polymorphism linkage map for the silkworm, Bombyx mori, based on bacterial artificial chromosome end sequences.

We have developed a linkage map for the silkworm Bombyx mori based on single nucleotide polymorphisms (SNPs) between strains p50T and C108T initially found on regions corresponding to the end sequences of bacterial artificial chromosome (BAC) clones. Using 190 segregants from a backcross of a p50T female x an F1 (p50T x C108T) male, we analyzed segregation patterns of 534 SNPs between p50T and C108T, detected among 3840 PCR amplicons, each associated with a p50T BAC end sequence. This enabled us to construct a linkage map composed of 534 SNP markers spanning 1305 cM in total length distributed over the expected 28 linkage groups. Of the 534 BACs whose ends harbored the SNPs used to construct the linkage map, 89 were associated with 107 different ESTs. Since each of the SNP markers is directly linked to a specific genomic BAC clone and to whole-genome sequence data, and some of them are also linked to EST data, the SNP linkage map will be a powerful tool for investigating silkworm genome properties, mutation mapping, and map-based cloning of genes of industrial and agricultural interest.

Expression map of a complete set of gustatory receptor genes in chemosensory organs of Bombyx mori.

Most lepidopteran species are herbivores, and interaction with host plants affects their gene expression and behavior as well as their genome evolution. Gustatory receptors (Grs) are expected to mediate host plant selection, feeding, oviposition and courtship behavior. However, due to their high diversity, sequence divergence and extremely low level of expression it has been difficult to identify precisely a complete set of Grs in Lepidoptera. By manual annotation and BAC sequencing, we improved annotation of 43 gene sequences compared with previously reported Grs in the most studied lepidopteran model, the silkworm, Bombyx mori, and identified 7 new tandem copies of BmGr30 on chromosome 7, bringing the total number of BmGrs to 76. Among these, we mapped 68 genes to chromosomes in a newly constructed chromosome distribution map and 8 genes to scaffolds; we also found new evidence for large clusters of BmGrs, especially from the bitter receptor family. RNA-seq analysis of diverse BmGr expression patterns in chemosensory organs of larvae and adults enabled us to draw a precise organ specific map of BmGr expression. Interestingly, most of the clustered genes were expressed in the same tissues and more than half of the genes were expressed in larval maxillae, larval thoracic legs and adult legs. For example, BmGr63 showed high expression levels in all organs in both larval and adult stages. By contrast, some genes showed expression limited to specific developmental stages or organs and tissues. BmGr19 was highly expressed in larval chemosensory organs (especially antennae and thoracic legs), the single exon genes BmGr53 and BmGr67 were expressed exclusively in larval tissues, the BmGr27-BmGr31 gene cluster on chr7 displayed a high expression level limited to adult legs and the candidate CO2 receptor BmGr2 was highly expressed in adult antennae, where few other Grs were expressed. Transcriptional analysis of the Grs in B. mori provides a valuable new reference for finding genes involved in plant-insect interactions in Lepidoptera and establishing correlations between these genes and vital insect behaviors like host plant selection and courtship for mating.

Genomic adaptation to polyphagy and insecticides in a major East Asian noctuid pest.

The tobacco cutworm, Spodoptera litura, is among the most widespread and destructive agricultural pests, feeding on over 100 crops throughout tropical and subtropical Asia. By genome sequencing, physical mapping and transcriptome analysis, we found that the gene families encoding receptors for bitter or toxic substances and detoxification enzymes, such as cytochrome P450, carboxylesterase and glutathione-S-transferase, were massively expanded in this polyphagous species, enabling its extraordinary ability to detect and detoxify many plant secondary compounds. Larval exposure to insecticidal toxins induced expression of detoxification genes, and knockdown of representative genes using short interfering RNA (siRNA) reduced larval survival, consistent with their contribution to the insect's natural pesticide tolerance. A population genetics study indicated that this species expanded throughout southeast Asia by migrating along a South India-South China-Japan axis, adapting to wide-ranging ecological conditions with diverse host plants and insecticides, surviving and adapting with the aid of its expanded detoxification systems. The findings of this study will enable the development of new pest management strategies for the control of major agricultural pests such as S. litura.

Multifaceted biological insights from a draft genome sequence of the tobacco hornworm moth, Manduca sexta.

Manduca sexta, known as the tobacco hornworm or Carolina sphinx moth, is a lepidopteran insect that is used extensively as a model system for research in insect biochemistry, physiology, neurobiology, development, and immunity. One important benefit of this species as an experimental model is its extremely large size, reaching more than 10 g in the larval stage. M. sexta larvae feed on solanaceous plants and thus must tolerate a substantial challenge from plant allelochemicals, including nicotine. We report the sequence and annotation of the M. sexta genome, and a survey of gene expression in various tissues and developmental stages. The Msex_1.0 genome assembly resulted in a total genome size of 419.4 Mbp. Repetitive sequences accounted for 25.8% of the assembled genome. The official gene set is comprised of 15,451 protein-coding genes, of which 2498 were manually curated. Extensive RNA-seq data from many tissues and developmental stages were used to improve gene models and for insights into gene expression patterns. Genome wide synteny analysis indicated a high level of macrosynteny in the Lepidoptera. Annotation and analyses were carried out for gene families involved in a wide spectrum of biological processes, including apoptosis, vacuole sorting, growth and development, structures of exoskeleton, egg shells, and muscle, vision, chemosensation, ion channels, signal transduction, neuropeptide signaling, neurotransmitter synthesis and transport, nicotine tolerance, lipid metabolism, and immunity. This genome sequence, annotation, and analysis provide an important new resource from a well-studied model insect species and will facilitate further biochemical and mechanistic experimental studies of many biological systems in insects.

Development of genetic sexing strains in Lepidoptera: from traditional to transgenic approaches.

The sterile insect technique (SIT) is currently being used for the control of many agricultural pests, including some lepidopteran species. The SIT relies on the rearing and release of large numbers of genetically sterile insects into a wild population. The holokinetic chromosomes of Lepidoptera respond differently to radiation than do species where there is a localized centromere. This difference has enabled a variation of the SIT to be developed for Lepidoptera where a substerilizing dose of radiation is given to the insects before their release with the result that a certain level of sterility is inherited by the F1 offspring. The development of genetic sexing strains for fruit flies, enabling the release of males only, has resulted in enormous economic benefits in the mass rearing and has increased the efficiency of the field operations severalfold. This article outlines Mendelian approaches that are currently available to separate large numbers of males and females efficiently for different lepidopteran species and describes their difficulties and constraints. Successful transgenesis in several lepidopteran species opens up new possibilities to develop genetic sexing strains. The proposal to develop genetic sexing strains described in this article takes advantage of the fact that in Lepidoptera, the female is the heterogametic sex, with most species having aWZ sex chromosome pair, whereas the males are ZZ. This means that if a conditional lethal gene can be inserted into the W chromosome, then all females should die after the application of the restrictive condition. The assumptions made to accommodate this model are discussed, and the advantages to be gained for control programs are elucidated.

Construction, complete sequence, and annotation of a BAC contig covering the silkworm chorion locus.

The silkmoth chorion was studied extensively by F.C. Kafatos' group for almost 40 years. However, the complete structure of the chorion locus was not obtained in the genome sequence of Bombyx mori published in 2008 due to repetitive sequences, resulting in gaps and an incomplete view of the locus. To obtain the complete sequence of the chorion locus, expressed sequence tags (ESTs) derived from follicular epithelium cells were used as probes to screen a bacterial artificial chromosome (BAC) library. Seven BACs were selected to construct a contig which covered the whole chorion locus. By Sanger sequencing, we successfully obtained complete sequences of the chorion locus spanning 871,711 base pairs on chromosome 2, where we annotated 127 chorion genes. The dataset reported here will recruit more researchers to revisit one of the oldest model systems which has been used to study developmentally regulated gene expression. It also provides insights into egg development and fertilization mechanisms and is relevant to applications related to improvements in breeding procedures and transgenesis.

A comprehensive analysis of the chorion locus in silkmoth.

Despite more than 40 years of intense study, essential features of the silkmoth chorion (eggshell) are still not fully understood. To determine the precise structure of the chorion locus, we performed extensive EST analysis, constructed a bacterial artificial chromosome (BAC) contig, and obtained a continuous genomic sequence of 871,711 base pairs. We annotated 127 chorion genes in two segments interrupted by a 164 kb region with 5 non-chorion genes, orthologs of which were on chorion bearing scaffolds in 4 ditrysian families. Detailed transcriptome analysis revealed expression throughout choriogenesis of most chorion genes originally categorized as "middle", and evidence for diverse regulatory mechanisms including cis-elements, alternative splicing and promoter utilization, and antisense RNA. Phylogenetic analysis revealed multigene family associations and faster evolution of early chorion genes and transcriptionally active pseudogenes. Proteomics analysis identified 99 chorion proteins in the eggshell and micropyle localization of 1 early and 6 Hc chorion proteins.

Chromatin-induced spindle assembly plays an important role in metaphase congression of silkworm holocentric chromosomes.

The kinetochore plays important roles in cell cycle progression. Interactions between chromosomes and spindle microtubules allow chromosomes to congress to the middle of the cell and to segregate the sister chromatids into daughter cells in mitosis. The chromosome passenger complex (CPC), composed of the Aurora B kinase and its regulatory subunits INCENP, Survivin, and Borealin, plays multiple roles in these chromosomal events. In the genome of the silkworm, Bombyx mori, which has holocentric chromosomes, the CPC components and their molecular interactions were highly conserved. In contrast to monocentric species, however, the silkworm CPC co-localized with the chromatin-driven spindles on the upper side of prometaphase chromosomes without forming bipolar mitotic spindles. Depletion of the CPC by RNAi arrested the cell cycle progression at prometaphase and disrupted the microtubule network of the chromatin-driven spindles. Interestingly, depletion of mitotic centromere-associated kinesin (MCAK) recovered formation of the microtubule network but did not overcome the cell cycle arrest at prometaphase. These results suggest that the CPC modulates the chromatin-induced spindle assembly and metaphase congression of silkworm holocentric chromosomes.

The Glanville fritillary genome retains an ancient karyotype and reveals selective chromosomal fusions in Lepidoptera.

Previous studies have reported that chromosome synteny in Lepidoptera has been well conserved, yet the number of haploid chromosomes varies widely from 5 to 223. Here we report the genome (393 Mb) of the Glanville fritillary butterfly (Melitaea cinxia; Nymphalidae), a widely recognized model species in metapopulation biology and eco-evolutionary research, which has the putative ancestral karyotype of n=31. Using a phylogenetic analyses of Nymphalidae and of other Lepidoptera, combined with orthologue-level comparisons of chromosomes, we conclude that the ancestral lepidopteran karyotype has been n=31 for at least 140 My. We show that fusion chromosomes have retained the ancestral chromosome segments and very few rearrangements have occurred across the fusion sites. The same, shortest ancestral chromosomes have independently participated in fusion events in species with smaller karyotypes. The short chromosomes have higher rearrangement rate than long ones. These characteristics highlight distinctive features of the evolutionary dynamics of butterflies and moths.