Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Genetic testing algorithm”

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 109 records · Page 6Linked to original sources
PubMed2025

Maternal-Foetal HLA-DQB1 Incompatibility Is Associated With Pregnancy-Induced Hypertensive Disorders in a Genetically Isolated Population.

In pregnancy, semi-allogenic foetal trophoblasts express a specific HLA profile mediating maternal leukocyte contact, crucial for placentation. Paradoxically, maternal immunomodulation requires foetal antigen recognition, especially involving certain HLA molecules. Pre-eclampsia, a severe hypertensive complication, has been linked to antigenic similarity. Previously, we showed no selection for HLA (in)compatibility in uncomplicated naturally conceived pregnancies. However, pre-eclamptic pregnancies were associated with increased total maternal-foetal HLA and HLA-C matching. These associations suggest a role for HLA mismatches in immune regulation leading to an uncomplicated pregnancy. To better understand HLA homozygosity in human reproduction, we aimed to determine if there is a preferential selection for HLA compatibility in a genetically isolated population, and its relation to hypertensive complications. A nested case-control study, comprising 125 uncomplicated pregnancies and 50 with hypertensive complications (29 with pregnancy-induced hypertension, 21 with pre-eclampsia) was conducted in a genetically isolated Dutch population (FROH 1.3-3.1). Maternal and foetal HLA-A, -B, -C, -DRB1, -DQA1, -DQB1 and maternal killer-cell immunoglobulin-like receptor (KIR) genotyping were performed. Maternal-foetal HLA (mis)match counts were compared to expected values from randomisation of paternal HLA haplotypes over maternal haplotypes of the foetuses. Mismatched CD4+ T cell epitopes presented by maternal HLA class II were predicted using the PIRCHE-II algorithm. In uncomplicated pregnancies, no difference was found between observed and expected maternal-foetal HLA (mis)matches. However, pregnancies with hypertensive complications showed significantly higher observed HLA-DQB1 mismatches, reflected in PIRCHE-II scores. No significant differences were found in KIR/HLA-C frequencies. Interpretation is limited by the small sample size and the grouping of distinct hypertensive disorders. Nonetheless, maternal-foetal HLA-DQB1 mismatch seems to play a role in the aetiology of hypertensive complications during pregnancy in this population.

Humans↗
PubMed1987

An exact test for Hardy-Weinberg and multiple alleles.

Algorithms for generating the exact distribution of a finite sample drawn from a population in Hardy-Weinberg equilibrium are given for multiple alleles. The finite sampling distribution is derived analogously to Fisher's 2 X 2 exact distribution and is equivalent to Levene's conditional finite sampling distribution for Hardy-Weinberg populations. The algorithms presented are fast computationally and allow for quick alternatives to standard methods requiring corrections and approximations. Computation time is on the order of a few seconds for three-allele examples and up to 2 minutes for four-allele examples on an IBM 3081 machine.

Algorithms↗
PubMed1987

Tertiary templates for proteins. Use of packing criteria in the enumeration of allowed sequences for different structural classes.

We assume that each class of protein has a core structure that is defined by internal residues, and that the external, solvent-contacting residues contribute to the stability of the structure, are of primary importance to function, but do not determine the architecture of the core portions of the polypeptide chain. An algorithm has been developed to supply a list of permitted sequences of internal residues compatible with a known core structure. This list is referred to as the tertiary template for that structure. In general the positions in the template are not sequentially adjacent and are distributed throughout the polypeptide chain. The template is derived using the fixed positions for the main-chain and beta-carbon atoms in the test structure and selected stereochemical rules. The focus of this paper is on the use of two packing criteria: avoidance of steric overlap and complete filling of available space. The program also notes potential polar group interactions and disulfide bonds as well as possible burial of formal charges. Central to the algorithm is the side-chain rotamer library. In an update of earlier studies by others, we show that 17 of the 20 amino acids (omitting Met, Lys and Arg) can be represented adequately by 67 side-chain rotamers. A list of chi angles and their standard deviations is given. The newer, high-resolution, refined structures in the Brookhaven Protein Data Bank show similar mean chi values, but have much smaller deviations than those of earlier studies. This suggests that a rotamer library may be a better structural approximation than was previously thought. In using packing constraints, it has been found essential to include all hydrogen atoms specifically. The "unified atom" representation is not adequate. The permitted rotamer sequences are severely restricted by the main-chain plus beta-carbon atoms of the test structure. Further restriction is introduced if the full set of atoms of the external residues are held fixed, the full-chain model. The space-filling requirement has a major role in restricting the template lists. The preliminary tests reported here make it appear likely that templates prepared from the currently known core structures will be able to discriminate between these structures. The templates should thus be useful in deciding whether a sequence of unknown tertiary structure fits any of the known core classes and, if a fit is found, how the sequence should be aligned in three dimensions to fit the core of that class.(ABSTRACT TRUNCATED AT 400 WORDS)

Algorithms↗
PubMedDate not supplied

Diagnosis of rare dementia syndromes: an algorithmic approach.

The etiology of dementia can be diagnosed in most patients using a standard clinical approach consisting of physical, neurologic, and mental status examinations, and laboratory testing, lumbar puncture, and neuroimaging. In some cases, however, the clinical presentation or historical data are unusual, or the results of the workup are inconclusive or atypical. A rare cause of dementia may then be present and a complicated evaluation may be necessary to identify the specific disease process. A potentially useful approach to the diagnosis of rare dementing disorders consists of a series of diagnostic algorithms. This approach utilizes results of neuroimaging studies to guide the evaluation through additional diagnostic steps such as specific enzymatic or immunologic assays or biopsy of extraneural tissues. The disorders potentially detected by these algorithms typically have unusual clinical features such as early age of onset, abnormal neurologic signs and symptoms early in the clinical course, early personality and mood changes, extrapyramidal or cerebellar signs and symptoms, seizures, peripheral neuropathy or myopathy, and extraneural abnormalities involving the dermatologic, cardiovascular, musculoskeletal, or ocular systems. Accurate diagnosis of these rare causes of dementia is important for medical and psychiatric management, prognosis, and genetic counseling.

Aged↗
PubMed1987

Recognition of characteristic patterns in sets of functionally equivalent DNA sequences.

An algorithm has been developed for the identification of unknown patterns which are distinctive for a set of short DNA sequences believed to be functionally equivalent. A pattern is defined as being a string, containing fully or partially specified nucleotides at each position of the string. The advantage of this 'vague' definition of the pattern is that it imposes minimum constraints on the characterization of patterns. A new feature of the approach developed here is that it allows a 'fair' simultaneous testing of patterns of all degrees of degeneracy. This analysis is based on an evaluation of inhomogeneity in the empirical occurrence distribution of any such pattern within a set of sequences. The use of the nonparametric kernel density estimation of Parzen allows one to assess small disturbances among the sequence alignments. The method also makes it possible to identify sequence subsets with different characteristic patterns. This algorithm was implemented in the analysis of patterns characteristic of sets of promoters, terminators and splice junction sequences. The results are compared with those obtained by other methods.

Algorithms↗
PubMed1992

Maternal serum screening for fetal Down syndrome in women less than 35 years of age using alpha-fetoprotein, hCG, and unconjugated estriol: a prospective 2-year study.

OBJECTIVE: To evaluate prospectively maternal serum screening with alpha-fetoprotein (AFP), hCG, and unconjugated estriol (uE3) as a screen for fetal Down syndrome. METHODS: Women less than 35 years of age were offered screening between 15-20 weeks' gestation. Screening results calculated by an algorithm to be equal to or greater than 1:274 (the risk of a 35-year-old for fetal Down syndrome at the second trimester) were considered positive. If gestational age was confirmed by ultrasonography, genetic counseling and amniocentesis were offered. RESULTS: In the first 2 years of our program, 9530 women were screened, of which 686 (7.2%) were found to be screen-positive. Ultrasonographic examination explained the abnormal values in 379 (4.0%). The remaining 307 (3.2%) received genetic counseling and 214 (2.2%) elected amniocentesis or CVS. Four cases of fetal Down syndrome and one de novo chromosomal marker were detected. In three additional cases of fetal Down syndrome, triple-analyte screening failed to identify the pregnancies to be at increased risk. None of the seven cases of fetal Down syndrome would have been detected through screening with maternal serum alpha-fetoprotein (MSAFP) and age alone. CONCLUSIONS: Measurement of MSAFP, hCG, and uE3 in women less than 35 years old is an effective screening test for fetal Down syndrome, with a sensitivity of 57% in our study and an amniocentesis rate (false-positive rate) of 3.2%.

Adult↗
PubMed1987

Impact of plasmids and genetic change on the numerical classification of staphylococci.

Newly isolated bacterial strains often contain extrachromosomal DNA as plasmid DNA. These accessory components of the DNA gene pool confer additional phenotypic properties on their host but, despite this, little attention has been paid to the impact of plasmid-mediated characters on bacterial classification. In the present study, the effect of antibiotic resistance plasmids on the classification of representative staphylococci was determined using numerical phenetic techniques. Over sixty percent of the eighty-one test strains contained one or more plasmids which varied in molecular weight from 1.4 to 36 Mdal. Antibiotic resistance phenotypes were eliminated from strains of S. aureus, S. chromogenes, S. cohnii, S. hyicus and S. xylosus, and from a laboratory isolate, to give sixteen derivative strains. Fourteen had lost one or more plasmids and two had deleted plasmids. In addition three further derivative strains were isolated which showed no plasmid loss but exhibited gross phenotypic changes. The test and derivative strains were the subject of numerical phenetic analyses based on seventy-eight unit characters. Data were examined using the simple matching, Jaccard and pattern coefficients and clustering achieved using the unweighted pair group method with arithmetic averages algorithm. Cluster composition was not markedly affected by the statistics used or by test error, estimated at 1.02%. Numerically circumscribed clusters and subclusters were equated with the established species S. aureus, S. chromogenes, S. cohnii, S. hyicus, S. lentus, S. intermedius, S. sciuri and S. xylosus. The sixteen derivative strains with either lost or delected plasmids were recovered in the same cluster or subcluster as their corresponding parent indicating that the removal of plasmid-expressed characters had little effect on the structure of the numerical classification. In contrast, two of the three strains of S. xylosus with genomically-derived phenotypic variation formed a cluster that separated from their parent strain at the 70% similarity level in the SSM, UPGMA analysis.

Animals↗
PubMedDate not supplied

Aligning amino acid sequences: comparison of commonly used methods.

We examined two extensive families of protein sequences using four different alignment schemes that employ various degrees of "weighting" in order to determine which approach is most sensitive in establishing relationships. All alignments used a similarity approach based on a general algorithm devised by Needleman and Wunsch. The approaches included a simple program, UM (unitary matrix), whereby only identities are scored; a scheme in which the genetic code is used as a basis for weighting (GC); another that employs a matrix based on structural similarity of amino acids taken together with the genetic basis of mutation (SG); and a fourth that uses the empirical log-odds matrix (LOM) developed by Dayhoff on the basis of observed amino acid replacements. The two sequence families examined were (a) nine different globins and (b) nine different tyrosine kinase-like proteins. It was assumed a priori that all members of a family share common ancestry. In cases where two sequences were more than 30% identical, alignments by all four methods were almost always the same. In cases where the percentage identity was less than 20%, however, there were often significant differences in the alignments. On the average, the Dayhoff LOM approach was the most effective in verifying distant relationships, as judged by an empirical "jumbling test." This was not universally the case, however, and in some instances the simple UM was actually as good or better. Trees constructed on the basis of the various alignments differed with regard to their limb lengths, but had essentially the same branching orders. We suggest some reasons for the different effectivenesses of the four approaches in the two different sequence settings, and offer some rules of thumb for assessing the significance of sequence relationships.

Amino Acid Sequence↗
PubMed2026

An end-to-end computational framework for "Record-seq" transcriptional recording data.

MOTIVATION: Record-seq captures cumulative transcriptional activity over time in engineered Escherichia coli by integrating cellular RNA-derived spacer sequences into clustered regularly interspaced short palindromic repeats (CRISPR) arrays, which are read out by sequencing. Unlike the approximately uniform transcript sampling of RNA-seq, Record-seq records biological signal as spacers sampled by the CRISPR spacer acquisition machinery. Consequently, standard RNA-seq analysis strategies are not directly applicable, limiting sensitivity and interpretability. Our previous pipeline addressed these challenges only partially, retained inherited RNA-seq assumptions, and had limited algorithmic efficiency. RESULTS: Here, we present an end-to-end computational framework for Record-seq data. To address the primary computational bottleneck of spacer sequence extraction, we implemented a wavefront alignment approach for efficient quasi-local pattern matching, achieving an approximately 30-fold speedup. We introduce transcription unit-based feature counting as an alternative to gene-body quantification to better represent prokaryotic transcription and increase statistical power by capturing signal from untranslated regions, which are spacer acquisition hotspots. For downstream analyses, we incorporate multiple normalization strategies and a nonparametric differential expression testing framework designed for sparse datasets. Further, we analyze spacer acquisition patterns and train sequence-based neural models that predict acquisition propensity from genomic sequence and annotations, providing a framework for assessing whether acquisition rules generalize as Record-seq is extended to new microbial hosts. AVAILABILITY AND IMPLEMENTATION: The primary analysis workflow, the recoRdseq package, acquisition modeling repository, and relevant data are all linked at https://github.com/plattlab/Record-seq-Framework. Acquisition models and training data are on Zenodo at https://doi.org/10.5281/zenodo.18891434.

Escherichia coli↗
PubMed2026

NLCD: A method to discover nonlinear causal relations among genes.

Distinguishing correlation from causation is a fundamental challenge in many scientific fields, including biology, especially when interventions like randomized controlled trials are infeasible and only observational data are available. Methods based on statistical tests of conditional independence within the Mendelian Randomization framework can detect causality between two observed variables that are each associated with a third instrumental variable. However, these methods for detecting causal relationships between traits (e.g., two gene expression or clinical traits associated with a genetic variant, all observed in the same population) often assume a linear relationship, thereby hindering the discovery of causal gene networks from genomics data. We have developed NLCD, a method for NonLinear Causal Discovery from genomics data based on nonlinear regression modeling and conditional feature importance scoring. NLCD uses these techniques to extend the statistical tests in an existing linear causal discovery method called the Causal Inference Test (CIT). We benchmarked NLCD against current state-of-the-art methods: CIT, Findr, and MRPC. On simulated datasets, NLCD performs comparably to most methods in detecting linear relations (Average AUPRC (Area Under the Precision-Recall Curve) of NLCD = 0.94, CIT = 0.94, Findr = 0.94, and MRPC = 0.99), and outperforms them in detecting nonlinear (sine and sawtooth type) relations between two genes (Average AUPRC of NLCD = 0.76, CIT = 0.60, Findr = 0.56, and MRPC = 0.73). When tested on a nonlinear subset of a yeast genomic dataset to recover known causal relations involving transcription factors, NLCD and CIT performed comparable to each other and slightly better than Findr and MRPC (Average AUPRC of NLCD = 0.82, CIT = 0.81, Findr = 0.71, and MRPC = 0.54). On application to a human genomic dataset, NLCD revealed active causal gene pairs (IRF1 → PSME1 and HLA-C → HLA-T) in the muscle tissue, and clarified the promises and challenges in discovering causal gene networks in tissues under in vivo human settings.

Humans↗
PubMed1990

Real-time multiprocessing of slit scan chromosome profiles.

The multiprocessor NERV and its application to slit scan flow cytometry is described. Up to 320 processors and 640 MBytes of RAM may be used in one VME crate, providing a computing power of less than or equal to 1300 MIPS. The multiprocessor is controlled by a host computer that provides a friendly user interface and comfortable program development tools. All hardware and software has been tested on a prototype NERV system with 5 processors. For a real-time classification/detection of normal and aberrant chromosomes, the centromeric index or the number of centromeres are computed or specifically labeled DNA sequences are detected. The program is partitioned into 60 tasks that can be executed concurrently. A total analysis time of less than 600 microseconds including system overhead will be achieved according to timing measurements which have been done for all individual tasks.

Algorithms↗
PubMed2025

Genetically distinct within-host subpopulations of hepatitis C virus persist after Direct-Acting Antiviral treatment failure.

Analysis of viral genetic data has previously revealed distinct within-host population structures in both untreated and interferon-treated chronic hepatitis C virus (HCV) infections. While multiple subpopulations persisted during the infection, each subpopulation was observed only intermittently. However, it was unknown whether similar patterns were also present after Direct-Acting Antiviral (DAA) treatment, where viral populations were often assumed to go through narrow bottlenecks. Here we tested for the maintenance of population structure after DAA treatment failure, and whether there were different evolutionary rates along distinct lineages where they were observed. We analysed whole-genome next-generation sequencing data generated from a randomised study using DAAs (the BOSON study). We focused on samples collected from patients (N=84) who did not achieve sustained virological response (i.e., treatment failure) and had sequenced virus from multiple timepoints. Given the short-read nature of the data, we used a number of methods to identify distinct within-host lineages including tracking concordance in intra-host nucleotide variant (iSNV) frequencies, applying sequenced-based and tree-based clustering algorithms to sliding windows along the genome, and haplotype reconstruction. Distinct viral subpopulations were maintained among a high proportion of individuals post DAA treatment failure. Using maximum likelihood modelling and model comparison, we found an overdispersion of viral evolutionary rates among individuals, and significant differences in evolutionary rates between lineages within individuals. These results suggest the virus is compartmentalised within individuals, with the varying evolutionary rates due to different viral replication rates and/or different selection pressures. We endorse lineage awareness in future analyses of HCV evolution and infections to avoid conflating patterns from distinct lineages, and to recognise the likely existence of unsampled subpopulations.

Humans↗
PubMed2026

Precision periodontology in clinical practice: bridging omics and clinical decision-making.

BACKGROUND: Precision periodontology integrates molecular diagnostics, genomics, and advanced imaging into clinical decision-making. Despite major advances in microbiome characterisation, host genetics, and inflammatory biomarkers, their translation into routine care remains limited. OBJECTIVES: To critically appraise current evidence on microbiome-based profiling, genetic and epigenetic markers, host-response biomarkers, and three-dimensional imaging in periodontology, and to propose a conceptual decision-support framework linking diagnostic outputs to potential therapeutic actions and future implementation research. MATERIALS AND METHODS: A narrative review searching PubMed/MEDLINE, Scopus, Embase, and the Cochrane Library (2010-2025) using terms related to precision periodontology, subgingival microbiome, periodontitis genetics and epigenetics, salivary and GCF biomarkers, aMMP-8, CBCT, risk assessment, and artificial intelligence. Priority was given to meta-analyses, systematic reviews, longitudinal studies, and guideline documents. RESULTS: Microbiological testing has defined but narrow indications; single-SNP genotyping has not demonstrated clinical utility commensurate with cost; aMMP-8 point-of-care testing is among the most extensively investigated host-response tools and may have adjunctive value in selected monitoring and peri-implant scenarios; however, current evidence remains insufficient to support routine diagnostic implementation. CBCT may directly influence surgical decision-making through defect morphology characterisation. AI-based models show promise but lack prospective clinical validation. These conclusions are consistent with the 20th EFP Workshop Consensus Report. CONCLUSIONS: Precision periodontology currently operates in addition to, rather than in replacement of, conventional staging and grading. We propose a conceptual decision-threshold framework for the selective consideration of molecular and advanced imaging tools when their additive contribution may meaningfully inform management. This framework should be regarded as a research-oriented decision-support model rather than a validated clinical algorithm. CLINICAL RELEVANCE: Clinicians are provided with a structured, evidence-based framework that identifies specific clinical scenarios where molecular diagnostics, host-response biomarkers, and three-dimensional imaging may meaningfully modify periodontal treatment decisions, supporting the operationalisation of precision approaches in daily practice.

Humans↗
PubMed1992

Biostatistical evaluation of blood group, HLA and DNA findings in Jeffreys' immigrant test-case using special kinship and DNA algorithms.

In the immigration case cited by Jeffreys et al. (1985), the biostatistical evaluation of blood group findings in 16 systems and of HLA-A,B findings for the mother and child, using a special kinship algorithm developed by Ihm and Hummel (1975), produced a probability of maternity of the Ghanaian-born putative mother of W = 6%; the probability of maternity for her sister was W = 94%. Using the DNA multilocus probes 33.15 and 33.6 and the bandsharing technique, the authors analysed band patterns from the putative mother and child as well as another 3 children of the same woman. It was concluded that the putative mother was the mother of all 4 children. An evaluation of the band patterns using the multi-di-allelism model and the kinship algorithm in accordance with the Essen-Möller principle produced: W = 99.99999999999999991%, or, if the "-----" constellations were not considered, W = 99.99998%.

Adult↗
PubMed1990

Interrelationships among major protistan groups based on a parsimony network of 5S rRNA sequences.

To test the validity of the maximum parsimony approach to discern protistan interrelationships, we have derived an optimal network of 16S-like rRNA sequences using our parsimony algorithm and compared it with those reported using the distance matrix method. We have also derived an optimal network topology of 50 5S rRNA sequences through an interactive search using our algorithm. In both these networks, the kinetoplastids and euglenoids form a linkage group with Dictyostelium emerging from its neighbourhood. The cryptophytes, dinoflagellates and chromophytes and green algae emerge as independent lines suggesting that plastids arose more than once during protistan evolution. The large 5S rRNA tree further indicates independent origins of mesozoa and metazoa; kinetoplastids and ciliates; and diphyletic origin of fungi. Comparatively close positions of charales and land plants, chytrids and Zygomycetes, Physarum and amoeba, and red algae and green algae are also seen in this network.

Animals↗
PubMed1989

A novel method for promoter search enhanced by function-specific subgrouping of promoters--developed and tested on E.coli system.

A new method for evaluating some complex characteristics of the primary structure of E.coli promoters is proposed. The method, of nonparametric statistical significance, selects important conserved single-base positions in combination with 2-base coupling relations of identity and complementarity. The extended consensus of promoter characteristics thus obtained was used to scan unknown sequences for similarity with E.coli promoters. In terms of this method, a complete set of 244 E.coli promoters was shown to be structurally inconsistent. The set was then broken down into functionally homogeneous subsets of promoters to enhance the selectivity of the search for E.coli-specific promoter sequences, with a high significance level being attained.

Algorithms↗
PubMed2026

Zone equalisation normalisation for improved alignment of epigenetic signal.

MOTIVATION: High-throughput genomic technologies have transformed our understanding of biological systems, yet direct comparison and visualisation of these complex datasets remains challenging. Existing normalisation methods often fail to align genomic signal across samples due to sensitivity to sequencing depth differences and localised high-signal artefacts, leading to inconsistent replicate behaviour and increased downstream variability. RESULTS: We introduce Zone Equalisation Normalisation (ZEN), a novel approach designed to improve cross-sample signal alignment of genomic data. ZEN rescales genomic signal based on variance estimated within biologically enriched regions, reducing the influence of extreme outliers while preserving underlying biological structure. Using a diverse collection of data and our new genome-wide benchmarking approach, we reveal that ZEN improves biological and technical replicate alignment across the majority of tested conditions and experimental platforms. We further show that this improved signal comparability is associated with fewer differential accessibility calls between technical replicates and a more conservative set of biological differences. Together, these results demonstrate that ZEN provides a complementary framework to improve the accuracy and reliability of genomic data analysis and that normalisation choice can affect downstream analyses and biological interpretation. AVAILABILITY AND IMPLEMENTATION: ZEN is available as an open-source Python package via conda and PyPI. Source code, documentation, tutorials, and code to reproduce the analyses are available at https://github.com/Genome-Function-Initiative-Oxford/Zone-Equalisation-Normalisation and Zenodo (https://doi.org/10.5281/zenodo.21067751).

Epigenesis, Genetic↗
PubMed2026

Episode clustering in phylogenetic networks.

MOTIVATION: The classical duplication episode clustering (EC) model introduced by Guigó et al. in the 1990s provides a foundational approach for inferring genomic duplication events crucial to understanding genome evolution. This model clusters single gene duplications from a collection of gene trees at locations in the species tree to minimize the total number of such locations, called duplication episodes. However, it does not capture reticulate evolutionary histories. RESULTS: Here, we introduce NetEC, a novel extension of this problem to phylogenetic networks. To solve NetEC, we first develop a polynomial-time dynamic programming (DP) algorithm for testing whether a given set of network nodes can serve as episode locations. We then propose a main inference algorithm that utilizes this DP component to optimize the episode count; while the feasibility test runs in polynomial time, the full optimization has exponential worst-case complexity, and an optional heuristic mode is provided for larger instances. We also propose an extended episode analysis procedure that identifies additional genomic duplication candidates below reticulation nodes, complementing the main algorithm by resolving potential upward clustering of duplications induced by reticulation. We evaluate our method on simulated data and on an empirical Pandanales dataset comprising over 29 000 gene trees, demonstrating exact and accurate inference of genomic duplication events even in the presence of multiple reticulations. AVAILABILITY AND IMPLEMENTATION: All experiments were conducted using the NetEC tool (https://github.com/ppgorecki/netec), with all input data, scripts, and parameter settings for reproduction available in the same repository.

Phylogeny↗