Holocentric chromosomes

Cytogenetics of insects in the era of chromosome-level genome assemblies

Over the past few years, a revolution has occurred in cytogenetics, driven by the emergence and spread of methods for obtaining high-quality chromosome-level genome assemblies. In fact, this has led to a new tool for studying chromosomes and chromosomal rearrangements, and this tool is thousands of times more powerful than light microscopy. This tool has revolutionized the cytogenetics of many groups of insects for which previously karyotype information, if available at all, was limited to the chromosome number. Even more impressive are the achievements of the genomic approach for studying the general patterns of chromosome organization and evolution in insects. Thus, it has been shown that rapid transformations of chromosomal numbers, which are often found in the order Lepidoptera, are most often carried out in the most parsimonious way, as a result of simple fusions and fissions of chromosomes. It has been established that these fusions and fissions are not random and occur independently in different phylogenetic lineages due to the reuse of the same ancestral chromosomal breakpoints. It has been shown that the tendency for chromosome fissions is correlated with the presence in chromosomes of the so-called interstitial telomeres, i. e. telomere-like structures located not at the ends of chromosomes, but inside them. It has been revealed that, in most insects, telomeric DNA is not just a set of short repeats, but a very long sequence consisting of (TTAGG)n (or other telomeric motifs), regularly and specifically interrupted by retrotransposons, and the telomeric motifs are diverse in terms of their length and nucleotide composition. The number of high-quality chromosome-level genome assemblies available for insects in the GenBank database is growing exponentially and now exceeds a thousand species. Therefore, the exceptional prospects for using genomic data for karyotype analysis are beyond doubt.

Towards identification of a holocentromere marker in the lepidopteran model Spodoptera frugiperda

Some insects have holocentric chromosomes, with multiple kinetochores rather than a single centromere. They also lack the CENP-A and CENP-C proteins, suggesting a kinetochore assembly process different from that of monocentric chromosomes. The homolog of CENP-T was recently shown to bind silent chromatin and to play a key role in kinetochore assembly in Bombyx mori, but its role in other insects with holocentric chromosomes is unknown. We identified kinetochore genes and analyzed their expression in Spodoptera frugiperda. We silenced the kinetochore genes cenp-L, cenp-S, cenp-X and ndc80 and searched for chromosome segregation defects in Sf9 cells. All kinetochore genes except cenp-S were more strongly expressed in gonadal than in somatic tissues. Immunofluorescence microscopy and RT-qPCR demonstrated the effective silencing of the target genes by transfection with dsRNA. In Sf9 cells depleted of CENP-L and NDC80, immunofluorescence microscopy revealed increases in mitotic index and in the proportion of cells with unaligned chromosomes or multipolar spindles. The depletion of CENP-S and CENP-X had no effect on mitotic index and no division defects were observed. This suggests that CENP-L and NDC80 play key roles in chromosome segregation, whereas the functions of CENP-S and CENP-X remain unknown. We have begun to characterize the kinetochore proteins (CENP-L, CENP-S, CENP-X, NDC80), a prerequisite for holocentromere identification in S. frugiperda. This study also provides the first information about the role, in Lepidoptera, of CENP-L, a protein essential to the structure of the constitutive centromere-associated network in species with monocentric chromosomes.

Genome Assembly of a Living Fossil, the Atlantic Horseshoe Crab Limulus polyphemus, Reveals Lineage-Specific Whole-Genome Duplications, Transposable Element-Based Centromeres, and a ZW Sex Chromosome System

Horseshoe crabs, considered living fossils with a stable morphotype spanning ∼445 million years, are evolutionarily, ecologically, and biomedically important species experiencing rapid population decline. Of the four extant species of horseshoe crabs, the Atlantic horseshoe crab, Limulus polyphemus, has become an essential component of the modern medicine toolkit. Here, we present the first chromosome-level genome assembly, and the most contiguous and complete assembly to date, for L. polyphemus using nanopore long-read sequencing and chromatin conformation analysis. We find support for three horseshoe crab-specific whole-genome duplications, but none shared with Arachnopulmonata (spiders and scorpions). Moreover, we discovered tandem duplicates of endotoxin detection pathway components Factors C and G, identify candidate centromeres consisting of Gypsy retroelements, and classify the ZW sex chromosome system for this species and a sister taxon, Carcinoscorpius rotundicauda. Finally, we revealed this species has been experiencing a steep population decline over the last 5 million years, highlighting the need for international conservation interventions and fisheries-based management for this critical species.

Centromere drive may propel the evolution of chromosome and genome size in plants

BACKGROUND

Genome size is influenced by natural selection and genetic drift acting on variations from polyploidy and repetitive DNA sequences. We hypothesized that centromere drive, where centromeres compete for inclusion in the functional gamete during meiosis, may also affect genome and chromosome size. This competition occurs in asymmetric meiosis, where only one of the four meiotic products becomes a gamete. If centromere drive influences chromosome size evolution, it may also impact post-polyploid diploidization, where a polyploid genome is restructured to function more like a diploid through chromosomal rearrangements, including fusions. We tested if plant lineages with asymmetric meiosis exhibit faster chromosome size evolution compared to those with only symmetric meiosis, which lack centromere drive as all four meiotic products become gametes. We also examined if positive selection on centromeric histone H3 (CENH3), a protein that can suppress centromere drive, is more frequent in these asymmetric lineages.

METHODS

We analysed plant groups with different meiotic modes: asymmetric in gymnosperms and angiosperms, and symmetric in bryophytes, lycophytes and ferns. We selected species based on available CENH3 gene sequences and chromosome size data. Using Ornstein-Uhlenbeck evolutionary models and phylogenetic regressions, we assessed the rates of chromosome size evolution and the frequency of positive selection on CENH3 in these clades.

RESULTS

Our analyses showed that clades with asymmetric meiosis have a higher frequency of positive selection on CENH3 and increased rates of chromosome size evolution compared to symmetric clades.

CONCLUSIONS

Our findings support the hypothesis that centromere drive accelerates chromosome and genome size evolution, potentially also influencing the process of post-polyploid diploidization. We propose a model which in a single framework helps explain the stability of chromosome size in symmetric lineages (bryophytes, lycophytes and ferns) and its variability in asymmetric lineages (gymnosperms and angiosperms), providing a foundation for future research in plant genome evolution.

De Novo Genome Assembly and Annotation for the Synanthropic Webbing Clothes Moth (Tineola bisselliella): A Globally Distributed, Economically Important Pest

Tineola bisselliella, the webbing clothes moth, is an economically important, globally distributed synanthropic pest species and member of the basal moth lineage Tineidae. These moths are facultatively keratinophagous, and their larvae can cause extensive damage, particularly to clothing, textiles, and museum specimens. Despite the economic and phylogenetic importance of T. bisselliella, there is a lack of quality genomic resources for this, or for other species within the Tineidae family. The T. bisselliella genome assembly presented here consists of 30 pseudochromosomes (29 autosomes and 1 Z chromosome) produced using synteny alignment of a preliminary contig-level assembly (256 contigs) to a closely related species, Tinea pellionella. The resulting final pseudochromosome-level assembly is 243.630 Mb and has an N50 length of 8.708 Mb. The assembly is highly contiguous and has similar or improved quality compared to other available Tineidae genomes, with 93.1% (91.8% single copy and 1.3% duplicated) of lepidopteran orthologs complete and present. Annotation of the pseudochromosome-level genome assembly with the transcriptome we produced ultimately yielded 11,259 annotated genes. Synteny alignments between the T. bisselliella genome assembly and other Tineidae genomes revealed evidence for numerous small rearrangements with high synteny conservation. In contrast, a synteny alignment performed between T. bisselliella and Melitaea cinxia, which is thought to have retained the ancestral karyotype (n = 31), revealed a fusion of the ancestral autosome 30 and Z chromosome that led to a reduction in T. bisselliella karyotype size. The reference quality annotated genome for T. bisselliella presented here will advance our understanding of the evolution of the lepidopteran karyotype by providing a chromosome-level genome for this basal moth lineage and provide future insights into the mechanisms underlying keratin digestion in T. bisselliella.

Holocentric repeat landscapes: From micro-evolutionary patterns to macro-evolutionary associations with karyotype evolution

Repetitive elements can cause large-scale chromosomal rearrangements, for example through ectopic recombination, potentially promoting reproductive isolation and speciation. Species with holocentric chromosomes, that lack a localized centromere, might be more likely to retain chromosomal rearrangements that lead to karyotype changes such as fusions and fissions. This is because chromosome segregation during cell division should be less affected than in organisms with a localized centromere. The relationships between repetitive elements and chromosomal rearrangements and how they may translate to patterns of speciation in holocentric organisms are though poorly understood. Here, we use a reference-free approach based on low-coverage short-read sequencing data to characterize the repeat landscape of two independently evolved holocentric groups: Erebia butterflies and Carex sedges. We consider both micro- and macro-evolutionary scales to investigate the repeat landscape differentiation between Erebia populations and the association between repeats and karyotype changes in a phylogenetic framework for both Erebia and Carex. At a micro-evolutionary scale, we found population differentiation in repeat landscape that increases with overall intraspecific genetic differentiation among four Erebia species. At a macro-evolutionary scale, we found indications for an association between repetitive elements and karyotype changes along both Erebia and Carex phylogenies. Altogether, our results suggest that repetitive elements are associated with the level of population differentiation and chromosomal rearrangements in holocentric clades and therefore likely play a role in adaptation and potentially species diversification.

The genome of Eleocharis vivipara elucidates the genetics of C3-C4 photosynthetic plasticity and karyotype evolution in the Cyperaceae

Eleocharis vivipara, an amphibious sedge in the Cyperaceae family, has several remarkable properties, most notably its alternate use of C3 photosynthesis underwater and C4 photosynthesis on land. However, the absence of genomic data has hindered its utility for evolutionary and genetic research. Here, we present a high-quality genome for E. vivipara, representing the first chromosome-level genome for the Eleocharis genus, with an approximate size of 965.22 Mb mainly distributed across 10 chromosomes. Its Hi-C pattern, chromosome clustering results, and one-to-one genome synteny across two subgroups indicates a tetraploid structure with chromosome count 2n = 4x = 20. Phylogenetic analysis suggests that E. vivipara diverged from Cyperus esculentus approximately 32.96 million years ago (Mya), and underwent a whole-genome duplication (WGD) about 3.5 Mya. Numerous fusion and fission events were identified between the chromosomes of E. vivipara and its close relatives. We demonstrate that E. vivipara has holocentromeres, a chromosomal feature which can maintain the stability of such chromosomal rearrangements. Experimental transplantation and cross-section studies showed its terrestrial culms developed C4 Kranz anatomy with increased number of chloroplasts in the bundle sheath (BS) cells. Gene expression and weighted gene co-expression network analysis (WGCNA) showed overall elevated expression of core genes associated with the C4 pathway, and significant enrichment of genes related to modified culm anatomy and photosynthesis efficiency. We found evidence of mixed nicotinamide adenine dinucleotide - malic enzyme and phosphoenolpyruvate carboxykinase type C4 photosynthesis in E. vivipara, and hypothesize that the evolution of C4 photosynthesis predates the WGD event. The mixed type is dominated by subgenome A and supplemented by subgenome B. Collectively, our findings not only shed light on the evolution of E. vivipara and karyotype within the Cyperaceae family, but also provide valuable insights into the transition between C3 and C4 photosynthesis, offering promising avenues for crop improvement and breeding.

The genome assembly of Carex breviculmis provides evidence for its phylogenetic localization and environmental adaptation

BACKGROUND AND AIMS

Carex breviculmis is a perennial herb with good resistance and is widely used for forage production and turf management. It is important in ecology, environmental protection and biodiversity conservation, but faces several challenges due to human activities. However, the absence of genome sequences has limited basic research and the improvement of wild plants.

METHODS

We annotated the genome of C. breviculmis and conducted a systematic analysis to explore its resistance to harsh environments. We also conducted a comparative analysis of Achnatherum splendens, which is similarly tolerant to harsh environments.

KEY RESULTS

The assembled the genome comprises 469.01 Mb, revealing 37 372 genes with a BUSCO completeness score of 99.0 %. The genome has 52.03 % repetitive sequences, primarily influenced by recent LTR insertions that have contributed to its expansion. Phylogenetic analysis suggested that C. breviculmis diverged from C. littledalei ~6.61 million years ago. Investigation of repetitive sequences and expanded gene families highlighted a rapid expansion of tandem duplicate genes, particularly in areas related to sugar metabolism, synthesis of various amino acids, and phenylpropanoid biosynthesis. Additionally, our analysis identified crucial genes involved in secondary metabolic pathways, such as glycolysis, phenylpropanoid biosynthesis and amino acid metabolism, which have undergone positive selection. We reconstructed the sucrose metabolic pathway and identified significant gene expansions, including 16 invertase, 9 sucrose phosphate synthase and 12 sucrose synthase genes associated with sucrose metabolism, which showed varying levels of expansion.

CONCLUSIONS

The expansion of these genes, coupled with subsequent positive selection, contributed to the ability of C. breviculmis to adapt to environmental stressors. This study lays the foundation for future research on the evolution of Carex plants, their environmental adaptations, and potential genetic breeding.

Condensin I folds the Caenorhabditis elegans genome

The structural maintenance of chromosome (SMC) complexes-cohesin and condensins-are crucial for chromosome separation and compaction during cell division. During the interphase, mammalian cohesins additionally fold the genome into loops and domains. Here we show that, in Caenorhabditis elegans, a species with holocentric chromosomes, condensin I is the primary, long-range loop extruder. The loss of condensin I and its X-specific variant, condensin IDC, leads to genome-wide decompaction, chromosome mixing and disappearance of X-specific topologically associating domains, while reinforcing fine-scale epigenomic compartments. In addition, condensin I/IDC inactivation led to the upregulation of X-linked genes and unveiled nuclear bodies grouping together binding sites for the X-targeting loading complex of condensin IDC. C. elegans condensin I/IDC thus uniquely organizes holocentric interphase chromosomes, akin to cohesin in mammals, as well as regulates X-chromosome gene expression.

Diversity of genome size and chromosome number in homothallic and heterothallic strains of the Closterium peracerosum-strigosum-littorale complex (Desmidiales, Zygnematophyceae, Streptophyta)

The evolutionary transitions of mating systems between outcrossing and self-fertilization are often suggested to associate with the cytological and genomic changes, but the empirical reports are limited in multicellular organisms. Here we used the unicellular zygnematophycean algae, the Closterium peracerosum-strigosum-littorale (C. psl.) complex, to address whether genomic properties such as genome sizes and chromosome numbers are associated with mating system transitions between homothallism (self-fertility) and heterothallism (self-sterility). Phylogenetic analysis revealed the polyphyly of homothallic strains, suggesting multiple transitions between homothallism and heterothallism in the C. psl. complex. Flow cytometry analysis identified a more than 2-fold genome size variation, ranging from 0.53 to 1.42 Gbp, which was positively correlated with chromosome number variation between strains. Although we did not find consistent trends in genome size change and mating system transitions, the mean chromosome sizes tend to be smaller in homothallic strains than in their relative heterothallic strains. This result suggests that homothallic strains possibly have more fragmented chromosomes, which is consistent with the argument that self-fertilizing populations may tolerate more chromosomal rearrangements.

Chromosomal evolution in Cryptangieae Benth. (Cyperaceae): Evidence of holocentrism and pseudomonads

Cryptangieae has recently been revised based on morphology and molecular phylogeny, but cytogenetic data is still scarce. We conducted this study with the aim of investigating the occurrence of holocentric chromosomes and pseudomonads, as well as understanding the mode of chromosomal evolution in the tribe. We performed analyses of meiotic behavior, chromosome counts, and reconstruction of the ancestral state for the haploid number. We present novel cytogenetic data for eight potentially holocentric species: Cryptangium verticillatum, Krenakia junciforme, K. minarum, Lagenocarpus bracteosus, L. griseus, L. inversus, L. rigidus, and L. tenuifolius. Meiotic abnormalities were observed, with parallel spindles being particularly noteworthy. Intra-specific variations in chromosome number were not found, which may indicate an efficient genetic control for the elimination of abnormal nuclei. The inferred ancestral haploid number was n = 16, with dysploidy being the main evolutionary mechanism. At least five chromosomal fissions occurred in Krenakia (n = 21), followed by a further ascending dysploidy event in Lagenocarpus (n = 17). As proposed for Cyperaceae, it is possible that cladogenesis events in Cryptangieae were marked by numerical and structural chromosomal changes.

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.

Evidence of Nonrandom Patterns of Functional Chromosome Organization in Danaus plexippus

Our understanding on the interplay between gene functionality and gene arrangement at different chromosome scales relies on a few Diptera and the honeybee, species with quality reference genome assemblies, accurate gene annotations, and abundant transcriptome data. Using recently generated 'omic resources in the monarch butterfly Danaus plexippus, a species with many more and smaller chromosomes relative to Drosophila species and the honeybee, we examined the organization of genes preferentially expressed at broadly defined developmental stages (larva, pupa, adult males, and adult females) at both fine and whole-chromosome scales. We found that developmental stage-regulated genes do not form more clusters, but do form larger clusters, than expected by chance, a pattern consistent across the gene categories examined. Notably, out of the 30 chromosomes in the monarch genome, 12 of them, plus the fraction of the chromosome Z that corresponds to the ancestral Z in other Lepidoptera, were found enriched for developmental stage-regulated genes. These two levels of nonrandom gene organization are not independent as enriched chromosomes for developmental stage-regulated genes tend to harbor disproportionately large clusters of these genes. Further, although paralogous genes were overrepresented in gene clusters, their presence is not enough to explain two-thirds of the documented cases of whole-chromosome enrichment. The composition of the largest clusters often included paralogs from more than one multigene family as well as unrelated single-copy genes. Our results reveal intriguing patterns at the whole-chromosome scale in D. plexippus while shedding light on the interplay between gene expression and chromosome organization beyond Diptera and Hymenoptera.

Phylloxera and Aphids Show Distinct Features of Genome Evolution Despite Similar Reproductive Modes

Genomes of aphids (family Aphididae) show several unusual evolutionary patterns. In particular, within the XO sex determination system of aphids, the X chromosome exhibits a lower rate of interchromosomal rearrangements, fewer highly expressed genes, and faster evolution at nonsynonymous sites compared with the autosomes. In contrast, other hemipteran lineages have similar rates of interchromosomal rearrangement for autosomes and X chromosomes. One possible explanation for these differences is the aphid's life cycle of cyclical parthenogenesis, where multiple asexual generations alternate with 1 sexual generation. If true, we should see similar features in the genomes of Phylloxeridae, an outgroup of aphids which also undergoes cyclical parthenogenesis. To investigate this, we generated a chromosome-level assembly for the grape phylloxera, an agriculturally important species of Phylloxeridae, and identified its single X chromosome. We then performed synteny analysis using the phylloxerid genome and 30 high-quality genomes of aphids and other hemipteran species. Unexpectedly, we found that the phylloxera does not share aphids' patterns of chromosome evolution. By estimating interchromosomal rearrangement rates on an absolute time scale, we found that rates are elevated for aphid autosomes compared with their X chromosomes, but this pattern does not extend to the phylloxera branch. Potentially, the conservation of X chromosome gene content is due to selection on XO males that appear in the sexual generation. We also examined gene duplication patterns across Hemiptera and uncovered horizontal gene transfer events contributing to phylloxera evolution.

Unzipped chromosome-level genomes reveal allopolyploid nematode origin pattern as unreduced gamete hybridization

The formation and consequences of polyploidization in animals with clonal reproduction remain largely unknown. Clade I root-knot nematodes (RKNs), characterized by parthenogenesis and allopolyploidy, show a widespread geographical distribution and extensive agricultural destruction. Here, we generated 4 unzipped polyploid RKN genomes and identified a putative novel alternative telomeric element. Then we reconstructed 4 chromosome-level assemblies and resolved their genome structures as AAB for triploid and AABB for tetraploid. The phylogeny of subgenomes revealed polyploid RKN origin patterns as hybridization between haploid and unreduced gametes. We also observed extensive chromosomal fusions and homologous gene expression decrease after polyploidization, which might offset the disadvantages of clonal reproduction and increase fitness in polyploid RKNs. Our results reveal a rare pathway of polyploidization in parthenogenic polyploid animals and provide a large number of high-precision genetic resources that could be used for RKN prevention and control.

The Impact of Chromosomal Rearrangements in Speciation: From Micro- to Macroevolution

Chromosomal rearrangements (CRs) have been known since almost the beginning of genetics. While an important role for CRs in speciation has been suggested, evidence primarily stems from theoretical and empirical studies focusing on the microevolutionary level (i.e., on taxon pairs where speciation is often incomplete). Although the role of CRs in eukaryotic speciation at a macroevolutionary level has been supported by associations between species diversity and rates of evolution of CRs across phylogenies, these findings are limited to a restricted range of CRs and taxa. Now that more broadly applicable and precise CR detection approaches have become available, we address the challenges in filling some of the conceptual and empirical gaps between micro- and macroevolutionary studies on the role of CRs in speciation. We synthesize what is known about the macroevolutionary impact of CRs and suggest new research avenues to overcome the pitfalls of previous studies to gain a more comprehensive understanding of the evolutionary significance of CRs in speciation across the tree of life.

Phylloxera and aphids show distinct features of genome evolution despite similar reproductive modes

Genomes of aphids (family Aphididae) show several unusual evolutionary patterns. In particular, within the XO sex determination system of aphids, the X chromosome exhibits a lower rate of interchromosomal rearrangements, fewer highly expressed genes, and faster evolution at nonsynonymous sites compared to the autosomes. In contrast, other hemipteran lineages have similar rates of interchromosomal rearrangement for autosomes and X chromosomes. One possible explanation for these differences is the aphid's life cycle of cyclical parthenogenesis, where multiple asexual generations alternate with one sexual generation. If true, we should see similar features in the genomes of Phylloxeridae, an outgroup of aphids which also undergoes cyclical parthenogenesis. To investigate this, we generated a chromosome-level assembly for the grape phylloxera, an agriculturally important species of Phylloxeridae, and identified its single X chromosome. We then performed synteny analysis using the phylloxerid genome and 30 high-quality genomes of aphids and other hemipteran species. Unexpectedly, we found that the phylloxera does not share aphids' patterns of chromosome evolution. By estimating interchromosomal rearrangement rates on an absolute time scale, we found that rates are elevated for aphid autosomes compared to their X chromosomes, but this pattern does not extend to the phylloxera branch. Potentially, the conservation of X chromosome gene content is due to selection on XO males that appear in the sexual generation. We also examined gene duplication patterns across Hemiptera and uncovered horizontal gene transfer events contributing to phylloxera evolution.

Cosegregation of recombinant chromatids maintains genome-wide heterozygosity in an asexual nematode

In asexual animals, female meiosis is modified to produce diploid oocytes. If meiosis still involves recombination, this is expected to lead to a rapid loss of heterozygosity, with adverse effects on fitness. Many asexuals, however, have a heterozygous genome, the underlying mechanisms being most often unknown. Cytological and population genomic analyses in the nematode Mesorhabditis belari revealed another case of recombining asexual being highly heterozygous genome-wide. We demonstrated that heterozygosity is maintained despite recombination because the recombinant chromatids of each chromosome pair cosegregate during the unique meiotic division. A theoretical model confirmed that this segregation bias is necessary to account for the observed pattern and likely to evolve under a wide range of conditions. Our study uncovers an unexpected type of non-Mendelian genetic inheritance involving cosegregation of recombinant chromatids.

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.

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.

Evolution of satDNAs on holocentric chromosomes: insights from hemipteran insects of the genus Mahanarva

Satellite DNAs (satDNAs) constitute one of the main components of eukaryote genomes and are involved in chromosomal organization and diversification. Although largely studied, little information was gathered about their evolution on holocentric species, i.e., diffuse centromeres, which, due to differences in repeat organization, could result in different evolutionary patterns. Here, we combined bioinformatics and cytogenetic approaches to evaluate the evolution of the satellitomes in Mahanarva holocentric insects. In two species, de novo identification revealed a high number of satDNAs, 110 and 113, with an extreme monomer length range of 18-4228 bp. The overall abundance of satDNAs was observed to be 6.67% in M. quadripunctata and 1.98% in M. spectabilis, with different abundances for the shared satDNAs. Chromosomal mapping of the most abundant repeats of M. quadripunctata and M. spectabilis on other Mahanarva reinforced the dynamic nature of satDNAs. Variable patterns of chromosomal distribution for the satDNAs were noticed, with the occurrence of clusters on distinct numbers of chromosomes and at different positions and the occurrence of scattered signals or nonclustered satDNAs. Altogether, our data demonstrated the high dynamism of satDNAs in Mahanarva with the involvement of this genomic fraction in chromosome diversification of the genus. The general characteristics and patterns of evolution of satDNAs are similar to those observed on monocentric chromosomes, suggesting that the differential organization of genome compartments observed on holocentric chromosomes compared with monocentric chromosomes does not have a large impact on the evolution of satDNAs. Analysis of the satellitomes of other holocentric species in a comparative manner will shed light on this issue.

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.

Kinetic Activity of Chromosomes and Expression of Recombination Genes in Achiasmatic Meiosis of Tityus (Archaeotityus) Scorpions

Several species of Tityus (Scorpiones, Buthidae) present multi-chromosomal meiotic associations and failures in the synaptic process, originated from reciprocal translocations. Holocentric chromosomes and achiasmatic meiosis in males are present in all members of this genus. In the present study, we investigated synapse dynamics, transcriptional silencing by γH2AX, and meiotic microtubule association in bivalents and a quadrivalent of the scorpion Tityus maranhensis. Additionally, we performed RT-PCR to verify the expression of mismatch repair enzymes involved in crossing-over formation in Tityus silvestris gonads. The quadrivalent association in T. maranhensis showed delay in the synaptic process and long asynaptic regions during pachytene. In this species, γH2AX was recorded only at the chromosome ends during early stages of prophase I; in metaphase I, bivalents and quadrivalents of T. maranhensis exhibited binding to microtubules along their entire length, while in metaphase II/anaphase II transition, spindle fibers interacted only with telomeric regions. Regarding T. silvestris, genes involved in the recombination process were transcribed in ovaries, testes and embryos, without significant difference between these tissues. The expression of these genes during T. silvestris achiasmatic meiosis is discussed in the present study. The absence of meiotic inactivation by γH2AX and holo/telokinetic behavior of the chromosomes are important factors for the maintenance of the quadrivalent in T. maranhensis and the normal continuation of the meiotic cycle in this species.

Evolution of eukaryotic centromeres by drive and suppression of selfish genetic elements

Despite the universal requirement for faithful chromosome segregation, eukaryotic centromeres are rapidly evolving. It is hypothesized that rapid centromere evolution represents an evolutionary arms race between selfish genetic elements that drive, or propagate at the expense of organismal fitness, and mechanisms that suppress fitness costs. Selfish centromere DNA achieves preferential inheritance in female meiosis by recruiting more effector proteins that alter spindle microtubule interaction dynamics. Parallel pathways for effector recruitment are adaptively evolved to suppress functional differences between centromeres. Opportunities to drive are not limited to female meiosis, and selfish transposons, plasmids and B chromosomes also benefit by maximizing their inheritance. Rapid evolution of selfish genetic elements can diversify suppressor mechanisms in different species that may cause hybrid incompatibility.

Celebrating Mendel, McClintock, and Darlington: On end-to-end chromosome fusions and nested chromosome fusions

The evolution of eukaryotic genomes is accompanied by fluctuations in chromosome number, reflecting cycles of chromosome number increase (polyploidy and centric fissions) and decrease (chromosome fusions). Although all chromosome fusions result from DNA recombination between two or more nonhomologous chromosomes, several mechanisms of descending dysploidy are exploited by eukaryotes to reduce their chromosome number. Genome sequencing and comparative genomics have accelerated the identification of inter-genome chromosome collinearity and gross chromosomal rearrangements and have shown that end-to-end chromosome fusions (EEFs) and nested chromosome fusions (NCFs) may have played a more important role in the evolution of eukaryotic karyotypes than previously thought. The present review aims to summarize the limited knowledge on the origin, frequency, and evolutionary implications of EEF and NCF events in eukaryotes and especially in land plants. The interactions between nonhomologous chromosomes in interphase nuclei and chromosome (mis)pairing during meiosis are examined for their potential importance in the origin of EEFs and NCFs. The remaining open questions that need to be addressed are discussed.

Evolution of sexual systems, sex chromosomes and sex-linked gene transcription in flatworms and roundworms

Many species with separate male and female individuals (termed ‘gonochorism’ in animals) have sex-linked genome regions. Here, we investigate evolutionary changes when genome regions become completely sex-linked, by analyses of multiple species of flatworms (Platyhelminthes; among which schistosomes recently evolved gonochorism from ancestral hermaphroditism), and roundworms (Nematoda) which have undergone independent translocations of different autosomes. Although neither the evolution of gonochorism nor translocations fusing ancestrally autosomal regions to sex chromosomes causes inevitable loss of recombination, we document that formerly recombining regions show genomic signatures of recombination suppression in both taxa, and become strongly genetically degenerated, with a loss of most genes. Comparisons with hermaphroditic flatworm transcriptomes show masculinisation and some defeminisation in schistosome gonad gene expression. We also find evidence that evolution of sex-linkage in nematodes is accompanied by transcriptional changes and dosage compensation. Our analyses also identify sex-linked genes that could assist future research aimed at controlling some of these important parasites.

Transitions between hermaphroditic and separate sexes are relatively understudied in animals compared to pants. Here, Wang et al. reconstruct the evolution of separate sexes in the flatworms and complex changes of sex chromosomes in the roundworms.

Wolffia, a minimalist plant and synthetic biology chassis

A highly simplified species for genome engineering would facilitate rational design of a synthetic plant. A candidate species is the aquatic, non-grass monocot wolffia (Wolffia australiana) in the Lemnaceae family. Commonly known as watermeal, wolffia is a rootless ball of several thousand cells the size of a pinhead and the fastest growing plant known on Earth. Its extreme morphological reduction is coupled to transposon-mediated streamlining of its transcriptome, which represents a core set of nonredundant protein coding genes. Despite its body plan and transcriptome being highly specialized for continuous growth, wolffia retains cell types relevant to higher plants. Systems level studies with this species could enable the creation of a defined biological chassis for synthetic plant construction.

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.

Hybridization fuelled diversification in Spialia butterflies

The importance of hybridization and introgression is well documented in the evolution of plants but, in insects, their role is not fully understood. Given the fact that insects are the most diverse group of organisms, assessing the impact of reticulation events on their evolution may be key to comprehend the emergence of such remarkable diversity. Here, we used an insect model, the Spialia butterflies, to gather genomic evidence of hybridization as a promoter of novel diversity. By using double‐digest RADseq (ddRADseq), we explored the phylogenetic relationships between Spialia orbifer, S. rosae and S. sertorius, and documented two independent events of interspecific gene flow. Our data support that the Iberian endemism S. rosae probably received genetic material from S. orbifer in both mitochondrial and nuclear DNA, which could have contributed to a shift in the ecological preferences of S. rosae. We also show that admixture between S. sertorius and S. orbifer probably occurred in Italy. As a result, the admixed Sicilian populations of S. orbifer are differentiated from the rest of populations both genetically and morphologically, and display signatures of reproductive character displacement in the male genitalia. Additionally, our analyses indicated that genetic material from S. orbifer is present in S. sertorius along the Italian Peninsula. Our findings add to the view that hybridization is a pervasive phenomenon in nature and in butterflies in particular, with important consequences for evolution due to the emergence of novel phenotypes.

Karyotype asymmetry in Cuscuta L. subgenus Pachystigma reflects its repeat DNA composition

Cuscuta is a cytogenetically diverse genus, with karyotypes varying 18-fold in chromosome number and 127-fold in genome size. Each of its four subgenera also presents particular chromosomal features, such as bimodal karyotypes in Pachystigma. We used low coverage sequencing of the Cuscuta nitida genome (subgenus Pachystigma), as well as chromosome banding and molecular cytogenetics of three subgenus representatives, to understand the origin of bimodal karyotypes. All three species, C. nitida, C. africana (2n = 28) and C. angulata (2n = 30), showed heterochromatic bands mainly in the largest chromosome pairs. Eighteen satellite DNAs were identified in C. nitida genome, two showing similarity to mobile elements. The most abundant were present at the largest pairs, as well as the highly abundant ribosomal DNAs. The most abundant Ty1/Copia and Ty3/Gypsy elements were also highly enriched in the largest pairs, except for the Ty3/Gypsy CRM, which also labelled the pericentromeric regions of the smallest chromosomes. This accumulation of repetitive DNA in the larger pairs indicates that these sequences are largely responsible for the formation of bimodal karyotypes in the subgenus Pachystigma. The repetitive DNA fraction is directly linked to karyotype evolution in Cuscuta.

Multiscale analysis of the randomization limits of the chromosomal gene organization between Lepidoptera and Diptera

How chromosome gene organization and gene content evolve among distantly related and structurally malleable genomes remains unresolved. This is particularly the case when considering different insect orders. We have compared the highly contiguous genome assemblies of the lepidopteran Danaus plexippus and the dipteran Drosophila melanogaster, which shared a common ancestor around 290 Ma. The gene content of 23 out of 30 D. plexippus chromosomes was significantly associated with one or two of the six chromosomal elements of the Drosophila genome, denoting common ancestry. Despite the phylogenetic distance, 9.6% of the 1-to-1 orthologues still reside within the same ancestral genome neighbourhood. Furthermore, the comparison D. plexippus-Bombyx mori indicated that the rates of chromosome repatterning are lower in Lepidoptera than in Diptera, although still within the same order of magnitude. Concordantly, 14 developmental gene clusters showed a higher tendency to retain full or partial clustering in D. plexippus, further supporting that the physical association between the SuperHox and NK clusters existed in the ancestral bilaterian. Our results illuminate the scope and limits of the evolution of the gene organization and content of the ancestral chromosomes to the Lepidoptera and Diptera while helping reconstruct portions of the genome in their most recent common ancestor.

The Evolution of Cytogenetic Traits in Cuscuta (Convolvulaceae), the Genus With the Most Diverse Chromosomes in Angiosperms

Karyotypes are characterized by traits such as chromosome number, which can change through whole-genome duplication and dysploidy. In the parasitic plant genus Cuscuta (Convolvulaceae), chromosome numbers vary more than 18-fold. In addition, species of this group show the highest diversity in terms of genome size among angiosperms, as well as a wide variation in the number and distribution of 5S and 35S ribosomal DNA (rDNA) sites. To understand its karyotypic evolution, ancestral character state reconstructions were performed for chromosome number, genome size, and position of 5S and 35S rDNA sites. Previous cytogenetic data were reviewed and complemented with original chromosome counts, genome size estimates, and rDNA distribution assessed via fluorescence in situ hybridization (FISH), for two, seven, and 10 species, respectively. Starting from an ancestral chromosome number of x = 15, duplications were inferred as the prevalent evolutionary process. However, in holocentric clade (subgenus Cuscuta), dysploidy was identified as the main evolutionary mechanism, typical of holocentric karyotypes. The ancestral genome size of Cuscuta was inferred as approximately 1C = 12 Gbp, with an average genome size of 1C = 2.8 Gbp. This indicates an expansion of the genome size relative to other Convolvulaceae, which may be linked to the parasitic lifestyle of Cuscuta. Finally, the position of rDNA sites varied mostly in species with multiple sites in the same karyotype. This feature may be related to the amplification of rDNA sites in association to other repeats present in the heterochromatin. The data suggest that different mechanisms acted in different subgenera, generating the exceptional diversity of karyotypes in Cuscuta.

Aphids and Ants, Mutualistic Species, Share a Mariner Element with an Unusual Location on Aphid Chromosomes

Aphids (Hemiptera, Aphididae) are small phytophagous insects. The aim of this study was to determine if the mariner elements found in the ant genomes are also present in Aphis fabae and Aphis hederae genomes and the possible existence of horizontal transfer events. Aphids maintain a relationship of mutualism with the ants. The close contact between these insects could favour horizontal transfer events of transposable elements. Myrmar mariner element isolated from Myrmica ruginodis and Tapinoma ibericum ants have also been found in the two Aphis species: A. fabae and A. hederae (Afabmar-Mr and Ahedmar-Mr elements). Besides, Afabmar-Mr could be an active transposon. Myrmar-like elements are also present in other insect species as well as in one Crustacean species. The phylogenetic study carried out with all Myrmar-like elements suggests the existence of horizontal transfer. Most aphids have 2n = 8 with a XX-X0 sex determination system. Their complicated life cycle is mostly parthenogenetic with sexual individuals only in autumn. The production of X0 males, originated by XX females which produce only spermatozoa with one X chromosome, must necessarily occur through specialized cytogenetic and molecular mechanisms which are not entirely known. In both aphid species, the mariner elements are located on all chromosomes, including the X chromosomes. However, on the two X chromosomes, no positive signals are detected in their small DAPI-negative telomere regions. The rDNA sites are located, as in the majority of Aphids species, on one of the telomere regions of each X chromosome. The hybridization patterns obtained by double FISH demonstrate that Afabmar-Mr and Ahedmar-Mr elements do not hybridize at the rDNA sites of their host species. Possible causes for the absence of these transposons in the rDNA genes are discussed, probably related with the X chromosome biology.

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

Chromosome Fusion Affects Genetic Diversity and Evolutionary Turnover of Functional Loci but Consistently Depends on Chromosome Size

Major changes in chromosome number and structure are linked to a series of evolutionary phenomena, including intrinsic barriers to gene flow or suppression of recombination due to chromosomal rearrangements. However, chromosome rearrangements can also affect the fundamental dynamics of molecular evolution within populations by changing relationships between linked loci and altering rates of recombination. Here, we build chromosome-level assembly Eueides isabella and, together with a recent chromosome-level assembly of Dryas iulia, examine the evolutionary consequences of multiple chromosome fusions in Heliconius butterflies. These assemblies pinpoint fusion points on 10 of the 20 autosomal chromosomes and reveal striking differences in the characteristics of fused and unfused chromosomes. The ten smallest autosomes in D. iulia and E. isabella, which have each fused to a longer chromosome in Heliconius, have higher repeat and GC content, and longer introns than predicted by their chromosome length. When fused, these characteristics change to become more in line with chromosome length. The fusions also led to reduced diversity, which likely reflects increased background selection and selection against introgression between diverging populations, following a reduction in per-base recombination rate. We further show that chromosome size and fusion impact turnover rates of functional loci at a macroevolutionary scale. Together these results provide further evidence that chromosome fusion in Heliconius likely had dramatic effects on population level processes shaping rates of neutral and adaptive divergence. These effects may have impacted patterns of diversification in Heliconius, a classic example of an adaptive radiation.

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.

Cell Division: Chromatin Dynamics Shape Insect Holocentromeres

Centromeres, the chromosomal loci that ensure chromosome segregation by directing kinetochore assembly, are typically marked by the histone CENP-A. A study in CENP-A-deficient insects finds that virtually any chromosomal region with low nucleosome turnover can assemble kinetochores, highlighting the extraordinary plasticity of holocentromeres.