Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “whole-genome duplication”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources
PubMed2026

Revisiting the genome assembly of Lupinus species reveals differential diploidization after a shared whole-genome duplication.

Accurate genome assemblies are essential for comparative genomics, yet Hi-C-guided scaffolding can introduce structural errors that misrepresent chromosome architecture and bias evolutionary inferences. Here, we identified pervasive scaffolding errors-including artificial fusions, internal inversions, and incomplete contig mounting-in 2 previously published Lupinus genomes (L. cosentinii and L. digitatus) using a segmentation method based on long terminal repeat (LTR) retrotransposon density. We reassembled both genomes, producing chromosome-level references of 472.7 Mb (16 chromosomes) and 427.2 Mb (21 chromosomes), with BUSCO completeness >98.5%. Synteny validation and reapplication of LTR profiling confirmed that all prior errors were resolved. Using these corrected genomes together with 4 additional Lupinus species and 2 outgroup legumes, we investigated postpolyploid evolution. Synonymous substitution rate (Ks) analysis revealed a genus-specific whole-genome duplication (WGD) event (Ks = 0.17) shared by all 6 Lupinus species. The proportion of WGD-derived genes varied markedly, from 60% in L. digitatus to only 36% in L. mutabilis, indicating differential diploidization. While all species retained a core set of WGD duplicates enriched in cytoskeleton organization, ion transport, and defense responses, each exhibited lineage-specific functional trajectories: cell wall modification in L. cosentinii and L. digitatus, nitrogen metabolism in L. albus and L. angustifolius, flower development in L. luteus, and stress/lipid metabolism in L. mutabilis. Our corrected assemblies provide optimal references for Lupinus comparative genomics, and our findings demonstrate that a shared WGD event can lead to both conserved and highly divergent postpolyploid fates, likely underpinning adaptive diversification within the genus.

Lupinus↗
PubMed2026

Haplotype-resolved genome of Forsythia suspensa reveals the reticulate evolution in Oleaceae and a novel gene cluster regulating stamen development.

The olive family (Oleaceae) comprises numerous species of economic, horticultural, and medicinal importance. Despite its significance, the evolutionary history of this complex family remains enigmatic. Here, we generated a high-quality haplotype-resolved genome of Forsythia suspensa, a distylous species that occupies a key phylogenetic position in Oleaceae. The 2 haplotypes exhibit significant allelic divergence with potential allele-specific regulation. We reconstructed the polyploidization history of Oleaceae by confirming and precisely dating a shared whole-genome triplication and an independent whole-genome duplication event. We revealed a complex reticulate evolution that gave rise to the tribe Oleeae: an initial hybridization between Forsythieae (♂) and Jasmineae (♀), a subsequent backcrossing event, and a final whole-genome duplication. We identified a novel tandemly duplicated pectin methylesterase inhibitor gene cluster that regulates filament length and pollen size via restricting cell elongation in the long-styled morph. Dosage augmentation via stepwise cluster formation (0.99 to 3.83 Mya) may contribute to maintaining stamen traits of the long-styled morph. These FsPMEIs are co-expressed with many cell wall-related genes, suggesting a functional link in cell wall modification. Our study reveals the reticulate evolution in Oleaceae and a novel gene cluster controlling stamen development in F. suspensa and provides valuable haplotype-resolved genomic resources for heterostylous species, offering novel framework and molecular pathways to understand plant adaptive evolution.

Forsythia↗
PubMed2026

Genome-wide cyclin gene evolution in Arabidopsis and Brassica reveals polyploidization-driven duplication and flowering-time associations.

Cyclin genes are plant cell cycle regulators that play essential roles in growth, development, and reproduction. However, the evolutionary dynamics and genomic organization of cyclin genes across the Brassicaceae family remain poorly understood, particularly in the context of allotetraploid genome evolution. Here, we investigated the diversity, expansion mechanisms, and potential functional diversification of cyclin genes across ten Brassicaceae genomes, including four Arabidopsis and six Brassica species. A total of 1087 cyclin genes representing 23 cyclin types were identified. Comparative genomic analyses revealed that cyclin gene expansion was strongly influenced by polyploidization in Brassica species, with 1845 duplication events involving 1063 genes. Whole-genome duplication was the predominant mechanism driving expansion, while both inter- and intra-genomic duplications contributed to gene retention in tetraploid Brassica species, with the highest duplication frequency observed in Brassica juncea. Across genomes, 120 physical gene clusters were identified, including homogeneous and heterogeneous types. Ortholog analysis between progenitor and allotetraploid species identified 852 orthologous pairs involving 366 genes, indicating extensive conservation following allotetraploid formation. Phylogenetic analysis resolved cyclins into three major clades, while expression-based clustering in Brassica napus grouped genes into four major clusters, suggesting functional diversification. Integration of pan-genomic and flowering-time QTL analyses further identified two cyclin genes, Bna21cycA2 and Bna113cycD4, which contain amino acid polymorphisms and represent putative candidate variations potentially associated with flowering-time variation across multiple genomes. These findings provide new insights into the evolutionary expansion, retention, and potential functional divergence of cyclin genes in Brassicaceae and highlight candidate loci for future functional studies and crop improvement.

Evolution, Molecular↗
PubMed2026

Ancient polyploidization waves as evolutionary shields for angiosperms.

Chen et al. identified 132 whole-genome duplications (WGDs) clustered around environmental crises. We highlight how, over longer evolutionary timescales, ancient WGDs convergently retained MADS-box, MYB, WRKY and HSF transcription factors, building stress-adaptation networks. These insights guide climate-resilient crop improvement through comparative genomics and CRISPR engineering.

MADS-box↗
PubMed2026

Algae-to-host horizontal gene transfer in Paramecium bursaria is associated with host adaptation during endosymbiosis.

Paramecium bursaria maintains a stable endosymbiosis with green algae, yet the evolutionary consequences of this association remain unclear. Here, we screened the host genome for algal-derived horizontally transferred genes (HTGs) using a lineage-aware workflow designed to detect horizontal gene transfer (HGT) between two defined lineages. We identified 16 candidate HTGs, including four putative newly transferred genes and 12 homologous transferred genes, most of which were functionally associated with redox homeostasis and metabolism. Five HTGs showed symbiosis-dependent expression. RNAi knockdown of GH32s and SATs reduced host proliferation, total cell area, and motility, while GH32s knockdown also reduced endosymbiont load. Duplication patterns suggest that most transfers may have occurred after the P. bursaria lineage diverged from the sampled Paramecium species but before its lineage-specific whole-genome duplication (WGD). The HTGs also showed host-associated shifts in GC content and gene length, while representative HTGs retained conserved domains and functional motifs. Together, our results support algae-to-host HGT in P. bursaria and suggest that some transferred genes may contribute to metabolic integration during endosymbiosis.

Gene Transfer, Horizontal↗
PubMed2026

Genome-wide identification, structural characterization, and evolutionary analysis of growth-related gene families in African catfish (Clarias gariepinus).

The somatotropic axis encompassing growth hormone (GH), insulin-like growth factor (IGF), myostatin (MSTN), and prolactin (PRL) signalling cascades is the master regulator of somatic growth, metabolism, and development in vertebrates. African catfish (Clarias gariepinus), a commercially pivotal aquaculture species, now possesses a chromosome-level reference genome (CGAR_prim_01v2); however, a systematic, genome-wide characterization spanning all five interconnected growth-related gene families has not previously been undertaken in this species. Here, we identified and characterized 15 growth-related genes spanning gh1, ghra, ghrb, Igf1, Igf2a, Igf2b, igf1ra, Igf1rb, Igf2r, Mstna, Mstnb, prl, prlra, prlrb, and smtlb distributed across 13 chromosomes. Complete one-to-one orthology with zebrafish confirmed strong dosage-balance conservation across >120 million years of teleost divergence. Physicochemical analysis resolved a clear biochemical dichotomy between compact, basic secreted ligands (19.88-45.81 kDa; pI up to 10.02) and large, acidic, heavily glycosylated membrane receptors (56.82-270.80 kDa; pI 4.85-5.97). Phylogenetic analysis confirmed 3R whole-genome duplication origins for all paralog pairs, while synteny analysis revealed a disruption of the ancestral gh1-prl chromosomal block in C. gariepinus, a finding that warrants further comparative and functional investigation. This genomic atlas provides the sequence and structural information including exon-intron boundaries, domain architecture, and chromosomal coordinates needed as a prerequisite for future marker-assisted selection and CRISPR-based myostatin-editing efforts in African catfish aquaculture, though translation into applied breeding outcomes will require subsequent functional and expression studies.

Animals↗
PubMed2026

Reconstruction of ancestral plant genomes for inter-crop translational research.

We present Ancestral Genome Reconstruction (AGR), an exploratory framework for the automated inference of "paleogenomes" from large-scale comparative datasets. By analyzing 84 extant angiosperm species, we reconstructed 10 key ancestral angiosperm genomes millions of years old. These reconstructed ancestors were instrumental in (1) estimating when angiosperms emerged, when major botanical families originated, and when shared ancestral whole-genome duplication events occurred; and (2) tracing the evolutionary trajectories of ancestral chromosomes and genes, especially those that may have driven the emergence of key life-history traits (e.g., woody vs. herbaceous, aquatic vs. terrestrial, C3 vs. C4, and symbiotic root-nodulating vs. non-nodulating species). We demonstrated that these paleogenomes serve as tractable backbones for inter-crop translational research. Through an open-access web tool, OrthoViewer, we identified orthologs that have retained the same ancestral genomic context, favoring the identification of genes associated with "phenologs"- orthologous genes across species driving analogous phenotypes, traits, or processes-exemplified by FUWA for yield components, FLC for flowering time, and DDM1 for DNA methylation. Taken together, this study provides a testable paleogenomic workflow, opening novel avenues for integrating evolutionary genomics data into modern climate-smart crop breeding and supporting the agroecological transition.

Genome, Plant↗
PubMed2026

From neuropeptides to toxins: illuminating the origins of venom complexity in cone snails.

New genes and gene functions are key drivers of evolutionary innovation. Venomous animals, such as cone snails, provide striking examples of gene innovation, yet the mechanisms by which toxins arise remain poorly understood. Using the Conus textile genome, we uncover how neuropeptide genes were recruited into the venom and neofunctionalized as doppelgänger toxins. We identify over 20 independent recruitment events that evolved dynamically across the Conus lineage. Rather than arising from ohnologs of a whole-genome duplication event ∼200 mya, these toxins evolved through diverse mechanisms, including exon shuffling, alternative splicing, and ectopic recombination, often facilitated by lineage-specific transposable elements. Our findings reveal a dynamic interplay between genome architecture and molecular innovation, offering broad insight into the evolution of complex gene repertoires in venoms and beyond.

Animals↗
PubMed2026

Mining the sHSP20 (small heat-shock protein) gene family in finger millet (Eleusine coracana (L.) Gaertn.): structural, evolutionary and predicted abiotic-stress-responsive insights.

Small heat-shock proteins (sHSPs, the HSP20 family) are ATP-independent molecular chaperones that hold partially unfolded substrates and protect the proteome during heat and other abiotic stresses; every member is defined by a conserved &#x3b1;-crystallin domain (ACD). Finger millet (Eleusine coracana) is a climate-resilient, calcium-rich allotetraploid cereal of the semi-arid tropics whose HSP20 repertoire had not been catalogued. The present study is an entirely computational (in silico) analysis of the chromosome-scale reference genome of finger millet (NCBI GenBank assembly GCA_032690845.1, cultivar KNE 796-S). Mining the predicted proteome with the ACD profile (Pfam PF00011) and confirming every candidate by NCBI CD-search recovered 76 non-redundant ACD-bearing HSP20 genes (EcHSP20-1-EcHSP20-76). Based on phylogeny and TargetP-predicted localization, the members were classified into ten subfamilies: seven cytosolic/nuclear classes (C-I to C-VII, 60 members) together with chloroplastic (11), mitochondrial (3) and endoplasmic-reticulum (2) groups. The proteins ranged from 110 to 355 amino acids (12.1-39.2&#xa0;kDa) with theoretical pI of 4.85-9.69. The 76 loci were distributed over 14 of the 18 chromosomes and were conspicuously absent from chromosomes 8&#xa0;A, 8B, 9&#xa0;A and 9B, with pronounced clustering on chromosomes 1, 2, 3 and 6. Duplication analysis detected 149 paralogous pairs (49 homoeologous, 80 segmental/dispersed and 18 tandem); 147 of 148 pairs for which substitution rates could be calculated returned Ka/Ks&#x2009;<&#x2009;1 (mean 0.20), indicating strong purifying selection consistent with retention after whole-genome/allopolyploid duplication. Promoter analysis (PlantCARE) revealed enrichment of abscisic-acid-responsive (ABRE), MYB/MYC drought-related, STRE, DRE, low-temperature (LTR) and methyl-jasmonate/salicylic-acid elements, whereas canonical heat-shock elements (HSE) were not recovered. Expression profiling against a public drought transcriptome (SRP081350) showed that about half of the genes (39 of 76) are transcribed in leaf tissue, the expressed fraction being dominated by the cytosolic class C-I. This first finger-millet HSP20 catalogue provides a verified, reproducible framework and nominates computationally predicted candidate genes for future functional work on thermotolerance in cereals.

Allotetraploid↗
PubMed2025

Comparative Pharmacological Analysis of Mitotic Inhibitors Using Isogenic Ploidy Series of HAP1 Cells.

Drastic changes in chromosome number and cellular contents upon ploidy alterations profoundly affect the stability of mitotic regulation in different biological and pathological processes. Isogenic ploidy series of somatic cell lines are useful for studying the effects of ploidy differences on mitotic regulation at cellular and molecular levels. This chapter describes experimental procedures using isogenic human HAP1 cell lines that cover haploid, diploid, and tetraploid states. We first describe methods to establish and maintain these isogenic HAP1 ploidy series using a flow cytometer. We then describe a procedure of comparative pharmacological assay for analyzing ploidy-dependent changes in the functionality of the mitotic spindle components.

Humans↗
PubMed2026

Establishment of induced pluripotent stem cell line TRNDi045-A-38 carrying homozygous DOK7-related Congenital Myasthenia patient-mutation knock-in variant from parental KOLF2.1J.

DOK7-related Congenital Myasthenic Syndrome (CMS) is a rare genetic neuromuscular junction disorder. This is one of the most common of the recessive forms of CMS, often presenting with more static proximal weakness (hence also referred to as limb girdle CMS). Whole-genome sequencing of affected patients implicates frameshift duplication mutations in DOK7 as drivers of impaired neuromuscular-junction signaling. In this study, we generated a human induced pluripotent stem cell (hiPSC) line TRNDi045-A-38 from the KOLF2.1J reference line, engineered to carry homozygous DOK7 c.1124_1127dupTGCC mutation knock-in using CRISPR/Cas9. This iPSC line could be used for in vitro disease modeling to study disease pathophysiology and for therapeutic development.

Humans↗
PubMed2026

Mechanisms of cefiderocol resistance in carbapenem-resistant Acinetobacter baumannii: a Swiss 2023-2025 collection.

OBJECTIVES: The numbers of infections caused by carbapenem-resistant Acinetobacter baumannii (CRAB) are increasing globally and present a significant burden on healthcare systems. This study describes the CRAB isolates received at the Swiss National Reference Centre for Emerging Antibiotic Resistance (NARA) over a 3-year period, from January 2022 to December 2025, and aimed to characterize the prevalence and mechanisms of FDC resistance. METHODS: Two-hundred and thirty-four non-duplicate CRAB isolates were submitted to NARA over the study period from hospitals and laboratories across Switzerland. Susceptibility testing was performed by disk diffusion and broth microdilution, according to EUCAST methodology. Whole-genome sequencing was performed on 11 isolates. ADC alleles were cloned into vector pVRL1 and transformed into Escherichia coli Top10. RESULTS: All isolates exhibited resistance to the carbapenems, and most were resistant to cephalosporins. Most isolates harboured an acquired class D carbapenemase, most frequently OXA-23 (181/234; 77.4%). One quarter of isolates were resistant to cefiderocol (FDC), exhibiting MICs ranging from 4->32 mg/L. Whole genome sequencing analyses, performed on 11 FDC-resistant isolates, identified that FDC resistance was due a combination of mechanisms including NDM and PER-production, mutations within the iron transporters, piuA and pirA, and the overexpression of ADC variants. CONCLUSIONS: This study showed that OXA-23 was the dominant mechanism of carbapenem-resistance in CRAB in Switzerland. Almost one quarter of CRAB isolates were resistant to "last resort" antimicrobial, FDC. The mechanisms of FDC resistance identified in this study emphasise that resistance to this antimicrobial is often complex and multifactorial, requiring high-resolution methods, including WGS, to identify.

Acinetobacter baumannii↗
PubMed2026

Efficacy of the NMIC-150 system in identifying extended-spectrum beta-lactamases in clinical isolates.

Extended-spectrum beta-lactamases (ESBLs) are significant contributors to the growing global crisis of antimicrobial resistance. This study evaluated the performance of the NMIC-150 System for susceptibility testing of third-generation cephalosporins (3GCs) and assessed whether ceftazidime-avibactam and aztreonam-avibactam could identify ESBL-producing carbapenem-resistant Enterobacterales (CREs). A total of 278 non-duplicate clinical isolates (Klebsiella pneumoniae, E. coli, and Proteus mirabilis) were analyzed. Antimicrobial susceptibility was determined using reference broth microdilution (BMD) and the NMIC-150 System. ESBL production was defined as an &#x2265;eight-fold reduction in the minimum inhibitory concentration (MIC) of 3GCs in the presence of clavulanic acid, according to CLSI criteria. Whole-genome sequencing was performed to characterize ESBL and carbapenemase genes among 3GC-resistant isolates. A Random Forest model was used to predict ESBL-producing isolates based on MIC values. The NMIC-150 System demonstrated over 90% categorical and essential agreement with BMD for ceftazidime and ceftriaxone, along with robust predictive performance via Random Forest analysis. These findings suggest that the NMIC-150 System is a reliable platform for 3GC susceptibility testing and that an &#x2265;eight-fold MIC reduction with ceftazidime-avibactam or aztreonam-avibactam may serve as a phenotypic indicator of ESBL production in CRE isolates. In conclusion, the NMIC-150 System shows potential for routine antimicrobial resistance surveillance and may facilitate the rapid identification of ESBL-producing CREs in clinical settings.

Microbial Sensitivity Tests↗
PubMed2026

Global Diffusion of IncC Plasmid Harboring blaNDM-1in the High-Risk Escherichia coli ST131 Clone.

AIMS: The global expansion of quinolone-resistant Escherichia coli (QR-EC) is increasingly associated with &#x3b2;-lactam resistance and mobile genetic elements that facilitate resistance dissemination. This study investigated the molecular mechanisms underlying fluoroquinolone and &#x3b2;-lactam resistance in clinical QR-EC isolates and explored the plasmid type associated. METHODS AND RESULTS: A total of 123 non-duplicate QR-EC clinical isolates responsible mainly for gastrointestinal colonization were collected between 2019 and 2021. Plasmid-mediated quinolone resistance (PMQR), extended-spectrum &#x3b2;-lactamase (ESBL), and carbapenemase genes were screened by PCR. Mutations in the quinolone resistance-determining regions (QRDR) of gyrA and parC were analyzed using sequencing and mismatch amplification mutation assay PCR (MAMA-PCR). Selected isolates underwent multilocus sequence typing (MLST). Whole-genome sequencing (WGS) of a representative extensively drug-resistant strain carrying multiple quinolone resistance determinants, ESBL genes, and carbapenemase genes, was performed. PMQR genes were prevalent among QR-EC, dominated by aac(6')-Ib-cr (60.9% of isolates). ESBL genes were identified in 93.5% of isolates, predominantly blaCTX-M (95.7%). Among ertapenem-resistant isolates (QCR-EC) (n=18), blaNDM-1 and blaOXA-48 were detected in 13 and 11 isolates, respectively. QRDR mutations were highly frequent, particularly gyrA83 (98.4%) and parC80 (30.9%). Major QCR-EC genotypes belonged to sequence types ST167 (n=2), ST1196, ST469, and ST410. High-risk E. coli ST131 clone harboring IncC plasmid encoding blaNDM-1 was described for the first time in Africa, following its emergence, in two continents, Asia and America. Despite the very rare description of these strains worldwide, their description in three continents sign their global diffusion. CONCLUSIONS: This finding highlights the ongoing spread of carbapenem resistance and underscores the urgent need for strengthened genomic surveillance.

Escherichia coli↗
PubMed2026

Genetic Screening of Colombian Patients With Early-Onset Parkinson Disease.

BACKGROUND AND OBJECTIVES: Early-onset Parkinson disease (EOPD), defined as symptom onset before 50 years of age, accounts for approximately 10% of patients and is suggested to have a greater genetic component than typical late-onset forms of the disease. Recessive variants in PRKN, PINK1, and DJ-1, are the most common genetic cause of EOPD, however, most studies are in patients of white ancestry. This study aims to analyze genetic variants in PRKN, PINK1, and DJ-1 in Colombian patients to help address the gap in EOPD genetic research of South American populations. METHODS: We analyzed 43 unrelated patients with EOPD using Sanger sequencing for the PRKN, PINK1, and DJ-1 genes and employed multiplex ligation-dependent probe amplification to detect copy number variants. Additionally, long-read whole-genome sequencing was conducted on 3 unresolved patients with age at onset before 30 years of age (long-read sequencing [LRS] patient A-C). RESULTS: We identified known pathogenic single-nucleotide variants and copy number variants in the PRKN gene accounting for 2 patients' disease (4.6% of patients). We observed 2 pathogenic variants in PRKN (c.155delA; p.N52Mfs*29 and c.1083+1G>A) in patient 1, who reported an age at onset of 16 years. We further detected a homozygous duplication of PRKN exons 5-6 in an additional patient, age at onset of 18 years. DISCUSSION: Our study helps characterize genetic contributors to EOPD in Colombian patients, demonstrating genetic forms (PRKN, PINK1, and DJ-1) are rare. Our results highlight a need to include diverse populations in research to improve genetic understanding of disease.

Journal Article↗
PubMed2026

Genetic mutations driving ciprofloxacin resistance in laboratory-evolved Salmonella Typhimurium.

Ciprofloxacin resistance in Salmonella Typhimurium is a significant public health concern, and the mechanisms by which the resistance evolves are poorly defined. Here, by serial passaging under antibiotic selection, we isolated ciprofloxacin-resistant S. Typhimurium mutants and subjected them to whole-genome sequencing to reveal the major mutations associated with resistance. The Low CipR mutant acquired four chromosomal mutations in ramR, icdA, lipB, and gyrA, and the High CipR mutant gained additional mutations in gyrB, yaiC, and corA. Functional characterization determined that mutations in ramR resulted in efflux pump upregulation, while disruptions in the TCA cycle caused by mutations in icdA and lipB led to metabolic alterations. These changes indirectly enhanced resistance by increasing the expression of the global regulator MarA and reducing OmpF-dependent membrane permeability. Despite the observation of the G105A substitution in GyrA, enzymatic assays confirmed the failure to support resistance to ciprofloxacin, possibly because the structural alteration remained minimal. GyrB488-489dup was associated with maintained supercoiling under ciprofloxacin and enhanced fluoroquinolone resistance, suggesting a major role in resistance evolution. Other mutations in yaiC impaired biofilm and, in corA, intracellular accumulation of magnesium, possibly stabilizing the bacterial cell envelope under antibiotic pressure. The findings provide novel explanations for the multifaceted mechanisms leading to ciprofloxacin resistance in Salmonella and suggest targets to combat antimicrobial resistance.IMPORTANCEAntibiotic resistance in Salmonella Typhimurium is an increasing public health concern, yet the genetic changes that allow bacteria to become resistant are not fully understood. In this study, we evolved ciprofloxacin-resistant Salmonella in the laboratory and identified the mutations that arise during resistance development. We found that resistance does not result from a single change but from multiple adaptations affecting drug efflux, metabolism, and the antibiotic target. Some mutations increased the activity of pumps that remove antibiotics from the cell, while others altered bacterial metabolism and reduced membrane permeability, making it harder for the drug to enter. A duplication in the DNA gyrase subunit GyrB played a particularly important role in maintaining DNA function under antibiotic stress. Together, these results reveal how diverse genetic changes cooperate to generate ciprofloxacin resistance and provide insights that may help guide strategies to combat drug-resistant Salmonella infections.

DNA gyrase↗
PubMed2026

A genome-wide coverage-based pipeline for the identification of host-derived candidate DNA biomarkers from cell-free blood.

We have created a new data-analysis pipeline for the discovery of host-specific candidate DNA biomarkers derived from sequencing data of cell-free blood. Unlike approaches that rely on specific molecular or genetic signatures, our method leverages the coverage distribution of cell-free DNA sequences mapped to a reference genome, applying statistical analyses to identify informative short genomic regions for biomarker discovery. The pipeline is applicable to diverse diseases and can be used to analyze cell-free DNA sequences from plasma or serum to identify candidate biomarkers that are characteristic of disease states in mammals. Core functionalities were developed in Java and integrated with open-source software tools for the preprocessing of raw sequencing data, complemented by Python scripts for the machine-learning analysis and statistical validation. The pipeline is designed for HPC use and users can access the pipeline through a Galaxy workflow, which offers a user-friendly web interface for input selection prior to execution and analysis progress monitoring. Performance tests, carried out using duplicate sets of COVID-19 samples and controls, showed linear scalability of execution time with an increasing dataset size, as well as a substantial reduction in execution time through parallelized computation, whereby each HPC node is used to process the data of one chromosome. Further statistical tests confirmed the quality of the pipeline's results by showing that the set of identified candidate biomarkers remained stable across varying dataset sizes.

Biomarkers↗