Search

The joint analysis of two datasets X and Y that describe the same phenomena (e.g. the cellular state), but measure disjoint sets of variables (e.g. mRNA vs. protein levels) is currently challenging. Traditional methods typically analyze single interaction patterns such as variance or covariance. However, problem-tailored external knowledge may contain multiple different information about the interaction between the measured variables. We introduce MIASA, a holistic framework for the joint analysis of multiple different variables. It consists of assembling multiple different information such as similarity vs. association, expressed in terms of interaction-scores or distances, for subsequent clustering/classification. In addition, our framework includes a novel qualitative Euclidean embedding method (qEE-Transition) which enables using Euclidean-distance/vector-based clustering/classification methods on datasets that have a non-Euclidean-based interaction structure. As an alternative to conventional optimization-based multidimensional scaling methods which are prone to uncertainties, our qEE-Transition generates a new vector representation for each element of the dataset union X∪Y in a common Euclidean space while strictly preserving the original ordering of the assembled interaction-distances. To demonstrate our work, we applied the framework to three types of simulated datasets: samples from families of distributions, samples from correlated random variables, and time-courses of statistical moments for three different types of stochastic two-gene interaction models. We then compared different clustering methods with vs. without the qEE-Transition. For all examples, we found that the qEE-Transition followed by Ward clustering had superior performance compared to non-agglomerative clustering methods but had a varied performance against ultrametric-based agglomerative methods. We also tested the qEE-Transition followed by supervised and unsupervised machine learning methods and found promising results, however, more work is needed for optimal parametrization of these methods. As a future perspective, our framework points to the importance of more developments and validation of distance-distribution models aiming to capture multiple-complex interactions between different variables. [ABSTRACT FROM AUTHOR]

Objective: Due to the small size of the murine knee joint, extracting the chondrocyte transcriptome from articular cartilage (AC) is a major technical challenge. In this study, we demonstrate a new pragmatic approach of combining bulk RNA-sequencing (RNA-seq) and single cell (sc)RNA-seq to address this problem.Design: We propose a new cutting strategy for the murine femur which produces three segments with a predictable mixed cell population, where one segment contains AC and growth plate (GP) chondrocytes, another GP chondrocytes, and the last segment only bone and bone marrow. We analysed the bulk RNA-seq of the different segments to find distinct genes between the segments. The segment containing AC chondrocytes was digested and analysed via scRNA-seq.Results: Differential expression analysis using bulk RNA-seq identified 350 candidate chondrocyte gene in the AC segment. Gene set enrichment analysis of these genes revealed biological processes related- and non-related to chondrocytes, including, cartilage development (adj. P-value: 3.45E-17) and endochondral bone growth (adj. P-value 1.22E-4), respectively. ScRNA-seq of the AC segment found a cluster of 131 cells containing mainly chondrocytes. This cluster had 759 differentially expressed genes which enriched for extracellular matrix organisation (adj. P-value 7.76E-40) and other joint development processes. The intersection of the gene sets of bulk- and scRNA-seq contained 75 genes.Conclusions: Based on our results, we conclude that the combination of the two RNA-seq methods is necessary to precisely delineate the chondrocyte transcriptome and to study the disease phenotypes of chondrocytes in murine OA models in the future. [ABSTRACT FROM AUTHOR]

Molecular simulations of ligand–receptor interactions are a computational challenge, especially when their association- ('on'-rate) and dissociation- ('off'-rate) mechanisms are working on vastly differing timescales. One way of tackling this multiscale problem is to compute the free-energy landscapes, where molecular dynamics (MD) trajectories are used to only produce certain statistical ensembles. The approach allows for deriving the transition rates between energy states as a function of the height of the activation-energy barriers. In this article, we derive the association rates of the opioids fentanyl and N-(3-fluoro-1-phenethylpiperidin-4-yl)-N-phenyl propionamide (NFEPP) in a μ-opioid receptor by combining the free-energy landscape approach with the square-root-approximation method (SQRA), which is a particularly robust version of Markov modelling. The novelty of this work is that we derive the association rates as a function of the pH level using only an ensemble of MD simulations. We also verify our MD-derived insights by reproducing the in vitro study performed by the Stein Lab. [ABSTRACT FROM AUTHOR]

Human T-lymphotropic virus type-1 (HTLV-1) persists within hosts via infectious spread (de novo infection) and mitotic spread (infected cell proliferation), creating a population structure of multiple clones (infected cell populations with identical genomic proviral integration sites). The relative contributions of infectious and mitotic spread to HTLV-1 persistence are unknown, and will determine the efficacy of different approaches to treatment. The prevailing view is that infectious spread is negligible in HTLV-1 persistence beyond early infection. However, in light of recent high-throughput data on the abundance of HTLV-1 clones, and recent estimates of HTLV-1 clonal diversity that are substantially higher than previously thought (typically between 104 and 105 HTLV-1+ T cell clones in the body of an asymptomatic carrier or patient with HTLV-1-associated myelopathy/tropical spastic paraparesis), ongoing infectious spread during chronic infection remains possible. We estimate the ratio of infectious to mitotic spread using a hybrid model of deterministic and stochastic processes, fitted to previously published HTLV-1 clonal diversity estimates. We investigate the robustness of our estimates using three alternative estimators. We find that, contrary to previous belief, infectious spread persists during chronic infection, even after HTLV-1 proviral load has reached its set point, and we estimate that between 100 and 200 new HTLV-1 clones are created and killed every day. We find broad agreement between all estimators. The risk of HTLV-1-associated malignancy and inflammatory disease is strongly correlated with proviral load, which in turn is correlated with the number of HTLV-1-infected clones, which are created by de novo infection. Our results therefore imply that suppression of de novo infection may reduce the risk of malignant transformation. Author summary: There is no effective antiretroviral treatment for infection with Human T-lymphotropic virus type-1 (HTLV-1), which causes a range of inflammatory diseases and the aggressive malignancy Adult T-cell Leukaemia/Lymphoma (ATL) in approximately 10% of infected people. Within hosts the virus spreads via infectious spread (de novo infection) and mitotic spread (infected cell division). The relative contributions of each mechanism are unknown, and have major implications for drug development and clinical management of infection. We estimate the ratio of infectious to mitotic spread during the infection's chronic phase using three methods. Each method indicates infectious spread at low but persistent levels after proviral load has reached set point, contrary to the prevailing view that infectious spread features in early infection only. Risk of disease in HTLV-1 infection is known to increase with proviral load, via mutations accrued from repeated infected cell division. Our analyses suggest that ongoing infectious spread may provide an additional mechanism whereby chronic infection becomes malignant. Further, because antiretroviral drugs against Human Immunodeficiency Virus (HIV) inhibit HTLV-1 infectious spread, they may reduce the risk of HTLV-1 malignancy. [ABSTRACT FROM AUTHOR]

Muscle fibre cross-sectional area (CSA) is an important biomedical measure used to determine the structural composition of skeletal muscle, and it is relevant for tackling research questions in many different fields of research. To date, time consuming and tedious manual delineation of muscle fibres is often used to determine the CSA. Few methods are able to automatically detect muscle fibres in muscle fibre cross-sections to quantify CSA due to challenges posed by variation of brightness and noise in the staining images. In this paper, we introduce the supervised learningcomputer vision combined pipeline (SLCV), a robust semi-automatic pipeline for muscle fibre detection, which combines supervised learning (SL) with computer vision (CV). SLCV is adaptable to different staining methods and is quickly and intuitively tunable by the user. We are the first to perform an error analysis with respect to cell count and area, based on which we compare SLCV to the best purely CV-based pipeline in order to identify the contribution of SL and CV steps to muscle fibre detection. Our results obtained on 27 fluorescence-stained cross-sectional images of varying staining quality suggest that combining SL and CV performs significantly better than both SL-based and CV-basedmethods with regards to both the cell separation- and the area reconstruction error. Furthermore, applying SLCV to our test set images yielded fibre detection results of very high quality, with average sensitivity values of 0.93 or higher on different cluster sizes and an average Dice similarity coefficient of 0.9778. [ABSTRACT FROM AUTHOR]

Osteoarthritis (OA) is the most common cause of disability in ageing societies, with no effective therapies available to date. Two preclinical models are widely used to validate novel OA interventions (MCL-MM and DMM). Our aim is to discern disease dynamics in these models to provide a clear timeline in which various pathological changes occur. OA was surgically induced in mice by destabilisation of the medial meniscus. Analysis of OA progression revealed that the intensity and duration of chondrocyte loss and cartilage lesion formation were significantly different in MCL-MM vs DMM. Firstly, apoptosis was seen prior to week two and was narrowly restricted to the weight bearing area. Four weeks post injury the magnitude of apoptosis led to a 40–60% reduction of chondrocytes in the non-calcified zone. Secondly, the progression of cell loss preceded the structural changes of the cartilage spatio-temporally. Lastly, while proteoglycan loss was similar in both models, collagen type II degradation only occurred more prominently in MCL-MM. Dynamics of chondrocyte loss and lesion formation in preclinical models has important implications for validating new therapeutic strategies. Our work could be helpful in assessing the feasibility and expected response of the DMM- and the MCL-MM models to chondrocyte mediated therapies. [ABSTRACT FROM AUTHOR]

Abstract: It is well known that many realistic mathematical models of biological and chemical systems, such as enzyme cascades and gene regulatory networks, need to include stochasticity. These systems can be described as Markov processes and are modelled using the Chemical Master Equation (CME). The CME is a differential-difference equation (continuous in time and discrete in the state space) for the probability of certain state at a given time. The state space is the population count of species in the system. A successful method for computing the CME is the Finite State Projection Method (FSP). In this paper we will give the mathematical background of the FSP and propose a new addition to the FSP which guaranties our approximation to have optimal order. [ABSTRACT FROM AUTHOR]

The reaction counts chemical master equation (CME) is a high-dimensional variant of the classical population counts CME. In the reaction counts CME setting, we count the reactions which have fired over time rather than monitoring the population state over time. Since a reaction either fires or not, the reaction counts CME transitions are only forward stepping. Typically there are more reactions in a system than species, this results in the reaction counts CME being higher in dimension, but simpler in dynamics. In this work, we revisit the reaction counts CME framework and its key theoretical results. Then we will extend the theory by exploiting the reactions counts' forward stepping feature, by decomposing the state space into independent continuous-time Markov chains (CTMC). We extend the reaction counts CME theory to derive analytical forms and estimates for the CTMC decomposition of the CME. This new theory gives new insights into solving hitting times-, rare events-, and a priori domain construction problems. [ABSTRACT FROM AUTHOR]

Background: 5q-Spinal muscular atrophy (SMA) is now one of the 5% treatable rare diseases worldwide. As disease-modifying therapies alter disease progression and patient phenotypes, paediatricians and consulting disciplines face new unknowns in their treatment decisions. Conclusions made from historical patient data sets are now mostly limited, and new approaches are needed to ensure our continued best standard-of-care practices for this exceptional patient group. Here, we present a data-driven machine learning approach to a rare disease data set to predict spinal muscular atrophy (SMA)-associated scoliosis., Methods: We collected data from 84 genetically confirmed 5q-SMA patients who have received novel SMA therapies. We performed expert domain knowledge-directed feature engineering, correlation and predictive power score (PPS) analyses for feature selection. To test the predictive performance of the selected features, we trained a Random Forest Classifier and evaluated model performance using standard metrics., Results: The SMA data set consisted of 1304 visits and over 360 variables. We performed feature engineering for variables related to 'interventions', 'devices', 'orthosis', 'ventilation', 'muscle contractures' and 'motor milestones'. Through correlation and PPS analysis paired with expert domain knowledge feature selection, we identified relevant features for scoliosis prediction in SMA that included disease progression markers: Hammersmith Functional Motor Scale Expanded 'HFMSE' (PPS = 0.27) and 6-Minute Walk Test '6MWT' scores (PPS = 0.44), 'age' (PPS = 0.41) and 'weight' (PPS = 0.49), 'contractures' (PPS = 0.17), the use of 'assistive devices' (PPS = 0.39, 'ventilation' (PPS = 0.16) and the presence of 'gastric tubes' (PPS = 0.35) in SMA patients. These features were validated using expert domain knowledge and used to train a Random Forest Classifier with an observed accuracy of 0.82 and an average receiver operating characteristic (ROC) area of 0.87., Conclusion: The introduction of disease-modifying SMA therapies, followed by the implementation of SMA in newborn screenings, has presented physicians with never-seen patients. We used feature engineering tools to overcome one of the main challenges when using data-driven approaches in rare disease data sets. Through predictive modelling of this data, we defined disease progression markers, which are easily assessed during patient visits and can help anticipate scoliosis onset. This highlights the importance of progressive features in the drug-induced revolution of this rare disease and further supports the ongoing efforts to update the SMA classification. We advocate for the consistent documentation of relevant progression markers, which will serve as a basis for data-driven models that physicians can use to update their best standard-of-care practices., (© 2024 The Author(s). Journal of Cachexia, Sarcopenia and Muscle published by Wiley Periodicals LLC.)

Deep convolutional neural networks (DCNNs) are routinely used for image segmentation of biomedical data sets to obtain quantitative measurements of cellular structures like tissues. These cellular structures often contain gaps in their boundaries, leading to poor segmentation performance when using DCNNs like the U-Net. The gaps can usually be corrected by post-hoc computer vision (CV) steps, which are specific to the data set and require a disproportionate amount of work. As DCNNs are Universal Function Approximators, it is conceivable that the corrections should be obsolete by selecting the appropriate architecture for the DCNN. In this article, we present a novel theoretical framework for the gap-filling problem in DCNNs that allows the selection of architecture to circumvent the CV steps. Combining information-theoretic measures of the data set with a fundamental property of DCNNs, the size of their receptive field, allows us to formulate statements about the solvability of the gap-filling problem independent of the specifics of model training. In particular, we obtain mathematical proof showing that the maximum proficiency of filling a gap by a DCNN is achieved if its receptive field is larger than the gap length. We then demonstrate the consequence of this result using numerical experiments on a synthetic and real data set and compare the gap-filling ability of the ubiquitous U-Net architecture with variable depths. Our code is available at https://github.com/ai-biology/dcnn-gap-filling.