2025 Hyper Recent •CC0 1.0 Universal

This work is dedicated to the public domain. No rights reserved.

Access Preprint From Server
February 1st, 2025
Version: 1
Inria
bioinformatics
biorxiv

Predicting coarse-grained representations of biogeochemical cycles from metabarcoding data

Belcour, A.Open in Google Scholar•Megy, L.Open in Google Scholar•Stephant, S.Open in Google Scholar•Michel, C.Open in Google Scholar•Rad, S.Open in Google Scholar•Bombach, P.Open in Google Scholar•Dopffel, N.Open in Google Scholar•de Jong, H.Open in Google Scholar•Ropers, D.Open in Google Scholar

Motivation: Taxonomic analysis of environmental microbial communities is now routinely performed thanks to advances in DNA sequencing. Determining the role of these communities in global biogeochemical cycles requires the identification of their metabolic functions, such as hydrogen oxidation, sulfur reduction, and carbon fixation. These functions can be directly inferred from metagenomics data, but in many environmental applications metabarcoding is still the method of choice. The reconstruction of metabolic functions from metabarcoding data and their integration into coarse-grained representations of geobiochemical cycles remains a difficult bioinformatics problem today. Results: We developed a pipeline, called Tabigecy, which exploits taxonomic affiliations to predict metabolic functions constituting biogeochemical cycles. In a first step, Tabigecy uses the tool EsMeCaTa to predict consensus proteomes from input affiliations. To optimise this process, we generated a precomputed database containing information about 2,404 taxa from UniProt. The consensus proteomes are searched using bigecyhmm, a newly developed Python package relying on Hidden Markov Models to identify key enzymes involved in metabolic function of biogeochemical cycles. The metabolic functions are then projected on coarse-grained representation of the cycles. We applied Tabigecy to two salt cavern datasets and validated its predictions with microbial activity and hydrochemistry measurements performed on the samples. The results highlight the utility of the approach to investigate the impact of microbial communities on geobiochemical processes. Availability: The Tabigecy pipeline is available at https://github.com/ArnaudBelcour/tabigecy. The Python package bigecyhmm and the precomputed EsMeCaTa database are also separately available at https://github.com/ArnaudBelcour/bigecyhmm and https://doi.org/10.5281/zenodo.13354073, respectively.

Similar Papers

biorxiv
Sun Feb 02 2025
GreedyMini: Generating low-density DNA minimizers
Minimizers is the most popular k-mer selection scheme in algorithms and data structures analyzing high-throughput sequencing (HTS) data. In a minimizers scheme, the smallest k-mer by some predefined order is selected as the representative of a sequence window containing w consecutive k-mers, which results in overlapping windows often selecting the same k-mer. Minimizers that achieve the lowest fre...
Golan, S.
•
Tziony, I.
•
Kraus, M.
•
Orenstein, Y.
•
Shur, A.
biorxiv
Sat Feb 01 2025
VIVIDHA: Variant Prediction and Visualization Interface for Dynamic High-throughput Analysis
Large scale genome sequencing projects have produced huge datasets that pose challenges of high processing times especially for variant calling, a significant downstream analysis step. Efficient utilization of computational resources for accurate variant prediction in a timely manner is possible using Hadoop MapReduce framework. We have developed VIVIDHA (Variant Prediction and Visualization Inter...
Barasiya, T.
•
Bharti, N.
•
Gadhari, R.
•
Bayaskar, A.
...•
Joshi, R.
biorxiv
Sat Feb 01 2025
Potentials and limitations in the application of Convolutional Neural Networks for mosquito species identification using wing images
1. This study addresses the pressing global health burden of mosquito-borne diseases by investigating the application of Convolutional Neural Networks (CNNs) for mosquito species identification using wing images. Conventional identification methods are hampered by the need for significant expertise and resources, while CNNs offer a promising alternative. Our research aimed to develop a reliable an...
Nolte, K.
•
Baumbach, J.
•
Lins, C.
•
Lohmann, J. J. G.
...•
Luehken, R.
biorxiv
Sat Feb 01 2025
AnchoRNA: Full virus genome alignments through conserved anchor regions
Motivation: Multiple sequence alignment of full viral genomes can be challenging due to factors such as long sequences, large insertions/deletions (spanning several 100 nucleotides), high number of sequences, sequence divergence, and high computational complexity, as e.g., in the context of secondary structure prediction. Standard alignment methods often face these issues, especially when processi...
Eulenfeld, T.
•
Triebel, S.
•
Marz, M.
biorxiv
Sat Feb 01 2025
What is a differentially expressed gene?
The concept of \'Differentially Expressed Genes\' (DEGs) is central to RNA-Seq studies, yet their identification suffers from reproducibility issues. This is largely a consequence of the inherent biological and technical variation that cannot be captured with small numbers of replicates. When thresholds for p-values and log fold changes are introduced, this variability can propagate an incomplete ...
Hoerbst, F.
•
Sidhu, G. S.
•
Tomkins, M.
•
Morris, R. J.
biorxiv
Sat Feb 01 2025
SCGclust: Single Cell Graph clustering using graphautoencoders integrating SNVs and CNAs
Intra-tumor heterogeneity (ITH) is a compounding factor for cancer prognosis and treatment. Single-cell DNA sequencing (scDNA-seq) provides cellular resolution of the variations in a cell and has been widely used to study cancer progression and responses to drug and treatment. While the low coverage scDNA-seq technologies typically provides a large number of cells, accurate cell clustering is esse...
Potu, T.
•
Hu, Y.
•
Khan, R.
•
Dharani, S.
...•
Mallory, X.
biorxiv
Sat Feb 01 2025
Benchmarking gene embeddings from sequence, expression, network, and text models for functional prediction tasks
Gene embeddings have emerged as transformative tools in computational biology, enabling the efficient translation of complex biological datasets into compact vector representations. This study presents a comprehensive benchmark by evaluating 38 classic and state-of-the-art gene embedding methods across a spectrum of functional prediction tasks. These embeddings, derived from data sources such as a...
Zhong, J.
•
Li, L.
•
Dannenfelser, R.
•
Yao, V.
biorxiv
Sat Feb 01 2025
Beyond contacts: the important role of the support region to distinguish stable and transient protein interfaces
Protein-protein interactions (PPIs) are fundamental to virtually all cellular processes; however, elucidating the principles that govern protein association into complexes remains a significant challenge. In this study, we present a comparative analysis of stable and transient protein interfaces, offering a detailed perspective on the interactions that form them. Moving beyond the traditional focu...
Miclot, T.
•
Timr, S.
biorxiv
Sat Feb 01 2025
Evaluation of sequencing reads at scale using rdeval
Motivation: Large sequencing data sets are produced and deposited into public archives at unprecedented rates. The availability of tools that can reliably and efficiently generate and store sequencing read summary statistics has become critical. Results: As part of the effort by the Vertebrate Genomes Project (VGP) to generate high-quality reference genomes at scale, we sought to address the commu...
Formenti, G.
•
Koo, B.
•
Sollitto, M.
•
Balacco, J.
...•
Jarvis, E. D.