Widespread genetic heterogeneity of human ribosomal RNA genes

Assessing Human Ribosomal DNA Variation and Its Association With Phenotypic Outcomes

Although genome-scale analyses have provided insights into the connection between genetic variability and complex human phenotypes, much trait variation is still not fully understood. Genetic variation within repetitive elements, such as the multi-copy, multi-locus ribosomal DNA (rDNA), has emerged as a potential contributor to trait variation. Whereas rDNA was long believed to be largely uniform within a species, recent studies have revealed substantial variability in the locus, both within and across individuals. This variation, which takes the form of copy number, structural arrangement, and sequence differences, has been found to be associated with human phenotypes. This review summarizes what is currently known about human rDNA variation, its causes, and its association with phenotypic outcomes, highlighting the technical challenges the field faces and the solutions proposed to address them. Finally, we suggest experimental approaches that can help clarify the elusive mechanisms underlying the phenotypic consequences of rDNA variation.

Small nucleolar RNAs: the hidden precursors of cancer ribosomes

Ribosomes are heterogeneous in terms of their constituent proteins, structural RNAs and ribosomal RNA (rRNA) modifications, resulting in diverse potential translatomes. rRNA modifications, guided by small nucleolar RNAs (snoRNAs), enable fine-tuning of ribosome function and translation profiles. Recent studies have begun linking dysregulation of snoRNAs, via rRNA modifications, to tumourigenesis. Deciphering the specific contributions of individual rRNA modifications to cancer hallmarks and identifying snoRNAs with oncogenic potential could lead to novel therapeutic strategies. These strategies might target snoRNAs or exploit the dependence of cancer cells on specific rRNA modification sites, potentially disrupting aberrant ribosomal translation programs and hindering tumour growth. This review discusses current evidence and challenges in linking changes in snoRNA expression to rRNA modification and cancer biology.This article is part of the discussion meeting issue 'Ribosome diversity and its impact on protein synthesis, development and disease'.

Impacts of ribosomal RNA sequence variation on gene expression and phenotype

Since the framing of the Central Dogma, it has been speculated that physically distinct ribosomes within cells may influence gene expression and cellular physiology. While heterogeneity in ribosome composition has been reported in bacteria, protozoans, fungi, zebrafish, mice and humans, its functional implications remain actively debated. Here, we review recent evidence demonstrating that expression of conserved variant ribosomal DNA (rDNA) alleles in bacteria, mice and humans renders their actively translating ribosome pool intrinsically heterogeneous at the level of ribosomal RNA (rRNA). In this context, we discuss reports that nutrient limitation-induced stress in Escherichia coli leads to changes in variant rRNA allele expression, programmatically altering transcription and cellular phenotype. We highlight that cells expressing ribosomes from distinct operons exhibit distinct drug sensitivities, which can be recapitulated in vitro and potentially rationalized by subtle perturbations in ribosome structure or in their dynamic properties. Finally, we discuss evidence that differential expression of variant rDNA alleles results in different populations of ribosome subtypes within mammalian tissues. These findings motivate further research into the impacts of rRNA heterogeneities on ribosomal function and predict that strategies targeting distinct ribosome subtypes may hold therapeutic potential.This article is part of the discussion meeting issue 'Ribosome diversity and its impact on protein synthesis, development and disease'.

High-throughput approaches for the identification of ribosome heterogeneity

Recent advances in the fields of RNA translation and ribosome biology have demonstrated the heterogeneous nature of ribosomes. This manifests not only across different cellular contexts but also within the same cell. Such variations in ribosomal composition, be it in ribosomal RNAs or proteins, can significantly influence cellular processes and responses by altering the mRNAs being translated or the dynamics of ribosomes during the translation process. Therefore, identifying this heterogeneity is crucial for unravelling the complexity of gene expression across different fields of biology. Here we provide an overview of recent advances in high-throughput techniques for identifying ribosomal heterogeneity. We cover methodologies for probing both rRNA and protein components of the ribosome and encompass the most recent next-generation sequencing and computational analyses, as well as a diverse array of mass spectrometry techniques.This article is part of the discussion meeting issue 'Ribosome diversity and its impact on protein synthesis, development and disease'.

Germline sequence variation within the ribosomal DNA is associated with human complex traits

The ribosome is one of the core macromolecules in the cell. The ribosomal RNAs (rRNA), which are essential components of the ribosome, are coded by the multi-copy ribosomal DNA (rDNA). Despite its highly conserved function, the rDNA displays substantial variation within all species analysed to date. This variation comprises both inter-individual differences in total copy number (CN) as well as inter- and intragenomic sequence variation in the form of single nucleotide variants (SNV) and insertions/deletions (INDELs) across rDNA copies. Whether germline variation of rDNA sequence associates with phenotypic traits in humans is, to date, unknown. Here, using the UK Biobank whole genome sequencing data, we first derive a high confidence list of rDNA-associated SNVs and INDELs that we validate in multiple ways. Using this list, we show that specific rDNA variants associate with several human traits. In particular, traits associated with body size appear enriched in variants within the Expansion Segment 15L region in the 28S rRNA. The strength of these associations does not diminish when accounting for the total rDNA CN of each individual. Our work represents the first large-scale association analysis of human traits with germline sequence variation in the rDNA, a source of human complex trait-relevant genetic variation that has thus far been largely ignored.

High-resolution reconstruction of a C. elegans ribosome sheds light on evolutionary dynamics and tissue specificity

Caenorhabditis elegans is an important model organism for human health and disease, with foundational contributions to the understanding of gene expression and tissue patterning in animals. An invaluable tool in modern gene expression research is the presence of a high-resolution ribosome structure, though no such structure exists for C. elegans Here, we present a high-resolution single-particle cryogenic electron microscopy (cryo-EM) reconstruction and molecular model of a C. elegans ribosome, revealing a significantly streamlined animal ribosome. Many facets of ribosome structure are conserved in C. elegans, including overall ribosomal architecture and the mechanism of cycloheximide, whereas other facets, such as expansion segments and eL28, are rapidly evolving. We identify uL5 and uL23 as two instances of tissue-specific ribosomal protein paralog expression conserved in Caenorhabditis, suggesting that C. elegans ribosomes vary across tissues. The C. elegans ribosome structure will provide a basis for future structural, biochemical, and genetic studies of translation in this important animal system.

Epigenetic control and inheritance of rDNA arrays

Ribosomal RNA (rRNA) genes exist in multiple copies arranged in tandem arrays known as ribosomal DNA (rDNA). The total number of gene copies is variable, and the mechanisms buffering this copy number variation remain unresolved. We surveyed the number, distribution, and activity of rDNA arrays at the level of individual chromosomes across multiple human and primate genomes. Each individual possessed a unique fingerprint of copy number distribution and activity of rDNA arrays. In some cases, entire rDNA arrays were transcriptionally silent. Silent rDNA arrays showed reduced association with the nucleolus and decreased interchromosomal interactions, indicating that the nucleolar organizer function of rDNA depends on transcriptional activity. Methyl-sequencing of flow-sorted chromosomes, combined with long read sequencing, showed epigenetic modification of rDNA promoter and coding region by DNA methylation. Silent arrays were in a closed chromatin state, as indicated by the accessibility profiles derived from Fiber-seq. Removing DNA methylation restored the transcriptional activity of silent arrays. Array activity status remained stable through the iPS cell re-programming. Family trio analysis demonstrated that the inactive rDNA haplotype can be traced to one of the parental genomes, suggesting that the epigenetic state of rDNA arrays may be heritable. We propose that the dosage of rRNA genes is epigenetically regulated by DNA methylation, and these methylation patterns specify nucleolar organizer function and can propagate transgenerationally.

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.

Ribosomes unraveled: The path from variant to impact

In this issue of Cell Genomics, Rothschild et al.1 reveal how ribosomal RNA diversity impacts ribosome structure and its implications for health and disease. Their innovative methodologies uncover distinct ribosome subtypes with significant structural variations and expression patterns. This work reveals connections to tissue-specific biology and cancer, positing new research avenues.

Origin and maintenance of large ribosomal RNA gene repeat size in mammals

The genes encoding ribosomal RNA are highly conserved across life and in almost all eukaryotes are present in large tandem repeat arrays called the rDNA. rDNA repeat unit size is conserved across most eukaryotes but has expanded dramatically in mammals, principally through the expansion of the intergenic spacer region that separates adjacent rRNA coding regions. Here, we used long-read sequence data from representatives of the major amniote lineages to determine where in amniote evolution rDNA unit size increased. We find that amniote rDNA unit sizes fall into two narrow size classes: "normal" (∼11-20 kb) in all amniotes except monotreme, marsupial, and eutherian mammals, which have "large" (∼35-45 kb) sizes. We confirm that increases in intergenic spacer length explain much of this mammalian size increase. However, in stark contrast to the uniformity of mammalian rDNA unit size, mammalian intergenic spacers differ greatly in sequence. These results suggest a large increase in intergenic spacer size occurred in a mammalian ancestor and has been maintained despite substantial sequence changes over the course of mammalian evolution. This points to a previously unrecognized constraint on the length of the intergenic spacer, a region that was thought to be largely neutral. We finish by speculating on possible causes of this constraint.

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.

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.

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.

Construction and validation of customized genomes for human and mouse ribosomal DNA mapping

rRNAs are transcribed from ribosomal DNA (rDNA) repeats, the most intensively transcribed loci in the genome. Due to their repetitive nature, there is a lack of genome assemblies suitable for rDNA mapping, creating a vacuum in our understanding of how the most abundant RNA in the cell is regulated. Our recent work revealed binding of numerous mammalian transcription and chromatin factors to rDNA. Several of these factors were known to play critical roles in development, tissue function, and malignancy, but their potential roles in rDNA regulation remained unexplored. This demonstrated the blind spot into which rDNA has fallen in genetic and epigenetic studies and highlighted an unmet need for public rDNA-optimized genome assemblies. Here, we customized five human and mouse assemblies-hg19 (GRCh37), hg38 (GRCh38), hs1 (T2T-CHM13), mm10 (GRCm38), and mm39 (GRCm39)-to render them suitable for rDNA mapping. The standard builds of these genomes contain numerous fragmented or repetitive rDNA loci. We identified and masked all rDNA-like regions, added a single rDNA reference sequence of the appropriate species as a ∼45 kb chromosome designated "chromosome R," and created annotation files to aid visualization of rDNA features in browser tracks. We validated these customized genomes for mapping of known rDNA-binding proteins and present a simple workflow for mapping chromatin immunoprecipitation-sequencing datasets. Customized genome assemblies, annotation files, positive and negative control tracks, and Snapgene files of standard rDNA reference sequences have been deposited to GitHub. These resources make rDNA mapping and visualization more readily accessible to a broad audience.

Hexavalent chromium [Cr(VI)]-induced ribosomal DNA copy number variation and DNA damage responses and their associations with nucleolar protein HRAS in humans and cells

The carcinogenicity of hexavalent chromium [Cr(VI)] and its compounds has been widely recognized, yet the mechanism of genetic damage is still not fully understood. The ribosomal DNA (rDNA) copy number is recently considered a potential marker of cancer-associated stress. To investigate the roles of rDNA copy number variation (CNV) in DNA damage responses (DDRs) induced by Cr(VI) and the potential mechanism from nucleolar protein HRAS, a cross-sectional study in Cr(Ⅵ)-exposed workers and an in vitro experiment using HeLa cells were conducted. Our results showed increased levels of rDNA CNV, DDRs, and HRAS expression in Cr(VI)-exposed workers. Generalized linear regression analyses showed that Cr(VI) exposure was significantly positively associated with increased levels of rDNA CNV, DDRs, and HRAS expression in Cr(VI)-exposed workers. Moreover, there were pairwise associations between rDNA CNV, DDRs, and HRAS levels. Mediation analyses found that rDNA CNV significantly mediated the association between Cr(VI) exposure and DDRs. The in vitro experiments further confirmed that Cr(VI) treatment induced increased levels of rDNA CNV, DDRs, and HRAS expression in HeLa cells. Cr(VI)-induced rDNA CNV, ATM activation, and apoptosis damage were then strongly enhanced by HRAS depletion with siRNA in vitro, suggesting the important role of HRAS in CNV and DDRs caused by Cr(VI). The combined results of the human and cell line studies indicated that Cr(VI) exposure might enhance rDNA CNV by regulation of HRAS expression, which leads to Cr(VI)-induced genetic damage.

rDNA Transcription in Developmental Diseases and Stem Cells

As the first and rate-limiting step in ribosome biogenesis, rDNA transcription undergoes significant dynamic changes during cell pluripotency alteration. Over the past decades, rDNA activity has demonstrated dynamic changes, but most people view it as passive compliance with cellular needs. The evidence for rDNA transcriptional activity determining stem cell pluripotency is growing as research advances, resulting in the arrest of embryonic development and impairment of stem cell lines stemness by rDNA transcription inhibition. The exact mechanism by which rDNA activation influences pluripotency remains unknown. The first objective of this opinion article is to describe rDNA changes in the pathological and physiological course of life, including developmental diseases, tumor genesis, and stem cell differentiation. After that, we propose three hypotheses regarding rDNA regulation of pluripotency: 1) Specialized ribosomes synthesized from rDNA variant, 2) Nucleolar stress induced by the drop of rDNA transcription, 3) Interchromosomal interactions between rDNA and other genes. The pluripotency regulatory center is expected to focus strongly on rDNA. A small molecule inhibitor of rDNA is used to treat tumors caused by abnormal pluripotency activation. By understanding how rDNA regulates pluripotency, we hope to treat developmental diseases and safely apply somatic cell reprogramming in clinical settings.

Regulation of RNA Polymerase I Stability and Function

RNA polymerase I is a highly processive enzyme with fast initiation and elongation rates. The structure of Pol I, with its in-built RNA cleavage ability and incorporation of subunits homologous to transcription factors, enables it to quickly and efficiently synthesize the enormous amount of rRNA required for ribosome biogenesis. Each step of Pol I transcription is carefully controlled. However, cancers have highjacked these control points to switch the enzyme, and its transcription, on permanently. While this provides an exceptional benefit to cancer cells, it also creates a potential cancer therapeutic vulnerability. We review the current research on the regulation of Pol I transcription, and we discuss chemical biology efforts to develop new targeted agents against this process. Lastly, we highlight challenges that have arisen from the introduction of agents with promiscuous mechanisms of action and provide examples of agents with specificity and selectivity against Pol I.

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.