Introduction: Ribosomal DNA (rDNA) loci have been widely used for identification of allopolyploids and hybrids, although few of these studies employed high-throughput sequencing data. Here we use graph clustering implemented in the RepeatExplorer (RE) pipeline to analyze homoeologous 5S rDNA arrays at the genomic level searching for hybridogenic origin of species. Data were obtained from more than 80 plant species, including several well-defined allopolyploids and homoploid hybrids of different evolutionary ages and from widely dispersed taxonomic groups. Results: (i) Diploids show simple circular-shaped graphs of their 5S rDNA clusters. In contrast, most allopolyploids and other interspecific hybrids exhibit more complex graphs composed of two or more interconnected loops representing intergenic spacers (IGS). (ii) There was a relationship between graph complexity and locus numbers. (iii) The sequences and lengths of the 5S rDNA units reconstituted in silico from k-mers were congruent with those experimentally determined. (iv) Three-genomic comparative cluster analysis of reads from allopolyploids and progenitor diploids allowed identification of homoeologous 5S rRNA gene families even in relatively ancient (c. 1 Myr) Gossypium and Brachypodium allopolyploids which already exhibit uniparental partial loss of rDNA repeats. (v) Finally, species harboring introgressed genomes exhibit exceptionally complex graph structures. Conclusion: We found that the cluster graph shapes and graph parameters (k-mer coverage scores and connected component index) well-reflect the organization and intragenomic homogeneity of 5S rDNA repeats. We propose that the analysis of 5S rDNA cluster graphs computed by the RE pipeline together with the cytogenetic analysis might be a reliable approach for the determination of the hybrid or allopolyploid plant species parentage and may also be useful for detecting historical introgression events.

Nucleolar dominance (ND), initially described as ‘differential amphiplasty’, is a phenomenon observed in some plant and animal allopolyploids and hybrids in which the selective suppression of the activity of 35S rRNA gene loci that have been inherited from one of the two or more ancestral genomes occurs. Although more than 80 years have passed since the discovery of ND, there is still a significant lack in our understanding of the mechanisms that determine this phenomenon. Here, we aimed to investigate the epigenetic status of 35S rRNA gene loci in the monocotyledonous Brachypodium hybridum, which is an allotetraploid that has resulted from a cross between B. distachyon and B. stacei. We revealed that the repressed B. stacei-inherited rDNA loci are characterised by a high level of DNA methylation. The global hypomethylation of B. hybridum nuclear DNA induced by 5-azacytidine, however, seems to be insufficient for the transcriptional reactivation of these loci, which indicates that factors other than DNA methylation are behind the suppression of B. stacei-originated loci. We also showed that the transcriptionally active and silenced fractions of rRNA genes that had been inherited from B. distachyon occupy different domains within the chromocentres adjacent to the nucleolus, depending on their epigenetic status.

Nucleolar dominance (ND) is selective epigenetic silencing of 35-48S rDNA loci. In allopolyploids, it is frequently manifested at the cytogenetic level by the inactivation of nucleolar organiser region(s) (NORs) inherited from one or several evolutionary ancestors. Grasses are ecologically and economically one of the most important land plant groups, which have frequently evolved through hybridisation and polyploidisation events. Here we review common and unique features of ND phenomena in this monocot family from cytogenetic, molecular, and genomic perspectives. We highlight recent advances achieved by using an allotetraploid model grass, Brachypodium hybridum, where ND commonly occurs at a population level, and we cover modern genomic approaches that decipher structural features of core arrays of NORs.

Summary: Nucleolar dominance (ND) consists of the reversible silencing of 35S/45S rDNA loci inherited from one of the ancestors of an allopolyploid. The molecular mechanisms by which one ancestral rDNA set is selected for silencing remain unclear. We applied a combination of molecular (Southern blot hybridization and reverse‐transcription cleaved amplified polymorphic sequence analysis), genomic (analysis of variants) and cytogenetic (fluorescence in situ hybridization) approaches to study the structure, expression and epigenetic landscape of 35S rDNA in an allotetraploid grass that exhibits ND, Brachypodium hybridum (genome composition DDSS), and its putative progenitors, Brachypodium distachyon (DD) and Brachypodium stacei (SS). In progenitor genomes, B. stacei showed a higher intragenomic heterogeneity of rDNA compared with B. distachyon. In all studied accessions of B. hybridum, there was a reduction in the copy number of S homoeologues, which was accompanied by their inactive transcriptional status. The involvement of DNA methylation in CG and CHG contexts in the silencing of the S‐genome rDNA loci was revealed. In the B. hybridum allotetraploid, ND is stabilized towards the D‐genome units, irrespective of the polyphyletic origin of the species, and does not seem to be influenced by homoeologous 35S rDNA ratios and developmental stage.

Sequential immunolocalisation of 5-methylcytosine (5-MeC) and fluorescence in situ hybridisation with chromosome-specific BAC clones were performed on Brachypodium distachyon mitotic metaphase chromosomes to determine specific DNA methylation patterns of each chromosome in the complement. In the majority of cells examined, chromosomes Bd4 and Bd5, which bear the loci of 5S and 35S ribosomal DNA, respectively, had characteristic 5-MeC patterns. In contrast, the distribution of 5-MeC along the metacentric chromosome pairs Bd1, Bd2 and Bd3 was more variable. There were numerous differences in distribution of methylated sites between homologous chromosomes as well as between chromosome arms. Some chromosome sites, such as pericentromeric regions, were highly methylated in all chromosomes. Additionally, the influence of a hypomethylating agent, 5-azacytidine, on B. distachyon chromosome methylation patterns was confirmed. It was found that some chromosome pairs underwent demethylation more easily than others, but there was no apparent regularity in demethylation of particular chromosome segments.

Cytogenetics constitutes a branch of genetics that is focused on the cellular components, especially chromosomes, in relation to heredity and genome structure, function and evolution. The use of modern cytogenetic approaches and the latest microscopes with image acquisition and processing systems enables the simultaneous two- or three-dimensional, multicolour visualisation of both single-copy and highly-repetitive sequences in the plant genome. The data that is gathered using the cytogenetic methods in the phylogenetic background enable tracing the evolution of the plant genome that involve changes in: (i) genome sizes; (ii) chromosome numbers and morphology; (iii) the content of repetitive sequences and (iv) ploidy level. Modern cytogenetic approaches such as FISH using chromosome- and genome-specific probes have been widely used in studies of the evolution of diploids and the consequences of polyploidy. Nowadays, modern cytogenetics complements analyses in other fields of cell biology and constitutes the linkage between genetics, molecular biology and genomics.