PubMed Health⌕ Search

Biomedical subjects

Manuel Rueda

Publications and source records attributed to Manuel Rueda.

9 recordsLinked to original sources

CBIcall: a configuration-driven framework for variant calling in large sequencing cohorts.

MOTIVATION: Variant calling for next-generation sequencing (NGS) data relies on a diverse ecosystem of tools and workflows. Large-scale collaborative studies increasingly adopt federated analysis, where each institution processes sensitive data locally using standardized pipelines. Deploying identical pipelines across multiple centers remains challenging because heterogeneous software environments and computing policies can cause workflow divergence and inconsistent results. RESULTS: We developed CBIcall, a workflow backend-flexible, configuration-driven framework that runs standardized variant-calling pipelines from raw FASTQ files to analysis-ready VCFs. Users define each analysis in a single YAML parameters file, which CBIcall resolves against a controlled workflow registry and resource catalog. The execution driver validates parameters and checks compatibility among pipelines, analysis modes, workflow backends, genome builds, tool versions, and resource bundles. CBIcall supports reproducibility auditing by comparing executions using recorded provenance and output fingerprints. CBIcall dispatches validated workflows natively through Bash, Cromwell, Nextflow and Snakemake backends and provides production-ready pipelines for germline WES, WGS (single-sample or cohort joint genotyping following GATK Best Practices), and mitochondrial DNA analysis. We evaluated analytical performance using public benchmark datasets and validated reproducibility across four computing environments. We further deployed CBIcall in the EU HEREDITARY project, where it processed 1102 samples with both WES and mtDNA pipelines on an institutional HPC system, supporting its suitability for reproducible cohort-scale genomic analyses. AVAILABILITY AND IMPLEMENTATION: CBIcall is open source (GPLv3) and distributed with ready-to-run pipelines; full dependency and installation documentation is available at https://github.com/CNAG-Biomedical-Informatics/cbicall.

Journal Article↗

ClarID: A Human-Readable and Compact Identifier Specification for Biomedical Metadata Integration.

BACKGROUND: In biomedical research, subjects and biospecimens are commonly tracked using simple IDs or UUIDs, which guarantee uniqueness but convey no embedded semantic information. Contextual metadata (such as tissue type, diagnosis, or assay) is often stored separately, making integration, cohort selection, and downstream analysis cumbersome. While structured barcoding systems exist in large consortia (e.g., TCGA, GTEx) or domain-specific contexts (e.g., SPREC, GOLD), no unified, extensible framework currently spans both subjects and biosamples in a human- and machine-readable way. METHODS: We developed ClarID, a domain-agnostic specification that supports two identifier formats: (i) a human-readable form (e.g., 'CNAG_Test-HomSap-00001-LIV-TUM-RNA-C22.0-TRT-P1W' that encodes key metadata such as project, species, subject_id, tissue, assay, disease, timepoint and duration (from that event); and (ii) a compact version named 'stub' (e.g., 'CT01001LTR0N401T1W') optimized for filenames, pipelines, and labeling.ClarID is implemented through an open-source command-line tool, ClarID-Tools, which processes tabular metadata files (CSV/TSV) and uses a YAML-based codebook to generate, decode, and validate identifiers, as well as to create and read QR codes. The tool supports bulk and single-sample processing and allows easy integration with institutional workflows. RESULTS: To demonstrate ClarID's utility, we applied it to datasets from the Genomic Data Commons (GDC), generating interpretable identifiers for more than 113,000 clinical records (subjects) and 4,255 biospecimen records. All materials, including pre-processing scripts, input and encoded data, are publicly available and fully reproducible via the accompanying GitHub repository and Google Colab. CONCLUSIONS: ClarID fills a critical gap between opaque accession numbers and rich metadata schemas by embedding key context directly into structured identifiers. It enhances traceability, facilitates downstream analysis, and remains adaptable to project-specific needs through a configurable codebook. The accompanying ClarID-Tools software is freely available, together with full documentation and reproducible pipelines, at https://github.com/CNAG-Biomedical-Informatics/clarid-tools.

Biosample identifiers↗

A consensus view of protein dynamics.

The dynamics of proteins in aqueous solution has been investigated through a massive approach based on "state of the art" molecular dynamics simulations performed for all protein metafolds using the four most popular force fields (OPLS, CHARMM, AMBER, and GROMOS). A detailed analysis of the massive database of trajectories (>1.5 terabytes of data obtained using approximately 50 years of CPU) allowed us to obtain a robust-consensus picture of protein dynamics in aqueous solution.

Computational Biology↗

G-quadruplexes can maintain their structure in the gas phase.

Several very extended (0.5-1 micros) molecular dynamics (MD) simulations of parallel and antiparallel G-quadruplex DNA strongly suggest that in the presence of suitable cations the quadruplex not only remains stable in the gas phase, but also displays a structure that closely resembles that found in extended (25-ns long) trajectories in aqueous solution. In the absence of the crucial cations, the trajectories become unstable and in general the quadruplex structure is lost. To our knowledge, this is the first physiologically relevant structure of DNA for which very large MD simulations suggest that the structure in water and in the gas phase are indistinguishable.

Computer Simulation↗

Data mining of molecular dynamics trajectories of nucleic acids.

Analysis, storage, and transfer of molecular dynamic trajectories are becoming the bottleneck of computer simulations. In this paper we discuss different approaches for data mining and data processing of huge trajectory files generated from molecular dynamic simulations of nucleic acids.

Computer Simulation↗

Nature of minor-groove binders-DNA complexes in the gas phase.

The structure of noncovalent complexes of DNA duplex with minor groove binders (mG-binders) has been analyzed by state of the art molecular dynamics (MD) simulations. More than 3.3 micros of MD trajectories (including 4 x 0.5 micros trajectories) were collected for the Dickerson's dodecamer bound to DAPI, Hoechst 33258, and Netropsin. Comparison of these trajectories with control simulations in water allowed us to determine that the extreme dehydration and partial neutralization occurring during electrospray experiments does not produce the disruption of the DNA:mG-binder complexes or the dissociation of the two strands of the duplex. Irrespective of the drug and the simulation conditions the mG-binders remains bound to the DNA near the preferential binding position in aqueous conditions. Large distortions appear in the two DNA strands, which maintain however a memory of the original DNA duplex structure in water, and a general helical-like conformation.

Binding Sites↗

Exploring the counterion atmosphere around DNA: what can be learned from molecular dynamics simulations?

The counterion distribution around a DNA dodecamer (5'-CGCGAATTCGCG-3') is analyzed using both standard and novel techniques based on state of the art molecular dynamics simulations. Specifically, we have explored the population of Na(+) in the minor groove of DNA duplex, and whether or not a string of Na(+) can replace the spine of hydration in the narrow AATT minor groove. The results suggest that the insertion of Na(+) in the minor groove is a very rare event, but that when once the ion finds specific sites deep inside the groove it can reside there for very long periods of time. According to our simulation the presence of Na(+) inside the groove does not have a dramatic influence in the structure or dynamics of the duplex DNA. The ability of current MD simulations to obtain equilibrated pictures of the counterion atmosphere around DNA is critically discussed.

Binding Sites↗

The structure and dynamics of DNA in the gas phase.

The impact of the transfer from water to the vacuum in the duplex DNA has been explored by using long molecular dynamic simulations. In opposition to chemical intuition, it is found that vaporization of DNA, even at high temperatures, does not lead to a total disruption of the double helix. Rather, the DNA duplex preserves gross structural, energetic, and dynamic features of the conformation of the double helix in aqueous solution. Thus, the two strands remain bound, the global structure has a slight helicity, and the total number of DNA-DNA interactions is not dramatically different from that found in solution. The results provide detailed structural and dynamic information useful to complement current mass spectroscopy studies of DNA structures.

Computer Simulation↗