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.

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.

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.

Boom-Bust Turnovers of Megabase-Sized Centromeric DNA in Solanum Species: Rapid Evolution of DNA Sequences Associated with Centromeres.

Centromeres are composed of long arrays of satellite repeats in most multicellular eukaryotes investigated to date. The satellite repeat-based centromeres are believed to have evolved from "neocentromeres" that originally contained only single- or low-copy sequences. However, the emergence and evolution of the satellite repeats in centromeres has been elusive. Potato (Solanum tuberosum) provides a model system for studying centromere evolution because each of its 12 centromeres contains distinct DNA sequences, allowing comparative analysis of homoeologous centromeres from related species. We conducted genome-wide analysis of the centromeric sequences in Solanum verrucosum, a wild species closely related to potato. Unambiguous homoeologous centromeric sequences were detected in only a single centromere (Cen9) between the two species. Four centromeres (Cen2, Cen4, Cen7, and Cen10) in S. verrucosum contained distinct satellite repeats that were amplified from retrotransposon-related sequences. Strikingly, the same four centromeres in potato contain either different satellite repeats (Cen2 and Cen7) or exclusively single- and low-copy sequences (Cen4 and Cen10). Our sequence comparison of five homoeologous centromeres in two Solanum species reveals rapid divergence of centromeric sequences among closely related species. We propose that centromeric satellite repeats undergo boom-bust cycles before a favorable repeat is fixed in the population.

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.

Repeat Composition of CenH3-chromatin and H3K9me2-marked heterochromatin in Sugar Beet (Beta vulgaris).

BACKGROUND: Sugar beet (Beta vulgaris) is an important crop of temperate climate zones, which provides nearly 30 % of the world's annual sugar needs. From the total genome size of 758 Mb, only 567 Mb were incorporated in the recently published genome sequence, due to the fact that regions with high repetitive DNA contents (e.g. satellite DNAs) are only partially included. Therefore, to fill these gaps and to gain information about the repeat composition of centromeres and heterochromatic regions, we performed chromatin immunoprecipitation followed by sequencing (ChIP-Seq) using antibodies against the centromere-specific histone H3 variant of sugar beet (CenH3) and the heterochromatic mark of dimethylated lysine 9 of histone H3 (H3K9me2). RESULTS: ChIP-Seq analysis revealed that active centromeres containing CenH3 consist of the satellite pBV and the Ty3-gypsy retrotransposon Beetle7, while heterochromatin marked by H3K9me2 exhibits heterogeneity in repeat composition. H3K9me2 was mainly associated with the satellite family pEV, the Ty1-copia retrotransposon family Cotzilla and the DNA transposon superfamily of the En/Spm type. In members of the section Beta within the genus Beta, immunostaining using the CenH3 antibody was successful, indicating that orthologous CenH3 proteins are present in closely related species within this section. CONCLUSIONS: The identification of repetitive genome portions by ChIP-Seq experiments complemented the sugar beet reference sequence by providing insights into the repeat composition of poorly characterized CenH3-chromatin and H3K9me2-heterochromatin. Therefore, our work provides the basis for future research and application concerning the sugar beet centromere and repeat-rich heterochromatic regions characterized by the presence of H3K9me2.

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.

Diverse retrotransposon families and an AT-rich satellite DNA revealed in giant genomes of Fritillaria lilies.

BACKGROUND AND AIMS: The genus Fritillaria (Liliaceae) comprises species with extremely large genomes (1C = 30 000-127 000 Mb) and a bicontinental distribution. Most North American species (subgenus Liliorhiza) differ from Eurasian Fritillaria species by their distinct phylogenetic position and increased amounts of heterochromatin. This study examined the contribution of major repetitive elements to the genome obesity found in Fritillaria and identified repeats contributing to the heterochromatin arrays in Liliorhiza species. METHODS: Two Fritillaria species of similar genome size were selected for detailed analysis, one from each phylogeographical clade: F. affinis (1C = 45·6 pg, North America) and F. imperialis (1C = 43·0 pg, Eurasia). Fosmid libraries were constructed from their genomic DNAs and used for identification, sequence characterization, quantification and chromosome localization of clones containing highly repeated sequences. KEY RESULTS AND CONCLUSIONS: Repeats corresponding to 6·7 and 4·7 % of the F. affinis and F. imperialis genome, respectively, were identified. Chromoviruses and the Tat lineage of Ty3/gypsy group long terminal repeat retrotransposons were identified as the predominant components of the highly repeated fractions in the F. affinis and F. imperialis genomes, respectively. In addition, a heterogeneous, extremely AT-rich satellite repeat was isolated from F. affinis. The FriSAT1 repeat localized in heterochromatic bands makes up approx. 26 % of the F. affinis genome and substantial genomic fractions in several other Liliorhiza species. However, no evidence of a relationship between heterochromatin content and genome size variation was observed. Also, this study was unable to reveal any predominant repeats which tracked the increasing/decreasing trends of genome size evolution in Fritillaria. Instead, the giant Fritillaria genomes seem to be composed of many diversified families of transposable elements. We hypothesize that the genome obesity may be partly determined by the failure of removal mechanisms to counterbalance effectively the retrotransposon amplification.

Complex sequence organization of heterochromatin in the holocentric plant Cuscuta europaea elucidated by the computational analysis of nanopore reads.

Repeat-rich regions of higher plant genomes are usually associated with constitutive heterochromatin, a specific type of chromatin that forms tightly packed nuclear chromocenters and chromosome bands. There is a large body of cytogenetic evidence that these chromosome regions are often composed of tandemly organized satellite DNA. However, comparatively little is known about the sequence arrangement within heterochromatic regions, which are difficult to assemble due to their repeated nature. Here, we explore long-range sequence organization of heterochromatin regions containing the major satellite repeat CUS-TR24 in the holocentric plant Cuscuta europaea. Using a combination of ultra-long read sequencing with assembly-free sequence analysis, we reveal the complex structure of these loci, which are composed of short arrays of CUS-TR24 interrupted frequently by emerging simple sequence repeats and targeted insertions of a specific lineage of LINE retrotransposons. These data suggest that the organization of satellite repeats constituting heterochromatic chromosome bands can be more complex than previously envisioned, and demonstrate that heterochromatin organization can be efficiently investigated without the need for genome assembly.

Repeat-sequence turnover shifts fundamentally in species with large genomes.

Given the 2,400-fold range of genome sizes (0.06-148.9 Gbp (gigabase pair)) of seed plants (angiosperms and gymnosperms) with a broadly similar gene content (amounting to approximately 0.03 Gbp), the repeat-sequence content of the genome might be expected to increase with genome size, resulting in the largest genomes consisting almost entirely of repetitive sequences. Here we test this prediction, using the same bioinformatic approach for 101 species to ensure consistency in what constitutes a repeat. We reveal a fundamental change in repeat turnover in genomes above around 10 Gbp, such that species with the largest genomes are only about 55% repetitive. Given that genome size influences many plant traits, habits and life strategies, this fundamental shift in repeat dynamics is likely to affect the evolutionary trajectory of species lineages.

Survey of repetitive sequences in Silene latifolia with respect to their distribution on sex chromosomes.

We carried out a global survey of all major types of transposable elements in Silene latifolia, a model species with sex chromosomes that are in the early stages of their evolution. A shotgun genomic library was screened with genomic DNA to isolate and characterize the most abundant elements. We found that the most common types of elements were the subtelomeric tandem repeat X-43.1 and Gypsy retrotransposons, followed by Copia retrotransposons and LINE non-LTR elements. SINE elements and DNA transposons were less abundant. We also amplified transposable elements with degenerate primers and used them to screen the library. The localization of elements by FISH revealed that most of the Copia elements were accumulated on the Y chromosome. Surprisingly, one type of Gypsy element, which was similar to Ogre elements known from legumes, was almost absent on the Y chromosome but otherwise uniformly distributed on all chromosomes. Other types of elements were ubiquitous on all chromosomes. Moreover, we isolated and characterized two new tandem repeats. One of them, STAR-C, was localized at the centromeres of all chromosomes except the Y chromosome, where it was present on the p-arm. Its variant, STAR-Y, carrying a small deletion, was specifically localized on the q-arm of the Y chromosome. The second tandem repeat, TR1, co-localized with the 45S rDNA cluster in the subtelomeres of five pairs of autosomes. FISH analysis of other Silene species revealed that some elements (e.g., Ogre-like elements) are confined to the section Elisanthe while others (e.g. Copia or Athila-like elements) are present also in more distant species. Similarly, the centromeric satellite STAR-C was conserved in the genus Silene whereas the subtelomeric satellite X-43.1 was specific for Elisanthe section. Altogether, our data provide an overview of the repetitive sequences in Silene latifolia and revealed that genomic distribution and evolutionary dynamics differ among various repetitive elements. The unique pattern of repeat distribution is found on the Y chromosome, where some elements are accumulated while other elements are conspicuously absent, which probably reflects different forces shaping the Y chromosome.

Mitotic Spindle Attachment to the Holocentric Chromosomes of Cuscuta europaea Does Not Correlate With the Distribution of CENH3 Chromatin.

The centromere is the region on a chromosome where the kinetochore assembles and spindle microtubules attach during mitosis and meiosis. In the vast majority of eukaryotes, the centromere position is determined epigenetically by the presence of the centromere-specific histone H3 variant CENH3. In species with monocentric chromosomes, CENH3 is confined to a single chromosomal region corresponding to the primary constriction on metaphase chromosomes. By contrast, in holocentrics, CENH3 (and thus centromere activity) is distributed along the entire chromosome length. Here, we report a unique pattern of CENH3 distribution in the holocentric plant Cuscuta europaea. This species expressed two major variants of CENH3, both of which were deposited into one to three discrete regions per chromosome, whereas the rest of the chromatin appeared to be devoid of CENH3. The two CENH3 variants fully co-localized, and their immunodetection signals overlapped with the positions of DAPI-positive heterochromatic bands containing the highly amplified satellite repeat CUS-TR24. This CENH3 distribution pattern contrasted with the distribution of the mitotic spindle microtubules, which attached at uniform density along the entire chromosome length. This distribution of spindle attachment sites proves the holocentric nature of C. europaea chromosomes and also suggests that, in this species, CENH3 either lost its function or acts in parallel to an additional CENH3-free mechanism of kinetochore positioning.

Survey of extrachromosomal circular DNA derived from plant satellite repeats.

BACKGROUND: Satellite repeats represent one of the most dynamic components of higher plant genomes, undergoing rapid evolutionary changes of their nucleotide sequences and abundance in a genome. However, the exact molecular mechanisms driving these changes and their eventual regulation are mostly unknown. It has been proposed that amplification and homogenization of satellite DNA could be facilitated by extrachromosomal circular DNA (eccDNA) molecules originated by recombination-based excision from satellite repeat arrays. While the models including eccDNA are attractive for their potential to explain rapid turnover of satellite DNA, the existence of satellite repeat-derived eccDNA has not yet been systematically studied in a wider range of plant genomes. RESULTS: We performed a survey of eccDNA corresponding to nine different families and three subfamilies of satellite repeats in ten species from various genera of higher plants (Arabidopsis, Oryza, Pisum, Secale, Triticum and Vicia). The repeats selected for this study differed in their monomer length, abundance, and chromosomal localization in individual species. Using two-dimensional agarose gel electrophoresis followed by Southern blotting, eccDNA molecules corresponding to all examined satellites were detected. EccDNA occurred in the form of nicked circles ranging from hundreds to over eight thousand nucleotides in size. Within this range the circular molecules occurred preferentially in discrete size intervals corresponding to multiples of monomer or higher-order repeat lengths. CONCLUSION: This work demonstrated that satellite repeat-derived eccDNA is common in plant genomes and thus it can be seriously considered as a potential intermediate in processes driving satellite repeat evolution. The observed size distribution of circular molecules suggests that they are most likely generated by molecular mechanisms based on homologous recombination requiring long stretches of sequence similarity.

Satellite DNA in Vicia faba is characterized by remarkable diversity in its sequence composition, association with centromeres, and replication timing.

Satellite DNA, a class of repetitive sequences forming long arrays of tandemly repeated units, represents substantial portions of many plant genomes yet remains poorly characterized due to various methodological obstacles. Here we show that the genome of the field bean (Vicia faba, 2n = 12), a long-established model for cytogenetic studies in plants, contains a diverse set of satellite repeats, most of which remained concealed until their present investigation. Using next-generation sequencing combined with novel bioinformatics tools, we reconstructed consensus sequences of 23 novel satellite repeats representing 0.008-2.700% of the genome and mapped their distribution on chromosomes. We found that in addition to typical satellites with monomers hundreds of nucleotides long, V. faba contains a large number of satellite repeats with unusually long monomers (687-2033 bp), which are predominantly localized in pericentromeric regions. Using chromatin immunoprecipitation with CenH3 antibody, we revealed an extraordinary diversity of centromeric satellites, consisting of seven repeats with chromosome-specific distribution. We also found that in spite of their different nucleotide sequences, all centromeric repeats are replicated during mid-S phase, while most other satellites are replicated in the first part of late S phase, followed by a single family of FokI repeats representing the latest replicating chromatin.

Significant expansion of Vicia pannonica genome size mediated by amplification of a single type of giant retroelement.

Amplification and eventual elimination of dispersed repeats, especially those of the retroelement origin, account for most of the profound size variability observed among plant genomes. In most higher plants investigated so far, differential accumulation of various families of elements contributes to these differences. Here we report the identification of giant Ty3/gypsy-like retrotransposons from the legume plant Vicia pannonica, which alone make up approximately 38% of the genome of this species. These retrotransposons have structural features of the Ogre elements previously identified in the genomes of pea and Medicago. These features include extreme size (25 kb), the presence of an extra ORF upstream of the gag-pol region, and a putative intron dividing the prot and rt coding sequences. The Ogre elements are evenly dispersed on V. pannonica chromosomes except for terminal regions containing satellite repeats, their individual copies show extraordinary sequence similarity, and at least part of them are transcriptionally active, which suggests their recent amplification. Similar elements were also detected in several other Vicia species but in most cases in significantly lower numbers. However, there was no obvious correlation of the abundance of Ogre sequences with the genome size of these species.