BACKGROUND: Polyploidisation often results in genome rearrangements that may involve changes in both the single-copy sequences and the repetitive genome fraction. In this study, we performed a comprehensive comparative analysis of repetitive DNA, with a particular focus on ribosomal DNA (rDNA), in Brachypodium hybridum (2n = 4x = 30, subgenome composition DDSS), an allotetraploid resulting from a natural cross between two diploid species that resemble the modern B. distachyon (2n = 10; DD) and B. stacei (2n = 20; SS). Taking advantage of the recurrent origin of B. hybridum, we investigated two genotypes, Bhyb26 and ABR113, differing markedly in their evolutionary age (1.4 and 0.14 Mya, respectively) and which resulted from opposite cross directions. To identify the origin of rDNA loci we employed cytogenetic and molecular methods (FISH, gCAPS and Southern hybridisation), phylogenetic and genomic approaches. RESULTS: Unlike the general maintenance of doubled gene dosage in B. hybridum, the rRNA genes showed a remarkable tendency towards diploidisation at both locus and unit levels. While the partial elimination of 35S rDNA units occurred in the younger ABR113 lineage, unidirectional elimination of the entire locus was observed in the older Bhyb26 lineage. Additionally, a novel 5S rDNA family was amplified in Bhyb26 replacing the parental units. The 35S and 5S rDNA units were preferentially eliminated from the S- and D-subgenome, respectively. Thus, in the more ancient B. hybridum lineage, Bhyb26, 5S and 35S rRNA genes are likely expressed from different subgenomes, highlighting the complexity of polyploid regulatory networks. CONCLUSION: Comparative analyses between two B. hybridum lineages of distinct evolutionary ages revealed that although the recent lineage ABR113 exhibited an additive pattern of rDNA loci distribution, the ancient lineage Bhyb26 demonstrated a pronounced tendency toward diploidisation manifested by the reduction in the number of both 35S and 5S loci. In conclusion, the age of the allopolyploid appears to be a decisive factor in rDNA turnover in B. hybridum.

• Keywords
• Brachypodium hybridum, 35S rDNA IGS, FISH, 5S rDNA NTS, nrITS, rDNA loci,
• MeSH
• Brachypodium * genetics   MeSH
• Phylogeny * MeSH
• Genetic Variation MeSH
• Genome, Plant MeSH
• Genes, rRNA genetics   MeSH
• Evolution, Molecular * MeSH
• Polyploidy * MeSH
• DNA, Ribosomal genetics   MeSH
• RNA, Ribosomal genetics   MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• DNA, Ribosomal MeSH
• RNA, Ribosomal MeSH

Despite the widely accepted involvement of DNA methylation in the regulation of rDNA transcription, the relative participation of different cytosine methylation pathways is currently described only for a few model plants. Using PacBio, Bisulfite, and RNA sequencing; PCR; Southern hybridizations; and FISH, the epigenetic consequences of rDNA copy number variation were estimated in two T. porrifolius lineages, por1 and por2, the latter with more than twice the rDNA copy numbers distributed approximately equally between NORs on chromosomes A and D. The lower rDNA content in por1 correlated with significantly reduced (>90%) sizes of both D-NORs. Moreover, two (L and S) prominent rDNA variants, differing in the repetitive organization of intergenic spacers, were detected in por2, while only the S-rDNA variant was detected in por1. Transcriptional activity of S-rDNA in por1 was associated with secondary constriction of both A-NORs. In contrast, silencing of S-rDNA in por2 was accompanied by condensation of A-NORs, secondary constriction on D-NORs, and L-rDNA transcriptional activity, suggesting (i) bidirectional nucleolar dominance and (ii) association of S-rDNAs with A-NORs and L-rDNAs with D-NORs in T. porrifolius. Each S- and L-rDNA array was formed of several sub-variants differentiating both genetically (specific SNPs) and epigenetically (transcriptional efficiency and cytosine methylation). The most significant correlations between rDNA silencing and methylation were detected for symmetric CWG motifs followed by CG motifs. No correlations were detected for external cytosine in CCGs or asymmetric CHHs, where methylation was rather position-dependent, particularly for AT-rich variants. We conclude that variations in rDNA copy numbers in plant diploids can be accompanied by prompt epigenetic responses to maintain an appropriate number of active rDNAs. The methylation dynamics of CWGs are likely to be the most responsible for regulating silent and active rDNA states.

• Keywords
• 35S rDNA copy number variations, CCGs and CHHs, CWGs, Tragopogon porrifolius ssp. porrifolius, bidirectional nucleolar dominance, methylation dynamics of CGs, transcriptional silencing/activation,
• MeSH
• Chromosomes, Plant genetics   MeSH
• Cytosine * metabolism   MeSH
• Epigenesis, Genetic MeSH
• Transcription, Genetic MeSH
• DNA Methylation * MeSH
• Gene Expression Regulation, Plant MeSH
• DNA, Ribosomal * genetics   MeSH
• Gene Silencing * MeSH
• DNA Copy Number Variations MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• Cytosine * MeSH
• DNA, Ribosomal * MeSH

Chamaecrista is a Pantropical legume genus of the tribe Cassieae, which includes six other genera. In contrast to most of the other Cassieae genera, Chamaecrista shows significant variability in chromosome number (from 2n = 14 to 2n = 56), with small and morphologically similar chromosomes. Here, we performed a new cytomolecular analysis on chromosome number, genome size, and rDNA site distribution in a molecular phylogenetic perspective to interpret the karyotype trends of Chamaecrista and other two genera of Cassieae, seeking to understand their systematics and evolution. Our phylogenetic analysis revealed that Chamaecrista is monophyletic and can be divided into four major clades corresponding to the four sections of the genus. Chromosome numbers ranged from 2n = 14, 16 (section Chamaecrista) to 2n = 28 (sections Absus, Apoucouita, and Baseophyllum). The number of 5S and 35S rDNA sites varied between one and three pairs per karyotype, distributed on different chromosomes or in synteny, with no obvious phylogenetic significance. Our data allowed us to propose x = 7 as the basic chromosome number of Cassieae, which was changed by polyploidy generating x = 14 (sections Absus, Apoucouita, and Baseophyllum) and by ascending dysploidy to x = 8 (section Chamaecrista). The DNA content values supported this hypothesis, with the genomes of the putative tetraploids being larger than those of the putative diploids. We hypothesized that ascending dysploidy, polyploidy, and rDNA amplification/deamplification are the major events in the karyotypic diversification of Chamaecrista. The chromosomal marks characterized here may have cytotaxonomic potential in future studies.

• Keywords
• trnL-trnF, C-value, Genome size, ITS1–5.8S–ITS2, Karyotype evolution, Ribosomal DNA,
• MeSH
• Chamaecrista * genetics   MeSH
• Chromosomes, Plant genetics   MeSH
• Fabaceae * genetics   MeSH
• Phylogeny MeSH
• Genome, Plant MeSH
• Karyotype MeSH
• Polyploidy MeSH
• DNA, Ribosomal genetics   MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• DNA, Ribosomal MeSH

Genome or genomic dominance (GD) is a phenomenon observed in hybrids when one parental genome becomes dominant over the other. It is manifested by the replacement of chromatin of the submissive genome by that of the dominant genome and by biased gene expression. Nucleolar dominance (ND) - the functional expression of only one parental set of ribosomal genes in hybrids - is another example of an intragenomic competitive process which, however, concerns ribosomal DNA only. Although GD and ND are relatively well understood, the nature and extent of their potential interdependence is mostly unknown. Here, we ask whether hybrids showing GD also exhibit ND and, if so, whether the dominant genome is the same. To test this, we used hybrids between Festuca and Lolium grasses (Festulolium), and between two Festuca species in which GD has been observed (with Lolium as the dominant genome in Festulolium and F. pratensis in interspecific Festuca hybrids). Using amplicon sequencing of ITS1 and ITS2 of the 45S ribosomal DNA (rDNA) cluster and molecular cytogenetics, we studied the organization and expression of rDNA in leaf tissue in five hybrid combinations, four generations and 31 genotypes [F. pratensis × L. multiflorum (F1, F2, F3, BC1), L. multiflorum × F. pratensis (F1), L. multiflorum × F. glaucescens (F2), L. perenne × F. pratensis (F1), F. glaucescens × F. pratensis (F1)]. We have found that instant ND occurs in Festulolium, where expression of Lolium-type rDNA reached nearly 100% in all F1 hybrids and was maintained through subsequent generations. Therefore, ND and GD in Festulolium are manifested by the same dominant genome (Lolium). We also confirmed the concordance between GD and ND in an interspecific cross between two Festuca species.

• Keywords
• Festuca, Lolium, fluorescent in situ hybridization, genome dominance, genomic in situ hybridization, internal transcribed spacer, nucleolar dominance, ribosomal DNA,
• Publication type
• Journal Article MeSH

The classical model of concerted evolution states that hundreds to thousands of ribosomal DNA (rDNA) units undergo homogenization, making the multiple copies of the individual units more uniform across the genome than would be expected given mutation frequencies and gene redundancy. While the universality of this over 50-year-old model has been confirmed in a range of organisms, advanced high throughput sequencing techniques have also revealed that rDNA homogenization in many organisms is partial and, in rare cases, even apparently failing. The potential underpinning processes leading to unexpected intragenomic variation have been discussed in a number of studies, but a comprehensive understanding remains to be determined. In this work, we summarize information on variation or polymorphisms in rDNAs across a wide range of taxa amongst animals, fungi, plants, and protists. We discuss the definition and description of concerted evolution and describe whether incomplete concerted evolution of rDNAs predominantly affects coding or non-coding regions of rDNA units and if it leads to the formation of pseudogenes or not. We also discuss the factors contributing to rDNA variation, such as interspecific hybridization, meiotic cycles, rDNA expression status, genome size, and the activity of effector genes involved in genetic recombination, epigenetic modifications, and DNA editing. Finally, we argue that a combination of approaches is needed to target genetic and epigenetic phenomena influencing incomplete concerted evolution, to give a comprehensive understanding of the evolution and functional consequences of intragenomic variation in rDNA.

• MeSH
• Phylogeny MeSH
• Genetic Variation * MeSH
• Fungi genetics   MeSH
• Evolution, Molecular MeSH
• Mutation MeSH
• Polymorphism, Genetic * MeSH
• DNA, Ribosomal genetics   MeSH
• Animals MeSH
• Check Tag
• Animals MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Review MeSH
• Names of Substances
• DNA, Ribosomal MeSH

The history of rDNA research started almost 90 years ago when the geneticist, Barbara McClintock observed that in interphase nuclei of maize the nucleolus was formed in association with a specific region normally located near the end of a chromosome, which she called the nucleolar organizer region (NOR). Cytologists in the twentieth century recognized the nucleolus as a common structure in all eukaryotic cells, using both light and electron microscopy and biochemical and genetic studies identified ribosomes as the subcellular sites of protein synthesis. In the mid- to late 1960s, the synthesis of nuclear-encoded rRNA was the only system in multicellular organisms where transcripts of known function could be isolated, and their synthesis and processing could be studied. Cytogenetic observations of NOR regions with altered structure in plant interspecific hybrids and detailed knowledge of structure and function of rDNA were prerequisites for studies of nucleolar dominance, epistatic interactions of rDNA loci, and epigenetic silencing. In this article, we focus on the early rDNA research in plants, performed mainly at the dawn of molecular biology in the 60 to 80-ties of the last century which presented a prequel to the modern genomic era. We discuss - from a personal view - the topics such as synthesis of rRNA precursor (35S pre-rRNA in plants), processing, and the organization of 35S and 5S rDNA. Cloning and sequencing led to the observation that the transcribed and processed regions of the rRNA genes vary enormously, even between populations and species, in comparison with the more conserved regions coding for the mature rRNAs. Epigenetic phenomena and the impact of hybridization and allopolyploidy on rDNA expression and homogenization are discussed. This historical view of scientific progress and achievements sets the scene for the other articles highlighting the immense progress in rDNA research published in this special issue of Frontiers in Plant Science on "Molecular organization, evolution, and function of ribosomal DNA."

• Keywords
• epigenetics, hybridization, molecular evolution, nucleolar dominance, polyploidy, rDNA research history, rRNA precursor, rRNA processing,
• Publication type
• Journal Article MeSH
• Review MeSH

Nuclear ribosomal DNA (nrDNA) has displayed extraordinary dynamics during the evolution of plant species. However, the patterns and evolutionary significance of nrDNA array expansion or contraction are still relatively unknown. Moreover, only little is known of the fate of minority nrDNA copies acquired between species via horizontal transfer. The barley genus Hordeum (Poaceae) represents a good model for such a study, as species of section Stenostachys acquired nrDNA via horizontal transfer from at least five different panicoid genera, causing long-term co-existence of native (Hordeum-like) and non-native (panicoid) nrDNAs. Using quantitative PCR, we investigated copy number variation (CNV) of nrDNA in the diploid representatives of the genus Hordeum. We estimated the copy number of the foreign, as well as of the native ITS types (ribotypes), and followed the pattern of their CNV in relation to the genus' phylogeny, species' genomes size and the number of nrDNA loci. For the native ribotype, we encountered an almost 19-fold variation in the mean copy number among the taxa analysed, ranging from 1689 copies (per 2C content) in H. patagonicum subsp. mustersii to 31342 copies in H. murinum subsp. glaucum. The copy numbers did not correlate with any of the genus' phylogeny, the species' genome size or the number of nrDNA loci. The CNV was high within the recognised groups (up to 13.2 × in the American I-genome species) as well as between accessions of the same species (up to 4×). Foreign ribotypes represent only a small fraction of the total number of nrDNA copies. Their copy numbers ranged from single units to tens and rarely hundreds of copies. They amounted, on average, to between 0.1% (Setaria ribotype) and 1.9% (Euclasta ribotype) of total nrDNA. None of the foreign ribotypes showed significant differences with respect to phylogenetic groups recognised within the sect. Stenostachys. Overall, no correlation was found between copy numbers of native and foreign nrDNAs suggesting the sequestration and independent evolution of native and non-native nrDNA arrays. Therefore, foreign nrDNA in Hordeum likely poses a dead-end by-product of horizontal gene transfer events.

• Keywords
• copy number variation (CNV), fluorescent in situ hybridisation (FISH), horizontal gene transfer (HGT), internal transcribed spacer (ITS), nuclear ribosomal DNA (nrDNA), phylogeny, qPCR (quantitative PCR), xenolog,
• Publication type
• Journal Article MeSH

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.

• Keywords
• Brachypodium hybridum, 35S rDNA evolution, 35S rRNA gene expression, allopolyploidy, nucleolar dominance,
• MeSH
• Brachypodium genetics   metabolism   MeSH
• Chromosomes, Plant genetics   MeSH
• Genetic Loci genetics   MeSH
• Genome, Plant genetics   MeSH
• Genes, rRNA genetics   MeSH
• DNA Methylation genetics   MeSH
• Evolution, Molecular MeSH
• Polymorphism, Genetic genetics   MeSH
• Genes, Plant genetics   MeSH
• Blotting, Southern MeSH
• Tetraploidy * MeSH
• DNA Copy Number Variations genetics   MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH

Plant genomes vary greatly in composition and size mainly due to the diversity of repetitive DNAs and the inherent propensity for their amplification and removal from the host genome. Most studies addressing repeatome dynamics focus on model organisms, whereas few provide comprehensive investigations across the genomes of related taxa. Herein, we analyze the evolution of repeats of the 13 species in Melampodium sect. Melampodium, representing all but two of its diploid taxa, in a phylogenetic context. The investigated genomes range in size from 0.49 to 2.27 pg/1C (ca. 4.5-fold variation), despite having the same base chromosome number (x = 10) and very strong phylogenetic affinities. Phylogenetic analysis performed in BEAST and ancestral genome size reconstruction revealed mixed patterns of genome size increases and decreases across the group. High-throughput genome skimming and the RepeatExplorer pipeline were utilized to determine the repeat families responsible for the differences in observed genome sizes. Patterns of repeat evolution were found to be highly correlated with phylogenetic position, namely taxonomic series circumscription. Major differences found were in the abundances of the SIRE (Ty1-copia), Athila (Ty3-gypsy), and CACTA (DNA transposon) lineages. Additionally, several satellite DNA families were found to be highly group-specific, although their overall contribution to genome size variation was relatively small. Evolutionary changes in repetitive DNA composition and genome size were complex, with independent patterns of genome up- and downsizing throughout the evolution of the analyzed diploids. A model-based analysis of genome size and repetitive DNA composition revealed evidence for strong phylogenetic signal and differential evolutionary rates of major lineages of repeats in the diploid genomes.

• Keywords
• Bayesian analysis, Melampodium, ancestral state reconstruction, genome size, phylogenetics, repetitive DNA, tandem repeats, transposable elements,
• Publication type
• Journal Article MeSH

Allopolyploidy has played an important role in the evolution of the flowering plants. Genome mergers are often accompanied by significant and rapid alterations of genome size and structure via chromosomal rearrangements and altered dynamics of tandem and dispersed repetitive DNA families. Recent developments in sequencing technologies and bioinformatic methods allow for a comprehensive investigation of the repetitive component of plant genomes. Interpretation of evolutionary dynamics following allopolyploidization requires both the knowledge of parentage and the age of origin of an allopolyploid. Whereas parentage is typically inferred from cytogenetic and phylogenetic data, age inference is hampered by the reticulate nature of the phylogenetic relationships. Treating subgenomes of allopolyploids as if they belonged to different species (i.e., no recombination among subgenomes) and applying cross-bracing (i.e., putting a constraint on the age difference of nodes pertaining to the same event), we can infer the age of allopolyploids within the framework of the multispecies coalescent within BEAST2. Together with a comprehensive characterization of the repetitive DNA fraction using the RepeatExplorer pipeline, we apply the dating approach in a group of closely related allopolyploids and their progenitor species in the plant genus Melampodium (Asteraceae). We dated the origin of both the allotetraploid, Melampodium strigosum, and its two allohexaploid derivatives, Melampodium pringlei and Melampodium sericeum, which share both parentage and the direction of the cross, to the Pleistocene ($<$1.4 Ma). Thus, Pleistocene climatic fluctuations may have triggered formation of allopolyploids possibly in short intervals, contributing to difficulties in inferring the precise temporal order of allopolyploid species divergence of M. sericeum and M. pringlei. The relatively recent origin of the allopolyploids likely played a role in the near-absence of major changes in the repetitive fraction of the polyploids' genomes. The repetitive elements most affected by the postpolyploidization changes represented retrotransposons of the Ty1-copia lineage Maximus and, to a lesser extent, also Athila elements of Ty3-gypsy family.

• MeSH
• Asteraceae classification   genetics   MeSH
• DNA, Plant genetics   MeSH
• Phylogeny MeSH
• Genome, Plant genetics   MeSH
• Evolution, Molecular * MeSH
• Polyploidy MeSH
• Repetitive Sequences, Nucleic Acid genetics   MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• DNA, Plant MeSH