Temporal dynamics of faster neo-Z evolution in butterflies

The faster-Z/X hypothesis predicts that sex-linked genes should diverge faster than autosomal genes. However, studies across different lineages have shown mixed support for this effect. So far, most analyses have focused on old and well-differentiated sex chromosomes, but less is known about the divergence of more recently acquired neo-sex chromosomes. In Lepidoptera (moths and butterflies), Z-autosome fusions are frequent, but the evolutionary dynamics of neo-Z chromosomes have not been explored in detail. Here, we analyzed the faster-Z effect in Leptidea sinapis, a butterfly with three Z chromosomes. We show that the neo-Z chromosomes have been acquired stepwise, resulting in strata of differentiation and masculinization. While all Z chromosomes showed evidence of the faster-Z effect, selection for genes on the youngest neo-Z chromosome (Z3) appears to have been hampered by a largely intact, homologous neo-W chromosome. However, the intermediately aged neo-Z chromosome (Z2), which lacks W gametologs, showed fewer evolutionary constraints, resulting in particularly fast evolution. Our results therefore support that neo-sex chromosomes can constitute temporary hot-spots of adaptation and divergence. The underlying dynamics are likely causally linked to shifts in selective constraints, evolution of gene expression, and degeneration of W-linked gametologs which gradually expose Z-linked genes to selection.

Structure and Evolution of Ribosomal Genes of Insect Chromosomes

Currently, clusters of 45S and 5S ribosomal DNA (rDNA) have been studied in about 1000 and 100 species of the class Insecta, respectively. Although the number of insect species with known 45S rDNA clusters (also referred to as nucleolus-organizing regions, or NORs) constitutes less than 0.1 percent of the described members of this enormous group, certain conclusions can already be drawn. Since haploid karyotypes with single 45S and 5S rDNA clusters predominate in both basal and derived insect groups, this character state is apparently ancestral for the class Insecta in general. Nevertheless, the number, chromosomal location, and other characteristics of both 45S and 5S rDNA sites substantially vary across different species, and sometimes even within the same species. There are several main factors and molecular mechanisms that either maintain these parameters or alter them on the short-term and/or long-term scale. Chromosome structure (i.e., monocentric vs. holokinetic chromosomes), excessive numbers of rRNA gene copies per cluster, interactions with transposable elements, pseudogenization, and meiotic recombination are perhaps the most important among them.

Genomic evidence reveals three W-autosome fusions in Heliconius butterflies

Sex chromosomes are evolutionarily labile in many animals and sometimes fuse with autosomes, creating so-called neo-sex chromosomes. Fusions between sex chromosomes and autosomes have been proposed to reduce sexual conflict and to promote adaptation and reproductive isolation among species. Recently, advances in genomics have fuelled the discovery of such fusions across the tree of life. Here, we discovered multiple fusions leading to neo-sex chromosomes in the sapho subclade of the classical adaptive radiation of Heliconius butterflies. Heliconius butterflies generally have 21 chromosomes with very high synteny. However, the five Heliconius species in the sapho subclade show large variation in chromosome number ranging from 21 to 60. We find that the W chromosome is fused with chromosome 4 in all of them. Two sister species pairs show subsequent fusions between the W and chromosomes 9 or 14, respectively. These fusions between autosomes and sex chromosomes make Heliconius butterflies an ideal system for studying the role of neo-sex chromosomes in adaptive radiations and the degeneration of sex chromosomes over time. Our findings emphasize the capability of short-read resequencing to detect genomic signatures of fusion events between sex chromosomes and autosomes even when sex chromosomes are not explicitly assembled.

Lepidopteran Synteny Units reveal deep chromosomal conservation in butterflies and moths

DNA is compacted into individual particles or chromosomes that form the basic units of inheritance. However, different animals and plants have widely different numbers of chromosomes. This means that we cannot readily tell which chromosomes are related to which. Here, we describe a simple technique that looks at the similarity of genes on each chromosome and thus gives us a true picture of their homology or similarity through evolutionary time. We use this new system to look at the chromosomes of butterflies and moths or Lepidoptera. We term the associated synteny units, Lepidopteran Synteny Units (LSUs). Using a sample of butterfly and moth genomes from across evolutionary time, we show that LSUs form a simple and reliable method of tracing chromosomal homology back through time. Surprisingly, this technique reveals that butterfly and moth chromosomes show conserved blocks dating back to their sister group the Trichoptera. As Lepidoptera have holocentric chromosomes, it will be interesting to see if similar levels of synteny are shown in groups of animals with monocentric chromosomes. The ability to define homology via LSU analysis makes it considerably easier to approach many questions in chromosomal evolution.

Nascent evolution of recombination rate differences as a consequence of chromosomal rearrangements

Reshuffling of genetic variation occurs both by independent assortment of chromosomes and by homologous recombination. Such reshuffling can generate novel allele combinations and break linkage between advantageous and deleterious variants which increases both the potential and the efficacy of natural selection. Here we used high-density linkage maps to characterize global and regional recombination rate variation in two populations of the wood white butterfly (Leptidea sinapis) that differ considerably in their karyotype as a consequence of at least 27 chromosome fissions and fusions. The recombination data were compared to estimates of genetic diversity and measures of selection to assess the relationship between chromosomal rearrangements, crossing over, maintenance of genetic diversity and adaptation. Our data show that the recombination rate is influenced by both chromosome size and number, but that the difference in the number of crossovers between karyotypes is reduced as a consequence of a higher frequency of double crossovers in larger chromosomes. As expected from effects of selection on linked sites, we observed an overall positive association between recombination rate and genetic diversity in both populations. Our results also revealed a significant effect of chromosomal rearrangements on the rate of intergenic diversity change between populations, but limited effects on polymorphisms in coding sequence. We conclude that chromosomal rearrangements can have considerable effects on the recombination landscape and consequently influence both maintenance of genetic diversity and efficiency of selection in natural populations.

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 Role of Repetitive Sequences in Repatterning of Major Ribosomal DNA Clusters in Lepidoptera

Genes for major ribosomal RNAs (rDNA) are present in multiple copies mainly organized in tandem arrays. The number and position of rDNA loci can change dynamically and their repatterning is presumably driven by other repetitive sequences. We explored a peculiar rDNA organization in several representatives of Lepidoptera with either extremely large or numerous rDNA clusters. We combined molecular cytogenetics with analyses of second- and third-generation sequencing data to show that rDNA spreads as a transcription unit and reveal association between rDNA and various repeats. Furthermore, we performed comparative long read analyses among the species with derived rDNA distribution and moths with a single rDNA locus, which is considered ancestral. Our results suggest that satellite arrays, rather than mobile elements, facilitate homology-mediated spread of rDNA via either integration of extrachromosomal rDNA circles or ectopic recombination. The latter arguably better explains preferential spread of rDNA into terminal regions of lepidopteran chromosomes as efficiency of ectopic recombination depends on the proximity of homologous sequences to telomeres.

More hidden diversity in a cryptic species complex: a new subspecies of Leptideasinapis (Lepidoptera, Pieridae) from Northern Iran

A new subspecies of Leptideasinapis from Northern Iran, discovered by means of DNA barcoding, is described as Leptideasinapistabarestanassp. nov. The new subspecies is allopatric with respect to other populations of L.sinapis and is genetically distinct, appearing as a well-supported sister clade to all other populations in COI-based phylogenetic reconstructions. Details on karyotype, genitalia, ecology and behaviour for the new subspecies are given and a biogeographical speciation scenario is proposed.

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.

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.

On the Origin of Neo-Sex Chromosomes in the Neotropical Dragonflies Rhionaeschna bonariensis and R. planaltica (Aeshnidae, Odonata)

Odonata have holokinetic chromosomes. About 95% of species have an XX/X0 sex chromosome system, with heterogametic males. There are species with neo-XX/neo-XY sex chromosomes resulting from an X chromosome/autosome fusion. The genus Rhionaeschna includes 42 species found in the Americas. We analyzed the distribution of the nucleolar organizer region (NOR) using FISH with rDNA probes in Rhionaeschna bonariensis (n = 12 + neo-XY), R. planaltica (n = 7 + neo-XY), and Aeshna cyanea (n = 13 + X0). In R. bonariensis and A. cyanea, the NOR is located on a large pair of autosomes, which have a secondary constriction in the latter species. In R. planaltica, the NOR is located on the ancestral part of the neo-X chromosome. Meiotic analysis and FISH results in R. planaltica led to the conclusion that the neo-XY system arose by insertion of the ancestral X chromosome into an autosome. Genomic in situ hybridization, performed for the first time in Odonata, highlighted the entire neo-Y chromosome in meiosis of R. bonariensis, suggesting that it consists mainly of repetitive DNA. This feature and the terminal chiasma localization suggest an ancient origin of the neo-XY system. Our study provides new information on the origin and evolution of neo-sex chromosomes in Odonata, including new types of chromosomal rearrangements, NOR transposition, and heterochromatin accumulation.

Representing sex chromosomes in genome assemblies

Sex chromosomes have evolved hundreds of independent times across eukaryotes. As genome sequencing, assembly, and scaffolding techniques rapidly improve, it is now feasible to build fully phased sex chromosome assemblies. Despite technological advances enabling phased assembly of whole chromosomes, there are currently no standards for representing sex chromosomes when publicly releasing a genome. Furthermore, most computational analysis tools are unable to efficiently investigate their unique biology relative to autosomes. We discuss a diversity of sex chromosome systems and consider the challenges of representing sex chromosome pairs in genome assemblies. By addressing these issues now as technologies for full phasing of chromosomal assemblies are maturing, we can collectively ensure that future genome analysis toolkits can be broadly applied to all eukaryotes with diverse types of sex chromosome systems. Here we provide best practice guidelines for presenting a genome assembly that contains sex chromosomes. These guidelines can also be applied to other non-recombining genomic regions, such as S-loci in plants and mating-type loci in fungi and algae.

Multiple sex chromosomes in teleost fishes from a cytogenetic perspective: state of the art and future challenges

Despite decades of cytogenetic and genomic research of dynamic sex chromosome evolution in teleost fishes, multiple sex chromosomes have been largely neglected. In this review, we compiled available data on teleost multiple sex chromosomes, identified major trends in their evolution and suggest further trajectories in their investigation. In a compiled dataset of 440 verified records of fish sex chromosomes, we counted 75 multiple sex chromosome systems with 60 estimated independent origins. We showed that male-heterogametic systems created by Y-autosome fusion predominate and that multiple sex chromosomes are over-represented in the order Perciformes. We documented a striking difference in patterns of differentiation of sex chromosomes between male and female heterogamety and hypothesize that faster W sex chromosome differentiation may constrain sex chromosome turnover in female-heterogametic systems. We also found no significant association between the mechanism of multiple sex chromosome formation and percentage of uni-armed chromosomes in teleost karyotypes. Last but not least, we hypothesized that interaction between fish populations, which differ in their sex chromosomes, can drive the evolution of multiple sex chromosomes in fishes. This underlines the importance of broader inter-population sampling in studies of fish sex chromosomes. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Degenerated, Undifferentiated, Rearranged, Lost: High Variability of Sex Chromosomes in Geometridae (Lepidoptera) Identified by Sex Chromatin

Sex chromatin is a conspicuous body that occurs in polyploid nuclei of most lepidopteran females and consists of numerous copies of the W sex chromosome. It is also a cytogenetic tool used to rapidly assess the W chromosome presence in Lepidoptera. However, certain chromosomal features could disrupt the formation of sex chromatin and lead to the false conclusion that the W chromosome is absent in the respective species. Here we tested the sex chromatin presence in 50 species of Geometridae. In eight selected species with either missing, atypical, or normal sex chromatin patterns, we performed a detailed karyotype analysis by means of comparative genomic hybridization (CGH) and fluorescence in situ hybridization (FISH). The results showed a high diversity of W chromosomes and clarified the reasons for atypical sex chromatin, including the absence or poor differentiation of W, rearrangements leading to the neo-W emergence, possible association with the nucleolus, and the existence of multiple W chromosomes. In two species, we detected intraspecific variability in the sex chromatin status and sex chromosome constitution. We show that the sex chromatin is not a sufficient marker of the W chromosome presence, but it may be an excellent tool to pinpoint species with atypical sex chromosomes.

Flavors of Non-Random Meiotic Segregation of Autosomes and Sex Chromosomes

Segregation of chromosomes is a multistep process occurring both at mitosis and meiosis to ensure that daughter cells receive a complete set of genetic information. Critical components in the chromosome segregation include centromeres, kinetochores, components of sister chromatid and homologous chromosomes cohesion, microtubule organizing centres, and spindles. Based on the cytological work in the grasshopper Brachystola, it has been accepted for decades that segregation of homologs at meiosis is fundamentally random. This ensures that alleles on chromosomes have equal chance to be transmitted to progeny. At the same time mechanisms of meiotic drive and an increasing number of other examples of non-random segregation of autosomes and sex chromosomes provide insights into the underlying mechanisms of chromosome segregation but also question the textbook dogma of random chromosome segregation. Recent advances provide a better understanding of meiotic drive as a prominent force where cellular and chromosomal changes allow autosomes to bias their segregation. Less understood are mechanisms explaining observations that autosomal heteromorphism may cause biased segregation and regulate alternating segregation of multiple sex chromosome systems or translocation heterozygotes as an extreme case of non-random segregation. We speculate that molecular and cytological mechanisms of non-random segregation might be common in these cases and that there might be a continuous transition between random and non-random segregation which may play a role in the evolution of sexually antagonistic genes and sex chromosome evolution.

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.

The Role of Satellite DNAs in Genome Architecture and Sex Chromosome Evolution in Crambidae Moths

Tandem repeats are important parts of eukaryotic genomes being crucial e.g., for centromere and telomere function and chromatin modulation. In Lepidoptera, knowledge of tandem repeats is very limited despite the growing number of sequenced genomes. Here we introduce seven new satellite DNAs (satDNAs), which more than doubles the number of currently known lepidopteran satDNAs. The satDNAs were identified in genomes of three species of Crambidae moths, namely Ostrinia nubilalis, Cydalima perspectalis, and Diatraea postlineella, using graph-based computational pipeline RepeatExplorer. These repeats varied in their abundance and showed high variability within and between species, although some degree of conservation was noted. The satDNAs showed a scattered distribution, often on both autosomes and sex chromosomes, with the exception of both satellites in D. postlineella, in which the satDNAs were located at a single autosomal locus. Three satDNAs were abundant on the W chromosomes of O. nubilalis and C. perspectalis, thus contributing to their differentiation from the Z chromosomes. To provide background for the in situ localization of the satDNAs, we performed a detailed cytogenetic analysis of the karyotypes of all three species. This comparative analysis revealed differences in chromosome number, number and location of rDNA clusters, and molecular differentiation of sex chromosomes.

Meiotic analyses show adaptations to maintenance of fertility in X1Y1X2Y2X3Y3X4Y4X5Y5 system of amazon frog Leptodactylus pentadactylus (Laurenti, 1768)

Heterozygous chromosomal rearrangements can result in failures during the meiotic cycle and the apoptosis of germline, making carrier individuals infertile. The Amazon frog Leptodactylus pentadactylus has a meiotic multivalent, composed of 12 sex chromosomes. The mechanisms by which this multi-chromosome system maintains fertility in males of this species remain undetermined. In this study we investigated the meiotic behavior of this multivalent to understand how synapse, recombination and epigenetic modifications contribute to maintaining fertility and chromosomal sexual determination in this species. Our sample had 2n = 22, with a ring formed by ten chromosomes in meiosis, indicating a new system of sex determination for this species (X1Y1X2Y2X3Y3X4Y4X5Y5). Synapsis occurs in the homologous terminal portion of the chromosomes, while part of the heterologous interstitial regions performed synaptic adjustment. The multivalent center remains asynaptic until the end of pachytene, with interlocks, gaps and rich-chromatin in histone H2A phosphorylation at serine 139 (γH2AX), suggesting transcriptional silence. In late pachytene, paired regions show repair of double strand-breaks (DSBs) with RAD51 homolog 1 (Rad51). These findings suggest that Rad51 persistence creates positive feedback at the pachytene checkpoint, allowing meiosis I to progress normally. Additionally, histone H3 trimethylation at lysine 27 in the pericentromeric heterochromatin of this anuran can suppress recombination in this region, preventing failed chromosomal segregation. Taken together, these results indicate that these meiotic adaptations are required for maintenance of fertility in L. pentadactylus.

Karyotype evolution and flexible (conventional versus inverted) meiosis in insects with holocentric chromosomes: a case study based on Polyommatus butterflies

The Polyommatus butterflies have holocentric chromosomes, which are characterized by kinetic activity distributed along the entire chromosome length, and the highest range of haploid chromosome numbers (n) known within a single eukaryotic genus (from n = 10 to n = 226). Previous analyses have shown that these numbers most likely evolved gradually from an ancestral karyotype, in accordance with the Brownian motion model of chromosome change accumulation. Here we studied chromosome sets within a monophyletic group of previously non-karyotyped Polyommatus species. We demonstrate that these species have a limited interspecific chromosome number variation from n = 16 to n = 25, which is consistent with the Brownian motion model prediction. We also found intra- and interpopulation variation in the chromosome numbers. These findings support the model of karyotype evolution through the gradual accumulation of neutral or weakly underdominant rearrangements that can persist in the heterozygous state within a population. For Polyommatus poseidonides we report the phenomenon of flexible meiosis in which the chromosome multivalents are able to undergo either conventional or inverted meiosis within the same individual. We hypothesise that the ability to invert the order of the meiotic events may be adaptive and can facilitate proper chromosome segregation in chromosomal heterozygotes, thus promoting rapid karyotype evolution.

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.

Cytogenetic markers applied to cytotaxonomy in two soybean pests: Anticarsia gemmatalis (Hübner, 1818) and Chrysodeixis includens (Walker, 1858)

Anticarsia gemmatalis (Hübner, 1818) and Chrysodeixis includens (Walker, 1858) are species of Lepidoptera that cause great damages in the soybean plantations of Brazil. Despite the importance they have in this regard, there are no studies on the chromosomal organization of these species and recently, A. gemmatalis, which belonged to the Noctuidae family, was allocated to the Erebidae family. Therefore, the objective of this paper was to analyze, through conventional and molecular cytogenetic markers, both species of Lepidoptera. A 2n = 62 was observed, with ZZ/ZW sex chromosome system and holokinetic chromosomes for both species. There was homogeneity in the number of 18S rDNA sites for both species. However, variations in heterochromatin distribution were observed between both species. The cytogenetic analyses enabled separation of the species, corroborating the transference of A. gemmatalis, from the family Noctuidae to the family Erebidae, suggesting new cytotaxonomic characteristics.

Advances and Challenges of Using the Sterile Insect Technique for the Management of Pest Lepidoptera

Over the past 30 years, the sterile insect technique (SIT) has become a regular component of area-wide integrated pest management (AW-IPM) programs against several major agricultural pests and vectors of severe diseases. The SIT-based programs have been especially successful against dipteran pests. However, the SIT applicability for controlling lepidopteran pests has been challenging, mainly due to their high resistance to the ionizing radiation that is used to induce sterility. Nevertheless, the results of extensive research and currently operating SIT programs show that most problems with the implementation of SIT against pest Lepidoptera have been successfully resolved. Here, we summarize the cytogenetic peculiarities of Lepidoptera that should be considered in the development and application of SIT for a particular pest species. We also discuss the high resistance of Lepidoptera to ionizing radiation, and present the principle of derived technology based on inherited sterility (IS). Furthermore, we present successful SIT/IS applications against five major lepidopteran pests, and summarize the results of research on the quality control of reared and released insects, which is of great importance for their field performance. In the light of new research findings, we also discuss options for the development of genetic sexing strains, which is a challenge to further improve the applicability of SIT/IS against selected lepidopteran pests.

The determinants of genetic diversity in butterflies

Under the neutral theory, genetic diversity is expected to increase with population size. While comparative analyses have consistently failed to find strong relationships between census population size and genetic diversity, a recent study across animals identified a strong correlation between propagule size and genetic diversity, suggesting that r-strategists that produce many small offspring, have greater long-term population sizes. Here we compare genome-wide genetic diversity across 38 species of European butterflies (Papilionoidea), a group that shows little variation in reproductive strategy. We show that genetic diversity across butterflies varies over an order of magnitude and that this variation cannot be explained by differences in current abundance, propagule size, host or geographic range. Instead, neutral genetic diversity is negatively correlated with body size and positively with the length of the genetic map. This suggests that genetic diversity is determined both by differences in long-term population size and the effect of selection on linked sites.

Deciphering the Origin and Evolution of the X1X2Y System in Two Closely-Related Oplegnathus Species (Oplegnathidae and Centrarchiformes)

Oplegnathus fasciatus and O. punctatus (Teleostei: Centrarchiformes: Oplegnathidae), are commercially important rocky reef fishes, endemic to East Asia. Both species present an X1X2Y sex chromosome system. Here, we investigated the evolutionary forces behind the origin and differentiation of these sex chromosomes, with the aim to elucidate whether they had a single or convergent origin. To achieve this, conventional and molecular cytogenetic protocols, involving the mapping of repetitive DNA markers, comparative genomic hybridization (CGH), and whole chromosome painting (WCP) were applied. Both species presented similar 2n, karyotype structure and hybridization patterns of repetitive DNA classes. 5S rDNA loci, besides being placed on the autosomal pair 22, resided in the terminal region of the long arms of both X1 chromosomes in females, and on the X1 and Y chromosomes in males. Furthermore, WCP experiments with a probe derived from the Y chromosome of O. fasciatus (OFAS-Y) entirely painted the X1 and X2 chromosomes in females and the X1, X2, and Y chromosomes in males of both species. CGH failed to reveal any sign of sequence differentiation on the Y chromosome in both species, thereby suggesting the shared early stage of neo-Y chromosome differentiation. Altogether, the present findings confirmed the origin of the X1X2Y sex chromosomes via Y-autosome centric fusion and strongly suggested their common origin.

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

Chromosome evolution presents an enigma in the mega-diverse Lepidoptera. Most species exhibit constrained chromosome evolution with nearly identical haploid chromosome counts and chromosome-level gene collinearity among species more than 140 million years divergent. However, a few species possess radically inflated chromosomal counts due to extensive fission and fusion events. To address this enigma of constraint in the face of an exceptional ability to change, we investigated an unprecedented reorganization of the standard lepidopteran chromosome structure in the green-veined white butterfly (Pieris napi). We find that gene content in P. napi has been extensively rearranged in large collinear blocks, which until now have been masked by a haploid chromosome number close to the lepidopteran average. We observe that ancient chromosome ends have been maintained and collinear blocks are enriched for functionally related genes suggesting both a mechanism and a possible role for selection in determining the boundaries of these genome-wide rearrangements.

Genome-wide RAD sequencing to identify a sex-specific marker in Chinese giant salamander Andrias davidianus

Background

Chinese giant salamander Andrias davidianus is an endangered species. The success of artificial breeding provides a useful way to protect this species. However, the method to identify the sex and mechanism of sex determination were unclear which hinder the improvement of the artificial breeding. Detection of a sex specific marker provides an effective approach to identify genetic sex and investigate the sex determination mechanism.

Results

We used restriction-site-associated DNA (RAD) sequencing to isolate a sex-specific genetic marker in A. davidianus to expand knowledge of the sex determination mechanism. Four male and four female specimens were subjected to RAD sequencing, which generated 934,072,989 reads containing approximately 134.4 Gb of sequences. The first round of comparison of the assembled sequence against the opposite sex raw reads revealed 19,097 female and 17,994 male unmatched sequences. Subsequently, 19,097 female sequences were subjected to a BLAST search against male genomic data, which revealed 308 sequences unmapped to the male genome. One hundred of these were randomly selected and validated by PCR in five male and five female specimens, and four putative sex-specific sequences were produced. Further validation was performed by PCR in another 24 females and 24 males, and all female individuals exhibited the expected specific bands, while the males did not. To apply the sex-specific marker, three specimens reversed from genetic female to physiological male were found in a group exposed to elevated temperature, and 13 individuals reversed from genetic male to physiological female were obtained in a 17β-estradiol exposed group.

Conclusion

This is the first report of a sex-specific marker in A. davidianus and may have potential for elucidation of its sex determination mechanism and, hence, its conservation.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5771-5) contains supplementary material, which is available to authorized users.

Sex Chromosome Turnover in Moths of the Diverse Superfamily Gelechioidea

Sex chromosomes play a central role in genetics of speciation and their turnover was suggested to promote divergence. In vertebrates, sex chromosome-autosome fusions resulting in neo-sex chromosomes occur frequently in male heterogametic taxa (XX/XY), but are rare in groups with female heterogamety (WZ/ZZ). We examined sex chromosomes of seven pests of the diverse lepidopteran superfamily Gelechioidea and confirmed the presence of neo-sex chromosomes in their karyotypes. Two synteny blocks, which correspond to autosomes 7 (LG7) and 27 (LG27) in the ancestral lepidopteran karyotype exemplified by the linkage map of Biston betularia (Geometridae), were identified as sex-linked in the tomato leafminer, Tuta absoluta (Gelechiidae). Testing for sex-linkage performed in other species revealed that while LG7 fused to sex chromosomes in a common ancestor of all Gelechioidea, the second fusion between the resulting neo-sex chromosome and the other autosome is confined to the tribe Gnoreschemini (Gelechiinae). Our data accentuate an emerging pattern of high incidence of neo-sex chromosomes in Lepidoptera, the largest clade with WZ/ZZ sex chromosome system, which suggest that the paucity of neo-sex chromosomes is not an intrinsic feature of female heterogamety. Furthermore, LG7 contains one of the major clusters of UDP-glucosyltransferases, which are involved in the detoxification of plant secondary metabolites. Sex chromosome evolution in Gelechioidea thus supports an earlier hypothesis postulating that lepidopteran sex chromosome-autosome fusions can be driven by selection for association of Z-linked preference or host-independent isolation genes with larval performance and thus can contribute to ecological specialization and speciation of moths.

Versatility of multivalent orientation, inverted meiosis, and rescued fitness in holocentric chromosomal hybrids

Chromosomal rearrangements (e.g., fusions/fissions) have the potential to drive speciation. However, their accumulation in a population is generally viewed as unlikely, because chromosomal heterozygosity should lead to meiotic problems and aneuploid gametes. Canonical meiosis involves segregation of homologous chromosomes in meiosis I and sister chromatid segregation during meiosis II. In organisms with holocentric chromosomes, which are characterized by kinetic activity distributed along almost the entire chromosome length, this order may be inverted depending on their metaphase I orientation. Here we analyzed the evolutionary role of this intrinsic versatility of holocentric chromosomes, which is not available to monocentric ones, by studying F₁ to F₄ hybrids between two chromosomal races of the Wood White butterfly (Leptidea sinapis), separated by at least 24 chromosomal fusions/fissions. We found that these chromosomal rearrangements resulted in multiple meiotic multivalents, and, contrary to the theoretical prediction, the hybrids displayed relatively high reproductive fitness (42% of that of the control lines) and regular behavior of meiotic chromosomes. In the hybrids, we also discovered inverted meiosis, in which the first and critical stage of chromosome number reduction was replaced by the less risky stage of sister chromatid separation. We hypothesize that the ability to invert the order of the main meiotic events facilitates proper chromosome segregation and hence rescues fertility and viability in chromosomal hybrids, potentially promoting dynamic karyotype evolution and chromosomal speciation.

Sex Chromosomes of the Iconic Moth Abraxas grossulariata (Lepidoptera, Geometridae) and Its Congener A. sylvata

The magpie moth, Abraxas grossulariata, is an iconic species in which female heterogamety was discovered at the beginning of the 20th century. However, the sex chromosomes of this species have not yet been cytologically identified. We describe the sex chromosomes of A. grossulariata and its congener, A. sylvata. Although these species split only around 9.5 million years ago, and both species have the expected WZ/ZZ chromosomal system of sex determination and their sex chromosomes share the major ribosomal DNA (rDNA) representing the nucleolar organizer region (NOR), we found major differences between their karyotypes, including between their sex chromosomes. The species differ in chromosome number, which is 2n = 56 in A. grossularita and 2n = 58 in A. sylvata. In addition, A. grossularita autosomes exhibit massive autosomal blocks of heterochromatin, which is a very rare phenomenon in Lepidoptera, whereas the autosomes of A. sylvata are completely devoid of distinct heterochromatin. Their W chromosomes differ greatly. Although they are largely composed of female-specific DNA sequences, as shown by comparative genomic hybridization, cross-species W-chromosome painting revealed considerable sequence differences between them. The results suggest a relatively rapid molecular divergence of Abraxas W chromosomes by the independent spreading of female-specific repetitive sequences.

More sex chromosomes than autosomes in the Amazonian frog Leptodactylus pentadactylus

Heteromorphic sex chromosomes are common in eukaryotes and largely ubiquitous in birds and mammals. The largest number of multiple sex chromosomes in vertebrates known today is found in the monotreme platypus (Ornithorhynchus anatinus, 2n = 52) which exhibits precisely 10 sex chromosomes. Interestingly, fish, amphibians, and reptiles have sex determination mechanisms that do or do not involve morphologically differentiated sex chromosomes. Relatively few amphibian species carry heteromorphic sex chromosomes, and when present, they are frequently represented by only one pair, either XX:XY or ZZ:ZW types. Here, in contrast, with several evidences, from classical and molecular cytogenetic analyses, we found 12 sex chromosomes in a Brazilian population of the smoky jungle frog, designated as Leptodactylus pentadactylus Laurenti, 1768 (Leptodactylinae), which has a karyotype with 2n = 22 chromosomes. Males exhibited an astonishing stable ring-shaped meiotic chain composed of six X and six Y chromosomes. The number of sex chromosomes is larger than the number of autosomes found, and these data represent the largest number of multiple sex chromosomes ever found among vertebrate species. Additionally, sequence and karyotype variation data suggest that this species may represent a complex of species, in which the chromosomal rearrangements may possibly have played an important role in the evolution process.

Molecular phylogenetic reconstruction and localization of the (TTAGG)n telomeric repeats in the chromosomes of Acromyrmex striatus (Roger, 1863) suggests a lower ancestral karyotype for leafcutter ants (Hymenoptera)

Chromosome counts and karyotype characterization have proved to be important features of a genome. Chromosome changes during the diversification of ants might play an important role, given the diversity and success of Formicidae. Comparative karyotype analyses on ants have enriched and helped ant systematics. Among leafcutter ants, two major chromosome counts have been described, one frequent in Atta Fabricius, 1804 (2n = 22 in all Atta spp. whose karyotype is known) and the other frequent in Acromyrmex Mayr, 1865 (2n = 38 in the majority of species whose karyotype is known). The main exception is Acromyrmex striatus (Roger, 1863), which harbors a diploid chromosome set of 22. Here we describe the use of fluorescence in situ hybridization (FISH) with telomeric probes with (TTAGG)₆ repeats to describe the telomere composition of A. striatus and to recover potential interstitial non-telomeric signals that may reflect fusion events during the evolution of leafcutter lineage from 38 to 22 chromosomes. Further, we reconstruct the ancestral chromosome numbers of the leafcutter clade based on a recently proposed molecular phylogenetic hypothesis and phylogenomic tree. Distinct signals have been observed in both extremities on the telomere chromosomes of A. striatus. Non-telomeric signals have not been retrieved in our analysis. It could be supposed that the low-numbered karyotype indeed represents the ancestral chromosome number of leafcutters. The phylogenetic reconstruction also recovered a low chromosome number from the diverse approaches implemented, suggesting that n = 11 is the most likely ancestral karyotype of the leafcutter ants and is a plesiomorphic feature shared between A. striatus and Atta spp.

A new butterfly species from south Russia revealed through chromosomal and molecular analysis of the Polyommatus (Agrodiaetus) damonides complex (Lepidoptera, Lycaenidae)

Finding a new species is a rare event in easy-to-see and well-studied organisms like butterflies, especially if they inhabit well-explored areas such as the Western Palaearctic. However, even in this region, gaps in taxonomic knowledge still exist and here we report such a discovery. Using a combined analysis of chromosomal and molecular markers we demonstrate that Polyommatus blue populations from Daghestan (South Russia), previously identified as P. aserbeidschanus, represent in fact a new species which is described here as P. australorossicussp. n. We also show that the enigmatic Polyommatus damonides described as a form of Polyommatus damone and later considered as an entity similar to P. poseidon or P. ninae is conspecific with a taxon previously known as P. elbursicus. As a result of our study, we propose several taxonomic changes within the P. damonides species complex and suggest the following new combinations: P. damonides elbursicus Forster, 1956, comb. n. and P. damonides gilanensis Eckweiler, 2002, comb. n.

Rapid Increase in Genome Size as a Consequence of Transposable Element Hyperactivity in Wood-White (Leptidea) Butterflies

Characterizing and quantifying genome size variation among organisms and understanding if genome size evolves as a consequence of adaptive or stochastic processes have been long-standing goals in evolutionary biology. Here, we investigate genome size variation and association with transposable elements (TEs) across lepidopteran lineages using a novel genome assembly of the common wood-white (Leptidea sinapis) and population re-sequencing data from both L. sinapis and the closely related L. reali and L. juvernica together with 12 previously available lepidopteran genome assemblies. A phylogenetic analysis confirms established relationships among species, but identifies previously unknown intraspecific structure within Leptidea lineages. The genome assembly of L. sinapis is one of the largest of any lepidopteran taxon so far (643 Mb) and genome size is correlated with abundance of TEs, both in Lepidoptera in general and within Leptidea where L. juvernica from Kazakhstan has considerably larger genome size than any other Leptidea population. Specific TE subclasses have been active in different Lepidoptera lineages with a pronounced expansion of predominantly LINEs, DNA elements, and unclassified TEs in the Leptidea lineage after the split from other Pieridae. The rate of genome expansion in Leptidea in general has been in the range of four Mb/Million year (My), with an increase in a particular L. juvernica population to 72 Mb/My. The considerable differences in accumulation rates of specific TE classes in different lineages indicate that TE activity plays a major role in genome size evolution in butterflies and moths.

Evolutionary mechanisms of runaway chromosome number change in Agrodiaetus butterflies

Despite predictions of the classic, hybrid-sterility model of chromosomal speciation, some organisms demonstrate high rate of karyotype evolution. This rate is especially impressive in Agrodiaetus butterflies that rapidly evolved the greatest chromosome number diversity known in animal kingdom within a single subgenus. Here we analyzed karyotype evolution in Agrodiaetus using phylogenetic comparative methods. We found that chromosome numbers possess a strong phylogenetic signal. This disproves the chromosome megaevolution model that proposes multiple chromosome rearrangements to accumulate independently in each of closely related species. We found that Brownian motion gives a more adequate description of observed trait changes than Ornstein-Uhlenbeck model. This indicates that chromosome numbers evolve via random walk along branches of the phylogeny. We discovered a correlation between karyotype changes and phylogeny branch lengths. This gradual pattern is inconsistent with the hybrid-sterility model which, due to association of major chromosome changes with cladogenetic events, predicts a high degree of punctualism in karyotype evolution. Thus, low underdominace of chromosomal rearrangements and/or prevalence of the recombination-suppression model over the hybrid-sterility model of chromosome speciation are the most common engines of the runaway chromosome number change observed.

Evolution of genome size and genomic GC content in carnivorous holokinetics (Droseraceae)

BACKGROUND AND AIMS

Studies in the carnivorous family Lentibulariaceae in the last years resulted in the discovery of the smallest plant genomes and an unusual pattern of genomic GC content evolution. However, scarcity of genomic data in other carnivorous clades still prevents a generalization of the observed patterns. Here the aim was to fill this gap by mapping genome evolution in the second largest carnivorous family, Droseraceae, where this evolution may be affected by chromosomal holokinetism in Drosera METHODS: The genome size and genomic GC content of 71 Droseraceae species were measured by flow cytometry. A dated phylogeny was constructed, and the evolution of both genomic parameters and their relationship to species climatic niches were tested using phylogeny-based statistics.

KEY RESULTS

The 2C genome size of Droseraceae varied between 488 and 10 927 Mbp, and the GC content ranged between 37·1 and 44·7 %. The genome sizes and genomic GC content of carnivorous and holocentric species did not differ from those of their non-carnivorous and monocentric relatives. The genomic GC content positively correlated with genome size and annual temperature fluctuations. The genome size and chromosome numbers were inversely correlated in the Australian clade of Drosera CONCLUSIONS: Our results indicate that neither carnivory (nutrient scarcity) nor the holokinetism have a prominent effect on size and DNA base composition of Droseraceae genomes. However, the holokinetic drive seems to affect karyotype evolution in one of the major clades of Drosera Our survey confirmed that the evolution of GC content is tightly connected with the evolution of genome size and also with environmental conditions.

Flow cytometry may allow microscope-independent detection of holocentric chromosomes in plants

Two chromosomal structures, known as monocentric and holocentric chromosomes, have evolved in eukaryotes. Acentric fragments of monocentric chromosomes are unequally distributed to daughter cells and/or lost, while holocentric fragments are inherited normally. In monocentric species, unequal distribution should generate chimeras of cells with different nuclear DNA content. We investigated whether such differences in monocentric species are detectable by flow cytometry (FCM) as (i) a decreased nuclear DNA content and (ii) an increased coefficient of variance (CV) of the G1 peak after gamma radiation-induced fragmentation. We compared 13 monocentric and 9 holocentric plant species. Unexpectedly, monocentrics and holocentrics did not differ with respect to parameters (i) and (ii) in their response to gamma irradiation. However, we found that the proportion of G2 nuclei was highly elevated in monocentrics after irradiation, while holocentrics were negligibly affected. Therefore, we hypothesize that DNA-damaging agents induce cell cycle arrest leading to endopolyploidy only in monocentric and not (or to much lesser extent) in holocentric plants. While current microscope-dependent methods for holocentrism detection are unreliable for small and numerous chromosomes, which are common in holocentrics, FCM can use somatic nuclei. Thus, FCM may be a rapid and reliable method of high-throughput screening for holocentric candidates across plant phylogeny.