Search

Aggregating transcriptomics data across hospitals can increase sensitivity and robustness of differential expression analyses, yielding deeper clinical insights. As data exchange is often restricted by privacy legislation, meta-analyses are frequently employed to pool local results. However, the accuracy might drop if class labels are inhomogeneously distributed among cohorts. Flimma ( https://exbio.wzw.tum.de/flimma/ ) addresses this issue by implementing the state-of-the-art workflow limma voom in a federated manner, i.e., patient data never leaves its source site. Flimma results are identical to those generated by limma voom on aggregated datasets even in imbalanced scenarios where meta-analysis approaches fail., (© 2021. The Author(s).)

In network and systems medicine, active module identification methods (AMIMs) are widely used for discovering candidate molecular disease mechanisms. To this end, AMIMs combine network analysis algorithms with molecular profiling data, most commonly, by projecting gene expression data onto generic protein-protein interaction (PPI) networks. Although active module identification has led to various novel insights into complex diseases, there is increasing awareness in the field that the combination of gene expression data and PPI network is problematic because up-to-date PPI networks have a very small diameter and are subject to both technical and literature bias. In this paper, we report the results of an extensive study where we analyzed for the first time whether widely used AMIMs really benefit from using PPI networks. Our results clearly show that, except for the recently proposed AMIM DOMINO, the tested AMIMs do not produce biologically more meaningful candidate disease modules on widely used PPI networks than on random networks with the same node degrees. AMIMs hence mainly learn from the node degrees and mostly fail to exploit the biological knowledge encoded in the edges of the PPI networks. This has far-reaching consequences for the field of active module identification. In particular, we suggest that novel algorithms are needed which overcome the degree bias of most existing AMIMs and/or work with customized, context-specific networks instead of generic PPI networks., (© The Author(s) 2021. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

The axolotl, known for its remarkable regenerative abilities, is an excellent model for studying regenerative therapies. Nevertheless, the precise molecular mechanisms governing its regenerative potential remain uncertain. In this study, we collected samples from axolotls of different ages, including 8‐year‐old individuals and 8‐month‐old juveniles, obtaining their blastemas 10 days after amputation. Subsequently, we conducted a transcriptomic analysis comparing our samples to a set of previously published experiments. Our analysis unveiled a distinctive transcriptional response in the blastema, characterized by differential gene expression associated with processes such as bone and tissue remodeling, transcriptional regulation, angiogenesis, and intercellular communication. To gain deeper insights, we compared these findings with those from aged axolotls that showed no signs of regeneration 10 days after amputation. We identified four genes—FSTL1, ADAMTS17, GPX7, and CTHRC1—that showed higher expression in regenerating tissue compared to aged axolotls. Further scrutiny, including structural and homology analysis, revealed that these genes are conserved across vertebrate species. Our discoveries point to a group of proteins relevant to tissue regeneration, with their conservation in vertebrates suggesting critical roles in development. These findings also propose a novel gene set involved in axolotl regeneration, laying a promising foundation for future investigations across vertebrates. [ABSTRACT FROM AUTHOR]

Architectural proteins are essential epigenetic regulators that play a critical role in organizing chromatin and controlling gene expression. CTCF (CCCTC-binding factor) is a key architectural protein responsible for maintaining the intricate 3D structure of chromatin. Because of its multivalent properties and plasticity to bind various sequences, CTCF is similar to a Swiss knife for genome organization. Despite the importance of this protein, its mechanisms of action are not fully elucidated. It has been hypothesized that its versatility is achieved through interaction with multiple partners, forming a complex network that regulates chromatin folding within the nucleus. In this review, we delve into CTCF's interactions with other molecules involved in epigenetic processes, particularly histone and DNA demethylases, as well as several long non-coding RNAs (lncRNAs) that are able to recruit CTCF. Our review highlights the importance of CTCF partners to shed light on chromatin regulation and pave the way for future exploration of the mechanisms that enable the finely-tuned role of CTCF as a master regulator of chromatin. [ABSTRACT FROM AUTHOR]

Motivation Cancer is one of the leading causes of death worldwide. Despite significant improvements in prevention and treatment, mortality remains high for many cancer types. Hence, innovative methods that use molecular data to stratify patients and identify biomarkers are needed. Promising biomarkers can also be inferred from competing endogenous RNA (ceRNA) networks that capture the gene–miRNA gene regulatory landscape. Thus far, the role of these biomarkers could only be studied globally but not in a sample-specific manner. To mitigate this, we introduce spongEffects, a novel method that infers subnetworks (or modules) from ceRNA networks and calculates patient- or sample-specific scores related to their regulatory activity. Results We show how spongEffects can be used for downstream interpretation and machine learning tasks such as tumor classification and for identifying subtype-specific regulatory interactions. In a concrete example of breast cancer subtype classification, we prioritize modules impacting the biology of the different subtypes. In summary, spongEffects prioritizes ceRNA modules as biomarkers and offers insights into the miRNA regulatory landscape. Notably, these module scores can be inferred from gene expression data alone and can thus be applied to cohorts where miRNA expression information is lacking. Availability and implementation https://bioconductor.org/packages/devel/bioc/html/SPONGE.html. [ABSTRACT FROM AUTHOR]

Background Eukaryotic gene expression is controlled by cis -regulatory elements (CREs), including promoters and enhancers, which are bound by transcription factors (TFs). Differential expression of TFs and their binding affinity at putative CREs determine tissue- and developmental-specific transcriptional activity. Consolidating genomic datasets can offer further insights into the accessibility of CREs, TF activity, and, thus, gene regulation. However, the integration and analysis of multimodal datasets are hampered by considerable technical challenges. While methods for highlighting differential TF activity from combined chromatin state data (e.g. chromatin immunoprecipitation [ChIP], ATAC, or DNase sequencing) and RNA sequencing data exist, they do not offer convenient usability, have limited support for large-scale data processing, and provide only minimal functionality for visually interpreting results. Results We developed TF-Prioritizer, an automated pipeline that prioritizes condition-specific TFs from multimodal data and generates an interactive web report. We demonstrated its potential by identifying known TFs along with their target genes, as well as previously unreported TFs active in lactating mouse mammary glands. Additionally, we studied a variety of ENCODE datasets for cell lines K562 and MCF-7, including 12 histone modification ChIP sequencing as well as ATAC and DNase sequencing datasets, where we observe and discuss assay-specific differences. Conclusion TF-Prioritizer accepts ATAC, DNase, or ChIP sequencing and RNA sequencing data as input and identifies TFs with differential activity, thus offering an understanding of genome-wide gene regulation, potential pathogenesis, and therapeutic targets in biomedical research. [ABSTRACT FROM AUTHOR]

Motivation Unsupervised learning approaches are frequently used to stratify patients into clinically relevant subgroups and to identify biomarkers such as disease-associated genes. However, clustering and biclustering techniques are oblivious to the functional relationship of genes and are thus not ideally suited to pinpoint molecular mechanisms along with patient subgroups. Results We developed the network-constrained biclustering approach Biclustering Constrained by Networks (BiCoN) which (i) restricts biclusters to functionally related genes connected in molecular interaction networks and (ii) maximizes the difference in gene expression between two subgroups of patients. This allows BiCoN to simultaneously pinpoint molecular mechanisms responsible for the patient grouping. Network-constrained clustering of genes makes BiCoN more robust to noise and batch effects than typical clustering and biclustering methods. BiCoN can faithfully reproduce known disease subtypes as well as novel, clinically relevant patient subgroups, as we could demonstrate using breast and lung cancer datasets. In summary, BiCoN is a novel systems medicine tool that combines several heuristic optimization strategies for robust disease mechanism extraction. BiCoN is well-documented and freely available as a python package or a web interface. Availability and implementation PyPI package: https://pypi.org/project/bicon. Web interface https://exbio.wzw.tum.de/bicon. Supplementary information Supplementary data are available at Bioinformatics online. [ABSTRACT FROM AUTHOR]

COVID-19 is an infection caused by SARS-CoV-2 (Severe Acute Respiratory Syndrome coronavirus 2), which has caused a global outbreak. Current research efforts are focused on the understanding of the molecular mechanisms involved in SARS-CoV-2 infection in order to propose drug-based therapeutic options. Transcriptional changes due to epigenetic regulation are key host cell responses to viral infection and have been studied in SARS-CoV and MERS-CoV; however, such changes are not fully described for SARS-CoV-2. In this study, we analyzed multiple transcriptomes obtained from cell lines infected with MERS-CoV, SARS-CoV, and SARS-CoV-2, and from COVID-19 patient-derived samples. Using integrative analyses of gene co-expression networks and de-novo pathway enrichment, we characterize different gene modules and protein pathways enriched with Transcription Factors or Epifactors relevant for SARS-CoV-2 infection. We identified EP300, MOV10, RELA, and TRIM25 as top candidates, and more than 60 additional proteins involved in the epigenetic response during viral infection that has therapeutic potential. Our results show that targeting the epigenetic machinery could be a feasible alternative to treat COVID-19. [ABSTRACT FROM AUTHOR]

Objective: Here, we applied a multi-omics approach (i) to examine molecular pathways related to de- and remyelination in multiple sclerosis (MS) lesions; and (ii) to translate these findings to the CSF proteome in order to identify molecules that are differentially expressed among MS subtypes. Methods: To relate differentially expressed genes in MS lesions to de- and remyelination, we compared transcriptome of MS lesions to transcriptome of cuprizone (CPZ)-induced de- and remyelination. Protein products of the overlapping orthologous genes were measured within the CSF by quantitative proteomics, parallel reaction monitoring (PRM). Differentially regulated proteins were correlated with molecular markers of inflammation by using MesoScale multiplex immunoassay. Expression kinetics of differentially regulated orthologous genes and proteins were examined in the CPZ model. Results: In the demyelinated and remyelinated corpus callosum, we detected 1239 differentially expressed genes; 91 orthologues were also differentially expressed in MS lesions. Pathway analysis of these orthologues suggested that the TYROBP (DAP12)-TREM2 pathway, TNF-receptor 1, CYBA and the proteasome subunit PSMB9 were related to de- and remyelination. We designed 129 peptides representing 51 orthologous proteins, measured them by PRM in 97 individual CSF, and compared their levels between relapsing (n = 40) and progressive MS (n = 57). Four proteins were differentially regulated among relapsing and progressive MS: tyrosine protein kinase receptor UFO (UFO), TIMP-1, apolipoprotein C-II (APOC2), and beta-2-microglobulin (B2M). The orthologous genes/proteins in the mouse brain peaked during acute remyelination. UFO, TIMP-1 and B2M levels correlated inversely with inflammation in the CSF (IL-6, MCP-1/CCL2, TARC/CCL17). APOC2 showed positive correlation with IL-2, IL-16 and eotaxin-3/CCL26. Conclusions: Pathology-based multi-omics identified four CSF markers that were differentially expressed in MS subtypes. Upregulated TIMP-1, UFO and B2M orthologues in relapsing MS were associated with reduced inflammation and reflected reparatory processes, in contrast to the upregulated orthologue APOC2 in progressive MS that reflected changes in lipid metabolism associated with increased inflammation. [ABSTRACT FROM AUTHOR]

Identifying groups of similar objects is a popular first step in biomedical data analysis, but it is error-prone and impossible to perform manually. Many computational methods have been developed to tackle this problem. Here we assessed 13 well-known methods using 24 data sets ranging from gene expression to protein domains. Performance was judged on the basis of 13 common cluster validity indices. We developed a clustering analysis platform, ClustEval (http://clusteval.mpi-inf.mpg.de), to promote streamlined evaluation, comparison and reproducibility of clustering results in the future. This allowed us to objectively evaluate the performance of all tools on all data sets with up to 1,000 different parameter sets each, resulting in a total of more than 4 million calculated cluster validity indices. We observed that there was no universal best performer, but on the basis of this wide-ranging comparison we were able to develop a short guideline for biomedical clustering tasks. ClustEval allows biomedical researchers to pick the appropriate tool for their data type and allows method developers to compare their tool to the state of the art. [ABSTRACT FROM AUTHOR]

Background: Recently, the research community has seen an influx of data relating to transcriptional regulatory interactions of Corynebacteria, organisms that are highly relevant to fields of systems biology, biotechnology, and human medicine. Information derived from DNA microarray experiments, computational predictions, and literature has opened the way for the graph-based analysis, visualization, and reconstruction of transcriptional regulatory networks across entire organisms. The reference database for corynebacterial gene regulatory networks CoryneRegNet provides methods for data storage and data exchange in a well-structured manner. Additional information on the model organism Escherichia coli K12 obtained from RegulonDB has been integrated. Generally, gene regulatory networks can be visualized as graphs by drawing directed edges between nodes, where a node represents a gene and an edge corresponds to a typed regulatory interaction. Cytoscape is an open-source software project whose aim is to provide graph-based visualization and analysis for biological networks. Its architecture allows the development and integration of user-made plugins to enhance core functionalities. Results: We introduce two novel plugins for the Cytoscape environment designed to enhance in silico studies of procaryotic transcriptional regulatory networks. Our plugins leverage the information from the cornyebacterial reference database CoryneRegNet with the graph analysis capabilities of Cytoscape. CoryneRegNetLoader queries the CoryneRegNet database to extract a gene regulatory network represented as a directed graph. Additional information is stored as node/ edge attributes within the graph. COMA facilitates consistency checks for gene expression studies given a gene regulatory network. COMA tests whether all gene expression levels correlate properly with the given network topology. Conclusion: The plugins facilitate in silico studies of procaryotic transcriptional gene regulation, particularly in Corynebacteria and E. coli, by combining the knowledge from the corynebacterial reference database with the graph analysis capabilities of Cytoscape, which is one of the most-widely used tools for biological network analyses. [ABSTRACT FROM AUTHOR]

Background: Detailed information on DNA-binding transcription factors (the key players in the regulation of gene expression) and on transcriptional regulatory interactions of microorganisms deduced from literature-derived knowledge, computer predictions and global DNA microarray hybridization experiments, has opened the way for the genome-wide analysis of transcriptional regulatory networks. The large-scale reconstruction of these networks allows the in silico analysis of cell behavior in response to changing environmental conditions. We previously published CoryneRegNet, an ontology-based data warehouse of corynebacterial transcription factors and regulatory networks. Initially, it was designed to provide methods for the analysis and visualization of the gene regulatory network of Corynebacterium glutamicum. Results: Now we introduce CoryneRegNet release 4.0, which integrates data on the gene regulatory networks of 4 corynebacteria, 2 mycobacteria and the model organism Escherichia coli K12. As the previous versions, CoryneRegNet provides a web-based user interface to access the database content, to allow various queries, and to support the reconstruction, analysis and visualization of regulatory networks at different hierarchical levels. In this article, we present the further improved database content of CoryneRegNet along with novel analysis features. The network visualization feature GraphVis now allows the inter-species comparisons of reconstructed gene regulatory networks and the projection of gene expression levels onto that networks. Therefore, we added stimulon data directly into the database, but also provide Web Service access to the DNA microarray analysis platform EMMA. Additionally, CoryneRegNet now provides a SOAP based Web Service server, which can easily be consumed by other bioinformatics software systems. Stimulons (imported from the database, or uploaded by the user) can be analyzed in the context of known transcriptional regulatory networks to predict putative contradictions or further gene regulatory interactions. Furthermore, it integrates protein clusters by means of heuristically solving the weighted graph cluster editing problem. In addition, it provides Web Service based access to up to date gene annotation data from GenDB. Conclusion: The release 4.0 of CoryneRegNet is a comprehensive system for the integrated analysis of procaryotic gene regulatory networks. It is a versatile systems biology platform to support the efficient and large-scale analysis of transcriptional regulation of gene expression in microorganisms. It is publicly available at http://www.CoryneRegNet.DE. [ABSTRACT FROM AUTHOR]

Background: The introduction of high-throughput genome sequencing and post-genome analysis technologies, e.g. DNA microarray approaches, has created the potential to unravel and scrutinize complex gene-regulatory networks on a large scale. The discovery of transcriptional regulatory interactions has become a major topic in modern functional genomics. Results: To facilitate the analysis of gene-regulatory networks, we have developed CoryneCenter, a web-based resource for the systematic integration and analysis of genome, transcriptome, and gene regulatory information for prokaryotes, especially corynebacteria. For this purpose, we extended and combined the following systems into a common platform: (1) GenDB, an open source genome annotation system, (2) EMMA, a MAGE compliant application for high-throughput transcriptome data storage and analysis, and (3) CoryneRegNet, an ontologybased data warehouse designed to facilitate the reconstruction and analysis of gene regulatory interactions. We demonstrate the potential of CoryneCenter by means of an application example. Using microarray hybridization data, we compare the gene expression of Corynebacterium glutamicum under acetate and glucose feeding conditions: Known regulatory networks are confirmed, but moreover CoryneCenter points out additional regulatory interactions. Conclusion: CoryneCenter provides more than the sum of its parts. Its novel analysis and visualization features significantly simplify the process of obtaining new biological insights into complex regulatory systems. Although the platform currently focusses on corynebacteria, the integrated tools are by no means restricted to these species, and the presented approach offers a general strategy for the analysis and verification of gene regulatory networks. CoryneCenter provides freely accessible projects with the underlying genome annotation, gene expression, and gene regulation data. The system is publicly available at http://www.CoryneCenter.de. [ABSTRACT FROM AUTHOR]

Background: The application of DNA microarray technology in post-genomic analysis of bacterial genome sequences has allowed the generation of huge amounts of data related to regulatory networks. This data along with literature-derived knowledge on regulation of gene expression has opened the way for genome-wide reconstruction of transcriptional regulatory networks. These large-scale reconstructions can be converted into in silico models of bacterial cells that allow a systematic analysis of network behavior in response to changing environmental conditions. Description: CoryneRegNet was designed to facilitate the genome-wide reconstruction of transcriptional regulatory networks of corynebacteria relevant in biotechnology and human medicine. During the import and integration process of data derived from experimental studies or literature knowledge CoryneRegNet generates links to genome annotations, to identified transcription factors and to the corresponding cis-regulatory elements. CoryneRegNet is based on a multi-layered, hierarchical and modular concept of transcriptional regulation and was implemented by using the relational database management system MySQL and an ontology-based data structure. Reconstructed regulatory networks can be visualized by using the yFiles JAVA graph library. As an application example of CoryneRegNet, we have reconstructed the global transcriptional regulation of a cellular module involved in SOS and stress response of corynebacteria. Conclusion: CoryneRegNet is an ontology-based data warehouse that allows a pertinent data management of regulatory interactions along with the genome-scale reconstruction of transcriptional regulatory networks. These models can further be combined with metabolic networks to build integrated models of cellular function including both metabolism and its transcriptional regulation. [ABSTRACT FROM AUTHOR]

Simple Summary: Patients with locally advanced rectal cancer have been treated with chemoradiotherapy followed by surgery, which results in variable therapy response. To date, there is a lack of established predictive biomarkers to distinguish responsive from non-responsive patients. Therefore, patients can be overtreated resulting in unnecessary toxicity, or therapy changes can occur later when the cancer is more aggressive. We used pre-treatment biopsies to evaluate changes in DNA methylation that could predict the chemoradiotherapy response. Data from an external dataset were used to confirm our findings. We identified and validated a classifier, composed of three candidates, that was able to distinguish responders from non-responders. The genomic context of our biomarkers was explored, giving evidence that they play a role in regulating gene expression. The biomarkers herein described can be easily evaluated in the clinical practice and can help to guide rectal cancer patients' treatment by identifying responders and non-responders. The treatment for locally advanced rectal carcinomas (LARC) is based on neoadjuvant chemoradiotherapy (nCRT) and surgery, which results in pathological complete response (pCR) in up to 30% of patients. Since epigenetic changes may influence response to therapy, we aimed to identify DNA methylation markers predictive of pCR in LARC patients treated with nCRT. We used high-throughput DNA methylation analysis of 32 treatment-naïve LARC biopsies and five normal rectal tissues to explore the predictive value of differentially methylated (DM) CpGs. External validation was carried out with The Cancer Genome Atlas-Rectal Adenocarcinoma (TCGA-READ 99 cases). A classifier based on three-CpGs DM (linked to OBSL1, GPR1, and INSIG1 genes) was able to discriminate pCR from incomplete responders with high sensitivity and specificity. The methylation levels of the selected CpGs confirmed the predictive value of our classifier in 77 LARCs evaluated by bisulfite pyrosequencing. Evaluation of external datasets (TCGA-READ, GSE81006, GSE75546, and GSE39958) reproduced our results. As the three CpGs were mapped near to regulatory elements, we performed an integrative analysis in regions associated with predicted cis-regulatory elements. A positive and inverse correlation between DNA methylation and gene expression was found in two CpGs. We propose a novel predictive tool based on three CpGs potentially useful for pretreatment screening of LARC patients and guide the selection of treatment modality. [ABSTRACT FROM AUTHOR]

The heterogeneity of multiple sclerosis is reflected by dynamic changes of different lesion types in the brain white matter (WM). To identify potential drivers of this process, we RNA-sequenced 73 WM areas from patients with progressive MS (PMS) and 25 control WM. Lesion endophenotypes were described by a computational systems medicine analysis combined with RNAscope, immunohistochemistry, and immunofluorescence. The signature of the normal-appearing WM (NAWM) was more similar to control WM than to lesions: one of the six upregulated genes in NAWM was CD26/DPP4 expressed by microglia. Chronic active lesions that become prominent in PMS had a signature that were different from all other lesion types, and were differentiated from them by two clusters of 62 differentially expressed genes (DEGs). An upcoming MS biomarker, CHI3L1 was among the top ten upregulated genes in chronic active lesions expressed by astrocytes in the rim. TGFβ-R2 was the central hub in a remyelination-related protein interaction network, and was expressed there by astrocytes. We used de novo networks enriched by unique DEGs to determine lesion-specific pathway regulation, i.e. cellular trafficking and activation in active lesions; healing and immune responses in remyelinating lesions characterized by the most heterogeneous immunoglobulin gene expression; coagulation and ion balance in inactive lesions; and metabolic changes in chronic active lesions. Because we found inverse differential regulation of particular genes among different lesion types, our data emphasize that omics related to MS lesions should be interpreted in the context of lesion pathology. Our data indicate that the impact of molecular pathways is substantially changing as different lesions develop. This was also reflected by the high number of unique DEGs that were more common than shared signatures. A special microglia subset characterized by CD26 may play a role in early lesion development, while astrocyte-derived TGFβ-R2 and TGFβ pathways may be drivers of repair in contrast to chronic tissue damage. The highly specific mechanistic signature of chronic active lesions indicates that as these lesions develop in PMS, the molecular changes are substantially skewed: the unique mitochondrial/metabolic changes and specific downregulation of molecules involved in tissue repair may reflect a stage of exhaustion. [ABSTRACT FROM AUTHOR]