Cytogenetic features of rRNA genes across land plants: analysis of the Plant rDNA database

Karyotype Description and Comparative Chromosomal Mapping of 5S rDNA in 42 Species

This study was conducted to evaluate the 5S rDNA site number, position, and origin of signal pattern diversity in 42 plant species using fluorescence in situ hybridization. The species were selected based on the discovery of karyotype rearrangement, or because 5S rDNA had not yet been explored the species. The chromosome number varied from 14 to 160, and the chromosome length ranged from 0.63 to 6.88 μm, with 21 species having small chromosomes (<3 μm). The chromosome numbers of three species and the 5S rDNA loci of nineteen species are reported for the first time. Six 5S rDNA signal pattern types were identified. The 5S rDNA varied and was abundant in signal site numbers (2-18), positions (distal, proximal, outside of chromosome arms), and even in signal intensity. Variation in the numbers and locations of 5S rDNA was observed in 20 species, whereas an extensive stable number and location of 5S rDNA was found in 22 species. The potential origin of the signal pattern diversity was proposed and discussed. These data characterized the variability of 5S rDNA within the karyotypes of the 42 species that exhibited chromosomal rearrangements and provided anchor points for genetic physical maps.

Genome Variability in Artificial Allopolyploid Hybrids of Avena sativa L. and Avena macrostachya Balansa ex Coss. et Durieu Based on Marker Sequences of Satellite DNA and the ITS1-5.8S rDNA Region

Artificial hybrids between cultivated Avena species and wild Avena macrostachya that possess genes for resistance to biotic and abiotic stresses can be important for oat breeding. For the first time, a comprehensive study of genomes of artificial fertile hybrids Avena sativa × Avena macrostachya and their parental species was carried out based on the chromosome FISH mapping of satellite DNA sequences (satDNAs) and also analysis of intragenomic polymorphism in the 18S-ITS1-5.8S rDNA region, using NGS data. Chromosome distribution patterns of marker satDNAs allowed us to identify all chromosomes in the studied karyotypes, determine their subgenomic affiliation, and detect several chromosome rearrangements. Based on the obtained cytogenomic data, we revealed differences between two A. macrostachya subgenomes and demonstrated that only one of them was inherited in the studied octoploid hybrids. Ribotype analyses showed that the second major ribotype of A. macrostachya was species-specific and was not represented in rDNA pools of the octoploids, which could be related to the allopolyploid origin of this species. Our results indicate that the use of marker satDNAs in cytogenomic studies can provide important data on genomic relationships within Avena allopolyploid species and hybrids, and also expand the potential for interspecific crosses for breeding.

The Application of Fluorescence In Situ Hybridization in the Prescreening of Veronica Hybrids

Fluorescence in situ hybridization (FISH), a molecular cytogenetic technique that enables the visualization and identification of specific DNA sequences within chromosomes, has emerged as a pivotal tool in plant breeding programs, particularly in the case of Veronica species. Veronica, a genus with a complex reproductive system, often poses challenges in accurately identifying hybrids because of its tendency to hybridize, which leads to intricate genetic variation. This study focused on the use of FISH as a prescreening method to identify true hybrids in Veronica breeding programs. FISH analysis was first performed on the parents to identify their 45S and 5S rDNA signals, along with their respective chromosome numbers. The signals were then compared with those of the twenty progenies with reference to their supposed parents. Five true hybrids, seven self-pollinated progenies, and eight false hybrids were identified through FISH. The findings highlight the significance of FISH as a screening method that contributes significantly to the efficiency of Veronica breeding programs by ensuring the preservation of desired genetic traits and minimizing the inadvertent inclusion of misidentified hybrids. To conclude, this study underscores the vital role of FISH in enhancing the precision and success of breeding programs and opens new avenues for improved breeding strategies and crop development.

Epigenetic histone H3 phosphorylation marks discriminate between univalent- and bivalent-forming chromosomes during canina asymmetrical meiosis

BACKGROUND AND AIMS

Dogroses (Rosa sect. Caninae) are mostly pentaploid, bearing 2n = 5x = 35 chromosomes in somatic cells. They evolved a unique form of asymmetrical meiosis characterized by two types of chromosomes: (1) chromosomes forming bivalents and distributed in the normal sexual way; and (2) chromosomes occurring as univalents and transferred by a female gamete only. In the mature pollen of pentaploid species, seven bivalent-derived chromosomes are transmitted to offspring, and 21 unpaired univalent chromosomes are eliminated during microsporogenesis. To discriminate between bivalent- and univalent-forming chromosomes, we studied histone H3 phosphorylation patterns regulating meiotic chromosome condensation and segregation.

METHODS

We analysed histone modification patterns during male canina meiosis in two representative dogrose species, 5x Rosa canina and 5x Rosa rubiginosa, by immunohistochemical and molecular cytogenetics approaches. Immunostaining of meiotic cells included α-tubulin, histone H3 phosphorylation (H3S10p, H3S28p and H3T3p) and methylation (H3K4me3 and H3K27me3) marks. In addition, fluorescent in situ hybridization was carried out with an 18S rDNA probe.

KEY RESULTS

In the first meiotic division, univalent chromosomes underwent equational division into chromatids, while homologues in bivalents were segregated as regular dyads. In diakinesis, bivalent chromosomes displayed strong H3 phosphorylation signals in proximal regions, spreading to the rest of the chromosome. In contrast, in univalents, the H3 phosphorylation signals were weaker, occurring mostly outside proximal regions largely overlapping with the H3K4me3 signals. Reduced phosphorylation was associated with relative under-condensation of the univalent chromosomes, particularly at early diakinesis.

CONCLUSIONS

We hypothesize that the absence of pairing and/or recombination in univalent chromosomes negatively affects the histone H3 phosphorylation of their chromatin and perhaps the loading of meiotic-specific cohesins. This apparently destabilizes cohesion of sister chromatids, leading to their premature split in the first meiotic division.

Molecular and cytogenetic analyses in Geranium macrorrhizum L. wild Italian plants

Geranium macrorrhizum L. is a herbaceous species native to southern Europe and was introduced in central Europe and North America. It is also widely distributed in Italy. In this study, molecular and cytogenetic analyses were carried out on 22 wild plants, collected in central and southern Italy, compared with five cultivated plants, with the main purpose to identify those living near the Marmore waterfalls in central Italy, recently described as the new species Geranium lucarinii. Four barcoding markers (rbcL, matK, trnH-psbA intergenic spacer and internal transcribed spacer region) were sequenced and their variability among the plants was evaluated. Chromosome numbers were determined and 45S rDNA was physically mapped by fluorescence in situ hybridization. Moreover, genomic affinity between wild and cultivated plants was evaluated by genomic in situ hybridization. The results of this study supported that all the plants belong to G. macrorrhizum, including the Marmore population. Barcoding analyses showed a close similarity among the wild plants, and a differentiation, although not significant, between the wild plants on one hand and the cultivated plants on the other. Integrated studies focusing on morphological, genetic and ecological characterization of a larger number of wild populations would allow us to know the extent of the variability within the species.

The Dynamic Interplay Between Ribosomal DNA and Transposable Elements: A Perspective From Genomics and Cytogenetics

Although both are salient features of genomes, at first glance ribosomal DNAs and transposable elements are genetic elements with not much in common: whereas ribosomal DNAs are mainly viewed as housekeeping genes that uphold all prime genome functions, transposable elements are generally portrayed as selfish and disruptive. These opposing characteristics are also mirrored in other attributes: organization in tandem (ribosomal DNAs) versus organization in a dispersed manner (transposable elements); evolution in a concerted manner (ribosomal DNAs) versus evolution by diversification (transposable elements); and activity that prolongs genomic stability (ribosomal DNAs) versus activity that shortens it (transposable elements). Re-visiting relevant instances in which ribosomal DNA-transposable element interactions have been reported, we note that both repeat types share at least four structural and functional hallmarks: (1) they are repetitive DNAs that shape genomes in evolutionary timescales, (2) they exchange structural motifs and can enter co-evolution processes, (3) they are tightly controlled genomic stress sensors playing key roles in senescence/aging, and (4) they share common epigenetic marks such as DNA methylation and histone modification. Here, we give an overview of the structural, functional, and evolutionary characteristics of both ribosomal DNAs and transposable elements, discuss their roles and interactions, and highlight trends and future directions as we move forward in understanding ribosomal DNA-transposable element associations.

High plasticity of ribosomal DNA organization in budding yeast

In eukaryotic genomes, rDNA generally resides as a highly repetitive and dynamic structure, making it difficult to study. Here, a synthetic rDNA array on chromosome III in budding yeast was constructed to serve as the sole source of rRNA. Utilizing the loxPsym site within each rDNA repeat and the Cre recombinase, we were able to reduce the copy number to as few as eight copies. Additionally, we constructed strains with two or three rDNA arrays and found that the presence of multiple arrays did not affect the formation of a single nucleolus. Although alteration of the position and number of rDNA arrays did impact the three-dimensional genome structure, the additional rDNA arrays had no deleterious influence on cell growth or transcriptomes. Overall, this study sheds light on the high plasticity of rDNA organization and opens up opportunities for future rDNA engineering.

Genomic and cytogenetic analyses reveal satellite repeat signature in allotetraploid okra (Abelmoschus esculentus)

BACKGROUND

Satellite repeats are one of the most rapidly evolving components in eukaryotic genomes and play vital roles in genome regulation, genome evolution, and speciation. As a consequence, the composition, abundance and chromosome distribution of satellite repeats often exhibit variability across various species, genome, and even individual chromosomes. However, we know little about the satellite repeat evolution in allopolyploid genomes.

RESULTS

In this study, we investigated the satellite repeat signature in five okra (Abelmoschus esculentus) accessions using genomic and cytogenetic methods. In each of the five accessions, we identified eight satellite repeats, which exhibited a significant level of intraspecific conservation. Through fluorescence in situ hybridization (FISH) experiments, we observed that the satellite repeats generated multiple signals and exhibited variations in copy number across chromosomes. Intriguingly, we found that five satellite repeats were interspersed with centromeric retrotransposons, signifying their involvement in centromeric satellite repeat identity. We confirmed subgenome-biased amplification patterns of these satellite repeats through existing genome assemblies or dual-color FISH, indicating their distinct dynamic evolution in the allotetraploid okra subgenome. Moreover, we observed the presence of multiple chromosomes harboring the 35 S rDNA loci, alongside another chromosomal pair carrying the 5 S rDNA loci in okra using FISH assay. Remarkably, the intensity of 35 S rDNA hybridization signals varied among chromosomes, with the signals predominantly localized within regions of relatively weak DAPI staining, associated with GC-rich heterochromatin regions. Finally, we observed a similar localization pattern between 35 S rDNA and three satellite repeats with high GC content and confirmed their origin in the intergenic spacer region of the 35 S rDNA.

CONCLUSIONS

Our findings uncover a unique satellite repeat signature in the allotetraploid okra, contributing to our understanding of the composition, abundance, and chromosomal distribution of satellite repeats in allopolyploid genomes, further enriching our understanding of their evolutionary dynamics in complex allopolyploid genomes.

High-fidelity (repeat) consensus sequences from short reads using combined read clustering and assembly

BACKGROUND

Despite the many cheap and fast ways to generate genomic data, good and exact genome assembly is still a problem, with especially the repeats being vastly underrepresented and often misassembled. As short reads in low coverage are already sufficient to represent the repeat landscape of any given genome, many read cluster algorithms were brought forward that provide repeat identification and classification. But how can trustworthy, reliable and representative repeat consensuses be derived from unassembled genomes?

RESULTS

Here, we combine methods from repeat identification and genome assembly to derive these robust consensuses. We test several use cases, such as (1) consensus building from clustered short reads of non-model genomes, (2) from genome-wide amplification setups, and (3) specific repeat-centred questions, such as the linked vs. unlinked arrangement of ribosomal genes. In all our use cases, the derived consensuses are robust and representative. To evaluate overall performance, we compare our high-fidelity repeat consensuses to RepeatExplorer2-derived contigs and check, if they represent real transposable elements as found in long reads. Our results demonstrate that it is possible to generate useful, reliable and trustworthy consensuses from short reads by a combination from read cluster and genome assembly methods in an automatable way.

CONCLUSION

We anticipate that our workflow opens the way towards more efficient and less manual repeat characterization and annotation, benefitting all genome studies, but especially those of non-model organisms.

Cyto-molecular characterization of rDNA and chromatin composition in the NOR-associated satellite in Chestnut (Castanea spp.)

The American chestnut (Castanea dentata, 2n = 2x = 24), once known as the "King of the Appalachian Forest", was decimated by chestnut blight during the first half of the twentieth century by an invasive fungus (Cryphonectria parasitica). The Chinese chestnut (C. mollissima, 2n = 2x = 24), in contrast to American chestnut, is resistant to this blight. Efforts are being made to transfer this resistance to American chestnut through backcross breeding and genetic engineering. Both chestnut genomes have been genetically mapped and recently sequenced to facilitate gene discovery efforts aimed at assisting molecular breeding and genetic engineering. To complement and extend this genomic work, we analyzed the distribution and organization of their ribosomal DNAs (35S and 5S rDNA), and the chromatin composition of the nucleolus organizing region (NOR)-associated satellites. Using fluorescent in situ hybridization (FISH), we have identified two 35S (one major and one minor) and one 5S rDNA sites. The major 35S rDNA sites are terminal and sub-terminal in American and Chinese chestnuts, respectively, originating at the end of the short arm of the chromosome, extending through the secondary constriction and into the satellites. An additional 5S locus was identified in certain Chinese chestnut accessions, and it was linked distally to the major 35S site. The NOR-associated satellite in Chinese chestnut was found to comprise a proximal region packed with 35S rDNA and a distinct distal heterochromatic region. In contrast, the American chestnut satellite was relatively small and devoid of the distal heterochromatic region.

Abundance of distal repetitive DNA sequences in Capsicum L. (Solanaceae) chromosomes

Chili peppers (Solanaceae family) have great commercial value. They are commercialized in natura and used as spices and for ornamental and medicinal purposes. Although three whole genomes have been published, limited information about satellite DNA sequences, their composition, and genomic distribution has been provided. Here, we exploited the noncoding repetitive fraction, represented by satellite sequences, that tends to accumulate in blocks along chromosomes, especially near the chromosome ends of peppers. Two satellite DNA sequences were identified (CDR-1 and CDR-2), characterized and mapped in silico in three Capsicum genomes (C. annuum, C. chinense, and C. baccatum) using data from the published high-coverage sequencing and repeats finding bioinformatic tools. Localization using FISH in the chromosomes of these species and in two others (C. frutescens and C. chacoense), totaling five species, showed signals adjacent to the rDNA sites. A sequence comparison with existing Solanaceae repeats showed that CDR-1 and CDR-2 have different origins but without homology to rDNA sequences. Satellites occupied subterminal chromosomal regions, sometimes collocated with or adjacent to 35S rDNA sequences. Our results expand knowledge about the diversity of subterminal regions of Capsicum chromosomes, showing different amounts and distributions within and between karyotypes. In addition, these sequences may be useful for future phylogenetic studies.

Intragenomic rDNA variation - the product of concerted evolution, mutation, or something in between?

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.

The Role of Repetitive Sequences in Repatterning of Major Ribosomal DNA Clusters in Lepidoptera

Genes for major ribosomal RNAs (rDNA) are present in multiple copies mainly organized in tandem arrays. The number and position of rDNA loci can change dynamically and their repatterning is presumably driven by other repetitive sequences. We explored a peculiar rDNA organization in several representatives of Lepidoptera with either extremely large or numerous rDNA clusters. We combined molecular cytogenetics with analyses of second- and third-generation sequencing data to show that rDNA spreads as a transcription unit and reveal association between rDNA and various repeats. Furthermore, we performed comparative long read analyses among the species with derived rDNA distribution and moths with a single rDNA locus, which is considered ancestral. Our results suggest that satellite arrays, rather than mobile elements, facilitate homology-mediated spread of rDNA via either integration of extrachromosomal rDNA circles or ectopic recombination. The latter arguably better explains preferential spread of rDNA into terminal regions of lepidopteran chromosomes as efficiency of ectopic recombination depends on the proximity of homologous sequences to telomeres.

Developing population identification tool based on polymorphism of rDNA for traditional Chinese medicine: Artemisia annua L

BACKGROUND

Artemisia annua, a well-known traditional Chinese medicine, is the main source for production of artemisinin, an anti-malaria drug. A. annua is distributed globally, with great diversity of morphological characteristics and artemisinin contents. Diverse traits among A. annua populations impeded the stable production of artemisinin, which needs an efficient tool to identify strains and assess population genetic homogeneity.

PURPOSE

In this study, ribosomal DNA (rDNA), were characterized for A. annua for strains identification and population genetic homogeneity assessment.

METHODS

The ribosomal RNA (rRNA) genes were identified using cmscan and assembled using rDNA unit of LQ-9 as a reference. rDNA among Asteraceae species were compared performing with 45S rDNA. The rDNA copy number was calculated based on sequencing depth. The polymorphisms of rDNA sequences were identified with bam-readcount, and confirmed by Sanger sequencing and restriction enzyme experiment. The ITS2 amplicon sequencing was used to verify the stability of ITS2 haplotype analysis.

RESULTS

Different from other Asteraceae species, 45S and 5S linked-type rDNA was only found in Artemisia genus. Rich polymorphisms of copy number and sequence of rDNA were identified in A. annua population. The haplotype composition of internal transcribed spacer 2 (ITS2) region which had moderate sequence polymorphism and relative short size was significantly different among A. annua strains. A population discrimination method was developed based on ITS2 haplotype analysis with high-throughput sequencing.

CONCLUSION

This study provides comprehensive characteristics of rDNA and suggests that ITS2 haplotype analysis is ideal tool for A. annua strain identification and population genetic homogeneity assessment.

Single Copies of the 5S rRNA Inserted into 45S rDNA Intergenic Spacers in the Genomes of Nototheniidae (Perciformes, Actinopterygii)

In the vast majority of Animalia genomes, the 5S rRNA gene repeats are located on chromosomes outside of the 45S rDNA arrays of the nucleolar organiser (NOR). We analysed the genomic databases available and found that a 5S rDNA sequence is inserted into the intergenic spacer (IGS) between the 45S rDNA repeats in ten species of the family Nototheniidae (Perciformes, Actinopterigii). We call this sequence the NOR-5S rRNA gene. Along with Testudines and Crocodilia, this is the second case of a close association between four rRNA genes within one repetitive unit in deuterostomes. In both cases, NOR-5S is oriented opposite the 45S rDNA. None of the three nucleotide substitutions compared to the canonical 5S rRNA gene influenced the 5S rRNA secondary structure. In transcriptomes of the Patagonian toothfish, we only found NOR-5S rRNA reads in ovaries and early embryos, but not in testis or somatic tissues of adults. Thus, we consider the NOR-5S gene to be a maternal-type 5S rRNA template. The colocalization of the 5S and 45S ribosomal genes appears to be essential for the equimolar production of all four rRNAs in the species that show rDNA amplification during oogenesis. Most likely, the integration of 5S and NOR rRNA genes occurred prior to Nototheniidae lineage diversification.

High Chromosomal Reorganization and Presence of Microchromosomes in Chactidae Scorpions from the Brazilian Amazon

Scorpions are of particular interest in cytogenomic studies, as they can present a high incidence of chromosomal rearrangements heterozygous in natural populations. In this study, we cytogenetically analyzed four species of Chactidae. In Brotheas, 2n = 40 was observed in Brotheas silvestris, 2n = 48 in Brotheas paraensis, and 2n = 50 (cytotype A) or 2n = 52 (cytotype B) among populations of Brotheas amazonicus. Our results showed a bimodal karyotype in Neochactas parvulus, 2n = 54, with microchromosomes and a concentration of constitutive heterochromatin in macrochromosomes. The 45S rDNA is located in only one pair of the karyotype, with different heteromorphisms of clusters of this rDNA in the cytotype B of B. amazonicus, with NOR-bearing chromosomes involved in multi-chromosomal associations during meiosis I. The U2 snDNA was mapped in the interstitial region of distinct karyotype pairs of three Chactidae species. Our results indicate the possible formation of cryptic species in B. amazonicus; the different 45S rDNA configurations in the genome of this species may result from amplification and degeneration. We suggest that the bimodal karyotype in N. parvulus results from fusion/fission events and that the unequal distribution of repetitive DNAs between macro and microchromosomes contributes to the maintenance of its asymmetry.

The origin of cultivated mangosteen (Garcinia mangostana L. var. mangostana): Critical assessments and an evolutionary-ecological perspective

Mangosteen (Garcinia mangostana var. mangostana) is a popular tropical fruit, yet many aspects of its biology and evolutionary history are little known. Its origin remains contentious, although recent findings suggest G. mangostana L. var. malaccensis (Hook. f.) Nazre (synonym: G. malaccensis Hook. f.) as the sole progenitor. We review hypotheses on the origin of mangosteen and clarify points that have been affected by errors of fact and interpretation. The narrow focus and lack of detail in published results make their interpretation difficult. When possible, we support our interpretations with field observations and examination of herbarium specimens. We outline the main biological traits (e.g., dioecy, facultative apomixis, and polyploidy) of mangosteen and its wild relatives to infer traits that might have evolved during domestication of mangosteen. We find no clear indication that apomixis and polyploidy evolved during domestication. Polyploidy is known in the wild relatives, but apomixis has not yet been demonstrated. Also, we propose a testable new evolutionary-ecological framework that we call "Forest-Dusun Interface" to infer processes in the origin of mangosteen. Dusun (Malay) refers to subsistence orchards in this context. Lastly, we propose future studies to address identified knowledge gaps.

Polyploidy on islands - concerted evolution and gene loss amid chromosomal stasis

BACKGROUND AND AIMS

Polyploidy is an important process that often generates genomic diversity within lineages, but it can also cause changes that result in loss of genomic material. Island lineages, while often polyploid, typically show chromosomal stasis but have not been investigated in detail regarding smaller-scale gene loss. Our aim was to investigate post-polyploidization genome dynamics in a chromosomally stable lineage of Malvaceae endemic to New Zealand.

METHODS

We determined chromosome numbers and used fluorescence in situ hybridization to localize 18S and 5S rDNA. Gene sequencing of 18S rDNA, the internal transcribed spacers (ITS) with intervening 5.8S rDNA, and a low-copy nuclear gene, GBSSI-1, was undertaken to determine if gene loss occurred in the New Zealand lineage following polyploidy.

KEY RESULTS

The chromosome number for all species investigated was 2n = 42, with the first published report for the monotypic Australian genus Asterotrichion. The five species investigated all had two 5S rDNA signals localized interstitially on the long arm of one of the largest chromosome pairs. All species, except Plagianthus regius, had two 18S rDNA signals localized proximally on the short arm of one of the smallest chromosome pairs. Plagianthus regius had two additional 18S rDNA signals on a separate chromosome, giving a total of four. Sequencing of nuclear ribosomal 18S rDNA and the ITS cistron indicated loss of historical ribosomal repeats. Phylogenetic analysis of a low-copy nuclear gene, GBSSI-1, indicated that some lineages maintained three copies of the locus, while others have lost one or two copies.

CONCLUSIONS

Although island endemic lineages show chromosomal stasis, with no additional changes in chromosome number, they may undergo smaller-scale processes of gene loss and concerted evolution ultimately leading to further genome restructuring and downsizing.

Release 4.0 of the Plant rDNA Database: A Database on Plant Ribosomal DNA Loci Number, Their Position, and Organization: An Information Source for Comparative Cytogenetics

This paper presents the latest update to the Plant rDNA database (Release 4.0), a valuable resource for researchers in the field of plant cytogenetics. The database provides information on the number, position, and arrangement of ribosomal DNA loci in plants, including angiosperms, gymnosperms, bryophytes, and pteridophytes. The new release includes new data for 820 species coming from additional 173 papers. In the updated version of the Plant rDNA database, 4948 entries comprising 2760 organisms can be found. A brief guide on how to navigate the database and obtain the desired information is also provided. The regular updating of the database is important for ensuring the information it contains is accurate, up-to-date, and useful for the research community. The Plant rDNA database continues to be beneficial for phylogenetic and cytogenetic studies in a wide range of taxa including angiosperms, gymnosperms, and early diverging groups, such as bryophytes and lycophytes.

Analysis of 5S rDNA Genomic Organization Through the RepeatExplorer2 Pipeline: A Simplified Protocol

The ribosomal RNA genes (rDNA) are universal genome components with a housekeeping function, given the crucial role of ribosomal RNA in the synthesis of ribosomes and thus for life-on-Earth. Therefore, their genomic organization is of considerable interest for biologists, in general. Ribosomal RNA genes have also been largely used to establish phylogenetic relationships, and to identify allopolyploid or homoploid hybridization.Here, we demonstrate how high-throughput sequencing data, through graph clustering implemented in RepeatExplorer2 pipeline ( https://repeatexplorer-elixir.cerit-sc.cz/galaxy/ ), can be helpful to decipher the genomic organization of 5S rRNA genes. We show that the linear shapes of cluster graphs are reminiscent to the linked organization of 5S and 35S rDNA (L-type arrangement) while the circular graphs correspond to their separate arrangement (S-type). We further present a simplified protocol based on the paper by (Garcia et al., Front Plant Sci 11:41, 2020) about the use of graph clustering of 5S rDNA homoeologs (S-type) to identify hybridization events in the species history. We found that the graph complexity (i.e., graph circularity in this case) is related to ploidy and genome complexity, with diploids typically showing circular-shaped graphs while allopolyploids and other interspecific hybrids display more complex graphs, with usually two or more interconnected loops representing intergenic spacers. When a three-genomic comparative clustering analysis from a given hybrid (homoploid/allopolyploid) and its putative progenitor species (diploids) is performed, it is possible to identify the corresponding homoeologous 5S rRNA gene families, and to elucidate the contribution of each putative parental genome to the 5S rDNA pool of the hybrid. Thus, the analysis of 5S rDNA cluster graphs by RepeatExplorer, together with information coming from other sources (e.g., morphology, cytogenetics) is a complementary approach for the determination of allopolyploid or homoploid hybridization and even ancient introgression events.

The Use of Ribosomal DNA for Comparative Cytogenetics

Fluorescence in situ hybridization (FISH) with ribosomal DNA (rDNA) sequences provides excellent chromosome markers for comparative cytogenetic analyses, especially in non-model plant species. The tandem repeat nature of a sequence and the presence of a highly conserved genic region make rDNA sequences relatively easy to isolate and clone. In this chapter, we describe the use of rDNA as markers for comparative cytogenetics studies. Traditionally, cloned probes labeled with Nick-translation have been used to detect rDNA loci. Recently, pre-labeled oligonucleotides are also employed quite frequently to detect both 35S and 5S rDNA loci. Ribosomal DNA sequences, together with other DNA probes in FISH/GISH or with fluorochromes such as CMA₃ banding or silver staining, are very useful tools in comparative analyses of plant karyotypes.

Hybridization and diversity of aquatic macrophyte Sparganium L. (Typhaceae) as revealed by high-throughput nrDNA sequencing

Sparganium is an emergent aquatic macrophyte widely spread in temperate and subtropical zones. Taxa of this genus feature high phenotypic plasticity and can produce interspecific hybrids. By means of high-throughput sequencing of the internal transcribed spacer (ITS1) of 35S rDNA, the status of 15 Eurasian Sparganium species and subspecies was clarified and the role of hybridization events in the recent evolution of the genus was investigated. It has been shown that a number of species such as S. angustifolium, S. fallax and S. subglobosum have homogenized rDNA represented by one major ribotype. The rDNA of other taxa is represented by two or more major ribotypes. Species with high rDNA heterogeneity are apparently of hybrid origin. Based on the differences in rDNA patterns, intraspecific diversity was identified in S. probatovae and S. emersum. Thus, we have concluded that Sparganium has extensive interspecific hybridization at the subgenus level, and there may also be occasional hybridization between species from different subgenera.

Agro-Morphological and Cytogenetic Characterization of Colchicine-Induced Tetraploid Plants of Polemonium caeruleum L. (Polemoniaceae)

Polemonium caeruleum L. (Polemoniaceae) is a valuable medicinal herb with a wide spectrum of biological activities. Under natural conditions, the productivity of this species is rather low. In this study, colchicine-induced tetraploid plants (2n = 4x = 36) of P. caeruleum were obtained, and for the first time, their morphological and cytogenetic characterization was performed. In the tetraploid plants, raw material productivity and also the content of triterpene saponins were significantly higher than in the control diploids. The analysis of chromosome behavior at meiosis and FISH chromosome mapping of 45S and 5S rDNA generally demonstrated stability of both genomes in the tetraploid plants. Based on chromosome morphology and distribution patterns of the studied molecular cytogenetic markers, all chromosome pairs in karyotypes were identified, and chromosome karyograms and idiograms of P. caeruleum were constructed. The revealed specific microdiagnostic characteristics of P. caeruleum (strongly sinuous cells and anomocytic stomata of the leaf epidermis, and also glandular hairs along the veins) could be useful for raw material identification. In the obtained tetraploids, the predominance of large stomata on the lower leaf epidermis was determined. The studied tetraploids can be used in various breeding programs to obtain high-quality pharmaceutical raw materials of P. caeruleum.

Karyotype Analysis, Genomic and Fluorescence In Situ Hybridization (GISH and FISH) Reveal the Ploidy and Parental Origin of Chromosomes in Paeonia Itoh Hybrids

Itoh hybrids are intersectional hybrids in Paeonia L. with sect. Moutan and sect. Paeonia as paternal and maternal parents, respectively. Therefore, these hybrids have herbaceous stems with improved ornamental value introduced by the paternal parent. Although both of their parents are diploids, Itoh hybrids are triploids. Moreover, the parental origin of their chromosomes has not been extensively studied. This study systematically analyzed the genome size, ploidy, and karyotype of Itoh hybrids and compared them with their parental taxa. Although the monoploid genome size of Itoh hybrids was different, it was not significantly different from that of the parents. However, the size of varieties in the two parental taxa was significantly different from the wild species, probably due to genome rearrangements caused by artificial selection. Further karyotype analysis, correlation analysis, and hierarchical clustering could not identify the parental origin of chromosomes in Itoh hybrids. Verification through genomic and fluorescence in situ hybridization (GISH and FISH) suggested that for the three sets of chromosomes in Itoh hybrids, two were from the paternal parent, and one was from the maternal parent. One of the first two sets was from wild species, and the other from a cultivated variety. GISH could not label the chromosomes of cultivated peonies from the sect. Moutan, probably due to the huge and complex genomes compared with the wild species. Meanwhile, 5S rDNA-based FISH was first applied in Paeonia, which may be used for ploidy assessment. This work may give insights into the utilization of Itoh hybrid resources.

The Chromosome Number and rDNA Loci Evolution in Onobrychis (Fabaceae)

The evolution of chromosome number and ribosomal DNA (rDNA) loci number and localisation were studied in Onobrychis Mill. Diploid and tetraploid species, as well as two basic chromosome numbers, x = 7 and x = 8, were observed among analysed taxa. The chromosomal distribution of rDNA loci was presented here for the first time using fluorescence in situ hybridisation (FISH) with 5S and 35S rDNA probes. Onobrychis species showed a high polymorphism in the number and localisation of rDNA loci among diploids, whereas the rDNA loci pattern was very similar in polyploids. Phylogenetic relationships among the species, inferred from nrITS sequences, were used as a framework to reconstruct the patterns of basic chromosome number and rDNA loci evolution. Analysis of the evolution of the basic chromosome numbers allowed the inference of x = 8 as the ancestral number and the descending dysploidy and polyploidisation as the major mechanisms of the chromosome number evolution. Analyses of chromosomal patterns of rRNA gene loci in a phylogenetic context resulted in the reconstruction of one locus of 5S rDNA and one locus of 35S rDNA in the interstitial chromosomal position as the ancestral state in this genus.

Comparative Analyses of Complete Chloroplast Genomes and Karyotypes of Allotetraploid Iris koreana and Its Putative Diploid Parental Species (Iris Series Chinenses, Iridaceae)

The Iris series Chinenses in Korea comprises four species (I. minutoaurea, I. odaesanensis, I. koreana, and I. rossii), and the group includes some endangered species, owing to their high ornamental, economic, and conservation values. Among them, the putative allotetraploid, Iris koreana (2n = 4x = 50), is hypothesized to have originated from the hybridization of the diploids I. minutoaurea (2n = 2x = 22) and I. odaesanensis (2n = 2x = 28) based on morphological characters, chromosome numbers, and genome size additivity. Despite extensive morphological and molecular phylogenetical studies on the genus Iris, little is known about Korean irises in terms of their complete chloroplast (cp) genomes and molecular cytogenetics that involve rDNA loci evolution based on fluorescence in situ hybridization (FISH). This study reports comparative analyses of the karyotypes of the three Iris species (I. koreana, I. odaesanensis, and I. minutoaurea), with an emphasis on the 5S and 35S rDNA loci number and localization using FISH together with the genome size and chromosome number. Moreover, the cp genomes of the same individuals were sequenced and assembled for comparative analysis. The rDNA loci numbers, which were localized consistently at the same position in all species, and the chromosome numbers and genome size values of tetraploid Iris koreana (four 5S and 35S loci; 2n = 50; 1C = 7.35 pg) were additively compared to its putative diploid progenitors, I. minutoaurea (two 5S and 35S loci; 2n = 22; 1C = 3.71 pg) and I. odaesanensis (two 5S and 35S loci; 2n = 28; 1C = 3.68 pg). The chloroplast genomes were 152,259–155,145 bp in length, and exhibited a conserved quadripartite structure. The Iris cp genomes were highly conserved and similar to other Iridaceae cp genomes. Nucleotide diversity analysis indicated that all three species had similar levels of genetic variation, but the cp genomes of I. koreana and I. minutoaurea were more similar to each other than to I. odaesanensis. Positive selection was inferred for psbK and ycf2 genes of the three Iris species. Phylogenetic analyses consistently recovered I. odaesanensis as a sister to a clade containing I. koreana and I. minutoaurea. Although the phylogenetic relationship, rDNA loci number, and localization, together with the genome size and chromosome number of the three species, allowed for the inference of I. minutoaurea as a putative maternal taxon and I. odaesanensis as a paternal taxon, further analyses involving species-specific molecular cytogenetic markers and genomic in situ hybridization are required to interpret the mechanisms involved in the origin of the chromosomal variation in Iris series Chinenses. This study contributes towards the genomic and chromosomal evolution of the genus Iris.

Structure and Methylation of 35S rDNA in Allopolyploids Anemone multifida (2n = 4x = 32, BBDD) and Anemone baldensis (2n = 6x = 48, AABBDD) and Their Parental Species Show Evidence of Nucleolar Dominance

Transcriptional silencing of 35S rDNA loci inherited from one parental species is occurring relatively frequently in allopolyploids. However, molecular mechanisms by which it is selected for transcriptional silencing remain unclear. We applied NGS, silver staining and bisulfite sequencing to study the structure, expression and methylation landscape of 35S rDNA in two allopolyploids of common origin, allotetraploid Anemone multifida (2n = 4x = 32, genome composition BBDD) and allohexaploid A. baldensis (2n = 6x = 48, AABBDD), and their genome donors, A. sylvestris (2n = 16, AA), A. cylindrica (2n = 16, BB) and A. parviflora (2n = 16, DD). The size of the recovered 35S rDNA units varied from 10,489 bp in A. cylindrica to 12,084 bp in A. sylvestris. Anemone showed an organization typical of most ribosomal 35S rDNA composed of NTS, ETS, rRNA genes, TTS and TIS with structural features of plant IGS sequences and all functional elements needed for rRNA gene activity. The NTS was more variable than the ETS and consisted of SRs which are highly variable among Anemone. Five to six CpG-rich islands were found within the ETS. CpG island located adjacent to the transcription initiation site (TIS) was highly variable regarding the sequence size and methylation level and exhibited in most of the species lower levels of methylation than CpG islands located adjacent to the 18S rRNA gene. Our results uncover hypomethylation of A. sylvestris- and A. parviflora-derived 35S rDNA units in allopolyploids A. multifida and A. baldensis. Hypomethylation of A. parviflora-derived 35S rDNA was more prominent in A. baldensis than in A. multifida. We showed that A. baldensis underwent coupled A. sylvestris-derived 35S rDNA array expansion and A. parviflora-derived 35S rDNA copy number decrease that was accompanied by lower methylation level of A. sylvestris-derived 35S rDNA units in comparison to A. parviflora-derived 35S rDNA units. These observations suggest that in A. baldensis nucleolar dominance is directed toward A. sylvestris-derived chromosomes. This work broadens our current knowledge of the 35S rDNA organization in Anemone and provides evidence of the progenitor-specific 35S rDNA methylation in nucleolar dominance.

Evolutionary Dynamics of the Repeatome Explains Contrasting Differences in Genome Sizes and Hybrid and Polyploid Origins of Grass Loliinae Lineages

The repeatome is composed of diverse families of repetitive DNA that keep signatures on the historical events that shaped the evolution of their hosting species. The cold seasonal Loliinae subtribe includes worldwide distributed taxa, some of which are the most important forage and lawn species (fescues and ray-grasses). The Loliinae are prone to hybridization and polyploidization. It has been observed a striking two-fold difference in genome size between the broad-leaved (BL) and fine-leaved (FL) Loliinae diploids and a general trend of genome reduction of some high polyploids. We have used genome skimming data to uncover the composition, abundance, and potential phylogenetic signal of repetitive elements across 47 representatives of the main Loliinae lineages. Independent and comparative analyses of repetitive sequences and of 5S rDNA loci were performed for all taxa under study and for four evolutionary Loliinae groups [Loliinae, Broad-leaved (BL), Fine-leaved (FL), and Schedonorus lineages]. Our data showed that the proportion of the genome covered by the repeatome in the Loliinae species was relatively high (average ∼ 51.8%), ranging from high percentages in some diploids (68.7%) to low percentages in some high-polyploids (30.7%), and that changes in their genome sizes were likely caused by gains or losses in their repeat elements. Ty3-gypsy Retand and Ty1-copia Angela retrotransposons were the most frequent repeat families in the Loliinae although the relatively more conservative Angela repeats presented the highest correlation of repeat content with genome size variation and the highest phylogenetic signal of the whole repeatome. By contrast, Athila retrotransposons presented evidence of recent proliferations almost exclusively in the Lolium clade. The repeatome evolutionary networks showed an overall topological congruence with the nuclear 35S rDNA phylogeny and a geographic-based structure for some lineages. The evolution of the Loliinae repeatome suggests a plausible scenario of recurrent allopolyploidizations followed by diploidizations that generated the large genome sizes of BL diploids as well as large genomic rearrangements in highly hybridogenous lineages that caused massive repeatome and genome contractions in the Schedonorus and Aulaxyper polyploids. Our study has contributed to disentangling the impact of the repeatome dynamics on the genome diversification and evolution of the Loliinae grasses.

Uneven Levels of 5S and 45S rDNA Site Number and Loci Variations across Wild Chrysanthemum Accessions

Ribosomal DNA (rDNA) is an excellent cytogenetic marker owing to its tandem arrangement and high copy numbers. However, comparative studies have focused more on the number of rDNA site variations within the Chrysanthemum genus, and studies on the types of rDNA sites with the same experimental procedures at the species levels are lacking. To further explore the number and types of rDNA site variations, we combined related data to draw ideograms of the rDNA sites of Chrysanthemum accessions using oligonucleotide fluorescence in situ hybridization (Oligo-FISH). Latent variations (such as polymorphisms of 45S rDNA sites and co-localized 5S-45S rDNA) also occurred among the investigated accessions. Meanwhile, a significant correlation was observed between the number of 5S rDNA sites and chromosome number. Additionally, the clumped and concentrated geographical distribution of different ploidy Chrysanthemum accessions may significantly promote the karyotype evolution. Based on the results above, we identified the formation mechanism of rDNA variations. Furthermore, these findings may provide a reliable method to examine the sites and number of rDNA variations among Chrysanthemum and its related accessions and allow researchers to further understand the evolutionary and phylogenetic relationships of the Chrysanthemum genus.

New Insights into the Genomic Structure of Avena L.: Comparison of the Divergence of A-Genome and One C-Genome Oat Species

We used next-generation sequencing analysis of the 3′-part of 18S rDNA, ITS1, and a 5′-part of the 5.8S rDNA region to understand genetic variation among seven diploid A-genome Avena species. We used 4−49 accessions per species that represented the As genome (A. atlantica, A. hirtula, and wiestii), Ac genome (A. canariensis), Ad genome (A. damascena), Al genome (A. longiglumis), and Ap genome (A. prostrata). We also took into our analysis one C-genome species, A. clauda, which previously was found to be related to A-genome species. The sequences of 169 accessions revealed 156 haplotypes of which seven haplotypes were shared by two to five species. We found 16 ribotypes that consisted of a unique sequence with a characteristic pattern of single nucleotide polymorphisms and deletions. The number of ribotypes per species varied from one in A. longiglumis to four in A. wiestii. Although most ribotypes were species-specific, we found two ribotypes shared by three species (one for A. damascena, A. hirtula, and A. wiestii, and the second for A. longiglumis, A. atlantica, and A. wiestii), and a third ribotype shared between A. atlantica and A. wiestii. A characteristic feature of the A. clauda ribotype, a diploid C-genome species, is that two different families of ribotypes have been found in this species. Some of these ribotypes are characteristic of Cc-genome species, whereas others are closely related to As-genome ribotypes. This means that A. clauda can be a hybrid between As- and C-genome oats.

Genomic architecture of 5S rDNA cluster and its variations within and between species

BACKGROUND

Ribosomal DNAs (rDNAs) are arranged in purely tandem repeats, preventing them from being reliably assembled onto chromosomes during generation of genome assembly. The uncertainty of rDNA genomic structure presents a significant barrier for studying their function and evolution.

RESULTS

Here we generate ultra-long Oxford Nanopore Technologies (ONT) and short NGS reads to delineate the architecture and variation of the 5S rDNA cluster in the different strains of C. elegans and C. briggsae. We classify the individual rDNA's repeating units into 25 types based on the unique sequence variations in each unit of C. elegans (N2). We next perform assembly of the cluster by taking advantage of the long reads that carry these units, which led to an assembly of 5S rDNA cluster consisting of up to 167 consecutive 5S rDNA units in the N2 strain. The ordering and copy number of various rDNA units are consistent with the separation time between strains. Surprisingly, we observed a drastically reduced level of variation in the unit composition in the 5S rDNA cluster in the C. elegans CB4856 and C. briggsae AF16 strains than in the C. elegans N2 strain, suggesting that N2, a widely used reference strain, is likely to be defective in maintaining the 5S rDNA cluster stability compared with other wild isolates of C. elegans or C. briggsae.

CONCLUSIONS

The results demonstrate that Nanopore DNA sequencing reads are capable of generating assembly of highly repetitive sequences, and rDNA units are highly dynamic both within and between population(s) of the same species in terms of sequence and copy number. The detailed structure and variation of the 5S rDNA units within the rDNA cluster pave the way for functional and evolutionary studies.

Molecular and Cytogenetic Analysis of rDNA Evolution in Crepis Sensu Lato

Although Crepis was the first model plant group in which chromosomal changes were considered to play an important role in speciation, their chromosome structure and evolution have been barely investigated using molecular cytogenetic methods. The aim of the study was to provide a better understanding of the patterns and directions of Crepis chromosome evolution, using comparative analyses of rDNA loci number and localisation. The chromosome base number and chromosomal organisation of 5S and 35S rDNA loci were analysed in the phylogenetic background for 39 species of Crepis, which represent the evolutionary lineages of Crepis sensu stricto and Lagoseris, including Lapsana communis. The phylogenetic relationships among all the species were inferred from nrITS and newly obtained 5S rDNA NTS sequences. Despite high variations in rDNA loci chromosomal organisation, most species had a chromosome with both rDNA loci within the same (usually short) chromosomal arm. The comparative analyses revealed several independent rDNA loci number gains and loci repositioning that accompanied diversification and speciation in Crepis. Some of the changes in rDNA loci patterns were reconstructed for the same evolutionary lineages as descending dysploidy.

Distribution patterns of rDNA loci in the Schedonorus-Lolium complex (Poaceae)

The Schedonorus-Lolium complex of the subtribe Loliinae (Poaceae) includes several economically important forage and turf grasses. This complex encompasses Lolium Linnaeus, 1753, Festuca Linnaeus, 1753 subgenus Schedonorus (P. Beauvois, 1824) Petermann, 1849 and Micropyropsis Romero Zarco et Cabezudo, 1983. New FISH results of 5S and 18S-26S rDNA sequences are presented for three species and the results are interpreted in a review of distribution patterns of 5S and 18S-26S rDNA sequences among other species in the complex. Micropyropsistuberosa Romero Zarco et Cabezudo, 1983 (2n = 2x = 14) displayed a distribution pattern of rDNA sequences identical to that of F.pratensis Hudson, 1762, supporting a close phylogenetic relationship at the bottom of the phylogenetic tree. "Loliummultiflorum" Lamarck, 1779 accessions sourced from Morocco showed a different pattern from European L.multiflorum and could be a unique and previously uncharacterised taxon. North African Festucasimensis Hochstetter ex A. Richard, 1851 had a marker pattern consistent with allotetraploidy and uniparental loss of one 18S-26S rDNA locus. This allotetraploid has previously been suggested to have originated from a hybrid with Festucaglaucescens (Festucaarundinaceavar.glaucescens Boissier, 1844). However, the distribution patterns of the two rDNA sequences in this allotetraploid do not align with F.glaucescens, suggesting that its origin from this species is unlikely. Furthermore, comparisons with other higher alloploids in the complex indicate that F.simensis was a potential donor of two sub-genomes of allohexaploid Festucagigantea (Linnaeus) Villars, 1787. In the overall complex, the proximal locations of both rDNA markers were conserved among the diploid species. Two types of synteny of the two markers could, to a considerable extent, distinguish allo- and autogamous Lolium species. The ancestral parentage of the three Festuca allotetraploids has not yet been determined, but all three appear to have been sub-genome donors to the higher allopolypoids of sub-genus Schedonorus. Terminal locations of both the markers were absent from the diploids but were very frequently observed in the polyploids.

The Ribosomal DNA Loci of the Ancient Monocot Pistia stratiotes L. (Araceae) Contain Different Variants of the 35S and 5S Ribosomal RNA Gene Units

The freshwater plant water lettuce (Pistia stratiotes L.) grows in warm climatic zones and is used for phytoremediation and biomass production. P. stratiotes belongs to the Araceae, an ecologically and structurally diverse early monocot family, but the phylogenetic relationships among Araceae members are poorly understood. Ribosomal DNAs (rDNAs), including the 35S and 5S rDNA, encode the RNA components of ribosomes and are widely used in phylogenetic and evolutionary studies of various plant taxa. Here, we comprehensively characterized the chromosomal locations and molecular organization of 35S and 5S rDNA genes in water lettuce using karyological and molecular methods. Fluorescence in situ hybridization revealed a single location for the 35S and 5S rDNA loci, each on a different pair of the species' 28 chromosomes. Molecular cloning and nucleotide sequencing of 35S rDNA of P. stratiotes, the first representative Araceae sensu stricto in which such a study was performed, displayed typical structural characteristics. The full-length repeat showed high sequence conservation of the regions producing the 18S, 5.8S, and 25S rRNAs and divergence of the internal transcribed spacers ITS1 and ITS2 as well as the large intergenic spacer (IGS). Alignments of the deduced sequence of 18S rDNA with the sequences available for other Araceae and representatives of other clades were used for phylogenetic analysis. Examination of 11 IGS sequences revealed significant intra-genomic length variability due to variation in subrepeat number, with four types of units detected within the 35S rDNA locus of the P. stratiotes genome (estimated size 407 Mb/1C). Similarly, the 5S rDNA locus harbors gene units comprising a conserved 119-bp sequence encoding 5S rRNA and two types of non-transcribed spacer (NTS) sequences. Type I was classified into four subtypes, which apparently originated via progressive loss of subrepeats within the duplicated NTS region containing the 3' part of the 5S rRNA gene. The minor Type II NTS is shorter than Type I and differs in nucleotide composition. Some DNA clones containing two or three consecutive 5S rDNA repeats harbored 5S rDNA genes with different types of NTSs, confirming the mosaic composition of the 5S rDNA locus.

First discovered, long out of sight, finally visible: ribosomal DNA

With the advent of long-read sequencing, previously unresolvable genomic elements are being revisited in an effort to generate fully complete reference genomes. One such element is ribosomal DNA (rDNA), the highly conserved genomic region that encodes rRNAs. Genomic structure and content of the rDNA are variable in both prokarya and eukarya, posing interesting questions about the biology of rDNA. Here, we consider the types of variation observed in rDNA - including locus structure and number, copy number, and sequence variation - and their known phenotypic consequences. With recent advances in long-read sequencing technology, incorporating the full rDNA sequence into reference genomes is within reach. This knowledge will have important implications for understanding rDNA biology within the context of cell physiology and whole-organism phenotypes.

Genome structure and diversity among Cynanchum wilfordii accessions

BACKGROUND

Cynanchum wilfordii (Cw) and Cynanchum auriculatum (Ca) have long been used in traditional medicine and as functional food in Korea and China, respectively. They have diverse medicinal functions, and many studies have been conducted, including pharmaceutical efficiency and metabolites. Especially, Cw is regarded as the most famous medicinal herb in Korea due to its menopausal symptoms relieving effect. Despite the high demand for Cw in the market, both species are cultivated using wild resources with rare genomic information.

RESULTS

We collected 160 Cw germplasm from local areas of Korea and analyzed their morphological diversity. Five Cw and one Ca of them, which were morphologically diverse, were sequenced, and nuclear ribosomal DNA (nrDNA) and complete plastid genome (plastome) sequences were assembled and annotated. We investigated the genomic characteristics of Cw as well as the genetic diversity of plastomes and nrDNA of Cw and Ca. The Cw haploid nuclear genome was approximately 178 Mbp. Karyotyping revealed the juxtaposition of 45S and 5S nrDNA on one of 11 chromosomes. Plastome sequences revealed 1226 interspecies polymorphisms and 11 Cw intraspecies polymorphisms. The 160 Cw accessions were grouped into 21 haplotypes based on seven plastome markers and into 108 haplotypes based on seven nuclear markers. Nuclear genotypes did not coincide with plastome haplotypes that reflect the frequent natural outcrossing events.

CONCLUSIONS

Cw germplasm had a huge morphological diversity, and their wide range of genetic diversity was revealed through the investigation with 14 molecular markers. The morphological and genomic diversity, chromosome structure, and genome size provide fundamental genomic information for breeding of undomesticated Cw plants.

The Nuclear 35S rDNA World in Plant Systematics and Evolution: A Primer of Cautions and Common Misconceptions in Cytogenetic Studies

The ubiquitous presence of rRNA genes in nuclear, plastid, and mitochondrial genomes has provided an opportunity to use genomic markers to infer patterns of molecular and organismic evolution as well as to assess systematic issues throughout the tree of life. The number, size, location, and activity of the 35S rDNA cistrons in plant karyotypes have been used as conventional cytogenetic landmarks. Their scrutiny has been useful to infer patterns of chromosomal evolution and the data have been used as a proxy for assessing species discrimination, population differentiation and evolutionary relationships. The correct interpretation of rDNA markers in plant taxonomy and evolution is not free of drawbacks given the complexities derived from the lability of the genetic architecture, the diverse patterns of molecular change, and the fate and evolutionary dynamics of the rDNA units in hybrids and polyploid species. In addition, the terminology used by independent authors is somewhat vague, which often complicates comparisons. To date, no efforts have been reported addressing the potential problems and limitations involved in generating, utilizing, and interpreting the data from the 35S rDNA in cytogenetics. This review discusses the main technical and conceptual limitations of these rDNA markers obtained by cytological and karyological experimental work, in order to clarify biological and evolutionary inferences postulated in a systematic and phylogenetic context. Also, we provide clarification for some ambiguity and misconceptions in terminology usually found in published work that may help to improve the usage of the 35S ribosomal world in plant evolution.

The Evolution of Cytogenetic Traits in Cuscuta (Convolvulaceae), the Genus With the Most Diverse Chromosomes in Angiosperms

Karyotypes are characterized by traits such as chromosome number, which can change through whole-genome duplication and dysploidy. In the parasitic plant genus Cuscuta (Convolvulaceae), chromosome numbers vary more than 18-fold. In addition, species of this group show the highest diversity in terms of genome size among angiosperms, as well as a wide variation in the number and distribution of 5S and 35S ribosomal DNA (rDNA) sites. To understand its karyotypic evolution, ancestral character state reconstructions were performed for chromosome number, genome size, and position of 5S and 35S rDNA sites. Previous cytogenetic data were reviewed and complemented with original chromosome counts, genome size estimates, and rDNA distribution assessed via fluorescence in situ hybridization (FISH), for two, seven, and 10 species, respectively. Starting from an ancestral chromosome number of x = 15, duplications were inferred as the prevalent evolutionary process. However, in holocentric clade (subgenus Cuscuta), dysploidy was identified as the main evolutionary mechanism, typical of holocentric karyotypes. The ancestral genome size of Cuscuta was inferred as approximately 1C = 12 Gbp, with an average genome size of 1C = 2.8 Gbp. This indicates an expansion of the genome size relative to other Convolvulaceae, which may be linked to the parasitic lifestyle of Cuscuta. Finally, the position of rDNA sites varied mostly in species with multiple sites in the same karyotype. This feature may be related to the amplification of rDNA sites in association to other repeats present in the heterochromatin. The data suggest that different mechanisms acted in different subgenera, generating the exceptional diversity of karyotypes in Cuscuta.

Personal Perspectives on Plant Ribosomal RNA Genes Research: From Precursor-rRNA to Molecular Evolution

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."

45S rDNA Repeats of Turtles and Crocodiles Harbor a Functional 5S rRNA Gene Specifically Expressed in Oocytes

In most eukaryotic genomes, tandemly repeated copies of 5S rRNA genes are clustered outside the nucleolus organizer region (NOR), which normally encodes three other major rRNAs: 18S, 5.8S, and 28S. Our analysis of turtle rDNA sequences has revealed a 5S rDNA insertion into the NOR intergenic spacer in antisense orientation. The insertion (hereafter called NOR-5S rRNA gene) has a length of 119 bp and coexists with the canonical 5S rDNA clusters outside the NOR. Despite the ∼20% nucleotide difference between the two 5S gene sequences, their internal control regions for RNA polymerase III are similar. Using the turtle Trachemys scripta as a model species, we showed the NOR-5S rDNA specific expression in oocytes. This expression is concurrent with the NOR rDNA amplification during oocyte growth. We show that in vitellogenic oocytes, the NOR-5S rRNA prevails over the canonical 5S rRNA in the ribosomes, suggesting a role of modified ribosomes in oocyte-specific translation. The orders Testudines and Crocodilia seem to be the only taxa of vertebrates with such a peculiar rDNA organization. We speculate that the amplification of the 5S rRNA genes as a part of the NOR DNA during oogenesis provides a dosage balance between transcription of all the four ribosomal RNAs while producing a maternal pool of extra ribosomes. We further hypothesize that the NOR-5S rDNA insertion appeared in the Archelosauria clade during the Permian period and was lost later in the ancestors of Aves.

Analyses of the Updated "Animal rDNA Loci Database" with an Emphasis on Its New Features

We report on a major update to the animal rDNA loci database, which now contains cytogenetic information for 45S and 5S rDNA loci in more than 2600 and 1000 species, respectively. The data analyses show the following: (i) A high variability in 5S and 45S loci numbers, with both showing 50-fold or higher variability. However, karyotypes with an extremely high number of loci were rare, and medians generally converged to two 5S sites and two 45S rDNA sites per diploid genome. No relationship was observed between the number of 5S and 45S loci. (ii) The position of 45S rDNA on sex chromosomes was relatively frequent in some groups, particularly in arthropods (14% of karyotypes). Furthermore, 45S rDNA was almost exclusively located in microchromosomes when these were present (in birds and reptiles). (iii) The proportion of active NORs (positively stained with silver staining methods) progressively decreased with an increasing number of 45S rDNA loci, and karyotypes with more than 12 loci showed, on average, less than 40% of active loci. In conclusion, the updated version of the database provides some new insights into the organization of rRNA genes in chromosomes. We expect that its updated content will be useful for taxonomists, comparative cytogeneticists, and evolutionary biologists.

Variation in Ribosomal DNA in the Genus Trifolium (Fabaceae)

The genus Trifolium L. is characterized by basic chromosome numbers 8, 7, 6, and 5. We conducted a genus-wide study of ribosomal DNA (rDNA) structure variability in diploids and polyploids to gain insight into evolutionary history. We used fluorescent in situ hybridization to newly investigate rDNA variation by number and position in 30 Trifolium species. Evolutionary history among species was examined using 85 available sequences of internal transcribed spacer 1 (ITS1) of 35S rDNA. In diploid species with ancestral basic chromosome number (x = 8), one pair of 5S and 26S rDNA in separate or adjacent positions on a pair of chromosomes was prevalent. Genomes of species with reduced basic chromosome numbers were characterized by increased number of signals determined on one pair of chromosomes or all chromosomes. Increased number of signals was observed also in diploids Trifolium alpestre and Trifolium microcephalum and in polyploids. Sequence alignment revealed ITS1 sequences with mostly single nucleotide polymorphisms, and ITS1 diversity was greater in diploids with reduced basic chromosome numbers compared to diploids with ancestral basic chromosome number (x = 8) and polyploids. Our results suggest the presence of one 5S rDNA site and one 26S rDNA site as an ancestral state.

The relevance of chromosome fissions for major ribosomal DNA dispersion in hymenopteran insects

Ribosomal DNA (rDNA) loci are essential for cellular metabolism due to their participation in ribosome biogenesis. Although these genes have been widely cytogenetically mapped, the evolutionary mechanisms behind their variability in number and chromosomal location remain elusive, even in well-known biological groups, such as ants, bees and wasps (Insecta: Hymenoptera). To address this question in Hymenoptera and therefore advance the understanding of rDNA evolution in insects in general, we integrated molecular cytogenetic data, a phylogenomic framework, model-based predictions and genome sequencing. Hence, we assessed the main evolutionary trends shaping the chromosomal distribution of rDNA loci in Hymenoptera. We noticed the conservation of one site of rDNA per haploid genome, suggesting that a single 45S rDNA locus is the putative ancestral pattern for aculeate Hymenoptera. Moreover, our results highlighted a nonrandom distribution of rDNA in Hymenoptera karyotypes, as well as a lineage-specific preferential location. The proximal location of rDNA is favoured in species with multiple loci and in the two families of Hymenoptera that show the highest range of chromosome numbers: Formicidae and Vespidae. We propose that chromosome fissions have played a crucial role in the distribution pattern of rDNA loci through the evolutionary diversification of Hymenoptera. Moreover, our genomic analysis of two species, one with a single locus of rDNA and one with multiple loci, supported that loci multiplication is followed by sequence divergence. Our results provide detailed information about the number and chromosomal position of rDNA in Hymenoptera and, therefore, broaden our knowledge regarding rDNA evolutionary dynamics in insects.

Tissue and/or developmental stage specific methylation of nrDNA in Capsicum annuum

The nuclear ribosomal DNA (nrDNA) sequences are often used for phylogenetic analysis among organisms. Because DNA cytosine methylation and nucleolar dominancy are two common epigenetic mechanisms of nrDNA, we hypothesized that internal transcribed spacer 1 (ITS1), 5.8S rRNA and ITS2 of nrDNA sequences could be used as epigenetic biomarkers. Thus, this research was undertaken to study level and pattern of site-specific cytosine methylation of ITS1, 5.8S and ITS2 in nine tissues and/or developmental stage of pepper Capsicum annuum L. cultivar Demre Sivrisi. Tissues studied consisted of young and old roots at 30 and 90 days after sowing (das), mature dry seeds and seeds at 26 days of post anthesis (dpa), flowering buds at 1 day before flowering, pericarps at 3, 15 and 65 dpa. Levels and patterns of DNA cytosine methylation were identified at single base resolution using bisulfite conversion sequencing. Results of this study revealed that DNA cytosine level and pattern of ITS1, 5.8S and ITS2 were different in most tissues and/or developmental stages studied. In addition, methylation levels of CG, CHG and CHH contexts were also significantly different among the regions. Based on the findings of this study, it was concluded that high level of methylation of nrDNA sequences was relatively higher as observed in transposable element and promoter. On the other hand, its tissue-specific gene expression was effective as that of gene body and promoter methylation. Overall findings revealed that methylation levels of nrDNA could be used as biomarkers for tissue identification or age estimation in plants.

Karyotypic stability of Fragaria (strawberry) species revealed by cross-species chromosome painting

Chromosome karyotyping analysis is particularly useful in determining species relationships and the origin of polyploid species. Identification of individual chromosomes is the foundation for karyotype development. For Fragaria (strawberry) species, definitive identification of the individual chromosomes is extremely difficult because of their small size and similar shape. Here, we identified all chromosomes for 11 representative Fragaria species with different ploidy using a set of oligonucleotide-based probes developed in Fragaria vesca. Comprehensive molecular cytogenetic karyotypes were established based on the individually identified chromosomes. In addition, we used oligo probes to assign the 5S and 45S rDNA loci to specific chromosomes in 16 Fragaria species. We found that these Fragaria species maintained a remarkably conserved karyotype. No inter-chromosomal structural rearrangements at the cytological level were observed in any of the chromosomes among these species. Despite karyotypic stability and similarity, variations in the signal intensity of oligo probes were observed among the homologous chromosomes in several polyploid species. Moreover, most Fragaria species also showed differences in the distribution patterns of 45S and 5S rDNA. These data provide new insights into the origins of several polyploid Fragaria species.

Cytogenomics of Deschampsia P. Beauv. (Poaceae) Species Based on Sequence Analyses and FISH Mapping of CON/COM Satellite DNA Families

The genus Deschampsia P. Beauv. (Poaceae) involves a group of widespread polymorphic species, and many of them are highly tolerant to stressful environmental conditions. Genome diversity and chromosomal phylogeny within the genus are still insufficiently studied. Satellite DNAs, including CON/COM families, are the main components of the plant repeatome, which contribute to chromosome organization. For the first time, using PCR-based (Polymerase Chain Reaction) techniques and sequential BLAST (Basic Local Alignment Search Tool) and MSA (Multiple Sequence Alignment) analyses, we identified and classified CON/COM repeats in genomes of eleven Deschampsia accessions and three accessions from related genera. High homology of CON/COM sequences were revealed in the studied species though differences in single-nucleotide alteration profiles detected in homologous CON/COM regions indicated that they tended to diverge independently. The performed chromosome mapping of 45S rDNA, 5S rDNA, and CON/COM repeats in six Deschampsia species demonstrated interspecific variability in localization of these cytogenetic markers and facilitated the identification of different chromosomal rearrangements. Based on the obtained data, the studied Deschampsia species were distinguished into karyological groups, and MSA-based schematic trees were built, which could clarify the relationships within the genus. Our findings can be useful for further genetic and phylogenetic studies.

Origin and evolution of highly polymorphic rDNA sites in Alstroemeria longistaminea (Alstroemeriaceae) and related species

Alstroemeria (Alstroemeriaceae) displays a conserved and highly asymmetric karyotype, where most rDNA sites can be properly recognized by the size and morphology of the chromosomes. We analyzed the intraspecific variation of rDNA sites in A. longistaminea and compared with their distribution in other species (A. caryophyllaea and A. piauhyensis) and a representative of a sister genus, Bomarea edulis. All three species of Alstroemeria presented 2n = 16, and one to six B chromosomes were found in some individuals of A. longistaminea. There was a set of 12 conserved rDNA sites (four 5S and eight 35S) and up to 11 variable sites. B chromosomes were almost entirely covered by 35S signals, coupled with tiny 5S sites. Noteworthy, most rDNA sites found in A. caryophyllaea and A. piauhyensis were localized in chromosome positions similar to those in A. longistaminea, suggesting the existence of conserved hotspots for rDNA accumulation. Some of these hotspots were absent in Chilean Alstromeria as well in B. edulis. We propose that insertions of rDNA sequences on chromosomes do not occur randomly but rather on preferential sites or hotspots for insertions. The maintenance of these arrays, however, may be favored/constrained by different factors, resulting in stable or polymorphic sites.

Ancient Origin of Two 5S rDNA Families Dominating in the Genus Rosa and Their Behavior in the Canina-Type Meiosis

The genus Rosa comprises more than 100 woody species characterized by intensive hybridization, introgression, and an overall complex evolutionary history. Besides many diploid species (2n = 2x = 14) polyploids ranging from 3x to 10x are frequently found. Here we analyzed 5S ribosomal DNA in 19 species covering two subgenera and the major sections within subg. Rosa. In addition to diploids and polyploids with regular meiosis, we focused on 5x dogroses (Rosa sect. Caninae), which exhibit an asymmetric meiosis differentiating between bivalent- and univalent-forming chromosomes. Using genomic resources, we reconstructed 5S rDNA units to reveal their phylogenetic relationships. Additionally, we designed locus-specific probes derived from intergenic spacers (IGSs) and determined the position and number of 5S rDNA families on chromosomes. Two major 5S rDNA families (termed 5S_A and 5S_B, respectively) were found at variable ratios in both diploid and polyploid species including members of the early diverging subgenera, Rosa persica and Rosa minutifolia. Within subg. Rosa species of sect. Rosa amplified the 5S_A variant only, while taxa of other sections contained both variants at variable ratios. The 5S_B family was often co-localized with 35S rDNA at the nucleolar organizer regions (NOR) chromosomes, whereas the co-localization of the 5S_A family with NOR was only exceptionally observed. The allo-pentaploid dogroses showed a distinct distribution of 5S rDNA families between bivalent- and univalent-forming chromosomes. In conclusion, two divergent 5S rDNA families dominate rose genomes. Both gene families apparently arose in the early history of the genus, already 30 myrs ago, and apparently survived numerous speciation events thereafter. These observations are consistent with a relatively slow genome turnover in the Rosa genus.

The rDNA Loci-Intersections of Replication, Transcription, and Repair Pathways

Genes encoding ribosomal RNA (rDNA) are essential for cell survival and are particularly sensitive to factors leading to genomic instability. Their repetitive character makes them prone to inappropriate recombinational events arising from collision of transcriptional and replication machineries, resulting in unstable rDNA copy numbers. In this review, we summarize current knowledge on the structure and organization of rDNA, its role in sensing changes in the genome, and its linkage to aging. We also review recent findings on the main factors involved in chromatin assembly and DNA repair in the maintenance of rDNA stability in the model plants Arabidopsis thaliana and the moss Physcomitrella patens, providing a view across the plant evolutionary tree.