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.

[Drosophila rDNA Genes Shape the Stable Contacts with the Tlk Gene at the Expression Area of Small RNAs and Affect on Looped Domains inside the Gene].

In experiments on mouse and human cells it was demonstrated that rDNA plays an important role in epigenetic regulation of many genes. To identify and study rDNA-contacting genes in Drosophila we used the 4С (circular chromosome conformation capture) approach. We detected very stable contacts of rDNA genes within a 5-kb region inside the Tlk gene residing in X chromosome. This 5-kb region corresponds to small RNAs. After heat shock treatment both the amount of contacts, and the expression level of the gene were increased. Tlk and Rala are genes that share the same short bidirectional promoter but exhibit different expression levels. Around the region of rDNA contacts inside the Tlk gene, looped domains were formed. We conclude that rDNA contact-dependent epigenetic regulation is guided by small RNAs and that the contacts are involved in rearrangements of the looped domains.

[Contact Sites of rDNA Clusters with FANK1 Gene Correspond to Repressed Chromatin].

rDNA genes play an important role in epigenetic regulation and in differentiation of eukaryotic cells. Using the 4C (circular chromosome conformation capture) approach and model HEK293T cells, we analyzed the rDNA-contacting gene FANK1, using anchor located inside rDNA genes. At the 5' end of the FANK1 gene we detected frequent contacts with rDNA clusters. The contact sites coincide with the border where chromatin state changes and nucleosome positioning. The adjacent genes DHX32, BCCIP and UROS are located in the active chromatin and are transcribed, but do not contact with rDNA genes, while FANK1 gene is silenced, and is located in repressed chromatin. Heat shock treatment dramatically changes the pattern of rDNA contacts in the region and induces about 4-fold increase in activation of the FANK1 gene. We conclude that rDNA contacts may be involved in repression of the FANK1 gene.

Interchromosomal Contacts of rDNA Clusters with DUX Genes in Human Chromosome 4 Are Very Sensitive to Heat Shock Treatment.

In order to study the effects of heat shock treatment on the distribution of rDNA contacts at the region possessing DUX genes inside chromosome 4 we used 4C approach. Our data indicate that the treatment removes the frequent rDNA contacts in this region. The recent data on involvement of superenhancers that are decorated by broad H3K27ac marks in the phase separation mechanisms and the previous data demonstrating that these broad marks are the favorite sites of rDNA contacts taken together with our data on sensitivity of the contacts to the heat shock treatment suggest that the phase separation mechanisms are involved in the reversible rDNA-mediated regulation of gene expression via the contacts.

[Homeotic DUX4 Genes that Control Human Embryonic Development at the Two-Cell Stage Are Surrounded by Regions Contacting with rDNA Gene Clusters].

Many human genes that control human embryonic development and differentiation of human cells form chromosomal contact with rRNA gene clusters, which are involved in the epigenetic regulation of many genes. The sites of rRNA gene contact often fall on extended (up to 50 kb) regions containing a chromatin mark, H3K27ac histone, typical for superenhancers, as well as on pericentromeric and subtelomeric regions of chromosomes. We found that the DUX4 genes located in the subtelomeric region of human chromosome 4 are surrounded by regions that are often in contact with the rRNA genes. The 25 kb region of this chromosome, presented in version hg19 of the sequenced human genome, contains several copies of the DUX4 gene. The sites of rRNA gene contacts located around this region contain methylation sites as well as CTCF binding sites. It is assumed that the rRNA gene contacts are important in silencing these DUX4 gene copies.

[Mutation Frequencies in RNAi Targets in HIV-1 Genomes Obtained from Blood Plasma of Patients Receiving Anti-Retroviral Therapy].

Gene therapy for AIDS based on RNA interference (RNAi) is currently looked upon as a promising alternative to conventional antiretroviral chemotherapy. The high variability of HIV-1 is the main challenge in developing new approaches to AIDS therapy. To date, about 18 million HIV-1 infected individuals receive antiretroviral therapy worldwide. As of 2017, about 44% of individuals with AIDS received antiretroviral therapy in Russia. Since the RNAs used for efficient RNAi and the corresponding targets in the viral transcript should be perfectly complementary to each other, it is necessary to continuously monitor the nucleotide sequences of clinical HIV-1 isolates obtained from blood and cells of naïve patients and patients receiving antiretroviral therapy. Comprehensive analysis of the mutation frequencies in the viral genome is only possible with deep sequencing approaches. The present paper reports on an analysis of the mutation frequencies in six 100 bp genome regions in clinical HIV-1 isolates obtained from blood plasma of four Russian AIDS patients who have been receiving antiretroviral therapy for several years. These regions contain efficient RNAi targets. The average frequencies of all possible transversions and transitions within the RNAi targets and in their proximity have been estimated. It has been demonstrated that reverse transcriptase inhibition decreases the frequency of a number of reverse mutations. It has been found that mutations in RNAi targets are rarer (5-75 times lower than the mutation frequency for different nucleotide substitutions) than in the adjacent sequences. Our findings speak in favor of these conservative targets for developing new approaches to gene therapy of AIDS.

[Mutation Frequencies in HIV-1 Genome in Regions Containing Efficient RNAi Targets As Calculated from Ultra-Deep Sequencing Data].

HIV-1 is one of the most variable viruses. The development of gene therapy technology using RNAi for AIDS/HIV-1 treatment is a potential alternative for traditional anti-retroviral therapy. Anti-HIV-1 siRNA should aim to exploit the most conserved viral targets. Using the deep sequencing of potential RNAi targets in 100-nt HIV-1 genome fragments from the clinical HIV-1 subtype A isolates in Russia, we found that the frequencies of all possible transversions and transitions in certain RNAi targets are 3-38 times lower than in adjacent sequences. Therefore, these targets are conserved. We propose the development of these RNAi targets for AIDS/HIV-1 treatment. Deep sequencing also enables the detection of the characteristic mutational bias of RT during the replication of viral RNA.

[Mutation frequencies in HIV-1 subtype-A genome in regions containing efficient RNAi targets].

The development of gene-therapy technology using RNAi for AIDS/HIV-1 treatment is a prospective alternative to traditional anti-retroviral therapy. RNAi targets could be selected in HIV-1 transcripts and in CCR5 mRNA. Previously, we experimentally selected a number of efficient siRNAs that target HIV-1 RNAs. The viral genome mutates frequently, and RNAi strength is very sensitive, even for a single mismatches. That is why it is important to study nucleotide sequences of targets in clinical isolates of HIV-1. In the present study, we analyzed mutations in 6 of about 300-bp regions containing RNAi targets from HIV-1 subtype A isolates in Russia. Estimates of the mean frequencies of mutations in the targets were obtained and the frequencies of mutations in the different codon positions were compared. The frequencies of mutations in the vicinity of the targets and directly within the targets were also compared and have been shown to be approximately the same. The frequencies of indels in the chosen regions have been assessed. Their frequencies have proved to be two to three orders of magnitude less compared to that for mutations.

Mapping of genomic double-strand breaks by ligation of biotinylated oligonucleotides to forum domains: Analysis of the data obtained for human rDNA units.

DNA double-strand breaks (DSBs) are associated with different physiological and pathological processes in different organisms. To understand the role of DSBs in multiple cellular mechanisms, a robust method for genome-wide mapping of chromosomal breaks at one-nucleotide resolution is required. Many years ago, we detected large DNA fragments migrating from DNA-agarose plugs in pulsed-field gels, which we named 'forum domains' [1,2]. Recently, we developed a method for genome-wide mapping of DSBs that produces these 50-150 kb DNA domains using microarrays or 454 sequencing (Tchurikov et al., 2011; 2013). Now we have used Illumina sequencing to map DSBs in repetitive rDNA units in human HEK293T cells. Here we describe in detail the experimental design and bioinformatics analysis of the data deposited in the Gene Expression Omnibus with accession number GSE49302 and associated with the study published in the Journal of Molecular Cell Biology (Tchurikov et al., 2014).

Genome-wide mapping of hot spots of DNA double-strand breaks in human cells as a tool for epigenetic studies and cancer genomics.

Hot spots of DNA double-strand breaks (DSBs) are associated with coordinated expression of genes in chromosomal domains (Tchurikov et al., 2011 [1]; 2013). These 50-150-kb DNA domains (denoted "forum domains") can be visualized by separation of undigested chromosomal DNA in pulsed-field agarose gels (Tchurikov et al., 1988; 1992) and used for genome-wide mapping of the DSBs that produce them. Recently, we described nine hot spots of DSBs in human rDNA genes and observed that, in rDNA units, the hot spots coincide with CTCF binding sites and H3K4me3 marks (Tchurikov et al., 2014), suggesting a role for DSBs in active transcription. Here we have used Illumina sequencing to map DSBs in chromosomes of human HEK293T cells, and describe in detail the experimental design and bioinformatics analysis of the data deposited in the Gene Expression Omnibus with accession number GSE53811 and associated with the study published in DNA Research (Kravatsky et al., 2015). Our data indicate that H3K4me3 marks often coincide with hot spots of DSBs in HEK293T cells and that the mapping of these hot spots is important for cancer genomic studies.

Molecular analysis of transcription start sites of separate copies of the suffix short retroelement in genome of Drosophila.

In order to study TSS in the suffix element, we used total RNA isolated from ovaries of Drosophila melanogaster. Using a 5'-RACE System (Invitrogen) and 454 sequencing, we found the full-length suffix sense transcripts. However, most 5'-RACE reads (>70%) correspond to 5'-truncated transcripts lacking the first 33-39 nucleotides. The data may indicate that these RNAs are products of the RNAi-related silencing mech-anism producing small RNAs that are larger than piRNAs. The full-length suffix transcripts could arise either from the F element or from active separate copies of suffix.