CBP and RAD21 Proteins Bind at the Termini of Forum Domains in Human Chromosomes

Forum domains are 50-100-kb stretches of DNA delimited by the hotspots of double-strand breaks (DSBs). These domains possess coordinately expressed genes. However, molecular mechanisms of such regulation are not clear. It is assumed that the proteins specifically binding at the termini of domains can be involved in coordinated regulation of expression. In this study, we used the results of precise mapping of hotspots of DSBs and ChIP-Seq data for ten nuclear proteins in HEK293T cell line for a search of proteins specifically binding at forum-domain termini. We detected that two proteins, CBP and RAD24, which are known to be involved in epigenetic regulation of gene expression and formation of 3D chromosomal structures, bind at the termini. We assume that these proteins may be involved in coordinated expression of genes in forum domains.

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.

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.

A giant left internal iliac artery aneurysm in a patient with Loeys-Dietz syndrome

A 16-year-old female presented with left iliac fossa pain. In January 2021, she was admitted to her local hospital with severe lower abdominal pain and the pelvic ultrasound demonstrated a 13-cm left internal iliac artery dissecting aneurysm with its partial thrombosis. On examination, she had a high-arched palate, multiple skin stretch marks, flat feet and a soft systolic ejection murmur at the left 5th mid-clavicular line. She had a mildly tender abdomen in the left iliac fossa. Computed tomography angiography demonstrated a 12.2 cm × 10.4 cm × 12.5 cm left internal iliac artery aneurysm. During surgery, the aneurysm was incised and the proximal and distal orifices of the internal iliac artery were ligated. Genetic testing yielded 2 mutations in the SMAD3 gene characteristic for Loeys-Dietz syndrome.

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.

[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.

Individual effects of the copia and gypsy enhancer and insulator on chromatin marks, eRNA synthesis, and binding of insulator proteins in transfected genetic constructs

Enhancers and insulators are involved in the regulation of gene expression, but the basic underlying mechanisms of action of these elements are unknown. We analyzed the individual effects of the enhancer and the insulator from Drosophila mobile elements copia [enh(copia)] and gypsy using transfected genetic constructs in S2 cells. This system excludes the influence of genomic cis regulatory elements. The enhancer-induced synthesis of 350-1050-nt-long enhancer RNAs (eRNAs) and H3K4me3 and H3K18ac marks, mainly in the region located about 300bp downstream of the enhancer. Insertion of the insulator between the enhancer and the promoter reduced these effects. We also observed the binding of dCTCF to the enhancer and to gypsy insulator. Our data indicate that a single gypsy insulator interacts with both the enhancer and the promoter, while two copies of the gypsy insulator preferentially interact with each other. Our results suggest the formation of chromatin loops that are shaped by the enhancer and the insulator.

[Molecular Mechanisms of HIV-1 Genetic Diversity]

High genetic diversity of HIV-1 is the main factor behind the fact that HIV infection is widespread and difficult to treat. Although a limited number (or only one) of virus particles enters the blood upon infection, the particles are replicated in infected cells and rapidly produce new genetic variants that are resistant to the host immune system and antiretroviral drugs. This circumstance hampers the development of anti-HIV-1 vaccines and requires new antiretroviral drugs to be designed. The cause of the high variation of HIV-1 is related to the properties of its reverse transcriptase, which is error prone and often makes mistakes when transcribing virus RNA. Moreover, host APOBEC3-family proteins deaminate cytosines in the resulting minus strand DNA copy, leading to C/G-T/A transitions. The review considers several mechanisms that generate HIV-1 variants, including multiple recombination events between two different RNA copies colocated within one capsid. To understand the mechanisms of high genetic diversity of HIV-1 is essential for designing basically new approaches to treatment of HIV infection and 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.

[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.

Hot spots of DNA double-strand breaks in human rDNA units are produced in vivo

Endogenous hot spots of DNA double-strand breaks (DSBs) are tightly linked with transcription patterns and cancer genomics(1,2). There are nine hot spots of DSBs located in human rDNA units(3-6). Here we describe that the profiles of these hot spots coincide with the profiles of γ-H2AX or H2AX, strongly 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. In metaphase chromosomes, we observed that only some portion of rDNA clusters possess γ-H2AX foci and that all γ-H2AX foci co-localize with UBF-1 binding sites, which strongly suggests that only active rDNA units possess the hot spots of DSBs. Both γ-H2AX and UBF-1 are epigenetically inherited and thus indicate the rDNA units that were active in the previous cell cycle. These results have implications for diverse fields, including epigenetics and cancer genomics.

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.

Hot spots of DNA double-strand breaks and genomic contacts of human rDNA units are involved in epigenetic regulation

DNA double-strand breaks (DSBs) are involved in many cellular mechanisms, including replication, transcription, and genome rearrangements. The recent observation that hot spots of DSBs in human chromosomes delimit DNA domains that possess coordinately expressed genes suggests a strong relationship between the organization of transcription patterns and hot spots of DSBs. In this study, we performed mapping of hot spots of DSBs in a human 43-kb ribosomal DNA (rDNA) repeated unit. We observed that rDNA units corresponded to the most fragile sites in human chromosomes and that these units possessed at least nine specific regions containing clusters of extremely frequently occurring DSBs, which were located exclusively in non-coding intergenic spacer (IGS) regions. The hot spots of DSBs corresponded to only a specific subset of DNase-hypersensitive sites, and coincided with CTCF, PARP1, and HNRNPA2B1 binding sites, and H3K4me3 marks. Our rDNA-4C data indicate that the regions of IGS containing the hot spots of DSBs often form contacts with specific regions in different chromosomes, including the pericentromeric regions, as well as regions that are characterized by H3K27ac and H3K4me3 marks, CTCF binding sites, ChIA-PET and RIP signals, and high levels of DSBs. The data suggest a strong link between chromosome breakage and several different mechanisms of epigenetic regulation of gene expression.

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.

DNA double-strand breaks coupled with PARP1 and HNRNPA2B1 binding sites flank coordinately expressed domains in human chromosomes

Genome instability plays a key role in multiple biological processes and diseases, including cancer. Genome-wide mapping of DNA double-strand breaks (DSBs) is important for understanding both chromosomal architecture and specific chromosomal regions at DSBs. We developed a method for precise genome-wide mapping of blunt-ended DSBs in human chromosomes, and observed non-random fragmentation and DSB hot spots. These hot spots are scattered along chromosomes and delimit protected 50-250 kb DNA domains. We found that about 30% of the domains (denoted forum domains) possess coordinately expressed genes and that PARP1 and HNRNPA2B1 specifically bind DNA sequences at the forum domain termini. Thus, our data suggest a novel type of gene regulation: a coordinated transcription or silencing of gene clusters delimited by DSB hot spots as well as PARP1 and HNRNPa2B1 binding sites.

Ribonuclease binase apoptotic signature in leukemic Kasumi-1 cells

Cytotoxic exogenous RNases triggering apoptotic response in malignant cells have potential as anticancer drugs; surprisingly, detailed characterization of the RNase-induced apoptosis has not been conducted so far. Here we show that a cytotoxic RNase from Bacillus intermedius (binase) induces extrinsic and intrinsic apoptotic pathways in leukemic Kasumi-1 cells. The experiments were performed using TaqMan Array Human Apoptosis 96-well Plate for gene expression analysis, and flow cytometry. Cytometric studies demonstrated dissipation of the mitochondrial membrane potential, opening of mitochondrial permeability transition pores, activation of caspases, increase of intracellular Ca(2+) and decrease of reactive oxygen species levels. We found that expression of 62 apoptotic genes is up-regulated, including 16 genes that are highly up-regulated, and only one gene was found to be down-regulated. The highest, 16 fold increase of the expression level was observed for TNF gene. Highly up-regulated genes also include the non-canonical NF-κB signaling pathway and inflammatory caspases 1,4. The obtained results suggest that binase induces evolutionary acquired cellular response to a microbial agent and triggers unusual apoptosis pathway.

Both piRNA and siRNA pathways are silencing transcripts of the suffix element in the Drosophila melanogaster germline and somatic cells

In the Drosophila melanogaster germline, the piRNA pathway silences retrotransposons as well as other transcribed repetitive elements. Suffix is an unusual short retroelement that was identified both as an actively transcribed repetitive element and also as an element at the 3' ends of the Drosophila non-LTR F element. The copies of suffix that are F element-independent are far more actively transcribed than their counterparts on the F element. We studied the patterns of small RNAs targeting both strands of suffix in Drosophila ovaries using an RNase protection assay and the analysis of the corresponding RNA sequences from the libraries of total small RNAs. Our results indicate that suffix sense and antisense transcripts are targeted mainly by 23-29 nucleotides in length piRNAs and also by 21 nucleotides in length siRNAs. Suffix sense transcripts actively form longer RNA species, corresponding either to partial digestion products of the RNAi and Piwi pathways or to another RNA silencing mechanism. Both sense and antisense suffix transcripts accumulated in the ovaries of homozygous spn-E, piwi and aub mutants. These results provide evidence that suffix sense and antisense transcripts in the germ line and soma are targeted by both RNAi and Piwi pathways and that a Dicer-independent pathway of biogenesis of siRNAs could exist in Drosophila cells.

Oncogenic c-kit transcript is a target for binase

Mutational activation of c-Kit receptor tyrosine kinase is common in acute myelogenous leukemia (AML). One such activating point mutation is the N822K replacement in the c-Kit protein. Here we investigate the selective cytotoxic effect of binase--RNase from Bacillus intermedius--on FDC-P1-N822K cells. These cells were derived from myeloid progenitor FDC-P1 cells, in which ectopic expression of N822K c-kit gene induces interleukin-3 independent growth. In order to determine whether the sensitivity of these cells to binase is caused by the expression of c-kit oncogene, the cytotoxicity of the RNase was studied in the presence of selective inhibitor of mutated c-Kit imatinib (Gleevec). Inhibition of mutated c-Kit protein leads to the loss of cell sensitivity to the apoptotic effect of binase, while the latter still decreases the amount of cellular RNA. Using green fluorescent protein as an expression marker for the c-Kit oncoprotein, we demonstrate that the elimination of c-Kit is the key factor in selective cytotoxicity of binase. Quantitative RT-PCR with RNA samples isolated from the binase-treated FDC-P1-N822K cells shows that binase treatment results in 41% reduction in the amount of с-kit mRNA. This indicates that the transcript of the activated mutant c-kit is the target for toxic action of binase. Thus, the combination of inhibition of oncogenic protein with the destruction of its mRNA is a promising approach to eliminating malignant cells.

Genome-wide profiling of forum domains in Drosophila melanogaster

Forum domains are stretches of chromosomal DNA that are excised from eukaryotic chromosomes during their spontaneous non-random fragmentation. Most forum domains are 50-200 kb in length. We mapped forum domain termini using FISH on polytene chromosomes and we performed genome-wide mapping using a Drosophila melanogaster genomic tiling microarray consisting of overlapping 3 kb fragments. We found that forum termini very often correspond to regions of intercalary heterochromatin and regions of late replication in polytene chromosomes. We found that forum domains contain clusters of several or many genes. The largest forum domains correspond to the main clusters of homeotic genes inside BX-C and ANTP-C, cluster of histone genes and clusters of piRNAs. PRE/TRE and transcription factor binding sites often reside inside domains and do not overlap with forum domain termini. We also found that about 20% of forum domain termini correspond to small chromosomal regions where Ago1, Ago2, small RNAs and repressive chromatin structures are detected. Our results indicate that forum domains correspond to big multi-gene chromosomal units, some of which could be coordinately expressed. The data on the global mapping of forum domains revealed a strong correlation between fragmentation sites in chromosomes, particular sets of mobile elements and regions of intercalary heterochromatin.

Evidence for RNA synthesis in the intergenic region between enhancer and promoter and its inhibition by insulators in Drosophila melanogaster

Uncovering the nature of communication between enhancers, promoters and insulators is important for understanding the fundamental mechanisms that ensure appropriate gene expression levels. Here we describe an approach employing transient expression of genetic luciferase reporter gene constructs with quantitative RT–PCR analysis of transcription between an enhancer and Hsp70 promoter. We tested genetic constructs containing gypsy and/or Fab7 insulators in different orientations, and an enhancer from copia LTR-retroelement [(enh)copia]. A single gypsy or Fab7 insulator inserted between the promoter and enhancer in any polarity reduced enhancer action. A pair of insulators flanking the gene in any orientation exhibited increased insulation activity. We detected promoter-independent synthesis of non-coding RNA in the intergenic region of the constructs, which was induced by the enhancer in both directions and repressed by a single insulator or a pair of insulators. These results highlight the involvement of RNA-tracking mechanisms in the communications between enhancers and promoters, which are inhibited by insulators.

Suffix-specific RNAi leads to silencing of F element in Drosophila melanogaster

Separate conserved copies of suffix, a short interspersed Drosophila retroelement (SINE), and also divergent copies in the 3′ untranslated regions of the three genes, have already been described. Suffix has also been identified on the 3′ end of the Drosophila non-LTR F element, where it forms the last conserved domain of the reverse transcriptase (RT). In our current study, we show that the separate copies of suffix are far more actively transcribed than their counterparts on the F element. Transcripts from both strands of suffix are present in RNA preparations during all stages of Drosophila development, providing the potential for the formation of double-stranded RNA and the initiation of RNA interference (RNAi). Using in situ RNA hybridization analysis, we have detected the expression of both sense and antisense suffix transcripts in germinal cells. These sense and antisense transcripts are colocalized in the primary spermatocytes and in the cytoplasm of the nurse cells, suggesting that they form double-stranded RNA. We performed further analyses of suffix-specific small RNAs using northern blotting and SI nuclease protection assays. Among the total RNA preparations isolated from embryos, larvae, pupae and flies, suffix-specific small interfering RNAs (siRNAs) were detected only in pupae. In wild type ovaries, both the siRNAs and longer suffix-specific Piwi-interacting RNAs (piRNAs) were observed, whereas in ovaries of the Dicer-2 mutant, only piRNAs were detected. We further found by 3′ RACE that in pupae and ovaries, F element transcripts lacking the suffix sequence are also present. Our data provide direct evidence that suffix-specific RNAi leads to the silencing of the relative LINE (long interspersed element), F element, and suggests that SINE-specific RNA interference could potentially downregulate a set of genes possessing SINE stretches in their 5′ or 3′ non-coding regions. These data also suggest that double stranded RNAs possessing suffix are processed by both RNAi and an additional silencing mechanism.

Molecular Mechanisms of Epigenetics

The main epigenetic mechanisms in regulation of gene expression are discussed. The definition of epigenetics and its specific mechanisms including DNA methylation and gene imprinting, modifications of nucleosomal histones associated with silencing or activation of gene transcription, RNA interference, chromosomal silencing, and the role of mobile elements are discussed.

SuUR protein binds to the boundary regions separating forum domains in Drosophila melanogaster

Forum domains are 50-150 kb DNA fragments that are released during spontaneous fragmentation of chromosomes. They are separated by islands of putative heterochromatin boundary regions. The SuUR protein, which is involved in the control of chromosome organization, is localized exclusively in heterochromatin and often colocalizes on chromosomes with Polycomb group proteins. To test whether the SuUR protein is associated with boundary regions, we used gel retardation assays and found that the SuUR protein binds specifically to boundary regions and that boundary regions are under-replicated. These results suggest that the regular distribution of boundary regions in chromosomes may represent the dispersion of sites designed for chromosomal silencing.

Gene-specific silencing by expression of parallel complementary RNA in Escherichia coli

Gene-specific silencing refers to a phenomenon in which expression of an individual gene can be specifically repressed by different mechanisms on the levels of transcription, RNA splicing, transport, degradation in nuclei or cytoplasm, or blocking of translation. In different species gene-specific silencing was observed by expression or injections of antiparallel double-stranded RNA formed by a fragment of mRNA and antisense RNA. Here we show a potent and specific gene silencing in bacteria by expression of RNA, that is complementary in a parallel orientation to Escherichia coli lon mRNA. Moreover, the expression of parallel RNA is more effective at producing interference than expression of antisense RNA corresponding to the same mRNA region. Both effects of interference mediated either by parallel RNA or antiparallel RNA gradually decrease up to the 40th generation. Together with in vitro nuclease protection studies these results indicate that a parallel RNA duplex might be formed in vivo and both types of duplexes, antiparallel or parallel, can induce gene-specific silencing by similar mechanisms.

Forum domain in Drosophila melanogaster cut locus possesses looped domains inside

We have studied the relationship between chromosomal forum domains and looped domains in the cut locus of Drosophila melanogaster . Forum domains were earlier detected by separation in pulsed-field gels of 50-150 kb chromosomal DNA fragments obtained after spontaneous non-random degradation of chromosomes. We have localized the boundary region where cleavage sites are scattered between two forum domains in the regulatory region of the cut locus. We have sequenced a 13 kb region spanning few kilobases from distal domain, the boundary region and part of the proximal forum domain where several scaffold associated regions (SARs) were observed. We conclude that forum domains and looped domains are physically different types of domains and belong to different levels of organization in eukaryotic chromosomes. The boundary region between the neighboring forum domains in the cut locus possesses the Doc element insertion and a micro-satellite stretch and thus might remind a small island of heterochromatin and correspond to so-called intercalary heterochromatin that is known to be located in the 7B1-2 band where the major part of the cut locus is reside.

IncI1 plasmid R64 encodes the ArsR protein that alleviates type I restriction

The host-controlled EcoK restriction of unmodified phage lambda was five-fold alleviated in the wild-type Escherichia coli strain K12 carrying the R64 plasmid of the incompatibility group I1. The relevant gene was mapped between the origin of vegetative replication (rep, oriV) and the tet(r) gene about 60 kbp downstream from the origin of transfer, oriT. We cloned this gene inside the 613 bp long EcoRI-PstI fragment and sequenced it. Only one 351 bp long open reading frame (ORF) starting at 124 bp from the beginning of the insert was found in the sequence. Computer search in the current databases revealed that the putative protein is identical to the ArsR protein specified by the IncFI plasmid R773. ArsR is a repressor of the arsenical resistance (ars) operon, arsRDABC. There are no arsABC genes in the R64 plasmid since plasmid R64- (or pSR8)-mediated resistance of E. coli K12 cells to the arsenicals arsenate and arsenite was not detected. The gene arsR and the antirestriction genes ard (ardA and ardB) are non-homologous. However, comparison of the deduced amino acid sequence of ArsR with the ArdA and ArdB sequences revealed only one small region of similarity, a 9 amino acid motif found in different antirestriction proteins that is hypothesized to be an interaction site for antirestriction proteins with restriction endonucleases.

Burdock, a novel retrotransposon in Drosophila melanogaster, integrates into the coding region of the cut locus

The burdock element is known to be the 2.6-kb insertion into the same region of the cut locus in 12 independently obtained ct-lethal mutants. Here we have determined the complete sequences of this insertion and of the hot spot region. It was found that the burdock is a short retrotransposon with long terminal repeats and a single open reading frame (ORF). The polypeptide encoded by the burdock ORF contains two adjacent regions homologous to the gag and pol polyproteins of the gypsy mobile element. The burdock insertion interrupts the short ORF of the cut locus. The target site sequence of the burdock insertions is similar to the Drosophila topoisomerase II cleavage site.