Sequencing of the Arabidopsis NOR2 reveals its distinct organization and tissue-specific rRNA ribosomal variants

Diversity of ribosomes at the level of rRNA variation associated with human health and disease

With hundreds of copies of rDNA, it is unknown whether they possess sequence variations that form different types of ribosomes. Here, we developed an algorithm for long-read variant calling, termed RGA, which revealed that variations in human rDNA loci are predominantly insertion-deletion (indel) variants. We developed full-length rRNA sequencing (RIBO-RT) and in situ sequencing (SWITCH-seq), which showed that translating ribosomes possess variation in rRNA. Over 1,000 variants are lowly expressed. However, tens of variants are abundant and form distinct rRNA subtypes with different structures near indels as revealed by long-read rRNA structure probing coupled to dimethyl sulfate sequencing. rRNA subtypes show differential expression in endoderm/ectoderm-derived tissues, and in cancer, low-abundance rRNA variants can become highly expressed. Together, this study identifies the diversity of ribosomes at the level of rRNA variants, their chromosomal location, and unique structure as well as the association of ribosome variation with tissue-specific biology and cancer.

Transcriptional Silencing of 35S rDNA in Tragopogon porrifolius Correlates with Cytosine Methylation in Sequence-Specific Manner

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.

Diversity of ribosomes at the level of rRNA variation associated with human health and disease

Ribosomal DNA and RNA (rDNA and rRNA) sequences are usually discarded from sequencing analyses. But with hundreds of copies of rDNA genes it is unknown whether they possess sequence variations that form different types of ribosomes that affect human physiology and disease. Here, we developed an algorithm for variant-calling between paralog genes (termed RGA) and compared rDNA variations found in short- and long-read sequencing data from the 1,000 Genomes Project (1KGP) and Genome In A Bottle (GIAB). We additionally developed a novel protocol for long-read sequencing full-length rRNA (RIBO-RT) from actively translating ribosomes. Our analyses identified hundreds of rDNA variants, most of which, surprisingly, are short insertion-deletions (indels) and dozens of highly abundant rRNA variants that are incorporated into translationally active ribosomes. To visualize variant ribosomes at the single cell level, we developed an in-situ rRNA sequencing method (SWITCH-seq) which revealed that variants are co-expressed within individual cells. Strikingly, by analyzing rDNA, we found that variants assemble into distinct ribosome subtypes. We discovered that these subtypes acquire different rRNA structures by successfully employing dimethyl sulfate (DMS) probing of full length rRNA. With this atlas we investigated rRNA variation changes across human tissues and cancer types. This revealed tissue-specific rRNA subtype expression in endoderm/ectoderm-derived tissues. In cancer, low abundant rRNA variants can become highly expressed, which suggests the presence of cancer-specific ribosomes. Together, this study identifies and comprehensively characterizes the diversity of ribosomes at the level of rRNA variants which is dominated by indel variants, their chromosomal location and unique structure as well as the association of ribosome variation with tissue-specific biology and cancer.

Uniparental silencing of 5S rRNA genes in plant allopolyploids - insights from Cardamine (Brassicaceae)

While the phenomenon of uniparental silencing of 35S rDNA in interspecific hybrids and allopolyploids is well documented, there is a notable absence of information regarding whether such silencing extends to the 5S RNA component of ribosomes. To address this gap in knowledge, we analyzed the 5S and 35S rDNA expression in Cardamine (Brassicaceae) allopolyploids, namely C. × insueta (2n = 3x = 24, genome composition RRA), C. flexuosa (2n = 4x = 32, AAHH), and C. scutata (2n = 4x = 32, PPAA) which share a common diploid ancestor (AA). We employed high-throughput sequencing of transcriptomes and genomes and phylogenetic analyses of 5S rRNA variants. The genomic organization of rDNA was further scrutinized through clustering and fluorescence in situ hybridization. In the C. × insueta allotriploid, we observed uniparental dominant expression of 5S and 35S rDNA loci. In the C. flexuosa and C. scutata allotetraploids, the expression pattern differed, with the 35S rDNA being expressed from the A subgenome, whereas the 5S rDNA was expressed from the partner subgenome. Both C. flexuosa and C. scutata but not C. × insueta showed copy and locus number changes. We conclude that in stabilized allopolyploids, transcription of ribosomal RNA components occurs from different subgenomes. This phenomenon appears to result in the formation of chimeric ribosomes comprising rRNA molecules derived from distinct parental origins. We speculate that the interplay of epigenetic silencing and rDNA rearrangements introduces an additional layer of variation in multimolecule ribosomal complexes, potentially contributing to the evolutionary success of allopolyploids.

The interplay of homology-directed repair pathways in the repair of zebularine-induced DNA-protein crosslinks in Arabidopsis

DNA-protein crosslinks (DPCs) are highly toxic DNA lesions represented by proteins covalently bound to the DNA. Persisting DPCs interfere with fundamental genetic processes such as DNA replication and transcription. Cytidine analog zebularine (ZEB) has been shown to crosslink DNA METHYLTRANSFERASE1 (MET1). Recently, we uncovered a critical role of the SMC5/6-mediated SUMOylation in the repair of DPCs. In an ongoing genetic screen, we identified two additional candidates, HYPERSENSITIVE TO ZEBULARINE 2 and 3, that were mapped to REGULATOR OF TELOMERE ELONGATION 1 (RTEL1) and polymerase TEBICHI (TEB), respectively. By monitoring the growth of hze2 and hze3 plants in response to zebularine, we show the importance of homologous recombination (HR) factor RTEL1 and microhomology-mediated end-joining (MMEJ) polymerase TEB in the repair of MET1-DPCs. Moreover, genetic interaction and sensitivity assays showed the interdependency of SMC5/6 complex, HR, and MMEJ in the homology-directed repair of MET1-DPCs in Arabidopsis. Altogether, we provide evidence that MET1-DPC repair in plants is more complex than originally expected.

The MBD-ACD DNA methylation reader complex recruits MICRORCHIDIA6 to regulate ribosomal RNA gene expression in Arabidopsis

DNA methylation is an important epigenetic mark implicated in selective rRNA gene expression, but the DNA methylation readers and effectors remain largely unknown. Here, we report a protein complex that reads DNA methylation to regulate variant-specific 45S ribosomal RNA (rRNA) gene expression in Arabidopsis (Arabidopsis thaliana). The complex, consisting of METHYL-CpG-BINDING DOMAIN PROTEIN5 (MBD5), MBD6, ALPHA-CRYSTALLIN DOMAIN PROTEIN15.5 (ACD15.5), and ACD21.4, directly binds to 45S rDNA. While MBD5 and MBD6 function redundantly, ACD15.5 and ACD21.4 are indispensable for variant-specific rRNA gene expression. These 4 proteins undergo phase separation in vitro and in vivo and are interdependent for their phase separation. The α-crystallin domain of ACD15.5 and ACD21.4, which is essential for their function, enables phase separation of the complex, likely by mediating multivalent protein interactions. The effector MICRORCHIDIA6 directly interacts with ACD15.5 and ACD21.4, but not with MBD5 and MBD6, and is recruited to 45S rDNA by the MBD-ACD complex to regulate variant-specific 45S rRNA expression. Our study reveals a pathway in Arabidopsis through which certain 45S rRNA gene variants are silenced, while others are activated.

Low dose ribosomal DNA P-loop mutation affects development and enforces autophagy in Arabidopsis

Arabidopsis contains hundreds of ribosomal DNA copies organized within the nucleolar organizing regions (NORs) in chromosomes 2 and 4. There are four major types of variants of rDNA, VAR1-4, based on the polymorphisms of 3' external transcribed sequences. The variants are known to be differentially expressed during plant development. We created a mutant by the CRISPR-Cas9-mediated excision of ~ 25 nt from predominantly NOR4 ribosomal DNA copies, obtaining mosaic mutational events on ~ 5% of all rDNA copies. The excised region consists of P-loop and Helix-82 segments of 25S rRNA. The mutation led to allelic, dosage-dependent defects marked by lateral root inhibition, reduced size, and pointy leaves, all previously observed for defective ribosomal function. The mutation in NOR4 led to dosage compensation from the NOR2 copies by elevated expression of VAR1 in mutants and further associated single-nucleotide variants, thus, resulting in altered rRNA sub-population. Furthermore, the mutants exhibited rRNA maturation defects specifically in the minor pathway typified by 32S pre-rRNA accumulation. Density-gradient fractionation and subsequent RT-PCR of rRNA analyses revealed that mutated copies were not incorporated into the translating ribosomes. The mutants in addition displayed an elevated autophagic flux as shown by the autophagic marker GFP-ATG8e, likely related to ribophagy.

Sequence and epigenetic landscapes of active and silent nucleolus organizer regions in Arabidopsis

Arabidopsis thaliana has two ribosomal RNA (rRNA) gene loci, nucleolus organizer regions NOR2 and NOR4, whose complete sequences are missing in current genome assemblies. Ultralong DNA sequences assembled using an unconventional approach yielded ~5.5- and 3.9-Mbp sequences for NOR2 and NOR4 in the reference strain, Col-0. The distinct rRNA gene subtype compositions of the NORs enabled the positional mapping of their active and inactive regions, using RNA sequencing to identify subtype-specific transcripts and DNA sequencing to identify subtypes associated with flow-sorted nucleoli. Comparisons of wild-type and silencing-defective plants revealed that most rRNA gene activity occurs in the central region of NOR4, whereas most, but not all, genes of NOR2 are epigenetically silenced. Intervals of low CG and CHG methylation overlap regions where gene activity and gene subtype homogenization are high. Collectively, the data reveal the genetic and epigenetic landscapes underlying nucleolar dominance (differential NOR activity) and implicate transcription as a driver of rRNA gene concerted evolution.

New insights on the evolution of nucleolar dominance in newly resynthesized hexaploid wheat Triticum zhukovskyi

Nucleolar dominance (ND) is a widespread epigenetic phenomenon in hybridizations where nucleolus transcription fails at the nucleolus organizer region (NOR). However, the dynamics of NORs during the formation of Triticum zhukovskyi (GGAu Au Am Am ), another evolutionary branch of allohexaploid wheat, remains poorly understood. Here, we elucidated genetic and epigenetic changes occurring at the NOR loci within the Am , G, and D subgenomes during allopolyploidization by synthesizing hexaploid wheat GGAu Au Am Am and GGAu Au DD. In T. zhukovskyi, Au genome NORs from T. timopheevii (GGAu Au ) were lost, while the second incoming NORs from T. monococcum (Am Am ) were retained. Analysis of the synthesized T. zhukovskyi revealed that rRNA genes from the Am genome were silenced in F1 hybrids (GAu Am ) and remained inactive after genome doubling and subsequent self-pollinations. We observed increased DNA methylation accompanying the inactivation of NORs in the Am genome and found that silencing of NORs in the S1 generation could be reversed by a cytidine methylase inhibitor. Our findings provide insights into the ND process during the evolutionary period of T. zhukovskyi and highlight that inactive rDNA units may serve as a 'first reserve' in the form of R-loops, contributing to the successful evolution of T. zhukovskyi.

45S rDNA Diversity In Natura as One Step towards Ribosomal Heterogeneity in Arabidopsis thaliana

The keystone of ribosome biogenesis is the transcription of 45S rDNA. The Arabidopsis thaliana genome contains hundreds of 45S rDNA units; however, they are not all transcribed. Notably, 45S rDNA units contain insertions/deletions revealing the existence of heterogeneous rRNA genes and, likely, heterogeneous ribosomes for rRNAs. In order to obtain an overall picture of 45S rDNA diversity sustaining the synthesis of rRNAs and, subsequently, of ribosomes in natura, we took advantage of 320 new occurrences of Arabidopsis thaliana as a metapopulation named At66, sampled from 0 to 1900 m of altitude in the eastern Pyrenees in France. We found that the 45S rDNA copy number is very dynamic in natura and identified new genotypes for both 5' and 3' External Transcribed Spacers (ETS). Interestingly, the highest 5'ETS genotype diversity is found in altitude while the highest 3'ETS genotype diversity is found at sea level. Structural analysis of 45S rDNA also shows conservation in natura of specific 5'ETS and 3'ETS sequences/features required to control rDNA expression and the processing of rRNAs. In conclusion, At66 is a worthwhile natural laboratory, and unraveled 45S rDNA diversity represents an interesting starting material to select subsets for rDNA transcription and alter the rRNA composition of ribosomes both intra- and inter-site.

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.

Switch them off or not: selective rRNA gene repression in grasses

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

Structural variation among assembled genomes facilitates development of rapid and low-cost NOR-linked markers and NOR-telomere junction mapping in Arabidopsis

KEY MESSAGE

Genome-wide structural variants we identified and new NOR-linked markers we developed would be useful for future genome-wide association studies (GWAS), and for new gene/trait mapping purposes. Bioinformatic alignment of the assembled genomes of Col-0 and Sha ecotypes of Arabidopsis thaliana revealed ~ 13,000 genome-wide structural variants involving simple insertions or deletions and repeat contractions or expansions. Using some of these structural variants, we developed new, rapid, and low-cost PCR-based molecular markers that are genetically linked to the nucleolus organizer regions (NORs). A. thaliana has two NORs, one each on chromosome 2 (NOR2) and chromosome 4 (NOR4). Both NORs are ~ 4 Mb each, and hundreds of 45S ribosomal RNA (rRNA) genes are tandemly arrayed at these loci. Using previously characterized recombinant inbred lines (RILs) derived from Sha x Col-0  crosses, we validated the utility of the newly developed NOR-linked markers in genetically mapping rRNA genes and the associated telomeres to either NOR2 or NOR4. Lastly, we sequenced Sha genome using Oxford Nanopore Technology (ONT) and used the data to obtain sequences of NOR-telomere junctions, and with the help of RILs, we mapped them as new genetic markers to their respective NORs (NOR2-TEL2N and NOR4-TEL4N). The structural variants obtained from this study would serve as valuable data for genome-wide association studies (GWAS), and to rapidly design more genome-wide genetic (molecular) markers for new gene/trait mapping purposes.

Under the magnifying glass: The ups and downs of rDNA copy number

The ribosomal DNA (rDNA) in Drosophila is found as two additive clusters of individual 35 S cistrons. The multiplicity of rDNA is essential to assure proper translational demands, but the nature of the tandem arrays expose them to copy number variation within and between populations. Here, we discuss means by which a cell responds to insufficient rDNA copy number, including a historical view of rDNA magnification whose mechanism was inferred some 35 years ago. Recent work has revealed that multiple conditions may also result in rDNA loss, in response to which rDNA magnification may have evolved. We discuss potential models for the mechanism of magnification, and evaluate possible consequences of rDNA copy number variation.

Varying strength of selection contributes to the intragenomic diversity of rRNA genes

Ribosome biogenesis in eukaryotes is supported by hundreds of ribosomal RNA (rRNA) gene copies that are encoded in the ribosomal DNA (rDNA). The multiple copies of rRNA genes are thought to have low sequence diversity within one species. Here, we present species-wide rDNA sequence analysis in Saccharomyces cerevisiae that challenges this view. We show that rDNA copies in this yeast are heterogeneous, both among and within isolates, and that many variants avoided fixation or elimination over evolutionary time. The sequence diversity landscape across the rDNA shows clear functional stratification, suggesting different copy-number thresholds for selection that contribute to rDNA diversity. Notably, nucleotide variants in the most conserved rDNA regions are sufficiently deleterious to exhibit signatures of purifying selection even when present in only a small fraction of rRNA gene copies. Our results portray a complex evolutionary landscape that shapes rDNA sequence diversity within a single species and reveal unexpectedly strong purifying selection of multi-copy genes.

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.

Upon heat stress processing of ribosomal RNA precursors into mature rRNAs is compromised after cleavage at primary P site in Arabidopsis thaliana

Transcription and processing of 45S rRNAs in the nucleolus are keystones of ribosome biogenesis. While these processes are severely impacted by stress conditions in multiple species, primarily upon heat exposure, we lack information about the molecular mechanisms allowing sessile organisms without a temperature-control system, like plants, to cope with such circumstances. We show that heat stress disturbs nucleolar structure, inhibits pre-rRNA processing and provokes imbalanced ribosome profiles in Arabidopsis thaliana plants. Notably, the accuracy of transcription initiation and cleavage at the primary P site in the 5’ETS (5’ External Transcribed Spacer) are not affected but the levels of primary 45S and 35S transcripts are, respectively, increased and reduced. In contrast, precursors of 18S, 5.8S and 25S RNAs are rapidly undetectable upon heat stress. Remarkably, nucleolar structure, pre-rRNAs from major ITS1 processing pathway and ribosome profiles are restored after returning to optimal conditions, shedding light on the extreme plasticity of nucleolar functions in plant cells. Further genetic and molecular analysis to identify molecular clues implicated in these nucleolar responses indicate that cleavage rate at P site and nucleolin protein expression can act as a checkpoint control towards a productive pre-rRNA processing pathway.

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.

Fine structure and transcription dynamics of bread wheat ribosomal DNA loci deciphered by a multi-omics approach

Three out of four RNA components of ribosomes are encoded by 45S ribosomal DNA (rDNA) loci, which are organized as long head-to-tail tandem arrays of nearly identical units, spanning several megabases of sequence. Due to this structure, the rDNA loci are the major sources of gaps in genome assemblies, and gene copy number, sequence composition, and expression status of particular arrays remain elusive, especially in complex genomes harboring multiple loci. Here we conducted a multi-omics study to decipher the 45S rDNA loci in hexaploid bread wheat. Coupling chromosomal genomics with optical mapping, we reconstructed individual rDNA arrays, enabling locus-specific analyses of transcription activity and methylation status from RNA- and bisulfite-sequencing data. We estimated a total of 6,650 rDNA units in the bread wheat genome, with approximately 2,321, 3,910, 253, and 50 gene copies located in short arms of chromosomes 1B, 6B, 5D, and 1A, respectively. Only 1B and 6B loci contributed substantially to rRNA transcription at a roughly 2:1 ratio. The ratio varied among five tissues analyzed (embryo, coleoptile, root tip, primary leaf, mature leaf), being the highest (2.64:1) in mature leaf and lowest (1.72:1) in coleoptile. Cytosine methylation was considerably higher in CHG context in the silenced 5D locus as compared with the active 1B and 6B loci. In conclusion, a fine genomic organization and tissue-specific expression of rDNA loci were deciphered, for the first time, in a complex polyploid species. The results are discussed in the context of wheat evolution and transcription regulation.

Nucleolar Organizer Regions as Transcription-Based Scaffolds of Nucleolar Structure and Function

Eukaryotic genomes maintain multiple copies of ribosomal DNA gene repeats in tandem arrays to provide sufficient ribosomal RNAs to make ribosomes. These DNA repeats are the most highly transcribed regions of the genome, with dedicated transcriptional machinery to manage the enormous task of producing more than 50% of the total RNA in a proliferating cell. The arrays are called nucleolar organizer regions (NORs) and constitute the scaffold of the nucleolar compartment, where ribosome biogenesis occurs. Advances in molecular and cellular biology have brought great insights into how these arrays are transcribed and organized within genomes. Much of their biology is driven by their high transcription level, which has also driven the development of unique methods to understand rDNA gene activity, beginning with classic techniques such as silver staining and Miller spreads. However, the application of modern methodologies such as CRISPR gene editing, super-resolution microscopy, and long-read sequencing has enabled recent advances described herein, with many more discoveries possible soon. This chapter highlights what is known about NOR transcription and organization and the techniques applied historically and currently. Given the potential for NORs to impact organismal health and disease, as highlighted at the end of the chapter, the field must continue to develop and apply innovative analysis to understand genetic, epigenetic, and organizer properties of the ribosomal DNA repeats.

Characterization of the Different Levels of Variation in 45S rRNA Genes

Although commonly regarded as nearly identical in sequence, 45S ribosomal RNA genes represent a massive source of genetic variation at different levels. Unfortunately, due to their repetitive nature and the difficulty to assemble their multiple copies in tandem, these important genomic elements remain largely unexplored in sequencing projects. Here, I describe how to exploit next generation sequencing data to estimate their copy number in an organism and detect true polymorphic sites within and among individuals. Furthermore, for species that carry multiple 45S ribosomal RNA gene clusters, I show how to make use of experimental populations to assign some of these variants to their cluster of origin.

Integrating Wheat Nucleolus Structure and Function: Variation in the Wheat Ribosomal RNA and Protein Genes

We review the coordinated production and integration of the RNA (ribosomal RNA, rRNA) and protein (ribosomal protein, RP) components of wheat cytoplasmic ribosomes in response to changes in genetic constitution, biotic and abiotic stresses. The components examined are highly conserved and identified with reference to model systems such as human, Arabidopsis, and rice, but have sufficient levels of differences in their DNA and amino acid sequences to form fingerprints or gene haplotypes that provide new markers to associate with phenotype variation. Specifically, it is argued that populations of ribosomes within a cell can comprise distinct complements of rRNA and RPs to form units with unique functionalities. The unique functionalities of ribosome populations within a cell can become central in situations of stress where they may preferentially translate mRNAs coding for proteins better suited to contributing to survival of the cell. In model systems where this concept has been developed, the engagement of initiation factors and elongation factors to account for variation in the translation machinery of the cell in response to stresses provided the precedents. The polyploid nature of wheat adds extra variation at each step of the synthesis and assembly of the rRNAs and RPs which can, as a result, potentially enhance its response to changing environments and disease threats.

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

Licorice Germplasm Resources Identification Using DNA Barcodes Inner-Variants

Based on the gradual transformation from wild growth to artificial cultivation, the accurate authentication of licorice seeds contributes to the first committed step of its quality control and is pivotal to ensure the clinical efficacy of licorice. However, it is still challenging to obtain genetically stable licorice germplasm resources due to the multi-source, multi-heterozygous, polyploid, and hybrid characteristics of licorice seeds. Here, a new method for determining the heterozygosity of licorice seed mixture, based on the various sites, and finding the composition characteristics of licorice seed is preliminarily designed and proposed. Namely, high-throughput full-length multiple DNA barcodes(HFMD), based on ITS multi-copy variation exist, the full-length amplicons of ITS2, psbA-trnH and ITS are directly sequenced by rDNA through the next-generation sequence(NGS) and single-molecule real-time (SMRT) technologies. By comparing the three sequencing methods, our results proved that SMRT sequencing successfully identified the complete gradients of complex mixed samples with the best performance. Meanwhile, HFMD is a brilliant and feasible method for evaluating the heterozygosity of licorice seeds. It shows a perfect interpretation of DNA barcoding and can be applied in multi-base multi-heterozygous and polyploid species.

Mosaic Arrangement of the 5S rDNA in the Aquatic Plant Landoltia punctata (Lemnaceae)

Duckweeds are a group of monocotyledonous aquatic plants in the Araceae superfamily, represented by 37 species divided into five genera. Duckweeds are the fastest growing flowering plants and are distributed around the globe; moreover, these plants have multiple applications, including biomass production, wastewater remediation, and making pharmaceutical proteins. Dotted duckweed (Landoltia punctata), the sole species in genus Landoltia, is one of the most resilient duckweed species. The ribosomal DNA (rDNA) encodes the RNA components of ribosomes and represents a significant part of plant genomes but has not been comprehensively studied in duckweeds. Here, we characterized the 5S rDNA genes in L. punctata by cloning and sequencing 25 PCR fragments containing the 5S rDNA repeats. No length variation was detected in the 5S rDNA gene sequence, whereas the nontranscribed spacer (NTS) varied from 151 to 524 bp. The NTS variants were grouped into two major classes, which differed both in nucleotide sequence and the type and arrangement of the spacer subrepeats. The dominant class I NTS, with a characteristic 12-bp TC-rich sequence present in 3-18 copies, was classified into four subclasses, whereas the minor class II NTS, with shorter, 9-bp nucleotide repeats, was represented by two identical sequences. In addition to these diverse subrepeats, class I and class II NTSs differed in their representation of cis-elements and the patterns of predicted G-quadruplex structures, which may influence the transcription of the 5S rDNA. Similar to related duckweed species in the genus Spirodela, L. punctata has a relatively low rDNA copy number, but in contrast to Spirodela and the majority of other plants, the arrangement of the 5S rDNA units demonstrated an unusual, heterogeneous pattern in L. punctata, as revealed by analyzing clones containing double 5S rDNA neighboring units. Our findings may further stimulate the research on the evolution of the plant rDNA and discussion of the molecular forces driving homogenization of rDNA repeats in concerted evolution.

From Microscopy to Nanoscopy: Defining an Arabidopsis thaliana Meiotic Atlas at the Nanometer Scale

Visualization of meiotic chromosomes and the proteins involved in meiotic recombination have become essential to study meiosis in many systems including the model plant Arabidopsis thaliana. Recent advances in super-resolution technologies changed how microscopic images are acquired and analyzed. New technologies enable observation of cells and nuclei at a nanometer scale and hold great promise to the field since they allow observing complex meiotic molecular processes with unprecedented detail. Here, we provide an overview of classical and advanced sample preparation and microscopy techniques with an updated Arabidopsis meiotic atlas based on super-resolution microscopy. We review different techniques, focusing on stimulated emission depletion (STED) nanoscopy, to offer researchers guidance for selecting the optimal protocol and equipment to address their scientific question.

Targeted Enrichment of rRNA Gene Tandem Arrays for Ultra-Long Sequencing by Selective Restriction Endonuclease Digestion

Large regions of nearly identical repeats, such as the 45S ribosomal RNA (rRNA) genes of Nucleolus Organizer Regions (NORs), can account for major gaps in sequenced genomes. To assemble these regions, ultra-long sequencing reads that span multiple repeats have the potential to reveal sets of repeats that collectively have sufficient sequence variation to unambiguously define that interval and recognize overlapping reads. Because individual repetitive loci typically represent a small proportion of the genome, methods to enrich for the regions of interest are desirable. Here we describe a simple method that achieves greater than tenfold enrichment of Arabidopsis thaliana 45S rRNA gene sequences among ultra-long Oxford Nanopore Technology sequencing reads. This method employs agarose-embedded genomic DNA that is subjected to restriction endonucleases digestion using a cocktail of enzymes predicted to be non-cutters of rRNA genes. Most of the genome is digested into small fragments that diffuse out of the agar plugs, whereas rRNA gene arrays are retained. In principle, the approach can also be adapted for sequencing other repetitive loci for which gaps exist in a reference genome.

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.

Sequencing of the Arabidopsis NOR2 reveals its distinct organization and tissue-specific rRNA ribosomal variants

Despite vast differences between organisms, some characteristics of their genomes are conserved, such as the nucleolus organizing region (NOR). The NOR is constituted of multiple, highly repetitive rDNA genes, encoding the catalytic ribosomal core RNAs which are transcribed from 45S rDNA units. Their precise sequence information and organization remain uncharacterized. Here, using a combination of long- and short-read sequencing technologies we assemble contigs of the Arabidopsis NOR2 rDNA domain. We identify several expressed rRNA gene variants which are integrated into translating ribosomes in a tissue-specific manner. These findings support the concept of tissue specific ribosome subpopulations that differ in their rRNA composition and provide insights into the higher order organization of NOR2.

The nucleolus organizing region (NOR) consists of multiple, highly repetitive rDNA genes. Here Sims et al. use both long- and short-read sequencing to determine the organization and sequence of Arabidopsis NOR2 rDNA and show that different rRNA gene variants are integrated into translating ribosomes in a tissue-specific manner.

It Is Just a Matter of Time: Balancing Homologous Recombination and Non-homologous End Joining at the rDNA Locus During Meiosis

Ribosomal RNA genes (rDNAs) are located in large domains of hundreds of rDNA units organized in a head-to-tail manner. The proper and stable inheritance of rDNA clusters is of paramount importance for survival. Yet, these highly repetitive elements pose a potential risk to the genome since they can undergo non-allelic exchanges. Here, we review the current knowledge of the organization of the rDNA clusters in Arabidopsis thaliana and their stability during meiosis. Recent findings suggest that during meiosis, all rDNA loci are embedded within the nucleolus favoring non-homologous end joining (NHEJ) as a repair mechanism, while DNA repair via homologous recombination (HR) appears to be a rare event. We propose a model where (1) frequent meiotic NHEJ events generate abundant single nucleotide polymorphisms and insertions/deletions within the rDNA, resulting in a heterogeneous population of rDNA units and (2) rare HR events dynamically change rDNA unit numbers, only to be observed in large populations over many generations. Based on the latest efforts to delineate the entire rDNA sequence in A. thaliana, we discuss evidence supporting this model. The results compiled so far draw a surprising picture of rDNA sequence heterogeneity between individual units. Furthermore, rDNA cluster sizes have been recognized as relatively stable when observing less than 10 generations, yet emerged as major determinant of genome size variation between different A. thaliana ecotypes. The sequencing efforts also revealed that transcripts from the diverse rDNA units yield heterogenous ribosome populations with potential functional implications. These findings strongly motivate further research to understand the mechanisms that maintain the metastable state of rDNA loci.