The Evolution of Chromosome Numbers: Mechanistic Models and Experimental Approaches

Chromosome fusions shaped karyotype evolution and evolutionary relationships in the model family Brassicaceae

The ancestral crucifer karyotype and 22 conserved genomic blocks (CGBs) facilitate phylogenomic analyses in the Brassicaceae. Chromosomal rearrangements reshuffled CGBs of ancestral chromosomes during karyotype evolution. Here, we identify eight protochromosomes representing the common ancestral karyotype (ACBK) of the two Brassicoideae supertribes: Camelinodae (Lineage I) and Brassicodae (Lineage II). The characterization of multiple cascading fusion events allows us to infer evolutionary relationships based on these events. In the Camelinodae, the ACBK first evolved into the AKI genome, which remained conserved in the Cardamineae, whereas it was altered to tAKI by a reciprocal translocation that preceded the diversification of most Camelinodae tribes. The identified fusion breakpoints largely overlap with CGB boundaries, suggesting that CGBs are mainly disrupted by chromosome fusions. Our results demonstrate the stable inheritance of chromosome fusions and their importance for reconstructing evolutionary relationships. The chromosomal breakpoint approach provides a basis for ancestral state reconstruction based on chromosome-level genome assemblies.

Karyotype evolution and speciation in Orthoptera

Karyotype evolution might fuel speciation and can thereby contribute to species diversity. To test the hypothesis that speciation and karyotype change are linked, we estimated anagenetic and cladogenetic rates of karyotype evolution as well as speciation rates in Orthoptera. We compiled the male diploid chromosome number and the number of visible chromosome arms (the fundamental number) from published sources for 1,541 species. Chromosome-associated speciation rates were estimated by jointly modelling cladogenetic and anagenetic character evolution and the phylogenetic birth-death process in a Bayesian statistical framework using a subset of 516 species from 14 families. Our findings unveiled heterogeneity among orthopteran families in the pace of karyotype evolution and whether it was linked to speciation. In 6/14 clades, we found evidence supporting speciation-associated (cladogenetic) karyotype changes, while in 6/14 clades karyotype evolution was primarily anagenetic. The remaining clades (2/14) showed uncertainty in favour of either model. We further analyzed whether flightless phenotype, and thus less mobile species, showed higher rates of karyotype evolution. We showed that the flightless phenotype is associated with the rate of chromosome loss. The finding indicates contrasting patterns of karyotype evolution within specific orthopteran lineages, thus emphasizing substantial diversity in the pace of this evolutionary process. It also implies that substantial changes in chromosome number, arising from instances of chromosomal gains and losses, are recurring events in orthopterans that are associated with reproductive isolation and speciation, at least in some groups.

Chromosome-scale genome assembly of Phyllanthus emblica L. 'Yingyu'

Phyllanthus emblica L. is an edible plant with medicinal properties native to the dry-hot valley of Yunnan, China. Here, we report a de novo chromosome-scale genome of P. emblica wild type 'Yingyu'. 'Yingyu' is an octopoid plant with a total of 104 chromosomes. In total, we assembled and clustered 480 Mb of the genome and constructed 26 pseudochromosomes (haplotypes) of P. emblica wild type 'Yingyu' that encompass 97.9% of the genome and demonstrate to have relatively high integrity. We annotated 31,111 genes found in the genome of P. emblica. We screened 5 different tissues for searching the tissue-specific expression candidate genes. Four unknown function candidate genes were expressed at high levels in the flowers while genes relating to the biosynthesis of gibberellins and cellulose were specifically expressed in the fruits. The ascorbate biosynthesis-related genes were screened on P. emblica 'Yingyu' genome. The high expression level of 2 GDP-mannose epimerases and one L-galactono-1,4- lactone dehydrogenases in the fruit may be related to the activity of absorbate biosynthesis in the fruit. The chromosome-level genomic data for P. emblica we report will be important for the development of molecular markers to facilitate the selection of superior cultivars for processing and pharmaceuticals.

On the way to diploidization and unexpected ploidy in the grass Sporobolus section Spartina mesopolyploids

Plant history is characterized by cyclical whole genome duplication and diploidization with important biological and ecological consequences. Here, we explore the genome history of two related iconic polyploid grasses (Sporobolus alterniflorus and S. maritimus), involved in a well-known example of neopolyploid speciation. We report particular genome dynamics where an ancestral Sporobolus genome (n = 2x = 20) duplicated 9.6-24.4 million years ago (MYA), which was followed by descending dysploidy resulting in a genome with an unexpected base chromosome number (n = 15). This diploidized genome duplicated again 2.1-6.2 MYA to form a tetraploid lineage (2n = 4x = 60), thus reshuffling the ploidy of these species previously thought hexaploids. We also elucidate the mechanism accompanying the speciation between S. maritimus (2n = 60) and S. alterniflorus (2n = 62), resulting from chromosome restructuring, and identify key adaptive genes in the corresponding regions. This represents critical findings to decipher molecular mechanisms underlying species expansion, adaptation to environmental challenge and invasiveness.

Molecular markers and cytogenetics of Eleven O'Clock Portulaca umbraticola: a non-conventional edible ornamental crop

Portulaca umbraticola, commonly known as "Eleven o'clock", is a popular ornamental plant in Brazil, but its potential as a non-conventional food source remains underexplored. Assessing its genetic and cytogenetic diversity is crucial for breeding and selecting optimal accessions. In this study, we analyzed the genetic diversity of P. umbraticola using RAPD markers and chromosomal traits in 20 accessions. We also compared them with P. oleracea to identify potential interspecific hybrids. Accessions were collected from Brazil and cultivated in 3L pots with commercial substrate. Controlled crosses and self-fertilizations were conducted in a greenhouse, followed by cytogenetic analyses. All amplified bands showed polymorphism among 18 accessions, indicating 100% polymorphism. The species displayed considerable dissimilarity, with all 20 accessions showing 2n = 18 chromosomes, while P. oleracea had 2n = 52 chromosomes. P. umbraticola exhibited allogamous reproduction tendencies and self-incompatibility, precluding self-fertilization. Analysis suggests that crosses between accessions 12 and 16 are recommended due to their substantial genetic distance. In conclusion, P. umbraticola and P. oleracea differ in chromosome numbers, and P. umbraticola shows allogamous reproduction with self-incompatibility.

A phylogenetic perspective of chromosome evolution in Formicidae

Chromosomes, as carriers of genes, are the fundamental units of heredity, with the eukaryotic genome divided into multiple chromosomes. Each species typically has a consistent number of chromosomes within its lineage. Ants, however, display remarkable diversity in chromosome numbers, and previous studies have shown that this variation may correlate with ant diversity. As ants evolved, various karyotypes emerged, primarily through chromosomal fissions, leading to an increase in chromosome number and a decrease in chromosome size. In this study, we investigate chromosome evolution in ants from a phylogenetic perspective using ancestral reconstruction. Our analysis indicates that the most recent common ancestor of ants had an ancestral haploid chromosome number of 11, likely composed of biarmed chromosomes. The bimodal distribution of karyotypes and the trend toward increased chromosome numbers align with previous assumptions. However, both dysploidy and ploidy changes have been indicated as likely mechanisms of chromosome number evolution. Descending dysploidy occurs consistently throughout the phylogeny, while changes in ploidy are believed to occur occasionally within the subfamilies during genus diversification. We propose, based on our results and previous evidence (e.g., genome size in ants), that both fusions and fissions contribute equally to karyotype changes in Formicidae. Additionally, changes in ploidy should not be fully ignored, as they can occur across specific lineages.

Allopolyploidy enhances survival advantages for urban environments in the native plant genus Commelina

BACKGROUND AND AIMS

Urbanization-induced environmental changes affect the geographical distribution of natural plant species. This study focused on how polyploidization, a dynamic genome change, influences the survival and distribution of Commelina communis (Cc) and its subspecies C. communis f. ciliata (Ccfc), which have variable chromosome numbers (e.g. Cc, 2n = 88 for Cc; Ccfc, 2n = 46 for Ccfc). The aim was to investigate polyploidization effects on natural plant distribution in urban environments.

METHODS

The geographical distribution across urban-rural gradients was investigated at a total of 218 sites in Japan. Stomata size and density were measured and compared between Cc and Ccfc. Flow cytometry determined genome size and polyploidy. Chromosome karyotyping was performed using the genomic in situ hybridization (GISH) method.

KEY RESULTS

Urban areas were exclusively dominated by Cc, while Cc and Ccfc coexisted in rural areas. Cc had larger and fewer stomata and a genome size more than twice that of Ccfc. GISH results indicated that Cc possesses Ccfc and another unknown genome, suggesting allopolyploidy.

CONCLUSIONS

Our results show that the ploidy difference affects the geographical distribution, stomata traits and genome size between two distinct taxa in the genus Commelina, C. communis as a neo-tetraploid and C. communis f. ciliata, the diploid. Cc is an allopolyploid and is therefore not only polyploidy but also has an additional genome that provides new sets of genes and alleles, contributing to Cc having enhanced survival potentials in urban environments compared with Ccfc. This is the first investigation to clarify the distribution difference related to urban environments, the difference in stomata traits and genome size, and to study chromosome composition in Commelina species.

New estimates and synthesis of chromosome numbers, ploidy levels and genome size variation in Allium sect. Codonoprasum: advancing our understanding of the unresolved diversification and evolution of this section

BACKGROUND

The genus Allium is known for its high chromosomal variability, but most chromosome counts are based on a few individuals and genome size (GS) reports are limited in certain taxonomic groups. This is evident in the Allium sect. Codonoprasum, a species-rich (> 150 species) and taxonomically complex section with weak morphological differences between taxa, the presence of polyploidy and frequent misidentification of taxa. Consequently, a significant proportion of older karyological reports may be unreliable and GS data are lacking for the majority of species within the section. This study, using chromosome counting and flow cytometry (FCM), provides the first comprehensive and detailed insight into variation in chromosome number, polyploid frequency and distribution, and GS in section members, marking a step towards understanding the unresolved diversification and evolution of this group.

RESULTS

We analysed 1578 individuals from 316 populations of 25 taxa and reported DNA ploidy levels and their GS, with calibration from chromosome counts in 22 taxa. Five taxa had multiple ploidy levels. First estimates of GS were obtained for 16 taxa. A comprehensive review of chromosome number and DNA-ploidy levels in 129 taxa of the section revealed that all taxa have x = 8, except A. rupestre with two polyploid series (x = 8, descending dysploidy x = 7), unique for this section. Diploid taxa dominated (72.1%), while di- & polyploid (12.4%) and exclusively polyploid (15.5%) taxa were less common. Ploidy diversity showed that diploid taxa dominated in the eastern Mediterranean and decreased towards the west and north, whereas only polyploid cytotypes of di- & polyploid taxa or exclusively polyploid taxa dominated in northern and northwestern Europe. A 4.1-fold variation in GS was observed across 33 taxa analysed so far (2C = 22.3-92.1 pg), mainly due to polyploidy, with GS downsizing observed in taxa with multiple ploidy levels. Intra-sectional GS variation suggests evolutionary relationships, and intraspecific GS variation within some taxa may indicate taxonomic heterogeneity and/or historical migration patterns.

CONCLUSIONS

Our study showed advantages of FCM as an effective tool for detecting ploidy levels and determining GS within the section. GS could be an additional character in understanding evolution and phylogenetic relationships within the section.

Physiological and Molecular Mechanisms of Lepidopteran Insects: Genomic Insights and Applications of Genome Editing for Future Research

Lepidopteran insects are a major threat to global agriculture, causing significant crop losses and economic damage. Traditional pest control methods are becoming less effective due to the rapid evolution of insecticide resistance. This study explores the current status and genomic characteristics of 1315 Lepidopteran records, alongside an overview of relevant research, utilizing advanced functional genomics techniques, including RNA-seq and CRISPR/Cas9 gene-editing technologies to uncover the molecular mechanisms underlying insecticide resistance. Our genomic analysis revealed significant variability in genome size, assembly quality, and chromosome number, which may influence species' biology and resistance mechanisms. We identified key resistance-associated genes and pathways, including detoxification and metabolic pathways, which help these insects evade chemical control. By employing CRISPR/Cas9 gene-editing techniques, we directly manipulated resistance-associated genes to confirm their roles in resistance, demonstrating their potential for targeted interventions in pest management. These findings emphasize the value of integrating genomic data into the development of effective and sustainable pest control strategies, reducing reliance on chemical insecticides and promoting environmentally friendly integrated pest management (IPM) approaches. Our study highlights the critical role of functional genomics in IPM and its potential to provide long-term solutions to the growing challenge of Lepidopteran resistance.

Integrated genomic, transcriptomic, and metabolomic analyses of Ilex hylonoma provide insights into the triterpenoid saponin biosynthesis

Ilex is known for its rich content of secondary metabolites, particularly triterpenoid saponins. These compounds hold significant value in natural remedies and herbal medicine. However, the molecular mechanisms responsible for triterpenoid biosynthesis in plants of this genus remain largely unexplored. In this study, we successfully generated the first chromosome-scale genome of Ilex hylonoma. The assembly, comprising 20 anchored chromosomes, has an N50 contig size of 2.13 Mb and a scaffold size of 33.68 Mb. Comparative genome analyses with two other congeners with available chromosome-level genomes suggested that an end-to-end chromosome fusion event likely contributed to the reduction in chromosome number from n = 20 to n = 19 within this genus. By integrating transcriptomic and metabolomic data, we identified the gene expression patterns and metabolite profiles of I. hylonoma across three commonly utilized medicinal tissues. We subsequently pinpointed candidate genes involved in the regulation of triterpenoid saponin biosynthesis, including CYP450 genes, UGT genes, and associated transcription factors. Furthermore, yeast heterologous expression analysis revealed that ihyl08363 catalyzed the conversion of β-amyrin into oleanolic acid, while ihyl04303 catalyzed the C-2α hydroxylation of oleanolic acid to produce maslinic acid. This integrated analysis provides valuable insights into the biosynthesis of important triterpenoid saponins in medicinal Ilex plants.

Chromosome-level assembly of Lindenbergia philippensis and comparative genomic analyses shed light on genome evolution in Lamiales

Lamiales, comprising over 23,755 species across 24 families, stands as a highly diverse and prolific plant group, playing a significant role in the cultivation of horticultural, ornamental, and medicinal plant varieties. Whole-genome duplication (WGD) and its subsequent post-polyploid diploidization (PPD) process represent the most drastic type of karyotype evolution, injecting significant potential for promoting the diversity of this lineage. However, polyploidization histories, as well as genome and subgenome fractionation following WGD events in Lamiales species, are still not well investigated. In this study, we constructed a chromosome-level genome assembly of Lindenbergia philippensis (Orobanchaceae) and conducted comparative genomic analyses with 14 other Lamiales species. L. philippensis is positioned closest to the parasitic lineage within Orobanchaceae and has a conserved karyotype. Through a combination of Ks analysis and syntenic depth analysis, we reconstructed and validated polyploidization histories of Lamiales species. Our results indicated that Primulina huaijiensis underwent three rounds of diploidization events following the γ-WGT event, rather than two rounds as reported. Besides, we reconfirmed that most Lamiales species shared a common diploidization event (L-WGD). Subsequently, we constructed the Lamiales Ancestral Karyotype (LAK), comprising 11 proto-chromosomes, and elucidated its evolutionary trajectory, highlighting the highly flexible reshuffling of the Lamiales paleogenome. We identified biased fractionation of subgenomes following the L-WGD event across eight species, and highlighted the positive impacts of non-WGD genes on gene family expansion. This study provides novel genomic resources and insights into polyploidy and karyotype remodeling of Lamiales species, essential for advancing our understanding of species diversification and genome evolution.

Evolutionary patterns of variations in chromosome counts and genome sizes show positive correlations with taxonomic diversity in tropical gingers

PREMISE

Cytogenetic traits such as an organism's chromosome number and genome size are taxonomically critical as they are instrumental in defining angiosperm diversity. Variations in these traits can be traced to evolutionary processes such as polyploidization, although geographic variations across cytogenetic traits remain underexplored. In the pantropical monocot family Zingiberaceae (~1500 species), cytogenetic traits have been well documented; however, the role of these traits in shaping taxonomic diversity and biogeographic patterns of gingers is not known.

METHODS

A time-calibrated Bayesian phylogenetic tree was constructed for 290 taxa covering three of the four subfamilies in Zingiberaceae. We tested models of chromosome number and genome size evolution within the family and whether lineage age, taxonomic diversity, and distributional range explain the variations in the cytogenetic traits. Tests were carried out at two taxonomic ranks: within Zingiberaceae and within genus Hedychium using correlations, generalized linear models and phylogenetic least square models.

RESULTS

The most frequent changes in chromosome number within Zingiberaceae were noted to be demi-polyploidization and polyploidization (~57% of the time), followed by ascending dysploidy (~27%). The subfamily Zingiberoideae showed descending dysploidy at its base, while Alpinioideae showed polyploidization at its internal nodes. Although chromosome counts and genome sizes did not corroborate with each other, suggesting that they are not equivalent; higher chromosome number variations and higher genome size variations were associated with higher taxonomic diversity and wider biogeographic distribution.

CONCLUSIONS

Within Zingiberaceae, multiple incidences of polyploidization were discovered, and cytogenetic events appear to have reduced the genome sizes and increased taxonomic diversity, distributional ranges and invasiveness.

Unfurling an improved method for visualizing mitotic chromosomes in ferns

Premise

Cytotaxonomy employs chromosome visualization to study organismal relationships and evolution. Despite the critical value of cytogenetic data, cytotypes are lacking for many plant groups. Here, we present an improved approach for visualizing mitotic chromosomes in ferns, a key lineage of land plants, using the dividing cells of unfurling croziers (fiddleheads).

Methods and Results

Our modified mitotic chromosome preparation incorporates a brief pectinase-cellulase pretreatment, as well as colchicine fixation and the Feulgen reaction to improve the staining and separation of mitotic chromosomes. To demonstrate this easy and efficient assessment, we determined the sporophytic (2n) chromosome number for three fern species: Cheilanthes mollis (2n = 60), Cheilanthes hypoleuca (2n = 120), and Nephrolepis cordifolia (2n = 82).

Conclusions

The new method presented here improves visualizations of mitotic chromosomes from the dividing nuclei of young fern croziers. Fiddleheads are widely accessible in nature and in living collections worldwide, and this modified approach increases their suitability for fern cytotaxonomic studies.

A locus-dependent mixed inheritance in the segmental allohexaploid sweetpotato (Ipomoea batatas [L.] Lam)

Two interrelated aspects of the sweetpotato genome, its polyploid origin and inheritance type, remain uncertain. We recently proposed a segmental allohexaploid sweetpotato and thus sought to clarify its inheritance type by direct analyses of homoeolog segregations at selected single-copy loci. For such analyses, we developed a digital quantitative PCR genotyping method using one nondiscriminatory and three discriminatory probes for each selected locus to discriminate and quantify three homoeolog-differentiating variation types (homoeolog-types) in genomic DNA samples for genotype fitting and constructed a F2 population for segregation analyses. We confirmed inter-subgenomic distinctions of three identified homoeolog-types at each of five selected loci by their interspecific differentiations among 14 species in Ipomoea section batatas and genotyped the loci in 549 F2 lines, selected F1 progenies, and their founding parents. Segregation and genotype analyses revealed a locus-dependent mixed inheritance (disomic, polysomic, and intermediate types) of the homoeolog-types at 4 loci in the F2 population, displaying estimated disomic-inheritance frequencies of 0, 2.72%, 14.52%, and 36.92%, and probably in the F1 population too. There were also low-frequency non-hexaploid F1 and F2 genotypes that were probably derived from double-reduction recombination or partially unreduced gametes, and F2 genotypes of apparent aneuploids/dysploids with neopolyploid-like frequencies. Additional analyses of homoeolog-type genotypes at the 5 loci in 46 lines from various regions revealed locus-dependent selection biases, favoring genotypes having more of one homoeolog-type, i.e. more of di- or homogenized homoeolog-type composition, and one-direction ploidy trending among apparent aneuploids/dysploids. These inheritance features pointed to an evolving segmental allohexaploid sweetpotato impacted by selection biases.

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.

Plant chromosome engineering - past, present and future

Spontaneous chromosomal rearrangements (CRs) play an essential role in speciation, genome evolution and crop domestication. To be able to use the potential of CRs for breeding, plant chromosome engineering was initiated by fragmenting chromosomes by X-ray irradiation. With the rise of the CRISPR/Cas system, it became possible to induce double-strand breaks (DSBs) in a highly efficient manner at will at any chromosomal position. This has enabled a completely new level of predesigned chromosome engineering. The genetic linkage between specific genes can be broken by inducing chromosomal translocations. Natural inversions, which suppress genetic exchange, can be reverted for breeding. In addition, various approaches for constructing minichromosomes by downsizing regular standard A or supernumerary B chromosomes, which could serve as future vectors in plant biotechnology, have been developed. Recently, a functional synthetic centromere could be constructed. Also, different ways of genome haploidization have been set up, some based on centromere manipulations. In the future, we expect to see even more complex rearrangements, which can be combined with previously developed engineering technologies such as recombinases. Chromosome engineering might help to redefine genetic linkage groups, change the number of chromosomes, stack beneficial genes on mini cargo chromosomes, or set up genetic isolation to avoid outcrossing.

Integration of fluorescence in situ hybridization and chromosome-length genome assemblies revealed synteny map for guinea pig, naked mole-rat, and human

Descriptions of karyotypes of many animal species are currently available. In addition, there has been a significant increase in the number of sequenced genomes and an ever-improving quality of genome assembly. To close the gap between genomic and cytogenetic data we applied fluorescent in situ hybridization (FISH) and Hi-C technology to make the first full chromosome-level genome comparison of the guinea pig (Cavia porcellus), naked mole-rat (Heterocephalus glaber), and human. Comparative chromosome maps obtained by FISH with chromosome-specific probes link genomic scaffolds to individual chromosomes and orient them relative to centromeres and heterochromatic blocks. Hi-C assembly made it possible to close all gaps on the comparative maps and to reveal additional rearrangements that distinguish the karyotypes of the three species. As a result, we integrated the bioinformatic and cytogenetic data and adjusted the previous comparative maps and genome assemblies of the guinea pig, naked mole-rat, and human. Syntenic associations in the two hystricomorphs indicate features of their putative ancestral karyotype. We postulate that the two approaches applied in this study complement one another and provide complete information about the organization of these genomes at the chromosome level.

PloiDB: the plant ploidy database

See also the Commentary on this article by Spoelhof et al., 240: 909–911.

Dramatic difference in rate of chromosome number evolution among sundew (Drosera L., Droseraceae) lineages

Chromosome number change is a driver of speciation in eukaryotic organisms. Carnivorous sundews in the plant genus Drosera L. exhibit single chromosome number variation both among and within species, especially in the Australian Drosera subg. Ergaleium D.C., potentially linked to atypical centromeres that span much of the length of the chromosomes. We critically reviewed the literature on chromosome counts in Drosera, verified the taxonomy and quality of the original counts, and reconstructed dated phylogenies. We used the BiChrom model to test whether rates of single chromosome number increase and decrease, and chromosome number doubling differed between D. subg. Ergaleium and the other subgenera and between self-compatible and self-incompatible lineages. The best model for chromosome evolution among subgenera had equal rates of chromosome number doubling but higher rates of single chromosome number change in D. subg. Ergaleium than in the other subgenera. Contrary to expectation, self-incompatible lineages had a significantly higher rate of single chromosome loss than self-compatible lineages. We found no evidence for an association between differences in single chromosome number changes and diploidization after polyploidy or centromere type. This study presents an exemplar for critically examining published cytological data and rigorously testing factors that may impact the rates of chromosome number evolution.

Microchromosome fusions underpin convergent evolution of chameleon karyotypes

Evolutionary shifts in chromosome compositions (karyotypes) are major drivers of lineage and genomic diversification. Fusion of ancestral chromosomes is one hypothesized mechanism for the evolutionary reduction of the total chromosome number, a frequently implied karyotypic shift. Empirical tests of this hypothesis require model systems with variable karyotypes, known chromosome features, and a robust phylogeny. Here we used chameleons, diverse lizards with exceptionally variable karyotypes ($2n=20\text{-}62$), to test whether chromosomal fusions explain the repeated evolution of karyotypes with fewer chromosomes than ancestral karyotypes. Using a multidisciplinary approach including cytogenetic analyses and phylogenetic comparative methods, we found that a model of constant loss through time best explained chromosome evolution across the chameleon phylogeny. Next, we tested whether fusions of microchromosomes into macrochromosomes explained these evolutionary losses using generalized linear models. Multiple comparisons supported microchromosome fusions as the predominant agent of evolutionary loss. We further compared our results to various natural history traits and found no correlations. As such, we infer that the tendency of microchromosomes to fuse was a quality of the ancestral chameleon genome and that the genomic predisposition of ancestors is a more substantive predictor of chromosome change than the ecological, physiological, and biogeographical factors involved in their diversification.

Dicentric chromosome breakage in Drosophila melanogaster is influenced by pericentric heterochromatin and occurs in nonconserved hotspots

Chromosome breakage plays an important role in the evolution of karyotypes and can produce deleterious effects within a single individual, such as aneuploidy or cancer. Forces that influence how and where chromosomes break are not fully understood. In humans, breakage tends to occur in conserved hotspots called common fragile sites (CFS), especially during replication stress. By following the fate of dicentric chromosomes in Drosophila melanogaster, we find that breakage under tension also tends to occur in specific hotspots. Our experimental approach was to induce sister chromatid exchange in a ring chromosome to generate a dicentric chromosome with a double chromatid bridge. In the following cell division, the dicentric bridges may break. We analyzed the breakage patterns of 3 different ring-X chromosomes. These chromosomes differ by the amount and quality of heterochromatin they carry as well as their genealogical history. For all 3 chromosomes, breakage occurs preferentially in several hotspots. Surprisingly, we found that the hotspot locations are not conserved between the 3 chromosomes: each displays a unique array of breakage hotspots. The lack of hotspot conservation, along with a lack of response to aphidicolin, suggests that these breakage sites are not entirely analogous to CFS and may reveal new mechanisms of chromosome fragility. Additionally, the frequency of dicentric breakage and the durability of each chromosome's spindle attachment vary significantly between the 3 chromosomes and are correlated with the origin of the centromere and the amount of pericentric heterochromatin. We suggest that different centromere strengths could account for this.

Three recent sex chromosome-to-autosome fusions in a Drosophila virilis strain with high satellite DNA content

The karyotype, or number and arrangement of chromosomes, has varying levels of stability across both evolution and disease. Karyotype changes often originate from DNA breaks near the centromeres of chromosomes, which generally contain long arrays of tandem repeats or satellite DNA. Drosophila virilis possesses among the highest relative satellite abundances of studied species, with almost half its genome composed of three related 7 bp satellites. We discovered a strain of D. virilis that we infer recently underwent three independent chromosome fusion events involving the X and Y chromosomes, in addition to one subsequent fission event. Here, we isolate and characterize the four different karyotypes we discovered in this strain which we believe demonstrates remarkable genome instability. We discovered that one of the substrains with an X-autosome fusion has an X-to-Y chromosome nondisjunction rate 20 × higher than the D. virilis reference strain (21% vs 1%). Finally, we found an overall higher rate of DNA breakage in the substrain with higher satellite DNA compared to a genetically similar substrain with less satellite DNA. This suggests that satellite DNA abundance may play a role in the risk of genome instability. Overall, we introduce a novel system consisting of a single strain with four different karyotypes, which we believe will be useful for future studies of genome instability, centromere function, and sex chromosome evolution.

Chromosome-length genome assemblies and cytogenomic analyses of pangolins reveal remarkable chromosome counts and plasticity

We report the first chromosome-length genome assemblies for three species in the mammalian order Pholidota: the white-bellied, Chinese, and Sunda pangolins. Surprisingly, we observe extraordinary karyotypic plasticity within this order and, in female white-bellied pangolins, the largest number of chromosomes reported in a Laurasiatherian mammal: 2n = 114. We perform the first karyotype analysis of an African pangolin and report a Y-autosome fusion in white-bellied pangolins, resulting in 2n = 113 for males. We employ a novel strategy to confirm the fusion and identify the autosome involved by finding the pseudoautosomal region (PAR) in the female genome assembly and analyzing the 3D contact frequency between PAR sequences and the rest of the genome in male and female white-bellied pangolins. Analyses of genetic variability show that white-bellied pangolins have intermediate levels of genome-wide heterozygosity relative to Chinese and Sunda pangolins, consistent with two moderate declines of historical effective population size. Our results reveal a remarkable feature of pangolin genome biology and highlight the need for further studies of these unique and endangered mammals.

An improved reference genome for Trifolium subterraneum L. provides insight into molecular diversity and intra-specific phylogeny

Subterranean clover (Trifolium subterraneum L., Ts) is a geocarpic, self-fertile annual forage legume with a compact diploid genome (n = x = 8, 544 Mb/1C). Its resilience and climate adaptivity have made it an economically important species in Mediterranean and temperate zones. Using the cultivar Daliak, we generated higher resolution sequence data, created a new genome assembly TSUd_3.0, and conducted molecular diversity analysis for copy number variant (CNV) and single-nucleotide polymorphism (SNP) among 36 cultivars. TSUd_3.0 substantively improves prior genome assemblies with new Hi-C and long-read sequence data, covering 531 Mb, containing 41,979 annotated genes and generating a 94.4% BUSCO score. Comparative genomic analysis among select members of the tribe Trifolieae indicated TSUd 3.0 corrects six assembly-error inversion/duplications and confirmed phylogenetic relationships. Its synteny with T. pratense, T. repens, Medicago truncatula and Lotus japonicus genomes were assessed, with the more distantly related T. repens and M. truncatula showing higher levels of co-linearity with Ts than between Ts and its close relative T. pratense. Resequencing of 36 cultivars discovered 7,789,537 SNPs subsequently used for genomic diversity assessment and sequence-based clustering. Heterozygosity estimates ranged from 1% to 21% within the 36 cultivars and may be influenced by admixture. Phylogenetic analysis supported subspecific genetic structure, although it indicates four or five groups, rather than the three recognized subspecies. Furthermore, there were incidences where cultivars characterized as belonging to a particular subspecies clustered with another subspecies when using genomic data. These outcomes suggest that further investigation of Ts sub-specific classification using molecular and morpho-physiological data is needed to clarify these relationships. This upgraded reference genome, complemented with comprehensive sequence diversity analysis of 36 cultivars, provides a platform for future gene functional analysis of key traits, and genome-based breeding strategies for climate adaptation and agronomic performance. Pangenome analysis, more in-depth intra-specific phylogenomic analysis using the Ts core collection, and functional genetic and genomic studies are needed to further augment knowledge of Trifolium genomes.

Whole-Genome Analysis Reveals the Dynamic Evolution of Holocentric Chromosomes in Satyrine Butterflies

Butterfly chromosomes are holocentric, i.e., lacking a localized centromere. Potentially, this can lead to rapid karyotypic evolution through chromosome fissions and fusions, since fragmented chromosomes retain kinetic activity, while fused chromosomes are not dicentric. However, the actual mechanisms of butterfly genome evolution are poorly understood. Here, we analyzed chromosome-scale genome assemblies to identify structural rearrangements between karyotypes of satyrine butterfly species. For the species pair Erebia ligea-Maniola jurtina, sharing the ancestral diploid karyotype 2n = 56 + ZW, we demonstrate a high level of chromosomal macrosynteny and nine inversions separating these species. We show that the formation of a karyotype with a low number of chromosomes (2n = 36 + ZW) in Erebia aethiops was based on ten fusions, including one autosome-sex chromosome fusion, resulting in a neo-Z chromosome. We also detected inversions on the Z sex chromosome that were differentially fixed between the species. We conclude that chromosomal evolution is dynamic in the satyrines, even in the lineage that preserves the ancestral chromosome number. We hypothesize that the exceptional role of Z chromosomes in speciation may be further enhanced by inversions and sex chromosome-autosome fusions. We argue that not only fusions/fissions but also inversions are drivers of the holocentromere-mediated mode of chromosomal speciation.

High-density linkage maps and chromosome level genome assemblies unveil direction and frequency of extensive structural rearrangements in wood white butterflies (Leptidea spp.)

Karyotypes are generally conserved between closely related species and large chromosome rearrangements typically have negative fitness consequences in heterozygotes, potentially driving speciation. In the order Lepidoptera, most investigated species have the ancestral karyotype and gene synteny is often conserved across deep divergence, although examples of extensive genome reshuffling have recently been demonstrated. The genus Leptidea has an unusual level of chromosome variation and rearranged sex chromosomes, but the extent of restructuring across the rest of the genome is so far unknown. To explore the genomes of the wood white (Leptidea) species complex, we generated eight genome assemblies using a combination of 10X linked reads and HiC data, and improved them using linkage maps for two populations of the common wood white (L. sinapis) with distinct karyotypes. Synteny analysis revealed an extensive amount of rearrangements, both compared to the ancestral karyotype and between the Leptidea species, where only one of the three Z chromosomes was conserved across all comparisons. Most restructuring was explained by fissions and fusions, while translocations appear relatively rare. We further detected several examples of segregating rearrangement polymorphisms supporting a highly dynamic genome evolution in this clade. Fusion breakpoints were enriched for LINEs and LTR elements, which suggests that ectopic recombination might be an important driver in the formation of new chromosomes. Our results show that chromosome count alone may conceal the extent of genome restructuring and we propose that the amount of genome evolution in Lepidoptera might still be underestimated due to lack of taxonomic sampling.

South American Plant Chromosome Numbers Databases: The Information We Have and the Information We Lack on the Most Plant-Diverse Continent

Chromosome numbers have been used in plant taxonomy, and they are still fundamental for taxon delimitation and genome evolution studies. South America is one of the most diverse continents in terms of plant species and there is a considerable number of species not yet analyzed. Accumulated knowledge about plant chromosome numbers has been compiled from online databases, and here we present an overview. The CCDB is one of the largest plant cytological databases and includes data for around 18% of known vascular plants in the world. In this work, we review the information contained in CCDB and in three databases with exclusive information for South America. At present, the three existing databases comprise information on around 1800 plant taxa related to specific regions, countries, or biomes. Efforts are necessary to expand cytological knowledge and to collect all the available information in a plant chromosome database for this continent.

The Chromosome Counts Database (CCDB)

For decades, plant biologists have been interested in the determination and documentation of chromosome numbers for extant taxa. This central cytological character has been used as an important phylogenetic marker and as an indicator for major genomic events such as polyploidy and dysploidy. Due to their significance and the relative ease by which chromosome numbers can be obtained, chromosome numbers have been extensively recorded across the plant kingdom and documented in a wide variety of resources. This makes the collection process a wearing task, often leading to partial data retrieval. In 2015, the Chromosome Counts Database (CCDB) was assembled, being an online unified community resource. This database compiles dozens of different chromosome counts sources, of which a significant portion had been unavailable before in a digitized, searchable format. The vast amount of data assembled in CCDB has already enabled a large number of analyses to examine the evolution of different plant hierarchies, as well as the application of various follow-up analyses, such as ploidy-level inference using chromEvol. CCDB ( http://ccdb.tau.ac.il/ ) encourages data sharing among the botanical community and is expected to continue expanding as additional chromosome numbers are recorded.

A Chromosome-level assembly of the Japanese eel genome, insights into gene duplication and chromosomal reorganization

Japanese eels (Anguilla japonica) are commercially important species, harvested extensively for food. Currently, this and related species (American and European eels) are challenging to breed on a commercial basis. As a result, the wild stock is used for aquaculture. Moreover, climate change, habitat loss, water pollution, and altered ocean currents affect eel populations negatively. Accordingly, the International Union for Conservation of Nature lists Japanese eels as endangered and on its red list. Here we presented a high-quality genome assembly for Japanese eels and demonstrated that large chromosome reorganizations occurred in the events of third-round whole-genome duplications (3R-WRDs). Several chromosomal fusions and fissions have reduced the ancestral protochromosomal number of 25 to 19 in the Anguilla lineage. A phylogenetic analysis of the expanded gene families showed that the olfactory receptors (group δ and ζ genes) and voltage-gated Ca2+ channels expanded significantly. Both gene families are crucial for olfaction and neurophysiology. Additional tandem and proximal duplications occurred following 3R-WGD to acquire immune-related genes for an adaptive advantage against various pathogens. The Japanese eel assembly presented here can be used to study other Anguilla species relating to evolution and conservation.

Comprehensive analysis of 1R- and 2R-MYBs reveals novel genic and protein features, complex organisation, selective expansion and insights into evolutionary tendencies

Myeloblastosis (MYB) family, the largest plant transcription factor family, has been subcategorised based on the number and type of repeats in the MYB domain. In spite of several reports, evolution of MYB genes and repeats remains enigmatic. Brassicaceae members are endowed with complex genomes, including dysploidy because of its unique history with multiple rounds of polyploidisation, genomic fractionations and rearrangements. The present study is an attempt to gain insights into the complexities of MYB family diversity, understand impacts of genome evolution on gene families and develop an evolutionary framework to understand the origin of various subcategories of MYB gene family. We identified and analysed 1129 MYBs that included 1R-, 2R-, 3R- and atypical-MYBs across sixteen species representing protists, fungi, animals and plants and exclude MYB identified from Brassicaceae except Arabidopsis thaliana; in addition, a total of 1137 2R-MYB genes from six Brassicaceae species were also analysed. Comparative analysis revealed predominance of 1R-MYBs in protists, fungi, animals and lower plants. Phylogenetic reconstruction and analysis of selection pressure suggested ancestral nature of R1-type repeat containing 1R-MYBs that might have undergone intragenic duplication to form multi-repeat MYBs. Distinct differences in gene structure between 1R-MYB and 2R-MYBs were observed regarding intron number, the ratio of gene length to coding DNA sequence (CDS) length and the length of exons encoding the MYB domain. Conserved as well as novel and lineage-specific intron phases were identified. Analyses of physicochemical properties revealed drastic differences indicating functional diversification in MYBs. Phylogenetic reconstruction of 1R- and 2R-MYB genes revealed a shared structure-function relationship in clades which was supported when transcriptome data was analysed in silico. Comparative genomics to study distribution pattern and mapping of 2R-MYBs revealed congruency and greater degree of synteny and collinearity among closely related species. Micro-synteny analysis of genomic segments revealed high conservation of genes that are immediately flanking the surrounding tandemly organised 2R-MYBs along with instances of local duplication, reorganisations and genome fractionation. In summary, polyploidy, dysploidy, reshuffling and genome fractionation were found to cause loss or gain of 2R-MYB genes. The findings need to be supported with functional validation to understand gene structure-function relationship along the evolutionary lineage and adaptive strategies based on comparative functional genomics in plants.

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.

Born in the mitochondrion and raised in the nucleus: evolution of a novel tandem repeat family in Medicago polymorpha (Fabaceae)

Plant nuclear genomes harbor sequence elements derived from the organelles (mitochondrion and plastid) through intracellular gene transfer (IGT). Nuclear genomes also show a dramatic range of repeat content, suggesting that any sequence can be readily amplified. These two aspects of plant nuclear genomes are well recognized but have rarely been linked. Through investigation of 31 Medicago taxa we detected exceptionally high post-IGT amplification of mitochondrial (mt) DNA sequences containing rps10 in the nuclear genome of Medicago polymorpha and closely related species. The amplified sequences were characterized as tandem arrays of five distinct repeat motifs (2157, 1064, 987, 971, and 587 bp) that have diverged from the mt genome (mitogenome) in the M. polymorpha nuclear genome. The mt rps10-like arrays were identified in seven loci (six intergenic and one telomeric) of the nuclear chromosome assemblies and were the most abundant tandem repeat family, representing 1.6-3.0% of total genomic DNA, a value approximately three-fold greater than the entire mitogenome in M. polymorpha. Compared to a typical mt gene, the mt rps10-like sequence coverage level was 691.5-7198-fold higher in M. polymorpha and closely related species. In addition to the post-IGT amplification, our analysis identified the canonical telomeric repeat and the species-specific satellite arrays that are likely attributable to an ancestral chromosomal fusion in M. polymorpha. A possible relationship between chromosomal instability and the mt rps10-like tandem repeat family in the M. polymorpha clade is discussed.

Chromosome Number and Genome Size Evolution in Brasolia and Sobralia (Sobralieae, Orchidaceae)

Despite the clear circumscription of tribe Sobralieae (Orchidaceae), its internal relationships are still dubious. The recently delimited genus Brasolia, based on previous Sobralia species, is now assumed to be paraphyletic, with a third genus, Elleanthus, nested in it. The morphology of these three genera is significantly different, indicating the necessity of new data for a better genera delimitation. Though morphology and molecular data are available, cytogenetics data for Sobralieae is restricted to two Sobralia and one Elleanthus species. Aiming to evaluate the potential of cytogenetic data for Brasolia-Elleanthus-Sobralia genera delimitation, we present chromosome number and genome size data for 21 and 20 species, respectively, and used such data to infer the pattern of karyotype evolution in these genera. The analysis allowed us to infer x = 24 as the base chromosome number and genome size of average 1C-value of 5.0 pg for the common ancestor of Brasolia-Elleanthus-Sobralia. The recurrent descending dysploidy in Sobralieae and the punctual genome upsize suggest a recent diversification in Sobralieae but did not allow differing between Brasolia and Sobralia. However, the basal position of tribe Sobralieae in the subfamily Epidendroideae makes this tribe of interest to further studies clarifying the internal delimitation and pattern of karyotype evolution.

Tracing the Evolution of the Angiosperm Genome from the Cytogenetic Point of View

Cytogenetics constitutes a branch of genetics that is focused on the cellular components, especially chromosomes, in relation to heredity and genome structure, function and evolution. The use of modern cytogenetic approaches and the latest microscopes with image acquisition and processing systems enables the simultaneous two- or three-dimensional, multicolour visualisation of both single-copy and highly-repetitive sequences in the plant genome. The data that is gathered using the cytogenetic methods in the phylogenetic background enable tracing the evolution of the plant genome that involve changes in: (i) genome sizes; (ii) chromosome numbers and morphology; (iii) the content of repetitive sequences and (iv) ploidy level. Modern cytogenetic approaches such as FISH using chromosome- and genome-specific probes have been widely used in studies of the evolution of diploids and the consequences of polyploidy. Nowadays, modern cytogenetics complements analyses in other fields of cell biology and constitutes the linkage between genetics, molecular biology and genomics.

Chromosome restructuring and number change during the evolution of Morus notabilis and Morus alba

Mulberry (Morus spp.) is an economically important plant as the main food plant used for rearing domesticated silkworm and it has multiple uses in traditional Chinese medicine. Two basic chromosome numbers (Morus notabilis, n = 7, and Morus alba, n = 14) have been reported in the genus Morus, but the evolutionary history and relationship between them remain unclear. In the present study, a 335-Mb high-quality chromosome-scale genome was assembled for the wild mulberry species M. notabilis. Comparative genomic analyses indicated high chromosomal synteny between the 14 chromosomes of cultivated M. alba and the six chromosomes of wild M. notabilis. These results were successfully verified by fluorescence in situ hybridization. Chromosomal fission/fusion events played crucial roles in the chromosome restructuring process between M. notabilis and M. alba. The activity of the centromere was another key factor that ensured the stable inheritance of chromosomes. Our results also revealed that long terminal repeat retrotransposons were a major driver of the genome divergence and evolution of the mulberry genomes after they diverged from each other. This study provides important insights and a solid foundation for studying the evolution of mulberry, allowing the accelerated genetic improvement of cultivated mulberry species.

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.

Diploidy within a Haploid Genus of Entomopathogenic Fungi

Fungi in the genus Metarhizium are soil-borne plant-root endophytes and rhizosphere colonizers, but also potent insect pathogens with highly variable host ranges. These ascomycete fungi are predominantly asexually reproducing and ancestrally haploid, but two independent origins of persistent diploidy within the Coleoptera-infecting Metarhizium majus species complex are known and has been attributed to incomplete chromosomal segregation following meiosis during the sexual cycle. There is also evidence for infrequent sexual cycles in the locust-specific pathogenic fungus Metarhizium acridum (Hypocreales: Clavicipitaceae), which is an important entomopathogenic biocontrol agent used for the control of grasshoppers in agricultural systems as an alternative to chemical control. Here, we show that the genome of the M. acridum isolate ARSEF 324, which is formulated and commercially utilized is functionally diploid. We used single-molecule real-time sequencing technology to complete a high-quality assembly of ARSEF 324. K-mer frequencies, intragenomic collinearity between contigs and single nucleotide variant read depths across the genome revealed the first incidence of diploidy described within the species M. acridum. The haploid assembly of 44.7 Mb consisted of 20.8% repetitive elements, which is the highest proportion described of any Metarhizium species. The long-read diploid genome assembly sheds light on past research on this strain, such as unusual high UVB tolerance. The data presented here could fuel future investigation into the fitness landscape of fungi with infrequent sexual reproduction and aberrant ploidy levels, not least in the context of biocontrol agents.