Search

Abstract Mosquito-borne diseases are a major global health threat. Traditional morphological or molecular methods for identifying mosquito species often require specialized expertise or expensive laboratory equipment. The use of convolutional neural networks (CNNs) to identify mosquito species based on images may offer a promising alternative, but their practical implementation often remains limited. This study explores the applicability of CNNs in classifying mosquito species. It compares the efficacy of body and wing depictions across three image collection methods: a smartphone, macro-lens attached to a smartphone and a professional stereomicroscope. The study included 796 specimens of four morphologically similar Aedes species, Aedes aegypti, Ae. albopictus, Ae. koreicus and Ae. japonicus japonicus. The findings of this study indicate that CNN models demonstrate superior performance in wing-based classification 87.6% (95% CI: 84.2–91.0) compared to body-based classification 78.9% (95% CI: 77.7–80.0). Nevertheless, there are notable limitations of CNNs as they perform reliably across multiple devices only when trained specifically on those devices, resulting in an average decline of mean accuracy by 14%, even with extensive image augmentation. Additionally, we also estimate the required training data volume for effective classification, noting a reduced requirement for wing-based classification compared to body-based methods. Our study underscores the viability of both body and wing classification methods for mosquito species identification while emphasizing the need to address practical constraints in developing accessible classification systems. Graphical abstract

Abstract Genomic data holds huge potential for medical progress but requires strict safety measures due to its sensitive nature to comply with data protection laws. This conflict is especially pronounced in genome-wide association studies (GWAS) which rely on vast amounts of genomic data to improve medical diagnoses. To ensure both their benefits and sufficient data security, we propose a federated approach in combination with privacy-enhancing technologies utilising the findings from a systematic review on federated learning and legal regulations in general and applying these to GWAS.

Abstract Artificial intelligence (AI) has the potential to transform clinical practice and healthcare. Following impressive advancements in fields such as computer vision and medical imaging, AI is poised to drive changes in microbiome‐based healthcare while facing challenges specific to the field. This review describes the state‐of‐the‐art use of AI in microbiome‐related healthcare. It points out limitations across topics such as data handling, AI modelling and safeguarding patient privacy. Furthermore, we indicate how these current shortcomings could be overcome in the future and discuss the influence and opportunities of increasingly complex data on microbiome‐based healthcare.

Summary: The reproductive phase of plants is highly sensitive to ambient temperature stresses. To investigate sensitivity of female reproductive organs in grass crops during the pollination phase, we exposed the elongated stigma (silk) of maize to ambient environment at the silking stage. Moderate heat stress causes cell death of silk hair cells but did not affect early pollen tube growth inside the silk. Late pollen tube growth arrest was observed, leading to sterility. Heat stress causes elevated levels of reactive oxygen species (ROS) in silks, whose levels can be reduced by scavengers partly restoring pollen tube growth and fertility. A number of biological processes including hydrogen peroxide catabolic processes and bHLH transcription factor genes are downregulated by heat stress, while some NAC transcription factor genes are strongly upregulated. In conclusion, this study now provides a basis to select genes for engineering heat-stress-tolerant grass crops during the pollination phase.

Abstract Bulk RNA sequencing (RNA-seq) of blood is typically used for gene expression analysis in biomedical research but is still rarely used in clinical practice. In this study, we propose that RNA-seq should be considered a diagnostic tool, as it offers not only insights into aberrant gene expression and splicing but also delivers additional readouts on immune cell type composition as well as B-cell and T-cell receptor (BCR/TCR) repertoires. We demonstrate that RNA-seq offers insights into a patient’s immune status via integrative analysis of RNA-seq data from patients infected with various SARS-CoV-2 variants (in total 196 samples with up to 200 million reads sequencing depth). We compare the results of computational cell-type deconvolution methods (e.g., MCP-counter, xCell, EPIC, quanTIseq) to complete blood count data, the current gold standard in clinical practice. We observe varying levels of lymphocyte depletion and significant differences in neutrophil levels between SARS-CoV-2 variants. Additionally, we identify B and T cell receptor (BCR/TCR) sequences using the tools MiXCR and TRUST4 to show that—combined with sequence alignments and BLASTp—they could be used to classify a patient's disease. Finally, we investigated the sequencing depth required for such analyses and concluded that 10 million reads per sample is sufficient. In conclusion, our study reveals that computational cell-type deconvolution and BCR/TCR methods using bulk RNA-seq analyses can supplement missing CBC data and offer insights into immune responses, disease severity, and pathogen-specific immunity, all achievable with a sequencing depth of 10 million reads per sample.

Abstract High drug development costs and the limited number of new annual drug approvals increase the need for innovative approaches for drug effect prediction. Severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2), the cause of coronavirus disease 2019 (COVID‐19), led to a global pandemic with high morbidity and mortality. Although effective preventive measures exist, there are few effective treatments for hospitalized patients with SARS‐CoV‐2 infection. Drug repurposing and drug effect prediction are promising strategies that could shorten development time and reduce costs compared with de novo drug discovery. In this work, we present a machine learning framework to integrate a variety of target network features and physicochemical properties of compounds, and analyze their influence on the therapeutic effects for SARS‐CoV‐2 infection and on host cell cytotoxic effects. Random forest models trained on compounds with known experimental effects on SARS‐CoV‐2 infection and subsequent feature importance analysis based on Shapley values provided insights into the determinants of drug efficacy and cytotoxicity, which can be incorporated into novel drug discovery approaches. Given the complexity of molecular mechanisms of drug action and limited sample sizes, our models achieve a reasonable mean area under the receiver operating characteristic curve (ROC‐AUC) of 0.73 on an unseen validation set. To our knowledge, this is the first work to incorporate a combination of network and physicochemical features of compounds into a machine learning model to predict drug effects on SARS‐CoV‐2 infection. Our systems pharmacology–based machine learning framework can be used to classify other existing drugs for SARS‐CoV‐2 infection and can easily be adapted to drug effect prediction for future viral outbreaks.

Abstract Proteomics technologies, which include a diverse range of approaches such as mass spectrometry-based, array-based, and others, are key technologies for the identification of biomarkers and disease mechanisms, referred to as mechanotyping. Despite over 15,000 published studies in 2022 alone, leveraging publicly available proteomics data for biomarker identification, mechanotyping and drug target identification is not readily possible. Proteomic data addressing similar biological/biomedical questions are made available by multiple research groups in different locations using different model organisms. Furthermore, not only various organisms are employed but different assay systems, such as in vitro and in vivo systems, are used. Finally, even though proteomics data are deposited in public databases, such as ProteomeXchange, they are provided at different levels of detail. Thus, data integration is hampered by non-harmonized usage of identifiers when reviewing the literature or performing meta-analyses to consolidate existing publications into a joint picture. To address this problem, we present ProHarMeD, a tool for harmonizing and comparing proteomics data gathered in multiple studies and for the extraction of disease mechanisms and putative drug repurposing candidates. It is available as a website, Python library and R package. ProHarMeD facilitates ID and name conversions between protein and gene levels, or organisms via ortholog mapping, and provides detailed logs on the loss and gain of IDs after each step. The web tool further determines IDs shared by different studies, proposes potential disease mechanisms as well as drug repurposing candidates automatically, and visualizes these results interactively. We apply ProHarMeD to a set of four studies on bone regeneration. First, we demonstrate the benefit of ID harmonization which increases the number of shared genes between studies by 50%. Second, we identify a potential disease mechanism, with five corresponding drug targets, and the top 20 putative drug repurposing candidates, of which Fondaparinux, the candidate with the highest score, and multiple others are known to have an impact on bone regeneration. Hence, ProHarMeD allows users to harmonize multi-centric proteomics research data in meta-analyses, evaluates the success of the ID conversions and remappings, and finally, it closes the gaps between proteomics, disease mechanism mining and drug repurposing. It is publicly available at https://apps.cosy.bio/proharmed/ .

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.

Large-scale disease-association data are widely used for pathomechanism mining, even if disease definitions used for annotation are mostly phenotype-based. Here, the authors show that this bias can lead to a blurred view on disease mechanisms, highlighting the need for close-up studies based on molecular data for well-characterized patient cohorts.

Molecular interaction networks lay the foundation for studying how biological functions are controlled by the complex interplay of genes and proteins. Investigating perturbed processes using biological networks has been instrumental in uncovering mechanisms that underlie complex disease phenotypes. Rapid advances in omics technologies have prompted the generation of high-throughput datasets, enabling large-scale, network-based analyses. Consequently, various modeling techniques, including network enrichment, differential network extraction, and network inference, have proven to be useful for gaining new mechanistic insights. We provide an overview of recent network-based methods and their core ideas to facilitate the discovery of disease modules or candidate mechanisms. Knowledge generated from these computational efforts will benefit biomedical research, especially drug development and precision medicine. We further discuss current challenges and provide perspectives in the field, highlighting the need for more integrative and dynamic network approaches to model disease development and progression.

BackgroundMachine learning and artificial intelligence have shown promising results in many areas and are driven by the increasing amount of available data. However, these data are often distributed across different institutions and cannot be easily shared owing to strict privacy regulations. Federated learning (FL) allows the training of distributed machine learning models without sharing sensitive data. In addition, the implementation is time-consuming and requires advanced programming skills and complex technical infrastructures. ObjectiveVarious tools and frameworks have been developed to simplify the development of FL algorithms and provide the necessary technical infrastructure. Although there are many high-quality frameworks, most focus only on a single application case or method. To our knowledge, there are no generic frameworks, meaning that the existing solutions are restricted to a particular type of algorithm or application field. Furthermore, most of these frameworks provide an application programming interface that needs programming knowledge. There is no collection of ready-to-use FL algorithms that are extendable and allow users (eg, researchers) without programming knowledge to apply FL. A central FL platform for both FL algorithm developers and users does not exist. This study aimed to address this gap and make FL available to everyone by developing FeatureCloud, an all-in-one platform for FL in biomedicine and beyond. MethodsThe FeatureCloud platform consists of 3 main components: a global frontend, a global backend, and a local controller. Our platform uses a Docker to separate the local acting components of the platform from the sensitive data systems. We evaluated our platform using 4 different algorithms on 5 data sets for both accuracy and runtime. ResultsFeatureCloud removes the complexity of distributed systems for developers and end users by providing a comprehensive platform for executing multi-institutional FL analyses and implementing FL algorithms. Through its integrated artificial intelligence store, federated algorithms can easily be published and reused by the community. To secure sensitive raw data, FeatureCloud supports privacy-enhancing technologies to secure the shared local models and assures high standards in data privacy to comply with the strict General Data Protection Regulation. Our evaluation shows that applications developed in FeatureCloud can produce highly similar results compared with centralized approaches and scale well for an increasing number of participating sites. ConclusionsFeatureCloud provides a ready-to-use platform that integrates the development and execution of FL algorithms while reducing the complexity to a minimum and removing the hurdles of federated infrastructure. Thus, we believe that it has the potential to greatly increase the accessibility of privacy-preserving and distributed data analyses in biomedicine and beyond.

Mass cytometry (CyTOF) is a new technology that allows the investigation of protein expression at single cell level with high resolution. While several protocols are available to investigate leukocyte expression, platelet staining and analysis with CyTOF have been described only from whole blood. Moreover, available protocols do not allow sample storage but require fresh samples to be obtained, processed, and measured immediately. We provide a structured and reproducible method to stain platelets from platelet-rich plasma to study thrombocyte protein expression and reactivity using mass cytometry. With our method, it is possible to acquire a large number of events allowing deep bioinformatic investigation of platelet expression heterogeneity. Integrated in our protocol is also a previously established freezing protocol that allows the storage of stained samples and to delay their measurement. Finally, we provide a structured workflow using different platelet stimulators and a freely available bioinformatic pipeline to analyze platelet expression. Our protocol unlocks the potential of CyTOF analysis for studying platelet biology in health and disease.

Abstract Differential gene expression normalised to a single housekeeping (HK) is used to identify disease mechanisms and therapeutic targets. HK gene selection is often arbitrary, potentially introducing systematic error and discordant results. Here we examine these risks in a disease model of brain hypoxia. We first identified the eight most frequently used HK genes through a systematic review. However, we observe that in both ex-vivo and in vivo, their expression levels varied considerably between conditions. When applying these genes to normalise expression levels of the validated stroke target gene, inducible Nox4, we obtained opposing results. As an alternative tool for unbiased HK gene selection, software tools exist but are limited to individual datasets lacking genome-wide search capability and user-friendly interfaces. We, therefore, developed the HouseKeepR algorithm to rapidly analyse multiple gene expression datasets in a disease-specific manner and rank HK gene candidates according to stability in an unbiased manner. Using a panel of de novo top-ranked HK genes for brain hypoxia, but not single genes, Nox4 induction was consistently reproduced. Thus, differential gene expression analysis is best normalised against a HK gene panel selected in an unbiased manner. HouseKeepR is the first user-friendly, bias-free, and broadly applicable tool to automatically propose suitable HK genes in a tissue- and disease-dependent manner.

BackgroundThe collection, storage, and analysis of large data sets are relevant in many sectors. Especially in the medical field, the processing of patient data promises great progress in personalized health care. However, it is strictly regulated, such as by the General Data Protection Regulation (GDPR). These regulations mandate strict data security and data protection and, thus, create major challenges for collecting and using large data sets. Technologies such as federated learning (FL), especially paired with differential privacy (DP) and secure multiparty computation (SMPC), aim to solve these challenges. ObjectiveThis scoping review aimed to summarize the current discussion on the legal questions and concerns related to FL systems in medical research. We were particularly interested in whether and to what extent FL applications and training processes are compliant with the GDPR data protection law and whether the use of the aforementioned privacy-enhancing technologies (DP and SMPC) affects this legal compliance. We placed special emphasis on the consequences for medical research and development. MethodsWe performed a scoping review according to the PRISMA-ScR (Preferred Reporting Items for Systematic Reviews and Meta-Analyses extension for Scoping Reviews). We reviewed articles on Beck-Online, SSRN, ScienceDirect, arXiv, and Google Scholar published in German or English between 2016 and 2022. We examined 4 questions: whether local and global models are “personal data” as per the GDPR; what the “roles” as defined by the GDPR of various parties in FL are; who controls the data at various stages of the training process; and how, if at all, the use of privacy-enhancing technologies affects these findings. ResultsWe identified and summarized the findings of 56 relevant publications on FL. Local and likely also global models constitute personal data according to the GDPR. FL strengthens data protection but is still vulnerable to a number of attacks and the possibility of data leakage. These concerns can be successfully addressed through the privacy-enhancing technologies SMPC and DP. ConclusionsCombining FL with SMPC and DP is necessary to fulfill the legal data protection requirements (GDPR) in medical research dealing with personal data. Even though some technical and legal challenges remain, for example, the possibility of successful attacks on the system, combining FL with SMPC and DP creates enough security to satisfy the legal requirements of the GDPR. This combination thereby provides an attractive technical solution for health institutions willing to collaborate without exposing their data to risk. From a legal perspective, the combination provides enough built-in security measures to satisfy data protection requirements, and from a technical perspective, the combination provides secure systems with comparable performance with centralized machine learning applications.

Background: Blood–barrier (BBB) breakdown and active inflammation are hallmarks of relapsing multiple sclerosis (RMS), but the molecular events contributing to the development of new lesions are not well explored. Leaky endothelial junctions are associated with increased production of endothelial-derived extracellular microvesicles (EVs) and result in the entry of circulating immune cells into the brain. MRI with intravenous gadolinium (Gd) can visualize acute blood–barrier disruption as the initial event of the evolution of new lesions. Methods: Here, weekly MRI with Gd was combined with proteomics, multiplex immunoassay, and endothelial stress-optimized EV array to identify early markers related to BBB disruption. Five patients with RMS with no disease-modifying treatment were monitored weekly using high-resolution 3T MRI scanning with intravenous gadolinium (Gd) for 8 weeks. Patients were then divided into three groups (low, medium, or high MRI activity) defined by the number of new, total, and maximally enhancing Gd-enhancing lesions and the number of new FLAIR lesions. Plasma samples taken at each MRI were analyzed for protein biomarkers of inflammation by quantitative proteomics, and cytokines using multiplex immunoassays. EVs were characterized with an optimized endothelial stress EV array based on exosome surface protein markers for the detection of soluble secreted EVs. Results: Proteomics analysis of plasma yielded quantitative information on 208 proteins at each patient time point (n = 40). We observed the highest number of unique dysregulated proteins (DEPs) and the highest functional enrichment in the low vs. high MRI activity comparison. Complement activation and complement/coagulation cascade were also strongly overrepresented in the low vs. high MRI activity comparison. Activation of the alternative complement pathway, pathways of blood coagulation, extracellular matrix organization, and the regulation of TLR and IGF transport were unique for the low vs. high MRI activity comparison as well, with these pathways being overrepresented in the patient with high MRI activity. Principal component analysis indicated the individuality of plasma profiles in patients. IL-17 was upregulated at all time points during 8 weeks in patients with high vs. low MRI activity. Hierarchical clustering of soluble markers in the plasma indicated that all four MRI outcomes clustered together with IL-17, IL-12p70, and IL-1β. MRI outcomes also showed clustering with EV markers CD62E/P, MIC A/B, ICAM-1, and CD42A. The combined cluster of these cytokines, EV markers, and MRI outcomes clustered also with IL-12p40 and IL-7. All four MRI outcomes correlated positively with levels of IL-17 (p < 0.001, respectively), and EV-ICAM-1 (p < 0.0003, respectively). IL-1β levels positively correlated with the number of new Gd-enhancing lesions (p < 0.01), new FLAIR lesions (p < 0.001), and total number of Gd-enhancing lesions (p < 0.05). IL-6 levels positively correlated with the number of new FLAIR lesions (p < 0.05). Random Forests and linear mixed models identified IL-17, CCL17/TARC, CCL3/MIP-1α, and TNF-α as composite biomarkers predicting new lesion evolution. Conclusions: Combination of serial frequent MRI with proteome, neuroinflammation markers, and protein array data of EVs enabled assessment of temporal changes in inflammation and endothelial dysfunction in RMS related to the evolution of new and enhancing lesions. Particularly, the Th17 pathway and IL-1β clustered and correlated with new lesions and Gd enhancement, indicating their importance in BBB disruption and initiating acute brain inflammation in MS. In addition to the Th17 pathway, abundant protein changes between MRI activity groups suggested the role of EVs and the coagulation system along with innate immune responses including acute phase proteins, complement components, and neutrophil degranulation.

Abstract High-grade serous ovarian carcinoma (HGSC) is the most lethal gynecologic malignancy due to the lack of reliable biomarkers, effective treatment, and chemoresistance. Improving the diagnosis and the development of targeted therapies is still needed. The molecular pathomechanisms driving HGSC progression are not fully understood though crucial for effective diagnosis and identification of novel targeted therapy options. The oncogene CTCFL (BORIS), the paralog of CTCF, is a transcriptional factor highly expressed in ovarian cancer (but in rarely any other tissue in females) with cancer-specific characteristics and therapeutic potential. In this work, we seek to understand the regulatory functions of CTCFL to unravel new target genes with clinical relevance. We used in vitro models to evaluate the transcriptional changes due to the presence of CTCFL, followed by a selection of gene candidates using de novo network enrichment analysis. The resulting mechanistic candidates were further assessed regarding their prognostic potential and druggability. We show that CTCFL-driven genes are involved in cytoplasmic membrane functions; in particular, the PI3K-Akt initiators EGFR1 and VEGFA, as well as ITGB3 and ITGB6 are potential drug targets. Finally, we identified the CTCFL targets ACTBL2, MALT1 and PCDH7 as mechanistic biomarkers to predict survival in HGSC. Finally, we elucidated the value of CTCFL in combination with its targets as a prognostic marker profile for HGSC progression and as putative drug targets.

Abstract Meta-analysis has been established as an effective approach to combining summary statistics of several genome-wide association studies (GWAS). However, the accuracy of meta-analysis can be attenuated in the presence of cross-study heterogeneity. We present sPLINK, a hybrid federated and user-friendly tool, which performs privacy-aware GWAS on distributed datasets while preserving the accuracy of the results. sPLINK is robust against heterogeneous distributions of data across cohorts while meta-analysis considerably loses accuracy in such scenarios. sPLINK achieves practical runtime and acceptable network usage for chi-square and linear/logistic regression tests. sPLINK is available at https://exbio.wzw.tum.de/splink .

Many disease-causing genetic variants converge on common biological functions and pathways. Precisely how to incorporate pathway knowledge in genetic association studies is not yet clear, however. Previous approaches employ a two-step approach, in which a regular association test is first performed to identify variants associated with the disease phenotype, followed by a test for functional enrichment within the genes implicated by those variants. Here we introduce a concise one-step approach, Hierarchical Genetic Analysis (Higana), which directly computes phenotype associations against each function in the large hierarchy of biological functions documented by the Gene Ontology. Using this approach, we identify risk genes and functions for Chronic Obstructive Pulmonary Disease (COPD), highlighting microtubule transport, muscle adaptation, and nicotine receptor signaling pathways. Microtubule transport has not been previously linked to COPD, as it integrates genetic variants spread over numerous genes. All associations validate strongly in a second COPD cohort.

Abstract Alternative splicing (AS) is an important aspect of gene regulation. Nevertheless, its role in molecular processes and pathobiology is far from understood. A roadblock is that tools for the functional analysis of AS-set events are lacking. To mitigate this, we developed NEASE, a tool integrating pathways with structural annotations of protein-protein interactions to functionally characterize AS events. We show in four application cases how NEASE can identify pathways contributing to tissue identity and cell type development, and how it highlights splicing-related biomarkers. With a unique view on AS, NEASE generates unique and meaningful biological insights complementary to classical pathways analysis.

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

Infectious agents have been long considered to play a role in the pathogenesis of neurological diseases as part of the interaction between genetic susceptibility and the environment. The role of bacteria in CNS autoimmunity has also been highlighted by changes in the diversity of gut microbiota in patients with neurological diseases such as Parkinson’s disease, Alzheimer disease and multiple sclerosis, emphasizing the role of the gut-brain axis. We discuss the hypothesis of a brain microbiota, the BrainBiota: bacteria living in symbiosis with brain cells. Existence of various bacteria in the human brain is suggested by morphological evidence, presence of bacterial proteins, metabolites, transcripts and mucosal-associated invariant T cells. Based on our data, we discuss the hypothesis that these bacteria are an integral part of brain development and immune tolerance as well as directly linked to the gut microbiome. We further suggest that changes of the BrainBiota during brain diseases may be the consequence or cause of the chronic inflammation similarly to the gut microbiota.

There is an unmet need for adaptable tools allowing biomedical researchers to employ network-based drug repurposing approaches for their individual use cases. Here, the authors close this gap with NeDRex, an integrative and interactive platform.

Clinical time-to-event studies are dependent on large sample sizes, often not available at a single institution. However, this is countered by the fact that, particularly in the medical field, individual institutions are often legally unable to share their data, as medical data is subject to strong privacy protection due to its particular sensitivity. But the collection, and especially aggregation into centralized datasets, is also fraught with substantial legal risks and often outright unlawful. Existing solutions using federated learning have already demonstrated considerable potential as an alternative for central data collection. Unfortunately, current approaches are incomplete or not easily applicable in clinical studies owing to the complexity of federated infrastructures. This work presents privacy-aware and federated implementations of the most used time-to-event algorithms (survival curve, cumulative hazard rate, log-rank test, and Cox proportional hazards model) in clinical trials, based on a hybrid approach of federated learning, additive secret sharing, and differential privacy. On several benchmark datasets, we show that all algorithms produce highly similar, or in some cases, even identical results compared to traditional centralized time-to-event algorithms. Furthermore, we were able to reproduce the results of a previous clinical time-to-event study in various federated scenarios. All algorithms are accessible through the intuitive web-app Partea (https://partea.zbh.uni-hamburg.de), offering a graphical user interface for clinicians and non-computational researchers without programming knowledge. Partea removes the high infrastructural hurdles derived from existing federated learning approaches and removes the complexity of execution. Therefore, it is an easy-to-use alternative to central data collection, reducing bureaucratic efforts but also the legal risks associated with the processing of personal data to a minimum.

Glomerular disease due to podocyte malfunction is a major factor in the pathogenesis of chronic kidney disease. Identification of podocyte-specific signaling pathways is therefore a prerequisite to characterizing relevant disease pathways and developing novel treatment approaches. Here, we employed loss of function studies for EPB41L5 (Yurt) as a central podocyte gene to generate a cell type-specific disease model. Loss of Yurt in fly nephrocytes caused protein uptake and slit diaphragm defects. Transcriptomic and proteomic analysis of human EPB41L5 knockout podocytes demonstrated impaired mechanotransduction via the YAP/TAZ signaling pathway. Further analysis of specific inhibition of the YAP/TAZ-TEAD transcription factor complex by TEADi led to the identification of ARGHAP29 as an EPB41L5 and YAP/TAZ-dependently expressed podocyte RhoGAP. Knockdown of ARHGAP29 caused increased RhoA activation, defective lamellipodia formation, and increased maturation of integrin adhesion complexes, explaining similar phenotypes caused by loss of EPB41L5 and TEADi expression in podocytes. Detection of increased levels of ARHGAP29 in early disease stages of human glomerular disease implies a novel negative feedback loop for mechanotransductive RhoA—YAP/TAZ signaling in podocyte physiology and disease.

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

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 To identify markers in the CSF of multiple sclerosis (MS) subtypes, we used a two-step proteomic approach: (i) Discovery proteomics compared 169 pooled CSF from MS subtypes and inflammatory/degenerative CNS diseases (NMO spectrum and Alzheimer disease) and healthy controls. (ii) Next, 299 proteins selected by comprehensive statistics were quantified in 170 individual CSF samples. (iii) Genes of the identified proteins were also screened among transcripts in 73 MS brain lesions compared to 25 control brains. F-test based feature selection resulted in 8 proteins differentiating the MS subtypes, and secondary progressive (SP)MS was the most different also from controls. Genes of 7 out these 8 proteins were present in MS brain lesions: GOLM was significantly differentially expressed in active, chronic active, inactive and remyelinating lesions, FRZB in active and chronic active lesions, and SELENBP1 in inactive lesions. Volcano maps of normalized proteins in the different disease groups also indicated the highest amount of altered proteins in SPMS. Apolipoprotein C-I, apolipoprotein A-II, augurin, receptor-type tyrosine-protein phosphatase gamma, and trypsin-1 were upregulated in the CSF of MS subtypes compared to controls. This CSF profile and associated brain lesion spectrum highlight non-inflammatory mechanisms in differentiating CNS diseases and MS subtypes and the uniqueness of SPMS.

Abstract Novel coronavirus disease 2019 (COVID-19) is associated with a hypercoagulable state, characterized by abnormal coagulation parameters and by increased incidence of cardiovascular complications. With this study, we aimed to investigate the activation state and the expression of transmembrane proteins in platelets of hospitalized COVID-19 patients. We investigated transmembrane proteins expression with a customized mass cytometry panel of 21 antibodies. Platelets of 8 hospitalized COVID-19 patients not requiring intensive care support and without pre-existing conditions were compared to platelets of healthy controls (11 donors) with and without in vitro stimulation with thrombin receptor-activating peptide (TRAP). Mass cytometry of non-stimulated platelets detected an increased surface expression of activation markers P-Selectin (0.67 vs. 1.87 median signal intensity for controls vs. patients, p = 0.0015) and LAMP-3 (CD63, 0.37 vs. 0.81, p = 0.0004), the GPIIb/IIIa complex (4.58 vs. 5.03, p

Microorganisms including bacteria, fungi, viruses, protists and archaea live as communities in complex and contiguous environments. They engage in numerous inter- and intra- kingdom interactions which can be inferred from microbiome profiling data. In particular, network-based approaches have proven helpful in deciphering complex microbial interaction patterns. Here we give an overview of state-of-the-art methods to infer intra-kingdom interactions ranging from simple correlation- to complex conditional dependence-based methods. We highlight common biases encountered in microbial profiles and discuss mitigation strategies employed by different tools and their trade-off with increased computational complexity. Finally, we discuss current limitations that motivate further method development to infer inter-kingdom interactions and to robustly and comprehensively characterize microbial environments in the future.

Multiple sclerosis (MS) is an inflammatory demyelinating and degenerative disease of the central nervous system (CNS). Although inflammatory responses are efficiently treated, therapies for progression are scarce and suboptimal, and biomarkers to predict the disease course are insufficient. Cure or preventive measures for MS require knowledge of core pathological events at the site of the tissue damage. Novelties in systems biology have emerged and paved the way for a more fine-grained understanding of key pathological pathways within the CNS, but they have also raised questions still without answers. Here, we systemically review the power of tissue and single-cell/nucleus CNS omics and discuss major gaps of integration into the clinical practice. Systemic search identified 49 transcriptome and 11 proteome studies of the CNS from 1997 till October 2021. Pioneering molecular discoveries indicate that MS affects the whole brain and all resident cell types. Despite inconsistency of results, studies imply increase in transcripts/proteins of semaphorins, heat shock proteins, myelin proteins, apolipoproteins and HLAs. Different lesions are characterized by distinct astrocytic and microglial polarization, altered oligodendrogenesis, and changes in specific neuronal subtypes. In all white matter lesion types, CXCL12, SCD, CD163 are highly expressed, and STAT6- and TGFβ-signaling are increased. In the grey matter lesions, TNF-signaling seems to drive cell death, and especially CUX2-expressing neurons may be susceptible to neurodegeneration. The vast heterogeneity at both cellular and lesional levels may underlie the clinical heterogeneity of MS, and it may be more complex than the current disease phenotyping in the clinical practice. Systems biology has not solved the mystery of MS, but it has discovered multiple molecules and networks potentially contributing to the pathogenesis. However, these results are mostly descriptive; focused functional studies of the molecular changes may open up for a better interpretation. Guidelines for acceptable quality or awareness of results from low quality data, and standardized computational and biological pipelines may help to overcome limited tissue availability and the “snap shot” problem of omics. These may help in identifying core pathological events and point in directions for focus in clinical prevention.

Information developed to understand the molecular mechanisms of SARS-CoV-2 infection for predicting drug repurposing candidates is time-consuming to integrate and explore. Here, the authors develop an interactive online platform for virus-host interactome exploration and drug (target) identification.

Abstract Background Large numbers of autosomal sites are found differentially methylated in the aging genome. Due to analytical difficulties in dealing with sex differences in X-chromosome content and X-inactivation (XCI) in females, this has not been explored for the X chromosome. Methods Using data from middle age to elderly individuals (age 55+ years) from two Danish cohorts of monozygotic twins and the Scottish Lothian Birth Cohort 1921, we conducted an X-chromosome-wide analysis of age-associated DNA methylation patterns with consideration of stably inferred XCI status. Results Through analysing and comparing sex-specific X-linked DNA methylation changes over age late in life, we identified 123, 293 and 55 CpG sites significant (FDR

De novo pathway enrichment is a systems biology approach in which OMICS data are projected onto a molecular interaction network to identify subnetworks representing condition-specific functional modules and molecular pathways. Compared to classical pathway enrichment analysis methods, de novo pathway enrichment is not limited to predefined lists of pathways from (curated) databases and thus particularly suited for discovering novel disease mechanisms. While several tools have been proposed for pathway enrichment, the integration of de novo pathway enrichment in end-to-end OMICS analysis workflows in the R programming language is currently limited to a single tool. To close this gap, we have implemented an R package KeyPathwayMineR (KPM-R). The package extends the features and usability of existing versions of KeyPathwayMiner by leveraging the power, flexibility and versatility of R and by providing various novel functionalities for performing data preparation, visualization, and comparison. In addition, thanks to its interoperability with a plethora of existing R packages in e.g., Bioconductor, CRAN, and GitHub, KPM-R allows carrying out the initial preparation of the datasets and to meaningfully interpret the extracted subnetworks. To demonstrate the package’s potential, KPM-R was applied to bulk RNA-Seq data of nasopharyngeal swabs from SARS-CoV-2 infected individuals, and on single cell RNA-Seq data of aging mice tissue from the Tabula Muris Senis atlas.

Abstract To identify pathogenetic markers and potential drivers of different lesion types in the white matter (WM) of patients with progressive multiple sclerosis (PMS), we sequenced RNA from 73 different WM areas. Compared to 25 WM controls, 6713 out of 18,609 genes were significantly differentially expressed in MS tissues (FDR

Cognitive aging is one of the major problems worldwide, especially as people get older. This study aimed to perform global gene expression profiling of cognitive function to identify associated genes and pathways and a novel transcriptional regulatory network analysis to identify important regulons. We performed single transcript analysis on 400 monozygotic twins using an assumption-free generalized correlation coefficient (GCC), linear mixed-effect model (LME) and kinship model and identified six probes (one significant at the standard FDR < 0.05 while the other results were suggestive with 0.18 ≤ FDR ≤ 0.28). We combined the GCC and linear model results to cover diverse patterns of relationships, and meaningful and novel genes like APOBEC3G, H6PD, SLC45A1, GRIN3B, and PDE4D were detected. Our exploratory study showed the downregulation of all these genes with increasing cognitive function or vice versa except the SLC45A1 gene, which was upregulated with increasing cognitive function. Linear models found only H6PD and SLC45A1, the other genes were captured by GCC. Significant functional pathways (FDR < 3.95e-10) such as focal adhesion, ribosome, cysteine and methionine metabolism, Huntington's disease, eukaryotic translation elongation, nervous system development, influenza infection, metabolism of RNA, and cell cycle were identified. A total of five regulons (FDR< 1.3e-4) were enriched in a transcriptional regulatory analysis in which CTCF and REST were activated and SP3, SRF, and XBP1 were repressed regulons. The genome-wide transcription analysis using both assumption-free GCC and linear models identified important genes and biological pathways implicated in cognitive performance, cognitive aging, and neurological diseases. Also, the regulatory network analysis revealed significant activated and repressed regulons on cognitive function.

Background: Artificial Intelligence (AI) in healthcare has demonstrated high efficiency in academic research, while only few, and predominantly small, real-world AI applications exist in the preventive, diagnostic and therapeutic contexts. Our identification and analysis of success factors for the implementation of AI aims to close the gap between recent years' significant academic AI advancements and the comparably low level of practical application in healthcare.Methods: A literature and real life cases analysis was conducted in Scopus and OpacPlus as well as the Google advanced search database. The according search queries have been defined based on success factor categories for AI implementation derived from a prior World Health Organization survey about barriers of adoption of Big Data within 125 countries. The eligible publications and real life cases were identified through a catalog of in- and exclusion criteria focused on concrete AI application cases. These were then analyzed to deduct and discuss success factors that facilitate or inhibit a broad-scale implementation of AI in healthcare.Results: The analysis revealed three categories of success factors, namely (1) policy setting, (2) technological implementation, and (3) medical and economic impact measurement. For each of them a set of recommendations has been deducted: First, a risk adjusted policy frame is required that distinguishes between precautionary and permissionless principles, and differentiates among accountability, liability, and culpability. Second, a “privacy by design” centered technology infrastructure shall be applied that enables practical and legally compliant data access. Third, the medical and economic impact need to be quantified, e.g., through the measurement of quality-adjusted life years while applying the CHEERS and PRISMA reporting criteria.Conclusions: Private and public institutions can already today leverage AI implementation based on the identified results and thus drive the translation from scientific development to real world application. Additional success factors could include trust-building measures, data categorization guidelines, and risk level assessments and as the success factors are interlinked, future research should elaborate on their optimal interaction to utilize the full potential of AI in real world application.

Small RNAs (sRNAs) are one of the key players in the post-transcriptional regulation of bacterial gene expression. These molecules, together with transcription factors, form regulatory networks and greatly influence the bacterial regulatory landscape. Little is known concerning sRNAs and their influence on the regulatory machinery in the genus Corynebacterium, despite its medical, veterinary and biotechnological importance. Here, we expand corynebacterial regulatory knowledge by integrating sRNAs and their regulatory interactions into the transcriptional regulatory networks of six corynebacterial species, covering four human and animal pathogens, and integrate this data into the CoryneRegNet database. To this end, we predicted sRNAs to regulate 754 genes, including 206 transcription factors, in corynebacterial gene regulatory networks. Amongst them, the sRNA Cd-NCTC13129-sRNA-2 is predicted to directly regulate ydfH, which indirectly regulates 66 genes, including the global regulator glxR in C. diphtheriae. All of the sRNA-enriched regulatory networks of the genus Corynebacterium have been made publicly available in the newest release of CoryneRegNet(www.exbio.wzw.tum.de/coryneregnet/) to aid in providing valuable insights and to guide future experiments.

Abstract Background Ovarian carcinomas presenting homologous recombination deficiency (HRD), which is observed in about 50% of cases, are more sensitive to platinum and PARP inhibitor therapies. Although platinum resistant disease has a low chance to be responsive to platinum-based chemotherapy, a set of patients is retreated with platinum and some of them are responsive. In this study, we evaluated copy number alterations, HR gene mutations and HR deficiency scores in ovarian cancer patients with prolonged platinum sensitivity. Methods In this retrospective study (2005 to 2014), we selected 31 patients with platinum resistant ovarian cancer retreated with platinum therapy. Copy number alterations and HR scores were evaluated using the OncoScan® FFPE platform in 15 cases. The mutational profile of 24 genes was investigated by targeted-NGS. Results The median values of the four HRD scores were higher in responders (LOH = 15, LST = 28, tAI = 33, CS = 84) compared with non-responders (LOH = 7.5, LST = 17.5, tAI = 23, CS = 47). Patients with high LOH, LST, tAI and CS scores had better response rates, although these differences were not statistically significant. Response rate to platinum retreatment was 22% in patients with CCNE1 gains and 83.5% in patients with no CCNE1 gains (p = 0.041). Furthermore, response rate was 54.5% in patients with RB1 loss and 25% in patients without RB1 loss (p = 0.569). Patients with CCNE1 gains showed a worse progression free survival (PFS = 11.1 months vs 3.7 months; p = 0.008) and a shorter overall survival (OS = 39.3 months vs 7.1 months; p = 0.007) in comparison with patients with no CCNE1 gains. Patients with RB1 loss had better PFS (9.0 months vs 2.6 months; p = 0.093) and OS (27.4 months vs 3.6 months; p = 0.025) compared with cases with no RB1 loss. Four tumor samples were BRCA mutated and tumor mutations were not associated with response to treatment. Conclusions HR deficiency was found in 60% of our cases and HRD medium values were higher in responders than in non-responders. Despite the small number of patients tested, CCNE1 gain and RB1 loss discriminate patients with tumors extremely sensitive to platinum retreatment.

Abstract 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 Background A revolutionary diversion from classical vaccinology to reverse vaccinology approach has been observed in the last decade. The ever-increasing genomic and proteomic data has greatly facilitated the vaccine designing and development process. Reverse vaccinology is considered as a cost-effective and proficient approach to screen the entire pathogen genome. To look for broad-spectrum immunogenic targets and analysis of closely-related bacterial species, the assimilation of pangenome concept into reverse vaccinology approach is essential. The categories of species pangenome such as core, accessory, and unique genes sets can be analyzed for the identification of vaccine candidates through reverse vaccinology. Results We have designed an integrative computational pipeline term as “PanRV” that employs both the pangenome and reverse vaccinology approaches. PanRV comprises of four functional modules including i) Pangenome Estimation Module (PGM) ii) Reverse Vaccinology Module (RVM) iii) Functional Annotation Module (FAM) and iv) Antibiotic Resistance Association Module (ARM). The pipeline is tested by using genomic data from 301 genomes of Staphylococcus aureus and the results are verified by experimentally known antigenic data. Conclusion The proposed pipeline has proved to be the first comprehensive automated pipeline that can precisely identify putative vaccine candidates exploiting the microbial pangenome. PanRV is a Linux based package developed in JAVA language. An executable installer is provided for ease of installation along with a user manual at https://sourceforge.net/projects/panrv2/.

Abstract Background Multiple epigenome-wide association studies have been performed to identify DNA methylation patterns regulated by aging or correlated with risk of death. However, the inter-relatedness of the epigenetic basis of aging and mortality has not been well investigated. Methods Using genome-wide DNA methylation data from the Lothian Birth Cohorts, we conducted a genome-wide association analysis of all-cause mortality and compared this with age-associated methylation patterns reported on the same samples. Results Survival analysis using the Cox regression model identified 2552 CpG sites with genome-wide significance (false discovery rate

In 16S rRNA gene sequencing, certain bacterial genera were found to be underrepresented or even missing in taxonomic profiles when using unsuitable primer combinations, outdated reference databases, or inadequate pipeline settings. Concerning the last, quality thresholds as well as bioinformatic settings (i.e., clustering approach, analysis pipeline, and specific adjustments such as truncation) are responsible for a number of observed differences between studies.

Patient classification based on clinical and genomic data will further the goal of precision medicine. Interpretability is of particular relevance for models based on genomic data, where sample sizes are relatively small (in the hundreds), increasing overfitting risk netDx is a machine learning method to integrate multi-modal patient data and build a patient classifier. Patient data are converted into networks of patient similarity, which is intuitive to clinicians who also use patient similarity for medical diagnosis. Features passing selection are integrated, and new patients are assigned to the class with the greatest profile similarity. netDx has excellent performance, outperforming most machine-learning methods in binary cancer survival prediction. It handles missing data – a common problem in real-world data – without requiring imputation. netDx also has excellent interpretability, with native support to group genes into pathways for mechanistic insight into predictive features. The netDx Bioconductor package provides multiple workflows for users to build custom patient classifiers. It provides turnkey functions for one-step predictor generation from multi-modal data, including feature selection over multiple train/test data splits. Workflows offer versatility with custom feature design, choice of similarity metric; speed is improved by parallel execution. Built-in functions and examples allow users to compute model performance metrics such as AUROC, AUPR, and accuracy. netDx uses RCy3 to visualize top-scoring pathways and the final integrated patient network in Cytoscape. Advanced users can build more complex predictor designs with functional building blocks used in the default design. Finally, the netDx Bioconductor package provides a novel workflow for pathway-based patient classification from sparse genetic data.

Hypertension is the most important cause of death and disability in the elderly. In 9 out of 10 cases, the molecular cause, however, is unknown. One mechanistic hypothesis involves impaired endothelium-dependent vasodilation through reactive oxygen species (ROS) formation. Indeed, ROS forming NADPH oxidase (Nox) genes associate with hypertension, yet target validation has been negative. We re-investigate this association by molecular network analysis and identify NOX5, not present in rodents, as a sole neighbor to human vasodilatory endothelial nitric oxide (NO) signaling. In hypertensive patients, endothelial microparticles indeed contained higher levels of NOX5-but not NOX1, NOX2, or NOX4-with a bimodal distribution correlating with disease severity. Mechanistically, mice expressing human Nox5 in endothelial cells developed-upon aging-severe systolic hypertension and impaired endothelium-dependent vasodilation due to uncoupled NO synthase (NOS). We conclude that NOX5-induced uncoupling of endothelial NOS is a causal mechanism and theragnostic target of an age-related hypertension endotype. Nox5 knock-in (KI) mice represent the first mechanism-based animal model of hypertension.

Summary: Manipulating molecules that impact T cell receptor (TCR) or cytokine signaling, such as the protein tyrosine phosphatase non-receptor type 2 (PTPN2), has significant potential for advancing T cell-based immunotherapies. Nonetheless, it remains unclear how PTPN2 impacts the activation, survival, and memory formation of T cells. We find that PTPN2 deficiency renders cells in vivo and in vitro less dependent on survival-promoting cytokines, such as interleukin (IL)-2 and IL-15. Remarkably, briefly ex vivo-activated PTPN2-deficient T cells accumulate in 3- to 11-fold higher numbers following transfer into unmanipulated, antigen-free mice. Moreover, the absence of PTPN2 augments the survival of short-lived effector T cells and allows them to robustly re-expand upon secondary challenge. Importantly, we find no evidence for impaired effector function or memory formation. Mechanistically, PTPN2 deficiency causes broad changes in the expression and phosphorylation of T cell expansion and survival-associated proteins. Altogether, our data underline the therapeutic potential of targeting PTPN2 in T cell-based therapies to augment the number and survival capacity of antigen-specific T cells.

Abstract Background Copy number variations (CNVs) are large segments of the genome that are duplicated or deleted. Structural variations in the genome have been linked to many complex diseases. Similar to how genome-wide association studies (GWAS) have helped discover single-nucleotide polymorphisms linked to disease phenotypes, the extension of GWAS to CNVs has aided the discovery of structural variants associated with human traits and diseases. Results We present CoNVaQ, an easy-to-use web-based tool for CNV-based association studies. The web service allows users to upload two sets of CNV segments and search for genomic regions where the occurrence of CNVs is significantly associated with the phenotype. CoNVaQ provides two models: a simple statistical model using Fisher’s exact test and a novel query-based model matching regions to user-defined queries. For each region, the method computes a global q-value statistic by repeated permutation of samples among the populations. We demonstrate our platform by using it to analyze a data set of HPV-positive and HPV-negative penile cancer patients. Conclusions CoNVaQ provides a simple workflow for performing CNV-based association studies. It is made available as a web platform in order to provide a user-friendly workflow for biologists and clinicians to carry out CNV data analysis without installing any software. Through the web interface, users are also able to analyze their results to find overrepresented GO terms and pathways. In addition, our method is also available as a package for the R programming language. CoNVaQ is available at https://convaq.compbio.sdu.dk.

As the experimental discovery of microRNAs (miRNAs) is cumbersome, computational tools have been developed for the prediction of pre-miRNAs. Here the authors develop a framework to assess the performance of existing and novel pre-miRNA prediction tools and provide guidelines for selecting an appropriate approach for a given data set.

Computational biology: Evaluation of network-based pathway enrichment tools De novo pathway enrichment methods are essential to understand disease complexity. They can uncover disease-specific functional modules by integrating molecular interaction networks with expression profiles. However, how should researchers choose one method out of several? In this article, a group of scientists from Denmark and Germany presents the first attempt to quantitatively evaluate existing methods. This framework will help the biomedical community to find the appropriate tool(s) for their data. They created synthetic gold standards and simulated expression profiles to perform a systematic assessment of various tools. They observed that the choice of interaction network, parameter settings, preprocessing of expression data and statistical properties of the expression profiles influence the results to a large extent. The results reveal strengths and limitations of the individual methods and suggest using two or more tools to obtain comprehensive disease-modules.

Pre-operative 5-fluoracil-based chemoradiotherapy (nCRT) is the standard treatment for patients with locally advanced rectal cancer (LARC). Patients with pathological complete response (pCR–0% of tumor cells in the surgical specimen after nCRT) have better overall survival and lower risk of recurrence in comparison with incomplete responders (pIR). Predictive biomarkers to be used for new therapeutic strategies and capable of stratifying patients to avoid overtreatment are needed. We evaluated the genomic profiles of 33 pre-treatment LARC biopsies using SNP array and targeted-next generation sequencing (tNGS). Based on the large number of identified genomic alterations, we calculated the genomic instability index (GII) and three homologous recombination deficiency (HRD) scores, which have been reported as impaired DNA repair markers. We observed high GII in our LARC cases, which was confirmed in 165 rectal cancer cases from TCGA. Patients with pCR presented higher GII compared with pIR. Moreover, a negative correlation between GII and the fraction of tumor cells remaining after surgery was observed (ρ = –0.382, P = 0.02). High HRD scores were detected in 61% of LARC, of which 70% were incomplete responders. Using tNGS (105 cancer-related genes, 13 involved in HR and 5 in mismatch repair pathways), we identified 23% of cases with mutations in HR genes, mostly in pIR cases (86% of mutated cases). In agreement, the analysis of the TCGA dataset (N = 145) revealed 21% of tumors with mutations in HR genes. The HRD scores were shown to be predictive of better response to PARP-inhibitors and platinum-based chemotherapy in breast and ovarian cancer. Our results suggest that the same strategy could be applied in a set of LARC patients with HRD. In conclusion, we identified high genomic instability in LARC, which was related to alterations in the HR pathway, especially in pIR. These findings suggest that patients with impaired HRD would clinically benefit from PARP-inhibitors and platinum-based therapy.