Integrating host and microbiome biology using holo-omics

438 | Mol. Omics, 2024, 20, 438–452 This journal is © The Royal Society of Chemistry 2024 Cite this: Mol. Omics, 2024, 20, 438 Integrating host and microbiome biology using holo-omics Carl M. Kobel, a Jenny Merkesvik, b Idun Maria Tokvam Burgos, c Wanxin Lai, b Ove Øyås, a Phillip B. Pope, abd Torgeir R. Hvidsten b and Velma T. E. Aho * a Holo-omics is the use of omics data to study a host and its inherent microbiomes – a biological system known as a ‘‘holobiont’’. A microbiome that exists in such a space often encounters habitat stability and in return provides metabolic capacities that can benefit their host. Here we present an overview of beneficial host–microbiome systems and propose and discuss several methodological frameworks that can be used to investigate the intricacies of the many as yet undefined host–microbiome interactions that influence holobiont homeostasis. While this is an emerging field, we anticipate that ongoing methodological advancements will enhance the biological resolution that is necessary to improve our understanding of host–microbiome interplay to make meaningful interpretations and biotechnological applications. Introduction Overview and potential of holo-omics In many biological systems and environments, both the host and its resident microbiomes are considered as important contributors to the total function of the overall system. 1 To study the biology of a living system, scientists regularly perform ‘‘omics’’ analyses of the various biomacromolecules that constitute a living cell, such as DNA, RNA, proteins, and metabolites (Table 1, Fig. 1). 2 Analysis of these different ‘‘layers’’ can be further empowered when performed in combination, an approach referred to as ‘‘multi-omics’’. The concept of ‘‘holoomics’’ represents an additional thematic shift whereby instead of focusing on molecular details of an individual organism or an isolated reaction in a specific environment, we can consider the biology of a host–microbiome system as a single unit of action. This makes it possible to understand overarching phenomena in the holobiont. 3 In this context, the objective of a Faculty of Biosciences, Norwegian University of Life Sciences, Ås, Norway. E-mail: velma.tea.[email protected] b Faculty of Chemistry, Biotechnology and Food Science, Norwegian University of Life Sciences, Ås, Norway c Faculty of Natural Sciences, Norwegian University of Science and Technology, Trondheim, Norway d Centre for Microbiome Research, School of Biomedical Sciences, Queensland University of Technology (QUT), Translational Research Institute, Woolloongabba, Queensland, Australia Carl M. Kobel Carl is a bioinformatician and a PhD candidate in the MEMO group at NMBU, Norway. Carl’s perspective is that microbiomes are largely undervalued and that we should better understand the minute interactions within them. Carl adopts a big data inspired approach, enjoys tinkering with hardware, and building parallelizable bioinformatics pipelines to gain insights into large microbiome datasets. Jenny Merkesvik Jenny is a PhD candidate in the Bioinformatics and Applied Statistics group at NMBU. She is part of 3D’omics, a European Union Horizon 2020 project in which she contributes to increase our understanding of host-microbiome interplay through holoomics. Her work is motivated by the aim of improving animal and feed production, benefitting both the animals and the growing human population, in a sustainable way. Received 1st February 2024, Accepted 10th June 2024 DOI: 10.1039/d4mo00017j rsc.li/molomics Molecular Omics REVIEW Open Access Article. Published on 04 July 2024. Downloaded on 10/30/2025 9:29:00 AM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence. View Article Online View Journal | View Issue This journal is © The Royal Society of Chemistry 2024 Mol. Omics, 2024, 20, 438–452 | 439 holo-omics is to study biomacromolecules that constitute biological interactions between a host and its microbiome (Fig. 1). Acquiring a dataset to study host–microbiome interactions is a matter of applying various omics technologies to measure Table 1 Glossary Term Definition Habitat A defined ecological niche that provides environmental parameters that supports a set of organisms. HoloFrom Ancient Greek ‘ o ´loB:ho ´los, ‘‘whole’’. Holo-omics Research that analyses one or more functional layers of omics data from both host and microbiome. The terms holo-omics and hologenomics might be used interchangeably because most omics layers arise from genomic DNA. Holobiont An ecological unit consisting of a host and its resident, interacting micro-organisms. Host–microbiome interface Any surface where biological features from either host or microbiome can interact. Integrative analysis Overlapping or relating the biological factors between two molecular layers or host–microbiome sources. Metagenomics Techniques used to study the collective genomic reads from all organisms in an ecological niche. Multi-omics Research covering more than one omics layer representing one or multiple interacting organisms. Examples of the former include human multi-omics with measurements that only reflect human biology; and microbial multi-omics without taking the host into account. Omics The study of all biomolecules of a specific type. This review focuses on functional omics data, which can be defined as omics data that change over time and across conditions. Proteomics Using a bespoke database which is based on in silico translation of the genomic sequences, to match mass spectrometric spectra to measure the abundance of proteins in a sample. Transcriptomics Techniques used to study an organism’s transcriptome, i.e. the sum of all of its RNA transcripts. Untargeted metabolomics Using methods such as mass spectrometry (MS) or nuclear magnetic resonance (NMR) to measure the abundance of all the metabolites in a sample. Idun Maria Tokvam Burgos Idun Burgos is a PhD candidate in the Systems Biology group at the Norwegian University of Science and Technology. She holds a master’s degree in chemical engineering with a specialization in systems biology from the same university. Her PhD project entails studying bacterial communities through genome-scale metabolic networks, applying methods from bioinformatics, biochemical engineering, and systems biology. Wanxin Lai Wanxin Lai is a PhD candidate affiliated with both the Bioinformatics and Applied Statistics group and the Faculty of Chemistry, Biotechnology and Food Science at the Norwegian University of Life Sciences (NMBU). Her research focuses on integrating multi omics data through networkbased approaches to uncover underlying molecular mechanisms and improve the prediction of phenotypes of interest. Ove Øyås Ove Øyås is a researcher in the Microbial Ecology and MetaOmics group at the Norwegian University of Life Sciences. In 2019, he obtained a PhD in computational systems biology from ETH Zurich, where he developed a passion for understanding biology through model-based integration of omics data. Currently, Dr Øyås is working on multi-omics data analysis and modeling of the rumen gut microbiome as part of the EU-funded HoloRuminant project. Most of his research involves development of scalable computational methods that make it possible to answer new biological questions. Velma T. E. Aho Velma Aho, PhD, has been involved in microbiome research since 2013, starting with amplicon sequencing studies and progressing towards increasingly complex multiand holo-omic projects, constantly striving for a deeper understanding of the roles of microbiomes in mammalian hosts. After ten years of focus on gut microbiota in Parkinson’s disease at the universities of Helsinki and Luxembourg, Dr Aho is currently exploring the microbial community of the cattle rumen as part of the Microbial Ecology and Meta-Omics group at the Norwegian University of Life Sciences. Review Molecular Omics Open Access Article. Published on 04 July 2024. Downloaded on 10/30/2025 9:29:00 AM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence. View Article Online 440 | Mol. Omics, 2024, 20, 438–452 This journal is © The Royal Society of Chemistry 2024 the molecular features of both sides of the holobiont. While this data acquisition used to be the limiting step in such analyses, modern molecular biology tools are making this process more efficient and economical. Today’s primary technical bottlenecks are (1) overcoming microbial community complexity, which can contain thousands of different genomes compared to their singular defined host, and (2) the computational analysis of holo-omic data so that the biological processes of both the host and its microbiome can be integrated computationally, interpreted, and visualised. 7 For example, performing data integration across the host–microbiome interface requires correlating individual biological features across various omics layers, which often cannot be scaled to the typical size of holo-omic datasets and can also suffer due to insufficient statistical power. To meet this challenge a new family of computational tools is needed: they must be able to cluster biological features into modules and cross-correlate features across the host–microbiome boundary, capturing the signals that represent the hypothesised cooperation between the host and its microbiome. Symbiotic interactions in host-associated microbiomes are generally defined by the mutualistic, commensalistic or neutral effects shared between each organism, which depend on whether the benefit involved is one-way, two-way or lacking, respectively. 8 Additionally, there is a spectrum of harmfulneutral interactions within the microbiome and between certain microorganisms, and opportunistic pathogens, viruses, and phages might play a role in defining the dynamics of the microbiome. 9 Additional layers of intra-microbiome complexity should also be considered, particularly for the existence of networks of symbioses within a given microbiome which can be characterised in isolation, as with any other microbial environment. What distinguishes holo-omics is that the host variation is integrated together with any intra-microbiome relationships. Subsequently, holo-omics makes it possible to understand the intra-microbiome dynamics where a hostdirected interaction is imposed on the microbiome. For this review, we discuss in detail how actual holo-omic analyses can be performed computationally and present several frameworks to take the typically massive and complex holoomic datasets and integrate the signal between the host and its microbiome. We consider host–microbiome studies where the host is a multicellular organism like an animal, fungus, or plant that forms a large surface or boundary from which it can interact with the microbiome that typically consists of a community of single-celled microorganisms (bacteria, archaea, eukaryotes) and possibly viruses, with varying degrees of diversity (Table 2). For simplicity, we do not consider parasite interactions in this review but focus on the beneficial interactions in holobionts. Many known hosts are obligate symbionts, meaning the host is non-viable when the microbiome is absent. One example of an obligate holobiont is lichen, where a fungus and a community of cyanobacteria represent a complete holobiont. The fungus provides physical anchoring and nutrient assimilation whereas the cyanobacteria provide carbohydrates assimilated through photosynthesis. Additionally, these holobionts may house Alphaproteobacteria which work in conjunction to fix nitrogen for the lichen, which may otherwise be nutrientlimited. 30 On the other end of the spectrum of dependency are several types of insects, such as ants and caterpillars, which harbour few or no resident microorganisms that are unlikely to have a large impact on fitness. 31 Mammalian hosts tend to fall between these two extreme examples: they are viable when Fig. 1 Holo-omics is a specialised case of multi-omics where biological features are linked across a host–microbiome interface. (A) This interface is idealised along the horizontal axis labelled ‘‘holo-omic’’ as an epithelium with a large surface area where biochemical compounds can be exchanged in both directions. The vertical axis labelled ‘‘multi-omic’’ highlights that interactions can occur on multiple levels in terms of coding sequences and biochemical compounds. (B) Examples of molecular interactions across a host–microbiome interface. 4–6 Created with biorender.com. Molecular Omics Review Open Access Article. Published on 04 July 2024. Downloaded on 10/30/2025 9:29:00 AM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence. View Article Online This journal is © The Royal Society of Chemistry 2024 Mol. Omics, 2024, 20, 438–452 | 441 raised in a germ-free setting, but experimental results suggest various abnormalities in such animals, ranging from changes in the immune system to altered neurodevelopment and behavior. 32,33 Host–microbiome orchestration The holobiont represents an evolutionary shortcut where the host and microbiome partners together orchestrate a metabolic capacity 34,35 that otherwise would have had to develop using horizontal gene transfer and recombination via sexual reproduction within the genome of the host organism itself. 36 In the holobiont perspective, the host provides a habitat for its associated microbiomes with defined and stable ecological factors, such as the presence or gradients of substrates and environmental factors like oxygen, temperature, and H + concentration. 37 In return, the microbiome provides complex biochemicals that the host otherwise would not have been able to synthesise or assimilate. In this context, host-directed internal environmental factors provide the selective pressure that defines which microorganisms are ultimately present. 38 However, many microorganisms, mainly prokaryotes, utilise promiscuous mechanisms for horizontal gene transfer. This gives them the ability to collect mobile and novel genetic elements from diverse sources such as viruses, and alternative genealogies across domains of life. 39 Mechanisms that enable the rapid evolution of microorganisms facilitate their competitive metabolic potential to assimilate both energy and nutrients from a spectrum of ecological niches. A host that has co-evolved with its microbiome can leverage its microbiome-based metabolic potential flexibility to adapt and thrive in niches that the host would have been unlikely to enter on its own. The microbiomes of holobionts are per definition not mediated through the somatic genome of the host which means that the microbiota must have its own way of transmitting genetic material to offspring or between individuals in a population. This means that the composition of species present in a microbiome is subject to change over time as new species colonise and take over functions of others. 40 Host–microbiome co-evolution and adaptation is possible when new microbiota become part of the holobiont in a population of hosts, and are inherited vertically to offspring or between individuals in a host population. This can give rise to endemic microbiota species which are exclusively found as part of a holobiont. The microorganisms can adapt to their host and thus diverge from their ancestral population. Hosts and microbiota are able to co-adapt evolutionarily which means that they can each specialise and optimise their function in the holobiont system over generations. 41,42 Idealised biological frameworks of holo-omic models Studying holobiont systems using holo-omics generally requires a statistical or mechanistic framework that can capture signals or patterns in the data to infer interacting biomacromolecules or biological features across the host and its microbiome. 43 When analysing holo-omic data, its size and complexity usually means it must first be constrained by dimensionality reduction or compression, or by clustering into modular groups of co-abundant biological features. This is to make the computational analysis tractable and to simplify the interpretation of its function. Therefore, it is necessary that the methodological framework chosen to perform this data constraining is able to capture the hypothesised interaction between the host and microbiome. Most frameworks are statistical in the sense that they test whether there are significant differences between treatments or co-appearing groups, but suitable mechanistic models are increasingly available and used for data integration as well. 43 To integrate omics data, these mechanistic models should ideally account for the dynamics of all relevant genome-scale networks in the holobiont system, but scaling to systems of this size entails major computational challenges for dynamic models in particular. 44 Because of this, mechanistic omics integration studies have mainly used genome-scale metabolic models (GEMs), which capture the steady-state flows of metabolites through an organism’s network of biochemical reactions 45 and are available for a range of hosts and microorganisms. 46 By linking metabolic flows to interactions between host and microbiome, GEMs integrated with holo-omics can allow mechanistic investigation of holobiont systems. Dynamic modelling Table 2 Selected examples of host–microbiome systems and their characteristics in terms of symbiotic benefit, dependency, species richness, and services exchanged between host and microbiome. These definitions depend on the ecological circumstances in which each host–microbiome system was considered Holobiont system Symbiosis Microbiome richness Host -microbiome services Microbiome -host services Cattle rumen Mutualistic 8,10 8500–16 994 prokaryotic species, 11,12 52 alveolata, 13 12 fungi 14 Habitat, substrates 15 Catabolism of complex plant fibres, 15 anabolism of essential chemicals Mouse gut Mutualistic 16 828–1573 species 17,18 Habitat, substrates Catabolism of feed matter, anabolism of essential chemicals 16,19 Salmon gut Commensalistic 30–40 species (prokaryotes) 20 Habitat, substrates Unknown Plant root-soil Mutualistic, commensalistic 21 2799–271 940 species 22,23 Energy (sugars, fibres) 24 Nutrients, nitrogen, 25 stress resistance 21 Bee gut Mutualistic o10 species 26–28 Habitat, substrates Modulate social behaviour, 29 catabolism of carbohydrates 28 Review Molecular Omics Open Access Article. Published on 04 July 2024. Downloaded on 10/30/2025 9:29:00 AM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence. View Article Online 442 | Mol. Omics, 2024, 20, 438–452 This journal is © The Royal Society of Chemistry 2024 of genome-scale interaction networks is also becoming feasible thanks to algorithmic and computational advances, 47 butmostof themethodsthatwewilldiscussheretakeastatisticalapproach where they compare and compute significance between groups. Examples of recent publications with a holo-omics approach Since the rise of modern molecular biology tools that have facilitated holo-omic analyses, the number of publications focusing on host–microbiome interactions has been growing. For the purposes of this review, we are particularly interested in studies that include an integrative analysis of two or more omic datasets and discuss both the host and its associated microbiome. Recent holo-omic research articles provide examples of the different types of questions that can be approached from a holo-omics point of view, ranging from experiments with model organisms to comparative evolutionary studies. In the classic experimental end of the spectrum, two studies used a mouse model to address two ‘‘epidemics’’ faced by human medicine: opioid overuse 48 and obesity. 49 Both studies included host transcriptomics, microbial shotgun metagenomics, and untargeted metabolomics, the latter capturing a mix of molecules produced by the host and the microbiome. Their results suggested that the tested medications – morphine in the opioid study, the antidiabetic drug empagliflozin in the obesity study – had effects on the host and microbial layers. 48,49 Both studies further confirmed that there are correlations between different omic layers, offering the simplest kind of evidence for host–microbiome interactions. The opioid study also tested this experimentally by showing that morphine-induced changes in host gene expression vary depending on the presence of a microbiome. 48 In an example closer to traditional ecology, a study focusing on the gut of the termite Labiotermes labralis used metagenomics, metatranscriptomics, and host transcriptomics data to demonstrate that the host and the microbiome provide complementary sets of carbohydrate-active enzymes, enabling the holobiont to degrade a wide range of soil polysaccharides. 50 Finally, a study taking a holo-omics approach to evolution compared several antand termite-eating mammals, with findings that supported convergent evolution not only in host genomes, but also in microbiomes. 51 Specifically, the gut metagenomes of these mammals were enriched in enzymes that are necessary for subsisting on an insectivorous diet, such as chitinases and trehalases, compared to mammals with other types of diets. While the existing publications showcase the exciting opportunities offered by holo-omic research, many of them include only one omics layer for each side of the holobiont. Comprehensive, multi-layered integrative studies remain rare, partly due to financial limitations, but also to the challenges presented by bioinformatic and statistical analyses. State of the art in integrative models Considerations for holo-omics tools Although the cost of generating omics data has come down considerably in recent years, it is still a major undertaking to run controlled animal experiments to obtain matching samples from hosts and microbiomes. As a consequence, holo-omic studies typically tend to have small sample sizes. At the same time, the number of measured biological features (genes, proteins, metabolites) may reach millions, considering that complex microbial communities contain hundreds of species. Let us consider a hypothetical holo-omic study, where we have measured the host transcriptome of the liver in 100 cows (n= 100) and the meta-transcriptome of the rumen content in those same individuals (p= 20 000 host genes + average 3000 microbial genes 200 microbial species = 620 000 features). Let us further assume that the experiment is set up to measure methane emission, and that half of the cows were given a methane-inhibiting feed additive (treatment) that indeed reduced emissions. This dataset would pose a massive challenge for data analysis, and not primarily because it would require considerable computational resources to assemble and annotate Metagenome Assembled Genomes (MAGs) and estimate expression (read mapping). The main challenge is related to the large number of features compared to samples. Naively one would think that this dataset could be analysed using multivariateor machine learning-based prediction methods, where the predictive model could be queried for features or combination of features that contributed significantly to the prediction; ‘‘IF gene G on MAG5 is up AND host gene H is down THEN low methane’’. However, with this many features there will be an enormous number of feature combinations that could separate low and high emitting cows, and with only 100 examples (cows) to constrain them, we would never be able to discern real biological feature-combinations from spurious ones (Fig. 2). This phenomenon is referred to as overfitting and is a consequence of the curse of dimensionality: the number of Fig. 2 Illustrating a common problem in multi-omics and holo-omics where a low number of samples with a high number of features are linked into a low number of traits (methane). The underlying data is arbitrary and represents a single omics layer. Created with biorender.com. Molecular Omics Review Open Access Article. Published on 04 July 2024. Downloaded on 10/30/2025 9:29:00 AM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence. View Article Online This journal is © The Royal Society of Chemistry 2024 Mol. Omics, 2024, 20, 438–452 | 443 examples (cows) needed to identify the biologically meaningful features grows exponentially with the number of features. Methods that divide the aforementioned examples into training and test sets, such as cross validation, would be able to tell us that we are overfitting, but will not be able to solve the problem. Even testing one feature at a time is problematic, since multiple hypothesis testing would severely limit the statistical power and thus only identify features with very large and consistent differences (i.e. large effect sizes) between the two treatments. Luckily, omics features are by no means independent and can be grouped into modules of co-abundant genes, proteins, or metabolites, for instance by correlation. This and other so-called dimensionality reduction approaches typically result in a few dozen distinct modules that can be used as our new features to reveal connections to methane emission and also to hypothesise putative interactions between host and microbiome. A note of caution here is that methods for module finding that rely on computing a distance matrix would require extreme amounts of memory. An approach used for instance by weighted gene co-expression network analysis (WGCNA, a method discussed later in this review) is to first group the data into ‘‘blocks’’ using k-means clustering, find modules in each group, and then combine similar modules at the end. Integrating several omics datasets for a multi-omics approach can help us hone in on biologically meaningful patterns, if done carefully. Assuming that we added metabolomics data to the aforementioned cow example; simply concatenating the transcriptomics and metabolomics table would leave us with even more features (number of genes + number of metabolites). Instead, one could first identify genes and metabolites that are differentially abundant between ‘‘low’’ and ‘‘high’’ methane-emitting cows, and then select pathways that are enriched in both differential genes and differential metabolites. Such consensus integration methods use information about multiple types of molecules to constrain the number of possible biological interpretations. Although there are strong functional interdependencies between rumen microbes converting feed into fatty acids and the host animal metabolizing fatty acids in the liver to produce energy, there are also clear physical boundaries separating these features, meaning that we should consider omic data origins in our holo-omic analysis design. In the case of pathway analysis, for example, one needs to consider that a pathway operates within the confines of a cell of a single organism. More generally, most integration methods are designed for a single species, and thus cannot be applied directly in a holo-omics setting. Any pattern discovered in omics data with the aim of describing host-microbiota interactions must include biomacromolecules originating from both sides of the holobiont boundary. This might be accomplished by first applying a standard (multi-)omics analysis method and then filtering the results afterwards, e.g. selecting modules containing genes from both the host and the microbiota. However, integrating the host-microbiota constraint as an integral part of the data analysis method could drastically reduce the search space, help deal with the curse of dimensionality and force results to include features from the host that might otherwise drown in the sea of microbial features. The methods described below are selected because we find them especially promising for solving challenges related specifically to holo-omics data sets. Existing methodological frameworks and tools Dimensionality reduction. The genetic repertoire of the host and its microbiomes captured by holo-omic data introduces complexities such as data sparsity, sampling variation, ecological differences, and host-specific genetic makeup. Furthermore, distinguishing between free-living and host-associated entities adds another layer of complexity. Since the number of biological features always surpasses the number of observations in holo-omic studies, dimensionality reduction is crucial to create human-interpretable visualisations to explore hidden structures and patterns, and prevent model overfitting. 52 Supervised dimensionality reduction – such as partial least squares discriminant analysis – relies on class labels or response variables to guide the dimensionality reduction process. However, such methods struggle when sample sizes are much smaller than the number of features. On the other hand, unsupervised dimensionality reduction like including matrix factorization and neighbour graph methods, allow discovery of structures in the data without relying on class labels or response variables. 52 Methods that find a few dimensions that are likely to be intrinsic come in two flavours; methods that identify a subset of relevant original features (feature selection), and methods that create new features by combining the original features (feature extraction). Feature extraction methods such as principal component analysis (PCA) (Fig. 3) and single value decomposition utilise variation preservation techniques to extract new features – so-called principal components – that are linear combinations of the original features. Principal components are commonly used for visualising clustering patterns and interpreting sample separation. 53 Canonical correlation analysis (CCA) is a statistical technique akin to PCA in terms of finding a linear transformation of the original variables that consists of orthogonal vectors. 54 The objective of CCA is to summarise the linear relationship between two sets of variables by identifying linear combinations – called canonical variables – that maximise correlations based on pairs of loading vectors. Although CCA is not primarily designed for dimensionality reduction, it plays a crucial role in comprehending multivariate relationships by revealing the directions in which two sets of variables are most interdependent. Several extensions of CCA further enhance its applicability: (i) multiset CCAs analyse maximal correlations across multiple sets of omics data; (ii) sparse CCAs identify a subset of variables most relevant to the canonical variables by introducing sparsity constraints; (iii) regularised CCAs incorporate regularisation which is particularly beneficial when dealing with high-dimensional data or when variables are not wellcaptured by linear transformations; and (iv) partial least squares CCAs which focus on predicting one set of variables using another, thus combining aspects of partial least squares Review Molecular Omics Open Access Article. Published on 04 July 2024. Downloaded on 10/30/2025 9:29:00 AM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence. View Article Online 444 | Mol. Omics, 2024, 20, 438–452 This journal is © The Royal Society of Chemistry 2024 regression with CCA. 55 These extensions cater to diverse scenarios, offering flexibility to address specific challenges in multivariate analysis and canonical correlation. Principal coordinates analysis (PCoA) is a linear transformation method similar to PCA which incorporates multidimensional scaling, creating dissimilarity matrices to visualise sample relationships. 56 Unlike PCA, PCoA is not limited to Euclidean measures and has been shown to be useful for comparing beta-diversity in microbial contexts. Non-metric multidimensional scaling (nMDS) is popular for amplicon/ shotgun sequencing data, offering a rank-based approach that handles non-linear relationships and outliers effectively, albeit with potential distortions in global structures. 57–60 Non-linear methods like t-distributed stochastic neighbour embedding (t-SNE) and uniform manifold approximation and projection (UMAP) belong to the second type of dimensionality reduction, known as neighbour graph algorithms. 59–62 These methods emphasise preserving local structures, relying on graph layout algorithms to create probabilistic weighted graphs representing relationships between high-dimensional data points. UMAP and t-SNE differ primarily in their theoretical foundation for balancing the local and global structures. 53 While t-SNE results can vary between runs due to its stochastic nature and sensitivity to initialisation, UMAP, although also stochastic, tends to demonstrate more stability across runs. UMAP excels in preserving the global structure of the final projection while still capturing local relationships, it is hence a better choice for prediction tasks. 59,60 Nonetheless, it may struggle to distinguish closely nested clusters. It is crucial to note that all three non-linear methods are sensitive to initialisation, and it is recommended to employ the first two principal components from the linear approach as seeds for initiation. Users should implement these exploratory methods with caution, exploring various hyperparameters, running multiple projections for stability. When choosing a non-linear dimensionality reduction method, careful consideration of data scale, characteristics, and specific research goals is essential. 63 Matrix factorisation (NMF and MCFA). Aforementioned methods for dimensionality reduction by matrix factorisation – such as PCA – enable compression of large datasets into a smaller feature space, and may thus facilitate identification of important biological factors for the variation in the observed Fig. 3 Figurative summary of the methods discussed in this review. All can reduce inputs with many features to a smaller number of components in order to simplify interpretation of the underlying biological phenomena. PCA: Principal component analysis; MCFA: Multiset correlation and factor analysis; LDA: Latent dirichlet allocation. Created with biorender.com. Molecular Omics Review Open Access Article. Published on 04 July 2024. Downloaded on 10/30/2025 9:29:00 AM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence. View Article Online This journal is © The Royal Society of Chemistry 2024 Mol. Omics, 2024, 20, 438–452 | 445 data. This is particularly relevant for holo-omic studies utilising a matrix factorisation approach, in which we consider complex systems through assembling a variety of data types from both sides of the holobiont, adding to the already prevalent imbalance of few biological samples and high feature counts. Challenges arise when size and heterogeneity of the dataset increases, which calls for adaptations of these matrix factorisation methods when applied in holo-omics. Non-negative matrix factorisation (NMF) 64 is a method for dimensionality reduction that has been used both in several multi-omic studies and as a basis for additional tools for multiomic data integration and analysis. 65–69 NMF has the same foundation as PCA, essentially decomposing a large data matrix (D) consisting of feature values (p) across biological replicates (n) into a reduced set of (r) linear expressions. These expressions are represented by two matrices smaller than the original data; one with weights (W,pr) and one with the reduced feature components (F,rn) (Fig. 4A). In contrast to PCA, NMF requires the decomposition matrices to contain non-negative values only. This constraint causes the NMF-derived linear expressions to only consist of addends, thereby preventing cancellations between biological factors with opposing signs. NMF thus reflects the idea of assembling parts – analogous to the omic data layers – into a larger image representing the whole system. Simultaneously, the non-negativity constraint of NMF necessitates the compressed data to be seen as an approximation (E) of the real data rather than as an equality (=). 70 Our objective function for determining the decomposition matrices then becomes to minimise the difference between the real data (D) and the approximation (WF). This iterative approach may yield different solutions based on the initial weight and reduced component matrices, potentially affecting the outcome of the analysis. 71 Hence dimensionality reduction by NMF may be more in line with the analogy of assembling omic datasets to uncover interactions between layers of the complex system, although resulting in an approximated model with a potentially large residual difference caused by the lossy factorisation. Another approach to holo-omic dataset integration based on matrix factorisation is multiset correlation and factor analysis (MCFA) 72 (Fig. 3 and 4C). While also seeking to compress observed data (D) into matrices for weights (W) and reduced components (F), MCFA effectively divides the model into two parts. One set of decomposition matrices fit the so-called shared space (S), consisting of reduced features with implied importance across all the included omics layers. This shared space is determined through an extension of CCA called probabilistic CCA (pCCA), and it serves the same purpose as the general decomposition seen in NMF. Additional sets of decomposition matrices are then fitted for each individual omics layer through factor analysis, based on the residual between the read data (D) and the modelled shared space. These ‘‘private’’ aspects of the model reflect contributions from factors that are only perceived as important for observations in specific omics layers. The full model then combines the shared and private spaces to approximate the real data, determining the weight and feature matrices through an expectation maximisation algorithm, with the remainder (c) being quantified a third addend to complete the expression. By fitting the observed data to shared and private reduced features separately, the MCFA method may help distinguish between components with implied importance across all levels of the holobiont and those that only appear relevant for a particular omics layer. Additionally, introducing a private model layer for each omic may leave a smaller residual than had the model only covered components relevant for all included data layers. At the time of writing, MCFA has not been applied in a peer-reviewed study since its publication in August 2023, thus its versatility for holo-omic data integration has yet to be demonstrated. Fig. 4 Comparison of two methods for matrix factorisation; (A) and (B) non-negative matrix factorisation (NMF) and (C) and (D) multiset correlation and factor analysis (MCFA). Both methods reduce a full set of observed data d(columns of D) into linear expressions of reduced features f(columns of F) transformed by multiplication with weights (W). (A) and (C) In contrast to NMF, the MCFA method reduces the dataset into two spaces, either shared between all omics layers (S) or private to each one (P). (B) and (D) All features contribute to approximate the observed data for each shared omics layer, visualised in the same style as Fig. 4 in ref. 66. Review Molecular Omics Open Access Article. Published on 04 July 2024. Downloaded on 10/30/2025 9:29:00 AM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence. View Article Online 446 | Mol. Omics, 2024, 20, 438–452 This journal is © The Royal Society of Chemistry 2024 Network analysis Networks are graphs that represent complex relationships between interacting entities within a system. 73 The network is a ubiquitous concept in informatics that can represent many analogous systems like social interactions, flow of information, internet connections, and biological systems like genome-scale metabolic networks (GEMs), genomic co-occurrences, RNA regulation, protein–protein interactions, and metabolics-driven networks. 74 We use correlation networks as an example in this review, as they are suitable for holo-omic studies (Fig. 3). Multiomics is a more mature concept than holo-omics, hence network methods for the former study type are more developed. 4 We suggest extending these multi-omic tools by integrating data crossing the holo-omic boundary as if it were another omics layer. In general, network analyses handle high-dimensional data well and can provide more interpretable results – compared to other approaches – in the form of node and edge statistics. WGCNA is a popular framework for investigating associations between biological features within a single omics layer. 75 It calculates an adjacency matrix containing transformed, pairwise correlations between biological features such as genes, proteins, and metabolites. The adjacencies are transformed in order to obtain a scale-free network, in which features can be related to continuous and categorical external data like phenotypic traits or treatment groups. On the basis of these adjacencies, the topological overlap measure can be used together with hierarchical clustering to obtain a set of clusters where each biological feature becomes part of only one of the formulated clusters. In WGCNA terminology, these clusters are referred to as ‘‘modules’’ and are represented by their first principal component. This linear combination of biological features is referred to as an eigengene and is idealised to capture the most important variation of the module with limited noise. Since these modules are called without utilising information about treatments or traits, the method can be characterised as unsupervised. WGCNA can be extended to holo-omic data 76 by relating the modules across the host–microbiome boundary. WGCNA has been applied for both clustering and dimensionality reduction in several multi and holo-omic studies related to both plant 77 and animal biology. 76,78–80 One study concerning the gut microbiome in patients with insulin sensitivity or resistance 79 applied a range of node selection and dimensionality reduction methods on their data, and used WGCNA to find clusters of hydrophilic and lipid metabolites. These were later connected to other omics layers to identify clusters associated with metabolism of the gut microbiome between the groups of patients. Alternative clustering methods can also be employed for dimensionality reduction. A state-of-the-art example is the Leiden algorithm, 81 which is an optimisation-based form of clustering. The algorithm was used in a study of HIV patients in which they investigated health in relation to the microbiome of the patients. Specifically, they used the Leiden algorithm to detect clusters of microbiome-derived metabolites before integrating these features with other omics layers. 82 Similarly, a study of the SARS-CoV-2 used the Leiden algorithm to detect clusters of metabolites. 80 Transkingdom network analysis (TkNA) 83 for holoand multi-omics is a network-based method that detects biological features that differentiate treatment groups. TkNA is designed to handle a binary testing condition, such as ‘‘disease’’ and ‘‘control’’. The method consists of a comprehensive pipeline containing all the functions needed to transform normalised data into a network that can be readily visualised. TkNA creates a co-variation network and calculates node statistics like node degrees and bipartite betweenness centrality (BiBC). This approach emphasises that hub nodes with high BiBC and degree represent potential modules of the biological network. Additionally, TkNA interfaces with the Infomap 84 and Louvain 85 network clustering algorithms, which can aid in the interpretation of a biological network further. The size and complexity of networks created from holoomics datasets make them hard to interpret, hence it is necessary to find ways to categorise and structure the represented data. Clustering nodes and thus reducing the number of visual features to consider can help organise the network. This is exemplified in the aforementioned SARS-CoV-2 study where WGCNA was used to recognise clusters across omics layers. The cross-omic clusters correlating with disease severity revealed a relationship between host serum metabolites and microorganisms. 80 In gene set enrichment analysis (GSEA), a gene set usually represents a metabolic pathway that performs a specific biological function. By testing whether there is an enrichment of genes from a specific pathway in a network cluster or module, we can argue that this pathway is captured by the module, thereby drawing further conclusions about its activity by interpreting the module’s omics profile and association to other phenotypic metadata. GSEA can be applied on clusters that are defined using any clustering algorithm. An example is a study on the Atlantic salmon 76 where gene enrichment analysis was used to show that certain host RNA genes responded to long chain fatty acids in the feed. A similar method 86 for improving interpretability is network enrichment, in which functional information and network connectivity is integrated. Instead of testing for a significant difference between treatment groups like GSEA, network enrichment quantifies the differential representation among neighbours in the gene network. 87 A network can be interpreted by statistical concepts that describe crucial properties of the nodes and how they are connected. Degree is simply the number of neighbours of any node. The degree can be expressed relative to the node with the highest number of neighbours, hence degree centrality. Node betweenness describes how many of the pairwise node connections in the network pass through a specific node. If this betweenness measurement is high, the node represents a bottleneck and is indicated to have a potential regulatory effect. 83 The cluster coefficient of a node describes the number of edges between its neighbours in relation to the possible number of edges between these. Coreness considers the neighbourhood of a node as it describes whether a node is part of a Molecular Omics Review Open Access Article. Published on 04 July 2024. Downloaded on 10/30/2025 9:29:00 AM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence. View Article Online