Next Article in Journal
Arabidopsis ANAC102, Chloroplastic or Nucleocytosolic Localization?
Next Article in Special Issue
Maps of Constitutive-Heterochromatin Distribution for Four Martes Species (Mustelidae, Carnivora, Mammalia) Show the Formative Role of Macrosatellite Repeats in Interspecific Variation of Chromosome Structure
Previous Article in Journal
Genetic Modifiers of Non-Penetrance and RNA Expression Levels in PRPF31-Associated Retinitis Pigmentosa in a Danish Cohort
Previous Article in Special Issue
Occurrence of Sex Chromosomes in Fish of the Genus Ancistrus with a New Description of Multiple Sex Chromosomes in the Ecuadorian Endemic Ancistrus clementinae (Loricariidae)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Elena A. Pazhenkova
1,* and
Vladimir A. Lukhtanov
2,*
1
Department of Biology, Biotechnical Faculty, University of Ljubljana, Večna pot 111, 1000 Ljubljana, Slovenia
2
Department of Karyosystematics, Zoological Institute of Russian Academy of Sciences, Universitetskaya nab. 1, St. Petersburg 199034, Russia
*
Authors to whom correspondence should be addressed.
Submission received: 20 December 2022 / Revised: 1 February 2023 / Accepted: 7 February 2023 / Published: 8 February 2023
(This article belongs to the Special Issue Chromosome Evolution and Karyotype Analysis)

Abstract

:
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 ligeaManiola 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.

1. Introduction

In contrast to monocentric chromosomes, in which centromeres are restricted to single regions on each chromosome, holocentric chromosomes are characterized by long centromeric areas (i.e., holocentromeres) [1,2,3,4]. Holocentromeres are distributed along the poleward surface of metaphase chromosomes, such that during cell divisions, microtubules bind to the chromosomes along almost their entire length [5]. The most recent study of holocentric chromosomes in beak-shaped species (Rhynchospora, family Cyperaceae) showed that these chromosomes are in fact polycentric [4], which significantly modifies previous models of holocentromere organization and indicates the need for a more detailed study of their structures in different organisms.
One of the most remarkable features of organisms with holocentric chromosomes is their ability to tolerate the fragmentation and fusion of chromosomes. In monocentric organisms, all chromosomal rearrangements must comply with the persistence of a functional centromere [6]; otherwise, they are not viable. For example, in monocentric organisms, fragmentation causes the formation of centric and acentric fragments, and fusion can cause the formation of dicentric chromosomes. In cell divisions, acentric fragments do not segregate and are consequently lost [7], whereas dicentric chromosomes are mostly broken [8]. In contrast, in holocentric organisms, fragmented chromosomes retain kinetic activity, and fused chromosomes are not dicentric and do not break during cell division [3,7].
Another remarkable feature found in many groups of holocentric organisms [9,10,11,12,13,14] is inverted meiosis. In standard conventional meiosis, the first meiotic division is reductional, resulting in the segregation of homologous chromosomes, whereas the second meiotic division is equational, resulting in the separation of sister chromatids. In inverted meiosis, the opposite order of these main meiotic events occurs, and homologous chromosomes separate after sister chromatids [15,16,17,18]. For butterflies of the genus Leptidea (Lepidoptera, family Pieridae), it has been shown that inverted meiosis facilitates proper chromosome segregation and rescues fertility and viability in hybrids heterozygous for chromosomal fusions and fissions [19,20]. Both the tolerance to chromosomal fusions and fragmentations and chromosomal heterozygosity are expected to favor the fixation of novel chromosomal rearrangements [3,7]. Therefore, holocentromeres have been hypothesized to drive rapid karyotype evolution [7,21].
This hypothesis, despite its popularity, has not yet received either reliable confirmation or refutation. In particular, comparative phylogenetic analysis has shown that chromosome number seems to evolve at equal rates in holocentric and monocentric clades [22]. There are numerous examples that demonstrate both the unusually high and unexpectedly conservative evolution of holocentric karyotypes. Thus, in the holocentric insect order Lepidoptera, which includes more than 150,000 species, in the blue butterfly genus Polyommatus (Lepidoptera, family Lycaenidae), in less than 5 million years of its evolution [23], a fan of chromosomally differentiated species arose. In these species, the diploid numbers (2n) of chromosomes vary from 2n = 20 to 2n = 452 [24,25,26]. A comparable level of interspecific variability in chromosome numbers is known for the satyrine genus Erebia (Lepidoptera, Nymphalidae, Satyrinae) (the diploid numbers vary from 2n = 14 to 2n = 104) and the tribe Ithomiini (Lepidoptera, family Nymphalidae) (the diploid numbers vary from 2n = 10 to 2n = 240) [27,28,29]. It would seem that these examples support the hypothesis that holocentricity accelerates karyotypic evolution. However, the vast majority of karyotyped Lepidoptera species demonstrate conservatism in chromosome numbers. Most species of Lepidoptera have a haploid number (n) of chromosomes, n = 31 [28]. This fundamental genome feature is ancestral to the order Lepidoptera [30,31] and has been preserved in most families for more than 100 million years of their evolution [32,33].
However, the stability of chromosome numbers is weak evidence of the stability of karyotypes, since chromosomal rearrangements, both intrachromosomal (e.g., inversions) and interchromosomal (e.g., translocations), can lead to the total reorganization of the gene order without changing the number of chromosomes. The high interspecies variability in the number of chromosomes is a good indication of rapid chromosomal evolution. Despite this, a simple count of the number of chromosomes is an unreliable tool for the analysis of the patterns and mechanisms of rapid karyotypic evolution. The chromosome number itself usually does not provide information on how individual chromosomes and the whole karyotypes originated or how they are related to chromosome complements of other species [6].
In the era of classical cytogenetics in the 20th century, for methodological reasons, the study of Lepidoptera karyotypes was mainly limited to counting numbers of chromosomes, mostly in meiotic metaphase plates [27,28,29,30,31]. The analytical tools of the 21st century (BAC-FISH mapping and whole-genome analyses) reveal chromosomal rearrangements, but so far, few species of butterflies and moths have been studied to the extent that allows the comparison of karyotypes [34,35,36,37,38,39,40]. The results of these studies are contradictory. The presence of multiple translocations between species with similar numbers of chromosomes has been directly demonstrated for the butterfly family Pieridae [38]. Multiple chromosomal rearrangements were detected between species and populations of the genus Leptidea (Pieridae) with different chromosome numbers [41]. Despite this, genomic data for other butterflies and moths (Lepidoptera) show a high level of chromosomal synteny when comparing chromosomes of different non-closely related species [33,36,42,43,44,45,46]. For example, extensive conserved synteny of genes was found between the karyotypes of Manduca sexta and Bombyx mori, belonging to two different moth families (Sphingidae and Bombycidae). Only a few rearrangements, including three inversions, three translocations and two fission/fusion events, were estimated to have occurred after the divergence of these families [47]. At the same time, M. sexta shows a high level of intraspecific differentiation, with the Arizona population being differentiated from other populations by two large inversions [48].
In the species-rich group of satyrine butterflies (Lepidoptera, subfamily Satyrinae of the family Nymphalidae), the modal and probably ancestral karyotype is 2n = 58, as is found in many genera [28]. In the genus Erebia, which belongs to the subfamily Satyrinae, there is a clear tendency toward a decrease in the diploid number of chromosomes up to 2n = 14 [49], against the background of the preservation of the ancestral 2n = 58 in many species of the genus [28]. The chromosomal rearrangements that may be behind these superficially similar or highly altered satyrine karyotypes remain unknown. In our work, we conducted a comparative analysis of genome-wide assemblies for three pairs of satyrine species. The first pair is represented by the species E. ligea and M. jurtina and demonstrates the preservation of the ancestral haploid chromosome number n = 29. The second pair is represented by the species E. aethiops and M. jurtina, demonstrating different chromosome numbers in different genera. The third pair is represented by the species E. ligea and E. aethiops, demonstrating a significant decrease in chromosome number within the same genus.

2. Materials and Methods

2.1. Materials

Chromosome-level genome assemblies of M. jurtina (Linnaeus, 1758), E. ligea (Linnaeus, 1758) and E. aethiops (Esper, 1777) (Insecta, Lepidoptera, Nymphalidae, Satyrinae) generated by the Darwin Tree of Life Project [50,51,52,53] and freely available upon their deposition in the European Nucleotide Archive (ENA) ( s://www.darwintreeoflife.org/wp-content/uploads/2020/03/DToL-Open-Data-Release-Policy-1.pdf, accessed on 30 January 2023) were used for the analysis of chromosomal macrosynteny and collinearity and the detection of chromosomal rearrangements (Table 1 and Table S1).

2.2. Detection of Macrosynteny and Chromosomal Rearrangements

We detected macrosynteny in the studied species through the pairwise whole-genome alignment of chromosome-level assemblies using minimap2 with the “asm20” preset, allowing higher divergence between genomes [54]. The advantages of this analysis are that it not only reveals syntenic regions but also visualizes chromosomal changes in terms and pictures of classical cytogenetics, i.e., provides graphic representations of chromosomes and chromosomal rearrangements. For each pair of species, this method made it possible to obtain multiple nucleotide alignments with varying lengths (in our case, up to 52 kb). Each such alignment was interpreted as a syntenic chromosomal block. We discarded alignments with a mapping quality of less than 60 and a length of less than 2 kb and visualized them as pairwise genomic dot plots with the pafr R package [55]. Pairwise comparisons were inspected visually. Karyotypes were drawn and colored using the RIdeogram R package [56]. At present, the circular plot is the standard for presenting comparative genomic data [57]. Therefore, we also visualized the pairwise synteny comparison as a circular layout (circos plot) showing the correspondence of each block between species pairs with the circlize R package [58].

3. Results

3.1. Analysis of Macrosynteny and Collinearity

The pairwise comparisons between the species were performed using the minimap2 algorithm [54] and revealed multiple syntenic blocks of nucleotides; the length of the detected blocks varied from 2 kb to 52 kb (Table 2). Within each compared pair of chromosomes, these syntenic blocks mainly formed continuous sequences (shown as diagonals in the genomic dot plots, Figure 1, Figure 2 and Figure 3). In most cases, the length of these diagonals and the order of blocks in them corresponded to the length and order of blocks in the compared chromosomes. This indicates not only the macrosynteny of the compared chromosomes (and, consequently, their homology) but also their collinearity.
Pairwise comparisons were also performed using a circus plot analysis [57]. This approach revealed the same patterns of relationships between species and chromosomes (Figure 4) as the minimap2 algorithm [54]. In the circus plot analysis, chromosomes of E. ligea and M. jurtina were found to be syntenic. Chromosomes 1–8, 18 and neo-Z of E. aethiops were found to be formed as a result of interchromosomal fusions.

3.2. Chromosomal Rearrangements: Inversions and Fusions

In the species pair M. jurtina–E. ligea, in eight autosomes and in the Z sex chromosome, the pairwise comparison revealed regions in which the sequence of syntenic blocks had the reverse orientation, thus indicating the presence of inversions (Figure 1). The inversions were located in subtelomeric (chromosomes 2, 6 22, 25 and Z of E. ligea) or interstitial (chromosomes 11, 17, 23 and 26 of E. ligea) positions. We did not find any interchromosomal rearrangements (translocations) between M. jurtina and E. ligea.
The comparison of E. ligea with E. aethiops shows that these taxa are separated by nine simple autosomal fusions (two smaller chromosomes of E. ligea correspond to a larger chromosome of E. aethiops) and the autosome 12–sex chromosome Z fusion resulting in a neo-Z chromosome in E. aethiops (Figure 3 and Figure 5). We also detected a single small subterminal inversion separating the Z chromosomes of this species (Figure 5). This is a different inversion than that found in the M. jurtina–E. ligea pair, since it has an almost terminal position, while in the M. jurtina–E. ligea pair, the inversion in the Z chromosome has a more interstitial localization.
The comparison between M. jurtina and E. aethiops (Figure 2) revealed nine autosomal fusions and one autosome–sex chromosome Z fusion. Taking into account the comparison between E. ligea and E. aethiops (Figure 3), we can conclude that all of these chromosome fusions arose in the lineage leading to the species E. aethiops, after its separation from E. ligea.
Since no inversions have been identified between E. ligea and E. aethiops, except for one terminal Z chromosome inversion (Figure 3), we can expect that all eight autosomal inversions and one subterminal Z chromosome inversion between M. jurtina and E. ligea (Figure 1) accumulated during evolution before the separation of E. ligea and E. aethiops. Therefore, one would expect that the same nine inversions (eight autosomal and one subterminal Z chromosome) would be found in the comparison between M. jurtina and E. aethiops. In fact, in the M. jurtina–E. aethiops comparison, seven inversions were revealed (Figure 2), which were identical to the inversions in the M. jurtina and E. ligea pair (in chromosomes 4jurtina, 7jurtina, 8jurtina, 21jurtina, 23jurtina, 25jurtina and Zjurtina). In the comparison between M. jurtina and E. aethiops, the expected inversions in chromosomes 19jurtina and 21jurtina were not detected. A possible explanation is that the inversions in chromosomes 4jurtina, 7jurtina, 8jurtina, 21jurtina, 23jurtina, 25jurtina and Zjurtina were fixed in evolution, while the inversions in chromosomes 19jurtina and 21jurtina are not fixed and may even have been lost in E. aethiops.

4. Discussion

4.1. Karyotypes of Maniola and Erebia Butterflies

The karyotype of M. jurtina was first studied by Federley [59] by using a standard microscopic technique for butterflies from Finland. Federley established that in the first metaphase of meiosis, both males and females have 29 bivalents (n = 29). This count of the haploid number (n = 29) was then confirmed by the studies by Lorković for Zagreb (Croatia) and Fontainebleau (France) [60] and Bigger for England [61].
The karyotype of E. ligea was also previously studied by Federley [59]. For butterflies from Finland, Federley identified 29 chromosome units in oogenesis at the stage of the first metaphase of meiosis (n = 29) [59]. In a population from Hokkaido (Japan), known as the subspecies E. l. rishirizana, Saitoh and Abe [62] found 28 bivalents in male meiosis and 56 chromosomes in male mitosis, that is, one pair of chromosomes fewer than in Finland. The chromosome-scale genome assemblies confirm the haploid number n = 29 for both species [51,52]. They also provide the first information on female heterogamety and sex chromosome systems in these species. M. jurtina has 28 pairs of autosomes, as well as sex chromosomes ZZ in males and ZW in females. In E. ligea, the whole genome was assembled for the male only. The male chromosome set has 28 autosomes and a pair of ZZ sex chromosomes.
The haploid number n = 29 is predominant in the subfamily Satyrinae [28], to which the genera Maniola and Erebia belong. As a modal number, n = 29 occurs in the vast majority of genera of the Palearctic Satyrinae [28,63,64]. The Neotropical satyrines, particularly the basal species, also have n = 29 as a weak modal number [65]. The African satyrine taxa have a strong modal n = 28; however, n = 29 also occurs in a few genera [65,66,67,68,69]. Based on the available data, it can be assumed that n = 29 (and, less likely, n = 28) is the ancestral number of chromosomes for the subfamily Satyrinae as a whole. The age of the subfamily is about 65 million years [70], during which this karyotype was preserved. Thus, in general, we can speak about the stability of chromosome numbers as one of the main trends in the karyotypic evolution of the Satyrinae.
Against the background of this stability, the species-rich satyrine genus Erebia is very unusual. Within this genus, haploid chromosome numbers vary from n = 7 in Erebia aethiopellus [49] to n = 51–52 in E. iranica [71]. Accordingly, the expected range of variability in diploid numbers is from 2n = 14 to 2n = 104, although care must be taken when interpreting the number of chromosome elements in meiosis in terms of diploid numbers. This is due to the fact that in normal meiosis in butterflies, multivalents and univalents can occur in meiosis, for example, as a consequence of the Z0 sex chromosome determination mechanism [32] or the presence of multiple sex chromosomes [34,35], as well as a consequence of polymorphism for chromosomal fusions and fissions [41,72].
The karyotype of E. aethiops was first studied by Lorković [60], who, based on the study of a few metaphases in a single male from Croatia, established the haploid number of n = 21. The chromosome-scale genome assembly shows that, in Scotland, this species has 18 pairs of autosomes and a pair of sex chromosomes: ZZ in males and ZW in females [53]. Our analysis shows that the Z chromosome of E. aethiops resulted from the fusion of the ancestral Z with an autosome and is thus the neo-Z chromosome [2,73]. This neo-Z chromosome is the largest in the set, and the tiny W chromosome is the smallest element in the set. Genomic data (Supplementary Table S1) show that the W chromosome of E. aethiops is 12 times smaller than the Z chromosome and 5 times smaller than the smallest (18th) autosome.

4.2. Dynamic Evolution of Holocentric Chromosomes in Satyrine Butterflies

Due to the extremely high level of interspecific karyotypic differentiation, the genus Erebia has become a model for studying rapid karyotypic evolution and the role of chromosomal changes in speciation and post-speciation evolution [74,75]. In this regard, it can be compared with another group of butterflies, the sister genera Polyommatus and Lysandra (Lycaenidae), in which there is an even higher level of interspecific differences in chromosome numbers [18,24,25,26,76,77]. For both Erebia and Polyommatus, attempts were made to identify the patterns, trends and possible mechanisms of karyotypic evolution using comparative phylogenetic analysis. In both cases, chromosome numbers were used as traits (proxy features for karyotype) for the analyses [74,75,76,77]. This approach has a serious limitation: although the decrease and increase in the number of chromosomes depend on chromosome fusions and fissions, the same chromosome numbers do not mean that the karyotypes are identical. The chromosome number itself does not provide information on how the whole karyotype originated or how chromosome complements of different species are related [6]. Robust evolutionary analysis should be based on a library of chromosomal rearrangements that differentiate the species under comparison. In our research, we took the first step toward creating such a library.
The comparison between M. jurtina and E. ligea revealed a high level of macrosynteny in all 29 chromosomes. We did not find any interchromosomal rearrangements (for example, translocations) between these species. Moreover, all of these chromosomes are mostly collinear: that is, they retain the same order of the studied chromosomal blocks. The exceptions are nine small inversions (eight in autosomes and one in the Z chromosome) that accumulated between these taxa over 30 million years of their independent evolution [78]. Thus, in the karyotype evolution of the ManiolaErebia lineage, an element of conservatism is observed, that is, the absence of large interchromosomal rearrangements, resulting in the preservation of the number of chromosomes. This conservatism is combined with a more dynamic repatterning of the gene order due to inversions.
A high level of macrosynteny does not necessarily indicate a complete absence of structural chromosomal evolution. It is known that homologous chromosomes can be differentiated through the accumulation of microarrangements, such as insertions and deletions [45,79]. Our data show that, against the background of the conservation of macrosynteny and collinearity, the compared species differ greatly (from 402 to 506 Mb, Table 1) in the total genome size. Most likely, these differences arose as a result of microinsertions and microdeletions of DNA repeats, leading to karyotype divergence, but not changing the gene order.
The evolution of karyotypes was more rapid in the lineage E. ligeaE. aethiops. The two species split about 12 million years ago [80]. A relatively rapid decrease in the number of chromosomes in E. aethiops occurred due to simple fusions, and we did not find any other large intrachromosomal rearrangements in autosomes. The evolution of the sex chromosome was more complex and included, in addition to fusion, the Z chromosome inversion. This finding is highly compatible with previous observations that, in Lepidoptera, sex chromosomes show a more dynamic structural evolution than autosomes [81,82,83,84,85,86,87]. Sex chromosomes are known to play a special role in the formation of reproductive barriers between species [88,89]. They evolve faster than autosomes [90,91], resulting in hybrid sterility genes that are preferentially localized on sex chromosomes [88,89]. We hypothesize that this exceptional role of the Z chromosome may be enhanced through sex chromosome–autosome fusions and Z chromosome inversions, i.e., via rearrangements that are known to promote reproductive isolation between nascent and closely related species [92].

4.3. Inversions and Holocentric Model of Chromosomal Speciation

The standard model of holocentric chromosomal speciation [19,93] is based on two main ideas. The first idea is that the holocentric nature of chromosomes contributes to the emergence and subsequent fixation of chromosomal fusions and divisions [3,7,20,75]. The second idea is that chromosomal fusions and fissions facilitate the divergence of incipient species through hybrid sterility and/or recombination repression mechanisms [19,94].
Here, based on our findings, we argue that not only fusions and fissions but also inversions play a significant role in holocentric-chromosome-based speciation. As our data show, inversions are the most common macroscale rearrangements in satyrine butterflies. The role of inversions in the origin and maintenance of reproductive isolation through a recombination suppression mechanism is well documented [92,95]. However, in monocentric organisms, the fixation of widespread paracentric inversions can be difficult. In a heterozygote for paracentric inversion, crossing over within the inverted region leads to the formation of an acentric chromatid and a dicentric chromatid. Both recombinants face problems. The acentric chromatid may be lost, and the dicentric recombinant generates a dicentric bridge during anaphase [2,92]. This problem is absent in holocentric organisms, providing an additional opportunity for karyotype evolution and the formation of new species.

5. Conclusions

We conclude that chromosomal evolution is dynamic in satyrine butterflies, even in the lineage that preserves the ancestral haploid chromosome number, n = 29. We hypothesize that the exceptional role of Z chromosomes in speciation may be further enhanced by inversions and sex chromosome–autosome fusions. Based on our analysis, we argue that not only fusions and fissions but also inversions are drivers of the holocentromere-mediated mode of chromosomal speciation.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/genes14020437/s1. Table S1: Studied species and chromosomes.

Author Contributions

Conceptualization, E.A.P. and V.A.L.; methodology, V.A.L.; software, E.A.P.; formal analysis, E.A.P.; resources, V.A.L.; writing—original draft preparation, V.A.L.; writing—review and editing, V.A.L. and E.A.P.; visualization, V.A.L.; funding acquisition, V.A.L. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by the Russian Science Foundation Grant No. 19-14-00202 (Continuation) to the Zoological Institute of the Russian Academy of Sciences (St. Petersburg, Russian Federation).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All of the analyzed chromosome-level assemblies are available via the GenBank links provided (Table S1).

Acknowledgments

We thank Valentina Kuznetsova (Zoological Institute of the Russian Academy of Sciences) for numerous discussions on the structure and evolution of holocentric chromosomes in insects.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schrader, F. Notes on the mitotic behavior of long chromosomes. Cytologia 1935, 6, 422–430. [Google Scholar] [CrossRef]
  2. White, M.J.D. Animal Cytology and Evolution, 3rd ed.; Cambridge University Press: Cambridge, UK, 1973; p. 468. [Google Scholar]
  3. Mandrioli, M.; Manicardi, G.C. Holocentric chromosomes. PLoS Genet. 2020, 16, e1008918. [Google Scholar] [CrossRef]
  4. Hofstatter, P.G.; Thangavel, G.; Lux, T.; Neumann, P.; Vondrak, T.; Novak, P.; Zhang, M.; Costa, L.; Castellani, M.; Scott, A.; et al. Repeat-based holocentromeres influence genome architecture and karyotype evolution. Cell 2022, 185, 3153–3168.e18. [Google Scholar] [CrossRef] [PubMed]
  5. Melters, D.P.; Paliulis, L.V.; Korf, I.F.; Chan, S.W. Holocentric chromosomes: Convergent evolution, meiotic adaptations, and genomic analysis. Chromosome Res. 2012, 20, 579–593. [Google Scholar] [CrossRef]
  6. Mayrose, I.; Lysak, M.A. The evolution of chromosome numbers: Mechanistic models and experimental approaches. Genome Biol. Evol. 2022, 13, evaa220. [Google Scholar] [CrossRef]
  7. Jankowska, M.; Fuchs, J.; Klocke, E.; Fojtová, M.; Polanská, P.; Fajkus, J.; Schubert, V.; Houben, A. Holokinetic centromeres and efficient telomere healing enable rapid karyotype evolution. Chromosoma 2015, 124, 519–528. [Google Scholar] [CrossRef]
  8. Haber, J.E.; Thorburn, P.; Rogers, D. Meiotic and mitotic behavior of dicentric chromosomes in Saccharomyces cerevisiae. Genetics 1983, 106, 185–205. [Google Scholar] [CrossRef]
  9. Bongiorni, S.; Fiorenzo, P.; Pippoletti, D.; Prantera, G. Inverted meiosis and meiotic drive in mealybugs. Chromosoma 2004, 112, 331–341. [Google Scholar] [CrossRef]
  10. Cabral, G.; Marques, A.; Schubert, V.; Pedrosa-Harand, A.; Schlögelhofer, P. Chiasmatic and achiasmatic inverted meiosis of plants with holocentric chromosomes. Nat. Commun. 2014, 5, 5070. [Google Scholar] [CrossRef]
  11. Heckmann, S.; Jankowska, M.; Schubert, V.; Kumke, K.; Ma, W.; Houben, A. Alternative meiotic chromatid segregation in the holocentric plant Luzula elegans. Nat. Commun. 2014, 5, 4979. [Google Scholar] [CrossRef] [Green Version]
  12. Manicardi, G.C.; Mandrioli, M.; Blackman, R.L. The cytogenetic architecture of the aphid genome. Biol. Rev. Camb. Philos. Soc. 2015, 90, 112–125. [Google Scholar] [CrossRef]
  13. Golub, N.V.; Golub, V.B.; Kuznetsova, V.G. New data on karyotypes of lace bugs (Tingidae, Cimicomorpha, Hemiptera) with analysis of the 18S rDNA clusters distribution. Comp. Cytogenet. 2018, 12, 515–528. [Google Scholar] [CrossRef]
  14. Lukhtanov, V.A.; Dantchenko, A.V.; Khakimov, F.R.; Sharafutdinov, D.; Pazhenkova, E.A. Karyotype evolution and flexible (conventional versus inverted) meiosis in insects with holocentric chromosomes: A case study based on Polyommatus butterflies. Biol. J. Linn. Soc. 2020, 130, 683–699. [Google Scholar] [CrossRef]
  15. Bogdanov, Y.F. Inverted meiosis and its place in the evolution of sexual reproduction pathways. Genetika 2016, 52, 541–560. [Google Scholar] [CrossRef]
  16. Lenormand, T.; Engelstaedter, J.; Johnston, S.E.; Wijnker, E.; Haag, C.R. Evolutionary mysteries in meiosis. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2016, 371, 20160001. [Google Scholar] [CrossRef]
  17. Loidl, J. Conservation and variability of meiosis across the eukaryotes. Ann. Rev. Genet. 2016, 50, 293–316. [Google Scholar] [CrossRef]
  18. Archetti, M. Inverted meiosis and the evolution of sex by loss of complementation. J. Evol. Biol. 2020, 33, 460–467. [Google Scholar] [CrossRef]
  19. Lukhtanov, V.A.; Dincă, V.; Friberg, M.; Šíchová, J.; Olofsson, M.; Vila, R.; Marec, F.; Wiklund, C. Versatility of multivalent orientation, inverted meiosis, and rescued fitness in holocentric chromosomal hybrids. Proc. Natl. Acad. Sci. USA 2018, 115, E9610–E9619. [Google Scholar] [CrossRef]
  20. Lukhtanov, V.A.; Dincă, V.; Friberg, M.; Vila, R.; Wiklund, C. Incomplete sterility of chromosomal hybrids: Implications for karyotype evolution and homoploid hybrid speciation. Front. Genet. 2020, 11, 583827. [Google Scholar] [CrossRef]
  21. Bureš, P.; Zedek, F.; Marková, M. Holocentric chromosomes. In Plant Genome Diversity; Greilhuber, J., Dolezel, J., Wendel, J.F., Eds.; Springer: Vienna, Austria, 2013; Volume 2, pp. 187–208. [Google Scholar]
  22. Ruckman, S.N.; Jonika, M.M.; Casola, C.; Blackmon, H. Chromosome number evolves at equal rates in holocentric and monocentric clades. PLoS Genet. 2020, 16, e1009076. [Google Scholar] [CrossRef]
  23. Talavera, G.; Lukhtanov, V.A.; Pierce, N.E.; Vila, R. DNA barcodes combined with multilocus data of representative taxa can generate reliable higher-level phylogenies. Syst. Biol. 2022, 71, 382–395. [Google Scholar] [CrossRef]
  24. Lukhtanov, V.A.; Kandul, N.P.; Plotkin, J.B.; Dantchenko, A.V.; Haig, D.; Pierce, N.E. Reinforcement of pre-zygotic isolation and karyotype evolution in Agrodiaetus butterflies. Nature 2005, 436, 385–389. [Google Scholar] [CrossRef]
  25. Talavera, G.; Lukhtanov, V.A.; Rieppel, L.; Pierce, N.E.; Vila, R. In the shadow of phylogenetic uncertainty: The recent diversification of Lysandra butterflies through chromosomal change. Mol. Phylogenet. Evol. 2013, 69, 469–478. [Google Scholar] [CrossRef]
  26. Lukhtanov, V.A. The blue butterfly Polyommatus (Plebicula) atlanticus (Lepidoptera, Lycaenidae) holds the record of the highest number of chromosomes in the non-polyploid eukaryotic organisms. Comp. Cytogenet. 2015, 9, 683–690. [Google Scholar] [CrossRef]
  27. de Lesse, H. Spéciation et variation chromosomique chez les Lépidoptères Rhopalocères. Annls Sci. Nat. Zool. Sér. 1960, 2, 1–223. [Google Scholar]
  28. Robinson, R. Lepidoptera Genetics; Pergamon Press: Oxford, UK, 1971. [Google Scholar]
  29. Brown, K.S.; von Schoultz, B.; Suomalainen, E. Chromosome evolution in Neotropical Danainae and Ithomiinae (Lepidoptera). Hereditas 2004, 141, 216–236. [Google Scholar] [CrossRef]
  30. Beliajeff, N.K. Die Chromosomenkomplexe und ihre Beziehung zur Phylogenie bei den Lepidopteren. Z. Für Inductive Abstamm. Und Verer. 1930, 14, 369–399. [Google Scholar] [CrossRef]
  31. Suomalainen, E. Chromosome evolution in Lepidoptera. Chromosome Today 1969, 2, 132–138. [Google Scholar]
  32. Lukhtanov, V.A. Sex chromatin and sex chromosome systems in nonditrysian Lepidoptera (Insecta). J. Zool. Syst. Evol. Res. 2000, 38, 73–79. [Google Scholar] [CrossRef]
  33. Ahola, V.; Lehtonen, R.; Somervuo, P.; Salmela, L.; Koskinen, P.; Rastas, P.; Välimäki, N.; Paulin, L.; Kvist, J.; Wahlberg, N.; et al. The Glanville fritillary genome retains an ancient karyotype and reveals selective chromosomal fusions in Lepidoptera. Nat. Commun. 2014, 5, 4737. [Google Scholar] [CrossRef]
  34. Šíchová, J.; Voleníková, A.; Dincă, V.; Nguyen, P.; Vila, R.; Sahara, K.; Marec, F. Dynamic karyotype evolution and unique sex determination systems in Leptidea wood white butterflies. BMC Evol. Biol. 2015, 15, 89. [Google Scholar] [CrossRef] [PubMed]
  35. Šíchová, J.; Ohno, M.; Dincă, V.; Watanabe, M.; Sahara, K.; Marec, F. Fissions, fusions, and translocations shaped the karyotype and multiple sex chromosome constitution of the northeast-Asian wood white butterfly, Leptidea amurensis. Biol. J. Linn. Soc. Lond. 2016, 118, 457–471. [Google Scholar] [CrossRef] [Green Version]
  36. Sahara, K.; Yoshido, A.; Marec, F.; Fuková, I.; Zhang, H.B.; Wu, C.C.; Goldsmith, M.R.; Yasukochi, Y. Conserved synteny of genes between chromosome 15 of Bombyx mori and a chromosome of Manduca sexta shown by five-color BAC-FISH. Genome 2007, 50, 1061–1065. [Google Scholar] [CrossRef] [PubMed]
  37. Sahara, K.; Yoshido, A.; Shibata, F.; Fujikawa-Kojima, N.; Okabe, T.; Tanaka-Okuyama, M.; Yasukochi, Y. FISH identification of Helicoverpa armigera and Mamestra brassicae chromosomes by BAC and fosmid probes. Insect Biochem. Mol. Biol. 2013, 43, 644–653. [Google Scholar] [CrossRef]
  38. Hill, J.; Rastas, P.; Hornett, E.A.; Neethiraj, R.; Clark, N.; Morehouse, N.; de la Paz Celorio-Mancera, M.; Cols, J.C.; Dircksen, H.; Meslin, C.; et al. Unprecedented reorganization of holocentric chromosomes provides insights into the enigma of lepidopteran chromosome evolution. Sci. Adv. 2019, 5, eaau3648. [Google Scholar] [CrossRef]
  39. Cicconardi, F.; Lewis, J.J.; Martin, S.H.; Reed, R.D.; Danko, C.G.; Montgomery, S.H. Chromosome fusion affects genetic diversity and evolutionary turnover of functional loci but consistently depends on chromosome size. Mol. Biol. Evol. 2021, 38, 4449–4462. [Google Scholar] [CrossRef]
  40. Mackintosh, A.; Laetsch, D.R.; Baril, T.; Foster, R.; Dincă, V.; Vila, R.; Hayward, A.; Lohse, K. The genome sequence of the lesser marbled fritillary, Brenthis ino, and evidence for a segregating neo-Z chromosome. G3-Genes Genom. Genet. 2022, 12, jkac069. [Google Scholar] [CrossRef]
  41. Höök, L.; Näsvall, K.; Vila, R.; Wiklund, C.; Backström, N. High-density linkage maps and chromosome level genome assemblies unveil direction and frequency of extensive structural rearrangements in wood white butterflies (Leptidea spp.). Chromosome Res. 2023, 31, 2. [Google Scholar] [CrossRef]
  42. Pringle, E.G.; Baxter, S.W.; Webster, C.L.; Papanicolaou, A.; Lee, S.F.; Jiggins, C.D. Synteny and chromosome evolution in the Lepidoptera: Evidence from mapping in Heliconius melpomene. Genetics 2007, 177, 417–426. [Google Scholar] [CrossRef]
  43. Papa, R.; Morrison, C.M.; Walters, J.R.; Counterman, B.A.; Chen, R.; Halder, G.; Ferguson, L.; Chamberlain, N.; Ffrench-Constant, R.; Kapan, D.D.; et al. Highly conserved gene order and numerous novel repetitive elements in genomic regions linked to wing pattern variation in Heliconius butterflies. BMC Genom. 2008, 9, 345. [Google Scholar] [CrossRef]
  44. Beldade, P.; Saenko, S.V.; Pul, N.; Long, A.D. A gene-based linkage map for Bicyclus anynana butterflies allows for a comprehensive analysis of synteny with the lepidopteran reference genome. PLoS Genet. 2009, 5, e1000366. [Google Scholar] [CrossRef]
  45. d’Alençon, E.; Sezutsu, H.; Legeai, F.; Permal, E.; Bernard-Samain, S.; Gimenez, S.; Gagneur, C.; Cousserans, F.; Shimomura, M.; Brun-Barale, A.; et al. Extensive synteny conservation of holocentric chromosomes in Lepidoptera despite high rates of local genome rearrangements. Proc. Natl. Acad. Sci. USA 2010, 107, 7680–7685. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Davey, J.W.; Chouteau, M.; Barker, S.L.; Maroja, L.; Baxter, S.W.; Simpson, F.; Merrill, R.M.; Joron, M.; Mallet, J.; Dasmahapatra, K.K.; et al. Major improvements to the Heliconius melpomene genome assembly used to confirm 10 chromosome fusion events in 6 million years of butterfly evolution. G3-Genes Genom. Genet. 2016, 6, 695–708. [Google Scholar] [CrossRef]
  47. Yasukochi, Y.; Tanaka-Okuyama, M.; Shibata, F.; Yoshido, A.; Marec, F.; Wu, C.; Zhang, H.; Goldsmith, M.R.; Sahara, K. Extensive conserved synteny of genes between the karyotypes of Manduca sexta and Bombyx mori revealed by BAC-FISH mapping. PLoS ONE 2009, 4, e7465. [Google Scholar] [CrossRef]
  48. Mongue, A.J.; Kawahara, A.Y. Population differentiation and structural variation in the Manduca sexta genome across the United States. G3-Genes Genom. Genet. 2022, 12, jkac047. [Google Scholar] [CrossRef]
  49. de Lesse, H. Caractères et repartition en France d’Erebia aethiopellus Hoffmsg. et E. mnestra Hb. Alexanor 1959, 1, 72–81. [Google Scholar]
  50. The Darwin Tree of Life Project Consortium. Sequence locally, think globally: The Darwin tree of life project. Proc. Natl. Acad. Sci. USA 2022, 119, e2115642118. [Google Scholar] [CrossRef]
  51. Lohse, K.; Weir, J.; Darwin Tree of Life Barcoding Collective; Wellcome Sanger Institute Tree of Life Programme; Wellcome Sanger Institute Scientific Operations: DNA Pipelines Collective; Tree of Life Core Informatics Collective; Darwin Tree of Life Consortium. The genome sequence of the meadow brown, Maniola jurtina (Linnaeus, 1758). Wellcome Open. Res. 2021, 6, 296. [Google Scholar] [CrossRef]
  52. Lohse, K.; Hayward, A.; Laetsch, D.R.; Vila, R.; Lucek, K.; Wellcome Sanger Institute Tree of Life Programme; Wellcome Sanger Institute Scientific Operations: DNA Pipelines Collective; Tree of Life Core Informatics Collective; Darwin Tree of Life Consortium. The genome sequence of the Arran brown, Erebia ligea (Linnaeus, 1758) [version 1; peer review: Awaiting peer review]. Wellcome Open. Res. 2022, 7, 259. [Google Scholar] [CrossRef]
  53. Lohse, O.; Lohse, K.; Wellcome Sanger Institute Tree of Life programme; Wellcome Sanger Institute Scientific Operations: DNA Pipelines collective; Tree of Life Core Informatics c ollective; Darwin Tree of Life Consortium. The genome sequence of the scotch argus butterfly, Erebia aethiops (Esper, 1777). Wellcome Open. Res. 2022, 7, 217. [Google Scholar] [CrossRef]
  54. Li, H. Minimap2: Pairwise alignment for nucleotide sequences. Bioinformatics 2018, 34, 3094–3100. [Google Scholar] [CrossRef] [PubMed]
  55. Winter, D.; Lee, K.; Cox, M. 2020 pafr: Read, Manipulate and Visualize Pairwise mApping Format. Available online: https://dwinter.github.io/pafr/ (accessed on 18 December 2022).
  56. Hao, Z.; Lv, D.; Ge, Y.; Shi, J.; Weijers, D.; Yu, G.; Chen, J. RIdeogram: Drawing SVG graphics to visualize and map genome-wide data on the idiograms. PeerJ Comput. Sci. 2020, 6, e251. [Google Scholar] [CrossRef] [PubMed]
  57. Krzywinski, M.; Schein, J.; Birol, I.; Connors, J.; Gascoyne, R.; Horsman, D.; Jones, S.J.; Marra, M.A. Circos: An information aesthetic for comparative genomics. Genome Res. 2009, 19, 1639–1645. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Gu, Z.; Gu, L.; Eils, R.; Schlesner, M.; Brors, B. circlize implements and enhances circular visualization in R. Bioinformatics 2014, 30, 2811–2812. [Google Scholar] [CrossRef] [PubMed]
  59. Federley, H. Chromosomenzahlen finnländischer Lepidopteren I. Rhopalocera. Hereditas 1938, 24, 397–464. [Google Scholar] [CrossRef]
  60. Lorković, Z. Die Chromosomezahlen in der Spermatogenese der Tagfalter. Chromosoma 1941, 2, 155–191. [Google Scholar] [CrossRef]
  61. Bigger, T.R.L. Chromosome numbers of Lepidoptera. Part I. Ent. Gaz. 1960, 11, 149–152. [Google Scholar]
  62. Saitoh, K.; Abe, A. The chromosomes of Erebia ligea rishirizana (Nymphalidae, Satyrinae). Nota lepid. 1997, 20, 326–329. [Google Scholar]
  63. Larsen, T.B. Chromosome numbers and notes on testicular morphology of some Lebanese Rhopalocera (Insecta Lepidoptera). Entomol. Scand. 1975, 6, 253–260. [Google Scholar] [CrossRef]
  64. Lorković, Z. The butterfly chromosomes and their application in systematics and phylogeny. In Butterflies of Europe; Kudrna, O., Ed.; Aula Verlag: Wiesbaden, Germany, 1990; Volume 2, pp. 332–396. [Google Scholar]
  65. Brown, K.; Freitas, A.V.L.; von Schoultz, B.; Saura, A.O.; Saura, A. Chromosomal evolution of South American frugivorous butterflies in the Satyroid clade (Nymphalidae: Charaxinae, Morphinae and Satyrinae). Biol. J. Linn. Soc. 2007, 92, 467–481. [Google Scholar] [CrossRef]
  66. de Lesse, H. Formules chromosomiques de Lépidoptères Rhopalocères d’Afrique du Nord. Bull. Soc. Entomol. Fr. 1967, 72, 20–25. [Google Scholar] [CrossRef]
  67. de Lesse, H. Formules chromosomiques de Lépidoptères Rhopalocères d’Uganda et du Kenya. Ann. Soc. Entomol. Fr. 1968, 4, 581–599. [Google Scholar]
  68. de Lesse, H.; Condamin, M. Formules chromosomiques de quelques Lépidoptères Rhopalocères du Sénégal. Bull. IFAN Ser. A 1962, 24, 464–473. [Google Scholar]
  69. de Lesse, H.; Condamin, M. Formules chromosomiques de quelques Lépidoptères Rhopalocères du Sénégal et de Côte d’Ivoire. Bull. IFAN Ser. A 1965, 27, 1089–1094. [Google Scholar]
  70. Espeland, M.; Breinholt, J.; Willmott, K.R.; Warren, A.D.; Vila, R.; Toussaint, E.F.A.; Maunsell, S.C.; Aduse-Poku, K.; Talavera, G.; Eastwood, R.; et al. A comprehensive and dated phylogenomic analysis of butterflies. Curr. Biol. 2018, 28, 770–778. [Google Scholar] [CrossRef]
  71. Lukhtanov, V.; Sourakov, A.; Tikhonov, V.; Zakharov, E. Taxonomic rearrangement of the Erebia tyndarus species group (Lepidoptera, Nymphalidae, Satyrinae) based on an analysis of COI barcodes, morphology and geographic distribution. Folia Biol. 2019, 67, 149–157. [Google Scholar] [CrossRef]
  72. Lukhtanov, V.A.; Dinca, V.; Talavera, G.; Vila, R. Unprecedented within-species chromosome number cline in the Wood White butterfly Leptidea sinapis and its significance for karyotype evolution and speciation. BMC Evol. Biol. 2011, 11, 109. [Google Scholar] [CrossRef]
  73. Hejníčková, M.; Dalíková, M.; Potocký, P.; Tammaru, T.; Trehubenko, M.; Kubíčková, S.; Marec, F.; Zrzavá, M. Degenerated, undifferentiated, rearranged, lost: High variability of sex chromosomes in Geometridae (Lepidoptera) identified by sex chromatin. Cells 2021, 10, 2230. [Google Scholar] [CrossRef]
  74. Lucek, K. Evolutionary mechanisms of varying chromosome numbers in the radiation of Erebia butterflies. Genes 2018, 9, 166. [Google Scholar] [CrossRef]
  75. de Vos, J.M.; Augustijnen, H.; Batscher, L.; Lucek, K. Speciation through chromosomal fusion and fission in Lepidoptera. Philos. Trans. R. Soc. B Biol. Sci. 2020, 375, 20190539. [Google Scholar] [CrossRef]
  76. Kandul, N.P.; Lukhtanov, V.A.; Pierce, N.E. Karyotypic diversity and speciation in Agrodiaetus butterflies. Evolution 2007, 61, 546–559. [Google Scholar] [CrossRef] [PubMed]
  77. Vershinina, A.O.; Lukhtanov, V.A. Evolutionary mechanisms of runaway chromosome number change in Agrodiaetus butterflies. Sci. Rep. 2017, 7, 8199. [Google Scholar] [CrossRef] [PubMed]
  78. Peña, C.; Nylin, S.; Wahlberg, N. The radiation of Satyrini butterflies (Nymphalidae: Satyrinae): A challenge for phylogenetic methods. Zool. J. Linn. Soc. 2011, 161, 64–87. [Google Scholar] [CrossRef]
  79. Ruggieri, A.A.; Livraghi, L.; Lewis, J.J.; Evans, E.; Cicconardi, F.; Hebberecht, L.; Ortiz-Ruiz, Y.; Montgomery, S.H.; Ghezzi, A.; Rodriguez-Martinez, J.A.; et al. A butterfly pan-genome reveals a large amount of structural variation underlies the evolution of chromatin accessibility. Genome Res. 2022, 32, 1862–1875. [Google Scholar] [CrossRef]
  80. Peña, C.; Witthauer, H.; Klečková, I.; Fric, Z.; Wahlberg, N. Adaptive radiations in butterflies: Evolutionary history of the genus Erebia (Nymphalidae: Satyrinae). Biol. J. Linn. Soc. 2015, 116, 449–467. [Google Scholar] [CrossRef]
  81. Nilsson, N.O.; Löfstedt, C.; Dävring, L. Unusual sex-chromosome inheritance in six species of small ermine moths (Yponomeuta, Yponomeutidae, Lepidoptera). Hereditas 2008, 108, 259–265. [Google Scholar] [CrossRef]
  82. Nguyen, P.; Sýkorová, M.; Šichová, J.; Kůta, V.; Daliková, M.; Frydrychová, R.Č.; Neven, L.G.; Sahara, K.; Marec, F. Neo-sex chromosomes and adaptive potential in tortricid pests. Proc. Natl. Acad. Sci. USA 2013, 110, 6931–6936. [Google Scholar] [CrossRef]
  83. Šíchová, J.; Nguyen, P.; Dalíková, M.; Marec, F. Chromosomal evolution in tortricid moths: Conserved karyotypes with diverged features. PLoS ONE 2013, 8, e64520. [Google Scholar] [CrossRef] [Green Version]
  84. Van’t Hof, A.; Nguyen, P.; Dalíková, M.; Edmonds, N.; Marec, F.; Saccheri, I.J. Linkage map of the peppered moth, Biston betularia (Lepidoptera, Geometridae): A model of industrial melanism. Heredity 2013, 110, 283–295. [Google Scholar] [CrossRef]
  85. Dalíková, M.; Zrzavá, M.; Hladová, I.; Nguyen, P.; Šonský, I.; Flegrová, M.; Kubíčková, S.; Voleníková, A.; Kawahara, A.Y.; Peters, R.S.; et al. New insights into the evolution of the W chromosome in Lepidoptera. J. Hered. 2017, 108, 709–719. [Google Scholar] [CrossRef]
  86. Mongue, A.J.; Nguyen, P.; Voleníková, A.; Walters, J.R. Neo-sex chromosomes in the monarch butterfly, Danaus plexippus. G3-Genes Genom. Genet. 2017, 7, 3281–3294. [Google Scholar] [CrossRef] [PubMed]
  87. Pazhenkova, E.A.; Lukhtanov, V.A. Chromosomal conservatism vs. chromosomal megaevolution: Enigma of karyotypic evolution in Lepidoptera. BioRxiv 2022, 5, 494852. [Google Scholar] [CrossRef]
  88. Coyne, J.A.; Orr, H.A. Speciation; Sinauer Associates: Sunderland, MA, USA, 2004; p. 545. [Google Scholar]
  89. Presgraves, D.C. Sex chromosomes and speciation in Drosophila. Trends Genet. 2008, 24, 336–343. [Google Scholar] [CrossRef]
  90. Charlesworth, B.; Campos, J.L.; Jackson, B.C. Faster-X evolution: Theory and evidence from Drosophila. Mol. Ecol. 2018, 27, 3753–3771. [Google Scholar] [CrossRef] [PubMed]
  91. Mongue, A.J.; Hansen, M.E.; Walters, J.R. Support for faster and more adaptive Z chromosome evolution in two divergent lepidopteran lineages. Evolution 2022, 76, 332–345. [Google Scholar] [CrossRef] [PubMed]
  92. King, M. Species Evolution: The Role of Chromosomal Change; Cambridge University Press: Cambridge, UK, 1993; p. 360. [Google Scholar]
  93. Lucek, K.; Augustijnen, H.; Escudero, M. A holocentric twist to chromosomal speciation? Trends Ecol. Evol. 2022, 37, 655–662. [Google Scholar] [CrossRef] [PubMed]
  94. Faria, R.; Navarro, A. Chromosomal speciation revisited: Rearranging theory with pieces of evidence. Trends Ecol. Evol. 2010, 25, 660–669. [Google Scholar] [CrossRef] [PubMed]
  95. Kirkpatrick, M.; Barton, N. Chromosome inversions, local adaptation and speciation. Genetics 2006, 173, 419–434. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Genome of E. ligea (horizontal axis) plotted against genome of M. jurtina (vertical axis). It can be seen that chromosome 1ligea corresponds to chromosome 1jurtina, chromosome 2ligea corresponds to chromosome 4jurtina, chromosome 3ligea corresponds to chromosome 5jurtina, chromosome 4ligea corresponds to chromosome 3jurtina, and so on, and chromosome Zligea corresponds to chromosome Zjurtina. Each short stroke represents a single syntenic block (a supported alignment), positioned according to its position in the genomes of the compared species. All compared chromosomes are syntenic. Twenty pairs of chromosomes are collinear. Nine pairs of chromosomes show the presence of inversions (shown by red arrows).
Figure 1. Genome of E. ligea (horizontal axis) plotted against genome of M. jurtina (vertical axis). It can be seen that chromosome 1ligea corresponds to chromosome 1jurtina, chromosome 2ligea corresponds to chromosome 4jurtina, chromosome 3ligea corresponds to chromosome 5jurtina, chromosome 4ligea corresponds to chromosome 3jurtina, and so on, and chromosome Zligea corresponds to chromosome Zjurtina. Each short stroke represents a single syntenic block (a supported alignment), positioned according to its position in the genomes of the compared species. All compared chromosomes are syntenic. Twenty pairs of chromosomes are collinear. Nine pairs of chromosomes show the presence of inversions (shown by red arrows).
Genes 14 00437 g001
Figure 2. Genome of M. jurtina (horizontal axis) plotted against genome of E. aethiops (vertical axis). It can be seen that chromosome 1jurtina corresponds to chromosome 11aethiops, chromosome 2jurtina corresponds to chromosome 10aethiops, chromosome 3jurtina corresponds to chromosome 12aethiops, chromosome 4jurtina corresponds to chromosome 5aethiops, and so on, and chromosomes Z + 10ligea correspond to chromosome neo-Zaethiops. Each short stroke represents a single syntenic block (a supported alignment), positioned according to its position in the genomes of the compared species. All compared chromosomes are syntenic. Seven pairs of chromosomes show the presence of inversions (shown by red arrows).
Figure 2. Genome of M. jurtina (horizontal axis) plotted against genome of E. aethiops (vertical axis). It can be seen that chromosome 1jurtina corresponds to chromosome 11aethiops, chromosome 2jurtina corresponds to chromosome 10aethiops, chromosome 3jurtina corresponds to chromosome 12aethiops, chromosome 4jurtina corresponds to chromosome 5aethiops, and so on, and chromosomes Z + 10ligea correspond to chromosome neo-Zaethiops. Each short stroke represents a single syntenic block (a supported alignment), positioned according to its position in the genomes of the compared species. All compared chromosomes are syntenic. Seven pairs of chromosomes show the presence of inversions (shown by red arrows).
Genes 14 00437 g002
Figure 3. Genome of E. ligea (horizontal axis) plotted against genome of E. aethiops (vertical axis). It can be seen that chromosome 1ligea corresponds to chromosome 11aetiops, chromosomes 2 + 24ligea correspond to chromosome 5aethiops, chromosome 3ligea corresponds to chromosome 9aethiops, chromosome 4ligea corresponds to chromosome 12aethiops, and so on, and chromosomes 12 + Zligea correspond to chromosome neo-Zaethiops. Each short stroke represents a single syntenic block (a supported alignment), positioned according to its position in the genomes of the compared species. The red arrow indicates a small terminal inversion in the Z chromosome.
Figure 3. Genome of E. ligea (horizontal axis) plotted against genome of E. aethiops (vertical axis). It can be seen that chromosome 1ligea corresponds to chromosome 11aetiops, chromosomes 2 + 24ligea correspond to chromosome 5aethiops, chromosome 3ligea corresponds to chromosome 9aethiops, chromosome 4ligea corresponds to chromosome 12aethiops, and so on, and chromosomes 12 + Zligea correspond to chromosome neo-Zaethiops. Each short stroke represents a single syntenic block (a supported alignment), positioned according to its position in the genomes of the compared species. The red arrow indicates a small terminal inversion in the Z chromosome.
Genes 14 00437 g003
Figure 4. Circos plot showing synteny comparison between E. ligea and E. aehtiops (a), E. ligea and M. jurtina (b) and E. aethiops and M. jurtina (c). Each link corresponds to a single syntenic block (alignment) and is colored by E. ligea (a,b) and M. jurtina (c) chromosomes. Chromosomes 1–8, 18 and neo-Z of E. aethiops were formed as a result of interchromosomal fusions. Chromosomes of E. ligea and M. jurtina are syntenic.
Figure 4. Circos plot showing synteny comparison between E. ligea and E. aehtiops (a), E. ligea and M. jurtina (b) and E. aethiops and M. jurtina (c). Each link corresponds to a single syntenic block (alignment) and is colored by E. ligea (a,b) and M. jurtina (c) chromosomes. Chromosomes 1–8, 18 and neo-Z of E. aethiops were formed as a result of interchromosomal fusions. Chromosomes of E. ligea and M. jurtina are syntenic.
Genes 14 00437 g004
Figure 5. Chromosome comparison between E. ligea and E. aethiops. Karyotypes of E. ligea (n = 29) (lower row) and E. aethiops (n = 19) (upper row) are differentiated by 10 chromosomal fusions and a single terminal inversion in the Z chromosome. The serial numbers of chromosomes of each species are taken from GenBank (Supplementary Table S1) and correspond to their sizes.
Figure 5. Chromosome comparison between E. ligea and E. aethiops. Karyotypes of E. ligea (n = 29) (lower row) and E. aethiops (n = 19) (upper row) are differentiated by 10 chromosomal fusions and a single terminal inversion in the Z chromosome. The serial numbers of chromosomes of each species are taken from GenBank (Supplementary Table S1) and correspond to their sizes.
Genes 14 00437 g005
Table 1. GenBank IDs, karyotypes, sex chromosome systems and genome sizes in the studied species of the subfamily Satyrinae (according to GenBank data and published data [51,52,53]).
Table 1. GenBank IDs, karyotypes, sex chromosome systems and genome sizes in the studied species of the subfamily Satyrinae (according to GenBank data and published data [51,52,53]).
SpeciesAssembly IDDiploid Chromosome Number and Sex Chromosomes in FemalesDiploid Chromosome Number and Sex Chromosomes in MalesTotal Genome Size
M. jurtinailManJurt1.12n = 56 + ZW2n = 56 + ZZ402.0 Mb
E. ligeailEreLige1.2unstudied2n = 56 + ZZ506.4 Mb
E. aethiopsilEreAeth2.22n = 36 + ZW2n = 36 + ZZ473.4 Mb
Table 2. Properties of pairwise alignments (i.e., syntenic blocks) of the Satyrinae whole-genome assemblies after filtration.
Table 2. Properties of pairwise alignments (i.e., syntenic blocks) of the Satyrinae whole-genome assemblies after filtration.
ComparisonMaximum Block Length, bpSum of Block Lengths, bpNumber of Blocks
M. jurtina–E. ligea13,73932,722,76315,870
E. ligea–E. aethiops52,04567,667,01330,566
M. jurtina–E. aethiops15,30516,426,95216,357
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pazhenkova, E.A.; Lukhtanov, V.A. Whole-Genome Analysis Reveals the Dynamic Evolution of Holocentric Chromosomes in Satyrine Butterflies. Genes 2023, 14, 437. https://doi.org/10.3390/genes14020437

AMA Style

Pazhenkova EA, Lukhtanov VA. Whole-Genome Analysis Reveals the Dynamic Evolution of Holocentric Chromosomes in Satyrine Butterflies. Genes. 2023; 14(2):437. https://doi.org/10.3390/genes14020437

Chicago/Turabian Style

Pazhenkova, Elena A., and Vladimir A. Lukhtanov. 2023. "Whole-Genome Analysis Reveals the Dynamic Evolution of Holocentric Chromosomes in Satyrine Butterflies" Genes 14, no. 2: 437. https://doi.org/10.3390/genes14020437

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop