Unprecedented reorganization of holocentric chromosomes provides insights into the enigma of lepidopteran chromosome evolution

Chromosome-level genome assembly and comparative genomics shed light on Helicoverpa assulta ecology and pest management

BACKGROUND

The Oriental tobacco budworm, Helicoverpa assulta, a specialist herbivorous insect that exclusively feeds on plants of the Solanaceae family, causes considerable damage to crops, such as tobacco and hot pepper. The absence of a genome sequence for this species hinders further research on its pest management and ecological adaptation.

RESULTS

Here, we present a high-quality chromosome-level genome of a Korean strain of H. assulta (Pyeongchang strain, K18). The total assembly spans 424.4 Mb with an N50 length of 14.54 Mb and 37% GC content. The assembled genome (ASM2961881v1) comprises 31 chromosomes, similar to other congeneric generalist species including H. armigera and H. zea. In terms of genomic assembly quality, the complete BUSCOs and repeat content accounted for 98.3% and 33.01% of the genome, respectively. Based on this assembly, 19 485 protein-coding genes were predicted in the genome annotation. A comparative analysis was conducted using the identified number of protein-coding genes in H. armigera (24154) and H. zea (23696). Out of the 19 485 predicted genes, 137 genes in 15 orthogroups were found to have expanded significantly in H. assulta, while 149 genes in 95 orthogroups contracted rapidly.

CONCLUSION

This study revealed specific gene expansions and contractions in H. assulta compared to those in its close relatives, indicating potential adaptations related to its specialized feeding habits. Also, the comparative genome analysis provides valuable insights for the integrated pest management of H. assulta and other globally significant pests in the Heliothinae subfamily. © 2024 Society of Chemical Industry.

Genome organization: Raison d'être of ancestral linkage groups

The genomes of extant organisms contain conserved blocks of regions that can be traced back to ancient ancestors, yet the evolutionary pressures that maintained such genomic segments remain unclear. New research on a curious organism with two different genomes sheds light on why our genomes are organized as they are.

De Novo Genome Assembly and Annotation of Leptosia nina Provide New Insights into the Evolutionary Dynamics of Genes Involved in Host-Plant Adaptation of Pierinae Butterflies

In interactions between plants and herbivorous insects, the traits enabling phytophagous insects to overcome chemical defenses of their host plants have evolved multiple times. A prominent example of such adaptive key innovations in herbivorous insects is nitrile specifier proteins (NSPs) that enabled Pierinae butterflies to colonize Brassicales host plants that have a glucosinolate-myrosinase defense system. Although the evolutionary aspects of NSP-encoding genes have been studied in some Pierinae taxa (especially among Pieris butterflies), the ancestral evolutionary state of NSPs is unclear due to the limited genomic information available for species within Pierinae. Here, we generate a high-quality genome assembly and annotation of Leptosia nina, a member of a small tribe, Leptosiaini. L. nina uses as its main host Capparaceae plants, one of the ancestral hosts within Pierinae. By using ∼90-fold coverage of Oxford Nanopore long reads and Illumina short reads for subsequent polishing and error correction, we constructed a final genome assembly that consisted of 286 contigs with a total of 225.8 Mb and an N50 of 10.7 Mb. Genome annotation with transcriptome hints predicted 16,574 genes and covered 98.3% of BUSCO genes. A typical NSP gene is composed of three tandem domains found in Pierinae butterflies; unexpectedly, we found a new NSP-like gene in Pierinae composed of only two tandem domains. This newly found NSP-like gene in L. nina provides important insights into the evolutionary dynamics of domain and gene duplication events relating to host-plant adaptation in Pierinae butterflies.

Extreme genome scrambling in marine planktonic Oikopleura dioica cryptic species

Genome structural variations within species are rare. How selective constraints preserve gene order and chromosome structure is a central question in evolutionary biology that remains unsolved. Our sequencing of several genomes of the appendicularian tunicate Oikopleura dioica around the globe reveals extreme genome scrambling caused by thousands of chromosomal rearrangements, although showing no obvious morphological differences between these animals. The breakpoint accumulation rate is an order of magnitude higher than in ascidian tunicates, nematodes, Drosophila, or mammals. Chromosome arms and sex-specific regions appear to be the primary unit of macrosynteny conservation. At the microsyntenic level, scrambling did not preserve operon structures, suggesting an absence of selective pressure to maintain them. The uncoupling of the genome scrambling with morphological conservation in O. dioica suggests the presence of previously unnoticed cryptic species and provides a new biological system that challenges our previous vision of speciation in which similar animals always share similar genome structures.

A macroevolutionary role for chromosomal fusion and fission in Erebia butterflies

The impact of large-scale chromosomal rearrangements, such as fusions and fissions, on speciation is a long-standing conundrum. We assessed whether bursts of change in chromosome numbers resulting from chromosomal fusion or fission are related to increased speciation rates in Erebia, one of the most species-rich and karyotypically variable butterfly groups. We established a genome-based phylogeny and used state-dependent birth-death models to infer trajectories of karyotype evolution. We demonstrated that rates of anagenetic chromosomal changes (i.e., along phylogenetic branches) exceed cladogenetic changes (i.e., at speciation events), but, when cladogenetic changes occur, they are mostly associated with chromosomal fissions rather than fusions. We found that the relative importance of fusion and fission differs among Erebia clades of different ages and that especially in younger, more karyotypically diverse clades, speciation is more frequently associated with cladogenetic chromosomal changes. Overall, our results imply that chromosomal fusions and fissions have contrasting macroevolutionary roles and that large-scale chromosomal rearrangements are associated with bursts of species diversification.

Comparative genomics reveals the dynamics of chromosome evolution in Lepidoptera

Chromosomes are a central unit of genome organization. One-tenth of all described species on Earth are butterflies and moths, the Lepidoptera, which generally possess 31 chromosomes. However, some species display dramatic variation in chromosome number. Here we analyse 210 chromosomally complete lepidopteran genomes and show that the chromosomes of extant lepidopterans are derived from 32 ancestral linkage groups, which we term Merian elements. Merian elements have remained largely intact through 250 million years of evolution and diversification. Against this stable background, eight lineages have undergone extensive reorganization either through numerous fissions or a combination of fusion and fission events. Outside these lineages, fusions are rare and fissions are rarer still. Fusions often involve small, repeat-rich Merian elements and the sex-linked element. Our results reveal the constraints on genome architecture in Lepidoptera and provide a deeper understanding of chromosomal rearrangements in eukaryotic genome evolution.

Evidence of Nonrandom Patterns of Functional Chromosome Organization in Danaus plexippus

Our understanding on the interplay between gene functionality and gene arrangement at different chromosome scales relies on a few Diptera and the honeybee, species with quality reference genome assemblies, accurate gene annotations, and abundant transcriptome data. Using recently generated 'omic resources in the monarch butterfly Danaus plexippus, a species with many more and smaller chromosomes relative to Drosophila species and the honeybee, we examined the organization of genes preferentially expressed at broadly defined developmental stages (larva, pupa, adult males, and adult females) at both fine and whole-chromosome scales. We found that developmental stage-regulated genes do not form more clusters, but do form larger clusters, than expected by chance, a pattern consistent across the gene categories examined. Notably, out of the 30 chromosomes in the monarch genome, 12 of them, plus the fraction of the chromosome Z that corresponds to the ancestral Z in other Lepidoptera, were found enriched for developmental stage-regulated genes. These two levels of nonrandom gene organization are not independent as enriched chromosomes for developmental stage-regulated genes tend to harbor disproportionately large clusters of these genes. Further, although paralogous genes were overrepresented in gene clusters, their presence is not enough to explain two-thirds of the documented cases of whole-chromosome enrichment. The composition of the largest clusters often included paralogs from more than one multigene family as well as unrelated single-copy genes. Our results reveal intriguing patterns at the whole-chromosome scale in D. plexippus while shedding light on the interplay between gene expression and chromosome organization beyond Diptera and Hymenoptera.

Sex-biased gene content is associated with sex chromosome turnover in Danaini butterflies

Sex chromosomes play an outsized role in adaptation and speciation, and thus deserve particular attention in evolutionary genomics. In particular, fusions between sex chromosomes and autosomes can produce neo-sex chromosomes, which offer important insights into the evolutionary dynamics of sex chromosomes. Here, we investigate the evolutionary origin of the previously reported Danaus neo-sex chromosome within the tribe Danaini. We assembled and annotated genomes of Tirumala septentrionis (subtribe Danaina), Ideopsis similis (Amaurina), Idea leuconoe (Euploeina) and Lycorea halia (Itunina) and identified their Z-linked scaffolds. We found that the Danaus neo-sex chromosome resulting from the fusion between a Z chromosome and an autosome corresponding to the Melitaea cinxia chromosome (McChr) 21 arose in a common ancestor of Danaina, Amaurina and Euploina. We also identified two additional fusions as the W chromosome further fused with the synteny block McChr31 in I. similis and independent fusion occurred between ancestral Z chromosome and McChr12 in L. halia. We further tested a possible role of sexually antagonistic selection in sex chromosome turnover by analysing the genomic distribution of sex-biased genes in I. leuconoe and L. halia. The autosomes corresponding to McChr21 and McChr31 involved in the fusions are significantly enriched in female- and male-biased genes, respectively, which could have hypothetically facilitated fixation of the neo-sex chromosomes. This suggests a role of sexual antagonism in sex chromosome turnover in Lepidoptera. The neo-Z chromosomes of both I. leuconoe and L. halia appear fully compensated in somatic tissues, but the extent of dosage compensation for the ancestral Z varies across tissues and species.

Unraveling the complex evolutionary history of lepidopteran chromosomes through ancestral chromosome reconstruction and novel chromosome nomenclature

BACKGROUND

Lepidoptera is one of the most species-rich animal groups, with substantial karyotype variations among species due to chromosomal rearrangements. Knowledge of the evolutionary patterns of lepidopteran chromosomes still needs to be improved.

RESULTS

Here, we used chromosome-level genome assemblies of 185 lepidopteran insects to reconstruct an ancestral reference genome and proposed a new chromosome nomenclature. Thus, we renamed over 5000 extant chromosomes with this system, revealing the historical events of chromosomal rearrangements and their features. Additionally, our findings indicate that, compared with autosomes, the Z chromosome in Lepidoptera underwent a fast loss of conserved genes, rapid acquisition of lineage-specific genes, and a low rate of gene duplication. Moreover, we presented evidence that all available 67 W chromosomes originated from a common ancestor chromosome, with four neo-W chromosomes identified, including one generated by fusion with an autosome and three derived through horizontal gene transfer. We also detected nearly 4000 inter-chromosomal gene movement events. Notably, Geminin is transferred from the autosome to the Z chromosome. When located on the autosome, Geminin shows female-biased expression, but on the Z chromosome, it exhibits male-biased expression. This contributes to the sexual dimorphism of body size in silkworms.

CONCLUSIONS

Our study sheds light on the complex evolutionary history of lepidopteran chromosomes based on ancestral chromosome reconstruction and novel chromosome nomenclature.

Exploring chromosome evolution in 250 million year old groups of dragonflies and damselflies (Insecta:Odonata)

Using recently published chromosome-length genome assemblies of two damselfly species, Ischnura elegans and Platycnemis pennipes, and two dragonfly species, Pantala flavescens and Tanypteryx hageni, we demonstrate that the autosomes of Odonata have undergone few fission, fusion, or inversion events, despite 250 million years of separation. In the four genomes discussed here, our results show that all autosomes have a clear ortholog in the ancestral karyotype. Despite this clear chromosomal orthology, we demonstrate that different factors, including concentration of repeat dynamics, GC content, relative position on the chromosome, and the relative proportion of coding sequence all influence the density of syntenic blocks across chromosomes. However, these factors do not interact to influence synteny the same way in any two pairs of species, nor is any one factor retained in all four species. Furthermore, it was previously unknown whether the micro-chromosomes in Odonata are descended from one ancestral chromosome. Despite structural rearrangements, our evidence suggests that the micro-chromosomes in the sampled Odonata do indeed descend from an ancestral chromosome, and that the micro-chromosome in P. flavescens was lost through fusion with autosomes.

Chromosomal-level Genome Assembly of the Coffee Bee Hawk Moth Reveals the Evolution of Chromosomes and the Molecular Basis of Distinct Phenotypes

Cephonodes hylas, the coffee bee hawk moth is a hawk moth species with unique characteristics, such as larvae feeding on gardenia, overcoming the toxicity of its iridoid glycosides, diurnal adults, and transparent wings. Although C. hylas is a fascinating model for molecular biological research, genome sequence analysis-based genetic approaches to elucidate these peculiarities have not yet been undertaken. We successfully achieved de novo genome assembly at the chromosome level of C. hylas comparable to the Lepidoptera model organism, silkworm. Additionally, 16,854 protein-coding genes were annotated, and the constructed genome sequence and annotated genes were of the highest quality BUSCO completion compared to closely related species. Comparative genome analysis revealed the process of chromosomal evolution from the Bombycoidea ancestral (n = 31) genome and changes in turnover at the chromosome level associated with chromosomal fusion events, such as the rate of repetitive sequence insertion. These analyses were only possible because the genome was constructed at the chromosome level. Additionally, increased the nonsynonymous/synonymous rate (dN/dS) ratios were observed in multiple photoreceptor-related genes that were strongly associated with the acquisition of diurnal activity. Furthermore, tandemly duplicated expanded genes containing many digestive and other enzymes and larval midgut-specific expression were also confirmed. These genes may be involved in the metabolism of genipin, a toxin found in gardenias. Using the genome sequence of C. hylas determined at the chromosome level, we have successfully identified new insights into the chromosomal evolution of Bombycoidea, as well as the relationship between the genome sequence and its characteristic traits.

Chromosomal conservatism vs chromosomal megaevolution: enigma of karyotypic evolution in Lepidoptera

In the evolution of many organisms, periods of slow genome reorganization (= chromosomal conservatism) are interrupted by bursts of numerous chromosomal changes (= chromosomal megaevolution). Using comparative analysis of chromosome-level genome assemblies, we investigated these processes in blue butterflies (Lycaenidae). We demonstrate that the phase of chromosome number conservatism is characterized by the stability of most autosomes and dynamic evolution of the sex chromosome Z, resulting in multiple variants of NeoZ chromosomes due to autosome-sex chromosome fusions. In contrast during the phase of rapid chromosomal evolution, the explosive increase in chromosome number occurs mainly through simple chromosomal fissions. We show that chromosomal megaevolution is a highly non-random canalized process, and in two phylogenetically independent Lysandra lineages, the drastic parallel increase in number of fragmented chromosomes was achieved, at least partially, through reuse of the same ancestral chromosomal breakpoints. In species showing chromosome number doubling, we found no blocks of duplicated sequences or duplicated chromosomes, thus refuting the hypothesis of polyploidy. In the studied taxa, long blocks of interstitial telomere sequences (ITSs) consist of (TTAGG)_n arrays interspersed with telomere-specific retrotransposons. ITSs are sporadically present in rapidly evolving Lysandra karyotypes, but not in the species with ancestral chromosome number. Therefore, we hypothesize that the transposition of telomeric sequences may be triggers of the rapid chromosome number increase. Finally, we discuss the hypothetical genomic and population mechanisms of chromosomal megaevolution and argue that the disproportionally high evolutionary role of the Z sex chromosome can be additionally reinforced by sex chromosome-autosome fusions and Z-chromosome inversions.

The genetic basis of wing spots in Pieris canidia butterflies

Spots in pierid butterflies and eyespots in nymphalid butterflies are likely non-homologous wing colour pattern elements, yet they share a few features in common. Both develop black scales that depend on the function of the gene spalt, and both might have central signalling cells. This suggests that both pattern elements may be sharing common genetic circuitry. Hundreds of genes have already been associated with the development of nymphalid butterfly eyespot patterns, but the genetic basis of the simpler spot patterns on the wings of pierid butterflies has not been investigated. To facilitate studies of pierid wing patterns, we report a high-quality draft genome assembly for Pieris canidia, the Indian cabbage white. We then conducted transcriptomic analyses of pupal wing tissues sampled from the spot and non-spot regions of P. canidia at 3-6 h post-pupation. A total of 1352 genes were differentially regulated between wing tissues with and without the black spot, including spalt, Krüppel-like factor 10, genes from the Toll, Notch, TGF-β, and FGFR signalling pathways, and several genes involved in the melanin biosynthetic pathway. We identified 14 genes that are up-regulated in both pierid spots and nymphalid eyespots and propose that spots and eyespots share regulatory modules despite their likely independent origins.

Evidence for a single, ancient origin of a genus-wide alternative life history strategy

Understanding the evolutionary origins and factors maintaining alternative life history strategies (ALHS) within species is a major goal of evolutionary research. While alternative alleles causing discrete ALHS are expected to purge or fix over time, one-third of the ~90 species of Colias butterflies are polymorphic for a female-limited ALHS called Alba. Whether Alba arose once, evolved in parallel, or has been exchanged among taxa is currently unknown. Using comparative genome-wide association study (GWAS) and population genomic analyses, we placed the genetic basis of Alba in time-calibrated phylogenomic framework, revealing that Alba evolved once near the base of the genus and has been subsequently maintained via introgression and balancing selection. CRISPR-Cas9 mutagenesis was then used to verify a putative cis-regulatory region of Alba, which we identified using phylogenetic foot printing. We hypothesize that this cis-regulatory region acts as a modular enhancer for the induction of the Alba ALHS, which has likely facilitated its long evolutionary persistence.

Whole-Genome Analysis Reveals the Dynamic Evolution of Holocentric Chromosomes in Satyrine Butterflies

Butterfly chromosomes are holocentric, i.e., lacking a localized centromere. Potentially, this can lead to rapid karyotypic evolution through chromosome fissions and fusions, since fragmented chromosomes retain kinetic activity, while fused chromosomes are not dicentric. However, the actual mechanisms of butterfly genome evolution are poorly understood. Here, we analyzed chromosome-scale genome assemblies to identify structural rearrangements between karyotypes of satyrine butterfly species. For the species pair Erebia ligea-Maniola jurtina, sharing the ancestral diploid karyotype 2n = 56 + ZW, we demonstrate a high level of chromosomal macrosynteny and nine inversions separating these species. We show that the formation of a karyotype with a low number of chromosomes (2n = 36 + ZW) in Erebia aethiops was based on ten fusions, including one autosome-sex chromosome fusion, resulting in a neo-Z chromosome. We also detected inversions on the Z sex chromosome that were differentially fixed between the species. We conclude that chromosomal evolution is dynamic in the satyrines, even in the lineage that preserves the ancestral chromosome number. We hypothesize that the exceptional role of Z chromosomes in speciation may be further enhanced by inversions and sex chromosome-autosome fusions. We argue that not only fusions/fissions but also inversions are drivers of the holocentromere-mediated mode of chromosomal speciation.

Tempo and Mode of Genome Structure Evolution in Insects

The division of the genome into discrete chromosomes is a fundamental characteristic of eukaryotic life. Insect taxonomists' early adoption of cytogenetics has led to an incredible amount of data describing genome structure across insects. In this article, we synthesize data from thousands of species and use biologically realistic models to infer the tempo and mode of chromosome evolution among insect orders. Our results show that orders vary dramatically in the overall rate of chromosome number evolution (a proxy of genome structural stability) and the pattern of evolution (e.g., the balance between fusions and fissions). These findings have important implications for our understanding of likely modes of speciation and offer insight into the most informative clades for future genome sequencing.

Diversity, duplication, and genomic organization of homeobox genes in Lepidoptera

Homeobox genes encode transcription factors with essential roles in patterning and cell fate in developing animal embryos. Many homeobox genes, including Hox and NK genes, are arranged in gene clusters, a feature likely related to transcriptional control. Sparse taxon sampling and fragmentary genome assemblies mean that little is known about the dynamics of homeobox gene evolution across Lepidoptera or about how changes in homeobox gene number and organization relate to diversity in this large order of insects. Here we analyze an extensive data set of high-quality genomes to characterize the number and organization of all homeobox genes in 123 species of Lepidoptera from 23 taxonomic families. We find most Lepidoptera have around 100 homeobox loci, including an unusual Hox gene cluster in which the lab gene is repositioned and the ro gene is next to pb A topologically associating domain spans much of the gene cluster, suggesting deep regulatory conservation of the Hox cluster arrangement in this insect order. Most Lepidoptera have four Shx genes, divergent zen-derived loci, but these loci underwent dramatic duplication in several lineages, with some moths having over 165 homeobox loci in the Hox gene cluster; this expansion is associated with local LINE element density. In contrast, the NK gene cluster content is more stable, although there are differences in organization compared with other insects, as well as major rearrangements within butterflies. Our analysis represents the first description of homeobox gene content across the order Lepidoptera, exemplifying the potential of newly generated genome assemblies for understanding genome and gene family evolution.

High-density linkage maps and chromosome level genome assemblies unveil direction and frequency of extensive structural rearrangements in wood white butterflies (Leptidea spp.)

Karyotypes are generally conserved between closely related species and large chromosome rearrangements typically have negative fitness consequences in heterozygotes, potentially driving speciation. In the order Lepidoptera, most investigated species have the ancestral karyotype and gene synteny is often conserved across deep divergence, although examples of extensive genome reshuffling have recently been demonstrated. The genus Leptidea has an unusual level of chromosome variation and rearranged sex chromosomes, but the extent of restructuring across the rest of the genome is so far unknown. To explore the genomes of the wood white (Leptidea) species complex, we generated eight genome assemblies using a combination of 10X linked reads and HiC data, and improved them using linkage maps for two populations of the common wood white (L. sinapis) with distinct karyotypes. Synteny analysis revealed an extensive amount of rearrangements, both compared to the ancestral karyotype and between the Leptidea species, where only one of the three Z chromosomes was conserved across all comparisons. Most restructuring was explained by fissions and fusions, while translocations appear relatively rare. We further detected several examples of segregating rearrangement polymorphisms supporting a highly dynamic genome evolution in this clade. Fusion breakpoints were enriched for LINEs and LTR elements, which suggests that ectopic recombination might be an important driver in the formation of new chromosomes. Our results show that chromosome count alone may conceal the extent of genome restructuring and we propose that the amount of genome evolution in Lepidoptera might still be underestimated due to lack of taxonomic sampling.

Testing hypotheses of a coevolutionary key innovation reveals a complex suite of traits involved in defusing the mustard oil bomb

Coevolutionary interactions are responsible for much of the Earth's biodiversity, with key innovations driving speciation bursts on both sides of the interaction. One persistent question is whether macroevolutionary traits identified as key innovations accurately predict functional performance and selection dynamics within species, as this necessitates characterizing their function, investigating their fitness consequences, and exploring the selection dynamics acting upon them. Here, we used CRISPR-Cas9 mediating nonhomologous end joining (NHEJ) in the butterfly species Pieris brassicae to knock out and directly assess the function and fitness impacts of nitrile specifier protein (NSP) and major allergen (MA). These are two closely related genes that facilitate glucosinolate (GSL) detoxification capacity, which is a key innovation in mustard feeding Pierinae butterflies. We find NSP and MA are both required for survival on plants containing GSLs, with expression differences arising in response to variable GSL profiles, concordant with detoxification performance. Importantly, this concordance was only observed when using natural host plants, likely reflecting the complexity of how these enzymes interact with natural plant variation in GSLs and myrosinases. Finally, signatures of positive selection for NSP and MA were detected across Pieris species, consistent with these genes' importance in recent coevolutionary interactions. Thus, the war between these butterflies and their host plants involves more than the mere presence of chemical defenses and detoxification mechanisms, as their regulation and activation represent key components of complex interactions. We find that inclusion of these dynamics, in ecologically relevant assays, is necessary for coevolutionary insights in this system and likely others.

A step forward in the genome characterization of the sugarcane borer, Diatraea saccharalis: karyotype analysis, sex chromosome system and repetitive DNAs through a cytogenomic approach

Moths of the family Crambidae include a number of pests that cause economic losses to agricultural crops. Despite their economic importance, little is known about their genome architecture and chromosome evolution. Here, we characterized the chromosomes and repetitive DNA of the sugarcane borer Diatraea saccharalis using a combination of low-pass genome sequencing, bioinformatics, and cytogenetic methods, focusing on the sex chromosomes. Diploid chromosome numbers differed between the sexes, i.e., 2n = 33 in females and 2n = 34 in males. This difference was caused by the occurrence of a WZ₁Z₂ trivalent in female meiosis, indicating a multiple sex-chromosome system WZ₁Z₂/Z₁Z₁Z₂Z₂. A strong interstitial telomeric signal was observed on the W chromosome, indicating a fusion of the ancestral W chromosome with an autosome. Among repetitive DNAs, transposable elements (TEs) accounted for 39.18% (males) to 41.35% (females), while satDNAs accounted for only 0.214% (males) and 0.215% (females) of the genome. FISH mapping revealed different chromosomal organization of satDNAs, such as single localized clusters, spread repeats, and non-clustered repeats. Two TEs mapped by FISH were scattered. Although we found a slight enrichment of some satDNAs in the female genome, they were not differentially enriched on the W chromosome. However, we found enriched FISH signals for TEs on the W chromosome, suggesting their involvement in W chromosome degeneration and differentiation. These data shed light on karyotype and repetitive DNA dynamics due to multiple chromosome fusions in D. saccharalis, contribute to the understanding of genome structure in Lepidoptera and are important for future genomic studies.

A high-quality Buxus austro-yunnanensis (Buxales) genome provides new insights into karyotype evolution in early eudicots

Background

Eudicots are the most diverse group of flowering plants that compromise five well-defined lineages: core eudicots, Ranunculales, Proteales, Trochodendrales, and Buxales. However, the phylogenetic relationships between these five lineages and their chromosomal evolutions remain unclear, and a lack of high-quality genome analyses for Buxales has hindered many efforts to address this knowledge gap.

Results

Here, we present a high-quality chromosome-level genome of Buxus austro-yunnanensis (Buxales). Our phylogenomic analyses revealed that Buxales and Trochodendrales are genetically similar and classified as sisters. Additionally, both are sisters to the core eudicots, while Ranunculales was found to be the first lineage to diverge from these groups. Incomplete lineage sorting and hybridization were identified as the main contributors to phylogenetic discordance (34.33%) between the lineages. In fact, B. austro-yunnanensis underwent only one whole-genome duplication event, and collinear gene phylogeny analyses suggested that separate independent polyploidizations occurred in the five eudicot lineages. Using representative genomes from these five lineages, we reconstructed the ancestral eudicot karyotype (AEK) and generated a nearly gapless karyotype projection for each eudicot species. Within core eudicots, we recovered one common chromosome fusion event in asterids and malvids, respectively. Further, we also found that the previously reported fused AEKs in Aquilegia (Ranunculales) and Vitis (core eudicots) have different fusion positions, which indicates that these two species have different karyotype evolution histories.

Conclusions

Based on our phylogenomic and karyotype evolution analyses, we revealed the likely relationships and evolutionary histories of early eudicots. Ultimately, our study expands genomic resources for early-diverging eudicots.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01420-1.

A genome for Edith's checkerspot butterfly: an insect with complex host-adaptive suites and rapid evolutionary responses to environmental changes

Insects have been key players in the assessments of biodiversity impacts of anthropogenically driven environmental change, including the evolutionary and ecological impacts of climate change. Populations of Edith’s Checkerspot Butterfly (Euphydryas editha) adapt rapidly to diverse environmental conditions, with numerous high-impact studies documenting these dynamics over several decades. However, studies of the underlying genetic bases of these responses have been hampered by missing genomic resources, limiting the ability to connect genomic responses to environmental change. Using a combination of Oxford Nanopore long reads, haplotype merging, HiC scaffolding followed by Illumina polishing, we generated a highly contiguous and complete assembly (contigs n = 142, N50 = 21.2 Mb, total length = 607.8 Mb; BUSCOs n = 5,286, single copy complete = 97.8%, duplicated = 0.9%, fragmented = 0.3%, missing = 1.0%). A total of 98% of the assembled genome was placed into 31 chromosomes, which displayed large-scale synteny with other well-characterized lepidopteran genomes. The E. editha genome, annotation, and functional descriptions now fill a missing gap for one of the leading field-based ecological model systems in North America.

Celebrating Mendel, McClintock, and Darlington: On end-to-end chromosome fusions and nested chromosome fusions

The evolution of eukaryotic genomes is accompanied by fluctuations in chromosome number, reflecting cycles of chromosome number increase (polyploidy and centric fissions) and decrease (chromosome fusions). Although all chromosome fusions result from DNA recombination between two or more nonhomologous chromosomes, several mechanisms of descending dysploidy are exploited by eukaryotes to reduce their chromosome number. Genome sequencing and comparative genomics have accelerated the identification of inter-genome chromosome collinearity and gross chromosomal rearrangements and have shown that end-to-end chromosome fusions (EEFs) and nested chromosome fusions (NCFs) may have played a more important role in the evolution of eukaryotic karyotypes than previously thought. The present review aims to summarize the limited knowledge on the origin, frequency, and evolutionary implications of EEF and NCF events in eukaryotes and especially in land plants. The interactions between nonhomologous chromosomes in interphase nuclei and chromosome (mis)pairing during meiosis are examined for their potential importance in the origin of EEFs and NCFs. The remaining open questions that need to be addressed are discussed.

A chromosome-level genome assembly of Paracymoriza distinctalis (Lepidoptera: Crambidae: Acentropinae)

Paracymoriza distinctalis is a semiaquatic lepidopteran insect, which is of great value for studying the differentiation of the Pyraloidea superfamily. However, the understanding of heredity, evolution, and functional genomics of P. distinctalis are limited by few genome-wide resources. Here, we applied PacBio sequencing and the chromosome capture technique to assemble the first P. distinctalis genome from a single female individual. The genome size is 1.2 Gb with 32 chromosomes and the N50 is 38.91 Mb. Approximately 576.37 Mb, accounting for 48.93% of the genome, was identified as repeats. The genome comprises 39,003 protein-coding genes, 66.56% of which were functionally annotated. Comparative genomics analysis suggested that the common ancestor of P. distinctalis and Chilo suppressalis lived ~83.5 million years ago. This chromosome-level genome assembly work is not only conducive to the understanding of P. distinctalis, but also may promote the study of the genomes of other lepidopteran species.

Dissecting cricket genomes for the advancement of entomology and entomophagy

Significant advances in biophysical methods such as next-generation sequencing technologies have now opened the way to conduct evolutionary and applied research based on the genomic information of greatly diverse insects. Crickets belonging to Orthoptera (Insecta: Polyneoptera), one of the most flourishing groups of insects, have contributed to the development of multiple scientific fields including developmental biology and neuroscience and have been attractive targets in evolutionary ecology for their diverse ecological niches. In addition, crickets have recently gained recognition as food and feed. However, the genomic information underlying their biological basis and application research toward breeding is currently underrepresented. In this review, we summarize the progress of genomics of crickets. First, we outline the phylogenetic position of crickets in insects and then introduce recent studies on cricket genomics and transcriptomics in a variety of fields. Furthermore, we present findings from our analysis of polyneopteran genomes, with a particular focus on their large genome sizes, chromosome number, and repetitive sequences. Finally, how the cricket genome can be beneficial to the food industry is discussed. This review is expected to enhance greater recognition of how important the cricket genomes are to the multiple biological fields and how basic research based on cricket genome information can contribute to tackling global food security.

Improved chromosome-level genome assembly of the Glanville fritillary butterfly (Melitaea cinxia) integrating Pacific Biosciences long reads and a high-density linkage map

BACKGROUND

The Glanville fritillary (Melitaea cinxia) butterfly is a model system for metapopulation dynamics research in fragmented landscapes. Here, we provide a chromosome-level assembly of the butterfly's genome produced from Pacific Biosciences sequencing of a pool of males, combined with a linkage map from population crosses.

RESULTS

The final assembly size of 484 Mb is an increase of 94 Mb on the previously published genome. Estimation of the completeness of the genome with BUSCO indicates that the genome contains 92-94% of the BUSCO genes in complete and single copies. We predicted 14,810 genes using the MAKER pipeline and manually curated 1,232 of these gene models.

CONCLUSIONS

The genome and its annotated gene models are a valuable resource for future comparative genomics, molecular biology, transcriptome, and genetics studies on this species.

The draft genome of the Asian corn borer yields insights into ecological adaptation of a devastating maize pest

The Asian corn borer (ACB) is the most devastating pest on maize in the western Pacific region of Asia. Despite broad interests in insecticide resistance, seasonal adaptation, and larval color mimicry regarding the ACB system, lacking of reference genomic information and a powerful gene editing approach have hindered the in-depth studies of these aspects. Here we present a 455.7 Mb draft genome of ACB with 98.4% completeness. Comparative genomics analysis showed an evident expansion in gene families of gustatory receptors (105), which is related to polyphagous characteristics. Based on the comparative transcriptome analysis of resistant and susceptible ACB against Bt Cry1Ab toxin, we identified 26 genes related to Cry1Ab resistance. Additionally, transcriptomics of insects exposed to conditions of low temperature and diapause (LT) vs. room temperature and diapause (RT) provided insights into the genetic mechanisms of cold adaptation. We also successfully developed an efficient CRISPR/Cas9-based genome editing system and applied it to explore the role of color pattern genes in the ecological adaptation of ACB. Taken together, our study provides a fully annotated high-quality reference genome and efficient gene editing system to realize the potential of ACB as a study system to address important biological questions such as insecticide resistance, seasonal adaptation, and coloration. These valuable genomic resources will also benefit the development of novel strategies for maize pest management.

The genome sequences of the male and female green-veined white, Pieris napi (Linnaeus, 1758)

We present genome assemblies from a male and female Pieris napi (the green-veined white; Arthropoda; Insecta; Lepidoptera; Pieridae). The genome sequences of the male and female are 320 and 319 megabases in span, respectively. The majority of the assembly (99.79% of the male assembly, 99.88% of the female) is scaffolded into 24 autosomal pseudomolecules, with the Z sex chromosome assembled for the male and Z and W chromosomes assembled for the female. Gene annotation of the male assembly on Ensembl has identified 13,221 protein coding genes.

Comparative analysis of chromosome numbers and sex chromosome systems in Paraneoptera (Insecta)

This article is part (the 4th article) of the themed issue (a monograph) "Aberrant cytogenetic and reproductive patterns in the evolution of Paraneoptera". The purpose of this article is to consider chromosome structure and evolution, chromosome numbers and sex chromosome systems, which all together constitute the chromosomal basis of reproduction and are essential for reproductive success. We are based on our own observations and literature data available for all major lineages of Paraneoptera including Zoraptera (angel insects), Copeognatha (=Psocoptera; bark lice), Parasita (=Phthiraptera s. str; true lice), Thysanoptera (thrips), Homoptera (scale insects, aphids, jumping plant-lice, whiteflies, and true hoppers), Heteroptera (true bugs), and Coleorrhyncha (moss bugs). Terminology, nomenclature, classification, and the study methods are given in the first paper of the issue (Gavrilov-Zimin et al. 2021).

The Easter Egg Weevil (Pachyrhynchus) genome reveals syntenic patterns in Coleoptera across 200 million years of evolution

Patterns of genomic architecture across insects remain largely undocumented or decoupled from a broader phylogenetic context. For instance, it is unknown whether translocation rates differ between insect orders. We address broad scale patterns of genome architecture across Insecta by examining synteny in a phylogenetic framework from open-source insect genomes. To accomplish this, we add a chromosome level genome to a crucial lineage, Coleoptera. Our assembly of the Pachyrhynchus sulphureomaculatus genome is the first chromosome scale genome for the hyperdiverse Phytophaga lineage and currently the largest insect genome assembled to this scale. The genome is significantly larger than those of other weevils, and this increase in size is caused by repetitive elements. Our results also indicate that, among beetles, there are instances of long-lasting (>200 Ma) localization of genes to a particular chromosome with few translocation events. While some chromosomes have a paucity of translocations, intra-chromosomal synteny was almost absent, with gene order thoroughly shuffled along a chromosome. This large amount of reshuffling within chromosomes with few inter-chromosomal events contrasts with patterns seen in mammals in which the chromosomes tend to exchange larger blocks of material more readily. To place our findings in an evolutionary context, we compared syntenic patterns across Insecta in a phylogenetic framework. For the first time, we find that synteny decays at an exponential rate relative to phylogenetic distance. Additionally, there are significant differences in decay rates between insect orders, this pattern was not driven by Lepidoptera alone which has a substantially different rate.

A telomere-to-telomere assembly of Oscheius tipulae and the evolution of rhabditid nematode chromosomes

Eukaryotic chromosomes have phylogenetic persistence. In many taxa, each chromosome has a single functional centromere with essential roles in spindle attachment and segregation. Fusion and fission can generate chromosomes with no or multiple centromeres, leading to genome instability. Groups with holocentric chromosomes (where centromeric function is distributed along each chromosome) might be expected to show karyotypic instability. This is generally not the case, and in Caenorhabditis elegans, it has been proposed that the role of maintenance of a stable karyotype has been transferred to the meiotic pairing centers, which are found at one end of each chromosome. Here, we explore the phylogenetic stability of nematode chromosomes using a new telomere-to-telomere assembly of the rhabditine nematode Oscheius tipulae generated from nanopore long reads. The 60-Mb O. tipulae genome is resolved into six chromosomal molecules. We find the evidence of specific chromatin diminution at all telomeres. Comparing this chromosomal O. tipulae assembly with chromosomal assemblies of diverse rhabditid nematodes, we identify seven ancestral chromosomal elements (Nigon elements) and present a model for the evolution of nematode chromosomes through rearrangement and fusion of these elements. We identify frequent fusion events involving NigonX, the element associated with the rhabditid X chromosome, and thus sex chromosome-associated gene sets differ markedly between species. Despite the karyotypic stability, gene order within chromosomes defined by Nigon elements is not conserved. Our model for nematode chromosome evolution provides a platform for investigation of the tensions between local genome rearrangement and karyotypic evolution in generating extant genome architectures.

The genome assembly and annotation of the Apollo butterfly Parnassius apollo, a flagship species for conservation biology

Conservation genomics has made dramatic improvements over the past decade, leveraging the power of genomes to infer diverse parameters central to conservation management questions. However, much of this effort has focused upon vertebrate species, despite insects providing similar flagship status with the added benefit of smaller genomes, shorter generation times and extensive historical collections in museums. Here we present the genome of the Apollo butterfly (Parnassius apollo, Papilionidae), an iconic endangered butterfly, which like many species in this genus, needs conservation genomic attention yet lacks a genome. Using 68.7 Gb of long-read data (N50 = 15.2 kb) we assembled a 1.4 Gb genome for the Apollo butterfly, making this the largest sequenced Lepidopteran genome to date. The assembly was highly contiguous (N50 = 7.1 Mb) and complete (97% of Lepidopteran BUSCOs were single-copy and complete) and consisted of 1,707 contigs. Using RNAseq data and Arthropoda proteins, we annotated 28.3K genes. Alignment with the closest-related chromosome-level assembly, Papilio bianor, reveals a highly conserved chromosomal organization, albeit genome size is highly expanded in the Apollo butterfly, due primarily to a dramatic increase in repetitive element content. Using this alignment for superscaffolding places the P. apollo genome in to 31 chromosomal scaffolds, and together with our functional annotation, provides an essential resource for advancing conservation genomics in a flagship species for insect conservation.

Large-scale comparative analysis of cytogenetic markers across Lepidoptera

Fluorescence in situ hybridization (FISH) allows identification of particular chromosomes and their rearrangements. Using FISH with signal enhancement via antibody amplification and enzymatically catalysed reporter deposition, we evaluated applicability of universal cytogenetic markers, namely 18S and 5S rDNA genes, U1 and U2 snRNA genes, and histone H3 genes, in the study of the karyotype evolution in moths and butterflies. Major rDNA underwent rather erratic evolution, which does not always reflect chromosomal changes. In contrast, the hybridization pattern of histone H3 genes was well conserved, reflecting the stable organisation of lepidopteran genomes. Unlike 5S rDNA and U1 and U2 snRNA genes which we failed to detect, except for 5S rDNA in a few representatives of early diverging lepidopteran lineages. To explain the negative FISH results, we used quantitative PCR and Southern hybridization to estimate the copy number and organization of the studied genes in selected species. The results suggested that their detection was hampered by long spacers between the genes and/or their scattered distribution. Our results question homology of 5S rDNA and U1 and U2 snRNA loci in comparative studies. We recommend the use of histone H3 in studies of karyotype evolution.

Zebrafish Chromosome 14 Gene Differential Expression in the fmr1 h u2787 Model of Fragile X Syndrome

Zebrafish represent a valuable model for investigating the molecular and cellular basis of Fragile X syndrome (FXS). Reduced expression of the zebrafish FMR1 orthologous gene, fmr1, causes developmental and behavioural phenotypes related to FXS. Zebrafish homozygous for the hu2787 non-sense mutation allele of fmr1 are widely used to model FXS, although FXS-relevant phenotypes seen from morpholino antisense oligonucleotide (morpholino) suppression of fmr1 transcript translation were not observed when hu2787 was first described. The subsequent discovery of transcriptional adaptation (a form of genetic compensation), whereby mutations causing non-sense-mediated decay of transcripts can drive compensatory upregulation of homologous transcripts independent of protein feedback loops, suggested an explanation for the differences reported. We examined the whole-embryo transcriptome effects of homozygosity for fmr1 h u2787 at 2 days post fertilisation. We observed statistically significant changes in expression of a number of gene transcripts, but none from genes showing sequence homology to fmr1. Enrichment testing of differentially expressed genes implied effects on lysosome function and glycosphingolipid biosynthesis. The majority of the differentially expressed genes are located, like fmr1, on Chromosome 14. Quantitative PCR tests did not support that this was artefactual due to changes in relative chromosome abundance. Enrichment testing of the "leading edge" differentially expressed genes from Chromosome 14 revealed that their co-location on this chromosome may be associated with roles in brain development and function. The differential expression of functionally related genes due to mutation of fmr1, and located on the same chromosome as fmr1, is consistent with R.A. Fisher's assertion that the selective advantage of co-segregation of particular combinations of alleles of genes will favour, during evolution, chromosomal rearrangements that place them in linkage disequilibrium on the same chromosome. However, we cannot exclude that the apparent differential expression of genes on Chromosome 14 genes was, (if only in part), caused by differences between the expression of alleles of genes unrelated to the effects of the fmr1 h u2787 mutation and made manifest due to the limited, but non-zero, allelic diversity between the genotypes compared.

Evolution of the codling moth pheromone via an ancient gene duplication

BACKGROUND

Defining the origin of genetic novelty is central to our understanding of the evolution of novel traits. Diversification among fatty acid desaturase (FAD) genes has played a fundamental role in the introduction of structural variation in fatty acyl derivatives. Because of its central role in generating diversity in insect semiochemicals, the FAD gene family has become a model to study how gene family expansions can contribute to the evolution of lineage-specific innovations. Here we used the codling moth (Cydia pomonella) as a study system to decipher the proximate mechanism underlying the production of the ∆8∆10 signature structure of olethreutine moths. Biosynthesis of the codling moth sex pheromone, (E8,E10)-dodecadienol (codlemone), involves two consecutive desaturation steps, the first of which is unusual in that it generates an E9 unsaturation. The second step is also atypical: it generates a conjugated diene system from the E9 monoene C12 intermediate via 1,4-desaturation.

RESULTS

Here we describe the characterization of the FAD gene acting in codlemone biosynthesis. We identify 27 FAD genes corresponding to the various functional classes identified in insects and Lepidoptera. These genes are distributed across the C. pomonella genome in tandem arrays or isolated genes, indicating that the FAD repertoire consists of both ancient and recent duplications and expansions. Using transcriptomics, we show large divergence in expression domains: some genes appear ubiquitously expressed across tissue and developmental stages; others appear more restricted in their expression pattern. Functional assays using heterologous expression systems reveal that one gene, Cpo_CPRQ, which is prominently and exclusively expressed in the female pheromone gland, encodes an FAD that possesses both E9 and ∆8∆10 desaturation activities. Phylogenetically, Cpo_CPRQ clusters within the Lepidoptera-specific ∆10/∆11 clade of FADs, a classic reservoir of unusual desaturase activities in moths.

CONCLUSIONS

Our integrative approach shows that the evolution of the signature pheromone structure of olethreutine moths relied on a gene belonging to an ancient gene expansion. Members of other expanded FAD subfamilies do not appear to play a role in chemical communication. This advises for caution when postulating the consequences of lineage-specific expansions based on genomics alone.

The Genome of the Margined White Butterfly (Pieris macdunnoughii): Sex Chromosome Insights and the Power of Polishing with PoolSeq Data

We report a chromosome-level assembly for Pieris macdunnoughii, a North American butterfly whose involvement in an evolutionary trap imposed by an invasive Eurasian mustard has made it an emerging model system for studying maladaptation in plant-insect interactions. Assembled using nearly 100× coverage of Oxford Nanopore long reads, the contig-level assembly comprised 106 contigs totaling 316,549,294 bases, with an N50 of 5.2 Mb. We polished the assembly with PoolSeq Illumina short-read data, demonstrating for the first time the comparable performance of individual and pooled short reads as polishing data sets. Extensive synteny between the reported contig-level assembly and a published, chromosome-level assembly of the European butterfly Pieris napi allowed us to generate a pseudochromosomal assembly of 47 contigs, placing 91.1% of our 317 Mb genome into a chromosomal framework. Additionally, we found support for a Z chromosome arrangement in P. napi, showing that the fusion event leading to this rearrangement predates the split between European and North American lineages of Pieris butterflies. This genome assembly and its functional annotation lay the groundwork for future research into the genetic basis of adaptive and maladaptive egg-laying behavior by P. macdunnoughii, contributing to our understanding of the susceptibility and responses of insects to evolutionary traps.

Chromosome-Scale Genome Assemblies of Aphids Reveal Extensively Rearranged Autosomes and Long-Term Conservation of the X Chromosome

Chromosome rearrangements are arguably the most dramatic type of mutations, often leading to rapid evolution and speciation. However, chromosome dynamics have only been studied at the sequence level in a small number of model systems. In insects, Diptera and Lepidoptera have conserved genome structure at the scale of whole chromosomes or chromosome arms. Whether this reflects the diversity of insect genome evolution is questionable given that many species exhibit rapid karyotype evolution. Here, we investigate chromosome evolution in aphids-an important group of hemipteran plant pests-using newly generated chromosome-scale genome assemblies of the green peach aphid (Myzus persicae) and the pea aphid (Acyrthosiphon pisum), and a previously published assembly of the corn-leaf aphid (Rhopalosiphum maidis). We find that aphid autosomes have undergone dramatic reorganization over the last 30 My, to the extent that chromosome homology cannot be determined between aphids from the tribes Macrosiphini (Myzus persicae and Acyrthosiphon pisum) and Aphidini (Rhopalosiphum maidis). In contrast, gene content of the aphid sex (X) chromosome remained unchanged despite rapid sequence evolution, low gene expression, and high transposable element load. To test whether rapid evolution of genome structure is a hallmark of Hemiptera, we compared our aphid assemblies with chromosome-scale assemblies of two blood-feeding Hemiptera (Rhodnius prolixus and Triatoma rubrofasciata). Despite being more diverged, the blood-feeding hemipterans have conserved synteny. The exceptional rate of structural evolution of aphid autosomes renders them an important emerging model system for studying the role of large-scale genome rearrangements in evolution.

A Population Genomic Investigation of Immune Cell Diversity and Phagocytic Capacity in a Butterfly

Insects rely on their innate immune system to successfully mediate complex interactions with their internal microbiota, as well as the microbes present in the environment. Given the variation in microbes across habitats, the challenges to respond to them are likely to result in local adaptations in the immune system. Here we focus upon phagocytosis, a mechanism by which pathogens and foreign particles are engulfed in order to be contained, killed, and processed. We investigated the phenotypic and genetic variation related to phagocytosis in two allopatric populations of the butterfly Pieris napi. Populations were found to differ in their hemocyte composition and overall phagocytic capability, driven by the increased phagocytic propensity of each cell type. Yet, genes annotated to phagocytosis showed no large genomic signal of divergence. However, a gene set enrichment analysis on significantly divergent genes identified loci involved in glutamine metabolism, which recently have been linked to immune cell differentiation in mammals. Together these results suggest that heritable variation in phagocytic capacity arises via a quantitative trait architecture with variation in genes affecting the activation and/or differentiation of phagocytic cells, suggesting them as potential candidate genes underlying these phenotypic differences.

Physiological Tradeoffs of Immune Response Differs by Infection Type in Pieris napi

Understanding the tradeoffs that result from successful infection responses is central to understanding how life histories evolve. Gaining such insights, however, can be challenging, as they may be pathogen specific and confounded with experimental design. Here, we investigated whether infection from gram positive or negative bacteria results in different physiological tradeoffs, and whether these infections impact life history later in life (post-diapause development), in the butterfly Pieris napi. During the first 24 h after infection (3, 6, 12, and 24 h), after removing effects due to injection, larvae infected with Micrococcus luteus showed a strong suppression of all non-immunity related processes while several types of immune responses were upregulated. In contrast, this tradeoff between homeostasis and immune response was much less pronounced in Escherichia coli infections. These differences were also visible long after infection, via weight loss and slower development, as well as an increased mortality at higher infection levels during later stages of development. Individuals infected with M. luteus, compared to E. coli, had a higher mortality rate, and a lower pupal weight, developmental rate and adult weight. Further, males exhibited a more negative impact of infection than females. Thus, immune responses come at a cost even when the initial infection has been overcome, and these costs are likely to affect later life history parameters with fitness consequences.

Meiosis Progression and Recombination in Holocentric Plants: What Is Known?

Differently from the common monocentric organization of eukaryotic chromosomes, the so-called holocentric chromosomes present many centromeric regions along their length. This chromosomal organization can be found in animal and plant lineages, whose distribution suggests that it has evolved independently several times. Holocentric chromosomes present an advantage: even broken chromosome parts can be correctly segregated upon cell division. However, the evolution of holocentricity brought about consequences to nuclear processes and several adaptations are necessary to cope with this new organization. Centromeres of monocentric chromosomes are involved in a two-step cohesion release during meiosis. To deal with that holocentric lineages developed different adaptations, like the chromosome remodeling strategy in Caenorhabditis elegans or the inverted meiosis in plants. Furthermore, the frequency of recombination at or around centromeres is normally very low and the presence of centromeric regions throughout the entire length of the chromosomes could potentially pose a problem for recombination in holocentric organisms. However, meiotic recombination happens, with exceptions, in those lineages in spite of their holocentric organization suggesting that the role of centromere as recombination suppressor might be altered in these lineages. Most of the available information about adaptations to meiosis in holocentric organisms is derived from the animal model C. elegans. As holocentricity evolved independently in different lineages, adaptations observed in C. elegans probably do not apply to other lineages and very limited research is available for holocentric plants. Currently, we still lack a holocentric model for plants, but good candidates may be found among Cyperaceae, a large angiosperm family. Besides holocentricity, chiasmatic and achiasmatic inverted meiosis are found in the family. Here, we introduce the main concepts of meiotic constraints and adaptations with special focus in meiosis progression and recombination in holocentric plants. Finally, we present the main challenges and perspectives for future research in the field of chromosome biology and meiosis in holocentric plants.

A region of the sex chromosome associated with population differences in diapause induction contains highly divergent alleles at clock genes

Developmental plasticity describes the capacity of individuals with the same genotype to induce permanent change in a phenotype depending on a specific external input. One well‐studied example of adaptive developmental plasticity is the induction of facultative diapause in insects. Studies investigating the inheritance of diapause induction have suggested diverse genetic origins. However, only few studies have performed genome‐wide scans to identify genes affecting the induction decision. Here we compare two populations of the butterfly Pieris napi that differ in the propensity to enter diapause, and despite showing a low genome‐wide divergence, we identify a few genomic regions that show high divergence between populations. We then identified a single genomic region associated with diapause induction by genotyping diapausing and directly developing siblings from backcrosses of these populations. This region is located on the Z chromosome and contained three circadian clock genes, cycle, clock, and period. Additionally, period harbored the largest number of SNPs showing complete fixation between populations. We conclude that the heritable basis of between‐population variation in the plasticity that determines diapause induction resides on the Z chromosome, with the period gene being the prime candidate for the genetic basis of adaptive plasticity.

Chromosome number evolves at equal rates in holocentric and monocentric clades

Despite the fundamental role of centromeres two different types are observed across plants and animals. Monocentric chromosomes possess a single region that function as the centromere while in holocentric chromosomes centromere activity is spread across the entire chromosome. Proper segregation may fail in species with monocentric chromosomes after a fusion or fission, which may lead to chromosomes with no centromere or multiple centromeres. In contrast, species with holocentric chromosomes should still be able to safely segregate chromosomes after fusion or fission. This along with the observation of high chromosome number in some holocentric clades has led to the hypothesis that holocentricity leads to higher rates of chromosome number evolution. To test for differences in rates of chromosome number evolution between these systems, we analyzed data from 4,393 species of insects in a phylogenetic framework. We found that insect orders exhibit striking differences in rates of fissions, fusions, and polyploidy. However, across all insects we found no evidence that holocentric clades have higher rates of fissions, fusions, or polyploidy than monocentric clades. Our results suggest that holocentricity alone does not lead to higher rates of chromosome number changes. Instead, we suggest that other co-evolving traits must explain striking differences between clades.

Holocentric chromosomes

Holocentric chromosomes possess multiple kinetochores along their length rather than the single centromere typical of other chromosomes [1]. They have been described for the first time in cytogenetic experiments dating from 1935 and, since this first observation, the term holocentric chromosome has referred to chromosomes that: i. lack the primary constriction corresponding to centromere observed in monocentric chromosomes [2]; ii. possess multiple kinetochores dispersed along the chromosomal axis so that microtubules bind to chromosomes along their entire length and move broadside to the pole from the metaphase plate [3]. These chromosomes are also termed holokinetic, because, during cell division, chromatids move apart in parallel and do not form the classical V-shaped figures typical of monocentric chromosomes [4–6]. Holocentric chromosomes evolved several times during both animal and plant evolution and are currently reported in about eight hundred diverse species, including plants, insects, arachnids and nematodes [7,8]. As a consequence of their diffuse kinetochores, holocentric chromosomes may stabilize chromosomal fragments favouring karyotype rearrangements [9,10]. However, holocentric chromosome may also present limitations to crossing over causing a restriction of the number of chiasma in bivalents [11] and may cause a restructuring of meiotic divisions resulting in an inverted meiosis [12].

Speciation through chromosomal fusion and fission in Lepidoptera

Changes in chromosome numbers may strongly affect reproductive barriers, because individuals heterozygous for distinct karyotypes are typically expected to be at least partially sterile or to show reduced recombination. Therefore, several classic speciation models are based on chromosomal changes. One import mechanism generating variation in chromosome numbers is fusion and fission of existing chromosomes, which is particularly likely in species with holocentric chromosomes, i.e. chromosomes that lack a single centromere. Holocentric chromosomes evolved repeatedly across the tree of life, including in Lepidoptera. Although changes in chromosome numbers are hypothesized to be an important driver of the spectacular diversification of Lepidoptera, comparative studies across the order are lacking. We performed the first comprehensive literature survey of karyotypes for Lepidoptera species since the 1970s and tested if, and how, chromosomal variation might affect speciation. Even though a meta-analysis of karyological differences between closely related taxa did not reveal an effect on the degree of reproductive isolation, phylogenetic diversification rate analyses across the 16 best-covered genera indicated a strong, positive association of rates of chromosome number evolution and speciation. These findings suggest a macroevolutionary impact of varying chromosome numbers in Lepidoptera and likely apply to other taxonomic groups, especially to those with holocentric chromosomes. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Evolution of multiple sex-chromosomes associated with dynamic genome reshuffling in Leptidea wood-white butterflies

Sex-chromosome systems tend to be highly conserved and knowledge about their evolution typically comes from macroevolutionary inference. Rapidly evolving complex sex-chromosome systems represent a rare opportunity to study the mechanisms of sex-chromosome evolution at unprecedented resolution. Three cryptic species of wood-white butterflies—Leptidea juvernica, L. sinapis and L. reali—have each a unique set of multiple sex-chromosomes with 3–4 W and 3–4 Z chromosomes. Using a transcriptome-based microarray for comparative genomic hybridisation (CGH) and a library of bacterial artificial chromosome (BAC) clones, both developed in L. juvernica, we identified Z-linked Leptidea orthologs of Bombyx mori genes and mapped them by fluorescence in situ hybridisation (FISH) with BAC probes on multiple Z chromosomes. In all three species, we determined synteny blocks of autosomal origin and reconstructed the evolution of multiple sex-chromosomes. In addition, we identified W homologues of Z-linked orthologs and characterised their molecular differentiation. Our results suggest that the multiple sex-chromosome system evolved in a common ancestor as a result of dynamic genome reshuffling through repeated rearrangements between the sex chromosomes and autosomes, including translocations, fusions and fissions. Thus, the initial formation of neo-sex chromosomes could not have played a role in reproductive isolation between these Leptidea species. However, the subsequent species-specific fissions of several neo-sex chromosomes could have contributed to their reproductive isolation. Then, significantly increased numbers of Z-linked genes and independent neo-W chromosome degeneration could accelerate the accumulation of genetic incompatibilities between populations and promote their divergence resulting in speciation.

Chromosomal-level reference genome of Chinese peacock butterfly (Papilio bianor) based on third-generation DNA sequencing and Hi-C analysis

BACKGROUND

Papilio bianor Cramer, 1777 (commonly known as the Chinese peacock butterfly) (Insecta, Lepidoptera, Papilionidae) is a widely distributed swallowtail butterfly with a wide number of geographic populations ranging from the southeast of Russia to China, Japan, India, Vietnam, Myanmar, and Thailand. Its wing color consists of both pigmentary colored scales (black, reddish) and structural colored scales (iridescent blue or green dust). A high-quality reference genome of P. bianor is an important foundation for investigating iridescent color evolution, phylogeography, and the evolution of swallowtail butterflies.

FINDINGS

We obtained a chromosome-level de novo genome assembly of the highly heterozygous P. bianor using long Pacific Biosciences sequencing reads and high-throughput chromosome conformation capture technology. The final assembly is 421.52 Mb on 30 chromosomes (29 autosomes and 1 Z sex chromosome) with 13.12 Mb scaffold N50. In total, 15,375 protein-coding genes and 233.09 Mb of repetitive sequences were identified. Phylogenetic analyses indicated that P. bianor separated from a common ancestor of swallowtails ∼23.69-36.04 million years ago. Demographic history suggested that the population expansion of this species from the last interglacial period to the last glacial maximum possibly resulted from its decreased natural enemies and its adaptation to climate change during the glacial period.

CONCLUSIONS

We present a high-quality chromosome-level reference genome of P. bianor using long-read single-molecule sequencing and Hi-C-based chromatin interaction maps. Our results lay the foundation for exploring the genetic basis of special biological features of P. bianor and also provide a useful data source for comparative genomics and phylogenomics among butterflies and moths.