rDNA Clusters Make Contact with Genes that Are Involved in Differentiation and Cancer and Change Contacts after Heat Shock Treatment

Homeotic DUX4 Genes Shape Dynamic Inter-Chromosomal Contacts with Nucleoli in Human Cells

Nucleoli form interchromosomal contacts with genes controlling differentiation and carcinogenesis. DUX4 genes specify transcription factor possessing two homeodomains. Previously, using Circular Chromosome Conformation Capture (4С) approach on population of cells, it was demonstrated that DUX4 gene clusters form frequent contacts with nucleoli. It was found also that these contacts are almost completely abolished after heat shock treatment. 4C approach as all ligation-mediated methods is capable to detect rather close interactions between chromatin loops in nuclei. In order to independently confirm the formation and the frequency of the contacts in single cells we used FISH approach. Here, we show that DUX genes in single cells form stable contacts in all tested HEK293T cells. During heat shock, DUX4 genes reversibly move 1-3 µm away from the nuclei. We conclude that interchromosomal contacts formed by nucleoli are strong, dynamic, and reversible, providing both the initiation and maintenance of a differentiated state.

Preferential Co-Expression and Colocalization of rDNA-Contacting Genes with LincRNAs Suggest Their Involvement in Shaping Inter-Chromosomal Interactions with Nucleoli

Different developmental genes shape frequent dynamic inter-chromosomal contacts with rDNA units in human and Drosophila cells. In the course of differentiation, changes in these contacts occur, coupled with changes in the expression of hundreds of rDNA-contacting genes. The data suggest a possible role of nucleoli in the global regulation of gene expression. However, the mechanism behind the specificity of these inter-chromosomal contacts, which are rebuilt in every cell cycle, is not yet known. Here, we describe the strong association of rDNA-contacting genes with numerous long intergenic non-coding RNAs (lincRNAs) in HEK293T cells and in initial and differentiated K562 cells. We observed that up to 600 different lincRNAs were preferentially co-expressed with multiple overlapping sets of rDNA-contacting developmental genes, and there was a strong correlation between the genomic positions of rDNA-contacting genes and lincRNA mappings. These two findings suggest that lincRNAs might guide the corresponding developmental genes toward rDNA clusters. We conclude that the inter-chromosomal interactions of rDNA-contacting genes with nucleoli might be guided by lincRNAs, which might physically link particular genomic regions with rDNA clusters.

Induction of the Erythroid Differentiation of K562 Cells Is Coupled with Changes in the Inter-Chromosomal Contacts of rDNA Clusters

The expression of clusters of rDNA genes influences pluripotency; however, the underlying mechanisms are not yet known. These clusters shape inter-chromosomal contacts with numerous genes controlling differentiation in human and Drosophila cells. This suggests a possible role of these contacts in the formation of 3D chromosomal structures and the regulation of gene expression in development. However, it has not yet been demonstrated whether inter-chromosomal rDNA contacts are changed during differentiation. In this study, we used human leukemia K562 cells and induced their erythroid differentiation in order to study both the changes in rDNA contacts and the expression of genes. We observed that approximately 200 sets of rDNA-contacting genes are co-expressed in different combinations in both untreated and differentiated K562 cells. rDNA contacts are changed during differentiation and coupled with the upregulation of genes whose products are mainly located in the nucleus and are highly associated with DNA- and RNA-binding, along with the downregulation of genes whose products mainly reside in the cytoplasm or intra- or extracellular vesicles. The most downregulated gene is ID3, which is known as an inhibitor of differentiation, and thus should be switched off to allow for differentiation. Our data suggest that the differentiation of K562 cells leads to alterations in the inter-chromosomal contacts of rDNA clusters and 3D structures in particular chromosomal regions as well as to changes in the expression of genes located in the corresponding chromosomal domains. We conclude that approximately half of the rDNA-contacting genes are co-expressed in human cells and that rDNA clusters are involved in the global regulation of gene expression.

Genes Possessing the Most Frequent DNA DSBs Are Highly Associated with Development and Cancers, and Essentially Overlap with the rDNA-Contacting Genes

Double-strand DNA breakes (DSBs) are the most deleterious and widespread examples of DNA damage. They inevitably originate from endogenous mechanisms in the course of transcription, replication, and recombination, as well as from different exogenous factors. If not properly repaired, DSBs result in cell death or diseases. Genome-wide analysis of DSBs has revealed the numerous endogenous DSBs in human chromosomes. However, until now, it has not been clear what kind of genes are preferentially subjected to breakage. We performed a genetic and epigenetic analysis of the most frequent DSBs in HEK293T cells. Here, we show that they predominantly occur in the active genes controlling differentiation, development, and morphogenesis. These genes are highly associated with cancers and other diseases. About one-third of the genes possessing frequent DSBs correspond to rDNA-contacting genes. Our data suggest that a specific set of active genes controlling morphogenesis are the main targets of DNA breakage in human cells, although there is a specific set of silent genes controlling metabolism that also are enriched in DSBs. We detected this enrichment by different activators and repressors of transcription at DSB target sites, as well breakage at promoters. We propose that both active transcription and silencing of genes give a propensity for DNA breakage. These results have implications for medicine and gene therapy.

Fragments of rDNA Genes Scattered over the Human Genome Are Targets of Small RNAs

Small noncoding RNAs of different origins and classes play several roles in the regulation of gene expression. Here, we show that diverged and rearranged fragments of rDNA units are scattered throughout the human genome and that endogenous small noncoding RNAs are processed by the Microprocessor complex from specific regions of ribosomal RNAs shaping hairpins. These small RNAs correspond to particular sites inside the fragments of rDNA that mostly reside in intergenic regions or the introns of about 1500 genes. The targets of these small ribosomal RNAs (srRNAs) are characterized by a set of epigenetic marks, binding sites of Pol II, RAD21, CBP, and P300, DNase I hypersensitive sites, and by enrichment or depletion of active histone marks. In HEK293T cells, genes that are targeted by srRNAs (srRNA target genes) are involved in differentiation and development. srRNA target genes are enriched with more actively transcribed genes. Our data suggest that remnants of rDNA sequences and srRNAs may be involved in the upregulation or downregulation of a specific set of genes in human cells. These results have implications for diverse fields, including epigenetics and gene therapy.

The Role of rDNA Clusters in Global Epigenetic Gene Regulation

The regulation of gene expression has been studied for decades, but the underlying mechanisms are still not fully understood. As well as local and distant regulation, there are specific mechanisms of regulation during development and physiological modulation of gene activity in differentiated cells. Current research strongly supports a role for the 3D chromosomal structure in the regulation of gene expression. However, it is not known whether the genome structure reflects the formation of active or repressed chromosomal domains or if these structures play a primary role in the regulation of gene expression. During early development, heterochromatinization of ribosomal DNA (rDNA) is coupled with silencing or activation of the expression of different sets of genes. Although the mechanisms behind this type of regulation are not known, rDNA clusters shape frequent inter-chromosomal contacts with a large group of genes controlling development. This review aims to shed light on the involvement of clusters of ribosomal genes in the global regulation of gene expression. We also discuss the possible role of RNA-mediated and phase-separation mechanisms in the global regulation of gene expression by nucleoli.

Zooming in: PAGE-Northern Blot Helps to Analyze Anti-Sense Transcripts Originating from Human rIGS under Transcriptional Stress

Ribosomal intergenic spacer (rIGS), located between the 45S rRNA coding arrays in humans, is a deep, unexplored source of small and long non-coding RNA molecules transcribed in certain conditions to help a cell generate a stress response, pass through a differentiation state or fine tune the functioning of the nucleolus as a ribosome biogenesis center of the cell. Many of the non-coding transcripts originating from the rIGS are not characterized to date. Here, we confirm the transcriptional activity of the region laying a 2 kb upstream of the rRNA promoter, and demonstrate its altered expression under transcriptional stress, induced by a wide range of known transcription inhibitors. We managed to show an increased variability of anti-sense transcripts in alpha-amanitin treated cells by applying the low-molecular RNA fraction extracted from agarose gel to PAGE-northern. Also, the fractioning of RNA by size using agarose gel slices occurred, being applicable for determining the sizes of target transcripts via RT-PCR.

Chromosomal Translocations in NK-Cell Lymphomas Originate from Inter-Chromosomal Contacts of Active rDNA Clusters Possessing Hot Spots of DSBs

Simple Summary

There are nine DSB hot spots located in the non-transcribed spacer of human rDNA units. Circular chromosome conformation capture data indicate that the rDNA clusters often shape contact with a specific set of chromosomal regions containing genes controlling differentiation and cancer, and often possessing the DSB hot spots. The data suggest a mechanism for rDNA-mediated translocation, and some of them could lead to tumorigenesis. Here, we searched for translocations in which rDNA clusters are involved. WGS data of normal T cells and NK-cell lymphomas from the same individuals were used. We revealed numerous translocations in which rDNA units are involved. The sites of these translocations in normal T cells and in the lymphomas were mostly different, but occurred at about the same frequency in both cell types. We conclude that oncogenic translocations lead to dysregulation of a specific set of genes controlling development.

Abstract

Endogenous hot spots of DNA double-strand breaks (DSBs) are tightly linked with transcription patterns and cancer. There are nine hot spots of DSBs (denoted Pleiades) in human rDNA units that are located exclusively inside the intergenic spacer (IGS). Profiles of Pleiades coincide with the profiles of γ-H2AX, suggesting a high level of in vivo breakage inside rDNA genes. The data were confirmed by microscopic observation of the largest γ-H2AX foci inside nucleoli in interphase chromosomes. Circular chromosome conformation capture (4C) data indicate that the rDNA units often make contact with a specific set of chromosomal regions containing genes that are involved in differentiation and cancer. Interestingly, these regions also often possess hot spots of DSBs that provide the potential for Robertsonian and oncogenic translocations. In this study, we searched for translocations in which rDNA clusters are involved. The whole genome sequence (WGS) data of normal T cells and NK-cell lymphomas from the same individuals revealed numerous translocations in which Pleiades were involved. The sites of these translocations in normal T cells and in the lymphomas were mostly different, although there were also some common sites. The genes at translocations in normal cells and in lymphomas are associated with predominantly non-overlapping lists of genes that are depleted with silenced genes. Our data indicate that rDNA-mediated translocations occur at about the same frequency in the normal T cells and NK-lymphoma cells but differ at particular sites that correspond to open chromatin. We conclude that oncogenic translocations lead to dysregulation of a specific set of genes controlling development. In normal T cells and in NK cells, there are hot spots of translocations at sites possessing strong H3K27ac marks. The data indicate that Pleiades are involved in rDNA-mediated translocation.

Interchromosomal Contacts of rDNA Clusters in Three Human Cell Lines Are Associated with Silencing of Genes Controlling Morphogenesis

To study the rDNA contacts with genes in three human cell lines of different origin, we used 4C approach. Our data indicate that the same set of about five hundred genes frequently shape contacts with rDNA clusters in HEK293T, K652, and hESM01 cells. Gene ontology search suggests that the genes are involved in development and morphogenesis. Approximately one hundred of these genes are highly associated with silencing by H3K27me3 mark in different normal cells, including bronchial epithelial cells, keratinocytes, myoblasts, monocytes, endothelial cells, kidney epithelial cells, and some others. We conclude that the concerted silencing of specific group of rDNA-contacting genes controlling development occurs during differentiation. We assume that the phase separation mechanisms may be involved in the rDNA-mediated silencing of a set of genes via the contacts with inactive rDNA clusters.

Dynamics of Whole-Genome Contacts of Nucleoli in Drosophila Cells Suggests a Role for rDNA Genes in Global Epigenetic Regulation

Chromosomes are organized into 3D structures that are important for the regulation of gene expression and differentiation. Important role in formation of inter-chromosome contacts play rDNA clusters that make up nucleoli. In the course of differentiation, heterochromatization of rDNA units in mouse cells is coupled with the repression or activation of different genes. Furthermore, the nucleoli of human cells shape the direct contacts with genes that are involved in differentiation and cancer. Here, we identified and categorized the genes located in the regions where rDNA clusters make frequent contacts. Using a 4C approach, we demonstrate that in Drosophila S2 cells, rDNA clusters form contacts with genes that are involved in chromosome organization and differentiation. Heat shock treatment induces changes in the contacts between nucleoli and hundreds of genes controlling morphogenesis. We show that nucleoli form contacts with regions that are enriched with active or repressive histone marks and where small non-coding RNAs are mapped. These data indicate that rDNA contacts are involved in the repression and activation of gene expression and that rDNA clusters orchestrate large groups of Drosophila genes involved in differentiation.

Copy number variations of satellite III (1q12) and ribosomal repeats in health and schizophrenia

OBJECTIVE

Earlier we studied the copy number variations (CNVs) of ribosomal repeat (rDNA) and the satellite III fragment (1q12) (f-SatIII) in the cells of schizophrenia patients (SZ) and healthy controls (HC). In the present study we pursued two main objectives: (1) to confirm the increased rDNA and decreased f-SatIII content in the genomes of enlarged SZ and HC samples and (2) to compare the rDNA and f-SatIII content in the same DNA samples of SZ and HC individuals.

METHODS

We determined the rDNA CN and f-SatIII content in the genomes of leukocytes of 1770 subjects [HC (N = 814) and SZ (N = 956)]. Non-radioactive quantitative hybridization method (NQH) was applied for analysis of the various combinations of the two repeats sizes in SZ and HC groups.

RESULTS

f-SatIII in human leukocytes (N = 1556) varies between 5.7 and 44.7 pg/ng DNA. RDNA CN varies between 200 and 896 (N = 1770). SZ group significantly differ from the HC group by lower f-SatIII content and by rDNA abundance. The f-SatIII and rDNA CN are not randomly combined in the genome. Higher rDNA CN values are associated with higher f-SatIII index values in SZ and HC. The f-SatIII variation interval in SZ group increases significantly in the subgroup with the high rDNA CN index values (>300 copies).

CONCLUSION

Schizophrenia patients' genomes contain low number of f-SatIII copies corresponding with a large ribosomal repeats CN. A scheme is proposed to explain the low f-SatIII content in SZ group against the background of high rDNA CN.

1Q12 Loci Movement in the Interphase Nucleus Under the Action of ROS Is an Important Component of the Mechanism That Determines Copy Number Variation of Satellite III (1q12) in Health and Schizophrenia

Introduction: Genome repeat cluster sizes can affect the chromatin spatial configuration and function. Low-dose ionizing radiation (IR) induces an adaptive response (AR) in human cells. AR includes the change in chromatin spatial configuration that is necessary to change the expression profile of the genome in response to stress. The 1q12 heterochromatin loci movement from the periphery to the center of the nucleus is a marker of the chromatin configuration change. We hypothesized that a large 1q12 domain could affect chromatin movement, thereby inhibiting the AR. Materials and Methods: 2D fluorescent in situ hybridization (FISH) method was used for the satellite III fragment from the 1q12 region (f-SatIII) localization analysis in the interphase nuclei of healthy control (HC) lymphocytes, schizophrenia (SZ) patients, and in cultured mesenchymal stem cells (MSCs). The localization of the nucleolus was analyzed by the nucleolus Ag staining. The non-radioactive quantitative hybridization (NQH) technique was used for the f-SatIII fragment content in DNA analysis. Satellite III fragments transcription was analyzed by reverse transcriptase quantitative PCR (RT-qPCR). Results: Low-dose IR induces the small-area 1q12 domains movement from the periphery to the central regions of the nucleus in HC lymphocytes and MSCs. Simultaneously, nucleolus moves from the nucleus center toward the nuclear envelope. The nucleolus in that period increases. The distance between the 1q12 domain and the nucleolus in irradiated cells is significantly reduced. The large-area 1q12 domains do not move in response to stress. During prolonged cultivation, the irradiated cells with a large f-SatIII amount die, and the population is enriched with the cells with low f-SatIII content. IR induces satellite III transcription in HC lymphocytes. Intact SZ patients' lymphocytes have the same signs of nuclei activation as irradiated HC cells. Conclusion: When a cell population responds to stress, cells are selected according to the size of the 1q12 domain (the f-SatIII content). The low content of the f-SatIII repeat in SZ patients may be a consequence of the chronic oxidative stress and of a large copies number of the ribosomal repeats.