Repeat-based holocentromeres influence genome architecture and karyotype evolution.

The centromere represents a single region in most eukaryotic chromosomes. However, several plant and animal lineages assemble holocentromeres along the entire chromosome length. Here, we compare genome organization and evolution as a function of centromere type by assembling chromosome-scale holocentric genomes with repeat-based holocentromeres from three beak-sedge (Rhynchospora pubera, R. breviuscula, and R. tenuis) and their closest monocentric relative, Juncus effusus. We demonstrate that transition to holocentricity affected 3D genome architecture by redefining genomic compartments, while distributing centromere function to thousands of repeat-based centromere units genome-wide. We uncover a complex genome organization in R. pubera that hides its unexpected octoploidy and describe a marked reduction in chromosome number for R. tenuis, which has only two chromosomes. We show that chromosome fusions, facilitated by repeat-based holocentromeres, promoted karyotype evolution and diploidization. Our study thus sheds light on several important aspects of genome architecture and evolution influenced by centromere organization.

The giant diploid faba genome unlocks variation in a global protein crop.

Increasing the proportion of locally produced plant protein in currently meat-rich diets could substantially reduce greenhouse gas emissions and loss of biodiversity1. However, plant protein production is hampered by the lack of a cool-season legume equivalent to soybean in agronomic value2. Faba bean (Vicia faba L.) has a high yield potential and is well suited for cultivation in temperate regions, but genomic resources are scarce. Here, we report a high-quality chromosome-scale assembly of the faba bean genome and show that it has expanded to a massive 13 Gb in size through an imbalance between the rates of amplification and elimination of retrotransposons and satellite repeats. Genes and recombination events are evenly dispersed across chromosomes and the gene space is remarkably compact considering the genome size, although with substantial copy number variation driven by tandem duplication. Demonstrating practical application of the genome sequence, we develop a targeted genotyping assay and use high-resolution genome-wide association analysis to dissect the genetic basis of seed size and hilum colour. The resources presented constitute a genomics-based breeding platform for faba bean, enabling breeders and geneticists to accelerate the improvement of sustainable protein production across the Mediterranean, subtropical and northern temperate agroecological zones.

Assembly of the 81.6 Mb centromere of pea chromosome 6 elucidates the structure and evolution of metapolycentric chromosomes.

Centromeres in the legume genera Pisum and Lathyrus exhibit unique morphological characteristics, including extended primary constrictions and multiple separate domains of centromeric chromatin. These so-called metapolycentromeres resemble an intermediate form between monocentric and holocentric types, and therefore provide a great opportunity for studying the transitions between different types of centromere organizations. However, because of the exceedingly large and highly repetitive nature of metapolycentromeres, highly contiguous assemblies needed for these studies are lacking. Here, we report on the assembly and analysis of a 177.6 Mb region of pea (Pisum sativum) chromosome 6, including the 81.6 Mb centromere region (CEN6) and adjacent chromosome arms. Genes, DNA methylation profiles, and most of the repeats were uniformly distributed within the centromere, and their densities in CEN6 and chromosome arms were similar. The exception was an accumulation of satellite DNA in CEN6, where it formed multiple arrays up to 2 Mb in length. Centromeric chromatin, characterized by the presence of the CENH3 protein, was predominantly associated with arrays of three different satellite repeats; however, five other satellites present in CEN6 lacked CENH3. The presence of CENH3 chromatin was found to determine the spatial distribution of the respective satellites during the cell cycle. Finally, oligo-FISH painting experiments, performed using probes specifically designed to label the genomic regions corresponding to CEN6 in Pisum, Lathyrus, and Vicia species, revealed that metapolycentromeres evolved via the expansion of centromeric chromatin into neighboring chromosomal regions and the accumulation of novel satellite repeats. However, in some of these species, centromere evolution also involved chromosomal translocations and centromere repositioning.

Disruption of the standard kinetochore in holocentric Cuscuta species.

The segregation of chromosomes depends on the centromere. Most species are monocentric, with the centromere restricted to a single region per chromosome. In some organisms, the monocentric organization changed to holocentric, in which the centromere activity is distributed over the entire chromosome length. However, the causes and consequences of this transition are poorly understood. Here, we show that the transition in the genus Cuscuta was associated with dramatic changes in the kinetochore, a protein complex that mediates the attachment of chromosomes to microtubules. We found that in holocentric Cuscuta species, the KNL2 genes were lost; the CENP-C, KNL1, and ZWINT1 genes were truncated; the centromeric localization of CENH3, CENP-C, KNL1, MIS12, and NDC80 proteins was disrupted; and the spindle assembly checkpoint (SAC) degenerated. Our results demonstrate that holocentric Cuscuta species lost the ability to form a standard kinetochore and do not employ SAC to control the attachment of microtubules to chromosomes.

Celine, a long interspersed nuclear element retrotransposon, colonizes in the centromeres of poplar chromosomes.

Centromeres in most multicellular eukaryotes are composed of long arrays of repetitive DNA sequences. Interestingly, several transposable elements, including the well-known long terminal repeat (LTR) retrotransposon CRM (centromeric retrotransposon of maize), were found to be enriched in functional centromeres marked by the centromeric histone H3 (CENH3). Here we report a centromeric long interspersed nuclear element (LINE), Celine, in Populus species. Celine has colonized preferentially in the CENH3-associated chromatin of every poplar chromosome, with 84% of the Celine elements localized in the CENH3-binding domains. By contrast, only 51% of the CRM elements were bound to CENH3 domains in Populus trichocarpa. These results suggest different centromere targeting mechanisms employed by Celine and CRM elements. Nevertheless, the high target specificity seems to be detrimental to further amplification of the Celine elements, leading to a shorter life span and patchy distribution among plant species compared to the CRM elements. Using a phylogenetically guided approach we were able to identify Celine-like LINE elements in tea plant (Camellia sinensis) and green ash tree (Fraxinus pennsylvanica). The centromeric localization of these Celine-like LINEs was confirmed in both species. We demonstrate that the centromere targeting property of Celine-like LINEs is of primitive origin and has been conserved among distantly related plant species.

KNL1 and NDC80 represent new universal markers for the detection of functional centromeres in plants.

Centromere is the chromosomal site of kinetochore assembly and microtubule attachment for chromosome segregation. Given its importance, markers that allow specific labeling of centromeric chromatin throughout the cell cycle and across all chromosome types are sought for facilitating various centromere studies. Antibodies against the N-terminal region of CENH3 are commonly used for this purpose, since CENH3 is the near-universal marker of functional centromeres. However, because the N-terminal region of CENH3 is highly variable among plant species, antibodies directed against this region usually function only in a small group of closely related species. As a more versatile alternative, we present here antibodies targeted to the conserved domains of two outer kinetochore proteins, KNL1 and NDC80. Sequence comparison of these domains across more than 350 plant species revealed a high degree of conservation, particularly within a six amino acid motif, FFGPVS in KNL1, suggesting that both antibodies would function in a wide range of plant species. This assumption was confirmed by immunolabeling experiments in angiosperm (monocot and dicot) and gymnosperm species, including those with mono-, holo-, and meta-polycentric chromosomes. In addition to centromere labeling on condensed chromosomes during cell division, both antibodies detected the corresponding regions in the interphase nuclei of most species tested. These results demonstrated that KNL1 and NDC80 are better suited for immunolabeling centromeres than CENH3, because antibodies against these proteins offer incomparably greater versatility across different plant species which is particularly convenient for studying the organization and function of the centromere in non-model species.

Characterization of repeat arrays in ultra-long nanopore reads reveals frequent origin of satellite DNA from retrotransposon-derived tandem repeats.

Amplification of monomer sequences into long contiguous arrays is the main feature distinguishing satellite DNA from other tandem repeats, yet it is also the main obstacle in its investigation because these arrays are in principle difficult to assemble. Here we explore an alternative, assembly-free approach that utilizes ultra-long Oxford Nanopore reads to infer the length distribution of satellite repeat arrays, their association with other repeats and the prevailing sequence periodicities. Using the satellite DNA-rich legume plant Lathyrus sativus as a model, we demonstrated this approach by analyzing 11 major satellite repeats using a set of nanopore reads ranging from 30 to over 200 kb in length and representing 0.73× genome coverage. We found surprising differences between the analyzed repeats because only two of them were predominantly organized in long arrays typical for satellite DNA. The remaining nine satellites were found to be derived from short tandem arrays located within LTR-retrotransposons that occasionally expanded in length. While the corresponding LTR-retrotransposons were dispersed across the genome, this array expansion occurred mainly in the primary constrictions of the L. sativus chromosomes, which suggests that these genome regions are favourable for satellite DNA accumulation.

Extraordinary Sequence Diversity and Promiscuity of Centromeric Satellites in the Legume Tribe Fabeae.

Satellite repeats are major sequence constituents of centromeres in many plant and animal species. Within a species, a single family of satellite sequences typically occupies centromeres of all chromosomes and is absent from other parts of the genome. Due to their common origin, sequence similarities exist among the centromere-specific satellites in related species. Here, we report a remarkably different pattern of centromere evolution in the plant tribe Fabeae, which includes genera Pisum, Lathyrus, Vicia, and Lens. By immunoprecipitation of centromeric chromatin with CENH3 antibodies, we identified and characterized a large and diverse set of 64 families of centromeric satellites in 14 species. These families differed in their nucleotide sequence, monomer length (33-2,979 bp), and abundance in individual species. Most families were species-specific, and most species possessed multiple (2-12) satellites in their centromeres. Some of the repeats that were shared by several species exhibited promiscuous patterns of centromere association, being located within CENH3 chromatin in some species, but apart from the centromeres in others. Moreover, FISH experiments revealed that the same family could assume centromeric and noncentromeric positions even within a single species. Taken together, these findings suggest that Fabeae centromeres are not shaped by the coevolution of a single centromeric satellite with its interacting CENH3 proteins, as proposed by the centromere drive model. This conclusion is also supported by the absence of pervasive adaptive evolution of CENH3 sequences retrieved from Fabeae species.

Impact of parasitic lifestyle and different types of centromere organization on chromosome and genome evolution in the plant genus Cuscuta.

The parasitic genus Cuscuta (Convolvulaceae) is exceptional among plants with respect to centromere organization, including both monocentric and holocentric chromosomes, and substantial variation in genome size and chromosome number. We investigated 12 species representing the diversity of the genus in a phylogenetic context to reveal the molecular and evolutionary processes leading to diversification of their genomes. We measured genome sizes and investigated karyotypes and centromere organization using molecular cytogenetic techniques. We also performed low-pass whole genome sequencing and comparative analysis of repetitive DNA composition. A remarkable 102-fold variation in genome sizes (342-34 734 Mbp/1C) was detected for monocentric Cuscuta species, while genomes of holocentric species were of moderate sizes (533-1545 Mbp/1C). The genome size variation was primarily driven by the differential accumulation of LTR-retrotransposons and satellite DNA. The transition to holocentric chromosomes in the subgenus Cuscuta was associated with loss of histone H2A phosphorylation and elimination of centromeric retrotransposons. In addition, basic chromosome number of holocentric species (x = 7) was smaller than in monocentrics (x = 15 or 16). We demonstrated that the transition to holocentricity in Cuscuta was accompanied by significant changes in epigenetic marks, chromosome number and the repetitive DNA sequence composition.

Super-Resolution Microscopy Reveals Diversity of Plant Centromere Architecture.

Centromeres are essential for proper chromosome segregation to the daughter cells during mitosis and meiosis. Chromosomes of most eukaryotes studied so far have regional centromeres that form primary constrictions on metaphase chromosomes. These monocentric chromosomes vary from point centromeres to so-called "meta-polycentromeres", with multiple centromere domains in an extended primary constriction, as identified in Pisum and Lathyrus species. However, in various animal and plant lineages centromeres are distributed along almost the entire chromosome length. Therefore, they are called holocentromeres. In holocentric plants, centromere-specific proteins, at which spindle fibers usually attach, are arranged contiguously (line-like), in clusters along the chromosomes or in bands. Here, we summarize findings of ultrastructural investigations using immunolabeling with centromere-specific antibodies and super-resolution microscopy to demonstrate the structural diversity of plant centromeres. A classification of the different centromere types has been suggested based on the distribution of spindle attachment sites. Based on these findings we discuss the possible evolution and advantages of holocentricity, and potential strategies to segregate holocentric chromosomes correctly.

TAREAN: a computational tool for identification and characterization of satellite DNA from unassembled short reads.

Satellite DNA is one of the major classes of repetitive DNA, characterized by tandemly arranged repeat copies that form contiguous arrays up to megabases in length. This type of genomic organization makes satellite DNA difficult to assemble, which hampers characterization of satellite sequences by computational analysis of genomic contigs. Here, we present tandem repeat analyzer (TAREAN), a novel computational pipeline that circumvents this problem by detecting satellite repeats directly from unassembled short reads. The pipeline first employs graph-based sequence clustering to identify groups of reads that represent repetitive elements. Putative satellite repeats are subsequently detected by the presence of circular structures in their cluster graphs. Consensus sequences of repeat monomers are then reconstructed from the most frequent k-mers obtained by decomposing read sequences from corresponding clusters. The pipeline performance was successfully validated by analyzing low-pass genome sequencing data from five plant species where satellite DNA was previously experimentally characterized. Moreover, novel satellite repeats were predicted for the genome of Vicia faba and three of these repeats were verified by detecting their sequences on metaphase chromosomes using fluorescence in situ hybridization.

Holocentromeres in Rhynchospora are associated with genome-wide centromere-specific repeat arrays interspersed among euchromatin.

Holocentric chromosomes lack a primary constriction, in contrast to monocentrics. They form kinetochores distributed along almost the entire poleward surface of the chromatids, to which spindle fibers attach. No centromere-specific DNA sequence has been found for any holocentric organism studied so far. It was proposed that centromeric repeats, typical for many monocentric species, could not occur in holocentrics, most likely because of differences in the centromere organization. Here we show that the holokinetic centromeres of the Cyperaceae Rhynchospora pubera are highly enriched by a centromeric histone H3 variant-interacting centromere-specific satellite family designated "Tyba" and by centromeric retrotransposons (i.e., CRRh) occurring as genome-wide interspersed arrays. Centromeric arrays vary in length from 3 to 16 kb and are intermingled with gene-coding sequences and transposable elements. We show that holocentromeres of metaphase chromosomes are composed of multiple centromeric units rather than possessing a diffuse organization, thus favoring the polycentric model. A cell-cycle-dependent shuffling of multiple centromeric units results in the formation of functional (poly)centromeres during mitosis. The genome-wide distribution of centromeric repeat arrays interspersing the euchromatin provides a previously unidentified type of centromeric chromatin organization among eukaryotes. Thus, different types of holocentromeres exist in different species, namely with and without centromeric repetitive sequences.

Karyotypes and Distribution of Tandem Repeat Sequences in Brassica nigra Determined by Fluorescence in situ Hybridization.

Whole-genome shotgun reads were analyzed to determine the repeat sequence composition in the genome of black mustard, Brassica nigra (L.) Koch. The analysis showed that satellite DNA sequences are very abundant in the black mustard genome. The distribution pattern of 7 new tandem repeats (BnSAT13, BnSAT28, BnSAT68, BnSAT76, BnSAT114, BnSAT180, and BnSAT200) on black mustard chromosomes was visualized using fluorescence in situ hybridization (FISH). The FISH signals of BnSAT13 and BnSAT76 provided useful cytogenetic markers; their position and fluorescence intensity allowed for unambiguous identification of all 8 somatic metaphase chromosomes. A karyotype showing the location and fluorescence intensity of these tandem repeat sequences together with the position of rDNAs and centromeric retrotransposons of Brassica (CRB) was constructed. The establishment of the FISH-based karyotype in B. nigra provides valuable information that can be used in detailed analyses of B. nigra accessions and derived allopolyploid Brassica species containing the B genome.

In vitro analysis of circumferential joint replacement, including bilateral facet joint replacement with lateral lumber disc prosthesis: a parametric investigation of disc sizing.

PURPOSE: Lateral lumbar disc prosthesis (LLDP) is an innovative device used to restore motion in select patients through a lateral retroperitoneal approach. No in vitro biomechanical studies have been published. Further, the potential for in toto circumferential joint restoration when use of this anterior disc is combined with facet replacement remains unqualified but signifies a potentially interesting clinical direction. METHODS: Researchers conducted a biomechanical feasibility study of an LLDP designed to investigate parameters of disc sizing used with bilateral facet joint replacement in a cadaveric model. Tested constructs at L4-L5 included (1) intact, (2) LLDP, (3) LLDP + wide discectomy, (4) LLDP + bilateral facetectomy, and (5) LLDP + bilateral facet joint replacement (BFJR). Investigators tested instrumented constructs (2-5) with an LLDP at compact-fit and lax-fit heights and used raw data to perform statistical analysis by repeated measures analysis of variance (ANOVA), along with Student-Newman-Keuls post hoc analysis (p ≤ 0.05). RESULTS: Increased height of the LLDP resulted in significantly less motion compared with intact. Widening the discectomy while using lax-fit sizing led to motion similar to intact in flexion-extension. As expected, motion was greater with lax-fit height than with compact-fit height in all loading modes and constructs, as is noted with a widened discectomy. The L4-L5 center of rotation was maintained regardless of placement of the LLDP. CONCLUSIONS: After bilateral facetectomy, reconstruction of the three-joint complex achieved by combining the LLDP with BFJR may provide a viable alternative to current clinical treatment regimens.

Centromere and telomere sequence alterations reflect the rapid genome evolution within the carnivorous plant genus Genlisea.

Linear chromosomes of eukaryotic organisms invariably possess centromeres and telomeres to ensure proper chromosome segregation during nuclear divisions and to protect the chromosome ends from deterioration and fusion, respectively. While centromeric sequences may differ between species, with arrays of tandemly repeated sequences and retrotransposons being the most abundant sequence types in plant centromeres, telomeric sequences are usually highly conserved among plants and other organisms. The genome size of the carnivorous genus Genlisea (Lentibulariaceae) is highly variable. Here we study evolutionary sequence plasticity of these chromosomal domains at an intrageneric level. We show that Genlisea nigrocaulis (1C = 86 Mbp; 2n = 40) and G. hispidula (1C = 1550 Mbp; 2n = 40) differ as to their DNA composition at centromeres and telomeres. G. nigrocaulis and its close relative G. pygmaea revealed mainly 161 bp tandem repeats, while G. hispidula and its close relative G. subglabra displayed a combination of four retroelements at centromeric positions. G. nigrocaulis and G. pygmaea chromosome ends are characterized by the Arabidopsis-type telomeric repeats (TTTAGGG); G. hispidula and G. subglabra instead revealed two intermingled sequence variants (TTCAGG and TTTCAGG). These differences in centromeric and, surprisingly, also in telomeric DNA sequences, uncovered between groups with on average a > 9-fold genome size difference, emphasize the fast genome evolution within this genus. Such intrageneric evolutionary alteration of telomeric repeats with cytosine in the guanine-rich strand, not yet known for plants, might impact the epigenetic telomere chromatin modification.

Centromeres Off the Hook: Massive Changes in Centromere Size and Structure Following Duplication of CenH3 Gene in Fabeae Species.

In most eukaryotes, centromere is determined by the presence of the centromere-specific histone variant CenH3. Two types of chromosome morphology are generally recognized with respect to centromere organization. Monocentric chromosomes possess a single CenH3-containing domain in primary constriction, whereas holocentric chromosomes lack the primary constriction and display dispersed distribution of CenH3. Recently, metapolycentric chromosomes have been reported in Pisum sativum, representing an intermediate type of centromere organization characterized by multiple CenH3-containing domains distributed across large parts of chromosomes that still form a single constriction. In this work, we show that this type of centromere is also found in other Pisum and closely related Lathyrus species, whereas Vicia and Lens genera, which belong to the same legume tribe Fabeae, possess only monocentric chromosomes. We observed extensive variability in the size of primary constriction and the arrangement of CenH3 domains both between and within individual Pisum and Lathyrus species, with no obvious correlation to genome or chromosome size. Search for CenH3 gene sequences revealed two paralogous variants, CenH3-1 and CenH3-2, which originated from a duplication event in the common ancestor of Fabeae species. The CenH3-1 gene was subsequently lost or silenced in the lineage leading to Vicia and Lens, whereas both genes are retained in Pisum and Lathyrus. Both of these genes appear to have evolved under purifying selection and produce functional CenH3 proteins which are fully colocalized. The findings described here provide the first evidence for a highly dynamic centromere structure within a group of closely related species, challenging previous concepts of centromere evolution.

Analysis of the giant genomes of Fritillaria (Liliaceae) indicates that a lack of DNA removal characterizes extreme expansions in genome size.

Plants exhibit an extraordinary range of genome sizes, varying by > 2000-fold between the smallest and largest recorded values. In the absence of polyploidy, changes in the amount of repetitive DNA (transposable elements and tandem repeats) are primarily responsible for genome size differences between species. However, there is ongoing debate regarding the relative importance of amplification of repetitive DNA versus its deletion in governing genome size. Using data from 454 sequencing, we analysed the most repetitive fraction of some of the largest known genomes for diploid plant species, from members of Fritillaria. We revealed that genomic expansion has not resulted from the recent massive amplification of just a handful of repeat families, as shown in species with smaller genomes. Instead, the bulk of these immense genomes is composed of highly heterogeneous, relatively low-abundance repeat-derived DNA, supporting a scenario where amplified repeats continually accumulate due to infrequent DNA removal. Our results indicate that a lack of deletion and low turnover of repetitive DNA are major contributors to the evolution of extremely large genomes and show that their size cannot simply be accounted for by the activity of a small number of high-abundance repeat families.

Repeatless and repeat-based centromeres in potato: implications for centromere evolution.

Centromeres in most higher eukaryotes are composed of long arrays of satellite repeats. By contrast, most newly formed centromeres (neocentromeres) do not contain satellite repeats and instead include DNA sequences representative of the genome. An unknown question in centromere evolution is how satellite repeat-based centromeres evolve from neocentromeres. We conducted a genome-wide characterization of sequences associated with CENH3 nucleosomes in potato (Solanum tuberosum). Five potato centromeres (Cen4, Cen6, Cen10, Cen11, and Cen12) consisted primarily of single- or low-copy DNA sequences. No satellite repeats were identified in these five centromeres. At least one transcribed gene was associated with CENH3 nucleosomes. Thus, these five centromeres structurally resemble neocentromeres. By contrast, six potato centromeres (Cen1, Cen2, Cen3, Cen5, Cen7, and Cen8) contained megabase-sized satellite repeat arrays that are unique to individual centromeres. The satellite repeat arrays likely span the entire functional cores of these six centromeres. At least four of the centromeric repeats were amplified from retrotransposon-related sequences and were not detected in Solanum species closely related to potato. The presence of two distinct types of centromeres, coupled with the boom-and-bust cycles of centromeric satellite repeats in Solanum species, suggests that repeat-based centromeres can rapidly evolve from neocentromeres by de novo amplification and insertion of satellite repeats in the CENH3 domains.

Stretching the rules: monocentric chromosomes with multiple centromere domains.

The centromere is a functional chromosome domain that is essential for faithful chromosome segregation during cell division and that can be reliably identified by the presence of the centromere-specific histone H3 variant CenH3. In monocentric chromosomes, the centromere is characterized by a single CenH3-containing region within a morphologically distinct primary constriction. This region usually spans up to a few Mbp composed mainly of centromere-specific satellite DNA common to all chromosomes of a given species. In holocentric chromosomes, there is no primary constriction; the centromere is composed of many CenH3 loci distributed along the entire length of a chromosome. Using correlative fluorescence light microscopy and high-resolution electron microscopy, we show that pea (Pisum sativum) chromosomes exhibit remarkably long primary constrictions that contain 3-5 explicit CenH3-containing regions, a novelty in centromere organization. In addition, we estimate that the size of the chromosome segment delimited by two outermost domains varies between 69 Mbp and 107 Mbp, several factors larger than any known centromere length. These domains are almost entirely composed of repetitive DNA sequences belonging to 13 distinct families of satellite DNA and one family of centromeric retrotransposons, all of which are unevenly distributed among pea chromosomes. We present the centromeres of Pisum as novel "meta-polycentric" functional domains. Our results demonstrate that the organization and DNA composition of functional centromere domains can be far more complex than previously thought, do not require single repetitive elements, and do not require single centromere domains in order to segregate properly. Based on these findings, we propose Pisum as a useful model for investigation of centromere architecture and the still poorly understood role of repetitive DNA in centromere evolution, determination, and function.

Selfish supernumerary chromosome reveals its origin as a mosaic of host genome and organellar sequences.

Supernumerary B chromosomes are optional additions to the basic set of A chromosomes, and occur in all eukaryotic groups. They differ from the basic complement in morphology, pairing behavior, and inheritance and are not required for normal growth and development. The current view is that B chromosomes are parasitic elements comparable to selfish DNA, like transposons. In contrast to transposons, they are autonomously inherited independent of the host genome and have their own mechanisms of mitotic or meiotic drive. Although B chromosomes were first described a century ago, little is known about their origin and molecular makeup. The widely accepted view is that they are derived from fragments of A chromosomes and/or generated in response to interspecific hybridization. Through next-generation sequencing of sorted A and B chromosomes, we show that B chromosomes of rye are rich in gene-derived sequences, allowing us to trace their origin to fragments of A chromosomes, with the largest parts corresponding to rye chromosomes 3R and 7R. Compared with A chromosomes, B chromosomes were also found to accumulate large amounts of specific repeats and insertions of organellar DNA. The origin of rye B chromosomes occurred an estimated ∼1.1-1.3 Mya, overlapping in time with the onset of the genus Secale (1.7 Mya). We propose a comprehensive model of B chromosome evolution, including its origin by recombination of several A chromosomes followed by capturing of additional A-derived and organellar sequences and amplification of B-specific repeats.