Prediction of EVT6-NTRK3-Dependent Papillary Thyroid Cancer Using Minor Expression Profile.

Solid tumors resulting from oncogenic stimulation of neurotrophin receptors (TRK) by chimeric proteins are a group of rare tumors of various localization that respond to therapy with targeted drugs entrectinib and larotrectinib. The standard method for detecting chimeric TRK genes in tumor samples today is considered to be next generation sequencing with the determination of the prime structure of the chimeric transcripts. We hypothesized that expression of the chimeric tyrosine kinase proteins in tumors can determine the specific transcriptomic profile of tumor cells. We detected differentially expressed genes allowing distinguishing between TRK-dependent tumors papillary thyroid cancer (TC) from other molecular variants of tumors of this type. Using PCR with reverse transcription (RT-PCR), we identified 7 samples of papillary TC carrying a EVT6-NTRK3 rearrangement (7/215, 3.26%). Using machine learning and the data extracted from TCGA, we developed of a recognition function for predicting the presence of rearrangement in NTRK genes based on the expression of 10 key genes: AUTS2, DTNA, ERBB4, HDAC1, IGF1, KDR, NTRK1, PASK, PPP2R5B, and PRSS1. The recognition function was used to analyze the expression data of the above genes in 7 TRK-dependent and 10 TRK-independent thyroid tumors obtained by RT-PCR. On the test samples from TCGA, the sensitivity was 72.7%, the specificity - 99.6%. On our independent validation samples tested by RT-PCR, sensitivity was 100%, specificity - 70%. We proposed an mRNA profile of ten genes that can classify TC in relation to the presence of driver NTRK-chimeric TRK genes with acceptable sensitivity and specificity.

[THE BIOMARKERS FOR TIMELY DIAGNOSTICS OF COLORECTAL CANCER].

The colorectal cancer (CC) is one of the most widespread type of cancer all over the world. It is confirmed that the screening procedures intended for timely detection of CC and adenomatous polyps, significantly decrease mortality. The colonoscopy and analysis offeces for occult blood are widely applied as screening procedures. However, they have a number of shortcomings. The studies of the last decade revealed number of genetic and epigenetic markers potentially permitting revealing patients with CC at early stages of development of disease. The article analyzes CC-specific microRNA and their possible interactions with different transcriptional factors. These factors, being integrated into the structure of so called network s with direct signal propagation, ensure special stability of all regulatory system. The derangement of functioning of these networks quite often results in pathological alterations.

[Comparative analysis of the DXPas34 regulatory region in rodents].

Mouse X chromosome inactivation center contains the DXPas34 minisatellite locus which plays an important role in expression regulation of the Tsix and Xist genes, involved into female dosage compensation. Comparative analysis of the DXPas34 locus from mouse, rat, and four common vole species revealed similar organization of this region in the form of tandem repeat blocks. A search for functionally important elements in this locus showed that all the species examined carried the conservative motif monomers, which could be involved in regulation of X inactivation.

[Molecular genetic characterization of the regulatory region of the Xist gene in the common vole Microtus rossiaemeridionalis].

Two conserved regions were discovered as a result of interspecific comparison of the 5'-region of the Xist gene, which is the key gene in the process of X-chromosome inactivation in mammalian females. The first region corresponds to the minimal promoter, and the second spans between -480 bp and -400 bp from the start of Xist transcription. Footprinting experiments revealed protected regions corresponding to the potential binding sites for TBP, SP1, API, SRY, ER, and some other transcription factors. They also demonstrated the interaction with the minimal promoter of the human recombinant transcription factor SP1 in vitro and of the transcription factor CTCF in vivo. Experiments with reporter constructs showed that repressors of Xist transcription were located between -100 bp and -200 bp and between -300 bp and -400 bp and activators of Xist transcription were located between -200 bp and -300 bp and between -400 bp and -500 bp.

TRANSFAC and its module TRANSCompel: transcriptional gene regulation in eukaryotes.

The TRANSFAC database on transcription factors, their binding sites, nucleotide distribution matrices and regulated genes as well as the complementing database TRANSCompel on composite elements have been further enhanced on various levels. A new web interface with different search options and integrated versions of Match and Patch provides increased functionality for TRANSFAC. The list of databases which are linked to the common GENE table of TRANSFAC and TRANSCompel has been extended by: Ensembl, UniGene, EntrezGene, HumanPSD and TRANSPRO. Standard gene names from HGNC, MGI and RGD, are included for human, mouse and rat genes, respectively. With the help of InterProScan, Pfam, SMART and PROSITE domains are assigned automatically to the protein sequences of the transcription factors. TRANSCompel contains now, in addition to the COMPEL table, a separate table for detailed information on the experimental EVIDENCE on which the composite elements are based. Finally, for TRANSFAC, in respect of data growth, in particular the gain of Drosophila transcription factor binding sites (by courtesy of the Drosophila DNase I footprint database) and of Arabidopsis factors (by courtesy of DATF, Database of Arabidopsis Transcription Factors) has to be stressed. The here described public releases, TRANSFAC 7.0 and TRANSCompel 7.0, are accessible under http://www.gene-regulation.com/pub/databases.html.

Corrigendum to "Data on master regulators and transcription factor binding sites found by upstream analysis of multi-omics data on methotrexate resistance of colon cancer" [Data Brief. 10 (2016) 499-504].

[This corrects the article DOI: 10.1016/j.dib.2016.11.096.].

P-Match: transcription factor binding site search by combining patterns and weight matrices.

P-Match is a new tool for identifying transcription factor (TF) binding sites in DNA sequences. It combines pattern matching and weight matrix approaches thus providing higher accuracy of recognition than each of the methods alone. P-Match is closely interconnected with the TRANSFAC database. In particular, P-Match uses the matrix library as well as sets of aligned known TF-binding sites collected in TRANSFAC and therefore provides the possibility to search for a large variety of different TF binding sites. Using results of extensive tests of recognition accuracy, we selected three sets of optimized cut-off values that minimize either false negatives or false positives, or the sum of both errors. Comparison with the weight matrix approaches such as Matchtrade mark tool shows that P-Match generally provides superior recognition accuracy in the area of low false negative errors (high sensitivity). As familiar to the user of Matchtrade mark, P-Match also allows to save user-specific profiles that include selected subsets of matrices with corresponding TF-binding sites or user-defined cut-off values. Furthermore, a number of tissue-specific profiles are provided that were compiled by the TRANSFAC team. A public version of the P-Match tool is available at http://www.gene-regulation.com/cgi-bin/pub/programs/pmatch/bin/p-match.cgi.

MATCH: A tool for searching transcription factor binding sites in DNA sequences.

Match is a weight matrix-based tool for searching putative transcription factor binding sites in DNA sequences. Match is closely interconnected and distributed together with the TRANSFAC database. In particular, Match uses the matrix library collected in TRANSFAC and therefore provides the possibility to search for a great variety of different transcription factor binding sites. Several sets of optimised matrix cut-off values are built in the system to provide a variety of search modes of different stringency. The user may construct and save his/her specific user profiles which are selected subsets of matrices including default or user-defined cut-off values. Furthermore a number of tissue-specific profiles are provided that were compiled by the TRANSFAC team. A public version of the Match tool is available at: http://www.gene-regulation.com/pub/programs.html#match. The same program with a different web interface can be found at http://compel.bionet.nsc.ru/Match/Match.html. An advanced version of the tool called Match Professional is available at http://www.biobase.de.

TRANSFAC: transcriptional regulation, from patterns to profiles.

The TRANSFAC database on eukaryotic transcriptional regulation, comprising data on transcription factors, their target genes and regulatory binding sites, has been extended and further developed, both in number of entries and in the scope and structure of the collected data. Structured fields for expression patterns have been introduced for transcription factors from human and mouse, using the CYTOMER database on anatomical structures and developmental stages. The functionality of Match, a tool for matrix-based search of transcription factor binding sites, has been enhanced. For instance, the program now comes along with a number of tissue-(or state-)specific profiles and new profiles can be created and modified with Match Profiler. The GENE table was extended and gained in importance, containing amongst others links to LocusLink, RefSeq and OMIM now. Further, (direct) links between factor and target gene on one hand and between gene and encoded factor on the other hand were introduced. The TRANSFAC public release is available at http://www.gene-regulation.com. For yeast an additional release including the latest data was made available separately as TRANSFAC Saccharomyces Module (TSM) at http://transfac.gbf.de. For CYTOMER free download versions are available at http://www.biobase.de:8080/index.html.

Computer-assisted identification of cell cycle-related genes: new targets for E2F transcription factors.

The processes that take place during development and differentiation are directed through coordinated regulation of expression of a large number of genes. One such gene regulatory network provides cell cycle control in eukaryotic organisms. In this work, we have studied the structural features of the 5' regulatory regions of cell cycle-related genes. We developed a new method for identifying composite substructures (modules) in regulatory regions of genes consisting of a binding site for a key transcription factor and additional contextual motifs: potential targets for other transcription factors that may synergistically regulate gene transcription. Applying this method to cell cycle-related promoters, we created a program for context-specific identification of binding sites for transcription factors of the E2F family which are key regulators of the cell cycle. We found that E2F composite modules are found at a high frequency and in close proximity to the start of transcription in cell cycle-related promoters in comparison with other promoters. Using this information, we then searched for E2F sites in genomic sequences with the goal of identifying new genes which play important roles in controlling cell proliferation, differentiation and apoptosis. Using a chromatin immunoprecipitation assay, we then experimentally verified the binding of E2F in vivo to the promoters predicted by the computer-assisted methods. Our identification of new E2F target genes provides new insight into gene regulatory networks and provides a framework for continued analysis of the role of contextual promoter features in transcriptional regulation. The tools described are available at http://compel.bionet.nsc.ru/FunSite/SiteScan.html.

Eukaryotic mRNAs encoding abundant and scarce proteins are statistically dissimilar in many structural features.

It is well known that non-coding mRNA sequences are dissimilar in many structural features. For individual mRNAs correlations were found for some of these features and their translational efficiency. However, no systematic statistical analysis was undertaken to relate protein abundance and structural characteristics of mRNA encoding the given protein. We have demonstrated that structural and contextual features of eukaryotic mRNAs encoding high- and low-abundant proteins differ in the 5' untranslated regions (UTR). Statistically, 5' UTRs of low-expression mRNAs are longer, their guanine plus cytosine content is higher, they have a less optimal context of the translation initiation codons of the main open reading frames and contain more frequently upstream AUG than 5' UTRs of high-expression mRNAs. Apart from the differences in 5' UTRs, high-expression mRNAs contain stronger termination signals. Structural features of low- and high-expression mRNAs are likely to contribute to the yield of their protein products.

[Theoretical analysis of mechanisms of occurrence of DNA deletions in prokaryotic genomes based on direct repeats]. / Teoreticheskii analiz mekhanizmov vozniknoveniia deletsiii DNK v genomakh prokariot na osnove priamykh povtorov.

In the present work a mechanism of deletions emergence on the basis of complementary DNA regions mispairing of direct repeats has been investigated theoretically. A quantitative dependence of the rates of deletions emergence on such parameters of the flanking repeats as the nucleotide composition of repeats, the number of homology damages and the distance between repeated regions has been constructed. It has been proved, that using this relationship one can reliably evaluate the total rates of deletions emergence in the lacI gene sequence of E. coli according to the repeats arrangement in this gene.

[Statistical evidence for the correlation of DNA deletions in prokaryotic genomes with direct repeats]. / Statisticheskoe dokazatel'stvo sviazi deletsii DNK v genomakh prokariot s priamymi povtorami.

In the present work a computer analysis of deletion localization in the sequence of the E. coli lacI gene has been carried out by the statistical weight method. Reliable statistical correlation of the deletions location sites with the arrangement of the most perfect direct repeats revealing the shortest distance between repeated fragments has been shown. At the same time statistical analysis did not reveal reliable connection of deletions localization regions with the expected sites of gyrase recognition, sites and other recombination sites. A conclusion has been drawn, that the mechanism of deletions emergence on the basis of repeats appears to be predominant.

[Theoretical analysis of the DNA duplication mechanisms in the prokaryotic genomes on the basis of repeats]. / Teoreticheskii analiz mekhanizmov duplikatsii DNK v genomakh prokariot na osnove povtorov.

In the present work a theoretical analysis of the molecular mechanisms on duplications emergence in the genomes of prokaryotes on the basis of direct repeats has been carried out. The correlations obtained have shown, that the duplication rate depends on such parameters as the distance between repeated regions, repeats nucleotide composition and the number of homology damages in them. It has been revealed that the rate of duplications decreases more readily than the deletion rate upon the growth of the distance between the repeats. Such prevalence of deletions over duplications must lead to the elimination of various types of direct repeats from the prokaryotic genomes in the course of their evolution.

[Cloning universal probes for detecting mycoplasma contamination of cell cultures]. / Kolnirovanie universal'nykh zondov dlia vyiavleniia mikoplazmennogo zarazheniia kletochnykh kul'tur.

To obtain the universal polynucleotide hybridization probes for testing mycoplasmal contaminations in cell cultures, we have cloned several DNA fragments from the srRNA gene of Acholeplasma laidlawii. Before cloning, in order to exclude cross-hybridization of these probes with eukaryotic rRNA, the thermodynamic parameters of duplex formation between DNA complementary to mycoplasmal rRNA and eukaryotic rRNA had been studied. Using a set of computer methods, the region which forms weak heteroduplexes with eukaryotic srRNA was revealed. This region occupies positions 250 to 550 position of the mycoplasmal srRNA. Three different DNA fragments which include the region were generated in PCR, cloned in pUC18, and their hybridization characteristics were evaluated. In appropriate hybridization conditions the probes hybridize with all mycoplasmal RNAs studied without cross-hybridization with eukaryotic ribosomal RNA and DNA, and allow one to detect virtually any mycoplasmas (or any prokaryote) in cell cultures. Blot-hybridization of universal probes with mycoplasmal DNA digested by BsuRI allows one to identify the different species of mycoplasmas.

[Modeling dynamics of gene net, regulating the cell cycle in mammalian cells]. / Modelirovanie dinamiki gennykh setei, reguliruiushchikh kletochnyi tsikl v kletakh mlekopitaiushchikh.

The study of the molecular mechanisms determining cellular programs of proliferation, differentiation, and apoptosis is currently attracting much attention. Recent studies have demonstrated that the system of cell-cycle control based on the transcriptional regulation of the expression of specific genes is responsible for the transition between programs. These groups of functionally connected genes from so-called gene networks characterized by numerous feedbacks and a complex behavioral dynamics. Computer simulation methods have been applied to studying the dynamics of gene networks regulating the cell cycle of vertebrates. The data on the regulation of the key genes obtained from the CYCLE-TRRD database have been used as a basis to construct gene networks of different degrees of complexity controlling the G1/S transition, one of the most important stages of the cell cycle. The behavior dynamics of the model constructed has been analyzed. Two qualitatively different functional modes of the system has been obtained. It has also been shown that the transition between these modes depends on the duration of the proliferation signal. It has also been demonstrated that the additional feedback from factor E2F to genes c-fos and c-jun, which was predicted earlier based on the computer analysis of promoters, plays an important role in the transition of the cell to the S phase.

[GeneExpress: an integrator for databases and computer systems accessible by the Internet and intended for studying eukaryotic gene expression]. / Gene Express: integrator baz dannykh i komp'iuternykh sistem, dostupnykh po seti internet i prednaznachennykh dlia izucheniia ékspresii genov éukariot.

We have developed GeneExpress that is the WWW-oriented integrator for the databases and systems supporting the investigation of gene expression. The total number of the Web-based resources integrated is 30. The database GeneNet on molecular events forming gene networks was assigned its integrative core. To navigate all these WWW-available resources, the SRS, HTML, and Java viewers were developed, http:@wwwmgs.bionet.nsc.ru/systems/GeneExpress/.

[WWWMGS: an integrated server for molecular-genetic studies]. / WWWMGS: integrirovannyi server dlia molekuliarno-geneticheskikh issledovanii.

We report an integrative technology for molecular biology studies in the field of transcription regulation by using Internet. A set of databases, programs, and systems are included into WWWMGS Web server. For example, the use of TRRD database information for site prediction is described. Using this method, the computer system SeqAnn was developed. The system performs the "real time" searching for prediction of initiation transcription site position according to database information. WWWMGS is available at URL: http://wwwmgs.bionet.nsc.ru/.

[Convergent origin of repeats in genes coding for globular proteins. An analysis of the factors determining the presence of inverted and symmetrical repeats]. / Konvergentnoe proiskhozhdenie povtorov v genakh, kodiruiushchikh globuliarnye belki. Analiz faktorov, obuslavlivaiushchikh nalichie invertirovannykh i simmetrichnykh povtorov.

The factors, determining the presence of inverted and symmetrical repeats in genes coding for globular proteins, have been analysed. An interesting property of genetical code has been revealed in the analysis of symmetrical repeats: the pairs of symmetrical codons corresponded to pairs of amino acids with mostly similar physical-chemical parameters. This property may explain the presence of symmetrical repeats and palindromes only in genes coding for beta-structural proteins-polypeptides, where amino acids with similar physical-chemical properties occupy symmetrical positions. A stochastic model of evolution of polynucleotide sequences has been used for analysis of inverted repeats. The modelling demonstrated that only limiting of sequences (uneven frequencies of used codons) is enough for arising of nonrandom inverted repeats in genes.