Reconstructing the Evolution of Brachypodium Genomes Using Comparative Chromosome Painting

Celebrating Mendel, McClintock, and Darlington: On end-to-end chromosome fusions and nested chromosome fusions

The evolution of eukaryotic genomes is accompanied by fluctuations in chromosome number, reflecting cycles of chromosome number increase (polyploidy and centric fissions) and decrease (chromosome fusions). Although all chromosome fusions result from DNA recombination between two or more nonhomologous chromosomes, several mechanisms of descending dysploidy are exploited by eukaryotes to reduce their chromosome number. Genome sequencing and comparative genomics have accelerated the identification of inter-genome chromosome collinearity and gross chromosomal rearrangements and have shown that end-to-end chromosome fusions (EEFs) and nested chromosome fusions (NCFs) may have played a more important role in the evolution of eukaryotic karyotypes than previously thought. The present review aims to summarize the limited knowledge on the origin, frequency, and evolutionary implications of EEF and NCF events in eukaryotes and especially in land plants. The interactions between nonhomologous chromosomes in interphase nuclei and chromosome (mis)pairing during meiosis are examined for their potential importance in the origin of EEFs and NCFs. The remaining open questions that need to be addressed are discussed.

Enriched root bacterial microbiome in invaded vs native ranges of the model grass allotetraploid Brachypodium hybridum

Invasive species can shift the composition of key soil microbial groups, thus creating novel soil microbial communities. To better understand the biological drivers of invasion, we studied plant-microbial interactions in species of the Brachypodium distachyon complex, a model system for functional genomic studies of temperate grasses and bioenergy crops. While Brachypodium hybridum invasion in California is in an incipient stage, threatening natural and agricultural systems, its diploid progenitor species B. distachyon is not invasive in California. We investigated the root, soil, and rhizosphere bacterial composition of Brachypodium hybridum in both its native and invaded range, and of B. distachyon in the native range. We used high-throughput, amplicon sequencing to evaluate if the bacteria associated with these plants differ, and whether biotic controls may be driving B. hybridum invasion. Bacterial community composition of B. hybridum differed based on provenance (native or invaded range) for root, rhizosphere, and bulk soils, as did the abundance of dominant bacterial taxa. Bacteroidetes, Cyanobacteria and Bacillus spp. (species) were significantly more abundant in B. hybridum roots from the invaded range, whereas Proteobacteria, Firmicutes, Erwinia and Pseudomonas were more abundant in the native range roots. Brachypodium hybridum forms novel biotic interactions with a diverse suite of rhizosphere microbes from the invaded range, which may not exert a similar influence within its native range, ostensibly contributing to B. hybridum’s invasiveness. These associated plant microbiomes could inform future management approaches for B. hybridum in its invaded range and could be key to understanding, predicting, and preventing future plant invasions.

Chromosome and Genome Diversity in the Genus Trifolium (Fabaceae)

Trifolium L. is an economically important genus that is characterized by variable karyotypes relating to its ploidy level and basic chromosome numbers. The advent of genomic resources combined with molecular cytogenetics provides an opportunity to develop our understanding of plant genomes in general. Here, we summarize the current state of knowledge on Trifolium genomes and chromosomes and review methodologies using molecular markers that have contributed to Trifolium research. We discuss possible future applications of cytogenetic methods in research on the Trifolium genome and chromosomes.

Chromosome identification in oil palm (Elaeis guineensis) using in situ hybridization with massive pools of single copy oligonucleotides and transferability across Arecaceae species

Chromosome identification is essential for linking sequence and chromosomal maps, verifying sequence assemblies, showing structural variations and tracking inheritance or recombination of chromosomes and chromosomal segments during evolution and breeding programs. Unfortunately, identification of individual chromosomes and chromosome arms has been a major challenge for some economically important crop species with a near-continuous chromosome size range and similar morphology. Here, we developed oligonucleotide-based chromosome-specific probes that enabled us to establish a reference chromosome identification system for oil palm (Elaeis guineensis Jacq., 2n = 32). Massive oligonucleotide sequence pools were anchored to individual chromosome arms using dual and triple fluorescent in situ hybridization (EgOligoFISH). Three fluorescently tagged probe libraries were developed to contain, in total 52,506 gene-rich single-copy 47-mer oligonucleotides spanning each 0.2-0.5 Mb across strategically placed chromosome regions. They generated 19 distinct FISH signals and together with rDNA probes enabled identification of all 32 E. guineensis chromosome arms. The probes were able to identify individual homoeologous chromosome regions in the related Arecaceae palm species: American oil palm (Elaeis oleifera), date palm (Phoenix dactylifera) and coconut (Cocos nucifera) showing the comparative organization and concerted evolution of genomes in the Arecaceae. The oligonucleotide probes developed here provide a valuable approach to chromosome arm identification and allow tracking chromosome transfer in hybridization and breeding programs in oil palm, as well as comparative studies within Arecaceae.

Using carrot centromeric repeats to study karyotype relationships in the genus Daucus (Apiaceae)

Background

In the course of evolution, chromosomes undergo evolutionary changes; thus, karyotypes may differ considerably among groups of organisms, even within closely related taxa. The genus Daucus seems to be a promising model for exploring the dynamics of karyotype evolution. It comprises some 40 wild species and the cultivated carrot, a crop of great economic significance. However, Daucus species are very diverse morphologically and genetically, and despite extensive research, the taxonomic and phylogenetic relationships between them have still not been fully resolved. Although several molecular cytogenetic studies have been conducted to investigate the chromosomal structure and karyotype evolution of carrot and other Daucus species, detailed karyomorphological research has been limited to carrot and only a few wild species. Therefore, to better understand the karyotype relationships within Daucus, we (1) explored the chromosomal distribution of carrot centromeric repeats (CentDc) in 34 accessions of Daucus and related species by means of fluorescence in situ hybridization (FISH) and (2) performed detailed karyomorphological analysis in 16 of them.

Results

We determined the genomic organization of CentDc in 26 accessions of Daucus (belonging to both Daucus I and II subclades) and one accession of closely related species. The CentDc repeats were present in the centromeric regions of all chromosomes of 20 accessions (representing 11 taxa). In the other Daucus taxa, the number of chromosome pairs with CentDc signals varied depending on the species, yet their centromeric localization was conserved. In addition, precise chromosome measurements performed in 16 accessions showed the inter- and intraspecific karyological relationships among them.

Conclusions

The presence of the CentDc repeats in the genomes of taxa belonging to both Daucus subclades and one outgroup species indicated the ancestral status of the repeat. The results of our study provide useful information for further evolutionary, cytotaxonomic, and phylogenetic research on the genus Daucus and may contribute to a better understanding of the dynamic evolution of centromeric satellites in plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07853-2.

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

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

Ribosomal DNA localization on Lathyrus species chromosomes by FISH

BACKGROUND

Fluorescence In Situ Hybridization (FISH) played an essential role to locate the ribosomal RNA genes on the chromosomes that offered a new tool to study the chromosome structure and evolution in plant. The 45S and 5S rRNA genes are independent and localized at one or more loci per the chromosome complement, their positions along chromosomes offer useful markers for chromosome discriminations. In the current study FISH has been performed to locate 45S and 5S rRNA genes on the chromosomes of nine Lathyrus species belong to five different sections, all have chromosome number 2n=14, Lathyrus gorgoni Parl, Lathyrus hirsutus L., Lathyrus amphicarpos L., Lathyrus odoratus L., Lathyrus sphaericus Retz, Lathyrus incospicuus L, Lathyrus paranensis Burkart, Lathyrus nissolia L., and Lathyrus articulates L.

RESULTS

The revealed loci of 45S and 5S rDNA by FISH on metaphase chromosomes of the examined species were as follow: all of the studied species have one 45S rDNA locus and one 5S rDNA locus except L. odoratus L., L. amphicarpos L. and L. sphaericus Retz L. have two loci of 5S rDNA. Three out of the nine examined species have the loci of 45S and 5S rRNA genes on the opposite arms of the same chromosome (L. nissolia L., L. amphicarpos L., and L. incospicuus L.), while L. hirsutus L. has both loci on the same chromosome arm. The other five species showed the loci of the two types of rDNA on different chromosomes.

CONCLUSION

The detected 5S and 45S rDNA loci in Lathyrus could be used as chromosomal markers to discriminate the chromosome pairs of the examined species. FISH could discriminate only one chromosome pair out of the seven pairs in three species, in L. hirsutus L., L. nissolia L. and L. incospicuus L., and two chromosome pairs in five species, in L. paranensis Burkart, L. odoratus L., L. amphicarpos L., L. gorgoni Parl. and L. articulatus L., while it could discriminate three chromosome pairs in L. sphaericus Retz. these results could contribute into the physical genome mapping of Lathyrus species and the evolution of rDNA patterns by FISH in the coming studies in future.

The Evolution of Chromosome Numbers: Mechanistic Models and Experimental Approaches

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

Progressive refinement of the karyotyping of Brachypodium genomes

Brachypodium distachyon is a weedy grass species that is firmly established as a model for the comparative and functional genomics of temperate cereals and grasses. Its simple, nuclear genome of five chromosomes contrasts it with other relatives of the genus with different, and usually higher, basic chromosome numbers and ploidy levels. This variation in karyotypic structure affords the possibility of reconstructing evolutionary pathways that have shaped the genome structure of extant species. This Tansley insight documents how key refinements in molecular cytogenetic approaches, from simple fluorescence in situ hybridization to comparative chromosome barcoding, have enabled genome structure studies and yielded valuable information about the drivers of karyotypic reorganization and evolution in the model grass genus Brachypodium.

Repetitive sequences and structural chromosome alterations promote intraspecific variations in Zea mays L. karyotype

LTR-retrotransposons, knobs and structural chromosome alterations contribute to shape the structure and organization of the Zea mays karyotype. Our initial nuclear DNA content data of Z. mays accessions revealed an intraspecific variation (2 C = 2.00 pg to 2 C = 6.10 pg), suggesting differences in their karyotypes. We aimed to compare the karyotypes of three Z. mays accessions in search of the differences and similarities among them. Karyotype divergences were demonstrated among the accessions, despite their common chromosome number (2n = 20) and ancestral origin. Cytogenomic analyses showed that repetitive sequences and structural chromosome alterations play a significant role in promoting intraspecific nuclear DNA content variation. In addition, heterozygous terminal deletion in chromosome 3 was pointed out as a cause of lower nuclear 2 C value. Besides this, translocation was also observed in the short arm of chromosome 1. Differently, higher 2 C value was associated with the more abundant distribution of LTR-retrotransposons from the family Grande in the karyotype. Moreover, heteromorphism involving the number and position of the 180-bp knob sequence was found among the accessions. Taken together, we provide insights on the pivotal role played by repetitive sequences and structural chromosome alterations in shaping the karyotype of Z. mays.

Flexible chromosome painting based on multiplex PCR of oligonucleotides and its application for comparative chromosome analyses in Cucumis

Chromosome painting is a powerful technique for chromosome and genome studies. We developed a flexible chromosome painting technique based on multiplex PCR of a synthetic oligonucleotide (oligo) library in cucumber (Cucumis sativus L., 2n = 14). Each oligo in the library was associated with a universal as well as nested specific primers for amplification, which allow the generation of different probes from the same oligo library. We were also able to generate double-stranded labelled oligos, which produced much stronger signals than single-stranded labelled oligos, by amplification using fluorophore-conjugated primer pairs. Oligos covering cucumber chromosome 1 (Chr1) and chromosome 4 (Chr4) consisting of eight segments were synthesized in one library. Different oligo probes generated from the library painted the corresponding chromosomes/segments unambiguously, especially on pachytene chromosomes. This technique was then applied to study the homoeologous relationships among cucumber, C. hystrix and C. melo chromosomes based on cross-species chromosome painting using Chr4 probes. We demonstrated that the probe was feasible to detect interspecies chromosome homoeologous relationships and chromosomal rearrangement events. Based on its advantages and great convenience, we anticipate that this flexible oligo-painting technique has great potential for the studies of the structure, organization, and evolution of chromosomes in any species with a sequenced genome.

Comparatively Barcoded Chromosomes of Brachypodium Perennials Tell the Story of Their Karyotype Structure and Evolution

The Brachypodium genus is an informative model system for studying grass karyotype organization. Previous studies of a limited number of species and reference chromosomes have not provided a comprehensive picture of the enigmatic phylogenetic relationships in the genus. Comparative chromosome barcoding, which enables the reconstruction of the evolutionary history of individual chromosomes and their segments, allowed us to infer the relationships between putative ancestral karyotypes of extinct species and extant karyotypes of current species. We used over 80 chromosome-specific BAC (bacterial artificial chromosome) clones derived from five reference chromosomes of B. distachyon as probes against the karyotypes of twelve accessions representing five diploid and polyploid Brachypodium perennials. The results showed that descending dysploidy is common in Brachypodium and occurs primarily via nested chromosome fusions. Brachypodiumdistachyon was rejected as a putative ancestor for allotetraploid perennials and B. stacei for B. mexicanum. We propose two alternative models of perennial polyploid evolution involving either the incorporation of a putative x = 5 ancestral karyotype with different descending dysploidy patterns compared to B. distachyon chromosomes or hybridization of two x = 9 ancestors followed by genome doubling and descending dysploidy. Details of the karyotype structure and evolution in several Brachypodium perennials are revealed for the first time.

The story of promiscuous crucifers: origin and genome evolution of an invasive species, Cardamine occulta (Brassicaceae), and its relatives

BACKGROUND AND AIMS

Cardamine occulta (Brassicaceae) is an octoploid weedy species (2n = 8x = 64) originated in Eastern Asia. It has been introduced to other continents including Europe and considered to be an invasive species. Despite its wide distribution, the polyploid origin of C. occulta remained unexplored. The feasibility of comparative chromosome painting (CCP) in crucifers allowed us to elucidate the origin and genome evolution in Cardamine species. We aimed to investigate the genome structure of C. occulta in comparison with its tetraploid (2n = 4x = 32, C. kokaiensis and C. scutata) and octoploid (2n = 8x = 64, C. dentipetala) relatives.

METHODS

Genomic in situ hybridization (GISH) and large-scale CCP were applied to uncover the parental genomes and chromosome composition of the investigated Cardamine species.

KEY RESULTS

All investigated species descended from a common ancestral Cardamine genome (n = 8), structurally resembling the Ancestral Crucifer Karyotype (n = 8), but differentiated by a translocation between chromosomes AK6 and AK8. Allotetraploid C. scutata originated by hybridization between two diploid species, C. parviflora and C. amara (2n = 2x = 16). By contrast, C. kokaiensis has an autotetraploid origin from a parental genome related to C. parviflora. Interestingly, octoploid C. occulta probably originated through hybridization between the tetraploids C. scutata and C. kokaiensis. The octoploid genome of C. dentipetala probably originated from C. scutata via autopolyploidization. Except for five species-specific centromere repositionings and one pericentric inversion post-dating the polyploidization events, the parental subgenomes remained stable in the tetra- and octoploids.

CONCLUSIONS

Comparative genome structure, origin and evolutionary history was reconstructed in C. occulta and related species. For the first time, whole-genome cytogenomic maps were established for octoploid plants. Post-polyploid evolution in Asian Cardamine polyploids has not been associated with descending dysploidy and intergenomic rearrangements. The combination of different parental (sub)genomes adapted to distinct habitats provides an evolutionary advantage to newly formed polyploids by occupying new ecological niches.

Fluorescence in situ hybridization in plants: recent developments and future applications

Fluorescence in situ hybridization (FISH) was developed more than 30 years ago and has been the most paradigm-changing technique in cytogenetic research. FISH has been used to answer questions related to structure, mutation, and evolution of not only individual chromosomes but also entire genomes. FISH has served as an important tool for chromosome identification in many plant species. This review intends to summarize and discuss key technical development and applications of FISH in plants since 2006. The most significant recent advance of FISH is the development and application of probes based on synthetic oligonucleotides (oligos). Oligos specific to a repetitive DNA sequence, to a specific chromosomal region, or to an entire chromosome can be computationally identified, synthesized in parallel, and fluorescently labeled. Oligo probes designed from conserved DNA sequences from one species can be used among genetically related species, allowing comparative cytogenetic mapping of these species. The advances with synthetic oligo probes will significantly expand the applications of FISH especially in non-model plant species. Recent achievements and future applications of FISH and oligo-FISH are discussed.

The large genome size variation in the Hesperis clade was shaped by the prevalent proliferation of DNA repeats and rarer genome downsizing

BACKGROUND AND AIMS

Most crucifer species (Brassicaceae) have small nuclear genomes (mean 1C-value 617 Mb). The species with the largest genomes occur within the monophyletic Hesperis clade (Mandáková et al., Plant Physiology174: 2062-2071; also known as Clade E or Lineage III). Whereas most chromosome numbers in the clade are 6 or 7, monoploid genome sizes vary 16-fold (256-4264 Mb). To get an insight into genome size evolution in the Hesperis clade (~350 species in ~48 genera), we aimed to identify, quantify and localize in situ the repeats from which these genomes are built. We analysed nuclear repeatomes in seven species, covering the phylogenetic and genome size breadth of the clade, by low-pass whole-genome sequencing.

METHODS

Genome size was estimated by flow cytometry. Genomic DNA was sequenced on an Illumina sequencer and DNA repeats were identified and quantified using RepeatExplorer; the most abundant repeats were localized on chromosomes by fluorescence in situ hybridization. To evaluate the feasibility of bacterial artificial chromosome (BAC)-based comparative chromosome painting in Hesperis-clade species, BACs of arabidopsis were used as painting probes.

KEY RESULTS

Most biennial and perennial species of the Hesperis clade possess unusually large nuclear genomes due to the proliferation of long terminal repeat retrotransposons. The prevalent genome expansion was rarely, but repeatedly, counteracted by purging of transposable elements in ephemeral and annual species.

CONCLUSIONS

The most common ancestor of the Hesperis clade has experienced genome upsizing due to transposable element amplification. Further genome size increases, dominating diversification of all Hesperis-clade tribes, contrast with the overall stability of chromosome numbers. In some subclades and species genome downsizing occurred, presumably as an adaptive transition to an annual life cycle. The amplification versus purging of transposable elements and tandem repeats impacted the chromosomal architecture of the Hesperis-clade species.

Fluorescence in situ hybridization in plants: recent developments and future applications

Fluorescence in situ hybridization (FISH) was developed more than 30 years ago and has been the most paradigm-changing technique in cytogenetic research. FISH has been used to answer questions related to structure, mutation, and evolution of not only individual chromosomes but also entire genomes. FISH has served as an important tool for chromosome identification in many plant species. This review intends to summarize and discuss key technical development and applications of FISH in plants since 2006. The most significant recent advance of FISH is the development and application of probes based on synthetic oligonucleotides (oligos). Oligos specific to a repetitive DNA sequence, to a specific chromosomal region, or to an entire chromosome can be computationally identified, synthesized in parallel, and fluorescently labeled. Oligo probes designed from conserved DNA sequences from one species can be used among genetically related species, allowing comparative cytogenetic mapping of these species. The advances with synthetic oligo probes will significantly expand the applications of FISH especially in non-model plant species. Recent achievements and future applications of FISH and oligo-FISH are discussed.

Dissecting the chromosomal composition of mutagen-induced micronuclei in Brachypodium distachyon using multicolour FISH

BACKGROUND AND AIMS

Brachypodium distachyon (Brachypodium) is a model species for temperate cereals and other economically important grasses. Its favourable cytogenetic features and advanced molecular infrastructure make it a good model for understanding the mechanisms of instability of plant genomes after mutagenic treatment. The aim of this study was to qualitatively and quantitatively assess the composition and origin of micronuclei arising from genomic fracture, and to detect possible 'hot spots' for mutagen-induced DNA breaks.

METHODS

Seeds of Brachypodium were treated with maleic hydrazide (MH) or X-rays. The structure of mutagen-induced micronuclei was analysed in root-tip meristematic cells using multicolour fluorescence in situ hybridization (mcFISH) with various repetitive (5S rDNA, 25S rDNA, telomeric, centromeric) and low-repeat [small and large pools of bacterial artificial chromosome (BAC) clones specific for chromosome Bd1] DNA sequences.

KEY RESULTS

The majority of micronuclei derive from large, acentric fragments. X-rays caused more interstitial DNA breaks than MH. Double-strand breaks rarely occurred in distal chromosome regions. Bd1 contributed to the formation of more mutagen-induced micronuclei than expected from random chromosome involvement.

CONCLUSIONS

mcFISH with chromosome-specific BAC clones offers insight into micronuclei composition, in so far as it allows their origin and formation to be determined more specifically. A reliable assay for micronuclei composition is crucial for the development of modern genotoxicity tests using plant cells. The combination of mutagenic treatments and well-developed cytomolecular resources in Brachypodium make this model species very promising for plant mutagenesis research.

Homoeolog-specific activation of genes for heat acclimation in the allopolyploid grass Brachypodium hybridum

Background

Allopolyploid plants often show wider environmental tolerances than their ancestors; this is expected to be due to the merger of multiple distinct genomes with a fixed heterozygosity. The complex homoeologous gene expression could have been evolutionarily advantageous for the adaptation of allopolyploid plants. Despite multiple previous studies reporting homoeolog-specific gene expression in allopolyploid species, there are no clear examples of homoeolog-specific function in acclimation to a long-term stress condition.

Results

We found that the allopolyploid grass Brachypodium hybridum and its ancestor Brachypodium stacei show long-term heat stress tolerance, unlike its other ancestor, Brachypodium distachyon. To understand the physiological traits of B. hybridum, we compared the transcriptome of the 3 Brachypodium species grown under normal and heat stress conditions. We found that the expression patterns of approximately 26% and approximately 38% of the homoeolog groups in B. hybridum changed toward nonadditive expression and nonancestral expression, respectively, under normal condition. Moreover, we found that B. distachyon showed similar expression patterns between normal and heat stress conditions, whereas B. hybridum and B. stacei significantly altered their transcriptome in response to heat after 3 days of stress exposure, and homoeologs that were inherited from B. stacei may have contributed to the transcriptional stress response to heat in B. hybridum. After 15 days of heat exposure, B. hybridum and B. stacei maintained transcriptional states similar to those under normal conditions. These results suggest that an earlier response to heat that was specific to homoeologs originating from B. stacei contributed to cellular homeostasis under long-term heat stress in B. hybridum.

Conclusions

Our results provide insights into different regulatory events of the homoeo-transcriptome that are associated with stress acclimation in allopolyploid plants.

Hybridization-facilitated genome merger and repeated chromosome fusion after 8 million years

The small genus Ricotia (nine species, Brassicaceae) is confined to the eastern Mediterranean. By comparative chromosome painting and a dated multi-gene chloroplast phylogeny, we reconstructed the origin and subsequent evolution of Ricotia. The ancestral Ricotia genome originated through hybridization between two older genomes with n = 7 and n = 8 chromosomes, respectively, on the Turkish mainland during the Early Miocene (c. 17.8 million years ago, Ma). Since then, the allotetraploid (n = 15) genome has been altered by two independent descending dysploidies (DD) to n = 14 in Ricotia aucheri and the Tenuifolia clade (2 spp.). By the Late Miocene (c. 10 Ma), the latter clade started to evolve in the most diverse Ricotia core clade (6 spp.), the process preceded by a DD event to n = 13. It is noteworthy that this dysploidy was mediated by a unique chromosomal rearrangement, merging together the same two chromosomes as were merged during the origin of a fusion chromosome within the paternal n = 7 genome c. 20 Ma. This shows that within a time period of c. 8 Myr genome evolution can repeat itself and that structurally very similar chromosomes may originate repeatedly from the same ancestral chromosomes by different pathways (end-to-end translocation versus nested chromosome insertion).

Chromosome identification and reconstruction of evolutionary rearrangements in Brachypodium distachyon, B. stacei and B. hybridum

BACKGROUND AND AIMS

The Brachypodium genus represents a useful model system to study grass genome organization. Palaeogenomic analyses (e.g. Murat F, Armero A, Pont C, Klopp C, Salse J. 2017. Reconstructing the genome of the most recent common ancestor of flowering plants. Nature Genetics49: 490-496) have identified polyploidization and dysploidy as the prime mechanisms driving the diversity of plant karyotypes and nested chromosome fusions (NCFs) crucial for shaping grass chromosomes. This study compares the karyotype structure and evolution in B. distachyon (genome Bd), B. stacei (genome Bs) and in their putative allotetraploid B. hybridum (genomes BdBs).

METHODS

Brachypodium chromosomes were measured and identified using multicolour fluorescence in situ hybridization (mcFISH). For higher resolution, comparative chromosome barcoding was developed using sets of low-repeat, physically mapped B. distachyon-derived bacterial artificial chromosome (BAC) clones.

KEY RESULTS

All species had rather small chromosomes, and essentially all in the Bs genome were morphometrically indistinguishable. Seven BACs combined with two rDNA-based probes provided unambiguous and reproducible chromosome discrimination. Comparative chromosome barcoding revealed NCFs that contributed to the reduction in the x = 12 chromosome number that has been suggested for the intermediate ancestral grass karyotype. Chromosome Bd3 derives from two NCFs of three ancestral chromosomes (Os2, Os8, Os10). Chromosome Bs6 shows an ancient Os8/Os10 NCF, whilst Bs4 represents Os2 only. Chromosome Bd4 originated from a descending dysploidy that involves two NCFs of Os12, Os9 and Os11. The specific distribution of BACs along Bs9 and Bs5, in both B. stacei and B. hybridum, suggests a Bs genome-specific Robertsonian rearrangement.

CONCLUSIONS

mcFISH-based karyotyping identifies all chromosomes in Brachypodium annuals. Comparative chromosome barcoding reveals rearrangements responsible for the diverse organization of Bd and Bs genomes and provides new data regarding karyotype evolution since the split of the two diploids. The fact that no chromosome rearrangements were observed in B. hybridum compared with the karyotypes of its phylogenetic ancestors suggests prolonged genome stasis after the formation of the allotetraploid.

Brachypodium: A Monocot Grass Model Genus for Plant Biology

The genus Brachypodium represents a model system that is advancing our knowledge of the biology of grasses, including small grains, in the postgenomics era. The most widely used species, Brachypodium distachyon, is a C3 plant that is distributed worldwide. B. distachyon has a small genome, short life cycle, and small stature and is amenable to genetic transformation. Due to the intensive and thoughtful development of this grass as a model organism, it is well-suited for laboratory and field experimentation. The intent of this review is to introduce this model system genus and describe some key outcomes of nearly a decade of research since the first draft genome sequence of the flagship species, B. distachyon, was completed. We discuss characteristics and features of B. distachyon and its congeners that make the genus a valuable model system for studies in ecology, evolution, genetics, and genomics in the grasses, review current hot topics in Brachypodium research, and highlight the potential for future analysis using this system in the coming years.

Analysis of Brachypodium genomes with genome-wide optical maps

Brachypodium distachyon (n = 5) is a diploid and has been widely used as a genetic model. Brachypodium stacei (n = 10) and B. hybridum (n = 15) are species that are related to B. distachyon, leading to an hypothesis that they are part of a polyploid series based on x = 5. Several lines of evidence suggest that this hypothesis is incorrect and that the genomes of the three taxa may have evolved by a more complex process. We constructed an optical whole-genome BioNano genome (BNG) map for each species and did pairwise alignment of the BNG maps. The maps showed that B. distachyon and B. stacei are both diploid, in spite of B. stacei having twice as many chromosomes as B. distachyon, and that B. hybridum is an allopolyploid formed from hybridization between B. distachyon and B. stacei. This study also demonstrated the use of BNG maps in the detection and quantification of structural variants among the genomes.

Chromosome painting and its applications in cultivated and wild rice

BACKGROUND

The chromosome-specific probe is a fundamental tool of chromosome painting and has been commonly applied in mammalian species. The technology, however, has not been widely applied in plants due to a lack of methodologies for probe development. Identification and labeling of a large number of oligonucleotides (oligos) specific to a single chromosome offers us an opportunity to establish chromosome-specific probes in plants. However, never before has whole chromosome painting been performed in rice.

RESULTS

We developed a pooled chromosome 9-specific probe in rice, which contains 25,000 oligos based on the genome sequence of a japonica rice (Oryza sativa L., AA, 2n = 2× = 24). Chromosome 9 was easily identified in both japonica and indica rice using this chromosome 9-painting probe. The probe was also successfully used to identify and characterize chromosome 9 in additional lines of O. sativa, a translocation line, two new aneuploids associated with chromosome 9 and a wild rice (Oryza eichingeri A. Peter, CC, 2n = 2× = 24).

CONCLUSION

The study reveals that a pool of oligos specific to a chromosome is a useful tool for chromosome painting in rice.

Reconstructing the origins and the biogeography of species' genomes in the highly reticulate allopolyploid-rich model grass genus Brachypodium using minimum evolution, coalescence and maximum likelihood approaches

The identification of homeologous genomes and the biogeographical analyses of highly reticulate allopolyploid-rich groups face the challenge of incorrectly inferring the genomic origins and the biogeographical patterns of the polyploids from unreliable strictly bifurcating trees. Here we reconstruct a plausible evolutionary scenario of the diverging and merging genomes inherited by the diploid and allopolyploid species and cytotypes of the model grass genus Brachypodium. We have identified the ancestral Brachypodium genomes and inferred the paleogeographical ranges for potential hybridization events that originated its allopolyploid taxa. We also constructed a comprehensive phylogeny of Brachypodium from five nuclear and plastid genes using Species Tree Minimum Evolution allele grafting and Species Network analysis. The divergence ages of the lineages were estimated from a consensus maximum clade credibility tree using fossil calibrations, whereas ages of origin of the diploid and allopolyploid species were inferred from coalescence Bayesian methods. The biogeographical events of the genomes were reconstructed using a stratified Dispersal-Extinction-Colonization model with three temporal windows. Our combined Minimum Evolution-coalescence-Bayesian approach allowed us to infer the origins and the identities of the homeologous genomes of the Brachypodium allopolyploids, matching the expected ploidy levels of the hybrids. To date, the current extant progenitor genomes (species) are only known for B. hybridum. Putative ancestral homeologous genome have been inherited by B. mexicanum, ancestral and recent genomes by B. boissieri, and only recently evolved genomes by B. retusum and the core perennial clade allopolyploids (B. phoenicoides, B. pinnatum 4x, B. rupestre 4x). We dissected the complex spatio-temporal evolution of ancestral and recent genomes and have detected successive splitting, dispersal and merging events for dysploid homeologous genomes in diverse geographical scenarios that have led to the current extant taxa. Our data support Mid-Miocene splits of the Holarctic ancestral genomes that preceded the Late Miocene origins of Brachypodium ancestors of the modern diploid species. Successive divergences of the annual B. stacei and B. distachyon diploid genomes were implied to have occurred in the Mediterranean region during the Late Miocene-Pliocene. By contrast, a profusion of splits, range expansions and different genome mergings were inferred for the perennial diploid genomes in the Mediterranean and Eurasian regions, with sporadic colonizations and further mergings in other continents during the Quaternary. A reliable biogeographical scenario was obtained for the Brachypodium genomes and allopolyploids where homeologous genomes split from their respective diploid counterpart lineages in the same ancestral areas, showing similar or distinct dispersals. By contrast, the allopolyploid taxa remained in the same ancestral ranges after hybridization and genome doubling events. Our approach should have utility in deciphering the genomic composition and the historical biogeography of other allopolyploid-rich organismal groups, which are predominant in eukaryotes.

Methods for Cytogenetic Chromosome Barcoding and Chromosome Painting in Brachypodium distachyon and Its Relative Species

Brachypodium distachyon provides a particularly appealing object for molecular cytogenetic analysis due to its compact genome and low repetitive DNA content, as well as low (x = 5) basic number of chromosomes easily identifiable on the basis of their morphometric features. Some of these features, such as genome compactness, are shared by the other members of the genus, thus making them amenable for comparative cytogenetic mapping. Cytogenetic infrastructure established for B. distachyon was initially based on fluorescence in situ hybridization with various tandemly repeated sequences as probes. The molecular cytogenetic studies advanced greatly with the development of B. distachyon large DNA insert genomic libraries. These resources coupled with the access to the fully sequenced genome of B. distachyon enabled chromosome painting in monocots for the first time. This pioneering work was subsequently extended to other Brachypodium species, allowing insight into grass karyotype evolution. In this protocol we describe the methods of making somatic and meiotic chromosome preparations, probe labeling, FISH with BAC clones, a strategy for chromosome barcoding and chromosome painting in B. distachyon, and comparative chromosome painting in the other Brachypodium species.

Centromeric DNA characterization in the model grass Brachypodium distachyon provides insights on the evolution of the genus

Brachypodium distachyon is a well-established model monocot plant, and its small and compact genome has been used as an accurate reference for the much larger and often polyploid genomes of cereals such as Avena sativa (oats), Hordeum vulgare (barley) and Triticum aestivum (wheat). Centromeres are indispensable functional units of chromosomes and they play a core role in genome polyploidization events during evolution. As the Brachypodium genus contains about 20 species that differ significantly in terms of their basic chromosome numbers, genome size, ploidy levels and life strategies, studying their centromeres may provide important insight into the structure and evolution of the genome in this interesting and important genus. In this study, we isolated the centromeric DNA of the B. distachyon reference line Bd21 and characterized its composition via the chromatin immunoprecipitation of the nucleosomes that contain the centromere-specific histone CENH3. We revealed that the centromeres of Bd21 have the features of typical multicellular eukaryotic centromeres. Strikingly, these centromeres contain relatively few centromeric satellite DNAs; in particular, the centromere of chromosome 5 (Bd5) consists of only ~40 kb. Moreover, the centromeric retrotransposons in B. distachyon (CRBds) are evolutionarily young. These transposable elements are located both within and adjacent to the CENH3 binding domains, and have similar compositions. Moreover, based on the presence of CRBds in the centromeres, the species in this study can be grouped into two distinct lineages. This may provide new evidence regarding the phylogenetic relationships within the Brachypodium genus.

Chromosome Evolution in Connection with Repetitive Sequences and Epigenetics in Plants

Chromosome evolution is a fundamental aspect of evolutionary biology. The evolution of chromosome size, structure and shape, number, and the change in DNA composition suggest the high plasticity of nuclear genomes at the chromosomal level. Repetitive DNA sequences, which represent a conspicuous fraction of every eukaryotic genome, particularly in plants, are found to be tightly linked with plant chromosome evolution. Different classes of repetitive sequences have distinct distribution patterns on the chromosomes. Mounting evidence shows that repetitive sequences may play multiple generative roles in shaping the chromosome karyotypes in plants. Furthermore, recent development in our understanding of the repetitive sequences and plant chromosome evolution has elucidated the involvement of a spectrum of epigenetic modification. In this review, we focused on the recent evidence relating to the distribution pattern of repetitive sequences in plant chromosomes and highlighted their potential relevance to chromosome evolution in plants. We also discussed the possible connections between evolution and epigenetic alterations in chromosome structure and repatterning, such as heterochromatin formation, centromere function, and epigenetic-associated transposable element inactivation.

Exploiting repetitive sequences and BAC clones in Festuca pratensis karyotyping

The Festuca genus is thought to be the most numerous genus of the Poaceae family. One of the most agronomically important forage grasses, Festuca pratensis Huds. is treated as a model plant to study the molecular mechanisms associated with tolerance to winter stresses, including frost. However, the precise mapping of the genes governing stress tolerance in this species is difficult as its karyotype remains unrecognized. Only two F. pratensis chromosomes with 35S and 5S rDNA sequences can be easily identified, but its remaining chromosomes have not been distinguished to date. Here, two libraries derived from F. pratensis nuclear DNA with various contents of repetitive DNA sequences were used as sources of molecular probes for fluorescent in situ hybridisation (FISH), a BAC library and a library representing sequences most frequently present in the F. pratensis genome. Using FISH, six groups of DNA sequences were revealed in chromosomes on the basis of their signal position, including dispersed-like sequences, chromosome painting-like sequences, centromeric-like sequences, knob-like sequences, a group without hybridization signals, and single locus-like sequences. The last group was exploited to develop cytogenetic maps of diploid and tetraploid F. pratensis, which are presented here for the first time and provide a remarkable progress in karyotype characterization.

Diversity and association of phenotypic and metabolomic traits in the close model grasses Brachypodium distachyon, B. stacei and B. hybridum

BACKGROUND AND AIMS

Morphological traits in combination with metabolite fingerprinting were used to investigate inter- and intraspecies diversity within the model annual grasses Brachypodium distachyon, Brachypodium stacei and Brachypodium hybridum .

METHODS

Phenotypic variation of 15 morphological characters and 2219 nominal mass ( m / z ) signals generated using flow infusion electrospray ionization-mass spectrometry (FIE-MS) were evaluated in individuals from a total of 174 wild populations and six inbred lines, and 12 lines, of the three species, respectively. Basic statistics and multivariate principal component analysis and discriminant analysis were used to differentiate inter- and intraspecific variability of the two types of variable, and their association was assayed with the rcorr function.

KEY RESULTS

Basic statistics and analysis of variance detected eight phenotypic characters [(stomata) leaf guard cell length, pollen grain length, (plant) height, second leaf width, inflorescence length, number of spikelets per inflorescence, lemma length, awn length] and 434 tentatively annotated metabolite signals that significantly discriminated the three species. Three phenotypic traits (pollen grain length, spikelet length, number of flowers per inflorescence) might be genetically fixed. The three species showed different metabolomic profiles. Discriminant analysis significantly discriminated the three taxa with both morphometric and metabolome traits and the intraspecific phenotypic diversity within B. distachyon and B. stacei . The populations of B. hybridum were considerably less differentiated.

CONCLUSIONS

Highly explanatory metabolite signals together with morphological characters revealed concordant patterns of differentiation of the three taxa. Intraspecific phenotypic diversity was observed between northern and southern Iberian populations of B. distachyon and between eastern Mediterranean/south-western Asian and western Mediterranean populations of B. stacei . Significant association was found for pollen grain length and lemma length and ten and six metabolomic signals, respectively. These results would guide the selection of new germplasm lines of the three model grasses in ongoing genome-wide association studies.

Brachypodium distachyon - A Useful Model in the Qualification of Mutagen-Induced Micronuclei Using Multicolor FISH

Brachypodium distachyon (Brachypodium) is now intensively utilized as a model grass species in various biological studies. Its favorable cytological features create a unique foundation for a convenient system in mutagenesis, thereby potentially enabling the 'hot spots' and 'cold spots' of DNA damage in its genome to be analyzed. The aim of this study was to analyze the involvement of 5S rDNA, 25S rDNA, the Arabidopsis-type (TTTAGGG)n telomeric sequence and the Brachypodium-originated centromeric BAC clone CB33J12 in the micronuclei formation in Brachypodium root tip cells that were subjected to the chemical clastogenic agent maleic hydrazide (MH). To the best of our knowledge, this is the first use of a multicolor fluorescence in situ hybridization (mFISH) with four different DNA probes being used simultaneously to study plant mutagenesis. A quantitative analysis allowed ten types of micronuclei, which were characterized by the presence or absence of specific FISH signal(s), to be distinguished, thus enabling some specific rules governing the composition of the MH-induced micronuclei with the majority of them originating from the terminal regions of chromosomes, to be identified. The application of rDNA sequences as probes showed that 5S rDNA-bearing chromosomes are involved in micronuclei formation more frequently than the 25S rDNA-bearing chromosomes. These findings demonstrate the promising potential of Brachypodium to be a useful model organism to analyze the effects of various genotoxic agents on the plant nuclear genome stability, especially when the complex FISH-based and chromosome-specific approaches such as chromosome barcoding and chromosome painting will be applied in future studies.

Recreating Stable Brachypodium hybridum Allotetraploids by Uniting the Divergent Genomes of B. distachyon and B. stacei

Brachypodium hybridum (2n = 30) is a natural allopolyploid with highly divergent sub-genomes derived from two extant diploid species, B. distachyon (2n = 10) and B. stacei (2n = 20) that differ in chromosome evolution and number. We created synthetic B. hybridum allotetraploids by hybridizing various lines of B. distachyon and B. stacei. The initial amphihaploid F1 interspecific hybrids were obtained at low frequencies when B. distachyon was used as the maternal parent (0.15% or 0.245% depending on the line used) and were sterile. No hybrids were obtained from reciprocal crosses or when autotetraploids of the parental species were crossed. Colchicine treatment was used to double the genome of the F1 amphihaploid lines leading to allotetraploids. The genome-doubled F1 plants produced a few S1 (first selfed generation) seeds after self-pollination. S1 plants from one parental combination (Bd3-1×Bsta5) were fertile and gave rise to further generations whereas those of another parental combination (Bd21×ABR114) were sterile, illustrating the importance of the parental lineages crossed. The synthetic allotetraploids were stable and resembled the natural B. hybridum at the phenotypic, cytogenetic and genomic levels. The successful creation of synthetic B. hybridum offers the possibility to study changes in genome structure and regulation at the earliest stages of allopolyploid formation in comparison with the parental species and natural B. hybridum.

Molecular, genetic and evolutionary analysis of a paracentric inversion in Arabidopsis thaliana

Chromosomal inversions can provide windows onto the cytogenetic, molecular, evolutionary and demographic histories of a species. Here we investigate a paracentric 1.17-Mb inversion on chromosome 4 of Arabidopsis thaliana with nucleotide precision of its borders. The inversion is created by Vandal transposon activity, splitting an F-box and relocating a pericentric heterochromatin segment in juxtaposition with euchromatin without affecting the epigenetic landscape. Examination of the RegMap panel and the 1001 Arabidopsis genomes revealed more than 170 inversion accessions in Europe and North America. The SNP patterns revealed historical recombinations from which we infer diverse haplotype patterns, ancient introgression events and phylogenetic relationships. We find a robust association between the inversion and fecundity under drought. We also find linkage disequilibrium between the inverted region and the early flowering Col-FRIGIDA allele. Finally, SNP analysis elucidates the origin of the inversion to South-Eastern Europe approximately 5000 years ago and the FRI-Col allele to North-West Europe, and reveals the spreading of a single haplotype to North America during the 17th to 19th century. The 'American haplotype' was identified from several European localities, potentially due to return migration.

The arrangement of Brachypodium distachyon chromosomes in interphase nuclei

The spatial organization of chromatin within the interphase nucleus and the interactions between chromosome territories (CTs) are essential for various biological processes, such as DNA replication, transcription, and repair. However, detailed data about the CT arrangement in monocotyledonous plants are scarce. In this study, chromosome painting was used to analyse the distribution and associations of individual chromosomes in the 3-D preserved nuclei of Brachypodium distachyon root cells in order to determine the factors that may have an impact on the homologous CT arrangement. It was shown that the frequency of CT association is linked to the steric constraints imposed by the limited space within the nucleus and may depend on chromosome size and morphology as well as on the nuclear shape. Furthermore, in order to assess whether the distribution of interphase chromosomes is random or is subject to certain patterns, a comparison between the experimental data and the results of a computer simulation (ChroTeMo), which was based on a fully probabilistic distribution of the CTs, was performed. This comparison revealed that homologous chromosome arm CTs associate more often than if they were randomly arranged inside the interphase nucleus.

Chromosome Territory Modeller and Viewer

This paper presents ChroTeMo, a tool for chromosome territory modelling, accompanied by ChroTeVi–a chromosome territory visualisation software that uses the data obtained by ChroTeMo. These tools have been developed in order to complement the molecular cytogenetic research of interphase nucleus structure in a model grass Brachypodium distachyon. Although the modelling tool has been initially created for one particular species, it has universal application. The proposed version of ChroTeMo allows for generating a model of chromosome territory distribution in any given plant or animal species after setting the initial, species-specific parameters. ChroTeMo has been developed as a fully probabilistic modeller. Due to this feature, the comparison between the experimental data on the structure of a nucleus and the results obtained from ChroTeMo can indicate whether the distribution of chromosomes inside a nucleus is also fully probabilistic or is subjected to certain non-random patterns. The presented tools have been written in Python, so they are multiplatform, portable and easy to read. Moreover, if necessary they can be further developed by users writing their portions of code. The source code, documentation, and wiki, as well as the issue tracker and the list of related articles that use ChroTeMo and ChroTeVi, are accessible in a public repository at Github under GPL 3.0 license.

Chromosome identification for the carnivorous plant Genlisea margaretae

Genlisea margaretae, subgenus Genlisea, section Recurvatae (184 Mbp/1C), belongs to a plant genus with a 25-fold genome size difference and an extreme genome plasticity. Its 19 chromosome pairs could be distinguished individually by an approach combining optimized probe pooling and consecutive rounds of multicolor fluorescence in situ hybridization (mcFISH) with bacterial artificial chromosomes (BACs) selected for repeat-free inserts. Fifty-one BACs were assigned to 18 chromosome pairs. They provide a tool for future assignment of genomic sequence contigs to distinct chromosomes as well as for identification of homeologous chromosome regions in other species of the carnivorous Lentibulariaceae family, and potentially of chromosome rearrangements, in cases where more than one BAC per chromosome pair was identified.

Evolutionary and Taxonomic Implications of Variation in Nuclear Genome Size: Lesson from the Grass Genus Anthoxanthum (Poaceae)

The genus Anthoxanthum (sweet vernal grass, Poaceae) represents a taxonomically intricate polyploid complex with large phenotypic variation and its evolutionary relationships still poorly resolved. In order to get insight into the geographic distribution of ploidy levels and assess the taxonomic value of genome size data, we determined C- and Cx-values in 628 plants representing all currently recognized European species collected from 197 populations in 29 European countries. The flow cytometric estimates were supplemented by conventional chromosome counts. In addition to diploids, we found two low (rare 3x and common 4x) and one high (~16x-18x) polyploid levels. Mean holoploid genome sizes ranged from 5.52 pg in diploid A. alpinum to 44.75 pg in highly polyploid A. amarum, while the size of monoploid genomes ranged from 2.75 pg in tetraploid A. alpinum to 9.19 pg in diploid A. gracile. In contrast to Central and Northern Europe, which harboured only limited cytological variation, a much more complex pattern of genome sizes was revealed in the Mediterranean, particularly in Corsica. Eight taxonomic groups that partly corresponded to traditionally recognized species were delimited based on genome size values and phenotypic variation. Whereas our data supported the merger of A. aristatum and A. ovatum, eastern Mediterranean populations traditionally referred to as diploid A. odoratum were shown to be cytologically distinct, and may represent a new taxon. Autopolyploid origin was suggested for 4x A. alpinum. In contrast, 4x A. odoratum seems to be an allopolyploid, based on the amounts of nuclear DNA. Intraspecific variation in genome size was observed in all recognized species, the most striking example being the A. aristatum/ovatum complex. Altogether, our study showed that genome size can be a useful taxonomic marker in Anthoxathum to not only guide taxonomic decisions but also help resolve evolutionary relationships in this challenging grass genus.

Chromosome-Specific Painting in Cucumis Species Using Bulked Oligonucleotides

Chromosome-specific painting is a powerful technique in molecular cytogenetic and genome research. We developed an oligonucleotide (oligo)-based chromosome painting technique in cucumber (Cucumis sativus) that will be applicable in any plant species with a sequenced genome. Oligos specific to a single chromosome of cucumber were identified using a newly developed bioinformatic pipeline and then massively synthesized de novo in parallel. The synthesized oligos were amplified and labeled with biotin or digoxigenin for use in fluorescence in situ hybridization (FISH). We developed three different probes with each containing 23,000-27,000 oligos. These probes spanned 8.3-17 Mb of DNA on targeted cucumber chromosomes and had the densities of 1.5-3.2 oligos per kilobases. These probes produced FISH signals on a single cucumber chromosome and were used to paint homeologous chromosomes in other Cucumis species diverged from cucumber for up to 12 million years. The bulked oligo probes allowed us to track a single chromosome in early stages during meiosis. We were able to precisely map the pairing between cucumber chromosome 7 and chromosome 1 of Cucumis hystrix in a F1 hybrid. These two homeologous chromosomes paired in 71% of prophase I cells but only 25% of metaphase I cells, which may provide an explanation of the higher recombination rates compared to the chiasma frequencies between homeologous chromosomes reported in plant hybrids.