Identification of unique rectal cancer-specific subtypes.

BACKGROUND: Existing colorectal cancer subtyping methods were generated without much consideration of potential differences in expression profiles between colon and rectal tissues. Moreover, locally advanced rectal cancers at resection often have received neoadjuvant chemoradiotherapy which likely has a significant impact on gene expression. METHODS: We collected mRNA expression profiles for rectal and colon cancer samples (n = 2121). We observed that (i) Consensus Molecular Subtyping (CMS) had a different prognosis in treatment-naïve rectal vs. colon cancers, and (ii) that neoadjuvant chemoradiotherapy exposure produced a strong shift in CMS subtypes in rectal cancers. We therefore clustered 182 untreated rectal cancers to find rectal cancer-specific subtypes (RSSs). RESULTS: We identified three robust subtypes. We observed that RSS1 had better, and RSS2 had worse disease-free survival. RSS1 showed high expression of MYC target genes and low activity of angiogenesis genes. RSS2 exhibited low regulatory T cell abundance, strong EMT and angiogenesis signalling, and high activation of TGF-ß, NF-κB, and TNF-α signalling. RSS3 was characterised by the deactivation of EGFR, MAPK and WNT pathways. CONCLUSIONS: We conclude that RSS subtyping allows for more accurate prognosis predictions in rectal cancers than CMS subtyping and provides new insight into targetable disease pathways within these subtypes.

BioUML-towards a universal research platform.

BioUML (https://www.biouml.org)-is a web-based integrated platform for systems biology and data analysis. It supports visual modelling and construction of hierarchical biological models that allow us to construct the most complex modular models of blood pressure regulation, skeletal muscle metabolism, COVID-19 epidemiology. BioUML has been integrated with git repositories where users can store their models and other data. We have also expanded the capabilities of BioUML for data analysis and visualization of biomedical data: (i) any programs and Jupyter kernels can be plugged into the BioUML platform using Docker technology; (ii) BioUML is integrated with the Galaxy and Galaxy Tool Shed; (iii) BioUML provides two-way integration with R and Python (Jupyter notebooks): scripts can be executed on the BioUML web pages, and BioUML functions can be called from scripts; (iv) using plug-in architecture, specialized viewers and editors can be added. For example, powerful genome browsers as well as viewers for molecular 3D structure are integrated in this way; (v) BioUML supports data analyses using workflows (own format, Galaxy, CWL, BPMN, nextFlow). Using these capabilities, we have initiated a new branch of the BioUML development-u-science-a universal scientific platform that can be configured for specific research requirements.

Comparative analysis of deeply phenotyped GBM cohorts of 'short-term' and 'long-term' survivors.

BACKGROUND: Glioblastoma (GBM) is an aggressive brain cancer that typically results in death in the first 15 months after diagnosis. There have been limited advances in finding new treatments for GBM. In this study, we investigated molecular differences between patients with extremely short (≤ 9 months, Short term survivors, STS) and long survival (≥ 36 months, Long term survivors, LTS). METHODS: Patients were selected from an in-house cohort (GLIOTRAIN-cohort), using defined inclusion criteria (Karnofsky score > 70; age < 70 years old; Stupp protocol as first line treatment, IDH wild type), and a multi-omic analysis of LTS and STS GBM samples was performed. RESULTS: Transcriptomic analysis of tumour samples identified cilium gene signatures as enriched in LTS. Moreover, Immunohistochemical analysis confirmed the presence of cilia in the tumours of LTS. Notably, reverse phase protein array analysis (RPPA) demonstrated increased phosphorylated GAB1 (Y627), SRC (Y527), BCL2 (S70) and RAF (S338) protein expression in STS compared to LTS. Next, we identified 25 unique master regulators (MR) and 13 transcription factors (TFs) belonging to ontologies of integrin signalling and cell cycle to be upregulated in STS. CONCLUSION: Overall, comparison of STS and LTS GBM patients, identifies novel biomarkers and potential actionable therapeutic targets for the management of GBM.

GTRD: an integrated view of transcription regulation.

The Gene Transcription Regulation Database (GTRD; http://gtrd.biouml.org/) contains uniformly annotated and processed NGS data related to gene transcription regulation: ChIP-seq, ChIP-exo, DNase-seq, MNase-seq, ATAC-seq and RNA-seq. With the latest release, the database has reached a new level of data integration. All cell types (cell lines and tissues) presented in the GTRD were arranged into a dictionary and linked with different ontologies (BRENDA, Cell Ontology, Uberon, Cellosaurus and Experimental Factor Ontology) and with related experiments in specialized databases on transcription regulation (FANTOM5, ENCODE and GTEx). The updated version of the GTRD provides an integrated view of transcription regulation through a dedicated web interface with advanced browsing and search capabilities, an integrated genome browser, and table reports by cell types, transcription factors, and genes of interest.

BioUML: an integrated environment for systems biology and collaborative analysis of biomedical data.

BioUML (homepage: http://www.biouml.org, main public server: https://ict.biouml.org) is a web-based integrated environment (platform) for systems biology and the analysis of biomedical data generated by omics technologies. The BioUML vision is to provide a computational platform to build virtual cell, virtual physiological human and virtual patient. BioUML spans a comprehensive range of capabilities, including access to biological databases, powerful tools for systems biology (visual modelling, simulation, parameters fitting and analyses), a genome browser, scripting (R, JavaScript) and a workflow engine. Due to integration with the Galaxy platform and R/Bioconductor, BioUML provides powerful possibilities for the analyses of omics data. The plug-in-based architecture allows the user to add new functionalities using plug-ins. To facilitate a user focus on a particular task or database, we have developed several predefined perspectives that display only those web interface elements that are needed for a specific task. To support collaborative work on scientific projects, there is a central authentication and authorization system (https://bio-store.org). The diagram editor enables several remote users to simultaneously edit diagrams.

Signaling Pathways Potentially Responsible for Foam Cell Formation: Cholesterol Accumulation or Inflammatory Response-What is First?

Accumulation of lipid-laden (foam) cells in the arterial wall is known to be the earliest step in the pathogenesis of atherosclerosis. There is almost no doubt that atherogenic modified low-density lipoproteins (LDL) are the main sources of accumulating lipids in foam cells. Atherogenic modified LDL are taken up by arterial cells, such as macrophages, pericytes, and smooth muscle cells in an unregulated manner bypassing the LDL receptor. The present study was conducted to reveal possible common mechanisms in the interaction of macrophages with associates of modified LDL and non-lipid latex particles of a similar size. To determine regulatory pathways that are potentially responsible for cholesterol accumulation in human macrophages after the exposure to naturally occurring atherogenic or artificially modified LDL, we used transcriptome analysis. Previous studies of our group demonstrated that any type of LDL modification facilitates the self-association of lipoprotein particles. The size of such self-associates hinders their interaction with a specific LDL receptor. As a result, self-associates are taken up by nonspecific phagocytosis bypassing the LDL receptor. That is why we used latex beads as a stimulator of macrophage phagocytotic activity. We revealed at least 12 signaling pathways that were regulated by the interaction of macrophages with the multiple-modified atherogenic naturally occurring LDL and with latex beads in a similar manner. Therefore, modified LDL was shown to stimulate phagocytosis through the upregulation of certain genes. We have identified at least three genes (F2RL1, EIF2AK3, and IL15) encoding inflammatory molecules and associated with signaling pathways that were upregulated in response to the interaction of modified LDL with macrophages. Knockdown of two of these genes, EIF2AK3 and IL15, completely suppressed cholesterol accumulation in macrophages. Correspondingly, the upregulation of EIF2AK3 and IL15 promoted cholesterol accumulation. These data confirmed our hypothesis of the following chain of events in atherosclerosis: LDL particles undergo atherogenic modification; this is accompanied by the formation of self-associates; large LDL associates stimulate phagocytosis; as a result of phagocytosis stimulation, pro-inflammatory molecules are secreted; these molecules cause or at least contribute to the accumulation of intracellular cholesterol. This chain of events may explain the relationship between cholesterol accumulation and inflammation. The primary sequence of events in this chain is related to inflammatory response rather than cholesterol accumulation.

Role of Phagocytosis in the Pro-Inflammatory Response in LDL-Induced Foam Cell Formation; a Transcriptome Analysis.

Excessive accumulation of lipid inclusions in the arterial wall cells (foam cell formation) caused by modified low-density lipoprotein (LDL) is the earliest and most noticeable manifestation of atherosclerosis. The mechanisms of foam cell formation are not fully understood and can involve altered lipid uptake, impaired lipid metabolism, or both. Recently, we have identified the top 10 master regulators that were involved in the accumulation of cholesterol in cultured macrophages induced by the incubation with modified LDL. It was found that most of the identified master regulators were related to the regulation of the inflammatory immune response, but not to lipid metabolism. A possible explanation for this unexpected result is a stimulation of the phagocytic activity of macrophages by modified LDL particle associates that have a relatively large size. In the current study, we investigated gene regulation in macrophages using transcriptome analysis to test the hypothesis that the primary event occurring upon the interaction of modified LDL and macrophages is the stimulation of phagocytosis, which subsequently triggers the pro-inflammatory immune response. We identified genes that were up- or downregulated following the exposure of cultured cells to modified LDL or latex beads (inert phagocytosis stimulators). Most of the identified master regulators were involved in the innate immune response, and some of them were encoding major pro-inflammatory proteins. The obtained results indicated that pro-inflammatory response to phagocytosis stimulation precedes the accumulation of intracellular lipids and possibly contributes to the formation of foam cells. In this way, the currently recognized hypothesis that the accumulation of lipids triggers the pro-inflammatory response was not confirmed. Comparative analysis of master regulators revealed similarities in the genetic regulation of the interaction of macrophages with naturally occurring LDL and desialylated LDL. Oxidized and desialylated LDL affected a different spectrum of genes than naturally occurring LDL. These observations suggest that desialylation is the most important modification of LDL occurring in vivo. Thus, modified LDL caused the gene regulation characteristic of the stimulation of phagocytosis. Additionally, the knock-down effect of five master regulators, such as IL15, EIF2AK3, F2RL1, TSPYL2, and ANXA1, on intracellular lipid accumulation was tested. We knocked down these genes in primary macrophages derived from human monocytes. The addition of atherogenic naturally occurring LDL caused a significant accumulation of cholesterol in the control cells. The knock-down of the EIF2AK3 and IL15 genes completely prevented cholesterol accumulation in cultured macrophages. The knock-down of the ANXA1 gene caused a further decrease in cholesterol content in cultured macrophages. At the same time, knock-down of F2RL1 and TSPYL2 did not cause an effect. The results obtained allowed us to explain in which way the inflammatory response and the accumulation of cholesterol are related confirming our hypothesis of atherogenesis development based on the following viewpoints: LDL particles undergo atherogenic modifications that, in turn, accompanied by the formation of self-associates; large LDL associates stimulate phagocytosis; as a result of phagocytosis stimulation, pro-inflammatory molecules are secreted; these molecules cause or at least contribute to the accumulation of intracellular cholesterol. Therefore, it became obvious that the primary event in this sequence is not the accumulation of cholesterol but an inflammatory response.

Walking pathways with positive feedback loops reveal DNA methylation biomarkers of colorectal cancer.

BACKGROUND: The search for molecular biomarkers of early-onset colorectal cancer (CRC) is an important but still quite challenging and unsolved task. Detection of CpG methylation in human DNA obtained from blood or stool has been proposed as a promising approach to a noninvasive early diagnosis of CRC. Thousands of abnormally methylated CpG positions in CRC genomes are often located in non-coding parts of genes. Novel bioinformatic methods are thus urgently needed for multi-omics data analysis to reveal causative biomarkers with a potential driver role in early stages of cancer. METHODS: We have developed a method for finding potential causal relationships between epigenetic changes (DNA methylations) in gene regulatory regions that affect transcription factor binding sites (TFBS) and gene expression changes. This method also considers the topology of the involved signal transduction pathways and searches for positive feedback loops that may cause the carcinogenic aberrations in gene expression. We call this method "Walking pathways", since it searches for potential rewiring mechanisms in cancer pathways due to dynamic changes in the DNA methylation status of important gene regulatory regions ("epigenomic walking"). RESULTS: In this paper, we analysed an extensive collection of full genome gene-expression data (RNA-seq) and DNA methylation data of genomic CpG islands (using Illumina methylation arrays) generated from a sample of tumor and normal gut epithelial tissues of 300 patients with colorectal cancer (at different stages of the disease) (data generated in the EU-supported SysCol project). Identification of potential epigenetic biomarkers of DNA methylation was performed using the fully automatic multi-omics analysis web service "My Genome Enhancer" (MGE) (my-genome-enhancer.com). MGE uses the database on gene regulation TRANSFAC®, the signal transduction pathways database TRANSPATH®, and software that employs AI (artificial intelligence) methods for the analysis of cancer-specific enhancers. CONCLUSIONS: The identified biomarkers underwent experimental testing on an independent set of blood samples from patients with colorectal cancer. As a result, using advanced methods of statistics and machine learning, a minimum set of 6 biomarkers was selected, which together achieve the best cancer detection potential. The markers include hypermethylated positions in regulatory regions of the following genes: CALCA, ENO1, MYC, PDX1, TCF7, ZNF43.

GTRD: a database of transcription factor binding sites identified by ChIP-seq experiments.

GTRD-Gene Transcription Regulation Database (http://gtrd.biouml.org)-is a database of transcription factor binding sites (TFBSs) identified by ChIP-seq experiments for human and mouse. Raw ChIP-seq data were obtained from ENCODE and SRA and uniformly processed: (i) reads were aligned using Bowtie2; (ii) ChIP-seq peaks were called using peak callers MACS, SISSRs, GEM and PICS; (iii) peaks for the same factor and peak callers, but different experiment conditions (cell line, treatment, etc.), were merged into clusters; (iv) such clusters for different peak callers were merged into metaclusters that were considered as non-redundant sets of TFBSs. In addition to information on location in genome, the sets contain structured information about cell lines and experimental conditions extracted from descriptions of corresponding ChIP-seq experiments. A web interface to access GTRD was developed using the BioUML platform. It provides: (i) browsing and displaying information; (ii) advanced search possibilities, e.g. search of TFBSs near the specified gene or search of all genes potentially regulated by a specified transcription factor; (iii) integrated genome browser that provides visualization of the GTRD data: read alignments, peaks, clusters, metaclusters and information about gene structures from the Ensembl database and binding sites predicted using position weight matrices from the HOCOMOCO database.

Integrated Bio-Search: challenges and trends for the integration, search and comprehensive processing of biological information.

Many efforts exist to design and implement approaches and tools for data capture, integration and analysis in the life sciences. Challenges are not only the heterogeneity, size and distribution of information sources, but also the danger of producing too many solutions for the same problem. Methodological, technological, infrastructural and social aspects appear to be essential for the development of a new generation of best practices and tools. In this paper, we analyse and discuss these aspects from different perspectives, by extending some of the ideas that arose during the NETTAB 2012 Workshop, making reference especially to the European context. First, relevance of using data and software models for the management and analysis of biological data is stressed. Second, some of the most relevant community achievements of the recent years, which should be taken as a starting point for future efforts in this research domain, are presented. Third, some of the main outstanding issues, challenges and trends are analysed. The challenges related to the tendency to fund and create large scale international research infrastructures and public-private partnerships in order to address the complex challenges of data intensive science are especially discussed. The needs and opportunities of Genomic Computing (the integration, search and display of genomic information at a very specific level, e.g. at the level of a single DNA region) are then considered. In the current data and network-driven era, social aspects can become crucial bottlenecks. How these may best be tackled to unleash the technical abilities for effective data integration and validation efforts is then discussed. Especially the apparent lack of incentives for already overwhelmed researchers appears to be a limitation for sharing information and knowledge with other scientists. We point out as well how the bioinformatics market is growing at an unprecedented speed due to the impact that new powerful in silico analysis promises to have on better diagnosis, prognosis, drug discovery and treatment, towards personalized medicine. An open business model for bioinformatics, which appears to be able to reduce undue duplication of efforts and support the increased reuse of valuable data sets, tools and platforms, is finally discussed.

A discriminative approach for unsupervised clustering of DNA sequence motifs.

Algorithmic comparison of DNA sequence motifs is a problem in bioinformatics that has received increased attention during the last years. Its main applications concern characterization of potentially novel motifs and clustering of a motif collection in order to remove redundancy. Despite growing interest in motif clustering, the question which motif clusters to aim at has so far not been systematically addressed. Here we analyzed motif similarities in a comprehensive set of vertebrate transcription factor classes. For this we developed enhanced similarity scores by inclusion of the information coverage (IC) criterion, which evaluates the fraction of information an alignment covers in aligned motifs. A network-based method enabled us to identify motif clusters with high correspondence to DNA-binding domain phylogenies and prior experimental findings. Based on this analysis we derived a set of motif families representing distinct binding specificities. These motif families were used to train a classifier which was further integrated into a novel algorithm for unsupervised motif clustering. Application of the new algorithm demonstrated its superiority to previously published methods and its ability to reproduce entrained motif families. As a result, our work proposes a probabilistic approach to decide whether two motifs represent common or distinct binding specificities.

MatrixCatch--a novel tool for the recognition of composite regulatory elements in promoters.

BACKGROUND: Accurate recognition of regulatory elements in promoters is an essential prerequisite for understanding the mechanisms of gene regulation at the level of transcription. Composite regulatory elements represent a particular type of such transcriptional regulatory elements consisting of pairs of individual DNA motifs. In contrast to the present approach, most available recognition techniques are based purely on statistical evaluation of the occurrence of single motifs. Such methods are limited in application, since the accuracy of recognition is greatly dependent on the size and quality of the sequence dataset. Methods that exploit available knowledge and have broad applicability are evidently needed. RESULTS: We developed a novel method to identify composite regulatory elements in promoters using a library of known examples. In depth investigation of regularities encoded in known composite elements allowed us to introduce a new characteristic measure and to improve the specificity compared with other methods. Tests on an established benchmark and real genomic data show that our method outperforms other available methods based either on known examples or statistical evaluations. In addition to better recognition, a practical advantage of this method is first the ability to detect a high number of different types of composite elements, and second direct biological interpretation of the identified results. The program is available at http://gnaweb.helmholtz-hzi.de/cgi-bin/MCatch/MatrixCatch.pl and includes an option to extend the provided library by user supplied data. CONCLUSIONS: The novel algorithm for the identification of composite regulatory elements presented in this paper was proved to be superior to existing methods. Its application to tissue specific promoters identified several highly specific composite elements with relevance to their biological function. This approach together with other methods will further advance the understanding of transcriptional regulation of genes.

Epigenome-Wide Changes in the Cell Layers of the Vein Wall When Exposing the Venous Endothelium to Oscillatory Shear Stress.

Epigenomic changes in the venous cells exerted by oscillatory shear stress towards the endothelium may result in consolidation of gene expression alterations upon vein wall remodeling during varicose transformation. We aimed to reveal such epigenome-wide methylation changes. Primary culture cells were obtained from non-varicose vein segments left after surgery of 3 patients by growing the cells in selective media after magnetic immunosorting. Endothelial cells were either exposed to oscillatory shear stress or left at the static condition. Then, other cell types were treated with preconditioned media from the adjacent layer's cells. DNA isolated from the harvested cells was subjected to epigenome-wide study using Illumina microarrays followed by data analysis with GenomeStudio (Illumina), Excel (Microsoft), and Genome Enhancer (geneXplain) software packages. Differential (hypo-/hyper-) methylation was revealed for each cell layer's DNA. The most targetable master regulators controlling the activity of certain transcription factors regulating the genes near the differentially methylated sites appeared to be the following: (1) HGS, PDGFB, and AR for endothelial cells; (2) HGS, CDH2, SPRY2, SMAD2, ZFYVE9, and P2RY1 for smooth muscle cells; and (3) WWOX, F8, IGF2R, NFKB1, RELA, SOCS1, and FXN for fibroblasts. Some of the identified master regulators may serve as promising druggable targets for treating varicose veins in the future.

Complexation of histone deacetylase inhibitor belinostat to Cu(II) prevents premature metabolic inactivation in vitro and demonstrates potent anti-cancer activity in vitro and ex vivo in colon cancer.

PURPOSE: The histone deacetylase inhibitor (HDACi), belinostat, has had limited therapeutic impact in solid tumors, such as colon cancer, due to its poor metabolic stability. Here we evaluated a novel belinostat prodrug, copper-bis-belinostat (Cubisbel), in vitro and ex vivo, designed to overcome the pharmacokinetic challenges of belinostat. METHODS: The in vitro metabolism of each HDACi was evaluated in human liver microsomes (HLMs) using mass spectrometry. Next, the effect of belinostat and Cubisbel on cell growth, HDAC activity, apoptosis and cell cycle was assessed in three colon cancer cell lines. Gene expression alterations induced by both HDACis were determined using RNA-Seq, followed by in silico analysis to identify master regulators (MRs) of differentially expressed genes (DEGs). The effect of both HDACis on the viability of colon cancer patient-derived tumor organoids (PDTOs) was also examined. RESULTS: Belinostat and Cubisbel significantly reduced colon cancer cell growth mediated through HDAC inhibition and apoptosis induction. Interestingly, the in vitro half-life of Cubisbel was significantly longer than belinostat. Belinostat and its Cu derivative commonly dysregulated numerous signalling and metabolic pathways while genes downregulated by Cubisbel were potentially controlled by VEGFA, ERBB2 and DUSP2 MRs. Treatment of colon cancer PDTOs with the HDACis resulted in a significant reduction in cell viability and downregulation of stem cell and proliferation markers. CONCLUSIONS: Complexation of belinostat to Cu(II) does not alter the HDAC activity of belinostat, but instead significantly enhances its metabolic stability in vitro and targets anti-cancer pathways by perturbing key MRs in colon cancer. Complexation of HDACis to a metal ion might improve the efficacy of clinically used HDACis in patients with colon cancer.

Inhibition of glycolytic enzymes mediated by pharmacologically activated p53: targeting Warburg effect to fight cancer.

Unique sensitivity of tumor cells to the inhibition of glycolysis is a good target for anticancer therapy. Here, we demonstrate that the pharmacologically activated tumor suppressor p53 mediates the inhibition of glycolytic enzymes in cancer cells in vitro and in vivo. We showed that p53 binds to the promoters of metabolic genes and represses their expression, including glucose transporters SLC2A12 (GLUT12) and SLC2A1 (GLUT1). Furthermore, p53-mediated repression of transcription factors c-Myc and HIF1α, key drivers of ATP-generating pathways in tumors, contributed to ATP production block. Inhibition of c-Myc by p53 mediated the ablation of several glycolytic genes in normoxia, whereas in hypoxia down-regulation of HIF1α contributed to this effect. We identified Sp1 as a transcription cofactor cooperating with p53 in the ablation of metabolic genes. Using different approaches, we demonstrated that glycolysis block contributes to the robust induction of apoptosis by p53 in cancer cells. Taken together, our data suggest that tumor-specific reinstatement of p53 function targets the "Achilles heel" of cancer cells (i.e. their dependence on glycolysis), which could contribute to the tumor-selective killing of cancer cells by pharmacologically activated p53.

Mild mitochondrial uncoupling does not affect mitochondrial biogenesis but downregulates pyruvate carboxylase in adipocytes: role for triglyceride content reduction.

In adipocytes, mitochondrial uncoupling is known to trigger a triglyceride loss comparable with the one induced by TNFα, a proinflammatory cytokine. However, the impact of a mitochondrial uncoupling on the abundance/composition of mitochondria and its connection with triglyceride content in adipocytes is largely unknown. In this work, the effects of a mild mitochondrial uncoupling triggered by FCCP were investigated on the mitochondrial population of 3T3-L1 adipocytes by both quantitative and qualitative approaches. We found that mild mitochondrial uncoupling does not stimulate mitochondrial biogenesis in adipocytes but induces an adaptive cell response characterized by quantitative modifications of mitochondrial protein content. Superoxide anion radical level was increased in mitochondria of both TNFα- and FCCP-treated adipocytes, whereas mitochondrial DNA copy number was significantly higher only in TNFα-treated cells. Subproteomic analysis revealed that the abundance of pyruvate carboxylase was reduced significantly in mitochondria of TNFα- and FCCP-treated adipocytes. Functional study showed that overexpression of this major enzyme of lipid metabolism is able to prevent the triglyceride content reduction in adipocytes exposed to mitochondrial uncoupling or TNFα. These results suggest a new mechanism by which the effects of mitochondrial uncoupling might limit triglyceride accumulation in adipocytes.

Transcriptome Analysis Reveals Intrinsic Proinflammatory Signaling in Healthy African American Skin.

Differences in the morphology and physiology of darkly pigmented skin compared with those of lightly pigmented skin are well-recognized. There are also disparities in the prevalence and clinical features for many inflammatory skin diseases, including atopic dermatitis and psoriasis; however, the underlying mechanisms are largely unknown. We compared the baseline gene expression in full-thickness skin biopsies from healthy individuals self-reporting as African American (AA) or as White non-Hispanic (WNH). Extensively validated RNA-sequencing analysis identified 570 differentially expressed genes in AA skin, including Igs and their receptors such as FCER1G; proinflammatory genes such as TNFα and IL32; and epidermal differentiation cluster and keratin genes. Differentially expressed genes were functionally enriched for inflammatory responses, keratinization, and cornified envelope formation. RNA-sequencing analysis of three-dimensional human skin equivalents made from AA and WNH primary keratinocytes revealed 360 differentially expressed genes (some shared with skin) that were enriched by similar functions. AA human skin equivalents appeared more responsive to TNF-α proinflammatory effects. Finally, AA-specific differentially expressed genes in the skin and human skin equivalents significantly overlapped with molecular signatures of skin in patients with atopic dermatitis and psoriasis. Overall, these findings suggest the existence of intrinsic proinflammatory circuits in AA keratinocytes/skin that may account for disease disparities and will help to build a foundation for the development of targeted skin disease prevention.

Master Transcription Regulators and Transcription Factors Regulate Immune-Associated Differences Between Patients of African and European Ancestry With Colorectal Cancer.

Background and Aims: Individuals of African (AFR) ancestry have a higher incidence of colorectal cancer (CRC) than those of European (EUR) ancestry and exhibit significant health disparities. Previous studies have noted differences in the tumor microenvironment between AFR and EUR patients with CRC. However, the molecular regulatory processes that underpin these immune differences remain largely unknown. Methods: Multiomics analysis was carried out for 55 AFR and 456 EUR patients with microsatellite-stable CRC using The Cancer Genome Atlas. We evaluated the tumor microenvironment by using gene expression and methylation data, transcription factor, and master transcriptional regulator analysis to identify the cell signaling pathways mediating the observed phenotypic differences. Results: We demonstrate that downregulated genes in AFR patients with CRC showed enrichment for canonical pathways, including chemokine signaling. Moreover, evaluation of the tumor microenvironment showed that cytotoxic lymphocytes and neutrophil cell populations are significantly decreased in AFR compared with EUR patients, suggesting AFR patients have an attenuated immune response. We further demonstrate that molecules called "master transcriptional regulators" (MTRs) play a critical role in regulating the expression of genes impacting key immune processes through an intricate signal transduction network mediated by disease-associated transcription factors (TFs). Furthermore, a core set of these MTRs and TFs showed a positive correlation with levels of cytotoxic lymphocytes and neutrophils across both AFR and EUR patients with CRC, thus suggesting their role in driving the immune infiltrate differences between the two ancestral groups. Conclusion: Our study provides an insight into the intricate regulatory landscape of MTRs and TFs that orchestrate the differences in the tumor microenvironment between patients with CRC of AFR and EUR ancestry.

Human-chimpanzee promoter comparisons: property-conserved evolution?

Identification of different functional elements and their properties is a fundamental need in biomedical research and phylogenetic comparisons of a growing number of sequenced genomes form a solid basis for this task. Most available phylogenetic approaches are focused on searching for individual sequence alterations, responsible for the observed phenotype, or statistically evaluate observed mutations to infer general trends. However, being applied to close genomes such methods suffer from poor statistics of rare mutations and give only (at its best) coarse results concerning the potential functional importance of the nucleotide differences. However, quantifying the changes in physical properties of DNA allows to see the strength of introduced mutations and hence to classify them for further investigations. In this work we present the comparative sequence analysis of two evolutionarily close species-human and chimpanzee. In contrast to previous studies we evaluate changes in melting enthalpy of DNA rather than count nucleotide mismatches. We find that nucleotide mismatches in promoters were apparently introduced in a correlated manner during the course of evolution, so that, for example, the DNA property "melting enthalpy" was retained. Such property conservation of promoters is significantly different from nucleotide conservation, shows significant positional and functional biases, and seems to represent a novel feature of gene regulation.

IGFBP2 Is a Potential Master Regulator Driving the Dysregulated Gene Network Responsible for Short Survival in Glioblastoma Multiforme.

Only 2% of glioblastoma multiforme (GBM) patients respond to standard therapy and survive beyond 36 months (long-term survivors, LTS), while the majority survive less than 12 months (short-term survivors, STS). To understand the mechanism leading to poor survival, we analyzed publicly available datasets of 113 STS and 58 LTS. This analysis revealed 198 differentially expressed genes (DEGs) that characterize aggressive tumor growth and may be responsible for the poor prognosis. These genes belong largely to the Gene Ontology (GO) categories "epithelial-to-mesenchymal transition" and "response to hypoxia." In this article, we applied an upstream analysis approach that involves state-of-the-art promoter analysis and network analysis of the dysregulated genes potentially responsible for short survival in GBM. Binding sites for transcription factors (TFs) associated with GBM pathology like NANOG, NF-κB, REST, FRA-1, PPARG, and seven others were found enriched in the promoters of the dysregulated genes. We reconstructed the gene regulatory network with several positive feedback loops controlled by five master regulators [insulin-like growth factor binding protein 2 (IGFBP2), vascular endothelial growth factor A (VEGFA), VEGF165, platelet-derived growth factor A (PDGFA), adipocyte enhancer-binding protein (AEBP1), and oncostatin M (OSMR)], which can be proposed as biomarkers and as therapeutic targets for enhancing GBM prognosis. A critical analysis of this gene regulatory network gives insights into the mechanism of gene regulation by IGFBP2 via several TFs including the key molecule of GBM tumor invasiveness and progression, FRA-1. All the observations were validated in independent cohorts, and their impact on overall survival has been investigated.