DArTseq genotyping facilitates the transfer of "exotic" chromatin from a Secale cereale × S. strictum hybrid into wheat.

Cultivated and wild species of the genus rye (Secale) are important but underexploited gene sources for increasing the genetic diversity of bread wheat. Gene transfer is possible via bridge genetic materials derived from intergeneric hybrids. During this process, it is essential to precisely identify the rye chromatin in the wheat genetic background. In the present study, backcross generation BC2F8 from a cross between Triticum aestivum (Mv9kr1) and S. cereanum ('Kriszta,' a cultivar from the artificial hybrid of S. cereale and S. strictum) was screened using in-situ hybridization (GISH and FISH) and analyzed by DArTseq genotyping in order to select potentially agronomically useful genotypes for prebreeding purposes. Of the 329,267 high-quality short sequence reads generated, 27,822 SilicoDArT and 8,842 SNP markers specific to S. cereanum 1R-7R chromosomes were identified. Heatmaps of the marker densities along the 'Lo7' rye reference pseudomolecules revealed subtle differences between the FISH- and DArTseq-based results. This study demonstrates that the "exotic" rye chromatin of S. cereanum introgressed into wheat can be reliably identified by high-throughput DArTseq genotyping. The Mv9kr1-'Kriszta' addition and translocation lines presented here may serve as valuable prebreeding genetic materials for the development of stress-tolerant or disease-resistant wheat varieties.

Stressed? Break-induced replication comes to the rescue!

Break-induced replication (BIR) is a homologous recombination (HR) pathway that repairs one-ended DNA double-strand breaks (DSBs), which can result from replication fork collapse, telomere erosion, and other events. Eukaryotic BIR has been mainly investigated in yeast, where it is initiated by invasion of the broken DNA end into a homologous sequence, followed by extensive replication synthesis proceeding to the chromosome end. Multiple recent studies have described BIR in mammalian cells, the properties of which show many similarities to yeast BIR. While HR is considered as "error-free" mechanism, BIR is highly mutagenic and frequently leads to chromosomal rearrangements-genetic instabilities known to promote human disease. In addition, it is now recognized that BIR is highly stimulated by replication stress (RS), including RS constantly present in cancer cells, implicating BIR as a contributor to cancer genesis and progression. Here, we discuss the past and current findings related to the mechanism of BIR, the association of BIR with replication stress, and the destabilizing effects of BIR on the eukaryotic genome. Finally, we consider the potential for exploiting the BIR machinery to develop anti-cancer therapeutics.

Complex genomic rearrangements of the Y chromosome in a premature infant.

BACKGROUND: Chromoanagenesis is an umbrella term used to describe catastrophic "all at once" cellular events leading to the chaotic reconstruction of chromosomes. It is characterized by numerous rearrangements involving a small number of chromosomes/loci, copy number gains in combination with deletions, reconstruction of chromosomal fragments with improper order/orientation, and preserved heterozygosity in copy number neutral regions. Chromoanagesis is frequently described in association with cancer; however, it has also been described in the germline. The clinical features associated with constitutional chromoanagenesis are typically due to copy number changes and/or disruption of genes or regulatory regions. CASE PRESENTATION: We present an 8-year-old male patient with complex rearrangements of the Y chromosome including a ring Y chromosome, a derivative Y;21 chromosome, and a complex rearranged Y chromosome. These chromosomes were characterized by G-banded chromosome analysis, SNP microarray, interphase FISH, and metaphase FISH. The mechanism(s) by which these rearrangements occurred is unclear; however, it is evocative of chromoanagenesis. CONCLUSION: This case is a novel example of suspected germline chromoanagenesis leading to large copy number changes that are well-tolerated, possibly because only the sex chromosomes are affected.

Complete genome of the Medicago anthracnose fungus, Colletotrichum destructivum, reveals a mini-chromosome-like region within a core chromosome.

Colletotrichum destructivum (Cd) is a phytopathogenic fungus causing significant economic losses on forage legume crops (Medicago and Trifolium species) worldwide. To gain insights into the genetic basis of fungal virulence and host specificity, we sequenced the genome of an isolate from Medicago sativa using long-read (PacBio) technology. The resulting genome assembly has a total length of 51.7 Mb and comprises ten core chromosomes and two accessory chromosomes, all of which were sequenced from telomere to telomere. A total of 15, 631 gene models were predicted, including genes encoding potentially pathogenicity-related proteins such as candidate-secreted effectors (484), secondary metabolism key enzymes (110) and carbohydrate-active enzymes (619). Synteny analysis revealed extensive structural rearrangements in the genome of Cd relative to the closely related Brassicaceae pathogen, Colletotrichum higginsianum. In addition, a 1.2 Mb species-specific region was detected within the largest core chromosome of Cd that has all the characteristics of fungal accessory chromosomes (transposon-rich, gene-poor, distinct codon usage), providing evidence for exchange between these two genomic compartments. This region was also unique in having undergone extensive intra-chromosomal segmental duplications. Our findings provide insights into the evolution of accessory regions and possible mechanisms for generating genetic diversity in this asexual fungal pathogen.

Major chromosome rearrangements in intergeneric wheat × rye hybrids in compatible and incompatible crosses detected by GBS read coverage analysis.

The presence of incompatibility alleles in primary amphidiploids constitutes a reproductive barrier in newly synthesized wheat-rye hybrids. To overcome this barrier, the genome stabilization process includes large-scale chromosome rearrangements. In incompatible crosses resulting in fertile amphidiploids, the elimination of one of the incompatible alleles Eml-A1 or Eml-R1b can occur already in the somatic tissue of the wheat × rye hybrid embryo. We observed that the interaction of incompatible loci Eml-A1 of wheat and Eml-R1b of rye after overcoming embryo lethality leads to hybrid sterility in primary triticale. During subsequent seed reproductions (R1, R2 or R3) most of the chromosomes of A, B, D and R subgenomes undergo rearrangement or eliminations to increase the fertility of the amphidiploid by natural selection. Genotyping-by-sequencing (GBS) coverage analysis showed that improved fertility is associated with the elimination of entire and partial chromosomes carrying factors that either cause the disruption of plant development in hybrid plants or lead to the restoration of the euploid number of chromosomes (2n = 56) in the absence of one of the incompatible alleles. Highly fertile offspring obtained in compatible and incompatible crosses can be successfully adapted for the production of triticale pre-breeding stocks.

Comparative karyotype analysis of the red brocket deer (M. americana sensu lato and M. rufa) complex: evidence of drastic chromosomal evolution and implications on speciation process.

Chromosomal rearrangements are often associated with playing a role in the speciation process. However, the underlying mechanism that favors the genetic isolation associated with chromosomal changes remains elusive. In this sense, the genus Mazama is recognized by its high level of karyotype diversity among species with similar morphology. A cryptic species complex has been identified within the genus, with the red brocket deer (Mazama americana and Mazama rufa) being the most impressive example. The chromosome variation was clustered in cytotypes with diploid numbers ranging from 42 to 53 and was correlated with geographical location. We conducted an analysis of chromosome evolution of the red brocket deer complex using comparative chromosome painting and Bacterial Artificial Chromosome (BAC) clones among different cytotypes. The aim was to deepen our understanding of the karyotypic relationships within the red brocket, thereby elucidating the significant chromosome variation among closely related species. This underscores the significance of chromosome changes as a key evolutionary process shaping their genomes. The results revealed the presence of three distinct cytogenetic lineages characterized by significant karyotypic divergence, suggesting the existence of efficient post-zygotic barriers. Tandem fusions constitute the main mechanism driving karyotype evolution, following a few centric fusions, inversion X-autosomal fusions. The BAC mapping has improved our comprehension of the karyotypic relationships within the red brocket deer complex, prompting questions regarding the role of these changes in the speciation process. We propose the red brocket as a model group to investigate how chromosomal changes contribute to isolation and explore the implications of these changes in taxonomy and conservation.

Oligo-barcode illuminates holocentric karyotype evolution in Rhynchospora (Cyperaceae).

Holocentric karyotypes are assumed to rapidly evolve through chromosome fusions and fissions due to the diffuse nature of their centromeres. Here, we took advantage of the recent availability of a chromosome-scale reference genome for Rhynchospora breviuscula, a model species of this holocentric genus, and developed the first set of oligo-based barcode probes for a holocentric plant. These probes were applied to 13 additional species of the genus, aiming to investigate the evolutionary dynamics driving the karyotype evolution in Rhynchospora. The two sets of probes were composed of 27,392 (green) and 23,968 (magenta) oligonucleotides (45-nt long), and generated 15 distinct FISH signals as a unique barcode pattern for the identification of all five chromosome pairs of the R. breviuscula karyotype. Oligo-FISH comparative analyzes revealed different types of rearrangements, such as fusions, fissions, putative inversions and translocations, as well as genomic duplications among the analyzed species. Two rounds of whole genome duplication (WGD) were demonstrated in R. pubera, but both analyzed accessions differed in the complex chain of events that gave rise to its large, structurally diploidized karyotypes with 2n = 10 or 12. Considering the phylogenetic relationships and divergence time of the species, the specificity and synteny of the probes were maintained up to species with a divergence time of ~25 My. However, karyotype divergence in more distant species hindered chromosome mapping and the inference of specific events. This barcoding system is a powerful tool to study chromosomal variations and genomic evolution in holocentric chromosomes of Rhynchospora species.

Unraveling the Unusual Subgenomic Organization in the Neopolyploid Free-Living Flatworm Macrostomum lignano.

Whole genome duplication (WGD) is an evolutionary event resulting in a redundancy of genetic material. Different mechanisms of WGD, allo- or autopolyploidization, lead to distinct evolutionary trajectories of newly formed polyploids. Genome studies on such species are important for understanding the early stages of genome evolution. However, assembling neopolyploid is a challenging task due to the presence of 2 homologous (or homeologous) chromosome sets and therefore the existence of the extended paralogous regions in its genome. Post-WGD evolution of polyploids includes cytogenetic diploidization leading to the formation of species, whose polyploid origin might be hidden by disomic inheritance. Earlier we uncovered the hidden polyploid origin of the free-living flatworms of the genus Macrostomum (Macrostomum lignano, M. janickei, and M. mirumnovem). Cytogenetic diploidization in these species is accompanied by intensive chromosomal rearrangements including chromosomes fusions. In this study, we unravel the M. lignano genome organization through generation and sequencing of 2 sublines of the commonly used inbred line of M. lignano (called DV1) differing only in a copy number of the largest chromosome (MLI1). Using nontrivial assembly free comparative analysis of their genomes, we deciphered DNA sequences belonging to MLI1 and validated them by sequencing the pool of microdissected MLI1. Here we presented the uncommon mechanism of genome rediplodization of M. lignano, which consists of (i) presence of 3 subgenomes, which emerged via formation of large fused chromosomes and its variants, and (ii) sustaining their heterozygosity through inter- and intrachromosomal rearrangements.

Comparative analysis of repetitive DNA in dysploid and non-dysploid Phaseolus beans.

Structural karyotype changes result from ectopic recombination events frequently associated with repetitive DNA. Although most Phaseolus species present relatively stable karyotypes with 2n = 22 chromosomes, the karyotypes of species of the Leptostachyus group show high rates of structural rearrangements, including a nested chromosome fusion that led to the dysploid chromosome number of the group (2n = 20). We examined the roles of repetitive landscapes in the rearrangements of species of the Leptostachyus group using genome-skimming data to characterize the repeatome in a range of Phaseolus species and compared them to species of that group (P. leptostachyus and P. macvaughii). LTR retrotransposons, especially the Ty3/gypsy lineage Chromovirus, were the most abundant elements in the genomes. Differences in the abundance of Tekay, Retand, and SIRE elements between P. macvaughii and P. leptostachyus were reflected in their total amounts of Ty3/gypsy and Ty1/copia. The satellite DNA fraction was the most divergent among the species, varying both in abundance and distribution, even between P. leptostachyus and P. macvaughii. The rapid turnover of repeats in the Leptostachyus group may be associated with the several rearrangements observed.

Fixed Allele Differences Associated With the Centromere Reveal Chromosome Morphology and Rearrangements in a Reptile (Varanus acanthurus BOULENGER).

Chromosome rearrangements are often implicated with genomic divergence and are proposed to be associated with species evolution. Rearrangements alter the genomic structure and interfere with homologous recombination by isolating a portion of the genome. Integration of multiplatform next-generation DNA sequencing technologies has enabled putative identification of chromosome rearrangements in many taxa; however, integrating these data sets with cytogenetics is still uncommon beyond model genetic organisms. Therefore, to achieve the ultimate goal for the genomic classification of eukaryotic organisms, physical chromosome mapping remains critical. The ridge-tailed goannas (Varanus acanthurus BOULENGER) are a group of dwarf monitor lizards comprised of several species found throughout northern Australia. These lizards exhibit extreme divergence at both the genic and chromosomal levels. The chromosome polymorphisms are widespread extending across much of their distribution, raising the question if these polymorphisms are homologous within the V. acanthurus complex. We used a combined genomic and cytogenetic approach to test for homology across divergent populations with morphologically similar chromosome rearrangements. We showed that more than one chromosome pair was involved with the widespread rearrangements. This finding provides evidence to support de novo chromosome rearrangements have occurred within populations. These chromosome rearrangements are characterized by fixed allele differences originating in the vicinity of the centromeric region. We then compared this region with several other assembled genomes of reptiles, chicken, and the platypus. We demonstrated that the synteny of genes in Reptilia remains conserved despite centromere repositioning across these taxa.

Case report: Optical genome mapping revealed double rearrangements in a male undergoing preimplantation genetic testing.

Chromosome rearrangement is one of the main causes of abortion. In individuals with double chromosomal rearrangements, the abortion rate and the risk of producing abnormal chromosomal embryos are increased. In our study, preimplantation genetic testing for structural rearrangement (PGT-SR) was performed for a couple because of recurrent abortion and the karyotype of the male was 45, XY der (14; 15)(q10; q10). The PGT-SR result of the embryo in this in vitro fertilization (IVF) cycle showed microduplication and microdeletion at the terminals of chromosomes 3 and 11, respectively. Therefore, we speculated whether the couple might have a cryptic reciprocal translocation which was not detected by karyotyping. Then, optical genome mapping (OGM) was performed for this couple, and cryptic balanced chromosomal rearrangements were detected in the male. The OGM data were consistent with our hypothesis according to previous PGT results. Subsequently, this result was verified by fluorescence in situ hybridization (FISH) in metaphase. In conclusion, the male's karyotype was 45, XY, t(3; 11)(q28; p15.4), der(14; 15)(q10; q10). Compared with traditional karyotyping, chromosomal microarray, CNV-seq and FISH, OGM has significant advantages in detecting cryptic and balanced chromosomal rearrangements.

Festuca pratensis-like Subgenome Reassembly from a "Chromosomal Cocktail" in the Intergeneric Festulolium (Poaceae) Hybrid: A Rare Chromoanagenesis Event in Grasses.

Festuca and Lolium grass species are used for Festulolium hybrid variety production where they display trait complementarities. However, at the genome level, they show antagonisms and a broad scale of rearrangements. A rare case of an unstable hybrid, a donor plant manifesting pronounced variability of its clonal parts, was discovered in the F2 group of 682 plants of Lolium multiflorum × Festuca arundinacea (2n = 6x = 42). Five phenotypically distinct clonal plants were determined to be diploids, having only 14 chromosomes out of the 42 in the donor. GISH defined the diploids as having the basic genome from F. pratensis (2n = 2x = 14), one of the progenitors of F. arundinacea (2n = 6x = 42), with minor components from L. multiflorum and another subgenome, F. glaucescens. The 45S rDNA position on two chromosomes also corresponded to the variant of F. pratensis in the F. arundinacea parent. In the highly unbalanced donor genome, F. pratensis was the least represented, but the most involved in numerous recombinant chromosomes. Specifically, FISH highlighted 45S rDNA-containing clusters involved in the formation of unusual chromosomal associations in the donor plant, suggesting their active role in karyotype realignment. The results of this study show that F. pratensis chromosomes have a particular fundamental drive for restructuring, which prompts the disassembly/reassembly processes. The finding of F. pratensis "escaping" and rebuilding itself from the chaotic "chromosomal cocktail" of the donor plant points to a rare chromoanagenesis event and extends the view of plant genome plasticity.

Fusion genes aiding the diagnosis of soft tissue tumours of the oral cavity: From bench to bedside.

Soft tissue tumours (STT) are a heterogeneous group of benign, malignant, and intermediate/borderline mesenchymal tumours. In the oral and maxillofacial region, less than 3% of all lesions correspond to benign STT and <1% are sarcomas. Overlapping microscopic features may lead to quite challenging diagnostic processes. Translocations and fusion genes are frequent, and type-specific genetic alterations are detected in these tumours. The detection of such alterations by classic cytogenetic, FISH, RT-PCR or NGS can help to define the diagnosis. This narrative review aims to review fusion genes reported for STT that affect the oral cavity and their use in diagnostic molecular pathology. Basic concepts regarding mechanisms of fusion genes formation are presented to clarify this information for surgical pathologists. The chromosomal rearrangements and fusion genes of adipocytic, fibroblastic and myofibroblastic, vascular, pericytic, smooth muscle, skeletal muscle, chondro-osseous, and uncertain origin STT are summarised. The advance in molecular pathology techniques has led not only to a better understanding of the molecular pathogenesis of STT, but also to the development of helpful diagnostic tools. Therefore, it is important for the oral and head and neck pathologists to familiarise with the signature rearrangements and fusion genes for each tumour.

Insight into the Organization of the B10v3 Cucumber Genome by Integration of Biological and Bioinformatic Data.

The availability of a well-organized and annotated reference genome is essential for genome research and the analysis of re-sequencing approaches. The B10v3 cucumber (Cucumis sativus L.) reference genome has been sequenced and assembled into 8035 contigs, a small fraction of which have been mapped to individual chromosomes. Currently, bioinformatics methods based on comparative homology have made it possible to re-order the sequenced contigs by mapping them to the reference genomes. The B10v3 genome (North-European, Borszczagowski line) was rearranged against the genomes of cucumber 9930 ('Chinese Long' line) and Gy14 (North American line). Furthermore, a better insight into the organization of the B10v3 genome was obtained by integrating the data available in the literature on the assignment of contigs to chromosomes in the B10v3 genome with the results of the bioinformatic analysis. The combination of information on the markers used in the assembly of the B10v3 genome and the results of FISH and DArT-seq experiments confirmed the reliability of the in silico assignment. Approximately 98% of the protein-coding genes within the chromosomes were assigned and a significant proportion of the repetitive fragments in the sequenced B10v3 genome were identified using the RagTag programme. In addition, BLAST analyses provided comparative information between the B10v3 genome and the 9930 and Gy14 data sets. This revealed both similarities and differences in the functional proteins found between the coding sequences region in the genomes. This study contributes to better knowledge and understanding of cucumber genome line B10v3.

Chromosome Fissions and Fusions Act as Barriers to Gene Flow between Brenthis Fritillary Butterflies.

Chromosome rearrangements are thought to promote reproductive isolation between incipient species. However, it is unclear how often, and under what conditions, fission and fusion rearrangements act as barriers to gene flow. Here we investigate speciation between two largely sympatric fritillary butterflies, Brenthis daphne and Brenthis ino. We use a composite likelihood approach to infer the demographic history of these species from whole-genome sequence data. We then compare chromosome-level genome assemblies of individuals from each species and identify a total of nine chromosome fissions and fusions. Finally, we fit a demographic model where effective population sizes and effective migration rate vary across the genome, allowing us to quantify the effects of chromosome rearrangements on reproductive isolation. We show that chromosomes involved in rearrangements experienced less effective migration since the onset of species divergence and that genomic regions near rearrangement points have a further reduction in effective migration rate. Our results suggest that the evolution of multiple rearrangements in the B. daphne and B. ino populations, including alternative fusions of the same chromosomes, have resulted in a reduction in gene flow. Although fission and fusion of chromosomes are unlikely to be the only processes that have led to speciation between these butterflies, this study shows that these rearrangements can directly promote reproductive isolation and may be involved in speciation when karyotypes evolve quickly.

Double-strand break repair and mis-repair in 3D.

DNA double-strand breaks (DSBs) are lesions that arise frequently from exposure to damaging agents as well as from ongoing physiological DNA transactions. Mis-repair of DSBs leads to rearrangements and structural variations in chromosomes, including insertions, deletions, and translocations implicated in disease. The DNA damage response (DDR) limits pathologic mutations and large-scale chromosome rearrangements. DSB repair initiates in 2D at DNA lesions with the stepwise recruitment of repair proteins and local chromatin remodeling which facilitates break accessibility. More complex structures are then formed via protein assembly into nanodomains and via genome folding into chromatin loops. Subsequently, 3D reorganization of DSBs is guided by clustering forces which drive the assembly of repair domains harboring multiple lesions. These domains are further stabilized and insulated into condensates via liquid-liquid phase-separation. Here, we discuss the benefits and risks associated with this 3D reorganization of the broken genome.

Mechanisms of Rapid Karyotype Evolution in Mammals.

Chromosome reshuffling events are often a foundational mechanism by which speciation can occur, giving rise to highly derivative karyotypes even amongst closely related species. Yet, the features that distinguish lineages prone to such rapid chromosome evolution from those that maintain stable karyotypes across evolutionary time are still to be defined. In this review, we summarize lineages prone to rapid karyotypic evolution in the context of Simpson's rates of evolution-tachytelic, horotelic, and bradytelic-and outline the mechanisms proposed to contribute to chromosome rearrangements, their fixation, and their potential impact on speciation events. Furthermore, we discuss relevant genomic features that underpin chromosome variation, including patterns of fusions/fissions, centromere positioning, and epigenetic marks such as DNA methylation. Finally, in the era of telomere-to-telomere genomics, we discuss the value of gapless genome resources to the future of research focused on the plasticity of highly rearranged karyotypes.

Oat chromosome and genome evolution defined by widespread terminal intergenomic translocations in polyploids.

Structural chromosome rearrangements involving translocations, fusions and fissions lead to evolutionary variation between species and potentially reproductive isolation and variation in gene expression. While the wheats (Triticeae, Poaceae) and oats (Aveneae) all maintain a basic chromosome number of x=7, genomes of oats show frequent intergenomic translocations, in contrast to wheats where these translocations are relatively rare. We aimed to show genome structural diversity and genome relationships in tetraploid, hexaploid and octoploid Avena species and amphiploids, establishing patterns of intergenomic translocations across different oat taxa using fluorescence in situ hybridization (FISH) with four well-characterized repetitive DNA sequences: pAs120, AF226603, Ast-R171 and Ast-T116. In A. agadiriana (2n=4x=28), the selected probes hybridized to all chromosomes indicating that this species originated from one (autotetraploid) or closely related ancestors with the same genomes. Hexaploid amphiploids were confirmed as having the genomic composition AACCDD, while octoploid amphiploids showed three different genome compositions: AACCCCDD, AAAACCDD or AABBCCDD. The A, B, C, and D genomes of oats differ significantly in their involvement in non-centromeric, intercalary translocations. There was a predominance of distal intergenomic translocations from the C- into the D-genome chromosomes. Translocations from A- to C-, or D- to C-genome chromosomes were less frequent, proving that at least some of the translocations in oat polyploids are non-reciprocal. Rare translocations from A- to D-, D- to A- and C- to B-genome chromosomes were also visualized. The fundamental research has implications for exploiting genomic biodiversity in oat breeding through introgression from wild species potentially with contrasting chromosomal structures and hence deleterious segmental duplications or large deletions in amphiploid parental lines.

Chromosome Evolution of the Liolaemus monticola (Liolaemidae) Complex: Chromosomal and Molecular Aspects.

Chromosomal rearrangements can directly influence population differentiation and speciation. The Liolaemus monticola complex in Chile is a unique model consisting of several chromosome races arranged in a latitudinal sequence of increasing karyotype complexity from south to north. Here, we compared chromosomal and mitochondrial cytochrome b data from 15 localities across the northern geographic distribution of L. monticola. We expanded the distribution of the previously described Multiple Fissions race (re-described as MF2), in the Coastal range between the Aconcagua River and the Petorca River, and described a new Multiple Fissions 1 (MF1) race in the Andean range. Both races present centric fissions in pairs 1 and 2, as well as a pericentric inversion in one fission product of pair 2 that changes the NOR position. Additionally, we detected a new chromosomal race north of the Petorca River, the Northern Modified 2 (NM2) race, which is polymorphic for novel centric fissions in pairs 3 and 4. Our results increase the number of chromosomal races in L. monticola to seven, suggesting a complex evolutionary history of chromosomal rearrangements, population isolation by barriers, and hybridization. These results show the relevant role of chromosome mutations in evolution, especially for highly speciose groups such as Liolaemus lizards.

Mega-sized pericentromeric blocks of simple telomeric repeats and their variants reveal patterns of chromosome evolution in ancient Cycadales genomes.

Simple telomeric repeats composed of six to seven iterating nucleotide units are important sequences typically found at the ends of chromosomes. Here we analyzed their abundance and homogeneity in 42 gymnosperm (29 newly sequenced), 29 angiosperm (one newly sequenced), and eight bryophytes using bioinformatics, conventional cytogenetic and molecular biology approaches to explore their diversity across land plants. We found more than 10 000-fold variation in the amounts of telomeric repeats among the investigated taxa. Repeat abundance was positively correlated with increasing intragenomic sequence heterogeneity and occurrence at non-telomeric positions, but there was no correlation with genome size. The highest abundance/heterogeneity was found in the gymnosperm genus Cycas (Cycadaceae), in which megabase-sized blocks of telomeric repeats (i.e., billions of copies) were identified. Fluorescent in situ hybridization experiments using variant-specific probes revealed canonical Arabidopsis-type telomeric TTTAGGG repeats at chromosome ends, while pericentromeric blocks comprised at least four major telomeric variants with decreasing abundance: TTTAGGG>TTCAGGG >TTTAAGG>TTCAAGG. Such a diversity of repeats was not found in the sister cycad family Zamiaceae or in any other species analyzed. Using immunocytochemistry, we showed that the pericentromeric blocks of telomeric repeats overlapped with histone H3 serine 10 phosphorylation signals. We show that species of Cycas have amplified their telomeric repeats in centromeric and telomeric positions on telocentric chromosomes to extraordinary high levels. The ancestral chromosome number reconstruction suggests their occurrence is unlikely to be the product of ancient Robertsonian chromosome fusions. We speculate as to how the observed chromosome dynamics may be associated with the diversification of cycads.