Evolutionary Mechanisms of Varying Chromosome Numbers in the Radiation of Erebia Butterflies

A macroevolutionary role for chromosomal fusion and fission in Erebia butterflies

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

Investigating the evolutionary dynamics of diploid number variation in Ctenomys (Ctenomyidae, Rodentia)

Contrary to predictions from classical hybrid sterility models of chromosomal speciation, some organisms display high rates of karyotype variation. Ctenomys are the current mammals with the greatest interspecific and intraspecific chromosomal variation. A large number of species have been studied cytogenetically. The diploid numbers of chromosomes range from 2n = 10 to 2n = 70. Here, we analyzed karyotype evolution in Ctenomys using comparative phylogenetic methods. We found a strong phylogenetic signal with chromosome number. This refutes the chromosomal megaevolution model, which proposes the independent accumulation of multiple chromosomal rearrangements in each closely related species. We found that Brownian motion (BM) described the observed characteristic changes more thoroughly than the Ornstein-Uhlenbeck and Early-Burst models. This suggests that the evolution of chromosome numbers occurs by a random walk along phylogenetic clades. However, our data indicate that the BM model alone does not fully characterize the chromosomal evolution of Ctenomys.

Two chromosome-level genome assemblies of Rhodiola shed new light on genome evolution in rapid radiation and evolution of the biosynthetic pathway of salidroside

Rhodiola L. is a genus that has undergone rapid radiation in the mid-Miocene and may represent a typic case of adaptive radiation. Many species of Rhodiola have also been widely used as an important adaptogen in traditional medicines for centuries. However, a lack of high-quality chromosome-level genomes hinders in-depth study of its evolution and biosynthetic pathway of secondary metabolites. Here, we assembled two chromosome-level genomes for two Rhodiola species with different chromosome number and sexual system. The assembled genome size of R. chrysanthemifolia (2n = 14; hermaphrodite) and R. kirilowii (2n = 22; dioecious) were of 402.67 and 653.62 Mb, respectively, with approximately 57.60% and 69.22% of transposable elements (TEs). The size difference between the two genomes was mostly due to proliferation of long terminal repeat-retrotransposons (LTR-RTs) in the R. kirilowii genome. Comparative genomic analysis revealed possible gene families responsible for high-altitude adaptation of Rhodiola, including a homolog of plant cysteine oxidase 2 gene of Arabidopsis thaliana (AtPCO2), which is part of the core molecular reaction to hypoxia and contributes to the stability of Group VII ethylene response factors (ERF-VII). We found extensive chromosome fusion/fission events and structural variations between the two genomes, which might have facilitated the initial rapid radiation of Rhodiola. We also identified candidate genes in the biosynthetic pathway of salidroside. Overall, our results provide important insights into genome evolution in plant rapid radiations, and possible roles of chromosome fusion/fission and structure variation played in rapid speciation.

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.

Genome assembly of 3 Amazonian Morpho butterfly species reveals Z-chromosome rearrangements between closely related species living in sympatry

The genomic processes enabling speciation and species coexistence in sympatry are still largely unknown. Here we describe the whole-genome sequencing and assembly of 3 closely related species from the butterfly genus Morpho: Morpho achilles (Linnaeus, 1758), Morpho helenor (Cramer, 1776), and Morpho deidamia (Höbner, 1819). These large blue butterflies are emblematic species of the Amazonian rainforest. They live in sympatry in a wide range of their geographical distribution and display parallel diversification of dorsal wing color pattern, suggesting local mimicry. By sequencing, assembling, and annotating their genomes, we aim at uncovering prezygotic barriers preventing gene flow between these sympatric species. We found a genome size of  480 Mb for the 3 species and a chromosomal number ranging from 2n = 54 for M. deidamia to 2n = 56 for M. achilles and M. helenor. We also detected inversions on the sex chromosome Z that were differentially fixed between species, suggesting that chromosomal rearrangements may contribute to their reproductive isolation. The annotation of their genomes allowed us to recover in each species at least 12,000 protein-coding genes and to discover duplications of genes potentially involved in prezygotic isolation like genes controlling color discrimination (L-opsin). Altogether, the assembly and the annotation of these 3 new reference genomes open new research avenues into the genomic architecture of speciation and reinforcement in sympatry, establishing Morpho butterflies as a new eco-evolutionary model.

Highly divergent karyotypes and barcoding of the East African genus Gonatoxia Karsch (Orthoptera: Phaneropterinae)

East Africa is a hotspot of biodiversity of many orthopteran taxa, including bushcrickets. Gonatoxia Karsch, 1889 species are fully alate Phaneropterinae, which are perfectly adapted to the foliage of forests. We examined five species using combined cytogenetic and molecular data to determine the inter- and intraspecific genetic diversity. The variation in the diploid number of chromosomes in males ranged from 2n = 28 + X0 and 26 + X0 to 2n = 6 + X0. Fluorescence in situ hybridization showed from one to many 18S rDNA loci as well as interstitial sequences, especially in G. helleri. 18S rDNA loci coincided with active NOR and C-banding patterns. The isolation of populations of the species explains differences in the number of chromosomes (G. maculata), chromosomal polymorphism and chromosomal heterozygosity (G. helleri). Our molecular phylogeny based on the COI locus supported the monophyly of the genus Gonatoxia and separateness of the five examined species in accordance with their morphological features and chromosome numbers as well as the species' distribution.

Dynamics of sex chromosome evolution in a rapid radiation of cichlid fishes

Sex is a fundamental trait determined by environmental and/or genetic factors, including sex chromosomes. Sex chromosomes are studied in species scattered across the tree of life, yet little is known about tempo and mode of sex chromosome evolution among closely related species. Here, we examine sex chromosome evolution in the adaptive radiation of cichlid fishes in Lake Tanganyika. Through the analysis of male and female genomes from 244 cichlid taxa (189 described species with 5 represented with two local variants/populations; 50 undescribed species) and of 396 multitissue transcriptomes from 66 taxa, we identify signatures of sex chromosomes in 79 taxa, involving 12 linkage groups. We find that Tanganyikan cichlids have the highest rates of sex chromosome turnover and heterogamety transitions known to date. We show that sex chromosome recruitment is not at random. Moreover convergently emerged sex chromosomes in cichlids support the “limited options” hypothesis of sex chromosome evolution.

Diversification of the Balloon bushcrickets (Orthoptera, Hexacentrinae, Aerotegmina) in the East African mountains

East African mountains constitute a network of isolated habitat islands among dry savannah and are thus ideal for studying species diversification processes. This study elucidated the phylogenetic and phylogeographic relationships of all bushcricket species comprising the genus Aerotegmina. Our analysis indicated that large-scale climatic and topographic processes in Africa are likely to have driven speciation in this group, and revealed the cytogenetic traits of the species. Molecular phylogeny supported the monophyly of Aerotegmina and showed that the genus probably originated in the old Eastern Arc Mountains of Tanzania and Kenya. Two lineages were distinguished: small- and large-sized species with geographically distinct habitats. The underlying processes are thought to be eight dispersals, ten vicariance events, and one extinction event linked to repeated fragmentation of the African rainforest. Those processes, in conjunction with habitat change, probably also led to the spatial separation of the species into a northern clade with a diploid number of chromosomes 2n = 32 + X0 or 2n = 30 + neo-XY and a southern clade with a reduced number of chromosomes (2n = 28 + X0 or 24 + neo-X₁X₂Y). Karyotype analysis suggests that Aerotegmina is currently in the process of speciation.

Genome size versus geographic range size in birds

Why do some species occur in small, restricted areas, while others are distributed globally? Environmental heterogeneity increases with area and so does the number of species. Hence, diverse biotic and abiotic conditions across large ranges may lead to specific adaptations that are often linked to a species' genome size and chromosome number. Therefore, a positive association between genome size and geographic range is anticipated. Moreover, high cognitive ability in organisms would be favored by natural selection to cope with the dynamic conditions within large geographic ranges. Here, we tested these hypotheses in birds-the most mobile terrestrial vertebrates-and accounted for the effects of various confounding variables, such as body mass, relative brain mass, and geographic latitude. Using phylogenetic generalized least squares and phylogenetic confirmatory path analysis, we demonstrated that range size is positively associated with bird genome size but probably not with chromosome number. Moreover, relative brain mass had no effect on range size, whereas body mass had a possible weak and negative effect, and range size was larger at higher geographic latitudes. However, our models did not fully explain the overall variation in range size. Hence, natural selection may impose larger genomes in birds with larger geographic ranges, although there may be additional explanations for this phenomenon.

Chromosome number evolves at equal rates in holocentric and monocentric clades

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

Secondary contact zones of closely-related Erebia butterflies overlap with narrow phenotypic and parasitic clines

Zones of secondary contact between closely related taxa are a common legacy of the Quaternary ice ages. Despite their abundance, the factors that keep species apart and prevent hybridization are often unknown. Here, we study a very narrow contact zone between three closely related butterfly species of the Erebia tyndarus species complex. Using genomic data, we first determined whether gene flow occurs and then assessed whether it might be hampered by differences in chromosome number between some species. We found interspecific gene flow between sibling species that differ in karyotype by one chromosome. Conversely, only F1 hybrids occurred between two species that have the same karyotype, forming a steep genomic cline. In a second step, we fitted clines to phenotypic, ecological and parasitic data to identify the factors associated with the genetic cline. We found clines for phenotypic data and the prevalence of the endosymbiont parasite Wolbachia to overlap with the genetic cline, suggesting that they might be drivers for separating the two species. Overall, our results highlight that some gene flow is possible between closely related species despite different chromosome numbers, but that other barriers restrict such gene flow.

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'.

Rapid chromosomal evolution in the bush-cricket Gonatoxia helleri Hemp, 2016 (Orthoptera, Phaneropterinae)

Gonatoxia helleri Hemp, 2016 is one of the most widespread bush-crickets of the genus Gonatoxia Karsch, 1889 in East Africa. This species with seven large chromosomes (2n♂ = 7) differs from other representatives of the genus Gonatoxia drastically by its reduced chromosome number, the asymmetrical karyotype including karyomorphs rarely found in tettigoniids, as well as in irregularities in the course of meiosis. To better understand the origin of such an exceptional karyotype, chromosomes of 29 specimens from four populations/localities were studied using classical techniques, such as C-banding, silver impregnation, fluorochrome double staining and fluorescence in situ hybridization (FISH) technique with 18S rDNA and (TTAGG) _n telomeric probes. FISH showed many 18S rDNA loci as well as interstitial telomeric sequences, where chromosome morphology varied in these components in terms of quantity and distribution. The 18S rDNA loci coincided with active NORs and C-banding patterns. We suggest that a combination of Robertsonian rearrangements and/or multiple common tandem fusions involving the same chromosomes contributed to the formation of this karyotype/karyomorphs. The results are the first step towards a better understanding of chromosomal reorganization and evolution within the genus Gonatoxia. Low chromosome number, together with the incidence of chromosomal polymorphism that is higher in G. helleri than previously reported in bush-crickets, implies that this species can be a valuable new model for cytogenetic and speciation studies. Our findings suggest that chromosomal translocations lead to diversification and speciation in this species and could be the driving force of adaptive radiation.

Postglacial ecotype formation under outcrossing and self-fertilization in Arabidopsis lyrata

The formation of ecotypes has been invoked as an important driver of postglacial biodiversity, because many species colonized heterogeneous habitats and experienced divergent selection. Ecotype formation has been predominantly studied in outcrossing taxa, while far less attention has been paid to the implications of mating system shifts. Here, we addressed whether substrate-related ecotypes exist in selfing and outcrossing populations of Arabidopsis lyrata subsp. lyrata and whether the genomic footprint differs between mating systems. The North American subspecies colonized both rocky and sandy habitats during postglacial range expansion and shifted the mating system from predominantly outcrossing to predominantly selfing in a number of regions. We performed an association study on pooled whole-genome sequence data of 20 selfing or outcrossing populations, which suggested genes involved in adaptation to substrate. Motivated by enriched gene ontology terms, we compared root growth between plants from the two substrates in a common environment and found that plants originating from sand grew roots faster and produced more side roots, independent of mating system. Furthermore, single nucleotide polymorphisms associated with substrate-related ecotypes were more clustered among selfing populations. Our study provides evidence for substrate-related ecotypes in A. lyrata and divergence in the genomic footprint between mating systems. The latter is the likely result of selfing populations having experienced divergent selection on larger genomic regions due to higher genome-wide linkage disequilibrium.