Regularized Generalized Canonical Correlation Analysis: A Framework for Sequential Multiblock Component Methods

A framework for a brain-derived nosology of psychiatric disorders

Current psychiatric diagnoses are not defined by neurobiological measures which hinders the development of therapies targeting mechanisms underlying mental illness 1,2 . Research confined to diagnostic boundaries yields heterogeneous biological results, whereas transdiagnostic studies often investigate individual symptoms in isolation. There is currently no paradigm available to comprehensively investigate the relationship between different clinical symptoms, individual disorders, and the underlying neurobiological mechanisms. Here, we propose a framework that groups clinical symptoms derived from ICD-10/DSM-V according to shared brain mechanisms defined by brain structure, function, and connectivity. The reassembly of existing ICD-10/DSM-5 symptoms reveal six cross-diagnostic psychopathology scores related to mania symptoms, depressive symptoms, anxiety symptoms, stress symptoms, eating pathology, and fear symptoms. They were consistently associated with multimodal neuroimaging components in the training sample of young adults aged 23, the independent test sample aged 23, participants aged 14 and 19 years, and in psychiatric patients. The identification of symptom groups of mental illness robustly defined by precisely characterized brain mechanisms enables the development of a psychiatric nosology based upon quantifiable neurobiological measures. As the identified symptom groups align well with existing diagnostic categories, our framework is directly applicable to clinical research and patient care.

DNA methylation haplotype block signatures responding to Staphylococcus aureus subclinical mastitis and association with production and health traits

BACKGROUND

DNA methylation has been documented to play vital roles in diseases and biological processes. In bovine, little is known about the regulatory roles of DNA methylation alterations on production and health traits, including mastitis.

RESULTS

Here, we employed whole-genome DNA methylation sequencing to profile the DNA methylation patterns of milk somatic cells from sixteen cows with naturally occurring Staphylococcus aureus (S. aureus) subclinical mastitis and ten healthy control cows. We observed abundant DNA methylation alterations, including 3,356,456 differentially methylated cytosines and 153,783 differential methylation haplotype blocks (dMHBs). The DNA methylation in regulatory regions, including promoters, first exons and first introns, showed global significant negative correlations with gene expression status. We identified 6435 dMHBs located in the regulatory regions of differentially expressed genes and significantly correlated with their corresponding genes, revealing their potential effects on transcriptional activities. Genes harboring DNA methylation alterations were significantly enriched in multiple immune- and disease-related pathways, suggesting the involvement of DNA methylation in regulating host responses to S. aureus subclinical mastitis. In addition, we found nine discriminant signatures (differentiates cows with S. aureus subclinical mastitis from healthy cows) representing the majority of the DNA methylation variations related to S. aureus subclinical mastitis. Validation of seven dMHBs in 200 cows indicated significant associations with mammary gland health (SCC and SCS) and milk production performance (milk yield).

CONCLUSIONS

In conclusion, our findings revealed abundant DNA methylation alterations in milk somatic cells that may be involved in regulating mammary gland defense against S. aureus infection. Particularly noteworthy is the identification of seven dMHBs showing significant associations with mammary gland health, underscoring their potential as promising epigenetic biomarkers. Overall, our findings on DNA methylation alterations offer novel insights into the regulatory mechanisms of bovine subclinical mastitis, providing further avenues for the development of effective control measures.

Exploring Metabolic and Gut Microbiome Responses to Paraquat Administration in Male Wistar Rats: Implications for Oxidative Stress

In this study, we examined the metabolic and gut microbiome responses to paraquat (PQ) in male Wistar rats, focusing on oxidative stress effects. Rats received a single intraperitoneal injection of PQ at 15 and 30 mg/kg, and various oxidative stress parameters (i.e., MDA, SOD, ROS, 8-isoprostanes) were assessed after three days. To explore the omic profile, GC-qTOF and UHPLC-qTOF were performed to assess the plasma metabolome; 1H-NMR was used to assess the urine metabolome; and shotgun metagenomics sequencing was performed to study the gut microbiome. Our results revealed reductions in body weight and tissue changes, particularly in the liver, were observed, suggesting a systemic effect of PQ. Elevated lipid peroxidation and reactive oxygen species levels in the liver and plasma indicated the induction of oxidative stress. Metabolic profiling revealed changes in the tricarboxylic acid cycle, accumulation of ketone body, and altered levels of key metabolites, such as 3-hydroxybutyric acid and serine, suggesting intricate links between energy metabolism and redox reactions. Plasma metabolomic analysis revealed alterations in mitochondrial metabolism, nicotinamide metabolism, and tryptophan degradation. The gut microbiome showed shifts, with higher PQ doses influencing microbial populations (e.g., Escherichia coli and Akkermansia muciniphila) and metagenomic functions (pyruvate metabolism, fermentation, nucleotide and amino acid biosynthesis). Overall, this study provides comprehensive insights into the complex interplay between PQ exposure, metabolic responses, and gut microbiome dynamics. These findings enhance our understanding of the mechanisms behind oxidative stress-induced metabolic alterations and underscore the connections between xenobiotic exposure, gut microbiota, and host metabolism.

Logistic regression with sparse common and distinctive covariates

Having large sets of predictor variables from multiple sources concerning the same individuals is becoming increasingly common in behavioral research. On top of the variable selection problem, predicting a categorical outcome using such data gives rise to an additional challenge of identifying the processes at play underneath the predictors. These processes are of particular interest in the setting of multi-source data because they can either be associated individually with a single data source or jointly with multiple sources. Although many methods have addressed the classification problem in high dimensionality, the additional challenge of distinguishing such underlying predictor processes from multi-source data has not received sufficient attention. To this end, we propose the method of Sparse Common and Distinctive Covariates Logistic Regression (SCD-Cov-logR). The method is a multi-source extension of principal covariates regression that combines with generalized linear modeling framework to allow classification of a categorical outcome. In a simulation study, SCD-Cov-logR resulted in outperformance compared to related methods commonly used in behavioral sciences. We also demonstrate the practical usage of the method under an empirical dataset.

Multinomial Logistic Factor Regression for Multi-source Functional Block-wise Missing Data

Multi-source functional block-wise missing data arise more commonly in medical care recently with the rapid development of big data and medical technology, hence there is an urgent need to develop efficient dimension reduction to extract important information for classification under such data. However, most existing methods for classification problems consider high-dimensional data as covariates. In the paper, we propose a novel multinomial imputed-factor Logistic regression model with multi-source functional block-wise missing data as covariates. Our main contribution is to establishing two multinomial factor regression models by using the imputed multi-source functional principal component scores and imputed canonical scores as covariates, respectively, where the missing factors are imputed by both the conditional mean imputation and the multiple block-wise imputation approaches. Specifically, the univariate FPCA is carried out for the observable data of each data source firstly to obtain the univariate principal component scores and the eigenfunctions. Then, the block-wise missing univariate principal component scores instead of the block-wise missing functional data are imputed by the conditional mean imputation method and the multiple block-wise imputation method, respectively. After that, based on the imputed univariate factors, the multi-source principal component scores are constructed by using the relationship between the multi-source principal component scores and the univariate principal component scores; and at the same time, the canonical scores are obtained by the multiple-set canonial correlation analysis. Finally, the multinomial imputed-factor Logistic regression model is established with the multi-source principal component scores or the canonical scores as factors. Numerical simulations and real data analysis on ADNI data show the proposed method works well.

Predictive overfitting in immunological applications: Pitfalls and solutions

Overfitting describes the phenomenon where a highly predictive model on the training data generalizes poorly to future observations. It is a common concern when applying machine learning techniques to contemporary medical applications, such as predicting vaccination response and disease status in infectious disease or cancer studies. This review examines the causes of overfitting and offers strategies to counteract it, focusing on model complexity reduction, reliable model evaluation, and harnessing data diversity. Through discussion of the underlying mathematical models and illustrative examples using both synthetic data and published real datasets, our objective is to equip analysts and bioinformaticians with the knowledge and tools necessary to detect and mitigate overfitting in their research.

Multimodal data analysis of knee osteoarthritis assessment: factors selection for conservative care decision making

When assessing a patient with knee osteoarthritis (OA), a number of factors are considered to guide treatment plan, namely, demographic, radiographic, clinical, musculoskeletal, and biomechanical factors. The aim of this study is to identify which of these factors are the most related to each other to potentially better prioritize the modifiable factors to be addressed as they may influence treatment outcomes. We investigated a multimodal canonical correlation analysis to evaluate associations between these factors. The analysis was performed on 415 OA patients who were not candidates for knee arthroplasty, to identify factors that are associated to the patients' clinical conditions.

Missing data in multi-omics integration: Recent advances through artificial intelligence

Biological systems function through complex interactions between various 'omics (biomolecules), and a more complete understanding of these systems is only possible through an integrated, multi-omic perspective. This has presented the need for the development of integration approaches that are able to capture the complex, often non-linear, interactions that define these biological systems and are adapted to the challenges of combining the heterogenous data across 'omic views. A principal challenge to multi-omic integration is missing data because all biomolecules are not measured in all samples. Due to either cost, instrument sensitivity, or other experimental factors, data for a biological sample may be missing for one or more 'omic techologies. Recent methodological developments in artificial intelligence and statistical learning have greatly facilitated the analyses of multi-omics data, however many of these techniques assume access to completely observed data. A subset of these methods incorporate mechanisms for handling partially observed samples, and these methods are the focus of this review. We describe recently developed approaches, noting their primary use cases and highlighting each method's approach to handling missing data. We additionally provide an overview of the more traditional missing data workflows and their limitations; and we discuss potential avenues for further developments as well as how the missing data issue and its current solutions may generalize beyond the multi-omics context.

Multivariate associations between behavioural dimensions and white matter across children and adolescents with and without attention-deficit/hyperactivity disorder

BACKGROUND

Attention deficit/hyperactivity disorder (ADHD) is a heterogeneous neurodevelopmental disorder. Integrity of white matter microstructure plays a key role in the neural mechanism of ADHD presentations. However, the relationships between specific behavioural dimensions and white matter microstructure are less well known. This study aimed to identify associations between white matter and a broad set of clinical features across children and adolescent with and without ADHD using a data-driven multivariate approach.

METHOD

We recruited a total of 130 children (62 controls and 68 ADHD) and employed regularized generalized canonical correlation analysis to characterize the associations between white matter and a comprehensive set of clinical measures covering three domains, including symptom, cognition and behaviour. We further applied linear discriminant analysis to integrate these associations to explore potential developmental effects.

RESULTS

We delineated two brain-behaviour dimensional associations in each domain resulting a total of six multivariate patterns of white matter microstructural alterations linked to hyperactivity-impulsivity and mild affected; executive functions and working memory; externalizing behaviour and social withdrawal, respectively. Apart from executive function and externalizing behaviour sharing similar white matter patterns, all other dimensions linked to a specific pattern of white matter microstructural alterations. The multivariate dimensional association scores showed an overall increase and normalization with age in ADHD group while remained stable in controls.

CONCLUSIONS

We found multivariate neurobehavioral associations exist across ADHD and controls, which suggested that multiple white matter patterns underlie ADHD heterogeneity and provided neural bases for more precise diagnosis and individualized treatment.

Randomized singular value decomposition for integrative subtype analysis of 'omics data' using non-negative matrix factorization

Integration of multiple 'omics datasets for differentiating cancer subtypes is a powerful technic that leverages the consistent and complementary information across multi-omics data. Matrix factorization is a common technique used in integrative clustering for identifying latent subtype structure across multi-omics data. High dimensionality of the omics data and long computation time have been common challenges of clustering methods. In order to address the challenges, we propose randomized singular value decomposition (RSVD) for integrative clustering using Non-negative Matrix Factorization: intNMF-rsvd. The method utilizes RSVD to reduce the dimensionality by projecting the data into eigen vector space with user specified lower rank. Then, clustering analysis is carried out by estimating common basis matrix across the projected multi-omics datasets. The performance of the proposed method was assessed using the simulated datasets and compared with six state-of-the-art integrative clustering methods using real-life datasets from The Cancer Genome Atlas Study. intNMF-rsvd was found working efficiently and competitively as compared to standard intNMF and other multi-omics clustering methods. Most importantly, intNMF-rsvd can handle large number of features and significantly reduce the computation time. The identified subtypes can be utilized for further clinical association studies to understand the etiology of the disease.

A guide to multi-omics data collection and integration for translational medicine

The emerging high-throughput technologies have led to the shift in the design of translational medicine projects towards collecting multi-omics patient samples and, consequently, their integrated analysis. However, the complexity of integrating these datasets has triggered new questions regarding the appropriateness of the available computational methods. Currently, there is no clear consensus on the best combination of omics to include and the data integration methodologies required for their analysis. This article aims to guide the design of multi-omics studies in the field of translational medicine regarding the types of omics and the integration method to choose. We review articles that perform the integration of multiple omics measurements from patient samples. We identify five objectives in translational medicine applications: (i) detect disease-associated molecular patterns, (ii) subtype identification, (iii) diagnosis/prognosis, (iv) drug response prediction, and (v) understand regulatory processes. We describe common trends in the selection of omic types combined for different objectives and diseases. To guide the choice of data integration tools, we group them into the scientific objectives they aim to address. We describe the main computational methods adopted to achieve these objectives and present examples of tools. We compare tools based on how they deal with the computational challenges of data integration and comment on how they perform against predefined objective-specific evaluation criteria. Finally, we discuss examples of tools for downstream analysis and further extraction of novel insights from multi-omics datasets.

High-dimensional gene expression and morphology profiles of cells across 28,000 genetic and chemical perturbations

Cells can be perturbed by various chemical and genetic treatments and the impact on gene expression and morphology can be measured via transcriptomic profiling and image-based assays, respectively. The patterns observed in these high-dimensional profile data can power a dozen applications in drug discovery and basic biology research, but both types of profiles are rarely available for large-scale experiments. Here, we provide a collection of four datasets with both gene expression and morphological profile data useful for developing and testing multimodal methodologies. Roughly a thousand features are measured for each of the two data types, across more than 28,000 chemical and genetic perturbations. We define biological problems that use the shared and complementary information in these two data modalities, provide baseline analysis and evaluation metrics for multi-omic applications, and make the data resource publicly available ( https://broad.io/rosetta/ ).

Network-based integration of multi-omics data for clinical outcome prediction in neuroblastoma

Multi-omics data are increasingly being gathered for investigations of complex diseases such as cancer. However, high dimensionality, small sample size, and heterogeneity of different omics types pose huge challenges to integrated analysis. In this paper, we evaluate two network-based approaches for integration of multi-omics data in an application of clinical outcome prediction of neuroblastoma. We derive Patient Similarity Networks (PSN) as the first step for individual omics data by computing distances among patients from omics features. The fusion of different omics can be investigated in two ways: the network-level fusion is achieved using Similarity Network Fusion algorithm for fusing the PSNs derived for individual omics types; and the feature-level fusion is achieved by fusing the network features obtained from individual PSNs. We demonstrate our methods on two high-risk neuroblastoma datasets from SEQC project and TARGET project. We propose Deep Neural Network and Machine Learning methods with Recursive Feature Elimination as the predictor of survival status of neuroblastoma patients. Our results indicate that network-level fusion outperformed feature-level fusion for integration of different omics data whereas feature-level fusion is more suitable incorporating different feature types derived from same omics type. We conclude that the network-based methods are capable of handling heterogeneity and high dimensionality well in the integration of multi-omics.

Imbalances in TCA, Short Fatty Acids and One-Carbon Metabolisms as Important Features of Homeostatic Disruption Evidenced by a Multi-Omics Integrative Approach of LPS-Induced Chronic Inflammation in Male Wistar Rats

Chronic inflammation is an important risk factor in a broad variety of physical and mental disorders leading to highly prevalent non-communicable diseases (NCDs). However, there is a need for a deeper understanding of this condition and its progression to the disease state. For this reason, it is important to define metabolic pathways and complementary biomarkers associated with homeostatic disruption in chronic inflammation. To achieve that, male Wistar rats were subjected to intraperitoneal and intermittent injections with saline solution or increasing lipopolysaccharide (LPS) concentrations (0.5, 5 and 7.5 mg/kg) thrice a week for 31 days. Biochemical and inflammatory parameters were measured at the end of the study. To assess the omics profile, GC-qTOF and UHPLC-qTOF were performed to evaluate plasma metabolome; 1H-NMR was used to evaluate urine metabolome; additionally, shotgun metagenomics sequencing was carried out to characterize the cecum microbiome. The chronicity of inflammation in the study was evaluated by the monitoring of monocyte chemoattractant protein-1 (MCP-1) during the different weeks of the experimental process. At the end of the study, together with the increased levels of MCP-1, levels of interleukin-6 (IL-6), tumour necrosis factor alpha (TNF-α) and prostaglandin E2 (PGE2) along with 8-isoprostanes (an indicative of oxidative stress) were significantly increased (p-value < 0.05). The leading features implicated in the current model were tricarboxylic acid (TCA) cycle intermediates (i.e., alpha-ketoglutarate, aconitic acid, malic acid, fumaric acid and succinic acid); lipids such as specific cholesterol esters (ChoEs), lysophospholipids (LPCs) and phosphatidylcholines (PCs); and glycine, as well as N, N-dimethylglycine, which are related to one-carbon (1C) metabolism. These metabolites point towards mitochondrial metabolism through TCA cycle, β-oxidation of fatty acids and 1C metabolism as interconnected pathways that could reveal the metabolic effects of chronic inflammation induced by LPS administration. These results provide deeper knowledge concerning the impact of chronic inflammation on the disruption of metabolic homeostasis.

Advancing tools for human early lifecourse exposome research and translation (ATHLETE): Project overview

Early life stages are vulnerable to environmental hazards and present important windows of opportunity for lifelong disease prevention. This makes early life a relevant starting point for exposome studies. The Advancing Tools for Human Early Lifecourse Exposome Research and Translation (ATHLETE) project aims to develop a toolbox of exposome tools and a Europe-wide exposome cohort that will be used to systematically quantify the effects of a wide range of community- and individual-level environmental risk factors on mental, cardiometabolic, and respiratory health outcomes and associated biological pathways, longitudinally from early pregnancy through to adolescence. Exposome tool and data development include as follows: (1) a findable, accessible, interoperable, reusable (FAIR) data infrastructure for early life exposome cohort data, including 16 prospective birth cohorts in 11 European countries; (2) targeted and nontargeted approaches to measure a wide range of environmental exposures (urban, chemical, physical, behavioral, social); (3) advanced statistical and toxicological strategies to analyze complex multidimensional exposome data; (4) estimation of associations between the exposome and early organ development, health trajectories, and biological (metagenomic, metabolomic, epigenetic, aging, and stress) pathways; (5) intervention strategies to improve early life urban and chemical exposomes, co-produced with local communities; and (6) child health impacts and associated costs related to the exposome. Data, tools, and results will be assembled in an openly accessible toolbox, which will provide great opportunities for researchers, policymakers, and other stakeholders, beyond the duration of the project. ATHLETE's results will help to better understand and prevent health damage from environmental exposures and their mixtures from the earliest parts of the life course onward.

Domestication Explains Two-Thirds of Differential-Gene-Expression Variance between Domestic and Wild Animals; The Remaining One-Third Reflects Intraspecific and Interspecific Variation

Belyaev's concept of destabilizing selection during domestication was a major achievement in the XX century. Its practical value has been realized in commercial colors of the domesticated fox that never occur in the wild and has been confirmed in a wide variety of pet breeds. Many human disease models involving animals allow to test drugs before human testing. Perhaps this is why investigators doing transcriptomic profiling of domestic versus wild animals have searched for breed-specific patterns. Here we sequenced hypothalamic transcriptomes of tame and aggressive rats, identified their differentially expressed genes (DEGs), and, for the first time, applied principal component analysis to compare them with all the known DEGs of domestic versus wild animals that we could find. Two principal components, PC1 and PC2, respectively explained 67% and 33% of differential-gene-expression variance (hereinafter: log2 value) between domestic and wild animals. PC1 corresponded to multiple orthologous DEGs supported by homologs; these DEGs kept the log2 value sign from species to species and from tissue to tissue (i.e., a common domestication pattern). PC2 represented stand-alone homologous DEG pairs reversing the log2 value sign from one species to another and from tissue to tissue (i.e., representing intraspecific and interspecific variation).

Integration strategies of multi-omics data for machine learning analysis

Increased availability of high-throughput technologies has generated an ever-growing number of omics data that seek to portray many different but complementary biological layers including genomics, epigenomics, transcriptomics, proteomics, and metabolomics. New insight from these data have been obtained by machine learning algorithms that have produced diagnostic and classification biomarkers. Most biomarkers obtained to date however only include one omic measurement at a time and thus do not take full advantage of recent multi-omics experiments that now capture the entire complexity of biological systems.

Multi-omics data integration strategies are needed to combine the complementary knowledge brought by each omics layer. We have summarized the most recent data integration methods/ frameworks into five different integration strategies: early, mixed, intermediate, late and hierarchical. In this mini-review, we focus on challenges and existing multi-omics integration strategies by paying special attention to machine learning applications.

Application of a path-modelling approach for deciphering causality relationships between microbiota, volatile organic compounds and off-odour profiles during meat spoilage

Microbiological spoilage of meat is considered as a process which involves mainly bacterial metabolism leading to degradation of meat sensory qualities. Studying spoilage requires the collection of different types of experimental data encompassing microbiological, physicochemical and sensorial measurements. Within this framework, the objective herein was to carry out a multiblock path modelling workflow to decipher causality relationships between different types of spoilage-related responses: composition of microbiota, volatilome and off-odour profiles. Analyses were performed with the Path-ComDim approach on a large-scale dataset collected on fresh turkey sausages. This approach enabled to quantify the importance of causality relationships determined a priori between each type of responses as well as to identify important responses involved in spoilage, then to validate causality assumptions. Results were very promising: the data integration confirmed and quantified the causality between data blocks, exhibiting the dynamical nature of spoilage, mainly characterized by the evolution of off-odour profiles caused by the production of volatile organic compounds such as ethanol or ethyl acetate. This production was possibly associated with several bacterial species like Lactococcus piscium, Leuconostoc gelidum, Psychrobacter sp. or Latilactobacillus fuchuensis. Likewise, the production of acetoin and diacetyl in meat spoilage was highlighted. The Path-ComDim approach illustrated here with meat spoilage can be applied to other large-scale and heterogeneous datasets associated with pathway scenarios and represents a promising key tool for deciphering causality in complex biological phenomena.

Multiblock metabolomics: An approach to elucidate whole-body metabolism with multiblock principal component analysis

•

“Multiblock metabolomics” elucidates the global metabolic network in a whole body.

•

“Multiblock metabolomics” combines LC/MS-based metabolomics with multiblock PCA.

•

“Multiblock metabolomics” highlights and elicits organ-specific metabolism.

•

TGs with less unsaturated fatty acids were highly accumulated in the diabetic liver.

Principal component analysis (PCA) is a useful tool for omics analysis to identify underlying factors and visualize relationships between biomarkers. However, this approach is limited in addressing life complexity and further improvement is required. This study aimed to develop a new approach that combines mass spectrometry-based metabolomics with multiblock PCA to elucidate the whole-body global metabolic network, thereby generating comparable metabolite maps to clarify the metabolic relationships among several organs. To evaluate the newly developed method, Zucker diabetic fatty (ZDF) rats (n = 6) were used as type 2 diabetic models and Sprague Dawley (SD) rats (n = 6) as controls. Metabolites in the heart, kidney, and liver were analyzed by capillary electrophoresis and liquid chromatography mass spectrometry, respectively, and the detected metabolites were analyzed by multiblock PCA. More than 300 metabolites were detected in the heart, kidney, and liver. When the metabolites obtained from the three organs were analyzed with multiblock PCA, the score and loading maps obtained were highly synchronized and their metabolism patterns were visually comparable. A significant finding in this study was the different expression patterns in lipid metabolism among the three organs; notably triacylglycerols with polyunsaturated fatty acids or less unsaturated fatty acids showed specific accumulation patterns depending on the organs.

Multi-dimensional data integration algorithm based on random walk with restart

BACKGROUND

The accumulation of various multi-omics data and computational approaches for data integration can accelerate the development of precision medicine. However, the algorithm development for multi-omics data integration remains a pressing challenge.

RESULTS

Here, we propose a multi-omics data integration algorithm based on random walk with restart (RWR) on multiplex network. We call the resulting methodology Random Walk with Restart for multi-dimensional data Fusion (RWRF). RWRF uses similarity network of samples as the basis for integration. It constructs the similarity network for each data type and then connects corresponding samples of multiple similarity networks to create a multiplex sample network. By applying RWR on the multiplex network, RWRF uses stationary probability distribution to fuse similarity networks. We applied RWRF to The Cancer Genome Atlas (TCGA) data to identify subtypes in different cancer data sets. Three types of data (mRNA expression, DNA methylation, and microRNA expression data) are integrated and network clustering is conducted. Experiment results show that RWRF performs better than single data type analysis and previous integrative methods.

CONCLUSIONS

RWRF provides powerful support to users to decipher the cancer molecular subtypes, thus may benefit precision treatment of specific patients in clinical practice.

Multi-omic modelling of inflammatory bowel disease with regularized canonical correlation analysis

BACKGROUND

Personalized medicine requires finding relationships between variables that influence a patient's phenotype and predicting an outcome. Sparse generalized canonical correlation analysis identifies relationships between different groups of variables. This method requires establishing a model of the expected interaction between those variables. Describing these interactions is challenging when the relationship is unknown or when there is no pre-established hypothesis. Thus, our aim was to develop a method to find the relationships between microbiome and host transcriptome data and the relevant clinical variables in a complex disease, such as Crohn's disease.

RESULTS

We present here a method to identify interactions based on canonical correlation analysis. We show that the model is the most important factor to identify relationships between blocks using a dataset of Crohn's disease patients with longitudinal sampling. First the analysis was tested in two previously published datasets: a glioma and a Crohn's disease and ulcerative colitis dataset where we describe how to select the optimum parameters. Using such parameters, we analyzed our Crohn's disease data set. We selected the model with the highest inner average variance explained to identify relationships between transcriptome, gut microbiome and clinically relevant variables. Adding the clinically relevant variables improved the average variance explained by the model compared to multiple co-inertia analysis.

CONCLUSIONS

The methodology described herein provides a general framework for identifying interactions between sets of omic data and clinically relevant variables. Following this method, we found genes and microorganisms that were related to each other independently of the model, while others were specific to the model used. Thus, model selection proved crucial to finding the existing relationships in multi-omics datasets.

Similarity-driven multi-view embeddings from high-dimensional biomedical data

Diverse, high-dimensional modalities collected in large cohorts present new opportunities for the formulation and testing of integrative scientific hypotheses. Similarity-driven multi-view linear reconstruction (SiMLR) is an algorithm that exploits inter-modality relationships to transform large scientific datasets into smaller, more well-powered and interpretable low-dimensional spaces. SiMLR contributes an objective function for identifying joint signal, regularization based on sparse matrices representing prior within-modality relationships and an implementation that permits application to joint reduction of large data matrices. We demonstrate that SiMLR outperforms closely related methods on supervised learning problems in simulation data, a multi-omics cancer survival prediction dataset and multiple modality neuroimaging datasets. Taken together, this collection of results shows that SiMLR may be applied to joint signal estimation from disparate modalities and may yield practically useful results in a variety of application domains.

Benchmarking joint multi-omics dimensionality reduction approaches for the study of cancer

High-dimensional multi-omics data are now standard in biology. They can greatly enhance our understanding of biological systems when effectively integrated. To achieve proper integration, joint Dimensionality Reduction (jDR) methods are among the most efficient approaches. However, several jDR methods are available, urging the need for a comprehensive benchmark with practical guidelines. We perform a systematic evaluation of nine representative jDR methods using three complementary benchmarks. First, we evaluate their performances in retrieving ground-truth sample clustering from simulated multi-omics datasets. Second, we use TCGA cancer data to assess their strengths in predicting survival, clinical annotations and known pathways/biological processes. Finally, we assess their classification of multi-omics single-cell data. From these in-depth comparisons, we observe that intNMF performs best in clustering, while MCIA offers an effective behavior across many contexts. The code developed for this benchmark study is implemented in a Jupyter notebook—multi-omics mix (momix)—to foster reproducibility, and support users and future developers.

Advances in omics technology have resulted in the generation of multi-view data for cancer samples. Here, the authors compare dimensionality reduction techniques using simulated and TCGA data and identify the features of the methods with superior performance.

Global Least Squares Path Modeling: A Full-Information Alternative to Partial Least Squares Path Modeling

Partial least squares path modeling has been widely used for component-based structural equation modeling, where constructs are represented by weighted composites or components of observed variables. This approach remains a limited-information method that carries out two separate stages sequentially to estimate parameters (component weights, loadings, and path coefficients), indicating that it has no single optimization criterion for estimating the parameters at once. In general, limited-information methods are known to provide less efficient parameter estimates than full-information ones. To address this enduring issue, we propose a full-information method for partial least squares path modeling, termed global least squares path modeling, where a single least squares criterion is consistently minimized via a simple iterative algorithm to estimate all the parameters simultaneously. We evaluate the relative performance of the proposed method through the analyses of simulated and real data. We also show that from algorithmic perspectives, the proposed method can be seen as a block-wise special case of another full-information method for component-based structural equation modeling-generalized structured component analysis.

Mortality Risk Profiling of Staphylococcus aureus Bacteremia by Multi-omic Serum Analysis Reveals Early Predictive and Pathogenic Signatures

Staphylococcus aureus bacteremia (SaB) causes significant disease in humans, carrying mortality rates of ∼25%. The ability to rapidly predict SaB patient responses and guide personalized treatment regimens could reduce mortality. Here, we present a resource of SaB prognostic biomarkers. Integrating proteomic and metabolomic techniques enabled the identification of >10,000 features from >200 serum samples collected upon clinical presentation. We interrogated the complexity of serum using multiple computational strategies, which provided a comprehensive view of the early host response to infection. Our biomarkers exceed the predictive capabilities of those previously reported, particularly when used in combination. Last, we validated the biological contribution of mortality-associated pathways using a murine model of SaB. Our findings represent a starting point for the development of a prognostic test for identifying high-risk patients at a time early enough to trigger intensive monitoring and interventions.

Assessing reproducibility of matrix factorization methods in independent transcriptomes

Motivation

Matrix factorization (MF) methods are widely used in order to reduce dimensionality of transcriptomic datasets to the action of few hidden factors (metagenes). MF algorithms have never been compared based on the between-datasets reproducibility of their outputs in similar independent datasets. Lack of this knowledge might have a crucial impact when generalizing the predictions made in a study to others.

Results

We systematically test widely used MF methods on several transcriptomic datasets collected from the same cancer type (14 colorectal, 8 breast and 4 ovarian cancer transcriptomic datasets). Inspired by concepts of evolutionary bioinformatics, we design a novel framework based on Reciprocally Best Hit (RBH) graphs in order to benchmark the MF methods for their ability to produce generalizable components. We show that a particular protocol of application of independent component analysis (ICA), accompanied by a stabilization procedure, leads to a significant increase in the between-datasets reproducibility. Moreover, we show that the signals detected through this method are systematically more interpretable than those of other standard methods. We developed a user-friendly tool for performing the Stabilized ICA-based RBH meta-analysis. We apply this methodology to the study of colorectal cancer (CRC) for which 14 independent transcriptomic datasets can be collected. The resulting RBH graph maps the landscape of interconnected factors associated to biological processes or to technological artifacts. These factors can be used as clinical biomarkers or robust and tumor-type specific transcriptomic signatures of tumoral cells or tumoral microenvironment. Their intensities in different samples shed light on the mechanistic basis of CRC molecular subtyping.

Availability and implementation

The RBH construction tool is available from http://goo.gl/DzpwYp

Supplementary information

Supplementary data are available at Bioinformatics online.

Changes of the human skin microbiota upon chronic exposure to polycyclic aromatic hydrocarbon pollutants

BACKGROUND

Polycyclic aromatic hydrocarbons (PAHs) are of environmental and public health concerns and contribute to adverse skin attributes such as premature skin aging and pigmentary disorder. However, little information is available on the potential roles of chronic urban PAH pollutant exposure on the cutaneous microbiota. Given the roles of the skin microbiota have on healthy and undesirable skin phenotypes and the relationships between PAHs and skin properties, we hypothesize that exposure of PAHs may be associated with changes in the cutaneous microbiota. In this study, the skin microbiota of over two hundred Chinese individuals from two cities in China with varying exposure levels of PAHs were characterized by bacterial and fungal amplicon and shotgun metagenomics sequencing.

RESULTS

Skin site and city were strong parameters in changing microbial communities and their assembly processes. Reductions of bacterial-fungal microbial network structural integrity and stability were associated with skin conditions (acne and dandruff). Multivariate analysis revealed associations between abundances of Propionibacterium and Malassezia with host properties and pollutant exposure levels. Shannon diversity increase was correlated to exposure levels of PAHs in a dose-dependent manner. Shotgun metagenomics analysis of samples (n = 32) from individuals of the lowest and highest exposure levels of PAHs further highlighted associations between the PAHs quantified and decrease in abundances of skin commensals and increase in oral bacteria. Functional analysis identified associations between levels of PAHs and abundance of microbial genes of metabolic and other pathways with potential importance in host-microbe interactions as well as degradation of aromatic compounds.

CONCLUSIONS

The results in this study demonstrated the changes in composition and functional capacities of the cutaneous microbiota associated with chronic exposure levels of PAHs. Findings from this study will aid the development of strategies to harness the microbiota in protecting the skin against pollutants. Video Abstract.

Multiway generalized canonical correlation analysis

Regularized generalized canonical correlation analysis (RGCCA) is a general multiblock data analysis framework that encompasses several important multivariate analysis methods such as principal component analysis, partial least squares regression, and several versions of generalized canonical correlation analysis. In this article, we extend RGCCA to the case where at least one block has a tensor structure. This method is called multiway generalized canonical correlation analysis (MGCCA). Convergence properties of the MGCCA algorithm are studied, and computation of higher-level components are discussed. The usefulness of MGCCA is shown on simulation and on the analysis of a cognitive study in human infants using electroencephalography (EEG).

Correlation and association analyses in microbiome study integrating multiomics in health and disease

Correlation and association analyses are one of the most widely used statistical methods in research fields, including microbiome and integrative multiomics studies. Correlation and association have two implications: dependence and co-occurrence. Microbiome data are structured as phylogenetic tree and have several unique characteristics, including high dimensionality, compositionality, sparsity with excess zeros, and heterogeneity. These unique characteristics cause several statistical issues when analyzing microbiome data and integrating multiomics data, such as large p and small n, dependency, overdispersion, and zero-inflation. In microbiome research, on the one hand, classic correlation and association methods are still applied in real studies and used for the development of new methods; on the other hand, new methods have been developed to target statistical issues arising from unique characteristics of microbiome data. Here, we first provide a comprehensive view of classic and newly developed univariate correlation and association-based methods. We discuss the appropriateness and limitations of using classic methods and demonstrate how the newly developed methods mitigate the issues of microbiome data. Second, we emphasize that concepts of correlation and association analyses have been shifted by introducing network analysis, microbe-metabolite interactions, functional analysis, etc. Third, we introduce multivariate correlation and association-based methods, which are organized by the categories of exploratory, interpretive, and discriminatory analyses and classification methods. Fourth, we focus on the hypothesis testing of univariate and multivariate regression-based association methods, including alpha and beta diversities-based, count-based, and relative abundance (or compositional)-based association analyses. We demonstrate the characteristics and limitations of each approaches. Fifth, we introduce two specific microbiome-based methods: phylogenetic tree-based association analysis and testing for survival outcomes. Sixth, we provide an overall view of longitudinal methods in analysis of microbiome and omics data, which cover standard, static, regression-based time series methods, principal trend analysis, and newly developed univariate overdispersed and zero-inflated as well as multivariate distance/kernel-based longitudinal models. Finally, we comment on current association analysis and future direction of association analysis in microbiome and multiomics studies.

Sparse multiple co-Inertia analysis with application to integrative analysis of multi -Omics data

Background

Multiple co-inertia analysis (mCIA) is a multivariate analysis method that can assess relationships and trends in multiple datasets. Recently it has been used for integrative analysis of multiple high-dimensional -omics datasets. However, its estimated loading vectors are non-sparse, which presents challenges for identifying important features and interpreting analysis results. We propose two new mCIA methods: 1) a sparse mCIA method that produces sparse loading estimates and 2) a structured sparse mCIA method that further enables incorporation of structural information among variables such as those from functional genomics.

Results

Our extensive simulation studies demonstrate the superior performance of the sparse mCIA and structured sparse mCIA methods compared to the existing mCIA in terms of feature selection and estimation accuracy. Application to the integrative analysis of transcriptomics data and proteomics data from a cancer study identified biomarkers that are suggested in the literature related with cancer disease.

Conclusion

Proposed sparse mCIA achieves simultaneous model estimation and feature selection and yields analysis results that are more interpretable than the existing mCIA. Furthermore, proposed structured sparse mCIA can effectively incorporate prior network information among genes, resulting in improved feature selection and enhanced interpretability.

MOTA: Network-Based Multi-Omic Data Integration for Biomarker Discovery

The recent advancement of omic technologies provides researchers with the possibility to search for disease-associated biomarkers at the system level. The integrative analysis of data from a large number of molecules involved at various layers of the biological system offers a great opportunity to rank disease biomarker candidates. In this paper, we propose MOTA, a network-based method that uses data acquired at multiple layers to rank candidate disease biomarkers. The networks constructed by MOTA allow users to investigate the biological significance of the top-ranked biomarker candidates. We evaluated the performance of MOTA in ranking disease-associated molecules from three sets of multi-omic data representing three cohorts of hepatocellular carcinoma (HCC) cases and controls with liver cirrhosis. The results demonstrate that MOTA allows the identification of more top-ranked metabolite biomarker candidates that are shared by two different cohorts compared to traditional statistical methods. Moreover, the mRNA candidates top-ranked by MOTA comprise more cancer driver genes compared to those ranked by traditional differential expression methods.

Identification of miRNA-mRNA associations in hepatocellular carcinoma using hierarchical integrative model

Background

The established role miRNA-mRNA regulation of gene expression has in oncogenesis highlights the importance of integrating miRNA with downstream mRNA targets. These findings call for investigations aimed at identifying disease-associated miRNA-mRNA pairs. Hierarchical integrative models (HIM) offer the opportunity to uncover the relationships between disease and the levels of different molecules measured in multiple omic studies.

Methods

The HIM model we formulated for analysis of mRNA-seq and miRNA-seq data can be specified with two levels: (1) a mechanistic submodel relating mRNAs to miRNAs, and (2) a clinical submodel relating disease status to mRNA and miRNA, while accounting for the mechanistic relationships in the first level.

Results

mRNA-seq and miRNA-seq data were acquired by analysis of tumor and normal liver tissues from 30 patients with hepatocellular carcinoma (HCC). We analyzed the data using HIM and identified 157 significant miRNA-mRNA pairs in HCC. The majority of these molecules have already been independently identified as being either diagnostic, prognostic, or therapeutic biomarker candidates for HCC. These pairs appear to be involved in processes contributing to the pathogenesis of HCC involving inflammation, regulation of cell cycle, apoptosis, and metabolism. For further evaluation of our method, we analyzed miRNA-seq and mRNA-seq data from TCGA network. While some of the miRNA-mRNA pairs we identified by analyzing both our and TCGA data are previously reported in the literature and overlap in regulation and function, new pairs have been identified that may contribute to the discovery of novel targets.

Conclusion

The results strongly support the hypothesis that miRNAs are important regulators of mRNAs in HCC. Furthermore, these results emphasize the biological relevance of studying miRNA-mRNA pairs.

Penalized co-inertia analysis with applications to -omics data

MOTIVATION

Co-inertia analysis (CIA) is a multivariate statistical analysis method that can assess relationships and trends in two sets of data. Recently CIA has been used for an integrative analysis of multiple high-dimensional omics data. However, for classical CIA, all elements in the loading vectors are nonzero, presenting a challenge for the interpretation when analyzing omics data. For other multivariate statistical methods such as canonical correlation analysis (CCA), penalized least squares (PLS), various approaches have been proposed to produce sparse loading vectors via l1-penalization/constraint. We propose a novel CIA method that uses l1-penalization to induce sparsity in estimators of loading vectors. Our method simultaneously conducts model fitting and variable selection. Also, we propose another CIA method that incorporates structure/network information such as those from functional genomics, besides using sparsity penalty so that one can get biologically meaningful and interpretable results.

RESULTS

Extensive simulations demonstrate that our proposed penalized CIA methods achieve the best or close to the best performance compared to the existing CIA method in terms of feature selection and recovery of true loading vectors. Also, we apply our methods to the integrative analysis of gene expression data and protein abundance data from the NCI-60 cancer cell lines. Our analysis of the NCI-60 cancer cell line data reveals meaningful variables for cancer diseases and biologically meaningful results that are consistent with previous studies.

AVAILABILITY AND IMPLEMENTATION

Our algorithms are implemented as an R package which is freely available at: https://www.med.upenn.edu/long-lab/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Does exercise impact gut microbiota composition in men receiving androgen deprivation therapy for prostate cancer? A single-blinded, two-armed, randomised controlled trial

INTRODUCTION

A potential link exists between prostate cancer (PCa) disease and treatment and increased inflammatory levels from gut dysbiosis. This study aims to examine if exercise favourably alters gut microbiota in men receiving androgen deprivation therapy (ADT) for PCa. Specifically, this study will explore whether: (1) exercise improves the composition of gut microbiota and increases the abundance of bacteria associated with health promotion and (2) whether gut health correlates with favourable inflammatory status, bowel function, continence and nausea among patients participating in the exercise intervention.

METHODS AND ANALYSIS

A single-blinded, two-armed, randomised controlled trial will explore the influence of a 3-month exercise programme (3 days/week) for men with high-risk localised PCa receiving ADT. Sixty patients will be randomly assigned to either exercise intervention or usual care. The primary endpoint (gut health and function assessed via feacal samples) and secondary endpoints (self-reported quality of life via standardised questionnaires, blood biomarkers, body composition and physical fitness) will be measured at baseline and following the intervention. A variety of statistical methods will be used to understand the covariance between microbial diversity and metabolomics profile across time and intervention. An intention-to-treat approach will be utilised for the analyses with multiple imputations followed by a secondary sensitivity analysis to ensure data robustness using a complete cases approach.

ETHICS AND DISSEMINATION

Ethics approval was obtained from the Human Research Ethics Committee of Edith Cowan University (ID: 19827 NEWTON). Findings will be reported in peer-reviewed publications and scientific conferences in addition to working with national support groups to translate findings for the broader community. If exercise is shown to result in favourable changes in gut microbial diversity, composition and metabolic profile, and reduce gastrointestinal complications in PCa patients receiving ADT, this study will form the basis of a future phase III trial.

TRIAL REGISTRATION NUMBER

ANZCTR12618000280202.

Multi-block Analysis of Genomic Data Using Generalized Canonical Correlation Analysis

Recently, there have been many studies in medicine related to genetic analysis. Many genetic studies have been performed to find genes associated with complex diseases. To find out how genes are related to disease, we need to understand not only the simple relationship of genotypes but also the way they are related to phenotype. Multi-block data, which is a summation form of variable sets, is used for enhancing the analysis of the relationships of different blocks. By identifying relationships through a multi-block data form, we can understand the association between the blocks in comprehending the correlation between them. Several statistical analysis methods have been developed to understand the relationship between multi-block data. In this paper, we will use generalized canonical correlation methodology to analyze multi-block data from the Korean Association Resource project, which has a combination of single nucleotide polymorphism blocks, phenotype blocks, and disease blocks.

Strategies for integrated analysis in imaging genetics studies

Imaging Genetics (IG) integrates neuroimaging and genomic data from the same individual, deepening our knowledge of the biological mechanisms behind neurodevelopmental domains and neurological disorders. Although the literature on IG has exponentially grown over the past years, the majority of studies have mainly analyzed associations between candidate brain regions and individual genetic variants. However, this strategy is not designed to deal with the complexity of neurobiological mechanisms underlying behavioral and neurodevelopmental domains. Moreover, larger sample sizes and increased multidimensionality of this type of data represents a challenge for standardizing modeling procedures in IG research. This review provides a systematic update of the methods and strategies currently used in IG studies, and serves as an analytical framework for researchers working in this field. To complement the functionalities of the Neuroconductor framework, we also describe existing R packages that implement these methodologies. In addition, we present an overview of how these methodological approaches are applied in integrating neuroimaging and genetic data.