Mechanisms of chromosome number reduction in Arabidopsis thaliana and related Brassicaceae species

The structure, function, and evolution of plant centromeres

Centromeres are essential regions of eukaryotic chromosomes responsible for the formation of kinetochore complexes, which connect to spindle microtubules during cell division. Notably, although centromeres maintain a conserved function in chromosome segregation, the underlying DNA sequences are diverse both within and between species and are predominantly repetitive in nature. The repeat content of centromeres includes high-copy tandem repeats (satellites), and/or specific families of transposons. The functional region of the centromere is defined by loading of a specific histone 3 variant (CENH3), which nucleates the kinetochore and shows dynamic regulation. In many plants, the centromeres are composed of satellite repeat arrays that are densely DNA methylated and invaded by centrophilic retrotransposons. In some cases, the retrotransposons become the sites of CENH3 loading. We review the structure of plant centromeres, including monocentric, holocentric, and metapolycentric architectures, which vary in the number and distribution of kinetochore attachment sites along chromosomes. We discuss how variation in CENH3 loading can drive genome elimination during early cell divisions of plant embryogenesis. We review how epigenetic state may influence centromere identity and discuss evolutionary models that seek to explain the paradoxically rapid change of centromere sequences observed across species, including the potential roles of recombination. We outline putative modes of selection that could act within the centromeres, as well as the role of repeats in driving cycles of centromere evolution. Although our primary focus is on plant genomes, we draw comparisons with animal and fungal centromeres to derive a eukaryote-wide perspective of centromere structure and function.

Translocations and inversions: major chromosomal rearrangements during Vigna (Leguminosae) evolution

KEY MESSAGE

Inversions and translocations are the major chromosomal rearrangements involved in Vigna subgenera evolution, being Vigna vexillata the most divergent species. Centromeric repositioning seems to be frequent within the genus. Oligonucleotide-based fluorescence in situ hybridization (Oligo-FISH) provides a powerful chromosome identification system for inferring plant chromosomal evolution. Aiming to understand macrosynteny, chromosomal diversity, and the evolution of bean species from five Vigna subgenera, we constructed cytogenetic maps for eight taxa using oligo-FISH-based chromosome identification. We used oligopainting probes from chromosomes 2 and 3 of Phaseolus vulgaris L. and two barcode probes designed from V. unguiculata (L.) Walp. genome. Additionally, we analyzed genomic blocks among the Ancestral Phaseoleae Karyotype (APK), two V. unguiculata subspecies (V. subg. Vigna), and V. angularis (Willd.) Ohwi & Ohashi (V. subg. Ceratotropis). We observed macrosynteny for chromosomes 2, 3, 4, 6, 7, 8, 9, and 10 in all investigated taxa except for V. vexillata (L.) A. Rich (V. subg. Plectrotropis), in which only chromosomes 4, 7, and 9 were unambiguously identified. Collinearity breaks involved with chromosomes 2 and 3 were revealed. We identified minor differences in the painting pattern among the subgenera, in addition to multiple intra- and interblock inversions and intrachromosomal translocations. Other rearrangements included a pericentric inversion in chromosome 4 (V. subg. Vigna), a reciprocal translocation between chromosomes 1 and 5 (V. subg. Ceratotropis), a potential deletion in chromosome 11 of V. radiata (L.) Wilczek, as well as multiple intrablock inversions and centromere repositioning via genomic blocks. Our study allowed the visualization of karyotypic patterns in each subgenus, revealing important information for understanding intrageneric karyotypic evolution, and suggesting V. vexillata as the most karyotypically divergent species.

Cyto-molecular characterization of rDNA and chromatin composition in the NOR-associated satellite in Chestnut (Castanea spp.)

The American chestnut (Castanea dentata, 2n = 2x = 24), once known as the "King of the Appalachian Forest", was decimated by chestnut blight during the first half of the twentieth century by an invasive fungus (Cryphonectria parasitica). The Chinese chestnut (C. mollissima, 2n = 2x = 24), in contrast to American chestnut, is resistant to this blight. Efforts are being made to transfer this resistance to American chestnut through backcross breeding and genetic engineering. Both chestnut genomes have been genetically mapped and recently sequenced to facilitate gene discovery efforts aimed at assisting molecular breeding and genetic engineering. To complement and extend this genomic work, we analyzed the distribution and organization of their ribosomal DNAs (35S and 5S rDNA), and the chromatin composition of the nucleolus organizing region (NOR)-associated satellites. Using fluorescent in situ hybridization (FISH), we have identified two 35S (one major and one minor) and one 5S rDNA sites. The major 35S rDNA sites are terminal and sub-terminal in American and Chinese chestnuts, respectively, originating at the end of the short arm of the chromosome, extending through the secondary constriction and into the satellites. An additional 5S locus was identified in certain Chinese chestnut accessions, and it was linked distally to the major 35S site. The NOR-associated satellite in Chinese chestnut was found to comprise a proximal region packed with 35S rDNA and a distinct distal heterochromatic region. In contrast, the American chestnut satellite was relatively small and devoid of the distal heterochromatic region.

Comparative genomics of Cryptococcus and Kwoniella reveals pathogenesis evolution and contrasting karyotype dynamics via intercentromeric recombination or chromosome fusion

A large-scale comparative genomic analysis was conducted for the global human fungal pathogens within the Cryptococcus genus, compared to non-pathogenic Cryptococcus species, and related species from the sister genus Kwoniella. Chromosome-level genome assemblies were generated for multiple species of both genera, resulting in a dataset encompassing virtually all of their known diversity. Although Cryptococcus and Kwoniella have comparable genome sizes (about 19.2 and 22.9 Mb) and similar gene content, hinting at pre-adaptive pathogenic potential, our analysis found evidence in pathogenic Cryptococcus species of specific examples of gene gain (via horizontal gene transfer) and gene loss, which might represent evolutionary signatures of pathogenic development. Genome analysis also revealed a significant variation in chromosome number and structure between the two genera. By combining synteny analysis and experimental centromere validation, we found that most Cryptococcus species have 14 chromosomes, whereas most Kwoniella species have fewer (11, 8, 5 or even as few as 3). Reduced chromosome number in Kwoniella is associated with formation of giant chromosomes (up to 18 Mb) through repeated chromosome fusion events, each marked by a pericentric inversion and centromere loss. While similar chromosome inversion-fusion patterns were observed in all Kwoniella species with fewer than 14 chromosomes, no such pattern was detected in Cryptococcus. Instead, Cryptococcus species with less than 14 chromosomes, underwent chromosome reductions primarily through rearrangements associated with the loss of repeat-rich centromeres. Additionally, Cryptococcus genomes exhibited frequent interchromosomal translocations, including intercentromeric recombination facilitated by transposons shared between centromeres. Taken together, our findings advance our understanding of genomic changes possibly associated with pathogenicity in Cryptococcus and provide a foundation to elucidate mechanisms of centromere loss and chromosome fusion driving distinct karyotypes in closely related fungal species, including prominent global human pathogens.

Structural evidence for MADS-box type I family expansion seen in new assemblies of Arabidopsis arenosa and A. lyrata

Arabidopsis thaliana diverged from A. arenosa and A. lyrata at least 6 million years ago. The three species differ by genome-wide polymorphisms and morphological traits. The species are to a high degree reproductively isolated, but hybridization barriers are incomplete. A special type of hybridization barrier is based on the triploid endosperm of the seed, where embryo lethality is caused by endosperm failure to support the developing embryo. The MADS-box type I family of transcription factors is specifically expressed in the endosperm and has been proposed to play a role in endosperm-based hybridization barriers. The gene family is well known for its high evolutionary duplication rate, as well as being regulated by genomic imprinting. Here we address MADS-box type I gene family evolution and the role of type I genes in the context of hybridization. Using two de-novo assembled and annotated chromosome-level genomes of A. arenosa and A. lyrata ssp. petraea we analyzed the MADS-box type I gene family in Arabidopsis to predict orthologs, copy number, and structural genomic variation related to the type I loci. Our findings were compared to gene expression profiles sampled before and after the transition to endosperm cellularization in order to investigate the involvement of MADS-box type I loci in endosperm-based hybridization barriers. We observed substantial differences in type-I expression in the endosperm of A. arenosa and A. lyrata ssp. petraea, suggesting a genetic cause for the endosperm-based hybridization barrier between A. arenosa and A. lyrata ssp. petraea.

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.

High rates of structural rearrangements have shaped the chromosome evolution in dysploid Phaseolus beans

KEY MESSAGE

Karyotypes evolve through numerical and structural chromosome rearrangements. We show that Phaseolus leptostachyus, a wild bean, underwent a rapid genome reshuffling associated with the reduction from 11 to 10 chromosome pairs, but without whole genome duplication, the highest chromosome evolution rate known for plants. Plant karyotypes evolve through structural rearrangements often associated with polyploidy or dysploidy. The genus Phaseolus comprises ~ 90 species, five of them domesticated due to their nutritional relevance. Most of the species have 2n = 22 karyotypes and are highly syntenic, except for three dysploid karyotypes of species from the Leptostachyus group (2n = 20) that have accumulated several rearrangements. Here, we investigated the degrees of structural rearrangements among Leptostachyus and other Phaseolus groups by estimating their chromosomal evolution rates (CER). For this, we combined our oligo-FISH barcode system for beans and chromosome-specific painting probes for chromosomes 2 and 3, with rDNA and a centromeric probe to establish chromosome orthologies and identify structural rearrangements across nine Phaseolus species. We also integrated the detected rearrangements with a phylogenomic approach to estimate the CERs for each Phaseolus lineage. Our data allowed us to identify translocations, inversions, duplications and deletions, mostly in species belonging to the Leptostachyus group. Phaseolus leptostachyus showed the highest CER (12.31 rearrangements/My), a tenfold increase in contrast to the 2n = 22 species analysed. This is the highest rate known yet for plants, making it a model species for investigating the mechanisms behind rapid genome reshuffling in early species diversification.

Transition to Self-compatibility Associated With Dominant S-allele in a Diploid Siberian Progenitor of Allotetraploid Arabidopsis kamchatica Revealed by Arabidopsis lyrata Genomes

A transition to selfing can be beneficial when mating partners are scarce, for example, due to ploidy changes or at species range edges. Here, we explain how self-compatibility evolved in diploid Siberian Arabidopsis lyrata, and how it contributed to the establishment of allotetraploid Arabidopsis kamchatica. First, we provide chromosome-level genome assemblies for two self-fertilizing diploid A. lyrata accessions, one from North America and one from Siberia, including a fully assembled S-locus for the latter. We then propose a sequence of events leading to the loss of self-incompatibility in Siberian A. lyrata, date this independent transition to ∼90 Kya, and infer evolutionary relationships between Siberian and North American A. lyrata, showing an independent transition to selfing in Siberia. Finally, we provide evidence that this selfing Siberian A. lyrata lineage contributed to the formation of the allotetraploid A. kamchatica and propose that the selfing of the latter is mediated by the loss-of-function mutation in a dominant S-allele inherited from A. lyrata.

Short Interrupted Repeat Cassette (SIRC)-Novel Type of Repetitive DNA Element Found in Arabidopsis thaliana

Short interrupted repeat cassette (SIRC)-a novel DNA element found throughout the A. thaliana nuclear genome. SIRCs are represented by short direct repeats interrupted by diverse DNA sequences. The maxima of SIRC's distribution are located within pericentromeric regions. We suggest that originally SIRC was a special case of the complex internal structure of the miniature inverted repeat transposable element (MITE), and further MITE amplification, transposition, and loss of terminal inverted repeats gave rise to SIRC as an independent DNA element. SIRC sites were significantly enriched with several histone modifications associated with constitutive heterochromatin and mobile genetic elements. The majority of DNA-binding proteins, strongly associated with SIRC, are related to histone modifications for transcription repression. A part of SIRC was found to overlap highly inducible protein-coding genes, suggesting a possible regulatory role for these elements, yet their definitive functions need further investigation.

Genomic analyses provide insights into the polyploidization-driven herbicide adaptation in Leptochloa weeds

Polyploidy confers a selective advantage under stress conditions; however, whether polyploidization mediates enhanced herbicide adaptation remains largely unknown. Tetraploid Leptochloa chinensis is a notorious weed in the rice ecosystem, causing severe yield loss in rice. In China, L. chinensis has only one sister species, the diploid L. panicea, whose damage is rarely reported. To gain insights into the effects of polyploidization on herbicide adaptation, we first assembled a high-quality genome of L. panicea and identified genome structure variations with L. chinensis. Moreover, we identified herbicide-resistance genes specifically expanded in L. chinensis, which may confer a greater herbicide adaptability in L. chinensis. Analysis of gene retention and loss showed that five herbicide target-site genes and several herbicide nontarget-site resistance gene families were retained during polyploidization. Notably, we identified three pairs of polyploidization-retained genes including LcABCC8, LcCYP76C1 and LcCYP76C4 that may enhance herbicide resistance. More importantly, we found that both copies of LcCYP76C4 were under herbicide selection during the spread of L. chinensis in China. Furthermore, we identified another gene potentially involved in herbicide resistance, LcCYP709B2, which is also retained during polyploidization and under selection. This study provides insights into the genomic basis of the enhanced herbicide adaptability of Leptochloa weeds during polyploidization and provides guidance for the precise and efficient control of polyploidy weeds.

Genomic rearrangements and evolutionary changes in 3D chromatin topologies in the cotton tribe (Gossypieae)

BACKGROUND

Analysis of the relationship between chromosomal structural variation (synteny breaks) and 3D-chromatin architectural changes among closely related species has the potential to reveal causes and correlates between chromosomal change and chromatin remodeling. Of note, contrary to extensive studies in animal species, the pace and pattern of chromatin architectural changes following the speciation of plants remain unexplored; moreover, there is little exploration of the occurrence of synteny breaks in the context of multiple genome topological hierarchies within the same model species.

RESULTS

Here we used Hi-C and epigenomic analyses to characterize and compare the profiles of hierarchical chromatin architectural features in representative species of the cotton tribe (Gossypieae), including Gossypium arboreum, Gossypium raimondii, and Gossypioides kirkii, which differ with respect to chromosome rearrangements. We found that (i) overall chromatin architectural territories were preserved in Gossypioides and Gossypium, which was reflected in their similar intra-chromosomal contact patterns and spatial chromosomal distributions; (ii) the non-random preferential occurrence of synteny breaks in A compartment significantly associate with the B-to-A compartment switch in syntenic blocks flanking synteny breaks; (iii) synteny changes co-localize with open-chromatin boundaries of topologically associating domains, while TAD stabilization has a greater influence on regulating orthologous expression divergence than do rearrangements; and (iv) rearranged chromosome segments largely maintain ancestral in-cis interactions.

CONCLUSIONS

Our findings provide insights into the non-random occurrence of epigenomic remodeling relative to the genomic landscape and its evolutionary and functional connections to alterations of hierarchical chromatin architecture, on a known evolutionary timescale.

Chromosome Painting Using Chromosome-Specific BAC Clones

Chromosome painting (CP) refers to visualization of large chromosome regions, chromosome arms or entire chromosomes via fluorescence in situ hybridization (FISH) of chromosome-specific DNA sequences. For CP in crucifers (Brassicaceae), typically contigs of chromosome-specific bacterial artificial chromosomes (BAC) from Arabidopsis thaliana are applied as painting probes on chromosomes of A. thaliana or other species (comparative chromosome painting, CCP). CP/CCP enables to identify and trace particular chromosome regions and/or chromosomes throughout all mitotic and meiotic stages as well as corresponding interphase chromosome territories. However, extended pachytene chromosomes provide the highest resolution of CP/CCP. Fine-scale chromosome structure, structural chromosome rearrangements (such as inversions, translocations, centromere repositioning), and chromosome breakpoints can be investigated by CP/CCP. BAC DNA probes can be accompanied by other types of DNA probes, such as repetitive DNA, genomic DNA, or synthetic oligonucleotide probes. Here, we describe a robust step-by-step protocol of CP and CCP which proved to be efficient across the family Brassicaceae, but which is also applicable to other angiosperm families.

Reconstruction of karyotypic evolution in Saccharum spontaneum species by comparative oligo-FISH mapping

BACKGROUND

Karyotype dynamics driven by chromosomal rearrangements has long been considered as a fundamental question in the evolutionary genetics. Saccharum spontaneum, the most primitive and complex species in the genus Saccharum, has reportedly undergone at least two major chromosomal rearrangements, however, its karyotypic evolution remains unclear.

RESULTS

In this study, four representative accessions, i.e., hypothetical diploid sugarcane ancestor (sorghum, x = 10), Sa. spontaneum Np-X (x = 10, tetraploid), 2012-46 (x = 9, hexaploid) and AP85-441 (x = 8, tetraploid), were selected for karyotype evolution studies. A set of oligonucleotide (oligo)-based barcode probes was developed based on the sorghum genome, which allowed universal identification of all chromosomes from sorghum and Sa. spontaneum. By comparative FISH assays, we reconstructed the karyotype evolutionary history and discovered that although chromosomal rearrangements resulted in greater variation in relative lengths of some chromosomes, all chromosomes maintained a conserved metacentric structure. Additionally, we found that the barcode oligo probe was not applicable for chromosome identification in both Sa. robustum and Sa. officinarum species, suggesting that sorghum is more distantly related to Sa. robustum and Sa. officinarum compared with Sa. spontaneum species.

CONCLUSIONS

Our study demonstrated that the barcode oligo-FISH is an efficient tool for chromosome identification and karyotyping research, and expanded our understanding of the karyotypic and chromosomal evolution in the genus Saccharum.

Comparative cytogenomics reveals genome reshuffling and centromere repositioning in the legume tribe Phaseoleae

The tribe Phaseoleae includes several legume crops with assembled genomes. Comparative genomic studies have evidenced the preservation of large genomic blocks among legumes, although chromosome dynamics during Phaseoleae evolution has not been investigated. We conducted a comparative genomic analysis to define an informative genomic block (GB) system and to reconstruct the ancestral Phaseoleae karyotype (APK). We identified GBs based on the orthologous genes between Phaseolus vulgaris and Vigna unguiculata and searched for GBs in different genomes of the Phaseolinae (P. lunatus) and Glycininae (Amphicarpaea edgeworthii) subtribes and Spatholobus suberectus (sister to Phaseolinae and Glycininae), using Medicago truncatula as the outgroup. We also used oligo-FISH probes of two P. vulgaris chromosomes to paint the orthologous chromosomes of two non-sequenced Phaseolinae species. We inferred the APK as having n = 11 and 19 GBs (A to S), hypothesizing five chromosome fusions that reduced the ancestral legume karyotype to n = 11. We identified the rearrangements among the APK and the subtribes and species, with extensive centromere repositioning in Phaseolus. We also reconstructed the chromosome number reduction in S. suberectus. The development of the GB system and the proposed APK provide useful approaches for future comparative genomic analyses of legume species.

Chromosome Numbers and Genome Sizes of All 36 Duckweed Species (Lemnaceae)

Usually, chromosome sets (karyotypes) and genome sizes are rather stable for distinct species and therefore of diagnostic value for taxonomy. In combination with (cyto)genomics, both features provide essential cues for genome evolution and phylogenetic relationship studies within and between taxa above the species level. We present for the first time a survey on chromosome counts and genome size measurement for one or more accessions from all 36 duckweed species and discuss the evolutionary impact and peculiarities of both parameters in duckweeds.

Development and applications of a collection of single copy gene-based cytogenetic DNA markers in garden asparagus

Garden asparagus (Asparagus officinalis, 2n = 2x = 20 chromosomes) is an important dioecious vegetable crop and a model species for studying sex chromosome formation and evolution. However, few molecular cytogenetic studies on garden asparagus have been reported because of its small metaphase chromosomes, the scarcity of distinguished cytogenetic markers, and the high content of repetitive sequences. In this study, a set of single copy genes free of repetitive sequences with sizes ranging from 4.3 kb to 8.2 kb were screened and used as probes for fluorescence in situ hybridization (FISH) to identify individual chromosomes of garden asparagus. The chromosome-specific signal distribution patterns of these probes enabled the distinguishment of each pair of chromosomes. The sequence assembly and cytogenetic map were successfully integrated, and the results confirmed that the chromosome 1 representing the sex chromosome in the genome assembly is chromosome 5 in the karyotype analysis. The cytogenetic identification of the male-specific region of the Y chromosome (MSY) was implemented using a mixed probe derived from a number of MSY-specific single copy sequences. In addition, the chromosome orthologous relationship between garden asparagus (A1-A10, karyotypic analysis) and its hermaphrodite close relative, A. setaceus (B1-B10, karyotypic analysis), was analyzed using this collection of chromosome-specific cytological markers. The results showed that B3 is the ortholog of sex chromosome A5 and thus may represent the ancestral autosome of the current sex chromosome in garden asparagus. Chromosomes B5, B4, B1, B8, B7, and B9 are the orthologs of A2, A3, A4, A7, A8, and A10, respectively. The chromosome identification, cytogenetic recognition of MSY, and the orthologous relationship analysis between garden asparagus and A. setaceus are valuable for the further investigation of the sex chromosome emergence and evolutionary mechanism of garden asparagus and genome structure evolution in the Asparagus genus.

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.

Recent advancements in CRISPR/Cas technology for accelerated crop improvement

Precise genome engineering approaches could be perceived as a second paradigm for targeted trait improvement in crop plants, with the potential to overcome the constraints imposed by conventional CRISPR/Cas technology. The likelihood of reduced agricultural production due to highly turbulent climatic conditions increases as the global population expands. The second paradigm of stress-resilient crops with enhanced tolerance and increased productivity against various stresses is paramount to support global production and consumption equilibrium. Although traditional breeding approaches have substantially increased crop production and yield, effective strategies are anticipated to restore crop productivity even further in meeting the world's increasing food demands. CRISPR/Cas, which originated in prokaryotes, has surfaced as a coveted genome editing tool in recent decades, reshaping plant molecular biology in unprecedented ways and paving the way for engineering stress-tolerant crops. CRISPR/Cas is distinguished by its efficiency, high target specificity, and modularity, enables precise genetic modification of crop plants, allowing for the creation of allelic variations in the germplasm and the development of novel and more productive agricultural practices. Additionally, a slew of advanced biotechnologies premised on the CRISPR/Cas methodologies have augmented fundamental research and plant synthetic biology toolkits. Here, we describe gene editing tools, including CRISPR/Cas and its imitative tools, such as base and prime editing, multiplex genome editing, chromosome engineering followed by their implications in crop genetic improvement. Further, we comprehensively discuss the latest developments of CRISPR/Cas technology including CRISPR-mediated gene drive, tissue-specific genome editing, dCas9 mediated epigenetic modification and programmed self-elimination of transgenes in plants. Finally, we highlight the applicability and scope of advanced CRISPR-based techniques in crop genetic improvement.

Nonhomologous end joining as key to CRISPR/Cas-mediated plant chromosome engineering

Although clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein (Cas)-mediated gene editing has revolutionized biology and plant breeding, large-scale, heritable restructuring of plant chromosomes is still in its infancy. Duplications and inversions within a chromosome, and also translocations between chromosomes, can now be achieved. Subsequently, genetic linkages can be broken or can be newly created. Also, the order of genes on a chromosome can be changed. While natural chromosomal recombination occurs by homologous recombination during meiosis, CRISPR/Cas-mediated chromosomal rearrangements can be obtained best by harnessing nonhomologous end joining (NHEJ) pathways in somatic cells. NHEJ can be subdivided into the classical (cNHEJ) and alternative NHEJ (aNHEJ) pathways, which partially operate antagonistically. The cNHEJ pathway not only protects broken DNA ends from degradation but also suppresses the joining of previously unlinked broken ends. Hence, in the absence of cNHEJ, more inversions or translocations can be obtained which can be ascribed to the unrestricted use of the aNHEJ pathway for double-strand break (DSB) repair. In contrast to inversions or translocations, short tandem duplications can be produced by paired single-strand breaks via a Cas9 nickase. Interestingly, the cNHEJ pathway is essential for these kinds of duplications, whereas aNHEJ is required for patch insertions that can also be formed during DSB repair. As chromosome engineering has not only been accomplished in the model plant Arabidopsis (Arabidopsis thaliana) but also in the crop maize (Zea mays), we expect that this technology will soon transform the breeding process.

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.

The Ribosomal DNA Loci of the Ancient Monocot Pistia stratiotes L. (Araceae) Contain Different Variants of the 35S and 5S Ribosomal RNA Gene Units

The freshwater plant water lettuce (Pistia stratiotes L.) grows in warm climatic zones and is used for phytoremediation and biomass production. P. stratiotes belongs to the Araceae, an ecologically and structurally diverse early monocot family, but the phylogenetic relationships among Araceae members are poorly understood. Ribosomal DNAs (rDNAs), including the 35S and 5S rDNA, encode the RNA components of ribosomes and are widely used in phylogenetic and evolutionary studies of various plant taxa. Here, we comprehensively characterized the chromosomal locations and molecular organization of 35S and 5S rDNA genes in water lettuce using karyological and molecular methods. Fluorescence in situ hybridization revealed a single location for the 35S and 5S rDNA loci, each on a different pair of the species' 28 chromosomes. Molecular cloning and nucleotide sequencing of 35S rDNA of P. stratiotes, the first representative Araceae sensu stricto in which such a study was performed, displayed typical structural characteristics. The full-length repeat showed high sequence conservation of the regions producing the 18S, 5.8S, and 25S rRNAs and divergence of the internal transcribed spacers ITS1 and ITS2 as well as the large intergenic spacer (IGS). Alignments of the deduced sequence of 18S rDNA with the sequences available for other Araceae and representatives of other clades were used for phylogenetic analysis. Examination of 11 IGS sequences revealed significant intra-genomic length variability due to variation in subrepeat number, with four types of units detected within the 35S rDNA locus of the P. stratiotes genome (estimated size 407 Mb/1C). Similarly, the 5S rDNA locus harbors gene units comprising a conserved 119-bp sequence encoding 5S rRNA and two types of non-transcribed spacer (NTS) sequences. Type I was classified into four subtypes, which apparently originated via progressive loss of subrepeats within the duplicated NTS region containing the 3' part of the 5S rRNA gene. The minor Type II NTS is shorter than Type I and differs in nucleotide composition. Some DNA clones containing two or three consecutive 5S rDNA repeats harbored 5S rDNA genes with different types of NTSs, confirming the mosaic composition of the 5S rDNA locus.

Genome Size and Chromosome Number Evaluation of Astragalus L. sect. Hymenostegis Bunge (Fabaceae)

Astragalus section Hymenostegis is one of the important characteristic elements of thorn-cushion formations in the Irano-Turanian floristic region. In this paper, we examined the chromosome number of 17 species (15 new reports) and provide estimates of genome size for 62 individuals belonging to 38 taxa of A. sect. Hymenostegis, some species outside this section, plus two Oxytropis species. Based on chromosome counts 11 species were found to be diploid (2n = 16), four species tetraploid (2n = 32) and two taxa hexaploid (2n = 48). From genome size measurements on silica-gel dried material, three ploidy levels (2x, 4x and 6x) were inferred, with a majority of species being diploid. The 2C values reach from 2.07 pg in diploid Astragalus zohrabi to 7.16 pg in hexaploid A. rubrostriatus. We found indications that species might occur with different cytotypes. A phylogenetic framework using nrDNA ITS sequences was constructed to understand the evolution of ploidy changes and genome sizes. It showed that genome size values among the studied taxa differ only slightly within ploidy levels and are nearly constant within most species and groups of closely related taxa within the genus Astragalus. The results of this study show that there is a rather strong correlation between genome sizes and chromosome numbers in sect. Hymenostegis. The resolution of the ITS-based phylogenetic tree is too low to infer evolutionary or environmental correlations of genome size differences. Polyploidization seems to contribute to the high species number in Astragalus, however, in sect. Hymenostegis it is not the main driver of speciation.

Chromosome-specific painting unveils chromosomal fusions and distinct allopolyploid species in the Saccharum complex

Karyotypes provide key cytogenetic information on the phylogenetic relationships and evolutionary origins in related eukaryotic species. Despite our knowledge of the chromosome numbers of sugarcane and its wild relatives, the chromosome composition and evolution among the species in the Saccharum complex have been elusive owing to the complex polyploidy and the large numbers of chromosomes of these species. Oligonucleotide-based chromosome painting has become a powerful tool of cytogenetic studies especially for plant species with large numbers of chromosomes. We developed oligo-based chromosome painting probes for all 10 chromosomes in Saccharum officinarum (2n = 8x = 80). The 10 painting probes generated robust fluorescence in situ hybridization signals in all plant species within the Saccharum complex, including species in the genera Saccharum, Miscanthus, Narenga and Erianthus. We conducted comparative chromosome analysis using the same set of probes among species from four different genera within the Saccharum complex. Excitingly, we discovered several novel cytotypes and chromosome rearrangements in these species. We discovered that fusion from two different chromosomes is a common type of chromosome rearrangement associated with the species in the Saccharum complex. Such fusion events changed the basic chromosome number and resulted in distinct allopolyploids in the Saccharum complex.

A chromosome-level genome assembly of an alpine plant Crucihimalaya lasiocarpa provides insights into high-altitude adaptation

It remains largely unknown how plants adapt to high-altitude habitats. Crucihimalaya (Brassicaceae) is an alpine genus occurring in the Qinghai-Tibet Plateau characterized by cold temperatures and strong ultraviolet radiation. Here, we generated a chromosome-level genome for C. lasiocarpa with a total size of 255.8 Mb and a scaffold N50 size of 31.9 Mb. We first examined the karyotype origin of this species and found that the karyotype of five chromosomes resembled the ancestral karyotype of the Brassicaceae family, while the other three showed strong chromosomal structural variations. In combination with the rough genome sequence of another congener (C. himalaica), we found that the significantly expanded gene families and positively selected genes involved in alpine adaptation have occurred since the origin of this genus. Our new findings provide valuable information for the chromosomal karyotype evolution of Brassicaceae and investigations of high-altitude environment adaptation of the genus.

Plant Cytogenetics in the Micronuclei Investigation-The Past, Current Status, and Perspectives

Cytogenetic approaches play an essential role as a quick evaluation of the first genetic effects after mutagenic treatment. Although labor-intensive and time-consuming, they are essential for the analyses of cytotoxic and genotoxic effects in mutagenesis and environmental monitoring. Over the years, conventional cytogenetic analyses were a part of routine laboratory testing in plant genotoxicity. Among the methods that are used to study genotoxicity in plants, the micronucleus test particularly represents a significant force. Currently, cytogenetic techniques go beyond the simple detection of chromosome aberrations. The intensive development of molecular biology and the significantly improved microscopic visualization and evaluation methods constituted significant support to traditional cytogenetics. Over the past years, distinct approaches have allowed an understanding the mechanisms of formation, structure, and genetic activity of the micronuclei. Although there are many studies on this topic in humans and animals, knowledge in plants is significantly limited. This article provides a comprehensive overview of the current knowledge on micronuclei characteristics in plants. We pay particular attention to how the recent contemporary achievements have influenced the understanding of micronuclei in plant cells. Together with the current progress, we present the latest applications of the micronucleus test in mutagenesis and assess the state of the environment.

Karyotype and genome size variation in white-flowered Eranthis sect. Shibateranthis (Ranunculaceae)

Comparative karyomorphological analyses of six out of the eight white-flowered species of Eranthissect.Shibateranthis have been carried out. All studied specimens of E.byunsanensis, E.lobulata, E.pinnatifida, and E.stellata had a somatic chromosome number 2n = 16 with basic chromosome number x = 8. On the contrary, E.tanhoensis and E.sibirica had a basic chromosome number x = 7. The specimens of E.tanhoensis were diploid with 2n = 14, while the specimens of E.sibirica were polyploid with 2n = 42. Monoploid chromosome sets of the investigated diploid species had 4-5 metacentric chromosomes and 2-4 submetacentric/subtelocentric/acrocentric chromosomes. The highest level of interchromosomal asymmetry, estimated via CV_CL, was found in E.byunsanensis and E.pinnatifida. The highest levels of intrachromosomal asymmetry (M_CA) and heterogeneity in centromere position (CV_CI) were found in E.lobulata and E.byunsanensis, while E.sibirica had the most symmetric karyotype. A multivariate PCoA analysis of basic karyotype parameters (2n, x, THL, CV_CL, M_CA, and CV_CI) highlighted no overlap among species accessions, which was also confirmed by LDA. The average absolute monoploid DNA content (1Cx) of the 23 investigated samples of six Eranthis species varied from 9.26 ± 0.25 pg in E.sibirica to 15.93 ± 0.32 pg in E.stellata. Overall karyological affinity was highlighted between E.lobulata and E.stellata, on one side, and between E.byunsanensis and E.pinnatifida, on the other side. Interestingly, there was no significant correlation between total haploid (monoploid) chromosome length (THL) and 1Cx values in these species.

Chromosomal fusions facilitate adaptation to divergent environments in threespine stickleback

Chromosomal fusions are hypothesized to facilitate adaptation to divergent environments, both by bringing together previously unlinked adaptive alleles and by creating regions of low recombination that facilitate the linkage of adaptive alleles; but, there is little empirical evidence to support this hypothesis. Here, we address this knowledge gap by studying threespine stickleback (Gasterosteus aculeatus), in which ancestral marine fish have repeatedly adapted to freshwater across the northern hemisphere. By comparing the threespine and ninespine stickleback (Pungitius pungitius) genomes to a de novo assembly of the fourspine stickleback (Apeltes quadracus) and an outgroup species, we find two chromosomal fusion events involving the same chromosomes have occurred independently in the threespine and ninespine stickleback lineages. On the fused chromosomes in threespine stickleback, we find an enrichment of quantitative trait loci underlying traits that contribute to marine versus freshwater adaptation. By comparing whole-genome sequences of freshwater and marine threespine stickleback populations, we also find an enrichment of regions under divergent selection on these two fused chromosomes. There is elevated genetic diversity within regions under selection in the freshwater population, consistent with a simulation study showing that gene flow can increase diversity in genomic regions associated with local adaptation and our demographic models showing gene flow between the marine and freshwater populations. Integrating our results with previous studies, we propose that these fusions created regions of low recombination that enabled the formation of adaptative clusters, thereby facilitating freshwater adaptation in the face of recurrent gene flow between marine and freshwater threespine sticklebacks.

Genome structure and evolution in the cruciferous tribe Thlaspideae (Brassicaceae)

Whole-genome duplications (WGDs) and chromosome rearrangements (CRs) play the key role in driving the diversification and evolution of plant lineages. Although the direct link between WGDs and plant diversification is well documented, relatively few studies focus on the evolutionary significance of CRs. The cruciferous tribe Thlaspideae represents an ideal model system to address the role of large-scale chromosome alterations in genome evolution, as most Thlaspideae species share the same diploid chromosome number (2n = 2x = 14). Here we constructed the genome structure in 12 Thlaspideae species, including field pennycress (Thlaspi arvense) and garlic mustard (Alliaria petiolata). We detected and precisely characterized genus- and species-specific CRs, mostly pericentric inversions, as the main genome-diversifying drivers in the tribe. We reconstructed the structure of seven chromosomes of an ancestral Thlaspideae genome, identified evolutionary stable chromosomes versus chromosomes prone to CRs, estimated the rate of CRs, and uncovered an allohexaploid origin of garlic mustard from diploid taxa closely related to A. petiolata and Parlatoria cakiloidea. Furthermore, we performed detailed bioinformatic analysis of the Thlaspideae repeatomes, and identified repetitive elements applicable as unique species- and genus-specific barcodes and chromosome landmarks. This study deepens our general understanding of the evolutionary role of CRs, particularly pericentric inversions, in plant genome diversification, and provides a robust base for follow-up whole-genome sequencing efforts.

A genome sequence resource for the genus Passiflora, the genome of the wild diploid species Passiflora organensis

The genus Passiflora comprises a large group of plants popularly known as passionfruit, much appreciated for their exotic flowers and edible fruits. The species (∼500) are morphologically variable (e.g., growth habit, size, and color of flowers) and are adapted to distinct tropical ecosystems. In this study, we generated the genome of the wild diploid species Passiflora organensis Gardner by adopting a hybrid assembly approach. Passiflora organensis has a small genome of 259 Mbp and a heterozygosity rate of 81%, consistent with its reproductive system. Most of the genome sequences could be integrated into its chromosomes with cytogenomic markers (satellite DNA) as references. The repeated sequences accounted for 58.55% of the total DNA analyzed, and the Tekay lineage was the prevalent retrotransposon. In total, 25,327 coding genes were predicted. Passiflora organensis retains 5,609 singletons and 15,671 gene families. We focused on the genes potentially involved in the locus determining self-incompatibility and the MADS-box gene family, allowing us to infer expansions and contractions within specific subfamilies. Finally, we recovered the organellar DNA. Structural rearrangements and two mitoviruses, besides relics of other mobile elements, were found in the chloroplast and mt-DNA molecules, respectively. This study presents the first draft genome assembly of a wild Passiflora species, providing a valuable sequence resource for genomic and evolutionary studies on the genus, and support for breeding cropped passionfruit species.

In Silico Analysis of Fatty Acid Desaturases Structures in Camelina sativa, and Functional Evaluation of Csafad7 and Csafad8 on Seed Oil Formation and Seed Morphology

Fatty acid desaturases add a second bond into a single bond of carbon atoms in fatty acid chains, resulting in an unsaturated bond between the two carbons. They are classified into soluble and membrane-bound desaturases, according to their structure, subcellular location, and function. The orthologous genes in Camelina sativa were identified and analyzed, and a total of 62 desaturase genes were identified. It was revealed that they had the common fatty acid desaturase domain, which has evolved separately, and the proteins of the same family also originated from the same ancestry. A mix of conserved, gained, or lost intron structure was obvious. Besides, conserved histidine motifs were found in each family, and transmembrane domains were exclusively revealed in the membrane-bound desaturases. The expression profile analysis of C. sativa desaturases revealed an increase in young leaves, seeds, and flowers. C. sativa ω3-fatty acid desaturases CsaFAD7 and CsaDAF8 were cloned and the subcellular localization analysis showed their location in the chloroplast. They were transferred into Arabidopsis thaliana to obtain transgenic lines. It was revealed that the ω3-fatty acid desaturase could increase the C18:3 level at the expense of C18:2, but decreases in oil content and seed weight, and wrinkled phenotypes were observed in transgenic CsaFAD7 lines, while no significant change was observed in transgenic CsaFAD8 lines in comparison to the wild-type. These findings gave insights into the characteristics of desaturase genes, which could provide an excellent basis for further investigation for C. sativa improvement, and overexpression of ω3-fatty acid desaturases in seeds could be useful in genetic engineering strategies, which are aimed at modifying the fatty acid composition of seed oil.

Functional homoeologous alleles of CONSTANS contribute to seasonal crop type in rapeseed

Two CO paralogs in Brassica napus were confirmed and shown distinct expression pattern and function in promoting flowering and allelic variation s within BnaCO.A10 were found closely associated with ecotype divergence. CONSTANS (CO) is a key gene that responds to photoperiod and in Arabidopsis can promote flowering under long-day (LD) conditions. Brassica napus L. is a major oil crop and close relative of Arabidopsis, and arose via allopolyploidization from the diploids B. rapa (A genome) and B. oleracea (C genome). In this study, we confirmed that B. napus has two CO genes located on the A10 (BnaCO.A10) and C9 (BnaCO.C9) chromosomes. Significant differences in level and temporal pattern of transcription, as well as in protein function, of these homoeologous may have resulted from sequence variation in the promoter as well as in the coding region. Apart from two insertions of 527 bp and 2002 bp in the promoter of BnaCO.C9 that function as transcriptional enhancers, this gene is otherwise highly conserved in both promoter and coding region. However, BnaCO.A10 was classified into two haplotypes and transgene analysis in Arabidopsis and backcross analysis in rapeseed indicated that the winter-type haplotype had a greater effect in promoting flowering than the spring type. We discuss the contribution of CO alleles to species evolution, and for eco-geographic radiation following crop domestication, alongside scope for managing this locus in future breeding.

The evolution of the expansin gene family in Brassica species

Expansin gene (EXP) family plays important roles in plant growth and crop improvement. However, it has not been well studied in the Brassica genus that includes several important agricultural and horticultural crops. To get insight to the evolution and expansion of EXP family in Brassica, Brassica EXPs which are homologues of 35 known AtEXPs of Arabidopsis were comprehensively and systematically analyzed in the present study. In total, 340 Brassica EXPs were clustered into four groups that corresponded multiple alignment to four subfamilies of AtEXPs, with divergent conserved motifs and cis-acting elements among groups. To understand the expansion of EXP family, an integrated genomic block system was constructed among Arabidopsis and Brassica species based on 24 known ancestral karyotype blocks. Obvious gene loss, segmental duplication, tandem duplication and DNA sequence repeat events were found during the expansion of Brassica EXPs, of which the segmental duplication was possibly the major driving force. The divergence time was estimated in 1109 orthologs pairs of EXPs, revealing the divergence of Brassica EXPs from AtEXPs during ~30 MYA, and the divergence of EXPs among Brassica species during 13.50-17.94 MYA. Selective mode analysis revealed that the purifying selection was the major contributor to expansion of Brassica EXPs. This study provides new insights into the evolution and expansion of the EXP family in Brassica genus.

Descending Dysploidy and Bidirectional Changes in Genome Size Accompanied Crepis (Asteraceae) Evolution

The evolution of the karyotype and genome size was examined in species of Crepis sensu lato. The phylogenetic relationships, inferred from the plastid and nrITS DNA sequences, were used as a framework to infer the patterns of karyotype evolution. Five different base chromosome numbers (x = 3, 4, 5, 6, and 11) were observed. A phylogenetic analysis of the evolution of the chromosome numbers allowed the inference of x = 6 as the ancestral state and the descending dysploidy as the major direction of the chromosome base number evolution. The derived base chromosome numbers (x = 5, 4, and 3) were found to have originated independently and recurrently in the different lineages of the genus. A few independent events of increases in karyotype asymmetry were inferred to have accompanied the karyotype evolution in Crepis. The genome sizes of 33 Crepis species differed seven-fold and the ancestral genome size was reconstructed to be 1C = 3.44 pg. Both decreases and increases in the genome size were inferred to have occurred within and between the lineages. The data suggest that, in addition to dysploidy, the amplification/elimination of various repetitive DNAs was likely involved in the genome and taxa differentiation in the genus.

Reconstruction of ancestral karyotype illuminates chromosome evolution in the genus Cucumis

Karyotype dynamics driven by complex chromosome rearrangements constitute a fundamental issue in evolutionary genetics. The evolutionary events underlying karyotype diversity within plant genera, however, have rarely been reconstructed from a computed ancestral progenitor. Here, we developed a method to rapidly and accurately represent extant karyotypes with the genus, Cucumis, using highly customizable comparative oligo-painting (COP) allowing visualization of fine-scale genome structures of eight Cucumis species from both African-origin and Asian-origin clades. Based on COP data, an evolutionary framework containing a genus-level ancestral karyotype was reconstructed, allowing elucidation of the evolutionary events that account for the origin of these diverse genomes within Cucumis. Our results characterize the cryptic rearrangement hotspots on ancestral chromosomes, and demonstrate that the ancestral Cucumis karyotype (n = 12) evolved to extant Cucumis genomes by hybridizations and frequent lineage- and species-specific genome reshuffling. Relative to the African species, the Asian species, including melon (Cucumis melo, n = 12), Cucumis hystrix (n = 12) and cucumber (Cucumis sativus, n = 7), had highly shuffled genomes caused by large-scale inversions, centromere repositioning and chromothripsis-like rearrangement. The deduced reconstructed ancestral karyotype for the genus allowed us to propose evolutionary trajectories and specific events underlying the origin of these Cucumis species. Our findings highlight that the partitioned evolutionary plasticity of Cucumis karyotype is primarily located in the centromere-proximal regions marked by rearrangement hotspots, which can potentially serve as a reservoir for chromosome evolution due to their fragility.

Chasing breeding footprints through structural variations in Cucumis melo and wild relatives

Cucumis melo (melon or muskmelon) is an important crop in the family of the Cucurbitaceae. Melon is cross pollinated and domesticated at several locations throughout the breeding history, resulting in highly diverse genetic structure in the germplasm. Yet, the relations among the groups and cultivars are still incomplete. We shed light on the melonbreeding history, analyzing structural variations ranging from 50 bp up to 100 kb, identified from whole genome sequences of 100 selected melon accessions and wild relatives. Phylogenetic trees based on SV types completely resolve cultivars and wild accessions into two monophyletic groups and clustering of cultivars largely correlates with their geographic origin. Taking into account morphology, we found six mis-categorized cultivars. Unique inversions are more often shared between cultivars, carrying advantageous genes and do not directly originate from wild species. Approximately 60% of the inversion breaks carry a long poly A/T motif, and following observations in other plant species, suggest that inversions in melon likely resulted from meiotic recombination events. We show that resistance genes in the linkage V region are expanded in the cultivar genomes compared to wild relatives. Furthermore, particular agronomic traits such as fruit ripening, fragrance, and stress response are specifically selected for in the melon subspecies. These results represent distinctive footprints of selective breeding that shaped today's melon. The sequences and genomic relations between land races, wild relatives, and cultivars will serve the community to identify genetic diversity, optimize experimental designs, and enhance crop development.

BAC- and oligo-FISH mapping reveals chromosome evolution among Vigna angularis, V. unguiculata, and Phaseolus vulgaris

Cytogenomic resources have accelerated synteny and chromosome evolution studies in plant species, including legumes. Here, we established the first cytogenetic map of V. angularis (Va, subgenus Ceratotropis) and compared this new map with those of V. unguiculata (Vu, subgenus Vigna) and P. vulgaris (Pv) by BAC-FISH and oligopainting approaches. We mapped 19 Vu BACs and 35S rDNA probes to the 11 chromosome pairs of Va, Vu, and Pv. Vigna angularis shared a high degree of macrosynteny with Vu and Pv, with five conserved syntenic chromosomes. Additionally, we developed two oligo probes (Pv2 and Pv3) used to paint Vigna orthologous chromosomes. We confirmed two reciprocal translocations (chromosomes 2 and 3 and 1 and 8) that have occurred after the Vigna and Phaseolus divergence (~9.7 Mya). Besides, two inversions (2 and 4) and one translocation (1 and 5) have occurred after Vigna and Ceratotropis subgenera separation (~3.6 Mya). We also observed distinct oligopainting patterns for chromosomes 2 and 3 of Vigna species. Both Vigna species shared similar major rearrangements compared to Pv: one translocation (2 and 3) and one inversion (chromosome 3). The sequence synteny identified additional inversions and/or intrachromosomal translocations involving pericentromeric regions of both orthologous chromosomes. We propose chromosomes 2 and 3 as hotspots for chromosomal rearrangements and de novo centromere formation within and between Vigna and Phaseolus. Our BAC- and oligo-FISH mapping contributed to physically trace the chromosome evolution of Vigna and Phaseolus and its application in further studies of both genera.

Tempo and mode in karyotype evolution revealed by a probabilistic model incorporating both chromosome number and morphology

Karyotype, including the chromosome and arm numbers, is a fundamental genetic characteristic of all organisms and has long been used as a species-diagnostic character. Additionally, karyotype evolution plays an important role in divergent adaptation and speciation. Centric fusion and fission change chromosome numbers, whereas the intra-chromosomal movement of the centromere, such as pericentric inversion, changes arm numbers. A probabilistic model simultaneously incorporating both chromosome and arm numbers has not been established. Here, we built a probabilistic model of karyotype evolution based on the "karyograph", which treats karyotype evolution as a walk on the two-dimensional space representing the chromosome and arm numbers. This model enables analysis of the stationary distribution with a stable karyotype for any given parameter. After evaluating their performance using simulated data, we applied our model to two large taxonomic groups of fish, Eurypterygii and series Otophysi, to perform maximum likelihood estimation of the transition rates and reconstruct the evolutionary history of karyotypes. The two taxa significantly differed in the evolution of arm number. The inclusion of speciation and extinction rates demonstrated possibly high extinction rates in species with karyotypes other than the most typical karyotype in both groups. Finally, we made a model including polyploidization rates and applied it to a small plant group. Thus, the use of this probabilistic model can contribute to a better understanding of tempo and mode in karyotype evolution and its possible role in speciation and extinction.

CRISPR/Cas-mediated chromosome engineering: opening up a new avenue for plant breeding

The advent of powerful site-specific nucleases, particularly the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein (Cas) system, which enables precise genome manipulation, has revolutionized plant breeding. Until recently, the main focus of researchers has been to simply knock-in or knock-out single genes, or to induce single base changes, but constant improvements of this technology have enabled more ambitious applications that aim to improve plant productivity or other desirable traits. One long-standing aim has been the induction of targeted chromosomal rearrangements (crossovers, inversions, or translocations). The feasibility of this technique has the potential to transform plant breeding, because natural rearrangements, like inversions, for example, typically present obstacles to the breeding process. In this way, genetic linkages between traits could be altered to combine or separate favorable and deleterious genes, respectively. In this review, we discuss recent breakthroughs in the field of chromosome engineering in plants and their potential applications in the field of plant breeding. In the future, these approaches might be applicable in shaping plant chromosomes in a directed manner, based on plant breeding needs.

Limitation of current probe design for oligo-cross-FISH, exemplified by chromosome evolution studies in duckweeds

Duckweeds represent a small, free-floating aquatic family (Lemnaceae) of the monocot order Alismatales with the fastest growth rate among flowering plants. They comprise five genera (Spirodela, Landoltia, Lemna, Wolffiella, and Wolffia) varying in genome size and chromosome number. Spirodela polyrhiza had the first sequenced duckweed genome. Cytogenetic maps are available for both species of the genus Spirodela (S. polyrhiza and S. intermedia). However, elucidation of chromosome homeology and evolutionary chromosome rearrangements by cross-FISH using Spirodela BAC probes to species of other duckweed genera has not been successful so far. We investigated the potential of chromosome-specific oligo-FISH probes to address these topics. We designed oligo-FISH probes specific for one S. intermedia and one S. polyrhiza chromosome (Fig. 1a). Our results show that these oligo-probes cross-hybridize with the homeologous regions of the other congeneric species, but are not suitable to uncover chromosomal homeology across duckweeds genera. This is most likely due to too low sequence similarity between the investigated genera and/or too low probe density on the target genomes. Finally, we suggest genus-specific design of oligo-probes to elucidate chromosome evolution across duckweed genera.

Chromosomal Mapping of Tandem Repeats Revealed Massive Chromosomal Rearrangements and Insights Into Senna tora Dysploidy

Tandem repeats can occupy a large portion of plant genomes and can either cause or result from chromosomal rearrangements, which are important drivers of dysploidy-mediated karyotype evolution and speciation. To understand the contribution of tandem repeats in shaping the extant Senna tora dysploid karyotype, we analyzed the composition and abundance of tandem repeats in the S. tora genome and compared the chromosomal distribution of these repeats between S. tora and a closely related euploid, Senna occidentalis. Using a read clustering algorithm, we identified the major S. tora tandem repeats and visualized their chromosomal distribution by fluorescence in situ hybridization. We identified eight independent repeats covering ~85 Mb or ~12% of the S. tora genome. The unit lengths and copy numbers had ranges of 7-5,833 bp and 325-2.89 × 10⁶, respectively. Three short duplicated sequences were found in the 45S rDNA intergenic spacer, one of which was also detected at an extra-NOR locus. The canonical plant telomeric repeat (TTTAGGG)_n was also detected as very intense signals in numerous pericentromeric and interstitial loci. StoTR05_180, which showed subtelomeric distribution in Senna occidentalis, was predominantly pericentromeric in S. tora. The unusual chromosomal distribution of tandem repeats in S. tora not only enabled easy identification of individual chromosomes but also revealed the massive chromosomal rearrangements that have likely played important roles in shaping its dysploid karyotype.

Identification of Chromosomes and Chromosome Rearrangements in Crop Brassicas and Raphanus sativus: A Cytogenetic Toolkit Using Synthesized Massive Oligonucleotide Libraries

Crop brassicas include three diploid [Brassica rapa (AA; 2n = 2x = 16), B. nigra (BB; 2n = 2x = 18), and B. oleracea (CC; 2n = 2x = 20)] and three derived allotetraploid species. It is difficult to distinguish Brassica chromosomes as they are small and morphologically similar. We aimed to develop a genome-sequence based cytogenetic toolkit for reproducible identification of Brassica chromosomes and their structural variations. A bioinformatic pipeline was used to extract repeat-free sequences from the whole genome assembly of B. rapa. Identified sequences were subsequently used to develop four c. 47-mer oligonucleotide libraries comprising 27,100, 11,084, 9,291, and 16,312 oligonucleotides. We selected these oligonucleotides after removing repeats from 18 identified sites (500-1,000 kb) with 1,997-5,420 oligonucleotides localized at each site in B. rapa. For one set of probes, a new method for amplification or immortalization of the library is described. oligonucleotide probes produced specific and reproducible in situ hybridization patterns for all chromosomes belonging to A, B, C, and R (Raphanus sativus) genomes. The probes were able to identify structural changes between the genomes, including translocations, fusions, and deletions. Furthermore, the probes were able to identify a structural translocation between a pak choi and turnip cultivar of B. rapa. Overall, the comparative chromosomal mapping helps understand the role of chromosome structural changes during genome evolution and speciation in the family Brassicaceae. The probes can also be used to identify chromosomes in aneuploids such as addition lines used for gene mapping, and to track transfer of chromosomes in hybridization and breeding programs.

Independent recurrent evolution of MICRORNA genes converging onto similar non-canonical organisation across green plant lineages is driven by local and segmental duplication events in species, family and lineages

The relationship between evolutionary history, organisation and transcriptional regulation of genes are intrinsically linked. These have been well studied in canonically organised protein-coding genes but not of MIRNAs. In the present study, we investigated the non-canonical arrangement of MIRNAs across taxonomic boundaries from algae to angiosperms employing a combination of genome organization, phylogeny and synteny. We retrieved the complete dataset of MIRNA from twenty-five species to identify and classify based on organisational patterns. The median size of cluster was between 2-5 kb and between 1-20 % of all MIRNAs are organized in head-to-head (with bidirectional promoter), head-to-tail (tandem), and overlapping manner. Although majority of the clusters are composed of MIRNA homologs, 25% of all clusters comprises of non-homologous genes with a potential of generating functional and regulatory complexity. A comparison of phylogeny and organizational patterns revealed that multiple independent events, some of which are species-specific, and some ancient, in different lineages, are responsible for non-canonical organization. Detailed investigation of MIR395 family across the plants revealed a complex origin of non-canonical arrangement through ancient and recent, segmental and local duplications; analysis of MIR399 family revealed major expansion occurred prior to monocot-dicot split, with few lineage-specific events. Evolution of "convergent" organization pattern of non-canonical arrangement originating from independent loci through recurrent event highlights our poor understanding of evolutionary process of MIRNA genes. The present investigation thus paves way for comparative functional genomics to understand the role of non-canonical organization on transcriptional regulation and regulatory diversity in MIRNA gene families.

Multiple and independent rearrangements revealed by comparative cytogenetic mapping in the dysploid Leptostachyus group (Phaseolus L., Leguminosae)

Polyploidy and dysploidy have been reported as the main events in karyotype evolution of plants. In the genus Phaseolus L. (2n = 22), a small monophyletic group of three species, the Leptostachyus group, presents a dysploid karyotype with 2n = 20. It was shown in Phaseolus leptostachyus that the dysploidy was caused by a nested chromosome fusion (NCF) accompanied by several translocations, suggesting a high rate of karyotype evolution in the group. To verify if this karyotype restructuring was a single event or occurred progressively during the evolution of this group, we analysed P. macvaughii, sister to Phaseolus micranthus + P. leptostachyus. Twenty-four genomic clones of P. vulgaris previously mapped on P. leptostachyus, in addition to the 5S and 35S rDNA probes, were used for fluorescence in situ hybridization. Only a single rearrangement was common to the two species: the nested chromosome fusion (NCF) involving chromosomes 10 and 11. The translocation of chromosome 2 is not the same found in P. leptostachyus, and pericentric inversions in chromosomed 3 and 4 were exclusive of P. macvaughii. The other rearrangements observed in P. leptostachyus were not shared with this species, suggesting that they occurred after the separation of these lineages. The presence of private rearrangements indicates a progressive accumulation of karyotype changes in the Leptostachyus group instead of an instant genome-wide repatterning.

Chromosome-scale genome assembly for the duckweed Spirodela intermedia, integrating cytogenetic maps, PacBio and Oxford Nanopore libraries

Duckweeds are small, free-floating, morphologically highly reduced organisms belonging to the monocot order Alismatales. They display the most rapid growth among flowering plants, vary ~ 14-fold in genome size and comprise five genera. Spirodela is the phylogenetically oldest genus with only two mainly asexually propagating species: S. polyrhiza (2n = 40; 160 Mbp/1C) and S. intermedia (2n = 36; 160 Mbp/1C). This study combined comparative cytogenetics and de novo genome assembly based on PacBio, Illumina and Oxford Nanopore (ON) reads to obtain the first genome reference for S. intermedia and to compare its genomic features with those of the sister species S. polyrhiza. Both species' genomes revealed little more than 20,000 putative protein-coding genes, very low rDNA copy numbers and a low amount of repetitive sequences, mainly Ty3/gypsy retroelements. The detection of a few new small chromosome rearrangements between both Spirodela species refined the karyotype and the chromosomal sequence assignment for S. intermedia.

The Evolution of Chromosome Numbers: Mechanistic Models and Experimental Approaches

Chromosome numbers have been widely used to describe the most fundamental genomic attribute of an organism or a lineage. Although providing strong phylogenetic signal, chromosome numbers vary remarkably among eukaryotes at all levels of taxonomic resolution. Changes in chromosome numbers regularly serve as indication of major genomic events, most notably polyploidy and dysploidy. Here, we review recent advancements in our ability to make inferences regarding historical events that led to alterations in the number of chromosomes of a lineage. We first describe the mechanistic processes underlying changes in chromosome numbers, focusing on structural chromosomal rearrangements. Then, we focus on experimental procedures, encompassing comparative cytogenomics and genomics approaches, and on computational methodologies that are based on explicit models of chromosome-number evolution. Together, these tools offer valuable predictions regarding historical events that have changed chromosome numbers and genome structures, as well as their phylogenetic and temporal placements.

Rapid Birth or Death of Centromeres on Fragmented Chromosomes in Maize

Comparative genomics has revealed common occurrences in karyotype evolution such as chromosomal end-to-end fusions and insertions of one chromosome into another near the centromere, as well as many cases of de novo centromeres that generate positional polymorphisms. However, how rearrangements such as dicentrics and acentrics persist without being destroyed or lost remains unclear. Here, we sought experimental evidence for the frequency and timeframe for inactivation and de novo formation of centromeres in maize (Zea mays). The pollen from plants with supernumerary B chromosomes was gamma-irradiated and then applied to normal maize silks of a line without B chromosomes. In ∼8,000 first-generation seedlings, we found many B-A translocations, centromere expansions, and ring chromosomes. We also found many dicentric chromosomes, but a fraction of these show only a single primary constriction, which suggests inactivation of one centromere. Chromosomal fragments were found without canonical centromere sequences, revealing de novo centromere formation over unique sequences; these were validated by immunolocalization with Thr133-phosphorylated histone H2A, a marker of active centromeres, and chromatin immunoprecipitation-sequencing with the CENH3 antibody. These results illustrate the regular occurrence of centromere birth and death after chromosomal rearrangement during a narrow window of one to potentially only a few cell cycles for the rearranged chromosomes to be recognized in this experimental regime.