Holokinetic drive: centromere drive in chromosomes without centromeres

Do chromosome rearrangements fix by genetic drift or natural selection? Insights from Brenthis butterflies

Large-scale chromosome rearrangements, such as fissions and fusions, are a common feature of eukaryote evolution. They can have considerable influence on the evolution of populations, yet it remains unclear exactly how rearrangements become established and eventually fix. Rearrangements could fix by genetic drift if they are weakly deleterious or neutral, or they may instead be favoured by positive natural selection. Here, we compare genome assemblies of three closely related Brenthis butterfly species and characterize a complex history of fission and fusion rearrangements. An inferred demographic history of these species suggests that rearrangements became fixed in populations with large long-term effective size (N e), consistent with rearrangements being selectively neutral or only very weakly underdominant. Using a recently developed analytic framework for characterizing hard selective sweeps, we find that chromosome fusions are not enriched for evidence of past sweeps compared to other regions of the genome. Nonetheless, we do infer a strong and recent selective sweep around one chromosome fusion in the B. daphne genome. Our results suggest that rearrangements in these species likely have weak absolute fitness effects and fix by genetic drift. However, one putative selective sweep raises the possibility that natural selection may sometimes play a role in the fixation of chromosome fusions.

Meiotic drive against chromosome fusions in butterfly hybrids

Species frequently differ in the number and structure of chromosomes they harbor, but individuals that are heterozygous for chromosomal rearrangements may suffer from reduced fitness. Chromosomal rearrangements like fissions and fusions can hence serve as a mechanism for speciation between incipient lineages, but their evolution poses a paradox. How can rearrangements get fixed between populations if heterozygotes have reduced fitness? One solution is that this process predominantly occurs in small and isolated populations, where genetic drift can override natural selection. However, fixation is also more likely if a novel rearrangement is favored by a transmission bias, such as meiotic drive. Here, we investigate chromosomal transmission distortion in hybrids between two wood white (Leptidea sinapis) butterfly populations with extensive karyotype differences. Using data from two different crossing experiments, we uncover that there is a transmission bias favoring the ancestral chromosomal state for derived fusions, a result that shows that chromosome fusions actually can fix in populations despite being counteracted by meiotic drive. This means that meiotic drive not only can promote runaway chromosome number evolution and speciation, but also that it can be a conservative force acting against karyotypic change and the evolution of reproductive isolation. Based on our results, we suggest a mechanistic model for why chromosome fusion mutations may be opposed by meiotic drive and discuss factors contributing to karyotype evolution in Lepidoptera.

Chromosome size matters: genome evolution in the cyperid clade

BACKGROUND AND AIMS

While variation in genome size and chromosome numbers and their consequences are often investigated in plants, the biological relevance of variation in chromosome size remains poorly known. Here, we examine genome and mean chromosome size in the cyperid clade (families Cyperaceae, Juncaceae and Thurniaceae), which is the largest vascular plant lineage with predominantly holocentric chromosomes.

METHODS

We measured genome size in 436 species of cyperids using flow cytometry, and augment these data with previously published datasets. We then separately compared genome and mean chromosome sizes (2C/2n) amongst the major lineages of cyperids and analysed how these two genomic traits are associated with various environmental factors using phylogenetically informed methods.

KEY RESULTS

We show that cyperids have the smallest mean chromosome sizes recorded in seed plants, with a large divergence between the smallest and largest values. We found that cyperid species with smaller chromosomes have larger geographical distributions and that there is a strong inverse association between mean chromosome size and number across this lineage.

CONCLUSIONS

The distinct patterns in genome size and mean chromosome size across the cyperids might be explained by holokinetic drive. The numerous small chromosomes might function to increase genetic diversity in this lineage where crossovers are limited during meiosis.

Comparative Analysis of Transposable Elements and the Identification of Candidate Centromeric Elements in the Prunus Subgenus Cerasus and Its Relatives

The subgenus Cerasus and its relatives include many crucial economic drupe fruits and ornamental plants. Repetitive elements make up a large part of complex genomes, and some of them play an important role in gene regulation that can affect phenotypic variation. However, the variation in their genomes remains poorly understood. This work conducted a comprehensive repetitive sequence identification across the draft genomes of eight taxa of the genus Prunus, including four of the Prunus subgenus Cerasus (Prunus pseudocerasus, P. avium, P. yedoensis and P. × yedoensis) as well as congeneric species (Prunus salicina, P. armeniaca, P. dulcis and P. persica). Annotation results showed high proportions of transposable elements in their genomes, ranging from 52.28% (P. armeniaca) to 61.86% (P. pseudocerasus). The most notable differences in the contents of long terminal repeat retrotransposons (LTR-RTs) and tandem repeats (TRs) were confirmed with de novo identification based on the structure of each genome, which significantly contributed to their genome size variation, especially in P. avium and P.salicina. Sequence comparisons showed many similar LTR-RTs closely related to their phylogenetic relationships, and a highly similar monomer unit of the TR sequence was conserved among species. Additionally, the predicted centromere-associated sequence was located in centromeric regions with FISH in the 12 taxa of Prunus. It presented significantly different signal intensities, even within the diverse interindividual phenotypes for Prunus tomentosa. This study provides insight into the LTR-RT and TR variation within Prunus and increases our knowledge about its role in genome evolution.

Centromere size scales with genome size across Eukaryotes

Previous studies on grass species suggested that the total centromere size (sum of all centromere sizes in a cell) may be determined by the genome size, possibly because stable scaling is important for proper cell division. However, it is unclear whether this relationship is universal. Here we analyze the total centromere size using the CenH3-immunofluorescence area as a proxy in 130 taxa including plants, animals, fungi, and protists. We verified the reliability of our methodological approach by comparing our measurements with available ChIP-seq-based measurements of the size of CenH3-binding domains. Data based on these two independent methods showed the same positive relationship between the total centromere size and genome size. Our results demonstrate that the genome size is a strong predictor (R-squared = 0.964) of the total centromere size universally across Eukaryotes. We also show that this relationship is independent of phylogenetic relatedness and centromere type (monocentric, metapolycentric, and holocentric), implying a common mechanism maintaining stable total centromere size in Eukaryotes.

Holocentric Chromosomes Probably Do Not Prevent Centromere Drive in Cyperaceae

Centromere drive model describes an evolutionary process initiated by centromeric repeats expansion, which leads to the recruitment of excess kinetochore proteins and consequent preferential segregation of an expanded centromere to the egg during female asymmetric meiosis. In response to these selfish centromeres, the histone protein CenH3, which recruits kinetochore components, adaptively evolves to restore chromosomal parity and counter the detrimental effects of centromere drive. Holocentric chromosomes, whose kinetochores are assembled along entire chromosomes, have been hypothesized to prevent expanded centromeres from acquiring a selective advantage and initiating centromere drive. In such a case, CenH3 would be subjected to less frequent or no adaptive evolution. Using codon substitution models, we analyzed 36 CenH3 sequences from 35 species of the holocentric family Cyperaceae. We found 10 positively selected codons in the CenH3 gene [six codons in the N-terminus and four in the histone fold domain (HFD)] and six branches of its phylogeny along which the positive selection occurred. One of the positively selected codons was found in the centromere targeting domain (CATD) that directly interacts with DNA and its mutations may be important in centromere drive suppression. The frequency of these positive selection events was comparable to the frequency of positive selection in monocentric clades with asymmetric female meiosis. Taken together, these results suggest that preventing centromere drive is not the primary adaptive role of holocentric chromosomes, and their ability to suppress it likely depends on their kinetochore structure in meiosis.

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

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

Chromosome numbers of Carex (Cyperaceae) and their taxonomic implications

Counting chromosomes is the first step towards a better understanding of the karyotype evolution and the role of chromosome evolution in species diversification within Carex; however, the chromosome count is not known yet for numerous sedges. In this paper chromosome counts were performed for 23 Carex taxa from Armenia, Austria, the Czech Republic, and Poland. Chromosome numbers were determined for the first time in three species (Carex cilicica, 2n = 54; C. phyllostachys, 2n = 56; C. randalpina, 2n = 78), two subspecies (C. muricata subsp. ashokae, 2n = 58; C. nigra subsp. transcaucasica, 2n = 84) and two hybrids (C. ×decolorans, 2n = 74; C. ×walasii, 2n = 108). Among the taxa whose number of chromosomes had been known before, the largest difference was found in C. hartmaniorum (here 2n = 52) and C. aterrima subsp. medwedewii (here 2n = 52). A difference in the chromosome count was demonstrated for C. cilicica (2n = 54) versus the species of the section Aulocystis (2n = 30 to 40) and for C. tomentosa (2n = 48) versus the species of the section Acrocystis (2n = 18 to 38). The results of this study indicate that the position of C. cilicica in Aulocystis section may raise doubts. Attention was paid to the relationship between C. phyllostachys and taxa of the subgenus Carex section Gynobasidae.

Holocentric Karyotype Evolution in Rhynchospora Is Marked by Intense Numerical, Structural, and Genome Size Changes

Cyperaceae is a family of Monocotyledons comprised of species with holocentric chromosomes that are associated with intense dysploidy and polyploidy events. Within this family the genus Rhynchospora has recently become the focus of several studies that characterize the organization of the holocentric karyotype and genome structures. To broaden our understanding of genome evolution in this genus, representatives of Rhynchospora were studied to contrast chromosome features, C-CMA/DAPI band distribution and genome sizes. Here, we carried out a comparative analysis for 35 taxa of Rhynchospora, and generated new genome size estimates for 20 taxa. The DNA 2C-values varied up to 22-fold, from 2C = 0.51 pg to 11.32 pg, and chromosome numbers ranged from 2n = 4 to 61. At least 37% of our sampling exhibited 2n different from the basic number x = 5, and chromosome rearrangements were also observed. A large variation in C-CMA/DAPI band accumulation and distribution was observed as well. We show that genome variation in Rhynchospora is much larger than previously reported. Phylogenetic analysis showed that most taxa were grouped in clades corresponding to previously described taxonomic sections. Basic chromosome numbers are the same within every section, however, changes appeared in all the clades. Ancestral chromosome number reconstruction revealed n = 5 as the most likely ancestral complements, but n = 10 appears as a new possibility. Chromosome evolution models point to polyploidy as the major driver of chromosome evolution in Rhynchospora, followed by dysploidy. A negative correlation between chromosome size and diploid number open the discussion for holokinetic drive-based genome evolution. This study explores relationships between karyotype differentiation and genome size variation in Rhynchospora, and contrasts it against the phylogeny of this holocentric group.

Insights into the karyotype and genome evolution of haplogyne spiders indicate a polyploid origin of lineage with holokinetic chromosomes

Spiders are an ancient and extremely diverse animal order. They show a considerable diversity of genome sizes, karyotypes and sex chromosomes, which makes them promising models to analyse the evolution of these traits. Our study is focused on the evolution of the genome and chromosomes in haplogyne spiders with holokinetic chromosomes. Although holokinetic chromosomes in spiders were discovered a long time ago, information on their distribution and evolution in these arthropods is very limited. Here we show that holokinetic chromosomes are an autapomorphy of the superfamily Dysderoidea. According to our hypothesis, the karyotype of ancestral Dysderoidea comprised three autosome pairs and a single X chromosome. The subsequent evolution has frequently included inverted meiosis of the sex chromosome and an increase of 2n. We demonstrate that caponiids, a sister clade to Dysderoidea, have enormous genomes and high diploid and sex chromosome numbers. This pattern suggests a polyploid event in the ancestors of caponiids. Holokinetic chromosomes could have arisen by subsequent multiple chromosome fusions and a considerable reduction of the genome size. We propose that spider sex chromosomes probably do not pose a major barrier to polyploidy due to specific mechanisms that promote the integration of sex chromosome copies into the genome.

Unexpected patterns of segregation distortion at a selfish supergene in the fire ant Solenopsis invicta

BACKGROUND

The Sb supergene in the fire ant Solenopsis invicta determines the form of colony social organization, with colonies whose inhabitants bear the element containing multiple reproductive queens and colonies lacking it containing only a single queen. Several features of this supergene - including suppressed recombination, presence of deleterious mutations, association with a large centromere, and "green-beard" behavior - suggest that it may be a selfish genetic element that engages in transmission ratio distortion (TRD), defined as significant departures in progeny allele frequencies from Mendelian inheritance ratios. We tested this possibility by surveying segregation ratios in embryo progenies of 101 queens of the "polygyne" social form (3512 embryos) using three supergene-linked markers and twelve markers outside the supergene.

RESULTS

Significant departures from Mendelian ratios were observed at the supergene loci in 3-5 times more progenies than expected in the absence of TRD and than found, on average, among non-supergene loci. Also, supergene loci displayed the greatest mean deviations from Mendelian ratios among all study loci, although these typically were modest. A surprising feature of the observed inter-progeny variation in TRD was that significant deviations involved not only excesses of supergene alleles but also similarly frequent excesses of the alternate alleles on the homologous chromosome. As expected given the common occurrence of such "drive reversal" in this system, alleles associated with the supergene gain no consistent transmission advantage over their alternate alleles at the population level. Finally, we observed low levels of recombination and incomplete gametic disequilibrium across the supergene, including between adjacent markers within a single inversion.

CONCLUSIONS

Our data confirm the prediction that the Sb supergene is a selfish genetic element capable of biasing its own transmission during reproduction, yet counterselection for suppressor loci evidently has produced an evolutionary stalemate in TRD between the variant homologous haplotypes on the "social chromosome". Evidence implicates prezygotic segregation distortion as responsible for the TRD we document, with "true" meiotic drive the most likely mechanism. Low levels of recombination and incomplete gametic disequilibrium across the supergene suggest that selection does not preserve a single uniform supergene haplotype responsible for inducing polygyny.

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.

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.

CenH3 evolution reflects meiotic symmetry as predicted by the centromere drive model

The centromere drive model explaining rapid evolution of eukaryotic centromeres predicts higher frequency of positive selection acting on centromeric histone H3 (CenH3) in clades with asymmetric meiosis compared to the clades with only symmetric meiosis. However, despite the impression one might get from the literature, this key prediction of the centromere drive model has not only never been confirmed, but it has never been tested, because all the previous studies dealt only with the presence or absence instead of the frequency of positive selection. To provide evidence for or against different frequencies of positively selected CenH3 in asymmetrics and symmetrics, we have inferred the selective pressures acting on CenH3 in seventeen eukaryotic clades, including plants, animals, fungi, ciliates and apicomplexa, using codon-substitution models, and compared the inferred frequencies between asymmetrics and symmetrics in a quantitative manner. We have found that CenH3 has been evolving adaptively much more frequently in clades with asymmetric meiosis compared with clades displaying only symmetric meiosis which confirms the prediction of centromere drive model. Our findings indicate that the evolution of asymmetric meiosis required CenH3 to evolve adaptively more often to counterbalance the negative consequences of centromere drive.

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.

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.

The blue butterfly Polyommatus (Plebicula) atlanticus (Lepidoptera, Lycaenidae) holds the record of the highest number of chromosomes in the non-polyploid eukaryotic organisms

The blue butterfly species Polyommatus (Plebicula) atlanticus (Elwes, 1906) (Lepidoptera, Lycaenidae) is known to have a very high haploid number of chromosomes (n= circa 223). However, this approximate count made by Hugo de Lesse 45 years ago was based on analysis of a single meiotic I metaphase plate, not confirmed by study of diploid chromosome set and not documented by microphotographs. Here I demonstrate that (1) Polyommatus atlanticus is a diploid (non-polyploid) species, (2) its meiotic I chromosome complement includes at least 224-226 countable chromosome bodies, and (3) all (or nearly all) chromosome elements in meiotic I karyotype are represented by bivalents. I also provide the first data on the diploid karyotype and estimate the diploid chromosome number as 2n=ca448-452. Thus, Polyommatus atlanticus is confirmed to possess the highest chromosome number among all the non-polyploid eukaryotic organisms.

Dynamic karyotype evolution and unique sex determination systems in Leptidea wood white butterflies

Background

Chromosomal rearrangements have the potential to limit the rate and pattern of gene flow within and between species and thus play a direct role in promoting and maintaining speciation. Wood white butterflies of the genus Leptidea are excellent models to study the role of chromosome rearrangements in speciation because they show karyotype variability not only among but also within species. In this work, we investigated genome architecture of three cryptic Leptidea species (L. juvernica, L. sinapis and L. reali) by standard and molecular cytogenetic techniques in order to reveal causes of the karyotype variability.

Results

Chromosome numbers ranged from 2n = 85 to 91 in L. juvernica and 2n = 69 to 73 in L. sinapis (both from Czech populations) to 2n = 51 to 55 in L. reali (Spanish population). We observed significant differences in chromosome numbers and localization of cytogenetic markers (rDNA and H3 histone genes) within the offspring of individual females. Using FISH with the (TTAGG)_n telomeric probe we also documented the presence of multiple chromosome fusions and/or fissions and other complex rearrangements. Thus, the intraspecific karyotype variability is likely due to irregular chromosome segregation of multivalent meiotic configurations. The analysis of female meiotic chromosomes by GISH and CGH revealed multiple sex chromosomes: W₁W₂W₃Z₁Z₂Z₃Z₄ in L. juvernica, W₁W₂W₃Z₁Z₂Z₃ in L. sinapis and W₁W₂W₃W₄Z₁Z₂Z₃Z₄ in L. reali.

Conclusions

Our results suggest a dynamic karyotype evolution and point to the role of chromosomal rearrangements in the speciation of Leptidea butterflies. Moreover, our study revealed a curious sex determination system with 3–4 W and 3–4 Z chromosomes, which is unique in the Lepidoptera and which could also have played a role in the speciation process of the three Leptidea species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0375-4) contains supplementary material, which is available to authorized users.