Ribosomal DNA copy number loss and sequence variation in cancer

Hypermethylation at 45S rDNA promoter in cancers

The ribosomal genes (rDNA genes) encode 47S rRNA which accounts for up to 80% of all cellular RNA. At any given time, no more than 50% of rDNA genes are actively transcribed, and the other half is silent by forming heterochromatin structures through DNA methylation. In cancer cells, upregulation of ribosome biogenesis has been recognized as a hallmark feature, thus, the reduced methylation of rDNA promoter has been thought to support conformational changes of chromatin accessibility and the subsequent increase in rDNA transcription. However, an increase in the heterochromatin state through rDNA hypermethylation can be a protective mechanism teetering on the brink of a threshold where cancer cells rarely successfully proliferate. Hence, clarifying hypo- or hypermethylation of rDNA will unravel its additional cellular functions, including organization of genome architecture and regulation of gene expression, in response to growth signaling, cellular stressors, and carcinogenesis. Using the bisulfite-based quantitative real-time methylation-specific PCR (qMSP) method after ensuring unbiased amplification and complete bisulfite conversion of the minuscule DNA amount of 1 ng, we established that the rDNA promoter was significantly hypermethylated in 107 breast, 65 lung, and 135 colon tumour tissue samples (46.81%, 51.02% and 96.60%, respectively) as compared with their corresponding adjacent normal samples (26.84%, 38.26% and 77.52%, respectively; p < 0.0001). An excessive DNA input of 1 μg resulted in double-stranded rDNA remaining unconverted even after bisulfite conversion, hence the dramatic drop in the single-stranded DNA that strictly required for bisulfite conversion, and leading to an underestimation of rDNA promoter methylation, in other words, a faulty hypomethylation status of the rDNA promoter. Our results are in line with the hypothesis that an increase in rDNA methylation is a natural pathway protecting rDNA repeats that are extremely sensitive to DNA damage in cancer cells.

Insulin signaling regulates R2 retrotransposon expression to orchestrate transgenerational rDNA copy number maintenance

Preserving a large number of essential yet highly unstable ribosomal DNA (rDNA) repeats is critical for the germline to perpetuate the genome through generations. Spontaneous rDNA loss must be countered by rDNA copy number (CN) expansion. Germline rDNA CN expansion is best understood in Drosophila melanogaster, which relies on unequal sister chromatid exchange (USCE) initiated by DNA breaks at rDNA. The rDNA-specific retrotransposon R2 responsible for USCE-inducing DNA breaks is typically expressed only when rDNA CN is low to minimize the danger of DNA breaks; however, the underlying mechanism of R2 regulation remains unclear. Here we identify the insulin receptor (InR) as a major repressor of R2 expression, limiting unnecessary R2 activity. Through single-cell RNA sequencing, we find that male germline stem cells (GSCs), the major cell type that undergoes rDNA CN expansion, have reduced InR expression when rDNA CN is low. Reduced InR activity in turn leads to R2 expression and CN expansion. We further find that dietary manipulation alters R2 expression and rDNA CN expansion activity. This work reveals that the insulin pathway integrates rDNA CN surveying with environmental sensing, revealing a potential mechanism by which diet exerts heritable changes to genomic content.

The complete telomere-to-telomere sequence of a mouse genome

The current reference genome of Mus musculus, GRCm39, has major gaps in both euchromatic and heterochromatic regions associated with repetitive sequences. In this work, we have sequenced and assembled the telomere-to-telomere genome of mouse haploid embryonic stem cells. The results reveal more than 7.7% of previously uncovered sequences of the mouse genome, including ribosomal DNA arrays and pericentromeric and subtelomeric regions, as well as an additional 140 genes predicted to be protein-coding. This study helps to address knowledge gaps in the mouse genome.

Protection of the genome and the central exome by peripheral non-coding DNA against DNA damage in health, ageing and age-related diseases

DNA in eukaryotic genomes is under constant assault from both exogenous and endogenous sources, leading to DNA damage, which is considered a major molecular driver of ageing. Fortunately, the genome and the central exome are safeguarded against these attacks by abundant peripheral non-coding DNA. Non-coding DNA codes for small non-coding RNAs that inactivate foreign nucleic acids in the cytoplasm and physically blocks these attacks in the nucleus. Damage to non-coding DNA produced during such blockage is removed in the form of extrachromosomal circular DNA (eccDNA) through nucleic pore complexes. Consequently, non-coding DNA serves as a line of defence for the exome against DNA damage. The total amount of non-coding DNA/heterochromatin declines with age, resulting in a decrease in both physical blockage and eccDNA exclusion, and thus an increase in the accumulation of DNA damage in the nucleus during ageing and in age-related diseases. Here, we summarize recent evidence supporting a protective role of non-coding DNA in healthy and pathological states and argue that DNA damage is the proximate cause of ageing and age-related genetic diseases. Strategies aimed at strengthening the protective role of non-coding DNA/heterochromatin could potentially offer better systematic protection for the dynamic genome and the exome against diverse assaults, reduce the burden of DNA damage to the exome, and thus slow ageing, counteract age-related genetic diseases and promote a healthier life for individuals.

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.

Hyper-recombination in ribosomal DNA is driven by long-range resection-independent RAD51 accumulation

Ribosomal DNA (rDNA) encodes the ribosomal RNA genes and represents an intrinsically unstable genomic region. However, the underlying mechanisms and implications for genome integrity remain elusive. Here, we use Bloom syndrome (BS), a rare genetic disease characterized by DNA repair defects and hyper-unstable rDNA, as a model to investigate the mechanisms leading to rDNA instability. We find that in Bloom helicase (BLM) proficient cells, the homologous recombination (HR) pathway in rDNA resembles that in nuclear chromatin; it is initiated by resection, replication protein A (RPA) loading and BRCA2-dependent RAD51 filament formation. However, BLM deficiency compromises RPA-loading and BRCA1/2 recruitment to rDNA, but not RAD51 accumulation. RAD51 accumulates at rDNA despite depletion of long-range resection nucleases and rDNA damage results in micronuclei when BLM is absent. In summary, our findings indicate that rDNA is permissive to RAD51 accumulation in the absence of BLM, leading to micronucleation and potentially global genomic instability.

Copy number and sequence variation in rDNA of Daphnia pulex from natural populations: insights from whole-genome sequencing

Ribosomal DNA (rDNA) has a vital role in ribosome biogenesis as it contains the genes that encode ribosomal RNA (rRNA) separated by intergenic spacers (IGSs). The rRNA genes occur in hundreds to tens of thousands of copies per haploid genome in eukaryotes and are generally highly conserved with low variation within species. Due to the repetitive nature and large size of rDNA arrays, detecting intraindividual variation can be difficult. In this study, we use whole-genome sequences of 169 Daphnia pulex individuals from 10 natural populations to measure the copy number and sequence variation in rDNA. This revealed that variation in rDNA copy number between individuals spans an order of magnitude. We further observed a substantial level of sequence variation within individual genomes. As expected, single-nucleotide polymorphisms occurred in regions of lower functional constraint such as the IGS and expansion segments of the rRNA genes. The presence of strong linkage disequilibrium among variants facilitated identification of haplotypes within each population. Although there was evidence of recombination among haplotypes from different populations, it is insufficient to eliminate linkage disequilibrium within populations. Estimating copy number and haplotype diversity within individuals revealed that the level of intraindividual sequence variation is not strongly correlated with copy number. The observed patterns of variation highlight a complex evolutionary history of rDNA in D. pulex. Future research should explore the functional implications of rDNA copy number and sequence variation on organismal phenotypes.

Insights into ribosomal DNA dominance and magnification through characterization of isogenic deletion alleles

The major loci for the large primary ribosomal RNA (rRNA) genes (35S rRNAs) exist as hundreds to thousands of tandem repeats in all organisms and dozens to hundreds in Drosophila. The highly repetitive nature of the ribosomal DNA (rDNA) makes it intrinsically unstable, and many conditions arise from the reduction in or magnification of copy number, but the conditions under which it does so remain unknown. By targeted DNA damage to the rDNA of the Y chromosome, we created and investigated a series of rDNA alleles. We found that complete loss of rDNA leads to lethality after the completion of embryogenesis, blocking larval molting and metamorphosis. We find that the resident retrotransposons-R1 and R2-are regulated by active rDNA such that reduction in copy number derepresses these elements. Their expression is highest during the early first instar, when loss of rDNA is lethal. Regulation of R1 and R2 may be related to their structural arrangement within the rDNA, as we find they are clustered in the flanks of the nucleolus organizing region (NOR; the cytological appearance of the rDNA). We assessed the complex nucleolar dominance relationship between X- and Y-linked rDNA using a histone H3.3-GFP reporter construct and incorporation at the NOR and found that dominance is controlled by rDNA copy number as at high multiplicity the Y-linked array is dominant, but at low multiplicity the X-linked array becomes derepressed. Finally, we found that multiple conditions that disrupt nucleolar dominance lead to increased rDNA magnification, suggesting that the phenomena of dominance and magnification are related, and a single mechanism may underlie and unify these two longstanding observations in Drosophila.

Ribosomal DNA copy number is associated with body mass in humans and other mammals

Body mass results from a complex interplay between genetics and environment. Previous studies of the genetic contribution to body mass have excluded repetitive regions due to the technical limitations of platforms used for population scale studies. Here we apply genome-wide approaches, identifying an association between adult body mass and the copy number (CN) of 47S-ribosomal DNA (rDNA). rDNA codes for the 18 S, 5.8 S and 28 S ribosomal RNA (rRNA) components of the ribosome. In mammals, there are hundreds of copies of these genes. Inter-individual variation in the rDNA CN has not previously been associated with a mammalian phenotype. Here, we show that rDNA CN variation associates with post-pubertal growth rate in rats and body mass index in adult humans. rDNA CN is not associated with rRNA transcription rates in adult tissues, suggesting the mechanistic link occurs earlier in development. This aligns with the observation that the association emerges by early adulthood.

Ribosomal DNA copy number variation associates with hematological profiles and renal function in the UK Biobank

The phenotypic impact of genetic variation of repetitive features in the human genome is currently understudied. One such feature is the multi-copy 47S ribosomal DNA (rDNA) that codes for rRNA components of the ribosome. Here, we present an analysis of rDNA copy number (CN) variation in the UK Biobank (UKB). From the first release of UKB whole-genome sequencing (WGS) data, a discovery analysis in White British individuals reveals that rDNA CN associates with altered counts of specific blood cell subtypes, such as neutrophils, and with the estimated glomerular filtration rate, a marker of kidney function. Similar trends are observed in other ancestries. A range of analyses argue against reverse causality or common confounder effects, and all core results replicate in the second UKB WGS release. Our work demonstrates that rDNA CN is a genetic influence on trait variance in humans.

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.

Comparative analysis and classification of highly divergent mouse rDNA units based on their intergenic spacer (IGS) variability

Ribosomal DNA (rDNA) repeat units are organized into tandem clusters in eukaryotic cells. In mice, these clusters are located on at least eight chromosomes and show extensive variation in the number of repeats between mouse genomes. To analyze intra- and inter-genomic variation of mouse rDNA repeats, we selectively isolated 25 individual rDNA units using Transformation-Associated Recombination (TAR) cloning. Long-read sequencing and subsequent comparative sequence analysis revealed that each full-length unit comprises an intergenic spacer (IGS) and a ∼13.4 kb long transcribed region encoding the three rRNAs, but with substantial variability in rDNA unit size, ranging from ∼35 to ∼46 kb. Within the transcribed regions of rDNA units, we found 209 variants, 70 of which are in external transcribed spacers (ETSs); but the rDNA size differences are driven primarily by IGS size heterogeneity, due to indels containing repetitive elements and some functional signals such as enhancers. Further evolutionary analysis categorized rDNA units into distinct clusters with characteristic IGS lengths; numbers of enhancers; and presence/absence of two common SNPs in promoter regions, one of which is located within promoter (p)RNA and may influence pRNA folding stability. These characteristic features of IGSs also correlated significantly with 5'ETS variant patterns described previously and associated with differential expression of rDNA units. Our results suggest that variant rDNA units are differentially regulated and open a route to investigate the role of rDNA variation on nucleolar formation and possible associations with pathology.

The impact of ribosome biogenesis in cancer: from proliferation to metastasis

The dysregulation of ribosome biogenesis is a hallmark of cancer, facilitating the adaptation to altered translational demands essential for various aspects of tumor progression. This review explores the intricate interplay between ribosome biogenesis and cancer development, highlighting dynamic regulation orchestrated by key oncogenic signaling pathways. Recent studies reveal the multifaceted roles of ribosomes, extending beyond protein factories to include regulatory functions in mRNA translation. Dysregulated ribosome biogenesis not only hampers precise control of global protein production and proliferation but also influences processes such as the maintenance of stem cell-like properties and epithelial-mesenchymal transition, contributing to cancer progression. Interference with ribosome biogenesis, notably through RNA Pol I inhibition, elicits a stress response marked by nucleolar integrity loss, and subsequent G1-cell cycle arrest or cell death. These findings suggest that cancer cells may rely on heightened RNA Pol I transcription, rendering ribosomal RNA synthesis a potential therapeutic vulnerability. The review further explores targeting ribosome biogenesis vulnerabilities as a promising strategy to disrupt global ribosome production, presenting therapeutic opportunities for cancer treatment.

U2AF2-SNORA68 promotes triple-negative breast cancer stemness through the translocation of RPL23 from nucleoplasm to nucleolus and c-Myc expression

BACKGROUND

Small nucleolar RNAs (snoRNAs) play key roles in ribosome biosynthesis. However, the mechanism by which snoRNAs regulate cancer stemness remains to be fully elucidated.

METHODS

SNORA68 expression was evaluated in breast cancer tissues by in situ hybridization and qRT‒PCR. Proliferation, migration, apoptosis and stemness analyses were used to determine the role of SNORA68 in carcinogenesis and stemness maintenance. Mechanistically, RNA pull-down, RNA immunoprecipitation (RIP), cell fractionation and coimmunoprecipitation assays were conducted.

RESULTS

SNORA68 exhibited high expression in triple-negative breast cancer (TNBC) and was significantly correlated with tumor size (P = 0.048), ki-67 level (P = 0.037), and TNM stage (P = 0.015). The plasma SNORA68 concentration was significantly lower in patients who achieved clinical benefit. The SNORA68-high patients had significantly shorter disease-free survival (DFS) (P = 0.036). Functionally, SNORA68 was found to promote the cell stemness and carcinogenesis of TNBC in vitro and in vivo. Furthermore, elevated SNORA68 expression led to increased nucleolar RPL23 expression and retained RPL23 in the nucleolus by binding U2AF2. RPL23 in the nucleolus subsequently upregulated c-Myc expression. This pathway was validated using a xenograft model.

CONCLUSION

U2AF2-SNORA68 promotes TNBC stemness by retaining RPL23 in the nucleolus and increasing c-Myc expression, which provides new insight into the regulatory mechanism of stemness.

Functions of RNAi Pathways in Ribosomal RNA Regulation

Argonaute proteins, guided by small RNAs, play crucial roles in gene regulation and genome protection through RNA interference (RNAi)-related mechanisms. Ribosomal RNAs (rRNAs), encoded by repeated rDNA units, constitute the core of the ribosome being the most abundant cellular transcripts. rDNA clusters also serve as sources of small RNAs, which are loaded into Argonaute proteins and are able to regulate rDNA itself or affect other gene targets. In this review, we consider the impact of small RNA pathways, specifically siRNAs and piRNAs, on rRNA gene regulation. Data from diverse eukaryotic organisms suggest the potential involvement of small RNAs in various molecular processes related to the rDNA transcription and rRNA fate. Endogenous siRNAs are integral to the chromatin-based silencing of rDNA loci in plants and have been shown to repress rDNA transcription in animals. Small RNAs also play a role in maintaining the integrity of rDNA clusters and may function in the cellular response to rDNA damage. Studies on the impact of RNAi and small RNAs on rRNA provide vast opportunities for future exploration.

Insulin signaling regulates R2 retrotransposon expression to orchestrate transgenerational rDNA copy number maintenance

Preserving a large number of essential yet highly unstable ribosomal DNA (rDNA) repeats is critical for the germline to perpetuate the genome through generations. Spontaneous rDNA loss must be countered by rDNA copy number (CN) expansion. Germline rDNA CN expansion is best understood in Drosophila melanogaster, which relies on unequal sister chromatid exchange (USCE) initiated by DNA breaks at rDNA. The rDNA-specific retrotransposon R2 responsible for USCE-inducing DNA breaks is typically expressed only when rDNA CN is low to minimize the danger of DNA breaks; however, the underlying mechanism of R2 regulation remains unclear. Here we identify the insulin receptor (InR) as a major repressor of R2 expression, limiting unnecessary R2 activity. Through single-cell RNA sequencing we find that male germline stem cells (GSCs), the major cell type that undergoes rDNA CN expansion, have reduced InR expression when rDNA CN is low. Reduced InR activity in turn leads to R2 expression and CN expansion. We further find that dietary manipulation alters R2 expression and rDNA CN expansion activity. This work reveals that the insulin pathway integrates rDNA CN surveying with environmental sensing, revealing a potential mechanism by which diet exerts heritable changes to genomic content.

Regulation and function of R-loops at repetitive elements

R-loops are atypical, three-stranded nucleic acid structures that contain a stretch of RNA:DNA hybrids and an unpaired, single stranded DNA loop. R-loops are physiological relevant and can act as regulators of gene expression, chromatin structure, DNA damage repair and DNA replication. However, unscheduled and persistent R-loops are mutagenic and can mediate replication-transcription conflicts, leading to DNA damage and genome instability if left unchecked. Detailed transcriptome analysis unveiled that 85% of the human genome, including repetitive regions, hold transcriptional activity. This anticipates that R-loops management plays a central role for the regulation and integrity of genomes. This function is expected to have a particular relevance for repetitive sequences that make up to 75% of the human genome. Here, we review the impact of R-loops on the function and stability of repetitive regions such as centromeres, telomeres, rDNA arrays, transposable elements and triplet repeat expansions and discuss their relevance for associated pathological conditions.

Transcription factor UBF depletion in mouse cells results in downregulation of both downstream and upstream elements of the rRNA transcription network

Transcription/processing of the ribosomal RNA (rRNA) precursor, as part of ribosome biosynthesis, is intensively studied and characterized in eukaryotic cells. Here, we constructed shRNA-based mouse cell lines partially silenced for the Upstream Binding Factor UBF, the master regulator of rRNA transcription and organizer of open rDNA chromatin. Full Ubf silencing in vivo is not viable, and these new tools allow further characterization of rRNA transcription and its coordination with cellular signaling. shUBF cells display cell cycle G1 delay and reduced 47S rRNA precursor and 28S rRNA at baseline and serum-challenged conditions. Growth-related mTOR signaling is downregulated with the fractions of active phospho-S6 Kinase and pEIF4E translation initiation factor reduced, similar to phosphorylated cell cycle regulator retinoblastoma, pRB, positive regulator of UBF availability/rRNA transcription. Additionally, we find transcription-competent pUBF (Ser484) severely restricted and its interacting initiation factor RRN3 reduced and responsive to extracellular cues. Furthermore, fractional UBF occupancy on the rDNA unit is decreased in shUBF, and expression of major factors involved in different aspects of rRNA transcription is severely downregulated by UBF depletion. Finally, we observe reduced RNA Pol1 occupancy over rDNA promoter sequences and identified unexpected regulation of RNA Pol1 expression, relative to serum availability and under UBF silencing, suggesting that regulation of rRNA transcription may not be restricted to modulation of Pol1 promoter binding/elongation rate. Overall, this work reveals that UBF depletion has a critical downstream and upstream impact on the whole network orchestrating rRNA transcription in mammalian cells.

Senescence in yeast is associated with amplified linear fragments of chromosome XII rather than ribosomal DNA circle accumulation

The massive accumulation of extrachromosomal ribosomal DNA circles (ERCs) in yeast mother cells has been long cited as the primary driver of replicative ageing. ERCs arise through ribosomal DNA (rDNA) recombination, and a wealth of genetic data connects rDNA instability events giving rise to ERCs with shortened life span and other ageing pathologies. However, we understand little about the molecular effects of ERC accumulation. Here, we studied ageing in the presence and absence of ERCs, and unexpectedly found no evidence of gene expression differences that might indicate stress responses or metabolic feedback caused by ERCs. Neither did we observe any global change in the widespread disruption of gene expression that accompanies yeast ageing, altogether suggesting that ERCs are largely inert. Much of the differential gene expression that accompanies ageing in yeast was actually associated with markers of the senescence entry point (SEP), showing that senescence, rather than age, underlies these changes. Cells passed the SEP irrespective of ERCs, but we found the SEP to be associated with copy number amplification of a region of chromosome XII between the rDNA and the telomere (ChrXIIr) forming linear fragments up to approximately 1.8 Mb size, which arise in aged cells due to rDNA instability but through a different mechanism to ERCs. Therefore, although rDNA copy number increases dramatically with age due to ERC accumulation, our findings implicate ChrXIIr, rather than ERCs, as the primary driver of senescence during budding yeast ageing.

Nuclear envelope-remodeling events as models to assess the potential role of membranes on genome stability

The nuclear envelope (NE) encloses the genetic material and functions in chromatin organization and stability. In Saccharomyces cerevisiae, the NE is bound to the ribosomal DNA (rDNA), highly repeated and transcribed, thus prone to genetic instability. While tethering limits instability, it simultaneously triggers notable NE remodeling. We posit here that NE remodeling may contribute to genome integrity maintenance. The NE importance in genome expression, structure, and integrity is well recognized, yet studies mostly focus on peripheral proteins and nuclear pores, not on the membrane itself. We recently characterized a NE invagination drastically obliterating the rDNA, which we propose here as a model to probe if and how membranes play an active role in genome stability preservation.

N6-methyladenosine upregulates ribosome biogenesis in environmental carcinogenesis

Many trace metal pollutants in surface water, the atmosphere, and soil are carcinogenic, and ribosome biogenesis plays an important role in the carcinogenicity of heavy metals. However, the contradiction between upregulated ribosome biogenesis and decreased ribosomal DNA copy number in environmental carcinogenesis is not fully understood. Here, from a perspective of the most predominant and abundant RNA epigenetic modification, N⁶-methyladenosine (m⁶A), we explored the reason behind this contradiction at the post-transcriptional level using arsenite-induced skin carcinogenesis models both in vitro and in vivo. Based on the m⁶A microarray assay and a series of experiments, we found for the first time that the elevated m⁶A in arsenite-induced transformation is mainly enriched in the genes regulating ribosome biogenesis. m⁶A upregulates ribosome biogenesis post-transcriptionally by stabilizing ribosomal proteins and modulating non-coding RNAs targeting ribosomal RNAs and proteins, leading to arsenite-induced skin carcinogenesis. Using multi-omics analysis of human subjects and experimental validation, we identified an unconventional role of a well-known key proliferative signaling node AKT1 as a vital mediator between m⁶A and ribosome biogenesis in arsenic carcinogenesis. m⁶A activates AKT1 and transmits proliferative signals to ribosome biogenesis, exacerbating the upregulation of ribosome biogenesis in arsenite-transformed keratinocytes. Similarly, m⁶A promotes cell proliferation by upregulating ribosome biogenesis in cell transformation induced by carcinogenic heavy metals (chromium and nickel). Importantly, inhibiting m⁶A reduces ribosome biogenesis. Targeted inhibition of m⁶A-upregulated ribosome biogenesis effectively prevents cell transformation induced by trace metals (arsenic, chromium, and nickel). Our results reveal the mechanism of ribosome biogenesis upregulated by m⁶A in the carcinogenesis of trace metal pollutants. From the perspective of RNA epigenetics, our study improves our understanding of the contradiction between upregulated ribosome biogenesis and decreased ribosomal DNA copy number in the carcinogenesis of environmental carcinogens.

Substantial rDNA copy number reductions alter timing of development and produce variable tissue-specific phenotypes in C. elegans

The genes that encode ribosomal RNAs are present in several hundred copies in most eukaryotes. These vast arrays of repetitive ribosomal DNA (rDNA) have been implicated not just in ribosome biogenesis, but also aging, cancer, genome stability, and global gene expression. rDNA copy number is highly variable among and within species; this variability is thought to associate with traits relevant to human health and disease. Here we investigate the phenotypic consequences of multicellular life at the lower bounds of rDNA copy number. We use the model Caenorhabditis elegans, which has previously been found to complete embryogenesis using only maternally provided ribosomes. We find that individuals with rDNA copy number reduced to ∼5% of wild type are capable of further development with variable penetrance. Such individuals are sterile and exhibit severe morphological defects, particularly in post-embryonically dividing tissues such as germline and vulva. Developmental completion and fertility are supported by an rDNA copy number ∼10% of wild type, with substantially delayed development. Worms with rDNA copy number reduced to ∼33% of wild type display a subtle developmental timing defect that was absent in worms with higher copy numbers. Our results support the hypothesis that rDNA requirements vary across tissues and indicate that the minimum rDNA copy number for fertile adulthood is substantially less than the lowest naturally observed total copy number. The phenotype of individuals with severely reduced rDNA copy number is highly variable in penetrance and presentation, highlighting the need for continued investigation into the biological consequences of rDNA copy number variation.

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.

RNA Polymerases I and III in development and disease

Ribosomes are macromolecular machines that are globally required for the translation of all proteins in all cells. Ribosome biogenesis, which is essential for cell growth, proliferation and survival, commences with transcription of a variety of RNAs by RNA Polymerases I and III. RNA Polymerase I (Pol I) transcribes ribosomal RNA (rRNA), while RNA Polymerase III (Pol III) transcribes 5S ribosomal RNA and transfer RNAs (tRNA) in addition to a wide variety of small non-coding RNAs. Interestingly, despite their global importance, disruptions in Pol I and Pol III function result in tissue-specific developmental disorders, with craniofacial anomalies and leukodystrophy/neurodegenerative disease being among the most prevalent. Furthermore, pathogenic variants in genes encoding subunits shared between Pol I and Pol III give rise to distinct syndromes depending on whether Pol I or Pol III function is disrupted. In this review, we discuss the global roles of Pol I and III transcription, the consequences of disruptions in Pol I and III transcription, disorders arising from pathogenic variants in Pol I and Pol III subunits, and mechanisms underpinning their tissue-specific phenotypes.

Ribosomal DNA replication time coordinates completion of genome replication and anaphase in yeast

Timely completion of genome replication is a prerequisite for mitosis, genome integrity, and cell survival. A challenge to this timely completion comes from the need to replicate the hundreds of untranscribed copies of rDNA that organisms maintain in addition to the copies required for ribosome biogenesis. Replication of these rDNA arrays is relegated to late S phase despite their large size, repetitive nature, and essentiality. Here, we show that, in Saccharomyces cerevisiae, reducing the number of rDNA repeats leads to early rDNA replication, which results in delaying replication elsewhere in the genome. Moreover, cells with early-replicating rDNA arrays and delayed genome-wide replication aberrantly release the mitotic phosphatase Cdc14 from the nucleolus and enter anaphase prematurely. We propose that rDNA copy number determines the replication time of the rDNA locus and that the release of Cdc14 upon completion of rDNA replication is a signal for cell cycle progression.

Repurposing Dihydroartemisinin to Combat Oral Squamous Cell Carcinoma, Associated with Mitochondrial Dysfunction and Oxidative Stress

Oral squamous cell carcinoma (OSCC), with aggressive locoregional invasion, has a high rate of early recurrences and poor prognosis. Dihydroartemisinin (DHA), as a derivative of artemisinin, has been found to exert potent antitumor activity. Recent studies reported that DHA suppresses OSCC cell growth and viability through the regulation of reactive oxygen species (ROS) production and mitochondrial calcium uniporter. However, the mechanism underlying the action of DHA on OSCCs remains elusive. In the study, we observed that 159 genes were remarkably misregulated in primary OSCC tumors associated with DHA-inhibited pathways, supporting that OSCCs are susceptible to DHA treatment. Herein, our study showed that DHA exhibited promising effects to suppress OSCC cell growth and survival, and single-cell colony formation. Interestingly, the combination of DHA and cisplatin (CDDP) significantly reduced the toxicity of CDDP treatment alone on human normal oral cells (NOK). Moreover, DHA remarkably impaired mitochondrial structure and function, and triggered DNA damage and ROS generation, and activation of mitophagy. In addition, DHA induced leakage of cytochrome C and apoptosis-inducing factor (AIF) from mitochondria, elevated Bax/cleaved-caspase 3 expression levels and compromised Bcl2 protein expression. In the OSCC tumor-xenograft mice model, DHA remarkably suppressed tumor growth and induced apoptosis of OSCCs in vivo. Intriguingly, a selective mitophagy inhibitor Mdivi-1 could significantly reinforce the anticancer activity of DHA treatment. DHA and Mdivi-1 can synergistically suppress OSCC cell proliferation and survival. These data uncover a previously unappreciated contribution of the mitochondria-associated pathway to the antitumor activity of DHA on OSCCs. Our study shed light on a new aspect of a DHA-based therapeutic strategy to combat OSCC tumors.

Proteome diversification by mRNA translation in cancer

mRNA translation is a highly conserved and tightly controlled mechanism for protein synthesis and is well known to be altered by oncogenes to promote cancer development. This distorted mRNA translation is accompanied by the vulnerability of cancer to inhibitors of key mRNA translation components. Novel studies also suggest that these alternations could be utilized for immunotherapy. Ribosome heterogeneity and alternative responses to nutrient shortages, which aid cancer growth and spread, are proposed to elicit aberrant protein production but may also result in previously unidentified therapeutic targets, such as the presentation of cancer-specific peptides at the surface of cancer cells (neoepitopes). This review will assess the driving forces in tRNA and ribosome function that underlie proteome diversification due to alterations in mRNA translation in cancer cells.

Fragile sites, chromosomal lesions, tandem repeats, and disease

Expanded tandem repeat DNAs are associated with various unusual chromosomal lesions, despiralizations, multi-branched inter-chromosomal associations, and fragile sites. Fragile sites cytogenetically manifest as localized gaps or discontinuities in chromosome structure and are an important genetic, biological, and health-related phenomena. Common fragile sites (∼230), present in most individuals, are induced by aphidicolin and can be associated with cancer; of the 27 molecularly-mapped common sites, none are associated with a particular DNA sequence motif. Rare fragile sites ( ≳ 40 known), ≤ 5% of the population (may be as few as a single individual), can be associated with neurodevelopmental disease. All 10 molecularly-mapped folate-sensitive fragile sites, the largest category of rare fragile sites, are caused by gene-specific CGG/CCG tandem repeat expansions that are aberrantly CpG methylated and include FRAXA, FRAXE, FRAXF, FRA2A, FRA7A, FRA10A, FRA11A, FRA11B, FRA12A, and FRA16A. The minisatellite-associated rare fragile sites, FRA10B, FRA16B, can be induced by AT-rich DNA-ligands or nucleotide analogs. Despiralized lesions and multi-branched inter-chromosomal associations at the heterochromatic satellite repeats of chromosomes 1, 9, 16 are inducible by de-methylating agents like 5-azadeoxycytidine and can spontaneously arise in patients with ICF syndrome (Immunodeficiency Centromeric instability and Facial anomalies) with mutations in genes regulating DNA methylation. ICF individuals have hypomethylated satellites I-III, alpha-satellites, and subtelomeric repeats. Ribosomal repeats and subtelomeric D4Z4 megasatellites/macrosatellites, are associated with chromosome location, fragility, and disease. Telomere repeats can also assume fragile sites. Dietary deficiencies of folate or vitamin B12, or drug insults are associated with megaloblastic and/or pernicious anemia, that display chromosomes with fragile sites. The recent discovery of many new tandem repeat expansion loci, with varied repeat motifs, where motif lengths can range from mono-nucleotides to megabase units, could be the molecular cause of new fragile sites, or other chromosomal lesions. This review focuses on repeat-associated fragility, covering their induction, cytogenetics, epigenetics, cell type specificity, genetic instability (repeat instability, micronuclei, deletions/rearrangements, and sister chromatid exchange), unusual heritability, disease association, and penetrance. Understanding tandem repeat-associated chromosomal fragile sites provides insight to chromosome structure, genome packaging, genetic instability, and disease.

53BP1-mediated recruitment of RASSF1A to ribosomal DNA breaks promotes local ATM signaling

DNA lesions occur across the genome and constitute a threat to cell viability; however, damage at specific genomic loci has a relatively greater impact on overall genome stability. The ribosomal RNA gene repeats (rDNA) are emerging fragile sites. Recent progress in understanding how the rDNA damage response is organized has highlighted a key role of adaptor proteins. Here, we show that the scaffold tumor suppressor RASSF1A is recruited to rDNA breaks. RASSF1A recruitment to double-strand breaks is mediated by 53BP1 and depends on RASSF1A phosphorylation at Serine 131 by ATM kinase. Employing targeted rDNA damage, we uncover that RASSF1A recruitment promotes local ATM signaling. RASSF1A silencing, a common epigenetic event during malignant transformation, results in persistent breaks, rDNA copy number alterations and decreased cell viability. Overall, we identify a novel role for RASSF1A at rDNA break sites, provide mechanistic insight into how the DNA damage response is organized in a chromatin context, and provide further evidence for how silencing of the RASSF1A tumor suppressor contributes to genome instability.

A new method for determining ribosomal DNA copy number shows differences between Saccharomyces cerevisiae populations

Ribosomal DNA genes (rDNA) encode the major ribosomal RNAs and in eukaryotes typically form tandem repeat arrays. Species have characteristic rDNA copy numbers, but there is substantial intra-species variation in copy number that results from frequent rDNA recombination. Copy number differences can have phenotypic consequences, however difficulties in quantifying copy number mean we lack a comprehensive understanding of how copy number evolves and the consequences. Here we present a genomic sequence read approach to estimate rDNA copy number based on modal coverage to help overcome limitations with existing mean coverage-based approaches. We validated our method using Saccharomyces cerevisiae strains with known rDNA copy numbers. Application of our pipeline to a global sample of S. cerevisiae isolates showed that different populations have different rDNA copy numbers. Our results demonstrate the utility of the modal coverage method, and highlight the high level of rDNA copy number variation within and between populations.

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.

Widespread genetic heterogeneity of human ribosomal RNA genes

Polymorphism drives survival under stress and provides adaptability. Genetic polymorphism of ribosomal RNA (rRNA) genes derives from internal repeat variation of this multicopy gene, and from interindividual variation. A considerable amount of rRNA sequence heterogeneity has been proposed but has been challenging to estimate given the scarcity of accurate reference sequences. We identified four rDNA copies on chromosome 21 (GRCh38) with 99% similarity to recently introduced reference sequence KY962518.1. We customized a GATK bioinformatics pipeline using the four rDNA loci, spanning a total 145 kb, for variant calling and used high-coverage whole-genome sequencing (WGS) data from the 1000 Genomes Project to analyze variants in 2504 individuals from 26 populations. We identified a total of 3791 variant positions. The variants positioned nonrandomly on the rRNA gene. Invariant regions included the promoter, early 5' ETS, most of 18S, 5.8S, ITS1, and large areas of the intragenic spacer. A total of 470 variant positions were observed on 28S rRNA. The majority of the 28S rRNA variants were located on highly flexible human-expanded rRNA helical folds ES7L and ES27L, suggesting that these represent positions of diversity and are potentially under continuous evolution. Several variants were validated based on RNA-seq analyses. Population analyses showed remarkable ancestry-linked genetic variance and the presence of both high penetrance and frequent variants in the 5' ETS, ITS2, and 28S regions segregating according to the continental populations. These findings provide a genetic view of rRNA gene array heterogeneity and raise the need to functionally assess how the 28S rRNA variants affect ribosome functions.

Short arms of human acrocentric chromosomes and the completion of the human genome sequence

The complete, ungapped sequence of the short arms of human acrocentric chromosomes (SAACs) is still unknown almost 20 years after the near completion of the Human Genome Project. Yet these short arms of Chromosomes 13, 14, 15, 21, and 22 contain the ribosomal DNA (rDNA) genes, which are of paramount importance for human biology. The sequences of SAACs show an extensive variation in the copy number of the various repetitive elements, the full extent of which is currently unknown. In addition, the full spectrum of repeated sequences, their organization, and the low copy number functional elements are also unknown. The Telomere-to-Telomere (T2T) Project using mainly long-read sequence technology has recently completed the assembly of the genome from a hydatidiform mole, CHM13, and has thus established a baseline reference for further studies on the organization, variation, functional annotation, and impact in human disorders of all the previously unknown genomic segments, including the SAACs. The publication of the initial results of the T2T Project will update and improve the reference genome for a better understanding of the evolution and function of the human genome.

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.

Global genomic instability caused by reduced expression of DNA polymerase ε in yeast

SignificanceAlthough most studies of the genetic regulation of genome stability involve an analysis of mutations within the coding sequences of genes required for DNA replication or DNA repair, recent studies in yeast show that reduced levels of wild-type enzymes can also produce a mutator phenotype. By whole-genome sequencing and other methods, we find that reduced levels of the wild-type DNA polymerase ε in yeast greatly increase the rates of mitotic recombination, aneuploidy, and single-base mutations. The observed pattern of genome instability is different from those observed in yeast strains with reduced levels of the other replicative DNA polymerases, Pol α and Pol δ. These observations are relevant to our understanding of cancer and other diseases associated with genetic instability.

Sexual Crossing, Chromosome-Level Genome Sequences, and Comparative Genomic Analyses for the Medicinal Mushroom Taiwanofungus Camphoratus (Syn. Antrodia Cinnamomea, Antrodia Camphorata)

Taiwanofungus camphoratus mushrooms are a complementary and alternative medicine for hangovers, cancer, hypertension, obesity, diabetes, and inflammation. Though Taiwanofungus camphoratus has attracted considerable biotechnological and pharmacological attention, neither classical genetic nor genomic approaches have been properly established for it. We isolated four sexually competent monokaryons from two T. camphoratus dikaryons used for the commercial cultivation of orange-red (HC1) and milky-white (SN1) mushrooms, respectively. We also sequenced, annotated, and comparatively analyzed high-quality and chromosome-level genome sequences of these four monokaryons. These genomic resources represent a valuable basis for understanding the biology, evolution, and secondary metabolite biosynthesis of this economically important mushrooms. We demonstrate that T. camphoratus has a tetrapolar mating system and that HC1 and SN1 represent two intraspecies isolates displaying karyotypic variation. Compared with several edible mushroom model organisms, T. camphoratus underwent a significant contraction in the gene family and individual gene numbers, most notably for plant, fungal, and bacterial cell-wall-degrading enzymes, explaining why T. camphoratus mushrooms are rare in natural environments, are difficult and time-consuming to artificially cultivate, and are susceptible to fungal and bacterial infections. Our results lay the foundation for an in-depth T. camphoratus study, including precise genetic manipulation, improvements to mushroom fruiting, and synthetic biology applications for producing natural medicinal products. IMPORTANCETaiwanofungus camphoratus (Tc) is a basidiomycete fungus that causes brown heart rot of the aromatic tree Cinnamomum kanehirae. The Tc fruiting bodies have been used to treat hangovers, abdominal pain, diarrhea, hypertension, and other diseases first by aboriginal Taiwanese and later by people in many countries. To establish classical genetic and genomic approaches for this economically important medicinal mushroom, we first isolated and characterized four sexually competent monokaryons from two dikaryons wildly used for commercial production of Tc mushrooms. We applied PacBio single molecule, real-time sequencing technology to determine the near-completed genome sequences of four monokaryons. These telomere-to-telomere and gapless haploid genome sequences reveal all genomic variants needed to be studied and discovered, including centromeres, telomeres, retrotransposons, mating type loci, biosynthetic, and metabolic gene clusters. Substantial interspecies diversities are also discovered between Tc and several other mushroom model organisms, including Agrocybe aegerita, Coprinopsis cinerea, and Schizophyllum commune, and Ganoderma lucidum.

The Psychoemotional Stress-Induced Changes in the Abundance of SatIII (1q12) and Telomere Repeats, but Not Ribosomal DNA, in Human Leukocytes

INTRODUCTION. As shown earlier, copy number variations (CNV) in the human satellite III (1q12) fragment (f-SatIII) and the telomere repeat (TR) reflects the cell’s response to oxidative stress. The contents of f-SatIII and TR in schizophrenic (SZ) patients were found to be lower than in healthy controls (HC) in previous studies. The major question of this study was: ‘What are the f-SatIII and TR CNV dynamic changes in human leukocytes, depending on psychoemotional stress?’ MATERIALS AND METHODS. We chose a model of psychoemotional stress experienced by second-year medical students during their exams. Blood samples were taken in stressful conditions (exams) and in a control non-stressful period. Biotinylated probes were used for f-SatIII, rDNA, and TR quantitation in leukocyte DNA by non-radioactive quantitative hybridization in SZ patients (n = 97), HC (n = 97), and medical students (n = 17, n = 42). A flow cytometry analysis was used for the oxidative stress marker (NOX4, 8-oxodG, and γH2AX) detection in the lymphocytes of the three groups. RESULTS. Oxidative stress markers increased significantly in the students’ lymphocytes during psychoemotional stress. The TR and f-SatIII, but not the rDNA, contents significantly changed in the DNA isolated from human blood leukocytes. After a restoration period (post-examinational vacations), the f-SatIII content decreased, and the TR content increased. Changes in the blood cells of students during examinational stress were similar to those in SZ patients during an exacerbation of the disease. CONCLUSIONS. Psychoemotional stress in students during exams triggers a universal mechanism of oxidative stress. The oxidative stress causes significant changes in the f-SatIII and TR contents, while the ribosomal repeat content remains stable. A hypothesis is proposed to explain the quantitative polymorphisms of f-SatIII and TR contents under transient (e.g., students’ exams) or chronic (in SZ patients) stress. The changes in the f-SatIII and TR copy numbers are non-specific events, irrespective of the source of stress. Thus, our findings suggest that the psychoemotional stress, common in SZ patients and healthy students during exams, but not in a schizophrenia-specific event, was responsible for the changes in the repeat contents that we observed earlier in SZ patients.

Epigenetics of Ribosomal RNA Genes

This review is focused on biology of genes encoding ribosomal RNA (rRNA) in mammals. rRNA is a structural component of the most abundant cellular molecule, the ribosome. There are many copies of rRNA genes per genome that are under tight transcriptional control by epigenetic mechanisms serving to meet cellular needs in protein synthesis. Curiously, only a fraction of rRNA genes is used even in the fast-growing cells, raising a question why unused copies of these genes have not been lost during evolution. Two plausible explanations are discussed. First, there is evidence that besides their direct function in production of rRNA, ribosomal RNA genes are involved in regulation of many other genes in the nucleus by forming either temporary or persistent complexes with these genes. Second, it seems that rRNA genes also play a role in the maintenance of genome stability, where lower copy number of rRNA genes destabilizes the genome. These "additional" functions of rRNA genes make them recurrent candidate drivers of chronic human diseases and aging. Experimental support for the involvement of these genes in human diseases and potential mechanisms are also discussed.

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.

Human rDNA and Cancer

In human cells, each rDNA unit consists of the ~13 kb long ribosomal part and ~30 kb long intergenic spacer (IGS). The ribosomal part, transcribed by RNA polymerase I (pol I), includes genes coding for 18S, 5.8S, and 28S RNAs of the ribosomal particles, as well as their four transcribed spacers. Being highly repetitive, intensively transcribed, and abundantly methylated, rDNA is a very fragile site of the genome, with high risk of instability leading to cancer. Multiple small mutations, considerable expansion or contraction of the rDNA locus, and abnormally enhanced pol I transcription are usual symptoms of transformation. Recently it was found that both IGS and the ribosomal part of the locus contain many functional/potentially functional regions producing non-coding RNAs, which participate in the pol I activity regulation, stress reactions, and development of the malignant phenotype. Thus, there are solid reasons to believe that rDNA locus plays crucial role in carcinogenesis. In this review we discuss the data concerning the human rDNA and its closely associated factors as both targets and drivers of the pathways essential for carcinogenesis. We also examine whether variability in the structure of the locus may be blamed for the malignant transformation. Additionally, we consider the prospects of therapy focused on the activity of rDNA.

PCR amplicons identify widespread copy number variation in human centromeric arrays and instability in cancer

Centromeric α-satellite repeats represent ~6% of the human genome, but their length and repetitive nature make sequencing and analysis of those regions challenging. However, centromeres are essential for the stable propagation of chromosomes, so tools are urgently needed to monitor centromere copy number and how it influences chromosome transmission and genome stability. We developed and benchmarked droplet digital PCR (ddPCR) assays that measure copy number for five human centromeric arrays. We applied them to characterize natural variation in centromeric array size, analyzing normal tissue from 37 individuals from China and 39 individuals from the US and UK. Each chromosome-specific array varies in size up to 10-fold across individuals and up to 50-fold across chromosomes, indicating a unique complement of arrays in each individual. We also used the ddPCR assays to analyze centromere copy number in 76 matched tumor-normal samples across four cancer types, representing the most-comprehensive quantitative analysis of centromeric array stability in cancer to date. In contrast to stable transmission in cultured cells, centromeric arrays show gain and loss events in each of the cancer types, suggesting centromeric α-satellite DNA represents a new category of genome instability in cancer. Our methodology for measuring human centromeric-array copy number will advance research on centromeres and genome integrity in normal and disease states.

Concerted evolution of ribosomal DNA: Somatic peace amid germinal strife: Intranuclear and cellular selection maintain the quality of rRNA

Most eukaryotes possess many copies of rDNA. Organismal selection alone cannot maintain rRNA function because the effects of mutations in one rDNA are diluted by the presence of many other rDNAs. rRNA quality is maintained by processes that increase homogeneity of rRNA within, and heterogeneity among, germ cells thereby increasing the effectiveness of cellular selection on ribosomal function. A successful rDNA repeat will possess adaptations for spreading within tandem arrays by intranuclear selection. These adaptations reside in the non-coding regions of rDNA. Single-copy genes are predicted to manage processes of intranuclear and cellular selection in the germline to maintain the quality of rRNA expressed in somatic cells of future generations.

Ribosomal proteins and human diseases: molecular mechanisms and targeted therapy

Ribosome biogenesis and protein synthesis are fundamental rate-limiting steps for cell growth and proliferation. The ribosomal proteins (RPs), comprising the structural parts of the ribosome, are essential for ribosome assembly and function. In addition to their canonical ribosomal functions, multiple RPs have extra-ribosomal functions including activation of p53-dependent or p53-independent pathways in response to stress, resulting in cell cycle arrest and apoptosis. Defects in ribosome biogenesis, translation, and the functions of individual RPs, including mutations in RPs have been linked to a diverse range of human congenital disorders termed ribosomopathies. Ribosomopathies are characterized by tissue-specific phenotypic abnormalities and higher cancer risk later in life. Recent discoveries of somatic mutations in RPs in multiple tumor types reinforce the connections between ribosomal defects and cancer. In this article, we review the most recent advances in understanding the molecular consequences of RP mutations and ribosomal defects in ribosomopathies and cancer. We particularly discuss the molecular basis of the transition from hypo- to hyper-proliferation in ribosomopathies with elevated cancer risk, a paradox termed "Dameshek's riddle." Furthermore, we review the current treatments for ribosomopathies and prospective therapies targeting ribosomal defects. We also highlight recent advances in ribosome stress-based cancer therapeutics. Importantly, insights into the mechanisms of resistance to therapies targeting ribosome biogenesis bring new perspectives into the molecular basis of cancer susceptibility in ribosomopathies and new clinical implications for cancer therapy.

Environmentally induced ribosomal DNA (rDNA) instability in human cells and populations exposed to hexavalent chromium [Cr (VI)]

Hexavalent Chromium [Cr (VI)] is an established toxicant, carcinogen, and a significant source of public health concern. The multicopy ribosomal DNA (rDNA) array is mechanistically linked to aging and cancer, is the most evolutionarily conserved segment of the human genome, and gives origin to nucleolus, a nuclear organelle where ribosomes are assembled. Here we show that exposure to Cr (VI) induces instability in the rDNA, triggering cycles of rapid, specific, and transient amplification and contraction of the array in human cells. The dynamic of environmentally responsive rDNA copy number (CN) amplification and contraction occurs at doses to which millions of individuals are regularly exposed. Finally, analyses of human populations occupationally exposed to Cr (VI) indicate that environmental exposure history and drinking habits but not age shape extensive naturally occurring rDNA copy number variation. Our observations identify a novel pathway of response to hexavalent chromium exposure and raise the prospect that a suite of environmental determinants of rDNA copy number remain to be discovered.

Harnessing the Nucleolar DNA Damage Response in Cancer Therapy

The nucleoli are subdomains of the nucleus that form around actively transcribed ribosomal RNA (rRNA) genes. They serve as the site of rRNA synthesis and processing, and ribosome assembly. There are 400-600 copies of rRNA genes (rDNA) in human cells and their highly repetitive and transcribed nature poses a challenge for DNA repair and replication machineries. It is only in the last 7 years that the DNA damage response and processes of DNA repair at the rDNA repeats have been recognized to be unique and distinct from the classic response to DNA damage in the nucleoplasm. In the last decade, the nucleolus has also emerged as a central hub for coordinating responses to stress via sequestering tumor suppressors, DNA repair and cell cycle factors until they are required for their functional role in the nucleoplasm. In this review, we focus on features of the rDNA repeats that make them highly vulnerable to DNA damage and the mechanisms by which rDNA damage is repaired. We highlight the molecular consequences of rDNA damage including activation of the nucleolar DNA damage response, which is emerging as a unique response that can be exploited in anti-cancer therapy. In this review, we focus on CX-5461, a novel inhibitor of Pol I transcription that induces the nucleolar DNA damage response and is showing increasing promise in clinical investigations.

Expansion of rDNA and pericentromere satellite repeats in the genomes of bank voles Myodes glareolus exposed to environmental radionuclides

Altered copy number of certain highly repetitive regions of the genome, such as satellite DNA within heterochromatin and ribosomal RNA loci (rDNA), is hypothesized to help safeguard the genome against damage derived from external stressors. We quantified copy number of the 18S rDNA and a pericentromeric satellite DNA (Msat-160) in bank voles (Myodes glareolus) inhabiting the Chernobyl Exclusion Zone (CEZ), an area that is contaminated by radionuclides and where organisms are exposed to elevated levels of ionizing radiation. We found a significant increase in 18S rDNA and Msat-160 content in the genomes of bank voles from contaminated locations within the CEZ compared with animals from uncontaminated locations. Moreover, 18S rDNA and Msat-160 copy number were positively correlated in the genomes of bank voles from uncontaminated, but not in the genomes of animals inhabiting contaminated, areas. These results show the capacity for local-scale geographic variation in genome architecture and are consistent with the genomic safeguard hypothesis. Disruption of cellular processes related to genomic stability appears to be a hallmark effect in bank voles inhabiting areas contaminated by radionuclides.

Cytogenetic Evidence for Sex Chromosomes and Karyotype Evolution in Anguimorphan Lizards

Anguimorphan lizards are a morphologically variable group of squamate reptiles with a wide geographical distribution. In spite of their importance, they have been cytogenetically understudied. Here, we present the results of the cytogenetic examination of 23 species from five anguimorphan families (Anguidae, Helodermatidae, Shinisauridae, Varanidae and Xenosauridae). We applied both conventional (Giemsa staining and C-banding) and molecular cytogenetic methods (fluorescence in situ hybridization with probes for the telomeric motifs and rDNA loci, comparative genome hybridization), intending to describe the karyotypes of previously unstudied species, to uncover the sex determination mode, and to reveal the distribution of variability in cytogenetic characteristics among anguimorphan lizards. We documented that karyotypes are generally quite variable across anguimorphan lineages, with anguids being the most varying. However, the derived chromosome number of 2n = 40 exhibits a notable long-term evolutionary stasis in monitors. Differentiated ZZ/ZW sex chromosomes were documented in monitors and helodermatids, as well as in the anguids Abronia lythrochila, and preliminary also in Celestus warreni and Gerrhonotus liocephalus. Several other anguimorphan species have likely poorly differentiated sex chromosomes, which cannot be detected by the applied cytogenetic methods, although the presence of environmental sex determination cannot be excluded. In addition, we uncovered a rare case of spontaneous triploidy in a fully grown Varanus primordius.