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.

Genome-Wide Study of Colocalization between Genomic Stretches: A Method and Applications to the Regulation of Gene Expression.

In this paper, we describe a method for the study of colocalization effects between stretch-stretch and stretch-point genome tracks based on a set of indices varying within the (-1, +1) interval. The indices combine the distances between the centers of neighboring stretches and their lengths. The extreme boundaries of the interval correspond to the complete colocalization of the genome tracks or its complete absence. We also obtained the relevant criteria of statistical significance for such indices using the complete permutation test. The method is robust with respect to strongly inhomogeneous positioning and length distribution of the genome tracks. On the basis of this approach, we created command-line software, the Genome Track Colocalization Analyzer. The software was tested, compared with other available packages, and applied to particular problems related to gene expression. The package, Genome Track Colocalization Analyzer (GTCA), is freely available to the users. GTCA complements our previous software, the Genome Track Analyzer, intended for the search for pairwise correlations between point-like genome tracks (also freely available). The corresponding details are provided in Data Availability Statement at the end of the text.

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.

Ribonucleocapsid assembly/packaging signals in the genomes of the coronaviruses SARS-CoV and SARS-CoV-2: detection, comparison and implications for therapeutic targeting.

The genomic ssRNA of coronaviruses is packaged within a helical nucleocapsid. Due to transitional symmetry of a helix, weakly specific cooperative interaction between ssRNA and nucleocapsid proteins leads to the natural selection of specific quasi-periodic assembly/packaging signals in the related genomic sequence. Such signals coordinated with the nucleocapsid helical structure were detected and reconstructed in the genomes of the coronaviruses SARS-CoV and SARS-CoV-2. The main period of the signals for both viruses was about 54 nt, that implies 6.75 nt per N protein. The complete coverage of the ssRNA genome of length about 30,000 nt by the nucleocapsid would need 4.4 × 103 N proteins, that makes them the most abundant among the structural proteins. The repertoires of motifs for SARS-CoV and SARS-CoV-2 were divergent but nearly coincided for different isolates of SARS-CoV-2. We obtained the distributions of assembly/packaging signals over the genomes with nonoverlapping windows of width 432 nt. Finally, using the spectral entropy, we compared the load from point mutations and indels during virus age for SARS-CoV and SARS-CoV-2. We found the higher mutational load on SARS-CoV. In this sense, SARS-CoV-2 can be treated as a 'newborn' virus. These observations may be helpful in practical medical applications and are of basic interest. Communicated by Ramaswamy H. Sarma.

Evolving ribonucleocapsid assembly/packaging signals in the genomes of the human and animal coronaviruses: targeting, transmission and evolution.

A world-wide COVID-19 pandemic intensified strongly the studies of molecular mechanisms related to the coronaviruses. The origin of coronaviruses and the risks of human-to-human, animal-to-human and human-to-animal transmission of coronaviral infections can be understood only on a broader evolutionary level by detailed comparative studies. In this paper, we studied ribonucleocapsid assembly-packaging signals (RNAPS) in the genomes of all seven known pathogenic human coronaviruses, SARS-CoV, SARS-CoV-2, MERS-CoV, HCoV-OC43, HCoV-HKU1, HCoV-229E and HCoV-NL63 and compared them with RNAPS in the genomes of the related animal coronaviruses including SARS-Bat-CoV, MERS-Camel-CoV, MHV, Bat-CoV MOP1, TGEV and one of camel alphacoronaviruses. RNAPS in the genomes of coronaviruses were evolved due to weakly specific interactions between genomic RNA and N proteins in helical nucleocapsids. Combining transitional genome mapping and Jaccard correlation coefficients allows us to perform the analysis directly in terms of underlying motifs distributed over the genome. In all coronaviruses, RNAPS were distributed quasi-periodically over the genome with the period about 54 nt biased to 57 nt and to 51 nt for the genomes longer and shorter than that of SARS-CoV, respectively. The comparison with the experimentally verified packaging signals for MERS-CoV, MHV and TGEV proved that the distribution of particular motifs is strongly correlated with the packaging signals. We also found that many motifs were highly conserved in both characters and positioning on the genomes throughout the lineages that make them promising therapeutic targets. The mechanisms of encapsidation can affect the recombination and co-infection as well.Communicated by Ramaswamy H. Sarma.

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.

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

Human rDNA clusters form numerous contacts with different chromosomal regions as evidenced by chromosome conformation capture data. Heterochromatization of rDNA genes leads to heterochromatization in different chromosomal regions coupled with the activation of the transcription of genes related to differentiation. These data suggest a role for rDNA clusters in the regulation of many human genes. However, the genes that reside within the rDNA-contacting regions have not been identified. The purpose of this study was to detect and characterize the regions where rDNA clusters make frequent contacts and to identify and categorize genes located in these regions. We analyzed the regions that contact rDNA using 4C data and show that these regions are enriched with genes specifying transcription factors and non-coding RNAs involved in differentiation and development. The rDNA-contacting genes are involved in neuronal development and are associated with different cancers. Heat shock treatment led to dramatic changes in the pattern of rDNA-contacting sites, especially in the regions possessing long stretches of H3K27ac marks. Whole-genome analysis of rDNA-contacting sites revealed specific epigenetic marks and the transcription sites of 20-100 nt non-coding RNAs in these regions. The rDNA-contacting genes jointly regulate many genes that are involved in the control of transcription by RNA polymerase II and the development of neurons. Our data suggest a role for rDNA clusters in the differentiation of human cells and carcinogenesis.

Six Highly Conserved Targets of RNAi Revealed in HIV-1-Infected Patients from Russia Are Also Present in Many HIV-1 Strains Worldwide.

RNAi has been suggested for use in gene therapy of HIV/AIDS, but the main problem is that HIV-1 is highly variable and could escape attack from the small interfering RNAs (siRNAs) due to even single nucleotide substitutions in the potential targets. To exhaustively check the variability in selected RNA targets of HIV-1, we used ultra-deep sequencing of six regions of HIV-1 from the plasma of two independent cohorts of patients from Russia. Six RNAi targets were found that are invariable in 82%-97% of viruses in both cohorts and are located inside the domains specifying reverse transcriptase (RT), integrase, vpu, gp120, and p17. The analysis of mutation frequencies and their characteristics inside the targets suggests a likely role for APOBEC3G (apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like 3G, A3G) in G-to-A mutations and a predominant effect of RT biases in the detected variability of the virus. The lowest frequency of mutations was detected in the central part of all six targets. We also discovered that the identical RNAi targets are present in many HIV-1 strains from many countries and from all continents. The data are important for both the understanding of the patterns of HIV-1 mutability and properties of RT and for the development of gene therapy approaches using RNAi for the treatment of HIV/AIDS.

Analysis of Variability in HIV-1 Subtype A Strains in Russia Suggests a Combination of Deep Sequencing and Multitarget RNA Interference for Silencing of the Virus.

Any method for silencing the activity of the HIV-1 retrovirus should tackle the extremely high variability of HIV-1 sequences and mutational escape. We studied sequence variability in the vicinity of selected RNA interference (RNAi) targets from isolates of HIV-1 subtype A in Russia, and we propose that using artificial RNAi is a potential alternative to traditional antiretroviral therapy. We prove that using multiple RNAi targets overcomes the variability in HIV-1 isolates. The optimal number of targets critically depends on the conservation of the target sequences. The total number of targets that are conserved with a probability of 0.7-0.8 should exceed at least 2. Combining deep sequencing and multitarget RNAi may provide an efficient approach to cure HIV/AIDS.

Conserved sequences in the current strains of HIV-1 subtype A in Russia are effectively targeted by artificial RNAi in vitro.

Highly active antiretroviral therapy has greatly reduced the morbidity and mortality of AIDS. However, many of the antiretroviral drugs are toxic with long-term use, and all currently used anti-HIV agents generate drug-resistant mutants. Therefore, there is a great need for new approaches to AIDS therapy. RNAi is a powerful means of inhibiting HIV-1 production in human cells. We propose to use RNAi for gene therapy of HIV/AIDS. Previously we identified a number of new biologically active siRNAs targeting several moderately conserved regions in HIV-1 transcripts. Here we analyze the heterogeneity of nucleotide sequences in three RNAi targets in sequences encoding the reverse transcriptase and integrase domains of current isolates of HIV-1 subtype A in Russia. These data were used to generate genetic constructs expressing short hairpin RNAs 28-30-bp in length that could be processed in cells into siRNAs. After transfection of the constructs we observed siRNAs that efficiently attacked the selected targets. We expect that targeting several viral genes important for HIV-1 reproduction will help overcome the problem of viral adaptation and will prevent the appearance of RNAi escape mutants in current virus strains, an important feature of gene therapy of HIV/AIDS.

Multiplex determination of serological signatures in the sera of colorectal cancer patients using hydrogel biochips.

Colorectal cancer (CRC) is the third most common malignancy in industrialized countries. Despite the advances in diagnostics and development of new drugs, the 5-year survival remains only 60-65%. Our approach to early diagnostics of CRC is based on the determination of serological signatures with an array of hemispherical hydrogel cells containing immobilized proteins and oligosaccharides (glycochip). The compounds immobilized on the glycochip include tumor-associated glycans (SiaTn, Tn, TF, Le(C) , Le(Y) , SiaLe(A) , and Manß1-4GlcNAcß) and antibodies against human immunoglobulins IgG, IgA, and IgM. The glycochip detects antibodies against tumor-associated glycans in patients' sera. The simultaneous measurement of the levels of immunoglobulins enhances the diagnostic impact of the signatures. In this work, we found previously unreported increase in antibodies against oligosaccharide Manß1-4GlcNAcß in patients with CRC. In parallel with these experiments, we determined the levels of oncomarkers carcinoembryonic antigen (CEA), cancer antigen (CA) 19-9, CA 125, CA 15-3, human chorionic gonadotropin (HCG), and alpha-fetoprotein (AFP) using another gel-based biochip with immobilized antibodies (oncochip) developed earlier in our laboratory. In total, 69 samples from healthy donors, 33 from patients with colorectal carcinoma, and 27 from patients with inflammatory bowel diseases were studied. The use of combined signatures of antiglycan antibodies and oncomarkers provides much better predictive value than the conventional measurement of oncomarkers CEA and CA 19-9. Positive predictive value of CRC diagnoses using together glycochip and oncochip reached 95% with the sensitivity and specificity 88% and 98%, respectively. Thus, the combination of antibody profiling with detection of conventional oncomarkers proved to be a promising tool in diagnostics of CRC.

Genome-wide study of correlations between genomic features and their relationship with the regulation of gene expression.

The broad class of tasks in genetics and epigenetics can be reduced to the study of various features that are distributed over the genome (genome tracks). The rapid and efficient processing of the huge amount of data stored in the genome-scale databases cannot be achieved without the software packages based on the analytical criteria. However, strong inhomogeneity of genome tracks hampers the development of relevant statistics. We developed the criteria for the assessment of genome track inhomogeneity and correlations between two genome tracks. We also developed a software package, Genome Track Analyzer, based on this theory. The theory and software were tested on simulated data and were applied to the study of correlations between CpG islands and transcription start sites in the Homo sapiens genome, between profiles of protein-binding sites in chromosomes of Drosophila melanogaster, and between DNA double-strand breaks and histone marks in the H. sapiens genome. Significant correlations between transcription start sites on the forward and the reverse strands were observed in genomes of D. melanogaster, Caenorhabditis elegans, Mus musculus, H. sapiens, and Danio rerio. The observed correlations may be related to the regulation of gene expression in eukaryotes. Genome Track Analyzer is freely available at http://ancorr.eimb.ru/.

Computational methods of identification of pseudogenes based on functionality: entropy and GC content.

Spectral entropy and GC content analyses reveal comprehensive structural features of DNA sequences. To illustrate the significance of these features, we analyze the ß-esterase gene cluster, including the Est-6 gene and the ψEst-6 putative pseudogene, in seven species of the Drosophila melanogaster subgroup. The spectral entropies show distinctly lower structural ordering for ψEst-6 than for Est-6 in all species studied. However, entropy accumulation is not a completely random process for either gene and it shows to be nucleotide dependent. Furthermore, GC content in synonymous positions is uniformly higher in Est-6 than in ψEst-6, in agreement with the reduced GC content generally observed in pseudogenes and nonfunctional sequences. The observed differences in entropy and GC content reflect an evolutionary shift associated with the process of pseudogenization and subsequent functional divergence of ψEst-6 and Est-6 after the duplication event. The data obtained show the relevance and significance of entropy and GC content analyses for pseudogene identification and for the comparative study of gene-pseudogene evolution.

Structural attributes of nucleotide sequences in promoter regions of supercoiling-sensitive genes: how to relate microarray expression data with genomic sequences.

The level of supercoiling in the chromosome can affect gene expression. To clarify the basis of supercoiling sensitivity, we analyzed the structural features of nucleotide sequences in the vicinity of promoters for the genes with expression enhanced and decreased in response to loss of chromosomal supercoiling in Escherichia coli. Fourier analysis of promoter sequences for supercoiling-sensitive genes reveals the tendency in selection of sequences with helical periodicities close to 10nt for relaxation-induced genes and to 11nt for relaxation-repressed genes. The helical periodicities in the subsets of promoters recognized by RNA polymerase with different sigma factors were also studied. A special procedure was developed for the study of correlations between the intensities of periodicities in promoter sequences and the expression levels of corresponding genes. Significant correlations of expression with the AT content and with AT periodicities about 10, 11, and 50nt indicate their role in regulation of supercoiling-sensitive genes.

Separate production of single-stranded DNA is not necessary: circuit denaturation of double-stranded DNA followed by hybridization of single strands on oligonucleotide microchips.

An approach to circuit renaturation-hybridization of dsDNA on oligonucleotide microchips is described. A close circuit cycling device has been developed, and the feasibility of the proposed technique was demonstrated on two platforms. First, a commercial microchip for detection of rifampicin resistance in Mycobacterium tuberculosis was used. Hybridization of a 126 nt long single-stranded DNA (ssDNA) fragment of the rpoB gene according to manufacturer's protocol has been compared to hybridization of the same double-stranded DNA (dsDNA) fragment using the developed approach. Hybridization signals obtained by both methods were comparable in intensity and correlated closely. Second, a 22 nt long hairpin-forming oligonucleotide was designed and hybridized with a custom microchip containing probes complementary to both strands of the oligonucleotide. Conventional hybridization of this oligonucleotide did not yield any significant signals. Cleavage of the hairpin loop resulted in the formation of a 9 bp long intermolecular duplex. Hybridization of the duplex using the suggested technique yielded strong signals. The proposed approach allows analyzing target DNA in double-stranded form bypassing the preparation of single-stranded targets. Moreover, both complementary chains could be analyzed simultaneously, providing a reliable internal control. Being combined with fragmentation this method opens new possibilities in analyzing ssDNA with complex secondary structure.

Gel-based oligonucleotide microarray approach to analyze protein-ssDNA binding specificity.

Gel-based oligonucleotide microarray approach was developed for quantitative profiling of binding affinity of a protein to single-stranded DNA (ssDNA). To demonstrate additional capabilities of this method, we analyzed the binding specificity of ribonuclease (RNase) binase from Bacillus intermedius (EC 3.1.27.3) to ssDNA using generic hexamer oligodeoxyribonucleotide microchip. Single-stranded octamer oligonucleotides were immobilized within 3D hemispherical gel pads. The octanucleotides in individual pads 5'-{N}N(1)N(2)N(3)N(4)N(5)N(6){N}-3' consisted of a fixed hexamer motif N(1)N(2)N(3)N(4)N(5)N(6) in the middle and variable parts {N} at the ends, where {N} represent A, C, G and T in equal proportions. The chip has 4096 pads with a complete set of hexamer sequences. The affinity was determined by measuring dissociation of the RNase-ssDNA complexes with the temperature increasing from 0 degrees C to 50 degrees C in quasi-equilibrium conditions. RNase binase showed the highest sequence-specificity of binding to motifs 5'-NNG(A/T/C)GNN-3' with the order of preference: GAG > GTG > GCG. High specificity towards G(A/T/C)G triplets was also confirmed by measuring fluorescent anisotropy of complexes of binase with selected oligodeoxyribonucleotides in solution. The affinity of RNase binase to other 3-nt sequences was also ranked. These results demonstrate the applicability of the method and provide the ground for further investigations of nonenzymatic functions of RNases.

Entropy and GC Content in the beta-esterase gene cluster of the Drosophila melanogaster subgroup.

We perform spectral entropy and GC content analyses in the beta-esterase gene cluster, including the Est-6 gene and the psiEst-6 putative pseudogene, in seven species of the Drosophila melanogaster species subgroup. psiEst-6 combines features of functional and nonfunctional genes. The spectral entropies show distinctly lower structural ordering for psiEst-6 than for Est-6 in all species studied. Our observations agree with previous results for D. melanogaster and provide additional support to our hypothesis that after the duplication event Est-6 retained the esterase-coding function and its role during copulation, while psiEst-6 lost that function but now operates in conjunction with Est-6 as an intergene. Entropy accumulation is not a completely random process for either gene. Structural entropy is nucleotide dependent. The relative normalized deviations for structural entropy are higher for G than for C nucleotides. The entropy values are similar for Est-6 and psiEst-6 in the case of A and T but are lower for Est-6 in the case of G and C. The GC content in synonymous positions is uniformly higher in Est-6 than in psiEst-6, which agrees with the reduced GC content generally observed in pseudogenes and nonfunctional sequences. The observed differences in entropy and GC content reflect an evolutionary shift associated with the process of pseudogenization and subsequent functional divergence of psiEst-6 and Est-6 after the duplication event.