Quenching of G4-DNA intrinsic fluorescence by ligands

G-quadruplex (G4) structures formed by the guanine-rich DNA regions exhibit several distinctive optical properties, including UV absorption and circular dichroism spectra. Some G4 DNA possess intrinsic UV fluorescence whose origin is not completely clear to date. In this work, we study the effect of TMPyP4 and Methylene Blue on the intrinsic fluorescence of the dimeric G4 DNA structure formed by two d(G3T)4 sequences. We demonstrate that binding of the ligands results in quenching of the intrinsic fluorescence, although the conformation of the G4 DNA and its dimeric structure remain preserved. The binding sites of the ligands were suggested by the photoinduced oxidation of guanines and analysis of binding isoterms. We discuss how DNA-ligand complexes can affect the intrinsic fluorescence of G4 DNA.

Taq-Polymerase Stop Assay to Determine Target Selectivity of G4 Ligands in Native Promoter Sequences of MYC, TERT, and KIT Oncogenes

Computational and high-throughput experimental methods predict thousands of potential quadruplex sequences (PQSs) in the human genome. Often these PQSs contain more than four G-runs, which introduce additional uncertainty into the conformational polymorphism of the G4 DNA. G4-specific ligands, which are currently being actively developed as potential anticancer agents or tools for studying G4 structures in genomes, may preferentially bind to specific G4 structures over the others that can be potentially formed in the extended G-rich genomic region. We propose a simple technique that identifies the sequences that tend to form G4 in the presence of potassium ions or a specific ligand. Thermostable DNA Taq-polymerase stop assay can detect the preferential position of the G4 -ligand binging within a long PQS-rich genomic DNA fragment. This technique was tested for four G4 binders PDS, PhenDC3, Braco-19, and TMPyP4 at three promoter sequences of MYC, KIT, and TERT that contain several PQSs each. We demonstrate that the intensity of polymerase pausing reveals the preferential binding of a ligand to particular G4 structures within the promoter. However, the strength of the polymerase stop at a specific site does not always correlate with the ligand-induced thermodynamic stabilization of the corresponding G4 structure.

Genome Assembly and Sex-Determining Region of Male and Female Populus × sibirica

The genus Populus is presented by dioecious species, and it became a promising object to study the genetics of sex in plants. In this work, genomes of male and female Populus × sibirica individuals were sequenced for the first time. To achieve high-quality genome assemblies, we used Oxford Nanopore Technologies and Illumina platforms. A protocol for the isolation of long and pure DNA from young poplar leaves was developed, which enabled us to obtain 31 Gb (N50 = 21 kb) for the male poplar and 23 Gb (N50 = 24 kb) for the female one using the MinION sequencer. Genome assembly was performed with different tools, and Canu provided the most complete and accurate assemblies with a length of 818 Mb (N50 = 1.5 Mb) for the male poplar and 816 Mb (N50 = 0.5 Mb) for the female one. After polishing with Racon and Medaka (Nanopore reads) and then with POLCA (Illumina reads), assembly completeness was 98.45% (87.48% duplicated) for the male and 98.20% (76.77% duplicated) for the female according to BUSCO (benchmarking universal single-copy orthologs). A high proportion of duplicated BUSCO and the increased genome size (about 300 Mb above the expected) pointed at the separation of haplotypes in a large part of male and female genomes of P. × sibirica. Due to this, we were able to identify two haplotypes of the sex-determining region (SDR) in both assemblies; and one of these four SDR haplotypes, in the male genome, contained partial repeats of the ARR17 gene (Y haplotype), while the rest three did not (X haplotypes). The analysis of the male P. × sibirica SDR suggested that the Y haplotype originated from P. nigra, while the X haplotype is close to P. trichocarpa and P. balsamifera species. Moreover, we revealed a Populus-specific repeat that could be involved in translocation of the ARR17 gene or its part to the SDR of P. × sibirica and other Populus species. The obtained results expand our knowledge on SDR features in the genus Populus and poplar phylogeny.

Genome and Transcriptome Sequencing of Populus × sibirica Identified Sex-Associated Allele-Specific Expression of the CLC Gene

Transcriptome sequencing of leaves, catkin axes, and flowers from male and female trees of Populus × sibirica and genome sequencing of the same plants were performed for the first time. The availability of both genome and transcriptome sequencing data enabled the identification of allele-specific expression. Such an analysis was performed for genes from the sex-determining region (SDR). P. × sibirica is an intersectional hybrid between species from sections Aigeiros (Populus nigra) and Tacamahaca (Populus laurifolia, Populus suaveolens, or Populus × moskoviensis); therefore, a significant number of heterozygous polymorphisms were identified in the SDR that allowed us to distinguish between alleles. In the SDR, both allelic variants of the TCP (T-complex protein 1 subunit gamma), CLC (Chloride channel protein CLC-c), and MET1 (DNA-methyltransferase 1) genes were expressed in females, while in males, two allelic variants were expressed for TCP and MET1 but only one allelic variant prevailed for CLC. Targeted sequencing of TCP, CLC, and MET1 regions on a representative set of trees confirmed the sex-associated allele-specific expression of the CLC gene in generative and vegetative tissues of P. × sibirica. Our study brings new knowledge on sex-associated differences in Populus species.

Recurrent Potential G-Quadruplex Sequences in Archaeal Genomes

Evolutionary conservation or over-representation of the potential G-quadruplex sequences (PQS) in genomes are usually considered as a sign of the functional relevance of these sequences. However, uneven base distribution (GC-content) along the genome may along the genome may result in seeming abundance of PQSs over average in the genome. Apart from this, a number of other conserved functional signals that are encoded in the GC-rich genomic regions may inadvertently result in emergence of G-quadruplex compatible sequences. Here, we analyze the genomes of archaea focusing our search to repetitive PQS (rPQS) motifs within each organism. The probability of occurrence of several identical PQSs within a relatively short archaeal genome is low and, thus, the structure and genomic location of such rPQSs may become a direct indication of their functionality. We have found that the majority of the genomes of Methanomicrobiaceae family of archaea contained multiple copies of the interspersed highly similar PQSs. Short oligonucleotides corresponding to the rPQS formed the G-quadruplex (G4) structure in presence of potassium ions as demonstrated by circular dichroism (CD) and enzymatic probing. However, further analysis of the genomic context for the rPQS revealed a 10–12 nt cytosine-rich track adjacent to 3'-end of each rPQS. Synthetic DNA fragments that included the C-rich track tended to fold into alternative structures such as hairpin structure and antiparallel triplex that were in equilibrium with G4 structure depending on the presence of potassium ions in solution. Structural properties of the found repetitive sequences, their location in the genomes of archaea, and possible functions are discussed.

[Functional Hypermethylation of ALDH1L1, PLCL2, and PPP2R3A in Colon Cancer]

DNA hypermethylation and mutations are key mechanisms for the downregulation of tumor suppressor genes. NotI-microarrays allowed us to detect hypermethylation and/or deletions in 180 NotI sites associated with 188 genes of human chromosome 3, in 24 paired (tumor/normal) colon samples. The most frequent aberrations (in more than 20% of tumor samples) were detected in the promoter regions of 20 genes. Expression and promoter methylation of these genes were analyzed using the data for paired colon samples from The Cancer Genome Atlas project. Three genes - ALDH1L1, PLCL2, and PPP2R3A - revealed a more than two-fold average decrease in expression and a negative correlation between mRNA level and promoter hypermethylation. The expression of these three genes was then evaluated in 30 paired colon samples by quantitative PCR. Frequent (in more than 60% of cases) and significant (5-9-fold on average) mRNA level decrease was found for each of the genes in the tumor samples. The results indicate a suppressor role of the ALDH1L1, PLCL2, and PPP2R3A genes in colon cancer, as well as functional significance of hypermethylation in the downregulation of these genes.

Stable G-Quadruplex Structures of Oncogene Promoters Induce Potassium-Dependent Stops of Thermostable DNA Polymerase

Amplification of GC-rich regions of genomic DNA is hindered either by high stability of DNA double helix or as a result of alternative structure formation by a guanine-rich DNA strand. Such potential G-quadruplex (G4) sequences are fairly common in promoters of the human genome. The efficiency of PCR amplification of promoter sequences for several human oncogenes (MYC, NRAS, TERT, KRAS, KIT) was studied. We demonstrate that the efficiency of DNA polymerase is reduced in the presence of potassium ions. The primer-extension technique localized DNA polymerase stops at the 3'-ends of potential quadruplex sequences. The structural and thermodynamic properties of short G-rich oligonucleotides corresponding to the stops of DNA polymerase were analyzed. These oligonucleotides formed stable parallel G4 in the presence of potassium ions. Correlation between the stability of G4 structure and efficiency of DNA polymerase stops was revealed. The results provide a method for detecting new G4 structures in extended genomic sequences and also clarify the mechanism of inhibition of DNA polymerase in G-rich regions of DNA.

Pan-Cancer Analysis of TCGA Data Revealed Promising Reference Genes for qPCR Normalization

Quantitative PCR (qPCR) remains the most widely used technique for gene expression evaluation. Obtaining reliable data using this method requires reference genes (RGs) with stable mRNA level under experimental conditions. This issue is especially crucial in cancer studies because each tumor has a unique molecular portrait. The Cancer Genome Atlas (TCGA) project provides RNA-Seq data for thousands of samples corresponding to dozens of cancers and presents the basis for assessment of the suitability of genes as reference ones for qPCR data normalization. Using TCGA RNA-Seq data and previously developed CrossHub tool, we evaluated mRNA level of 32 traditionally used RGs in 12 cancer types, including those of lung, breast, prostate, kidney, and colon. We developed an 11-component scoring system for the assessment of gene expression stability. Among the 32 genes, PUM1 was one of the most stably expressed in the majority of examined cancers, whereas GAPDH, which is widely used as a RG, showed significant mRNA level alterations in more than a half of cases. For each of 12 cancer types, we suggested a pair of genes that are the most suitable for use as reference ones. These genes are characterized by high expression stability and absence of correlation between their mRNA levels. Next, the scoring system was expanded with several features of a gene: mutation rate, number of transcript isoforms and pseudogenes, participation in cancer-related processes on the basis of Gene Ontology, and mentions in PubMed-indexed articles. All the genes covered by RNA-Seq data in TCGA were analyzed using the expanded scoring system that allowed us to reveal novel promising RGs for each examined cancer type and identify several "universal" pan-cancer RG candidates, including SF3A1, CIAO1, and SFRS4. The choice of RGs is the basis for precise gene expression evaluation by qPCR. Here, we suggested optimal pairs of traditionally used RGs for 12 cancer types and identified novel promising RGs that demonstrate high expression stability and other features of reliable and convenient RGs (high expression level, low mutation rate, non-involvement in cancer-related processes, single transcript isoform, and absence of pseudogenes).

[Differential expression of an ensemble of the key genes involved in cell-cycle regulation in lung cancer]

Targeted cancer therapy directed at individual targets is often accompanied by the rapid development of drug resistance. The development of a new generation of antitumor drugs involves the search for many targets simultaneously to block or, conversely, restore their activity. In this regard, simultaneous analysis of gene expression in a complex network of interactions, primarily cell cycle control elements, is relevant for the search of specific molecular markers for the differential diagnosis of adenocarcinoma (ADC) and squamous cell lung cancer (SCC), as well as new targets for therapy. In this paper we performed an extended quantitative analysis of the expression of two suppressor genes, CTDSPL and its target RB1, as well as 84 genes of the main participants of the p16^(INK4A)-Cdk/cyclin D1-Rb and p53/p21^(Waf1) signaling pathways in the histological types of non-small-cell lung cancer (NSCLC), i.e., ADC and SCC, using the special panel of the Human Cell Cycle Regulation Panel. The expression profile of some genes shows the specificity to the histological type of NSCLC and the presence of metastases. The genes with a significantly increased expression that affect the activity of Rb (cyclins, cyclin-dependent kinases, their activators, inhibitors, etc.) can serve as potential targets for combined therapy of both ADC and SCC.

Deep Sequencing Revealed a CpG Methylation Pattern Associated With ALDH1L1 Suppression in Breast Cancer

Hypermethylation of promoter CpG islands is generally recognized epigenetic mechanism responsible for gene silencing in cancer. However, molecular details on how this epigenetic mark triggers the process of gene downregulation are still elusive. Here, we used deep bisulfite sequencing and qPCR analysis to investigate the pattern of CpG methylation of ALDH1L1 promoter region and its association with the gene expression level in 16 paired breast cancer (BC) samples of different clinical stages. Expression of ALDH1L1 gene was suppressed in all examined BC samples up to 200-fold, and average hypermethylation level of the promoter region correlated positively with ALDH1L1 downregulation. We determined the role of every individual CpG site within the ALDH1L1 promoter, including upstream untranscribed region, first untranslated exon, and the start of the first intron, in aberrant gene expression by correlation analysis. The search revealed CpG sites which methylation has the highest impact on intensity of gene transcription. The majority of such CpG sites are located in a compact region in the first intron of the ALDH1L1 gene. These results assist in unraveling of dynamic nature of CpG promoter hypermethylation as well as demonstrate the efficiency of deep bisulfite sequencing in search for novel epigenetic markers in cancer.

[Interaction of two tumor suppressors: Phosphatase CTDSPL and Rb protein]

Earlier we established that CTDSPL gene encoding small carboxy-terminal domain serine phosphatase can be considered a classical tumor suppressor gene. Besides, transfection of tumor cell line MCF-7 with CTDSPL led to the content decrease of inactive phosphorylated form of another tumor suppressor, retinoblastoma protein (Rb), and subsequently to cell cycle arrest at the G1/S boundary. This result implied that small phosphatase CTDSPL is able to specifically dephosphorylate and activate Rb protein. In order to add some fuel to this hypothesis, in the present work we studied the interaction of two tumor suppressors CTDSPL and Rb in vitro. GST pool-down assay revealed that CTDSPL is able to precipitate Rb protein from MCF-7 cell extracts, while surface plasmon resonance technique showed that interaction of the two proteins is direct. Results of this study reassert that phosphatase CTDSPL and Rb could be involved in the common mechanism of cell cycle regulation.

Light-induced oxidation of the telomeric G4 DNA in complex with Zn(II) tetracarboxymethyl porphyrin

Structure-specific ligands are convenient tools for the recognition, targeting or probing of non-canonical DNA structures. Porphyrin derivatives exhibit a preference for interaction with G-quadruplex (G4) structures over canonical duplex DNA and are able to cause photoinducible damage to nucleic acids. Here, we show that Zn(II) 5,10,15,20-tetrakis(N-carboxymethyl-4-pyridinium)porphyrin (ZnP1) interacts with different conformations of the telomeric sequence d(TAGGG(TTAGGG)₃) at submicromolar concentrations without any detectible disturbance of the particular fold. Among different folds, potassium (3+1) hybrid G4-structure. reveal the highest affinity to ZnP1. The pattern of guanine oxidation is specific for each telomeric DNA conformation and may serve as an additional tool for probing the G4 topology. The potassium (3+1) and parallel G4 conformations are more susceptible to light-induced oxidation than the sodium G4 conformation or double helix of the telomeric DNA. The major products of the guanine modifications are spiroiminodihydantoin (Sp) and 8-oxoguanine (8-oxoG). ZnP1-induced oxidation of guanines results in the structural rearrangement of parallel and (3+1) G4 conformations yielding an antiparallel-like G4 conformation. The mechanism of the observed light-induced conformational changes is discussed.

Identification of Novel Epigenetic Markers of Prostate Cancer by NotI-Microarray Analysis

A significant need for reliable and accurate cancer diagnostics and prognosis compels the search for novel biomarkers that would be able to discriminate between indolent and aggressive tumors at the early stages of disease. The aim of this work was identification of potential diagnostic biomarkers for characterization of different types of prostate tumors. NotI-microarrays with 180 clones associated with chromosome 3 genes/loci were applied to determine genetic and epigenetic alterations in 33 prostate tumors. For 88 clones, aberrations were detected in more than 10% of tumors. The major types of alterations were DNA methylation and/or deletions. Frequent methylation of the discovered loci was confirmed by bisulfite sequencing on selective sampling of genes: FGF12, GATA2, and LMCD1. Three genes (BHLHE40, BCL6, and ITGA9) were tested for expression level alterations using qPCR, and downregulation associated with hypermethylation was shown in the majority of tumors. Based on these data, we proposed the set of potential biomarkers for detection of prostate cancer and discrimination between prostate tumors with different malignancy and aggressiveness: BHLHE40, FOXP1, LOC285205, ITGA9, CTDSPL, FGF12, LOC440944/SETD5, VHL, CLCN2, OSBPL10/ZNF860, LMCD1, FAM19A4, CAND2, MAP4, KY, and LRRC58. Moreover, we probabilistically estimated putative functional relations between the genes within each set using the network enrichment analysis.

[G-quadruplex ligands: mechanisms of anticancer action and target binding]

Since the discovery of potential therapeutic value of quadruplex secondary nucleic acids structures, many compounds that stabilize these targets were found. Such progress became possible due to understanding of the structural features of G-quadruplexes. Quadruplex ligands selectively suppress the growth of tumor cells by indirect inhibition of the telomerase activity and/or attenuation of oncogenes' overexpression. Therapeutic effect demonstrated in vivo supports the prospect of such compounds for the development of the targeted anticancer drugs. This review reveals the significance of G-quadruplexes as therapeutic targets and focuses on biochemical properties of the low molecular weight quadruplex ligands.

2-aminopurine fluorescence: discrimination between specific and unspecific ligand binding to the kissing-loop dimer of the HIV-1 RNA

The fluorescent 2-aminopurine probe (2-AP) incorporated into the loop of 23-mer RNA hairpin of HIV-1 genome dimerization initiation site (DIS) was used for discrimination of specific and unspecific binding of paromomycin and spermine to the kissing loop dimer (KD) formed in solution. While both ligands stabilized the KD RNA structure, only paromomycin binding resulted in significant increase of 2-AP fluorescence. These observations suggest that the 2-AP fluorescent RNA construct might be useful for selecting ligands specifically binding the HIV-1 kissing loop RNA dimer.

A protein component at the heart of an RNA machine: the importance of protein l27 for the function of the bacterial ribosome

Deletion of the gene for protein L27 from the E. coli chromosome results in severe defects in cell growth. This deficiency is corrected by the expression of wild-type (wt) protein L27 from a plasmid. Examination of strains expressing L27 variants truncated at the N terminus reveals that the absence of as few as three amino acids leads to a decrease in growth rate, an impairment in peptidyl transferase activity, and a sharp decline in the labeling of L27 from the 3' end of a photoreactive tRNA at the ribosomal P site. These findings suggest that the flexible N-terminal sequence of L27, which protrudes onto the interface of the bacterial 50S subunit, can reach the peptidyl transferase active site and contribute to its function, possibly by helping to correctly position tRNA substrates at the catalytic site.

[Characteristic of the two-step RNA dimerization in avian sarcoma and leukosis viruses]

Dimerization of two copies of genomic RNA is a necessary step of retroviral replication. In the case of human immunodeficiency virus type 1 (HIV-1) the process is explored in many details. It is proved that conserved stem-loop structure is an essential element in RNA dimerization. Similar model of two-step dimerization mechanism can be considered for avian sarcoma and leukosis virus group (ASLV) in spite of the absence of homology between dimer initiation site (DIS) of ASLV and that of HIV-1. In this paper, short RNA fragments of two viruses: avian sarcoma virus CT-10 and avian leukosis virus HPRS-103 have been chosen in order to investigate the structural requirements of dimerization process and compare them to that of HIV-1. The rate of spontaneous transition from loose to tight dimer was studied as a function of stem length and temperature. Although both types of dimers were observed for both avian retroviruses chosen, fragments of CT-10 requires much higher RNA concentration to form loose dimer. In spite of identical sequence of the loops (5'-A-CUGCAG-3') avian sarcoma virus CT-10 RNA fragments dimerization was greatly impaired. The differences can be explained by deletion of adenine 271 in avian sarcoma virus CT-10 in the stem and by resulting shortening of the self-complementary loop.

[The role of structural elements in dimerization of RNA fragments in avian retroviruses]

The slipped loop structure, earlier identified as an unusual DNA structure, was found to be a possible element of the RNA folding. In order to experimentally test this suggestion, model oligoribonucleotides capable of forming the SLS were synthesized. Treatment of the oligoribonucleotides with nuclease S1 and RNases specific for single- and double-stranded RNA demonstrated the steric possibility of SLS formation. To determine the possible functional role of SLS-RNA, various naturally occurring RNAs were screened in silico. Among the most interesting findings were dimerization initiation sites of avian retroviral genomic RNAs. Analysis of RNA from 31 viruses showed that formation of the intermolecular SLS during RNA dimerization is theoretically possible, competing with the formation of an alternative hairpin structure. Identification of the secondary structure of selected RNA dimers employing nuclease digestion techniques as well as covariance analysis of the retroviral RNA dimerization initiation site sequences were used to show that the alternative conformation (loop-loop interaction of two hairpins) is the most preferred. Alternative structures and conformational transitions in RNA dimerization mechanisms in avian retroviruses are discussed.

[Why termination tRNAs have not been found? Because they are hidden in the large ribosomal RNA]

It is well known that protein synthesis in ribosomes on mRNA requires two kinds of tRNAs: initiation and elongation. The former initiates the process (formylmethionine tRNA in prokaryotes and special methionine tRNA in eukaryotes). The latter participates in the synthesis proper, recognizing the sense codons. The synthesis is assisted by special proteins: initiation, elongation, and termination factors. The termination factors are necessary to recognize stop codons (UAG, UGA, and UAA) and to release the complete protein chain from the elongation tRNA preceding a stop codon. No termination tRNA capable of recognizing stop codons by its anticodon is known. The termination factors are thought to do this. We discovered in the large ribosomal RNA two regions that, like tRNAs, contain the anticodon hairpin, but with triplets complementary to stop codons. By analogy, we called them termination tRNAs (Ter-tRNA1 and Ter-tRNA2), though they transport no amino acids, and suggested them to directly recognize stop codons. The termination factors only condition such a recognition, making it specific and reliable (of course, they fulfill the hydrolysis of the ester bond between the polypeptide and tRNA). A strong argument in favor of our hypothesis came from vertebrate mitochondria. They acquired two new stop codons, AGA and AGG (in the standard code, they are two out of six arginine codons). We revealed that the corresponding anticodons appear in Ter-tRNA1.

A pseudoknot-compatible universal site is located in the large ribosomal RNA in the peptidyltransferase center

The RNA secondary structure is not confined to a system of the hairpins and can contain pseudoknots as well as topologically equivalent slipped-loop structure (SLS) conformations. A specific primary structure that directs folding to the pseudoknot or SLS is called SL-palindrome (SLP). Using a computer program for searching the SLP in the genomic sequences, 419 primary structures of large ribosomal RNAs from different kingdoms (prokaryota, eukaryota, archaebacteria) as well as plastids and mitochondria were analyzed. A universal site was found in the peptidyltransferase center (PTC) capable of folding to a pseudoknot of 48 nucleotides in length. Phylogenetic conservation of its helices (concurrent replacements with no violation of base pairing, covariation) has been demonstrated. We suggest the reversible folding-unfolding of the pseudoknot for certain stages of the ribosome functioning.