Clustering algorithms in biomedical research: a review

Machine learning algorithms for FPGA Implementation in biomedical engineering applications: A review

Field Programmable Gate Arrays (FPGAs) are integrated circuits that can be configured by the user after manufacturing, making them suitable for customized hardware prototypes, a feature not available in general-purpose processors in Application Specific Integrated Circuits (ASIC). In this paper, we review the vast Machine Learning (ML) algorithms implemented on FPGAs to increase performance and capabilities in healthcare technology over 2001-2023. In particular, we focus on real-time ML algorithms targeted to FPGAs and hybrid System-on-a-chip (SoC) FPGA architectures for biomedical applications. We discuss how previous works have customized and optimized their ML algorithm and FPGA designs to address the putative embedded systems challenges of limited memory, hardware, and power resources while maintaining scalability to accommodate different network sizes and topologies. We provide a synthesis of articles implementing classifiers and regression algorithms, as they are significant algorithms that cover a wide range of ML algorithms used for biomedical applications. This article is written to inform the biomedical engineering and FPGA design communities to advance knowledge of FPGA-enabled ML accelerators for biomedical applications.

PWSC: a novel clustering method based on polynomial weight-adjusted sparse clustering for sparse biomedical data and its application in cancer subtyping

BACKGROUND

Clustering analysis is widely used to interpret biomedical data and uncover new knowledge and patterns. However, conventional clustering methods are not effective when dealing with sparse biomedical data. To overcome this limitation, we propose a hierarchical clustering method called polynomial weight-adjusted sparse clustering (PWSC).

RESULTS

The PWSC algorithm adjusts feature weights using a polynomial function, redefines the distances between samples, and performs hierarchical clustering analysis based on these adjusted distances. Additionally, we incorporate a consensus clustering approach to determine the optimal number of classifications. This consensus approach utilizes relative change in the cumulative distribution function to identify the best number of clusters, resulting in more stable clustering results. Leveraging the PWSC algorithm, we successfully classified a cohort of gastric cancer patients, enabling categorization of patients carrying different types of altered genes. Further evaluation using Entropy showed a significant improvement (p = 2.905e-05), while using the Calinski-Harabasz index demonstrates a remarkable 100% improvement in the quality of the best classification compared to conventional algorithms. Similarly, significantly increased entropy (p = 0.0336) and comparable CHI, were observed when classifying another colorectal cancer cohort with microbial abundance. The above attempts in cancer subtyping demonstrate that PWSC is highly applicable to different types of biomedical data. To facilitate its application, we have developed a user-friendly tool that implements the PWSC algorithm, which canbe accessed at http://pwsc.aiyimed.com/ .

CONCLUSIONS

PWSC addresses the limitations of conventional approaches when clustering sparse biomedical data. By adjusting feature weights and employing consensus clustering, we achieve improved clustering results compared to conventional methods. The PWSC algorithm provides a valuable tool for researchers in the field, enabling more accurate and stable clustering analysis. Its application can enhance our understanding of complex biological systems and contribute to advancements in various biomedical disciplines.

KMD clustering: robust general-purpose clustering of biological data

The noisy and high-dimensional nature of biological data has spawned advanced clustering algorithms that are tailored for specific biological datatypes. However, the performance of such methods varies greatly between datasets and they require post hoc tuning of cryptic hyperparameters. We present k minimal distance (KMD) clustering, a general-purpose method based on a generalization of single and average linkage hierarchical clustering. We introduce a generalized silhouette-like function to eliminate the cryptic hyperparameter k, and use sampling to enable application to million-object datasets. Rigorous comparisons to general and specialized clustering methods on simulated, mass cytometry and scRNA-seq datasets show consistent high performance of KMD clustering across all datasets.

Assessing the surveillance use of 2018 EFP/AAP classification of periodontitis: A validation study and clustering analysis

BACKGROUND

The performance of the 2018 European Federation of Periodontology/American Academy of Periodontology (EFP/AAP) classification of periodontitis for epidemiology surveillance purposes remains to be investigated. This study assessed the surveillance use of the 2018 EFP/AAP classification and its agreement with the unsupervised clustering method compared with the 2012 Centers for Disease Control and Prevention(CDC)/AAP case definition.

METHODS

Participants (n = 9424) in the National Health and Nutrition Examination Survey (NHANES) were staged by the 2018 EFP/AAP classification and classified into subgroups via k-medoids clustering. Concordance levels between periodontitis definitions and the clustering method were evaluated through the multiclass area under the receiver operating characteristic curve (multiclass AUC) among "periodontitis cases" and the general population, respectively. The multiclass AUC of the 2012 CDC/AAP definition versus clustering was used as a reference. The associations of periodontitis with chronic diseases were estimated using multivariable logistic regression.

RESULTS

All the participants were identified as "periodontitis cases" by the 2018 EFP/AAP classification, and the prevalence of stage III-IV was 30%. The optimal numbers of clusters were three and four. The 2012 CDC/AAP definition versus clustering yielded a multiclass AUC of 0.82 and 0.85 among the general population and "periodontitis cases," respectively. The multiclass AUC of the 2018 EFP/AAP classification versus clustering was 0.77 and 0.78 for different target populations. Similar patterns prevailed in associations with chronic diseases between the 2018 EFP/AAP classification and clustering.

CONCLUSIONS

The validity of the 2018 EFP/AAP classification was verified by the unsupervised clustering method, which performed better in distinguishing "periodontitis cases" than classifying the general population. For surveillance purposes, the 2012 CDC/AAP definition showed a higher agreement level with the clustering method than the 2018 EFP/AAP classification.

Two-Dimensional Clustering of Spectral Changes for the Interpretation of Raman Hyperspectra

Two-dimensional correlation spectroscopy (2D-COS) is a technique that permits the examination of synchronous and asynchronous changes present in hyperspectral data. It produces two-dimensional correlation coefficient maps that represent the mutually correlated changes occurring at all Raman wavenumbers during an implemented perturbation. To focus our analysis on clusters of wavenumbers that tend to change together, we apply a k-means clustering to the wavenumber profiles in the perturbation domain decomposition of the two-dimensional correlation coefficient map. These profiles (or trends) reflect peak intensity changes as a function of the perturbation. We then plot the co-occurrences of cluster members two-dimensionally in a manner analogous to a two-dimensional correlation coefficient map. Because wavenumber profiles are clustered based on their similarity, two-dimensional cluster member spectra reveal which Raman peaks change in a similar manner, rather than how much they are correlated. Furthermore, clustering produces a discrete partitioning of the wavenumbers, thus a two-dimensional cluster member spectrum exhibits a discrete presentation of related Raman peaks as opposed to the more continuous representations in a two-dimensional correlation coefficient map. We demonstrate first the basic principles of the technique with the aid of synthetic data. We then apply it to Raman spectra obtained from a polystyrene perchlorate model system followed by Raman spectra from mammalian cells fixed with different percentages of methanol. Both data sets were designed to produce differential changes in sample components. In both cases, all the peaks pertaining to a given component should then change in a similar manner. We observed that component-based profile clustering did occur for polystyrene and perchlorate in the model system and lipids, nucleic acids, and proteins in the mammalian cell example. This confirmed that the method can translate to "real world" samples. We contrast these results with two-dimensional correlation spectroscopy results. To supplement interpretation, we present the cluster-segmented mean spectrum of the hyperspectral data. Overall, this technique is expected to be a valuable adjunct to two-dimensional correlation spectroscopy to further facilitate hyperspectral data interpretation and analysis.

Seismic background level (SBL) growth can reveal slowly developing long-term eruption precursors

The accelerating growth of seismic unrest before eruptions has been observed at many volcanoes and utilized for eruption forecasts. However, there are still many eruptions for which no precursory unrest has been identified, even at well-monitored volcanoes. The recent eruptions of Shinmoe-dake, Japan, have been another negative example of this kind. Here we present seismological evidence that the eruption preparation had been ongoing at the shallow depths beneath Shinmoe-dake for several months to a year. We investigated the seismic background level (SBL) of eleven-year data recorded around the volcano, including two stations about 1 km from the eruptive crater. We searched for persistent weak signals, focusing on low-amplitude time windows recorded during quiet nighttime. Then the spectra of daily background noise were classified by clustering analysis. The SBL analysis successfully revealed very weak precursory tremors from more than several months before the eruption, and residual tremors to the end of the eruptive period. The precursory signals grew acceleratory in a similar way as is assumed in the material failure forecast method applied to eruption forecasts. However, their growth was significantly slower and longer compared to other cases reported in the literature. Such slow and quiet eruption preparations would not be captured by conventional seismological methods. We expect that long-term SBL analyses on proximal seismic data will help detect early precursors, even at seismically quiet volcanoes, and will also help towards judging the end of an eruptive period.

Qluster: An easy-to-implement generic workflow for robust clustering of health data

The exploration of heath data by clustering algorithms allows to better describe the populations of interest by seeking the sub-profiles that compose it. This therefore reinforces medical knowledge, whether it is about a disease or a targeted population in real life. Nevertheless, contrary to the so-called conventional biostatistical methods where numerous guidelines exist, the standardization of data science approaches in clinical research remains a little discussed subject. This results in a significant variability in the execution of data science projects, whether in terms of algorithms used, reliability and credibility of the designed approach. Taking the path of parsimonious and judicious choice of both algorithms and implementations at each stage, this article proposes Qluster, a practical workflow for performing clustering tasks. Indeed, this workflow makes a compromise between (1) genericity of applications (e.g. usable on small or big data, on continuous, categorical or mixed variables, on database of high-dimensionality or not), (2) ease of implementation (need for few packages, few algorithms, few parameters, ...), and (3) robustness (e.g. use of proven algorithms and robust packages, evaluation of the stability of clusters, management of noise and multicollinearity). This workflow can be easily automated and/or routinely applied on a wide range of clustering projects. It can be useful both for data scientists with little experience in the field to make data clustering easier and more robust, and for more experienced data scientists who are looking for a straightforward and reliable solution to routinely perform preliminary data mining. A synthesis of the literature on data clustering as well as the scientific rationale supporting the proposed workflow is also provided. Finally, a detailed application of the workflow on a concrete use case is provided, along with a practical discussion for data scientists. An implementation on the Dataiku platform is available upon request to the authors.

Subtypes of relapsing-remitting multiple sclerosis identified by network analysis

We used network analysis to identify subtypes of relapsing-remitting multiple sclerosis subjects based on their cumulative signs and symptoms. The electronic medical records of 113 subjects with relapsing-remitting multiple sclerosis were reviewed, signs and symptoms were mapped to classes in a neuro-ontology, and classes were collapsed into sixteen superclasses by subsumption. After normalization and vectorization of the data, bipartite (subject-feature) and unipartite (subject-subject) network graphs were created using NetworkX and visualized in Gephi. Degree and weighted degree were calculated for each node. Graphs were partitioned into communities using the modularity score. Feature maps visualized differences in features by community. Network analysis of the unipartite graph yielded a higher modularity score (0.49) than the bipartite graph (0.25). The bipartite network was partitioned into five communities which were named fatigue, behavioral, hypertonia/weakness, abnormal gait/sphincter, and sensory, based on feature characteristics. The unipartite network was partitioned into five communities which were named fatigue, pain, cognitive, sensory, and gait/weakness/hypertonia based on features. Although we did not identify pure subtypes (e.g., pure motor, pure sensory, etc.) in this cohort of multiple sclerosis subjects, we demonstrated that network analysis could partition these subjects into different subtype communities. Larger datasets and additional partitioning algorithms are needed to confirm these findings and elucidate their significance. This study contributes to the literature investigating subtypes of multiple sclerosis by combining feature reduction by subsumption with network analysis.

Bayesian profile regression to study the ecologic associations of correlated environmental exposures with excess mortality risk during the first year of the Covid-19 epidemic in lombardy, Italy

Many countries, including Italy, have experienced significant social and spatial inequalities in mortality during the Covid-19 pandemic. This study applies a multiple exposures framework to investigate how joint place-based factors influence spatial inequalities of excess mortality during the first year of the Covid -19 pandemic in the Lombardy region of Italy. For the Lombardy region, we integrated municipality-level data on all-cause mortality between 2015 and 2020 with 13 spatial covariates, including 5-year average concentrations of six air pollutants, the average temperature in 2020, and multiple socio-demographic factors, and health facilities per capita. Using the clustering algorithm Bayesian profile regression, we fit spatial covariates jointly to identify clusters of municipalities with similar exposure profiles and estimated associations between clusters and excess mortality in 2020. Cluster analysis resulted in 13 clusters. Controlling for spatial autocorrelation of excess mortality and health-protective agency, two clusters had significantly elevated excess mortality than the rest of Lombardy. Municipalities in these highest-risk clusters are in Bergamo, Brescia, and Cremona provinces. The highest risk cluster (C11) had the highest long-term particulate matter air pollution levels (PM2.5 and PM10) and significantly elevated NO2 and CO air pollutants, temperature, proportion ≤18 years, and male-to-female ratio. This cluster is significantly lower for income and ≥65 years. The other high-risk cluster, Cluster 10 (C10), is elevated significantly for ozone but significantly lower for other air pollutants. Covariates with elevated levels for C10 include proportion 65 years or older and a male-to-female ratio. Cluster 10 is significantly lower for income, temperature, per capita health facilities, ≤18 years, and population density. Our results suggest that joint built, natural, and socio-demographic factors influenced spatial inequalities of excess mortality in Lombardy in 2020. Studies must apply a multiple exposures framework to guide policy decisions addressing the complex and multi-dimensional nature of spatial inequalities of Covid-19-related mortality.

Data clustering: application and trends

Clustering has primarily been used as an analytical technique to group unlabeled data for extracting meaningful information. The fact that no clustering algorithm can solve all clustering problems has resulted in the development of several clustering algorithms with diverse applications. We review data clustering, intending to underscore recent applications in selected industrial sectors and other notable concepts. In this paper, we begin by highlighting clustering components and discussing classification terminologies. Furthermore, specific, and general applications of clustering are discussed. Notable concepts on clustering algorithms, emerging variants, measures of similarities/dissimilarities, issues surrounding clustering optimization, validation and data types are outlined. Suggestions are made to emphasize the continued interest in clustering techniques both by scholars and Industry practitioners. Key findings in this review show the size of data as a classification criterion and as data sizes for clustering become larger and varied, the determination of the optimal number of clusters will require new feature extracting methods, validation indices and clustering techniques. In addition, clustering techniques have found growing use in key industry sectors linked to the sustainable development goals such as manufacturing, transportation and logistics, energy, and healthcare, where the use of clustering is more integrated with other analytical techniques than a stand-alone clustering technique.

Analysis of Structured Data in Biomedicine Using Soft Computing Techniques and Computational Analysis

In the field of biomedicine, enormous data are generated in a structured and unstructured form every day. Soft computing techniques play a major role in the interpretation and classification of the data to make appropriate decisions for making policies. The field of medical science and biomedicine needs efficient soft computing-based methods which can process all kind of data such as structured data, categorical data, and unstructured data to generate meaningful outcome for decision-making. The soft-computing methods allow clustering of similar data, classification of data, predictions from big-data analysis, and decision-making on the basis of analysis of data. A novel method is proposed in the paper using soft-computing methods where clustering mechanisms and classification mechanisms are used to process the biomedicine data for productive outcomes. Fuzzy logic and C-means clustering are devised as a collaborative approach to analyze the biomedicine data by reducing the time and space complexity of the clustering solutions. This research work is considering categorical data, numeric data, and structured data for the interpretation of data to make further decisions. Timely decisions are very important especially in the field of biomedicine because human health and human lives are involved in this field and delays in decision-making may cause threats to human lives. The COVID-19 situation was a recent example where timely diagnosis and interpretations played significant roles in saving the lives of people. Therefore, this research work has attempted to use soft computing techniques for the successful clustering of similar medical data and for quicker interpretation of data to support the decision-making processes related to medical fields.

A novel analytical framework for risk stratification of real-world data using machine learning: A small cell lung cancer study

In recent studies, small cell lung cancer (SCLC) treatment guidelines based on Veterans' Administration Lung Study Group limited/extensive disease staging and resulted in broad and inseparable prognostic subgroups. Evidence suggests that the eight versions of tumor, node, and metastasis (TNM) staging can play an important role to address this issue. The aim of the present study was to improve the detection of prognostic subgroups from a real-word data (RWD) cohort of patients and analyze their patterns using a development pipeline with thoracic oncologists and machine learning methods. The method detected subgroups of patients informing unsupervised learning (partition around medoids) including the impact of covariates on prognosis (Cox regression and random survival forest). An analysis was carried out using patients with SCLC (n = 636) with stage IIIA-IVB according to TNM classification. The analysis yielded k = 7 compacted and well-separated clusters of patients. Performance status (Eastern Cooperative Oncology Group-Performance Status), lactate dehydrogenase, spreading of metastasis, cancer stage, and CRP were the baselines that characterized the subgroups. The selected clustering method outperformed standard clustering techniques, which were not capable of detecting meaningful subgroups. From the analysis of cluster treatment decisions, we showed the potential of future RWD applications to understand disease, develop individualized therapies, and improve healthcare decision making.

Identification of an immune-related gene prognostic index for predicting survival and immunotherapy efficacy in papillary renal cell carcinoma

Despite considerable progress has been made in the understanding of the genetics and molecular biology of renal cell carcinoma (RCC), therapeutic options of patients with papillary renal cell carcinoma (PRCC) are limited. Immunotherapy based on immune checkpoint inhibitors (ICIs) has become a hot point in researching new drug for tumor and been tested in a number of human clinical trials. In this study, an immune-related gene prognostic index (IRGPI) was developed and provided a comprehensive and systematic analysis of distinct phenotypic and molecular portraits in the recognition, surveillance, and prognosis of PRCC. The reliability of the IRGPI was evaluated using independent datasets from GEO database and the expression levels of the genes in the IRGPI detected by real-time PCR. Collectively, the currently established IRGPI could be used as a potential biomarker to evaluate the response and efficacy of immunotherapy in PRCC.

Deep multi-modal learning for joint linear representation of nonlinear dynamical systems

Dynamical systems pervasively seen in most real-life applications are complex and behave by following certain evolution rules or dynamical patterns, which are linear, non-linear, or stochastic. The underlying dynamics (or evolution rule) of such complex systems, if found, can be used for understanding the system behavior, and furthermore for system prediction and control. It is common to analyze the system's dynamics through observations in different modality approaches. For instance, to recognize patient deterioration in acute care, it usually relies on monitoring and analyzing vital signs and other observations, such as blood pressure, heart rate, respiration, and electroencephalography. These observations convey the information describing the same target system, but the dynamics is not able to be directly characterized due to high complexity of individual modality and maybe time-delay interactions among modalities. In this work, we suppose that the state behavior of a dynamical system follows an intrinsic dynamics shared among these modalities. We specifically propose a new deep auto-encoder framework using the Koopman operator theory to derive the joint linear dynamics for a target system in a space spanned by the intrinsic coordinates. The proposed method aims to reconstruct the original system states by learning the information provided among multiple modalities. Furthermore, with the derived intrinsic dynamics, our method is capable of restoring the missing observations within and across modalities, and used for predicting the future states of the system that follows the same evolution rule.

Cluster analysis of extracellular matrix biomarkers predicts the development of impaired systolic function within 1 year of acute myocardial infarction

The clinical utility of combining extracellular matrix (ECM) biomarkers to predict the development of impaired systolic function following acute myocardial infarction (AMI) remains largely undetermined. A combination of ELISA and multiplexing assays were performed to measure matrix metalloproteinase (MMP)-2, MMP-3, MMP-8, MMP-9, periostin, N-terminal type I procollagen (PINP) and tissue inhibitor of matrix metalloproteinase-1 (TIMP-1) in plasma samples from 120 AMI patients. All patients had an echocardiogram within 1 year of AMI, and were divided into impaired (n = 37, LVEF < 50%) and preserved (n = 83, LVEF ≥ 50%) systolic function groups. Exploratory factor analysis was performed on log-transformed biomarkers using principle axis analysis with Oblimin rotation. Cluster analysis was performed on log-transformed and normalised biomarkers using Ward’s method of minimum variance and the squared Euclidean distance metric. Upon univariate analysis, current smoking, prescription of ACE inhibitors at discharge, peak hsTnT > 610 ng/L (median), MMP-8 levels, Factor 1 scores and Cluster One assignment were predictive of impaired systolic function. Upon multivariate analysis, Cluster One assignment (odds ratio [95% CI], 2.74 [1.04–7.23], p = 0.04) remained an independent predictor of systolic dysfunction in combination with clinical variables. These observations support the usefulness of combining ECM biomarkers using cluster analysis for predicting the development of impaired systolic function in AMI patients.

Unsupervised Approach for Malignancy Assessment of Lung Nodules in Computed Tomography Scans Using Radiomic Features

Lung cancer is the leading cause of cancer death worldwide. Early low-dose computed tomography (CT) screening can decrease its mortality rate and computer-aided diagnoses systems may make these screenings more accessible. Radiomic features and supervised machine learning have traditionally been employed in these systems. Contrary to supervised methods, unsupervised learning techniques do not require large amounts of annotated data which are labor-intensive to gather and long training times. Therefore, recent approaches have used unsupervised methods, such as clustering, to improve the performance of supervised models. However, an analysis of purely unsupervised methods for malignancy prediction of lung nodules from CT images has not been performed. This work studies nodule malignancy in the LIDC-IDRI image collection of chest CT scans using established radiomic features and unsupervised learning methods based on k-Means, Spectral Clustering, and Gaussian Mixture clustering. All tested methods resulted in clusters of high homogeneity malignancy. Results suggest convex feature distributions and well-separated feature subspaces associated with different diagnoses. Furthermore, diagnosis uncertainty may be explained by common characteristics captured by radiomic features. The k-Means and Gaussian Mixture models are able to generalize to unseen data, achieving a balanced accuracy of 87.23% and 86.96% when inference was tested. These results motivate the usage of unsupervised approaches for malignancy prediction of lung nodules, such as cluster-then-label models. Clinical Relevance- Unsupervised clustering of radiomic features of lung nodules in chest CT scans can differentiate between malignant and benign cases and reflects experts' diagnosis uncertainty.

iFeatureOmega: an integrative platform for engineering, visualization and analysis of features from molecular sequences, structural and ligand data sets

The rapid accumulation of molecular data motivates development of innovative approaches to computationally characterize sequences, structures and functions of biological and chemical molecules in an efficient, accessible and accurate manner. Notwithstanding several computational tools that characterize protein or nucleic acids data, there are no one-stop computational toolkits that comprehensively characterize a wide range of biomolecules. We address this vital need by developing a holistic platform that generates features from sequence and structural data for a diverse collection of molecule types. Our freely available and easy-to-use iFeatureOmega platform generates, analyzes and visualizes 189 representations for biological sequences, structures and ligands. To the best of our knowledge, iFeatureOmega provides the largest scope when directly compared to the current solutions, in terms of the number of feature extraction and analysis approaches and coverage of different molecules. We release three versions of iFeatureOmega including a webserver, command line interface and graphical interface to satisfy needs of experienced bioinformaticians and less computer-savvy biologists and biochemists. With the assistance of iFeatureOmega, users can encode their molecular data into representations that facilitate construction of predictive models and analytical studies. We highlight benefits of iFeatureOmega based on three research applications, demonstrating how it can be used to accelerate and streamline research in bioinformatics, computational biology, and cheminformatics areas. The iFeatureOmega webserver is freely available at http://ifeatureomega.erc.monash.edu and the standalone versions can be downloaded from https://github.com/Superzchen/iFeatureOmega-GUI/ and https://github.com/Superzchen/iFeatureOmega-CLI/.

Over-optimistic evaluation and reporting of novel cluster algorithms: an illustrative study

When researchers publish new cluster algorithms, they usually demonstrate the strengths of their novel approaches by comparing the algorithms’ performance with existing competitors. However, such studies are likely to be optimistically biased towards the new algorithms, as the authors have a vested interest in presenting their method as favorably as possible in order to increase their chances of getting published. Therefore, the superior performance of newly introduced cluster algorithms is over-optimistic and might not be confirmed in independent benchmark studies performed by neutral and unbiased authors. This problem is known among many researchers, but so far, the different mechanisms leading to over-optimism in cluster algorithm evaluation have never been systematically studied and discussed. Researchers are thus often not aware of the full extent of the problem. We present an illustrative study to illuminate the mechanisms by which authors—consciously or unconsciously—paint their cluster algorithm’s performance in an over-optimistic light. Using the recently published cluster algorithm Rock as an example, we demonstrate how optimization of the used datasets or data characteristics, of the algorithm’s parameters and of the choice of the competing cluster algorithms leads to Rock’s performance appearing better than it actually is. Our study is thus a cautionary tale that illustrates how easy it can be for researchers to claim apparent “superiority” of a new cluster algorithm. This illuminates the vital importance of strategies for avoiding the problems of over-optimism (such as, e.g., neutral benchmark studies), which we also discuss in the article.

Latent variable and clustering methods in intersectionality research: systematic review of methods applications

PURPOSE

An intersectionality framework has been increasingly incorporated into quantitative study of health inequity, to incorporate social power in meaningful ways. Researchers have identified "person-centered" methods that cluster within-individual characteristics as appropriate to intersectionality. We aimed to review their use and match with theory.

METHODS

We conducted a multidisciplinary systematic review of English-language quantitative studies wherein authors explicitly stated an intersectional approach, and used clustering methods. We extracted study characteristics and applications of intersectionality.

RESULTS

782 studies with quantitative applications of intersectionality were identified, of which 16 were eligible: eight using latent class analysis, two latent profile analysis, and six clustering methods. Papers used cross-sectional data (100.0%) primarily had U.S. lead authors (68.8%) and were published within psychology, social sciences, and health journals. While 87.5% of papers defined intersectionality and 93.8% cited foundational authors, engagement with intersectionality method literature was more limited. Clustering variables were based on social identities/positions (e.g., gender), dimensions of identity (e.g., race centrality), or processes (e.g., stigma). Results most commonly included four classes/clusters (60.0%), which were frequently used in additional analyses. These described sociodemographic differences across classes/clusters, or used classes/clusters as an exposure variable to predict outcomes in regression analysis, structural equation modeling, mediation, or survival analysis. Author rationales for method choice included both theoretical/intersectional and statistical arguments.

CONCLUSION

Latent variable and clustering methods were used in varied ways in intersectional approaches, and reflected differing matches between theory and methods. We highlight situations in which these methods may be advantageous, and missed opportunities for additional uses.

Analysis of cancer omics data: a selective review of statistical techniques

Cancer is an omics disease. The development in high-throughput profiling has fundamentally changed cancer research and clinical practice. Compared with clinical, demographic and environmental data, the analysis of omics data-which has higher dimensionality, weaker signals and more complex distributional properties-is much more challenging. Developments in the literature are often 'scattered', with individual studies focused on one or a few closely related methods. The goal of this review is to assist cancer researchers with limited statistical expertise in establishing the 'overall framework' of cancer omics data analysis. To facilitate understanding, we mainly focus on intuition, concepts and key steps, and refer readers to the original publications for mathematical details. This review broadly covers unsupervised and supervised analysis, as well as individual-gene-based, gene-set-based and gene-network-based analysis. We also briefly discuss 'special topics' including interaction analysis, multi-datasets analysis and multi-omics analysis.

Automatic Quasi-Clique Merger Algorithm - a Hierarchical Clustering Based on Subgraph-Density

The Automatic Quasi-clique Merger algorithm is a new algorithm adapted from early work published under the name QCM (introduced by Ou and Zhang in 2007). The AQCM algorithm performs hierarchical clustering in any data set for which there is an associated similarity measure quantifying the similarity of any data i and data j. Importantly, the method exhibits two valuable performance properties: 1) the ability to automatically return either a larger or smaller number of clusters depending on the inherent properties of the data rather than on a parameter. 2) the ability to return a very large number of relatively small clusters automatically when such clusters are reasonably well defined in a data set. In this work we present the general idea of a quasi-clique agglomerative approach, provide the full details of the mathematical steps of the AQCM algorithm, and explain some of the motivation behind the new methodology. The main achievement of the new methodology is that the agglomerative process now unfolds adaptively according to the inherent structure unique to a given data set, and this happens without the time-costly parameter adjustment that drove the previous QCM algorithm. For this reason we call the new algorithm automatic. We provide a demonstration of the algorithm's performance at the task of community detection in a social media network of 22,900 nodes.

Machine learning and big data provide crucial insight for future biomaterials discovery and research

Machine learning have been widely adopted in a variety of fields including engineering, science, and medicine revolutionizing how data is collected, used, and stored. Their implementation has led to a drastic increase in the number of computational models for the prediction of various numerical, categorical, or association events given input variables. We aim to examine recent advances in the use of machine learning when applied to the biomaterial field. Specifically, quantitative structure properties relationships offer the unique ability to correlate microscale molecular descriptors to larger macroscale material properties. These new models can be broken down further into four categories: regression, classification, association, and clustering. We examine recent approaches and new uses of machine learning in the three major categories of biomaterials: metals, polymers, and ceramics for rapid property prediction and trend identification. While current research is promising, limitations in the form of lack of standardized reporting and available databases complicates the implementation of described models. Herein, we hope to provide a snapshot of the current state of the field and a beginner's guide to navigating the intersection of biomaterials research and machine learning. STATEMENT OF SIGNIFICANCE: Machine learning and its methods have found a variety of uses beyond the field of computer science but have largely been neglected by those in realm of biomaterials. Through the use of more computational methods, biomaterials development can be expediated while reducing the need for standard trial and error methods. Within, we introduce four basic models that readers can potentially apply to their current research as well as current applications within the field. Furthermore, we hope that this article may act as a "call to action" for readers to realize and address the current lack of implementation within the biomaterials field.

A combined biomarker approach for characterising extracellular matrix profiles in acute myocardial infarction

Extracellular matrix (ECM) biomarkers are useful for measuring underlying molecular activity associated with cardiac repair following acute myocardial infarction (AMI). The aim of this study was to conduct exploratory factor analysis (EFA) to examine the interrelationships between ECM biomarkers, and cluster analysis to identify if distinct ECM profiles could distinguish patient risk in AMI. Ten ECM biomarkers were measured from plasma in 140 AMI patients: MMP-2, -3, -8, -9, periostin, procollagen I N-Terminal propeptide, osteopontin, TGF-β1, TIMP-1 and -4. EFA grouped eight ECM biomarkers into a two-factor solution, which comprised three biomarkers in Factor 1 and five biomarkers in Factor 2. Notably, ECM biomarkers were not separated based on biological function. Cluster analysis grouped AMI patients into three distinct clusters. Cluster One (n = 54) had increased levels of MMP-8, MMP-9, and TGF-B1. Cluster Two (n = 43) had elevated levels of MMP-2, MMP-3, osteopontin, periostin and TIMP-1, and increased high-sensitivity troponin T and GRACE scores. Cluster Three (n = 43) had decreased levels of ECM biomarkers. Circulating ECM biomarkers demonstrated collinearity and entwined biological functions based on EFA analysis. Using cluster analysis, patients with similar clinical presentations could be separated into distinct ECM profiles that were associated with differential patient risk. Clinical significance remains to be determined.

VICTOR: A visual analytics web application for comparing cluster sets

Clustering is the process of grouping different data objects based on similar properties. Clustering has applications in various case studies from several fields such as graph theory, image analysis, pattern recognition, statistics and others. Nowadays, there are numerous algorithms and tools able to generate clustering results. However, different algorithms or parameterizations may produce quite dissimilar cluster sets. In this way, the user is often forced to manually filter and compare these results in order to decide which of them generate the ideal clusters. To automate this process, in this study, we present VICTOR, the first fully interactive and dependency-free visual analytics web application which allows the visual comparison of the results of various clustering algorithms. VICTOR can handle multiple cluster set results simultaneously and compare them using ten different metrics. Clustering results can be filtered and compared to each other with the use of data tables or interactive heatmaps, bar plots, correlation networks, sankey and circos plots. We demonstrate VICTOR's functionality using three examples. In the first case, we compare five different network clustering algorithms on a Yeast protein-protein interaction dataset whereas in the second example, we test four different parameters of the MCL clustering algorithm on the same dataset. Finally, as a third example, we compare four different meta-analyses with hierarchically clustered differentially expressed genes found to be involved in myocardial infarction. VICTOR is available at http://victor.pavlopouloslab.info or http://bib.fleming.gr:3838/VICTOR.

Comparing Statistical Tests for Differential Network Analysis of Gene Modules

Genes often work together to perform complex biological processes, and "networks" provide a versatile framework for representing the interactions between multiple genes. Differential network analysis (DiNA) quantifies how this network structure differs between two or more groups/phenotypes (e.g., disease subjects and healthy controls), with the goal of determining whether differences in network structure can help explain differences between phenotypes. In this paper, we focus on gene co-expression networks, although in principle, the methods studied can be used for DiNA for other types of features (e.g., metabolome, epigenome, microbiome, proteome, etc.). Three common applications of DiNA involve (1) testing whether the connections to a single gene differ between groups, (2) testing whether the connection between a pair of genes differs between groups, or (3) testing whether the connections within a "module" (a subset of 3 or more genes) differs between groups. This article focuses on the latter, as there is a lack of studies comparing statistical methods for identifying differentially co-expressed modules (DCMs). Through extensive simulations, we compare several previously proposed test statistics and a new p-norm difference test (PND). We demonstrate that the true positive rate of the proposed PND test is competitive with and often higher than the other methods, while controlling the false positive rate. The R package discoMod (differentially co-expressed modules) implements the proposed method and provides a full pipeline for identifying DCMs: clustering tools to derive gene modules, tests to identify DCMs, and methods for visualizing the results.

Augmented Two-Dimensional Correlation Spectroscopy for the Joint Analysis of Correlated Changes in Spectroscopic and Disparate Sources

Here, we present an augmented form of two-dimensional correlation spectroscopy, that integrates in a single format data from spectroscopic and multiple non-spectroscopic sources for analysis. The integration is affected by augmenting every spectrum in a hyperspectral data set with relevant non-spectroscopic data to permit two-dimensional correlation analysis^(2D-COS) of the ensemble of augmented spectra. A k-means clustering is then applied to the results of the perturbation domain decomposition to determine which Raman peaks cluster with any of the non-spectroscopic data. We introduce and explain the method with the aid of synthetic spectra and synthetic non-spectroscopic data. We then demonstrate this approach with data using Raman spectra from human embryonic stem cell aggregates undergoing directed differentiation toward pancreatic endocrine cells and parallel bioassays of hormone mRNA expression and C-peptide levels in spent medium. These pancreatic endocrine cells generally contain insulin or glucagon. Insulin has disulfide bonds that produce Raman scattering near 513 cm^-1, but no tryptophan. For insulin-positive cells, we found that the application of multisource correlation analysis revealed a high correlation between insulin mRNA and Raman scattering in the disulfide region. In contrast, glucagon has no disulfide bonds but does contain tryptophan. For glucagon-positive cells, we also observed a high correlation between glucagon mRNA and tryptophan Raman scattering (∼757 cm^-1). We conclude with a discussion of methods to enhance spectral resolution and its effects on the performance of multisource correlation analysis.

High-Throughput Methods in the Discovery and Study of Biomaterials and Materiobiology

The complex interaction of cells with biomaterials (i.e., materiobiology) plays an increasingly pivotal role in the development of novel implants, biomedical devices, and tissue engineering scaffolds to treat diseases, aid in the restoration of bodily functions, construct healthy tissues, or regenerate diseased ones. However, the conventional approaches are incapable of screening the huge amount of potential material parameter combinations to identify the optimal cell responses and involve a combination of serendipity and many series of trial-and-error experiments. For advanced tissue engineering and regenerative medicine, highly efficient and complex bioanalysis platforms are expected to explore the complex interaction of cells with biomaterials using combinatorial approaches that offer desired complex microenvironments during healing, development, and homeostasis. In this review, we first introduce materiobiology and its high-throughput screening (HTS). Then we present an in-depth of the recent progress of 2D/3D HTS platforms (i.e., gradient and microarray) in the principle, preparation, screening for materiobiology, and combination with other advanced technologies. The Compendium for Biomaterial Transcriptomics and high content imaging, computational simulations, and their translation toward commercial and clinical uses are highlighted. In the final section, current challenges and future perspectives are discussed. High-throughput experimentation within the field of materiobiology enables the elucidation of the relationships between biomaterial properties and biological behavior and thereby serves as a potential tool for accelerating the development of high-performance biomaterials.

Simulation-derived best practices for clustering clinical data

INTRODUCTION

Clustering analyses in clinical contexts hold promise to improve the understanding of patient phenotype and disease course in chronic and acute clinical medicine. However, work remains to ensure that solutions are rigorous, valid, and reproducible. In this paper, we evaluate best practices for dissimilarity matrix calculation and clustering on mixed-type, clinical data.

METHODS

We simulate clinical data to represent problems in clinical trials, cohort studies, and EHR data, including single-type datasets (binary, continuous, categorical) and 4 data mixtures. We test 5 single distance metrics (Jaccard, Hamming, Gower, Manhattan, Euclidean) and 3 mixed distance metrics (DAISY, Supersom, and Mercator) with 3 clustering algorithms (hierarchical (HC), k-medoids, self-organizing maps (SOM)). We quantitatively and visually validate by Adjusted Rand Index (ARI) and silhouette width (SW). We applied our best methods to two real-world data sets: (1) 21 features collected on 247 patients with chronic lymphocytic leukemia, and (2) 40 features collected on 6000 patients admitted to an intensive care unit.

RESULTS

HC outperformed k-medoids and SOM by ARI across data types. DAISY produced the highest mean ARI for mixed data types for all mixtures except unbalanced mixtures dominated by continuous data. Compared to other methods, DAISY with HC uncovered superior, separable clusters in both real-world data sets.

DISCUSSION

Selecting an appropriate mixed-type metric allows the investigator to obtain optimal separation of patient clusters and get maximum use of their data. Superior metrics for mixed-type data handle multiple data types using multiple, type-focused distances. Better subclassification of disease opens avenues for targeted treatments, precision medicine, clinical decision support, and improved patient outcomes.

Research of insomnia on traditional Chinese medicine diagnosis and treatment based on machine learning

Background

Insomnia as one of the dominant diseases of traditional Chinese medicine (TCM) has been extensively studied in recent years. To explore the novel approaches of research on TCM diagnosis and treatment, this paper presents a strategy for the research of insomnia based on machine learning.

Methods

First of all, 654 insomnia cases have been collected from an experienced doctor of TCM as sample data. Secondly, in the light of the characteristics of TCM diagnosis and treatment, the contents of research samples have been divided into four parts: the basic information, the four diagnostic methods, the treatment based on syndrome differentiation and the main prescription. And then, these four parts have been analyzed by three analysis methods, including frequency analysis, association rules and hierarchical cluster analysis. Finally, a comprehensive study of the whole four parts has been conducted by random forest.

Results

Researches of the above four parts revealed some essential connections. Simultaneously, based on the algorithm model established by the random forest, the accuracy of predicting the main prescription by the combinations of the four diagnostic methods and the treatment based on syndrome differentiation was 0.85. Furthermore, having been extracted features through applying the random forest, the syndrome differentiation of five zang-organs was proven to be the most significant parameter of the TCM diagnosis and treatment.

Conclusions

The results indicate that the machine learning methods are worthy of being adopted to study the dominant diseases of TCM for exploring the crucial rules of the diagnosis and treatment.

Multi-assignment clustering: Machine learning from a biological perspective

A common approach for analyzing large-scale molecular data is to cluster objects sharing similar characteristics. This assumes that genes with highly similar expression profiles are likely participating in a common molecular process. Biological systems are extremely complex and challenging to understand, with proteins having multiple functions that sometimes need to be activated or expressed in a time-dependent manner. Thus, the strategies applied for clustering of these molecules into groups are of key importance for translation of data to biologically interpretable findings. Here we implemented a multi-assignment clustering (MAsC) approach that allows molecules to be assigned to multiple clusters, rather than single ones as in commonly used clustering techniques. When applied to high-throughput transcriptomics data, MAsC increased power of the downstream pathway analysis and allowed identification of pathways with high biological relevance to the experimental setting and the biological systems studied. Multi-assignment clustering also reduced noise in the clustering partition by excluding genes with a low correlation to all of the resulting clusters. Together, these findings suggest that our methodology facilitates translation of large-scale molecular data into biological knowledge. The method is made available as an R package on GitLab (https://gitlab.com/wolftower/masc).

Transport for the Elderly: Activity Patterns, Mode Choices, and Spatiotemporal Constraints

The rapid aging of the population has posed significant challenges to society and raised new demand for transportation services. Understanding travel needs of the elderly is crucial to making effective strategies for accommodating their demand in many newly motorized cities in developing countries such as China. Using a Markov-chain-based mixture model, we identify two main activity patterns of the elderly: recreation-shopping-oriented (RS-oriented) pattern and schooling-drop-off/pick-up-oriented (SDP-oriented) pattern. Elderly people in the RS-oriented pattern enjoy a cozy life with much time spent on recreation and shopping activities, while those in the SDP-oriented pattern take responsibility of sending grandchildren to school and taking them back home. The RS-oriented elderly people are faced with spatial constraints to access the sparsely distributed recreational sights; however, the SDP-oriented group is subject to temporal constraints when making daily trips. These results would encourage policy makers to reconsider the role of transportation in aged people’s lives and better accommodate their demand through designing safer walking and cycling environment and improving the quality of transit services.

Hypercluster: a flexible tool for parallelized unsupervised clustering optimization

BACKGROUND

Unsupervised clustering is a common and exceptionally useful tool for large biological datasets. However, clustering requires upfront algorithm and hyperparameter selection, which can introduce bias into the final clustering labels. It is therefore advisable to obtain a range of clustering results from multiple models and hyperparameters, which can be cumbersome and slow.

RESULTS

We present hypercluster, a python package and SnakeMake pipeline for flexible and parallelized clustering evaluation and selection. Users can efficiently evaluate a huge range of clustering results from multiple models and hyperparameters to identify an optimal model.

CONCLUSIONS

Hypercluster improves ease of use, robustness and reproducibility for unsupervised clustering application for high throughput biology. Hypercluster is available on pip and bioconda; installation, documentation and example workflows can be found at: https://github.com/ruggleslab/hypercluster .

Evaluation of standard and semantically-augmented distance metrics for neurology patients

Background

Patient distances can be calculated based on signs and symptoms derived from an ontological hierarchy. There is controversy as to whether patient distance metrics that consider the semantic similarity between concepts can outperform standard patient distance metrics that are agnostic to concept similarity. The choice of distance metric can dominate the performance of classification or clustering algorithms. Our objective was to determine if semantically augmented distance metrics would outperform standard metrics on machine learning tasks.

Methods

We converted the neurological findings from 382 published neurology cases into sets of concepts with corresponding machine-readable codes. We calculated patient distances by four different metrics (cosine distance, a semantically augmented cosine distance, Jaccard distance, and a semantically augmented bipartite distance). Semantic augmentation for two of the metrics depended on concept similarities from a hierarchical neuro-ontology. For machine learning algorithms, we used the patient diagnosis as the ground truth label and patient findings as machine learning features. We assessed classification accuracy for four classifiers and cluster quality for two clustering algorithms for each of the distance metrics.

Results

Inter-patient distances were smaller when the distance metric was semantically augmented. Classification accuracy and cluster quality were not significantly different by distance metric.

Conclusion

Although semantic augmentation reduced inter-patient distances, we did not find improved classification accuracy or improved cluster quality with semantically augmented patient distance metrics when applied to a dataset of neurology patients. Further work is needed to assess the utility of semantically augmented patient distances.

Cluster analysis in 975 patients with current cough identifies a phenotype with several cough triggers, many background disorders, and low quality of life

BACKGROUND

Recognition of disorder phenotypes may help to estimate prognosis and to guide the clinical management. Current cough management guidelines classify patients according to the duration of the cough episode. However, this classification is not based on phenotype analyses. The present study aimed to identify cough phenotypes by clustering.

METHODS

An email survey among employed, working-age subjects identified 975 patients with current cough. All filled in a comprehensive 80-item questionnaire including the Leicester Cough Questionnaire. Phenotypes were identified utilizing K-means partitional clustering. A subgroup filled in a follow-up questionnaire 12 months later to investigate the possible differences in the prognosis between the phenotypes.

RESULTS

Two clusters were found. The cluster A included 608 patients (62.4% of the population) and the cluster B 367 patients (37.6%). The three most important variables to separate the clusters were the number of the triggers of cough (mean 2.63 (SD 2.22) vs. 6.95 (2.30), respectively, p < 0.001), the number of the cough background disorders (chronic rhinosinusitis, current asthma, gastroesophageal reflux disease, 0.29 (0.50) vs. 1.28 (0.75), respectively, p < 0.001), and the Leicester Cough Questionnaire physical domain (5.33 (0.76) vs. 4.25 (0.84), respectively, p < 0.001). There were significant interrelationships between these three variables (each p < 0.001). Duration of the episode was not among the most important variables to separate the clusters. At 12 months, 27.0% of the patients of the cluster A and 46.1% of the patients of the cluster B suffered from cough that had continued without interruptions from the first survey (p < 0.001).

CONCLUSIONS

Two cough phenotypes could be identified. Cluster A represents phenotype A, which includes the majority of patients and has a tendency to heal by itself. The authors propose that cluster B represents phenotype TBQ (Triggers, Background disorders, Quality of life impairment). Given the poor prognosis of this phenotype, it urges a prompt and comprehensive clinical evaluation regardless of the duration of the cough episode.

Potential quality improvement of stochastic optical localization nanoscopy images obtained by frame by frame localization algorithms

A data movie of stochastic optical localization nanoscopy contains spatial and temporal correlations, both providing information of emitter locations. The majority of localization algorithms in the literature estimate emitter locations by frame-by-frame localization (FFL), which exploit only the spatial correlation and leave the temporal correlation into the FFL nanoscopy images. The temporal correlation contained in the FFL images, if exploited, can improve the localization accuracy and the image quality. In this paper, we analyze the properties of the FFL images in terms of root mean square minimum distance (RMSMD) and root mean square error (RMSE). It is shown that RMSMD and RMSE can be potentially reduced by a maximum fold equal to the square root of the average number of activations per emitter. Analyzed and revealed are also several statistical properties of RMSMD and RMSE and their relationship with respect to a large number of data frames, bias and variance of localization errors, small localization errors, sample drift, and the worst FFL image. Numerical examples are taken and the results confirm the prediction of analysis. The ideas about how to develop an algorithm to exploit the temporal correlation of FFL images are also briefly discussed. The results suggest development of two kinds of localization algorithms: the algorithms that can exploit the temporal correlation of FFL images and the unbiased localization algorithms.

Toward a Better Understanding of Listeriosis Risk among Older Adults in the United States: Characterizing Dietary Patterns and the Sociodemographic and Economic Attributes of Consumers with These Patterns

ABSTRACT

Older adults are at higher risk of invasive listeriosis compared with the general population. Some foods are more likely than others to be contaminated with or to contain high levels of Listeria monocytogenes. The objectives of this study were to (i) determine dietary consumption patterns among older adults in the United States; (ii) evaluate sociodemographic and economic characteristics of older adults associated with each pattern; (iii) determine intake of foods associated with larger relative risk of listeriosis within these patterns; and (iv) rank these patterns based on risk. Data related to older adults (age 60 and older) participating in the cross-sectional National Health and Nutrition Examination Surveys (NHANES) 2009 to 2010, 2011 to 2012, and 2013 to 2014 (n = 4,967) were included in these analyses. Cluster analysis was used to define dietary patterns based on 24-h dietary recalls from day 1 and day 2. Mean intake of foods associated with higher risk of listeriosis was examined within each pattern, and analysis of variance with Dunnett's method of adjustment was used to evaluate significant differences in mean intake of foods. Patterns were ranked based on relative risk of listeriosis, using outbreak illness attribution data. Five distinct dietary patterns were identified. Patterns ranked at highest relative risk of listeriosis, based on U.S. outbreak illness attribution data, were characterized by relatively higher intakes of fruits, vegetables, and cheeses (∼13% respondents) or cereal, milk, and yogurt (∼14% respondents). Individuals consuming these dietary patterns differed in sex, race, food security, self-rated diet quality, and self-rated health. Cluster analysis, despite methodological limitations, provides new information on consumption, sociodemographic, and economic characteristics of subgroups within susceptible populations, which may be used to target educational messages.

HIGHLIGHTS

Unsupervised machine learning and prognostic factors of survival in chronic lymphocytic leukemia

OBJECTIVE

Unsupervised machine learning approaches hold promise for large-scale clinical data. However, the heterogeneity of clinical data raises new methodological challenges in feature selection, choosing a distance metric that captures biological meaning, and visualization. We hypothesized that clustering could discover prognostic groups from patients with chronic lymphocytic leukemia, a disease that provides biological validation through well-understood outcomes.

METHODS

To address this challenge, we applied k-medoids clustering with 10 distance metrics to 2 experiments ("A" and "B") with mixed clinical features collapsed to binary vectors and visualized with both multidimensional scaling and t-stochastic neighbor embedding. To assess prognostic utility, we performed survival analysis using a Cox proportional hazard model, log-rank test, and Kaplan-Meier curves.

RESULTS

In both experiments, survival analysis revealed a statistically significant association between clusters and survival outcomes (A: overall survival, P = .0164; B: time from diagnosis to treatment, P = .0039). Multidimensional scaling separated clusters along a gradient mirroring the order of overall survival. Longer survival was associated with mutated immunoglobulin heavy-chain variable region gene (IGHV) status, absent Zap 70 expression, female sex, and younger age.

CONCLUSIONS

This approach to mixed-type data handling and selection of distance metric captured well-understood, binary, prognostic markers in chronic lymphocytic leukemia (sex, IGHV mutation status, ZAP70 expression status) with high fidelity.

Artificial intelligence in multiparametric prostate cancer imaging with focus on deep-learning methods

Prostate cancer represents today the most typical example of a pathology whose diagnosis requires multiparametric imaging, a strategy where multiple imaging techniques are combined to reach an acceptable diagnostic performance. However, the reviewing, weighing and coupling of multiple images not only places additional burden on the radiologist, it also complicates the reviewing process. Prostate cancer imaging has therefore been an important target for the development of computer-aided diagnostic (CAD) tools. In this survey, we discuss the advances in CAD for prostate cancer over the last decades with special attention to the deep-learning techniques that have been designed in the last few years. Moreover, we elaborate and compare the methods employed to deliver the CAD output to the operator for further medical decision making.

Screening as a strategy to drive regenerative medicine research

In the field of regenerative medicine, optimization of the parameters leading to a desirable outcome remains a huge challenge. Examples include protocols for the guided differentiation of pluripotent cells towards specialized and functional cell types, phenotypic maintenance of primary cells in cell culture, or engineering of materials for improved tissue interaction with medical implants. This challenge originates from the enormous design space for biomaterials, chemical and biochemical compounds, and incomplete knowledge of the guiding biological principles. To tackle this challenge, high-throughput platforms allow screening of multiple perturbations in one experimental setup. In this review, we provide an overview of screening platforms that are used in regenerative medicine. We discuss their fabrication techniques, and in silico tools to analyze the extensive data sets typically generated by these platforms.

Unsupervised Classification of Multivariate Time Series Using VPCA and Fuzzy Clustering With Spatial Weighted Matrix Distance

Due to high dimensionality and multiple variables, unsupervised classification of multivariate time series (MTS) involves more challenging problems than those of univariate ones. Unlike the vectorization of a feature matrix in traditional clustering algorithms, an unsupervised pattern recognition scheme based on matrix data is proposed for MTS samples in this paper. To reduce the computational load and time consumption, a novel variable-based principal component analysis (VPCA) is first devised for the dimensionality reduction of MTS samples. Afterward, a spatial weighted matrix distance-based fuzzy clustering (SWMDFC) algorithm is proposed to directly group MTS samples into clusters as well as preserve the structure of the data matrix. The spatial weighted matrix distance (SWMD) integrates the spatial dimensionality difference of elements of data into the distance of MST pairs. In terms of the SWMD, the MTS samples are clustered without vectorization in the dimensionality-reduced feature matrix space. Finally, three open-access datasets are utilized for the validation of the proposed unsupervised classification scheme. The results show that the VPCA can capture more features of MTS data than principal component analysis (PCA) and 2-D PCA. Furthermore, the clustering performance of SWMDFC is superior to that of fuzzy c -means clustering algorithms based on the Euclidean distance or image Euclidean distance.

Ginkgolic acid inhibits the growth of renal cell carcinoma cells via inactivation of the EGFR signaling pathway

Renal cell carcinoma (RCC) is one of the most common urological malignancies occurring in adult human kidneys worldwide. Recent research on antitumor drugs has focused on plant extracts, a class of compounds that play critical roles in cancer treatment. The present study aimed to investigate the potential antitumor effect of ginkgolic acid (GA) in RCC. Transwell invasion assay, cell counting kit-8 assay and flow cytometry were used to measure cell migration, cell viability and apoptosis, respectively. A network pharmacology approach was applied to identify pathway information, combining molecular docking techniques to screen for key target information. In the present study, GA inhibited the viability and proliferation of RCC cells (786-O and A498), both in vitro and in vivo, via G₁ arrest. GA also reduced RCC cell invasion and migration. In addition, the epidermal growth factor receptor (EGFR) was identified as a critical target protein of GA, which significantly inactivated EGFR signaling in RCC (P<0.05). Collectively, the present study provided evidence that GA exerts its anticancer function by directly targeting the EGFR signaling pathway, revealing the potential of GA therapy for RCC.

Molecular subtyping of cancer: current status and moving toward clinical applications

Cancer is a collection of genetic diseases, with large phenotypic differences and genetic heterogeneity between different types of cancers and even within the same cancer type. Recent advances in genome-wide profiling provide an opportunity to investigate global molecular changes during the development and progression of cancer. Meanwhile, numerous statistical and machine learning algorithms have been designed for the processing and interpretation of high-throughput molecular data. Molecular subtyping studies have allowed the allocation of cancer into homogeneous groups that are considered to harbor similar molecular and clinical characteristics. Furthermore, this has helped researchers to identify both actionable targets for drug design as well as biomarkers for response prediction. In this review, we introduce five frequently applied techniques for generating molecular data, which are microarray, RNA sequencing, quantitative polymerase chain reaction, NanoString and tissue microarray. Commonly used molecular data for cancer subtyping and clinical applications are discussed. Next, we summarize a workflow for molecular subtyping of cancer, including data preprocessing, cluster analysis, supervised classification and subtype characterizations. Finally, we identify and describe four major challenges in the molecular subtyping of cancer that may preclude clinical implementation. We suggest that standardized methods should be established to help identify intrinsic subgroup signatures and build robust classifiers that pave the way toward stratified treatment of cancer patients.

Bioinformatics Analysis to Screen the Key Prognostic Genes in Tumor Microenvironment of Bladder Cancer

Bladder cancer (BLCA) is the fifth most common cancer and has the features of low survival rate and high morbidity and mortality. The Cancer Genome Atlas (TCGA) is a pool of global gene expression profile and contains huge amounts of cancer genomics data, which makes it possible to inquire the relationship between gene expression and prognosis of a series of malignant tumors including BLCA. Immune and stromal cells are two major components of tumor microenvironment (TME) which play an important role in judging the prognosis of tumor and influencing the progression of malignant, inflammatory, and metabolic disorders. In our study, we conducted a quantitative analysis of immune and stromal elements based on the ESTIMATE algorithm and thus divided BLCA cases into high and low groups. Then the differentially expressed genes closely related to tumor prognosis between groups were identified and had been shown to correlate with immune response and stromal alterations, which was further confirmed by functional enrichment analysis and protein-protein interaction networks. We validated those genes through BLCA dates downloaded from ArrayExpress and thus got the marker genes to predict prognosis of BLCA. Additionally, immune cell infiltration analysis explored the correlation between the verified genes and immune cells. In conclusion, we identified a series of TME-related genes that assess the prognosis and explored the interaction between TME and tumor prognosis to guide clinical individualized treatment.

Hierarchical Clustering Analyses of Plasma Proteins in Subjects With Cardiovascular Risk Factors Identify Informative Subsets Based on Differential Levels of Angiogenic and Inflammatory Biomarkers

Agglomerative hierarchical clustering analysis (HCA) is a commonly used unsupervised machine learning approach for identifying informative natural clusters of observations. HCA is performed by calculating a pairwise dissimilarity matrix and then clustering similar observations until all observations are grouped within a cluster. Verifying the empirical clusters produced by HCA is complex and not well studied in biomedical applications. Here, we demonstrate the comparability of a novel HCA technique with one that was used in previous biomedical applications while applying both techniques to plasma angiogenic (FGF, FLT, PIGF, Tie-2, VEGF, VEGF-D) and inflammatory (MMP1, MMP3, MMP9, IL8, TNFα) protein data to identify informative subsets of individuals. Study subjects were diagnosed with mild cognitive impairment due to cerebrovascular disease (MCI-CVD). Through comparison of the two HCA techniques, we were able to identify subsets of individuals, based on differences in VEGF (p < 0.001), MMP1 (p < 0.001), and IL8 (p < 0.001) levels. These profiles provide novel insights into angiogenic and inflammatory pathologies that may contribute to VCID.