Species coherence in the face of karyotype diversification in holocentric organisms: the case of a cytogenetically variable sedge (Carex scoparia, Cyperaceae)

Founder events and subsequent genetic bottlenecks underlie karyotype evolution in the Ibero-North African endemic Carex helodes

BACKGROUND AND AIMS

Despite chromosomal evolution being one of the major drivers of diversification in plants, we do not yet have a clear view of how new chromosome rearrangements become fixed within populations, which is a crucial step forward for understanding chromosomal speciation.

METHODS

In this study, we test the role of genetic drift in the establishment of new chromosomal variants in the context of hybrid dysfunction models of chromosomal speciation. We genotyped 178 individuals from seven populations (plus 25 seeds from one population) across the geographical range of Carex helodes (Cyperaceae). We also characterized karyotype geographical patterns of the species across its distribution range. For one of the populations, we performed a detailed study of the fine-scale, local spatial distribution of its individuals and their genotypes and karyotypes.

KEY RESULTS

Synergistically, phylogeographical and karyotypic evidence revealed two main genetic groups: southwestern Iberian Peninsula vs. northwestern African populations; and within Europe our results suggest a west-to-east expansion with signals of genetic bottlenecks. Additionally, we inferred a pattern of descending dysploidy, plausibly as a result of a west-to-east process of post-glacial colonization in Europe.

CONCLUSIONS

Our results give experimental support to the role of geographical isolation, drift and inbreeding in the establishment of new karyotypes, which is key in the speciation models of hybrid dysfunction.

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.

The holocentric chromosome microevolution: From phylogeographic patterns to genomic associations with environmental gradients

Geographic isolation and chromosome evolution are two of the major drivers of diversification in eukaryotes in general, and specifically, in plants. On one hand, range shifts induced by Pleistocene glacial oscillations deeply shaped the evolutionary trajectories of species in the Northern Hemisphere. On the other hand, karyotype variability within species or species complexes may have adaptive potential as different karyotypes may represent different recombination rates and linkage groups that may be associated with locally adapted genes or supergenes. Organisms with holocentric chromosomes are ideal to study the link between local adaptation and chromosome evolution, due to their high cytogenetic variability, especially when it seems to be related to environmental variation. Here, we integrate the study of the phylogeography, chromosomal evolution and ecological requirements of a plant species complex distributed in the Western Euro-Mediterranean region (Carex gr. laevigata, Cyperaceae). We aim to clarify the relative influence of these factors on population differentiation and ultimately on speciation. We obtained a well-resolved RADseq phylogeny that sheds light on the phylogeographic patterns of molecular and chromosome number variation, which are compatible with south-to-north postglacial migration. In addition, landscape genomics analyses identified candidate loci for local adaptation, and also strong significant associations between the karyotype and the environment. We conclude that karyotype distribution in C. gr. laevigata has been constrained by both range shift dynamics and local adaptation. Our study demonstrates that chromosome evolution may be responsible, at least partially, for microevolutionary patterns of population differentiation and adaptation in Carex.

Genomic hotspots of chromosome rearrangements explain conserved synteny despite high rates of chromosome evolution in a holocentric lineage

Holocentric organisms, unlike typical monocentric organisms, have kinetochore activity distributed along almost the whole length of the chromosome. Because of this, chromosome rearrangements through fission and fusion are more likely to become fixed in holocentric species, which may account for the extraordinary rates of chromosome evolution that many holocentric lineages exhibit. Long blocks of genome synteny have been reported in animals with holocentric chromosomes despite high rates of chromosome rearrangements. Nothing is known from plants, however, despite the fact that holocentricity appears to have played a key role in the diversification of one of the largest angiosperm genera, Carex (Cyperaceae). In the current study, we compared genomes of Carex species and a distantly related Cyperaceae species to characterize conserved and rearranged genome regions. Our analyses span divergence times ranging between 2 and 50 million years. We also compared a C. scoparia chromosome-level genome assembly with a linkage map of the same species to study rearrangements at a population level and suppression of recombination patterns. We found longer genome synteny blocks than expected under a null model of random rearrangement breakpoints, even between very distantly related species. We also found repetitive DNA to be non-randomly associated with holocentromeres and rearranged regions of the genome. The evidence of conserved synteny in sedges despite high rates of chromosome fission and fusion suggests that conserved genomic hotspots of chromosome evolution related to repetitive DNA shape the evolution of recombination, gene order and crossability in sedges. This finding may help explain why sedges are able to maintain species cohesion even in the face of high interspecific chromosome rearrangements.

Chromosome-scale genome assemblies and annotations for Poales species Carex cristatella, Carex scoparia, Juncus effusus, and Juncus inflexus

The majority of sequenced genomes in the monocots are from species belonging to Poaceae, which include many commercially important crops. Here, we expand the number of sequenced genomes from the monocots to include the genomes of 4 related cyperids: Carex cristatella and Carex scoparia from Cyperaceae and Juncus effusus and Juncus inflexus from Juncaceae. The high-quality, chromosome-scale genome sequences from these 4 cyperids were assembled by combining whole-genome shotgun sequencing of Nanopore long reads, Illumina short reads, and Hi-C sequencing data. Some members of the Cyperaceae and Juncaceae are known to possess holocentric chromosomes. We examined the repeat landscapes in our sequenced genomes to search for potential repeats associated with centromeres. Several large satellite repeat families, comprising 3.2–9.5% of our sequenced genomes, showed dispersed distribution of large satellite repeat clusters across all Carex chromosomes, with few instances of these repeats clustering in the same chromosomal regions. In contrast, most large Juncus satellite repeats were clustered in a single location on each chromosome, with sporadic instances of large satellite repeats throughout the Juncus genomes. Recognizable transposable elements account for about 20% of each of the 4 genome assemblies, with the Carex genomes containing more DNA transposons than retrotransposons while the converse is true for the Juncus genomes. These genome sequences and annotations will facilitate better comparative analysis within monocots.

A holocentric twist to chromosomal speciation?

Chromosomal rearrangements trigger speciation by acting as barriers to gene flow. However, the underlying theory was developed with monocentric chromosomes in mind. Holocentric chromosomes, lacking a centromeric region, have repeatedly evolved and account for a significant fraction of extant biodiversity. Because chromosomal rearrangements may be more likely retained in holocentric species, holocentricity could provide a twist to chromosomal speciation. Here, we discuss how the abundance of chromosome-scale genomes, combined with novel analytical tools, offer the opportunity to assess the impacts of chromosomal rearrangements on rates of speciation by outlining a phylogenetic framework that aligns with the two major lines of chromosomal speciation theory. We further highlight how holocentric species could help to test for causal roles of chromosomal rearrangements in speciation.

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

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

Holocentric chromosomes

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

RAD-seq linkage mapping and patterns of segregation distortion in sedges: meiosis as a driver of karyotypic evolution in organisms with holocentric chromosomes

Meiotic drive, the class of meiotic mechanisms that drive unequal segregation of alleles among gametes, may be an important force in karyotype evolution. Its role in holocentric organisms, whose chromosomes lack localized centromeres, is poorly understood. We crossed two individuals of Carex scoparia (Cyperaceae) with different chromosome numbers (2n = 33^II  = 66 × 2n = 32^II  = 64) to obtain F1 individuals, which we then self-pollinated to obtain second-generation (F2) crosses. RAD-seq was performed for 191 individuals (including the parents, five F1 individuals and 184 F2 individuals). Our F2 linkage map based on stringent editing of the RAD-seq data set yielded 32 linkage groups. In the final map, 865 loci were located on a linkage map of 3966.99 cM (linkage groups ranged from 24.39 to 193.31 cM in length and contained 5-51 loci each). Three linkage groups exhibit more loci under segregation distortion than expected by chance; within linkage groups, loci exhibiting segregation distortion are clustered. This finding implicates meiotic drive in the segregation of chromosome variants, suggesting that selection of chromosome variants in meiosis may contribute to the establishment and fixation of chromosome variants in Carex, which is renowned for high chromosomal and species diversity. This is an important finding as previous studies demonstrate that chromosome divergence may play a key role in differentiation and speciation in Carex.

Two closely related species differ in their regional genetic differentiation despite admixing

Regional genetic differentiation within species is often addressed in evolutionary ecology and conservation biology. Here, we address regional differentiation in two closely related hybridizing taxa, the perennial sedges Carex flava and C. viridula and their hybrid C. × subviridula in 37 populations in the north and centre of their distribution range in Europe (Estonia, Lowland (<1000 m a.s.l.) and Highland Switzerland) using 10 putative microsatellite loci. We ask whether regional differentiation was larger in the less common taxon C. viridula or whether, possibly due to hybridization, it was similar between taxa. Our results showed similar, low to moderate genetic diversity for the three studied taxa. In total, we found 12 regional species-specific alleles. Analysis of molecular variance (AMOVA), STRUCTURE and multidimensional scaling analysis showed regional structure in genetic variation, where intraspecific differentiation between regions was lower for C. flava (AMOVA: 6.84 %) than for C. viridula (20.77 %) or C. × subviridula (18.27 %) populations. Hybrids differed from the parental taxa in the two regions where they occurred, i.e. in Estonia and Lowland Switzerland. We conclude that C. flava and C. viridula clearly differ from each other genetically, that there is pronounced regional differentiation and that, despite hybridization, this regional differentiation is more pronounced in the less common taxon, C. viridula. We encourage future studies on hybridizing taxa to work with plant populations from more than one region.

Do holocentric chromosomes represent an evolutionary advantage? A study of paired analyses of diversification rates of lineages with holocentric chromosomes and their monocentric closest relatives

Despite most of the cytogenetic research is focused on monocentric chromosomes, chromosomes with kinetochoric activity localized in a single centromere, several studies have been centered on holocentric chromosomes which have diffuse kinetochoric activity along the chromosomes. The eukaryotic organisms that present this type of chromosomes have been relatively understudied despite they constitute rather diversified species lineages. On the one hand, holocentric chromosomes may present intrinsic benefits (chromosome mutations such as fissions and fusions are potentially neutral in holocentrics). On the other hand, they present restrictions to the spatial separation of the functions of recombination and segregation during meiotic divisions (functions that may interfere), separation that is found in monocentric chromosomes. In this study, we compare the diversification rates of all known holocentric lineages in animals and plants with their most related monocentric lineages in order to elucidate whether holocentric chromosomes constitute an evolutionary advantage in terms of diversification and species richness. The results showed that null hypothesis of equal mean diversification rates cannot be rejected, leading us to surmise that shifts in diversification rates between holocentric and monocentric lineages might be due to other factors, such as the idiosyncrasy of each lineage or the interplay of evolutionary selections with the benefits of having either monocentric or holocentric chromosomes.

Chromosomal rearrangements in holocentric organisms lead to reproductive isolation by hybrid dysfunction: The correlation between karyotype rearrangements and germination rates in sedges

PREMISE OF THE STUDY

Understanding the drivers of speciation is a central task of evolutionary biology. Chromosomal rearrangements are known to play an important role in species diversification, but the role of rearrangements of holocentric chromosomes-chromosomes without localized centromeres-is poorly understood.

METHODS

We made numerous artificial crosses between Carex scoparia individuals of different diploid chromosome numbers and, for comparison, between individuals of the same chromosome number. We studied chromosome pairing and chromosomal rearrangements in the F1 individuals using light microscopy. We then estimated germination rates as a function of geographic distance, genetic distance, chromosome number differences in parents, and pairing irregularities in F1 individuals, using generalized least squares to fit alternative regression models.

KEY RESULTS

The most informative predictors of germination rates in the F1 generation are chromosome number differences and minimum number of chromosome pairing irregularities in the F1 individuals. Genetic and geographic distances between parents are not significant predictors.

CONCLUSIONS

Holocentric chromosomal rearrangements play an important role in postzygotic reproductive isolation in Carex through F1 hybrid inviability and sterility. Hybrid dysfunction seems to be a suitable model for chromosomal speciation when there are several chromosomal rearrangements between parents. However, we have not tested the hypothesis that genome rearrangements may also play an important role in suppressing recombination between cytogenetically divergent populations.

Genome size stability despite high chromosome number variation in Carex gr. laevigata

PREMISE OF THE STUDY

In organisms with holocentric chromosomes like Carex species, chromosome number evolution has been hypothesized to be a result of fission, fusion, and/or translocation events. Negative, positive, or the absence of correlations have been found between chromosome number and genome size in Carex.

METHODS

Using the inferred diploid chromosome number and 80 genome size measurements from 26 individuals and 20 populations of Carex gr. laevigata, we tested the null hypothesis of chromosome number evolution by duplication and deletion of whole chromosomes.

KEY RESULTS

Our results show a significant positive correlation between genome size and chromosome number, but the slope of such correlation supports the hypothesis of proliferation and removal of repetitive DNA fragments to explain genome size variation rather than duplication and deletion of whole chromosomes.

CONCLUSIONS

Our results refine the theory of the holokinetic drive: this mechanism is proposed to facilitate repetitive DNA removal (or any segmental deletion) when smaller homologous chromosomes are preferentially inherited, or repetitive DNA proliferation (or any segmental duplication) when larger homologs are preferred. This study sheds light on how karyotype evolution plays an important role in the diversification of the species of the genus Carex.