Volume 402, Issue 8

Using results from genome-wide association studies for understanding complex traits is a current challenge. Here we review how genotype data can be used with different machine learning (ML) methods to predict phenotype occurrence and severity from genotype data. We discuss common feature encoding schemes and how studies handle the often small number of samples compared to the huge number of variants. We compare which ML methods are being applied, including recent results using deep neural networks. Further, we review the application of methods for feature explanation and interpretation.

The placenta is a highly vascularized and complex foetal organ that performs various tasks, crucial to a healthy pregnancy. Its dysfunction leads to complications such as stillbirth, preeclampsia, and intrauterine growth restriction. The specific cause of placental dysfunction remains unknown. Recently, the role of mitochondrial function and mitochondrial adaptations in the context of angiogenesis and placental dysfunction is getting more attention. The required energy for placental remodelling, nutrient transport, hormone synthesis, and the reactive oxygen species leads to oxidative stress, stemming from mitochondria. Mitochondria adapt to environmental changes and have been shown to adjust their oxygen and nutrient use to best support placental angiogenesis and foetal development. Angiogenesis is the process by which blood vessels form and is essential for the delivery of nutrients to the body. This process is regulated by different factors, pro-angiogenic factors and anti-angiogenic factors, such as sFlt-1. Increased circulating sFlt-1 levels have been linked to different preeclamptic phenotypes. One of many effects of increased sFlt-1 levels, is the dysregulation of mitochondrial function. This review covers mitochondrial adaptations during placentation, the importance of the anti-angiogenic factor sFlt-1in placental dysfunction and its role in the dysregulation of mitochondrial function.

Ischaemic heart disease is among the most frequent causes of death. Early detection of myocardial pathologies can increase the benefit of therapy and reduce the number of lethal cases. Presence of myocardial scar is an indicator for developing ischaemic heart disease and can be detected with high diagnostic precision by magnetic resonance imaging. However, magnetic resonance imaging scanners are expensive and of limited availability. It is known that presence of myocardial scar has an impact on the well-established, reasonably low cost, and almost ubiquitously available electrocardiogram. However, this impact is non-specific and often hard to detect by a physician. We present an artificial intelligence based approach — namely a deep learning model — for the prediction of myocardial scar based on an electrocardiogram and additional clinical parameters. The model was trained and evaluated by applying 6-fold cross-validation to a dataset of 12-lead electrocardiogram time series together with clinical parameters. The proposed model for predicting the presence of scar tissue achieved an area under the curve score, sensitivity, specificity, and accuracy of 0.89, 70.0, 84.3, and 78.0%, respectively. This promisingly high diagnostic precision of our electrocardiogram-based deep learning models for myocardial scar detection may support a novel, comprehensible screening method.

Reactive oxygen species are produced by a number of stimuli and can lead both to irreversible intracellular damage and signaling through reversible post-translational modification. It is unclear which factors contribute to the sensitivity of cysteines to redox modification. Here, we used statistical and machine learning methods to investigate the influence of different structural and sequence features on the modifiability of cysteines. We found several strong structural predictors for redox modification. Sensitive cysteines tend to be characterized by higher exposure, a lack of secondary structure elements, and a high number of positively charged amino acids in their close environment. Our results indicate that modified cysteines tend to occur close to other post-translational modifications, such as phosphorylated serines. We used these features to create models and predict the presence of redox-modifiable cysteines in human mitochondrial complex I as well as make novel predictions regarding redox-sensitive cysteines in proteins.

Matrix targeting sequences (MTSs) direct proteins from the cytosol into mitochondria. Efficient targeting often relies on internal matrix targeting-like sequences (iMTS-Ls) which share structural features with MTSs. Predicting iMTS-Ls was tedious and required multiple tools and webservices. We present iMLP, a deep learning approach for the prediction of iMTS-Ls in protein sequences. A recurrent neural network has been trained to predict iMTS-L propensity profiles for protein sequences of interest. The iMLP predictor considerably exceeds the speed of existing approaches. Expanding on our previous work on iMTS-L prediction, we now serve an intuitive iMLP webservice available at http://iMLP.bio.uni-kl.de and a stand-alone command line tool for power user in addition.

According to the amino acid composition of natural proteins, it could be expected that all possible sequences of three or four amino acids will occur at least once in large protein datasets purely by chance. However, in some species or cellular context, specific short amino acid motifs are missing due to unknown reasons. We describe these as Avoided Motifs, short amino acid combinations missing from biological sequences. Here we identify 209 human and 154 bacterial Avoided Motifs of length four amino acids, and discuss their possible functionality according to their presence in other species. Furthermore, we determine two Avoided Motifs of length three amino acids in human proteins specifically located in the cytoplasm, and two more in secreted proteins. Our results support the hypothesis that the characterization of Avoided Motifs in particular contexts can provide us with information about functional motifs, pointing to a new approach in the use of molecular sequences for the discovery of protein function.

In order to take advantage of the continuously increasing number of transcriptome studies, it is important to develop strategies that integrate multiple expression datasets addressing the same biological question to allow a robust analysis. Here, we propose a meta-analysis framework that integrates enriched pathways identified through the Gene Set Enrichment Analysis (GSEA) approach and calculates for each meta-pathway an empirical p -value. Validation of our approach on benchmark datasets showed comparable or even better performance than existing methods and an increase in robustness with increasing number of integrated datasets. We then applied the meta-analysis framework to 15 functional genomics datasets of physiological and pathological cardiac hypertrophy. Within these datasets we grouped expression sets measured at time points that represent the same hallmarks of heart tissue remodeling (‘aggregated time points’) and performed meta-analysis on the expression sets assigned to each aggregated time point. To facilitate biological interpretation, results were visualized as gene set enrichment networks. Here, our meta-analysis framework identified well-known biological mechanisms associated with pathological cardiac hypertrophy (e.g., cardiomyocyte apoptosis, cardiac contractile dysfunction, and alteration in energy metabolism). In addition, results highlighted novel, potentially cardioprotective mechanisms in physiological cardiac hypertrophy involving the down-regulation of immune cell response, which are worth further investigation.

Genome-wide CRISPR screens are becoming more widespread and allow the simultaneous interrogation of thousands of genomic regions. Although recent progress has been made in the analysis of CRISPR screens, it is still an open problem how to interpret CRISPR mutations in non-coding regions of the genome. Most of the tools concentrate on the interpretation of mutations introduced in gene coding regions. We introduce a computational pipeline that uses epigenomic information about regulatory elements for the interpretation of CRISPR mutations in non-coding regions. We illustrate our analysis protocol on the analysis of a genome-wide CRISPR screen in hTERT-RPE1 cells and reveal novel regulatory elements that mediate chemoresistance against doxorubicin in these cells. We infer links to established and to novel chemoresistance genes. Our analysis protocol is general and can be applied on any cell type and with different CRISPR enzymes.

A common approach to genome annotation involves the use of homology-based tools for the prediction of the functional role of proteins. The quality of functional annotations is dependent on the reference data used, as such, choosing the appropriate sources is crucial. Unfortunately, no single reference data source can be universally considered the gold standard, thus using multiple references could potentially increase annotation quality and coverage. However, this comes with challenges, particularly due to the introduction of redundant and exclusive annotations. Through text mining it is possible to identify highly similar functional descriptions, thus strengthening the confidence of the final protein functional annotation and providing a redundancy-free output. Here we present UniFunc, a text mining approach that is able to detect similar functional descriptions with high precision. UniFunc was built as a small module and can be independently used or integrated into protein function annotation pipelines. By removing the need to individually analyse and compare annotation results, UniFunc streamlines the complementary use of multiple reference datasets.

Human lymph nodes play a central part of immune defense against infection agents and tumor cells. Lymphoid follicles are compartments of the lymph node which are spherical, mainly filled with B cells. B cells are cellular components of the adaptive immune systems. In the course of a specific immune response, lymphoid follicles pass different morphological differentiation stages. The morphology and the spatial distribution of lymphoid follicles can be sometimes associated to a particular causative agent and development stage of a disease. We report our new approach for the automatic detection of follicular regions in histological whole slide images of tissue sections immuno-stained with actin. The method is divided in two phases: (1) shock filter-based detection of transition points and (2) segmentation of follicular regions. Follicular regions in 10 whole slide images were manually annotated by visual inspection, and sample surveys were conducted by an expert pathologist. The results of our method were validated by comparing with the manual annotation. On average, we could achieve a Zijbendos similarity index of 0.71, with a standard deviation of 0.07.

Biomedicines are complex biochemical formulations with multiple components that require extensive quality control during manufacturing and in subsequent batch testing. A proof-of-concept study has shown that an application of Raman spectroscopy can be beneficial for a classification of vaccines. However, the complexity of biomedicines introduces new challenges to spectroscopic methodology that require advanced experimental protocols. We further show the impact of analytical protocols on vaccine classification using R as an Open Source data analysis platform. In conclusion, we advocate for standardized and transparent experimental and analytical procedures and discuss current findings and open challenges.