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.

Equal impact of diffusion and DNA binding rates on the potential spatial distribution of nuclear factor κB transcription factor inside the nucleus.

There are two physical processes that influence the spatial distribution of transcription factor molecules entering the nucleus of a eukaryotic cell, the binding to genomic DNA and the diffusion throughout the nuclear volume. Comparison of the DNA-protein association rate constant and the protein diffusion constant may determine which one is the limiting factor. If the process is diffusion-limited, transcription factor molecules are captured by DNA before their even distribution in the nuclear volume. Otherwise, if the reaction rate is limiting, these molecules diffuse evenly and then find their binding sites. Using well-studied human NF-κB dimer as an example, we calculated its diffusion constant using the Debye-Smoluchowski equation. The value of diffusion constant was about 10(-15) cm(3)/s, and it was comparable to the NF-κB association rate constant for DNA binding known from previous studies. Thus, both diffusion and DNA binding play an equally important role in NF-κB spatial distribution. The importance of genome 3D-structure in gene expression regulation and possible dependence of gene expression on the local concentration of open chromatin can be hypothesized from our theoretical estimate.

3DTF: a web server for predicting transcription factor PWMs using 3D structure-based energy calculations.

We present the webserver 3D transcription factor (3DTF) to compute position-specific weight matrices (PWMs) of transcription factors using a knowledge-based statistical potential derived from crystallographic data on protein-DNA complexes. Analysis of available structures that can be used to construct PWMs shows that there are hundreds of 3D structures from which PWMs could be derived, as well as thousands of proteins homologous to these. Therefore, we created 3DTF, which delivers binding matrices given the experimental or modeled protein-DNA complex. The webserver can be used by biologists to derive novel PWMs for transcription factors lacking known binding sites and is freely accessible at http://www.gene-regulation.com/pub/programs/3dtf/.

[Master regulators associated with poor prognosis in glioblastoma multiforme]. / Master-reguliatory, sviazannye s plokhim prognozom pri glioblastome.

Glioblastoma multiforme (GBM) is a highly malignant brain tumor with average survival time of 15 months. Less than 2% of the patients survive beyond 36 months. To understand the molecular mechanism responsible for poor prognosis, we analyzed GBM samples of TCGA microarray (n=560) data. We have identified 720 genes that have a significant impact upon survival based on univariate cox regression. We applied the Genome Enhancer pipeline to analyze potential mechanisms of regulation of activity of these genes and to build gene regulatory networks. We identified 12 transcription factors enriched in the promoters of these genes including the key molecule of GBM - STAT3. We found that STAT3 had significant differential expression across extreme survivor groups (short-term survivors- survival 36 months) and also had a significant impact on survival. In the next step, we identified master regulators in the signal transduction network that regulate the activity of these transcription factors. Master regulators are filtered based on their differential expression across extreme survivors groups and impact on survival. This work validates our earlier report on master regulators IGFBP2, PDGFA, OSMR, and AEBP1 driving short survival. Additionally, we propose CD14, CD44, DUSP6, GRB10, IL1RAP, FGFR3, and POSTN as master regulators driving poor survival. These master regulators are proposed as promising therapeutic targets to counter poor prognosis in GBM. Finally, the algorithm has prioritized several drugs for the further study as potential remedies to conquer the aggressive forms of GBM and to extend survival of the patients.

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

Towards a European health research and innovation cloud (HRIC).

The European Union (EU) initiative on the Digital Transformation of Health and Care (Digicare) aims to provide the conditions necessary for building a secure, flexible, and decentralized digital health infrastructure. Creating a European Health Research and Innovation Cloud (HRIC) within this environment should enable data sharing and analysis for health research across the EU, in compliance with data protection legislation while preserving the full trust of the participants. Such a HRIC should learn from and build on existing data infrastructures, integrate best practices, and focus on the concrete needs of the community in terms of technologies, governance, management, regulation, and ethics requirements. Here, we describe the vision and expected benefits of digital data sharing in health research activities and present a roadmap that fosters the opportunities while answering the challenges of implementing a HRIC. For this, we put forward five specific recommendations and action points to ensure that a European HRIC: i) is built on established standards and guidelines, providing cloud technologies through an open and decentralized infrastructure; ii) is developed and certified to the highest standards of interoperability and data security that can be trusted by all stakeholders; iii) is supported by a robust ethical and legal framework that is compliant with the EU General Data Protection Regulation (GDPR); iv) establishes a proper environment for the training of new generations of data and medical scientists; and v) stimulates research and innovation in transnational collaborations through public and private initiatives and partnerships funded by the EU through Horizon 2020 and Horizon Europe.

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

Composite Module Analyst: identification of transcription factor binding site combinations using genetic algorithm.

Composite Module Analyst (CMA) is a novel software tool aiming to identify promoter-enhancer models based on the composition of transcription factor (TF) binding sites and their pairs. CMA is closely interconnected with the TRANSFAC database. In particular, CMA uses the positional weight matrix (PWM) library collected in TRANSFAC and therefore provides the possibility to search for a large variety of different TF binding sites. We model the structure of the long gene regulatory regions by a Boolean function that joins several local modules, each consisting of co-localized TF binding sites. Having as an input a set of co-regulated genes, CMA builds the promoter model and optimizes the parameters of the model automatically by applying a genetic-regression algorithm. We use a multicomponent fitness function of the algorithm which includes several statistical criteria in a weighted linear function. We show examples of successful application of CMA to a microarray data on transcription profiling of TNF-alpha stimulated primary human endothelial cells. The CMA web server is freely accessible at http://www.gene-regulation.com/pub/programs/cma/CMA.html. An advanced version of CMA is also a part of the commercial system ExPlaintrade mark (www.biobase.de) designed for causal analysis of gene expression data.

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.

c-Jun N-terminal kinase 2 promotes enterocyte survival and goblet cell differentiation in the inflamed intestine.

c-Jun N-terminal kinases (JNKs) contribute to immune signaling but their functional role during intestinal mucosal inflammation has remained ill defined. Using genetic mouse models, we characterized the role of JNK1 and JNK2 during homeostasis and acute colitis. Epithelial apoptosis, regeneration, differentiation, and barrier function were analyzed in intestinal epithelium-specific (ΔIEC) or complete JNK1 and bone marrow chimeric or complete JNK2 deficient mice as well as double-knockout animals (JNK1ΔIECJNK2-/-) during homeostasis and acute dextran sulfate sodium (DSS)-induced colitis. Results were confirmed using human HT-29 cells and wild-type or JNK2-deficient mouse intestinal organoid cultures. We show that nonhematopoietic JNK2 but not JNK1 expression confers protection from DSS-induced intestinal inflammation reducing epithelial barrier dysfunction and enterocyte apoptosis. JNK2 additionally enhanced Atonal homolog 1 expression, goblet cell and enteroendocrine cell differentiation, and mucus production under inflammatory conditions. Our results identify a protective role of epithelial JNK2 signaling to maintain mucosal barrier function, epithelial cell integrity, and mucus layer production in the event of inflammatory tissue damage.

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.

Detecting reliable non interacting proteins (NIPs) significantly enhancing the computational prediction of protein-protein interactions using machine learning methods.

Protein-protein interactions (PPIs) play a vital role in most biological processes. Hence their comprehension can promote a better understanding of the mechanisms underlying living systems. However, besides the cost and the time limitation involved in the detection of experimentally validated PPIs, the noise in the data is still an important issue to overcome. In the last decade several in silico PPI prediction methods using both structural and genomic information were developed for this purpose. Here we introduce a unique validation approach aimed to collect reliable non interacting proteins (NIPs). Thereafter the most relevant protein/protein-pair related features were selected. Finally, the prepared dataset was used for PPI classification, leveraging the prediction capabilities of well-established machine learning methods. Our best classification procedure displayed specificity and sensitivity values of 96.33% and 98.02%, respectively, surpassing the prediction capabilities of other methods, including those trained on gold standard datasets. We showed that the PPI/NIP predictive performances can be considerably improved by focusing on data preparation.

Prediction of protein-DNA interactions of transcription factors linking proteomics and transcriptomics data.

We compared positional weight matrix-based prediction methods for transcription factor (TF) binding sites using selected fraction of ChIP-seq data with the help of partial AUC measure (limited to false positive rate 0.1, that is the most relevant for the application of the TF search in the genome scale). Comparison of three prediction methods-additive, multiplicative and information-vector based (MATCH) showed an advantage of the MATCH method for majority of transcription factors tested. We demonstrated that application of TF site identifying methods can help to connect the proteomics and phosphoproteomics world of signaling networks to gene regulation and transcriptomics world.

Recognition of NFATp/AP-1 composite elements within genes induced upon the activation of immune cells.

Composite elements are regulatory modules of promoters or enhancers that consist of binding sites of two different but synergizing transcription factors. A well-studied example is nuclear factors of activated T-cell (NFAT) sites which are composite elements of a NFATp/c and an activating protein 1 (AP-1) binding site. We have developed a computational approach to identify potential NFAT target genes which (a) comprises an improved method to scan for individual NFAT composite elements; (b) considers positional effects relative to transcription start sites; and (c) involves cluster analysis of potential NFAT composite elements. All three steps progressively helpX?ed to discriminate T-cell-specific promoter sequences against other functional regions (coding and intronic sequences) of the same genes, against promoters of muscle-specific genes or against random sequences. Using this approach, we identified potential NFAT composite elements in promoters of cytokine genes and their receptors as well as in promoters of genes for AP-1 family members, Ca2+-binding proteins and some other components of the regulatory network operating in activated T-cells and other immune cells. The method developed can be adapted to characterize and identify other composite elements as well. The program for recognition NFAT composite elements is available through the World Wide Web (http://compel.bionet.nsc.ru/FunSite/CompelScan. html and http://transfac.gbf.de/dbsearch/funsitep/s _comp.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.