Amalgams: data-driven amalgamation for the dimensionality reduction of compositional data

Evaluation of bi-directional causal association between obstructive sleep apnoea syndrome and diabetic microangiopathy: a Mendelian randomization study

Background

The relationship between obstructive sleep apnea syndrome (OSAS) and diabetic microangiopathy remains controversial.

Objective

This study aimed to use bidirectional two-sample Mendelian Randomization (MR) to assess the causal relationship between OSAS and diabetic microangiopathy.

Methods

First, we used the Linkage Disequilibrium Score Regression(LDSC) analysis to assess the genetic correlation. Then, the bidirectional two-sample MR study was conducted in two stages: OSAS and lung function-related indicators (forced vital capacity (FVC) and forced expiratory volume in 1 s (FEV1)) were investigated as exposures, with diabetic microangiopathy as the outcome in the first stage, and genetic tools were used as proxy variables for OSAS and lung function-related measures in the second step. Genome-wide association study data came from the open GWAS database. We used Inverse-Variance Weighted (IVW), MR-Egger regression, Weighted median, Simple mode, and Weighted mode for effect estimation and pleiotropy testing. We also performed sensitivity analyses to test the robustness of the results. Furthermore, we performed multivariate and mediation MR analyses.

Results

In the LDSC analysis, We found a genetic correlation between OSAS, FVC, FEV 1, and diabetic microangiopathy. In the MR analysis, based on IVW analysis, genetically predicted OSAS was positively correlated with the incidence of diabetic retinopathy (DR), diabetic kidney disease (DKD), and diabetic neuropathy (DN). In the subgroup analysis of DR, there was a significant causal relationship between OSAS and background diabetic retinopathy (BDR) and proliferative diabetic retinopathy (PDR). The reverse MR did not show a correlation between the incidence of diabetic microangiopathy and OSAS. Reduced FVC had a potential causal relationship with increased incidence of DR and PDR. Reduced FEV1 had a potential causal relationship with the increased incidence of BDR, PDR, and DKD. Multivariate MR analysis showed that the association between OSAS and diabetic microangiopathy remained significant after adjusting for confounding factors. However, we did not find the significant mediating factors.

Conclusion

Our results suggest that OSAS may be a cause of the development of diabetic microangiopathy, and OSAS may also be associated with a high risk of diabetic microangiopathy, providing a reference for a better understanding of the prevention of diabetic microangiopathy.

Unveiling the Microbial Realm with VEBA 2.0: A modular bioinformatics suite for end-to-end genome-resolved prokaryotic, (micro)eukaryotic, and viral multi-omics from either short- or long-read sequencing

The microbiome is a complex community of microorganisms, encompassing prokaryotic (bacterial and archaeal), eukaryotic, and viral entities. This microbial ensemble plays a pivotal role in influencing the health and productivity of diverse ecosystems while shaping the web of life. However, many software suites developed to study microbiomes analyze only the prokaryotic community and provide limited to no support for viruses and microeukaryotes. Previously, we introduced the Viral Eukaryotic Bacterial Archaeal (VEBA) open-source software suite to address this critical gap in microbiome research by extending genome-resolved analysis beyond prokaryotes to encompass the understudied realms of eukaryotes and viruses. Here we present VEBA 2.0 with key updates including a comprehensive clustered microeukaryotic protein database, rapid genome/protein-level clustering, bioprospecting, non-coding/organelle gene modeling, genome-resolved taxonomic/pathway profiling, long-read support, and containerization. We demonstrate VEBA's versatile application through the analysis of diverse case studies including marine water, Siberian permafrost, and white-tailed deer lung tissues with the latter showcasing how to identify integrated viruses. VEBA represents a crucial advancement in microbiome research, offering a powerful and accessible platform that bridges the gap between genomics and biotechnological solutions.

Oral mucosal breaks trigger anti-citrullinated bacterial and human protein antibody responses in rheumatoid arthritis

Periodontal disease is more common in individuals with rheumatoid arthritis (RA) who have detectable anti-citrullinated protein antibodies (ACPAs), implicating oral mucosal inflammation in RA pathogenesis. Here, we performed paired analysis of human and bacterial transcriptomics in longitudinal blood samples from RA patients. We found that patients with RA and periodontal disease experienced repeated oral bacteremias associated with transcriptional signatures of ISG15+HLADRhi and CD48highS100A2pos monocytes, recently identified in inflamed RA synovia and blood of those with RA flares. The oral bacteria observed transiently in blood were broadly citrullinated in the mouth, and their in situ citrullinated epitopes were targeted by extensively somatically hypermutated ACPAs encoded by RA blood plasmablasts. Together, these results suggest that (i) periodontal disease results in repeated breaches of the oral mucosa that release citrullinated oral bacteria into circulation, which (ii) activate inflammatory monocyte subsets that are observed in inflamed RA synovia and blood of RA patients with flares and (iii) activate ACPA B cells, thereby promoting affinity maturation and epitope spreading to citrullinated human antigens.

The diagnostic potential and barriers of microbiome based therapeutics

High throughput technological innovations in the past decade have accelerated research into the trillions of commensal microbes in the gut. The 'omics' technologies used for microbiome analysis are constantly evolving, and large-scale datasets are being produced. Despite of the fact that much of the research is still in its early stages, specific microbial signatures have been associated with the promotion of cancer, as well as other diseases such as inflammatory bowel disease, neurogenerative diareses etc. It has been also reported that the diversity of the gut microbiome influences the safety and efficacy of medicines. The availability and declining sequencing costs has rendered the employment of RNA-based diagnostics more common in the microbiome field necessitating improved data-analytical techniques so as to fully exploit all the resulting rich biological datasets, while accounting for their unique characteristics, such as their compositional nature as well their heterogeneity and sparsity. As a result, the gut microbiome is increasingly being demonstrating as an important component of personalised medicine since it not only plays a role in inter-individual variability in health and disease, but it also represents a potentially modifiable entity or feature that may be addressed by treatments in a personalised way. In this context, machine learning and artificial intelligence-based methods may be able to unveil new insights into biomedical analyses through the generation of models that may be used to predict category labels, and continuous values. Furthermore, diagnostic aspects will add value in the identification of the non invasive markers in the critical diseases like cancer.

Principal microbial groups: compositional alternative to phylogenetic grouping of microbiome data

Statistical and machine learning techniques based on relative abundances have been used to predict health conditions and to identify microbial biomarkers. However, high dimensionality, sparsity and the compositional nature of microbiome data represent statistical challenges. On the other hand, the taxon grouping allows summarizing microbiome abundance with a coarser resolution in a lower dimension, but it presents new challenges when correlating taxa with a disease. In this work, we present a novel approach that groups Operational Taxonomical Units (OTUs) based only on relative abundances as an alternative to taxon grouping. The proposed procedure acknowledges the compositional data making use of principal balances. The identified groups are called Principal Microbial Groups (PMGs). The procedure reduces the need for user-defined aggregation of $\textrm{OTU}$s and offers the possibility of working with coarse group of $\textrm{OTU}$s, which are not present in a phylogenetic tree. PMGs can be used for two different goals: (1) as a dimensionality reduction method for compositional data, (2) as an aggregation procedure that provides an alternative to taxon grouping for construction of microbial balances afterward used for disease prediction. We illustrate the procedure with a cirrhosis study data. PMGs provide a coherent data analysis for the search of biomarkers in human microbiota. The source code and demo data for PMGs are available at: https://github.com/asliboyraz/PMGs.

Three approaches to supervised learning for compositional data with pairwise logratios

Logratios between pairs of compositional parts (pairwise logratios) are the easiest to interpret in compositional data analysis, and include the well-known additive logratios as particular cases. When the number of parts is large (sometimes even larger than the number of cases), some form of logratio selection is needed. In this article, we present three alternative stepwise supervised learning methods to select the pairwise logratios that best explain a dependent variable in a generalized linear model, each geared for a specific problem. The first method features unrestricted search, where any pairwise logratio can be selected. This method has a complex interpretation if some pairs of parts in the logratios overlap, but it leads to the most accurate predictions. The second method restricts parts to occur only once, which makes the corresponding logratios intuitively interpretable. The third method uses additive logratios, so that K-1 selected logratios involve a K-part subcomposition. Our approach allows logratios or non-compositional covariates to be forced into the models based on theoretical knowledge, and various stopping criteria are available based on information measures or statistical significance with the Bonferroni correction. We present an application on a dataset from a study predicting Crohn's disease.

Principal Amalgamation Analysis for Microbiome Data

In recent years microbiome studies have become increasingly prevalent and large-scale. Through high-throughput sequencing technologies and well-established analytical pipelines, relative abundance data of operational taxonomic units and their associated taxonomic structures are routinely produced. Since such data can be extremely sparse and high dimensional, there is often a genuine need for dimension reduction to facilitate data visualization and downstream statistical analysis. We propose Principal Amalgamation Analysis (PAA), a novel amalgamation-based and taxonomy-guided dimension reduction paradigm for microbiome data. Our approach aims to aggregate the compositions into a smaller number of principal compositions, guided by the available taxonomic structure, by minimizing a properly measured loss of information. The choice of the loss function is flexible and can be based on familiar diversity indices for preserving either within-sample or between-sample diversity in the data. To enable scalable computation, we develop a hierarchical PAA algorithm to trace the entire trajectory of successive simple amalgamations. Visualization tools including dendrogram, scree plot, and ordination plot are developed. The effectiveness of PAA is demonstrated using gut microbiome data from a preterm infant study and an HIV infection study.

The relationship between land cover and microbial community composition in European lakes

Microbes are essential for element cycling and ecosystem functioning. However, many questions central to understanding the role of microbes in ecology are still open. Here, we analyze the relationship between lake microbiomes and the lakes' land cover. By applying machine learning methods, we quantify the covariance between land cover categories and the microbial community composition recorded in the largest amplicon sequencing dataset of European lakes available to date. Our results show that the aggregation of environmental features or microbial taxa before analysis can obscure ecologically relevant patterns. We observe a comparatively high covariation of the lakes' microbial community with herbaceous and open spaces surrounding the lake; nevertheless, the microbial covariation with land cover categories is generally lower than the covariation with physico-chemical parameters. Combining land cover and physico-chemical bioindicators identified from the same amplicon sequencing dataset, we develop analytical data structures that facilitate insights into the ecology of the lake microbiome. Among these, a list of the environmental parameters sorted by the number of microbial bioindicators we have identified for them points towards apparent environmental drivers of the lake microbial community composition, such as the altitude, conductivity, and area covered herbaceous vegetation surrounding the lake. Furthermore, the response map, a similarity matrix calculated from the Jaccard similarity of the environmental parameters' lists of bioindicators, allows us to study the ecosystem's structure from the standpoint of the microbiome. More specifically, we identify multiple clusters of highly similar and possibly functionally linked ecological parameters, including one that highlights the importance of the calcium-bicarbonate equilibrium for lake ecology. Taken together, we demonstrate the use of machine learning approaches in studying the interplay between microbial diversity and environmental factors and introduce novel approaches to integrate environmental molecular diversity into monitoring and water quality assessments.

Approximation of a Microbiome Composition Shift by a Change in a Single Balance Between Two Groups of Taxa

The method proposed here extends the range of compositional methods providing interpretation of classical statistical tools applied to data converted to the ILR coordinates. It provides a strictly optimal solution in several special cases. The approach is universally applicable to compositional data of any nature, including microbiome data sets.

Environmental DNA metabarcoding for benthic monitoring: A review of sediment sampling and DNA extraction methods

Environmental DNA (eDNA) metabarcoding (parallel sequencing of DNA/RNA for identification of whole communities within a targeted group) is revolutionizing the field of aquatic biomonitoring. To date, most metabarcoding studies aiming to assess the ecological status of aquatic ecosystems have focused on water eDNA and macroinvertebrate bulk samples. However, the eDNA metabarcoding has also been applied to soft sediment samples, mainly for assessing microbial or meiofaunal biota. Compared to classical methodologies based on manual sorting and morphological identification of benthic taxa, eDNA metabarcoding offers potentially important advantages for assessing the environmental quality of sediments. The methods and protocols utilized for sediment eDNA metabarcoding can vary considerably among studies, and standardization efforts are needed to improve their robustness, comparability and use within regulatory frameworks. Here, we review the available information on eDNA metabarcoding applied to sediment samples, with a focus on sampling, preservation, and DNA extraction steps. We discuss challenges specific to sediment eDNA analysis, including the variety of different sources and states of eDNA and its persistence in the sediment. This paper aims to identify good-practice strategies and facilitate method harmonization for routine use of sediment eDNA in future benthic monitoring.

A prospective investigation into the association between the gut microbiome composition and cognitive performance among healthy young adults

Background

There is emerging evidence that the gut microbiome composition is associated with several human health outcomes, which include cognitive performance. However, only a few prospective epidemiological studies exist and none among young adults. Here we address the gap in the literature by investigating whether the gut microbiome composition is prospectively linked to fluid intelligence among healthy young adults.

Methods

Forty individuals (65% females, 26 years) from the DOrtmund Nutritional and Anthropometric Longitudinally Designed (DONALD) study provided a fecal sample for gut microbiome composition and subsequently (average of 166 days) completed a cognitive functioning test using the Cattell’s Culture Fair Intelligence Test, revised German version (CFT 20-R). The assessment of the gut microbiome at the genera level was by 16S rRNA V3-V4 Illumina sequencing. The relative abundance of 158 genera was summarized into bacterial communities using a novel data-driven dimension reduction, amalgamation. The fluid intelligence score was regressed on the relative abundance of the bacterial communities and adjusted for selected covariates.

Results

The 158 genera were amalgamated into 12 amalgams (bacterial communities), which were composed of 18, 6, 10, 14, 8, 10, 16, 13, 12, 12, 3, and 11 genera. Only the 14-genera bacterial community, named the “Ruminococcaceae- and Coriobacteriaceae-dominant community” was positively associated with fluid intelligence score (β = 7.8; 95% CI: 0.62, 15.65, P = 0.04).

Conclusion

Among healthy young adults, the abundance of a gut bacterial community was associated with fluid intelligence score. This study suggests that cognitive performance may potentially benefit from gut microbiome-based intervention.

Learning sparse log-ratios for high-throughput sequencing data

MOTIVATION

The automatic discovery of sparse biomarkers that are associated with an outcome of interest is a central goal of bioinformatics. In the context of high-throughput sequencing (HTS) data, and compositional data (CoDa) more generally, an important class of biomarkers are the log-ratios between the input variables. However, identifying predictive log-ratio biomarkers from HTS data is a combinatorial optimization problem, which is computationally challenging. Existing methods are slow to run and scale poorly with the dimension of the input, which has limited their application to low- and moderate-dimensional metagenomic datasets.

RESULTS

Building on recent advances from the field of deep learning, we present CoDaCoRe, a novel learning algorithm that identifies sparse, interpretable and predictive log-ratio biomarkers. Our algorithm exploits a continuous relaxation to approximate the underlying combinatorial optimization problem. This relaxation can then be optimized efficiently using the modern ML toolbox, in particular, gradient descent. As a result, CoDaCoRe runs several orders of magnitude faster than competing methods, all while achieving state-of-the-art performance in terms of predictive accuracy and sparsity. We verify the outperformance of CoDaCoRe across a wide range of microbiome, metabolite and microRNA benchmark datasets, as well as a particularly high-dimensional dataset that is outright computationally intractable for existing sparse log-ratio selection methods.

AVAILABILITY AND IMPLEMENTATION

The CoDaCoRe package is available at https://github.com/egr95/R-codacore. Code and instructions for reproducing our results are available at https://github.com/cunningham-lab/codacore.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

DisBalance: a platform to automatically build balance-based disease prediction models and discover microbial biomarkers from microbiome data

How best to utilize the microbial taxonomic abundances in regard to the prediction and explanation of human diseases remains appealing and challenging, and the relative nature of microbiome data necessitates a proper feature selection method to resolve the compositional problem. In this study, we developed an all-in-one platform to address a series of issues in microbiome-based human disease prediction and taxonomic biomarkers discovery. We prioritize the interpretation, runtime and classification accuracy of the distal discriminative balances analysis (DBA-distal) method in selecting a set of distal discriminative balances, and develop DisBalance, a comprehensive platform, to integrate and streamline the workflows of disease model building, disease risk prediction and disease-related biomarker discovery for microbiome-based binary classifications. DisBalance allows the de novo model-building and disease risk prediction in a very fast and convenient way. To facilitate the model-driven and knowledge-driven discoveries, DisBalance dedicates multiple strategies for the mining of microbial biomarkers. The independent validation of the models constructed by the DisBalance pipeline is performed on seven microbiome datasets from the original article of DBA-distal. The implementation of the DisBalance platform is demonstrated by a complete analysis of a shotgun metagenomic dataset of Ulcerative Colitis (UC). As a free and open-source, DisBlance can be accessed at http://lab.malab.cn/soft/DisBalance. The source code and demo data for Disbalance are available at https://github.com/yangfenglong/DisBalance.