Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “incomplete lineage sorting”

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.

17 recordsLinked to original sources
PubMed2026

Complex evolutionary history of Rosales mediated by extensive incomplete lineage sorting and hybridization.

The angiosperm order Rosales still represents a major challenge for phylogenetic reconstruction. Although its circumscription is now well-defined, phylogenetic relationships among families are still uncertain. Here, we used nuclear, plastid, and mitochondrial genomic data from 33 species representing all nine families to further clarify interfamilial relationships and the group's evolutionary history. We detected significant phylogenetic conflict among the three datasets. Further analyses at the nuclear level identified incomplete lineage sorting (ILS) as the main cause of unstable phylogenetic positions among families. The discordant placements of Rhamnaceae and Elaeagnaceae based on plastid and mitochondrial data are caused by ancient hybridization events, potentially involving differences in organellar inheritance. Our molecular dating confirms earlier suggestions that the ancient rapid diversification of the three Rosaceae subfamilies could be the main reason for the difficulties in resolving their phylogenetic relationships. Our findings provide new insights into the interfamilial relationships of Rosales and demonstrate that the evolutionary history of this order was shaped by ancient and rapid radiation as well as extensive ILS and reticulate evolution. They also suggest that previous attempts to clarify interfamilial relationships in this order were hampered by combining nuclear and organellar sequence data, leading to inconsistent topologies observed across earlier studies.

Phylogeny↗
PubMed2026

Phylotranscriptomics Allows Distinguishing Major Gene Flow Events from Incomplete Lineage Sorting in Rapidly Diversifying Mimetic Orchids (Genus Ophrys).

Ophrys orchids (or bee orchids) provide an outstanding example of a plant adaptive radiation. Over the last 5 million years, this genus has diversified into hundreds of taxa as a result of its unconventional pollination strategy, known as "sexual swindling". However, the rapid and substantial diversification of this genus, combined with its capacity for hybridization and large genome size, poses significant challenges in addressing its systematics. We used phylotranscriptomics as a genome complexity reduction technique to infer the phylogenetic relationships among Ophrys main lineages. More than seven thousand gene trees enabled us to determine the relative contributions of gene flow and incomplete lineage sorting (ILS) in Ophrys evolution. First, we propose a new phylogenetic hypothesis for the genus with an unprecedented resolution that largely confirms the relationships between the main Ophrys lineages, but also provides new insights within each subgenera. By combining phylogenetic network inference with introgression analyzes based on gene tree topologies and branch lengths, we then show that the numerous phylogenetic incongruences among gene tree topologies result from a pervasive background of ILS, over which stand out several well-supported, ancient and potentially adaptive gene flow events between lineages. These major gene flow events provide a new perspective on the evolution of the Ophrys genus and its pollination, questioning previous hypotheses inferred without considering its reticulate evolution, and providing a better understanding of discrepancies observed among previous phylogenetic studies of the genus.

Orchidaceae↗
PubMed2025

Ancient Introgression Explains Mitochondrial Genome Capture and Mitonuclear Discordance Among South American Collared Tropidurus Lizards.

Mitonuclear discordance-evolutionary discrepancies between mitochondrial and nuclear DNA phylogenies-can arise from various factors, including introgression, incomplete lineage sorting, recent or ancient demographic fluctuations, sex-biased dispersal asymmetries, among others. Understanding this phenomenon is crucial for accurately reconstructing evolutionary histories, as failing to account for discordance can lead to misinterpretations of species boundaries, phylogenetic relationships, and historical biogeographic patterns. We investigate the evolutionary drivers of mitonuclear discordance in the Tropidurus spinulosus species group, which contains nine species of lizards inhabiting open tropical and subtropical environments in South America. Using a combination of population genetic and phylogenomic approaches applied to mitochondrial and nuclear data, we identified different instances of gene flow that occurred in ancestral lineages of extant species. Our results point to a complex evolutionary history marked by prolonged isolation between species, demographic fluctuations, and potential episodes of secondary contact with genetic admixture. These conditions likely facilitated mitochondrial genome capture while diluting signals of nuclear introgression. Furthermore, we found no strong evidence supporting incomplete lineage sorting or natural selection as primary drivers of the observed mitonuclear discordance. Therefore, the unveiled patterns are most consistent with neutral demographic processes, coupled with ancient mitochondrial introgression, as the main factors underlying the mismatch between nuclear and mitochondrial phylogenies in this system. Future research could further explore the role of other demographic processes, such as asymmetric sex-biased dispersal, in shaping these complex evolutionary patterns.

Animals↗
PubMed2026

Beyond species trees: pervasive gene flow limits phylogenomic resolution in the diversification of Juniperus from the Qinghai-Tibet Plateau.

Understanding how lineages diversify despite persistent ancestral polymorphism and recurrent gene flow remains a central challenge in evolutionary biology. Juniperus distributed across the Qinghai-Tibet Plateau provide an ideal system for addressing this question because repeated geological uplift and climatic oscillations have likely promoted cycles of lineage divergence, range shifts, and secondary contact. Here, we combined approximately 1.08 million genome-wide SNPs from 164 individuals representing thirteen Juniperus lineages with phylogenomic datasets comprising 3,381 nuclear single-copy genes and nearly complete plastomes. We detected extensive phylogenomic discordance and cytonuclear incongruence across genomic datasets. Topology weighting, coalescent simulations, quartet-based tests, and analyses of gene flow and reticulation collectively support the interpretation that these patterns were shaped by the combined effects of prolonged incomplete lineage sorting and gene flow during lineage diversification. Ecological niche analyses further provide a spatial and climatic context in which environmentally similar lineages may have had greater opportunities for secondary contact during historical range shifts. Collectively, our results reveal that the evolutionary history of Qinghai-Tibet Plateau Juniperus is characterized by reticulate diversification rather than strictly bifurcating evolution, and demonstrate how genome-wide discordance can provide biological insights into the evolutionary processes underlying lineage diversification.

Gene Flow↗
PubMed2026

Introgression shapes the genomic conflict landscape of Malus, providing evidence for a reticulate backbone in a woody crop lineage.

Phylogenomic discordance is widespread across plants, but its evolutionary significance is often obscured when conflict is treated primarily as analytical noise rather than as evidence of underlying processes. In woody lineages in particular, incomplete lineage sorting, introgression, and genome duplication can interact over long timescales to produce complex genomic histories that are not adequately summarized by a strictly bifurcating tree. Here, we use Malus as a model woody genus to investigate how these processes structure conflict across a genus-scale, accession-based phylogenomic framework. Using broad taxon sampling, hundreds of nuclear loci, plastid genomes, and genome-wide SNP summaries, we reconstruct a robust nuclear backbone for sampled Malus lineages and evaluate where discordance is concentrated and which processes best explain it. Nuclear analyses resolve eight major clades, whereas conflict is non-random and localized to recurrent hotspots rather than evenly distributed across the tree. Cytonuclear discordance is similarly concentrated, especially around Clade H, represented by sampled accessions of M. tschonoskii, where localized plastid-nuclear disagreement is consistent with candidate plastid capture or organellar introgression. Multiple complementary analyses further indicate that the strongest conflict is not explained by ILS alone, but instead reflects lineage-structured introgression, while polyploid complexes represent additional localized sources of evolutionary complexity. Together, these results provide evidence for a reticulate genomic backbone in Malus and show how integrating nuclear, plastid, and genome-wide conflict analyses can help distinguish background discordance from process-specific signals in woody plant radiations. Several lineage-level reticulation hypotheses identified here should now be tested with broader population-level sampling and curated reference accessions.

Malus↗
PubMed2026

Cytonuclear conflict and reticulate evolution in the Morelloid clade (Solanum, Solanaceae): Insights from genome skimming and network Phylogenomics.

The Morelloid clade (black nightshades) is one of the most strongly supported clades within the megadiverse Solanum genus. It comprises 76 globally distributed, non-spiny herbaceous and suffrutescent species. While often erroneously considered poisonous weeds, several species are economically important as orphan crops. The clade is closely related to tomato and potato but, due to a lack of focused breeding efforts, remains a putative reservoir of genetic diversity for crop improvement. Despite this potential, we lack fundamental knowledge on the evolution of the Morelloid clade. The group includes polyploid species with unknown parental origins-likely reflecting reticulate processes such as hybridization, introgression, and associated backcrossing events. Prior analyses have been unable to disentangle these processes, leaving the mechanisms underlying reticulate evolution in the Morelloid clade poorly understood. Here, we use genome skimming to produce a well-supported maximum likelihood plastid phylogeny from complete circularized plastomes and a coalescent-based species tree from combined Angiosperms353 and conserved ortholog set nuclear markers. Our dataset, composed of previously published data and deep genome skimming from herbarium samples, spans 26 Morelloid species. To investigate phylogenetic discordance, we used a nuclear phylogenetic network, multispecies coalescent simulations, a fused rooted nuclear chloroplast tree, and quantification of nuclear gene tree concordance. We show that incongruence between nuclear and plastid trees is pervasive and cannot be explained by incomplete lineage sorting alone. Instead, our results demonstrate that events consistent with repeated chloroplast capture have shaped the reticulate evolutionary history of the clade, especially among African polyploid and Pan-American diploid lineages.

Phylogeny↗
PubMed2026

Genome-wide insights into the evolutionary and demographic history of the red alga Mazzaella laminarioides: Evidence for speciation with ancient migration along the southeast Pacific coast.

The mechanisms driving lineage divergence in red algae remain unexplored, despite the group's remarkable diversity and ancient evolutionary history. The red alga Mazzaella laminarioides, a Chilean intertidal species complex composed of three parapatric cryptic lineages (North, Center, South), offers a valuable system to evaluate these processes, as its life history combines severe dispersal limitation with a haploid-diploid cycle that may influence the emergence of reproductive barriers. We reconstructed its evolutionary history using whole-genome sequencing and nuclear genome assembly of representative individuals from each lineage. Phylogenomic analyses based on 1,507 single-copy orthologs recovered three deeply divergent lineages with limited nuclear discordance consistent with incomplete lineage sorting. For both splits, demographic modelling was most consistent with an Ancient Migration scenario, although support over strict isolation was moderate, suggesting that divergence may have begun with low asymmetric ancestral gene flow followed by subsequent loss of connectivity, demographic bottlenecks, and later population expansion. Coding sequence analyses revealed lineage-specific dN/dS heterogeneity; only one South-lineage locus passed FDR correction (metaxin-1, mitochondrial protein import), with two further South-lineage candidates in chlorophyll and heme biosynthesis falling below the FDR threshold. Together, these signals suggest that divergent selective pressures on energy acquisition may have contributed to divergence at the southern end of the distribution. These results add to the small but growing body of whole-genome data for red algae and, alongside recent macroalgal studies, suggest that ancestral connectivity could be a recurrent feature of lineage divergence even in marine organisms with extremely restricted dispersal.

Rhodophyta↗
PubMed2026

Phylogenomics and evolution of the Lauraceae based on targeted capture data.

The family Lauraceae, a hyper-diverse magnoliid family comprising approximately 63 genera and over 3,000 species, plays a key ecological role in tropical and subtropical forests. Yet deep relationships among its nine tribes remain unresolved, likely due to limited sampling and complex evolutionary processes such as incomplete lineage sorting (ILS) and gene flow. To address these challenges, we generated datasets of 255 single-copy nuclear genes and chloroplast genomes using a newly designed Lauraceae-specific probe set, achieving the most comprehensive genus-level sampling (84%) to date. Phylogenomic analyses reconstructed a robust nuclear tree, which resolved the Neocinnamomeae as sister to the Caryodaphnopsideae and revealed pronounced gene tree conflict and pervasive cytonuclear discordance. To investigate the evolutionary processes underlying these patterns, comprehensive analyses were conducted. The results indicate that conflicting nuclear gene trees reflect the combined effects of ILS, gene tree estimation error, and gene flow, with ILS dominating across the core Lauraceae, whereas cytonuclear discordance is primarily driven by extensive gene flow. Diversification analyses further indicate that episodes of rapid lineage accumulation coincide with major gene flow events, suggesting a potential role of gene flow in the diversification of Lauraceae. Overall, this study provides a robust nuclear phylogenomic framework for Lauraceae and demonstrates that gene flow had profound effects on its evolutionary history, shedding light on the contribution of gene flow to the diversification of hyper-diverse tropical plant lineages.

Cytonuclear discordance↗
PubMed2026

Nuclear single-copy orthologous genes as phylogenomic markers for resolving the closely related firefly genera Pteroptyx, Medeopteryx, and Trisinuata (Coleoptera: Lampyridae: Luciolinae).

Fireflies (Lampyridae) are bioluminescent beetles with broad ecological roles across temperate and tropical ecosystems, occupying diverse habitats including forests, wetlands, grasslands, mangroves, and riverine systems. The subfamily Luciolinae is primarily distributed across Asia and the Indo-Pacific. Phylogenetic relationships among three closely related Luciolinae genera - Medeopteryx, Pteroptyx, and Trisinuata - remain unresolved using mitochondrial genome data alone. This study used nuclear genome data to resolve relationships among these genera and identify a lighter-weight nuclear marker panel for expanding taxon sampling. Draft genomes were reconstructed for fifteen firefly species, eight from the focal genera, and analyzed with five published firefly genomes. Using BUSCO and OrthoFinder, 1,011 nuclear single-copy orthologs (SCOs) were identified for phylogenomic inference. Discordance between concatenation- and coalescence-based phylogenies indicated incomplete lineage sorting (ILS). The coalescence-based phylogeny recoveredPteroptyxas monophyletic and sister to a (Medeopteryx,Trisinuata) clade, with Trisinuata nested within a non-monophyletic Medeopteryx; however, quartet support at the base of Pteroptyx, particularly at Pt. valida, was low.Filtering for compositional homogeneity, clock-likeness, and species-tree concordance yielded 103 SCOs with a significantly higher proportion of parsimony-informative sites than non-selected loci, retaining the backbone topology with higher gene concordance support at scored clades, while ILS-driven discordance at Pt. valida persists - confirming that the reduced panel retains phylogenetic resolving power for future taxon sampling. These findings demonstrate a practical framework for using nuclear SCOs to resolve close phylogenetic relationships within Luciolinae. Future work should expand taxon sampling - especially forTrisinuata - alongside long-read assemblies, for a more robust phylogenomic framework.

Fireflies↗
PubMed2025

Hybridization as driving force for cryptic species diversity in the Caribbean coral genus Madracis.

Species boundaries in scleractinian corals remain highly elusive due to conflicting patterns between morphological and molecular phylogenies, often caused by morphological plasticity, occurrence of cryptic species, incomplete lineage sorting or introgressive hybridization. Here, we use an integrated systematics approach, which combines reduced representation genome sequencing (nextRAD), micro-morphometric characterization, SEM analyses and compilation of life history traits, to infer phylogenetic relationships among closely related species in the Caribbean coral genus Madracis. In total, we analyzed 235 Madracis specimens from Curaçao and Bermuda collected from 10-90 m depth. Sequence- and SNP-based analyses for 115 samples generated unprecedented species resolution in Madracis, greatly supporting the morphology-based taxonomy of the current, accepted Caribbean species M. senaria, M. decactis, M. formosa, M. carmabi and M. mirabilis (M. auretenra). The exception was M. pharensis, in which we found evidence for three separate lineages, and for which we found signatures of admixture and introgression. These three M. pharensis lineages showed distinct depth distributions (thus classified as shallow, deep and very deep) and were partially distinguishable on the basis of fine microstructural elements of the collumella, septa and coenosteum. Further taxonomic comparisons are needed to formalize these putative cryptic species. Overall, our integrated systematics approach further resolves species relationships in the Caribbean genus Madracis, supports the morphological descriptions for most of the recognized species, but also reveals the existence of cryptic diversity in groups marked by high admixture, thus suggesting hybridization as a driving force in coral species diversity.

Animals↗
PubMed2026

Genomic Footprints of Historical Introgression Between Ancient Lineages of Wild Oryza AA-Genome Species With Widely Separated Contemporary Distributions.

Phylogenetic incongruence is increasingly recognized as pervasive, yet the extent to which reticulate evolution occurs between groups separated by substantial geographical distances and deep phylogenetic divergence remains poorly characterized. In the Oryza AA-genome group-a model for plant speciation and domestication-the traditional bifurcation model posits that Australian Oryza meridionalis and African Oryza longistaminata occupy basal branches, distinct from the more recently diversified monophyletic clade comprising Asian and other African lineages, including major cultivars. However, recent evidence from endogenous viral sequences has hinted at unexpected genetic relatedness between African O. longistaminata and Asian Oryza sativa, which are geographically and phylogenetically distant. Here, we conducted a genome-wide survey across 11 Oryza species to systematically identify genomic regions exhibiting phylogenetic incongruence. Widespread phylogenetic discordance was observed, notably involving genomic segments in which O. longistaminata showed phylogenetic proximity to Asian species, contradicting their established deep divergence. To distinguish between introgression and incomplete lineage sorting, we performed four-taxon ABBA-BABA tests, which provided statistical support for introgression. Furthermore, divergence time estimates for these incongruent regions were younger than the species divergence times, suggesting historical introgression between the ancestors of lineages that are currently separated by vast geographical distances. Systematic assessments indicated that potential analytical artifacts, such as compositional bias and substitution saturation, were unlikely to explain the observations. These convergent lines of evidence suggest that ancient introgression had occurred between currently geographically separated and evolutionarily divergent Oryza lineages, leaving detectable footprints across their modern genomes.

Oryza↗
PubMed2025

Rapid Radiations Outweigh Reticulations During the Evolution of a 750-Million-Year-Old Lineage of Cyanobacteria.

Species are a fundamental unit of biodiversity. Yet, the existence of clear species boundaries among bacteria has long been a subject of debate. Here, we studied species boundaries in the context of the phylogenetic history of Nostoc, a widespread genus of photoautotrophic and nitrogen-fixing cyanobacteria that includes many lineages that form symbiotic associations with plants (e.g. cycads and bryophytes) and fungi (e.g. cyanolichens). We found that the evolution of Nostoc was characterized by eight rapid radiations, many of which were associated with major events in the evolution of plants. In addition, incomplete lineage sorting associated with these rapid radiations outweighed reticulations during Nostoc evolution. We then show that the pattern of diversification of Nostoc shapes the distribution of average nucleotide identities (ANIs) into a complex mosaic, wherein some closely related clades are clearly isolated from each other by gaps in genomic similarity, while others form a continuum where genomic species boundaries are expected. Nevertheless, recently diverged Nostoc lineages often form cohesive clades that are maintained by within-clade gene flow. Boundaries to homologous recombination between these cohesive clades persist even when the potential for gene flow is high, i.e. when closely related clades of Nostoc co-occur or are locally found in symbiotic associations with the same lichen-forming fungal species. Our results demonstrate that rapid radiations are major contributors to the complex speciation history of Nostoc. This underscores the need to consider evolutionary information beyond thresholds of genomic similarity to delimit biologically meaningful units of biodiversity for bacteria.

Phylogeny↗
PubMed2026

Phylogenomic evidence provides insights into the phylogeny and speciation patterns in the subgenus Caloscordum (Allium, Amaryllidaceae).

Phylogenomics with abundant informative sites offers a powerful means for elucidating complex diversification history. Here, we collected 22 samples from 18 populations representing all species of subgenus Caloscordum (Allium). Using transcriptome and whole-genome resequencing data, we generated 1755 low-copy nuclear genes and 81 plastid genes. By integrating morphological and phylogenomic evidence, we clarified the subgenus's complex evolutionary histories and speciation patterns. A total of 18 morphological characteristics were analysed, with a taxonomic framework established. Our analyses resolved robust species relationships despite detecting extensive phylogenetic discordances, which were attributed to incomplete lineage sorting (ILS) and hybridization. Specifically, our results suggest that A. inutile originated via rapid budding speciation from the widespread A. tubiflorum. This process likely coincided with mid-Pleistocene glacial-interglacial cycles and may have been reinforced by geographic isolation and ecological adaptation. In contrast, the sole tetraploid species, A. peikingense, was confirmed to be of hybrid origin, derived from A. neriniflorum and A. tubiflorum. This allopolyploidization event appears to have been facilitated by secondary contact between the parent species, which was likely associated with climatic oscillations within the 35° N-45° N arid belt. Overall, our findings elucidate the intricate speciation patterns within Caloscordum and highlight how the interplay of polyploidization, ecological isolation, and tectonic uplift-driven aridification has shaped plant diversity in East Asia.

Phylogeny↗
PubMed2025

Unraveling evolutionary pathways: allopolyploidization and introgression in polyploid Prunus (Rosaceae).

Allopolyploidization, resulting from hybridization and subsequent whole-genome duplication (WGD), is a fundamental mechanism driving evolutionary diversification across various lineages within the Tree of Life. The polyploid Prunus (Rosaceae), significant for its economic and agricultural value, provides an ideal model for investigating the evolutionary dynamics associated with allopolyploidy. In this study, we utilized deep genome skimming (DGS) data to demonstrate a comprehensive analytical framework for elucidating the underlying allopolyploidy that includes a newly adapted tool (DGS-Tree2GD) tailored explicitly for accurately detecting WGD events. Additionally, we introduced two methods to evaluate the contribution of incomplete lineage sorting (ILS) to lineage diversification. Phylogenomic discordance analyses revealed that allopolyploidization, rather than ILS, played a dominant role in the origin and dynamics of polyploid Prunus. Moreover, we inferred that the uplift of the Himalayas from the Middle to Late Miocene was a key driver in the rapid diversification of the Maddenia clade, an endemic group in East Asia. This geological event facilitated extensive hybridization and allopolyploidization, particularly the introgression between the Himalayas-Hengduan and Central-Eastern China clades. This case study demonstrates the robustness and efficacy of our analytical approach in precisely identifying WGD events and elucidating the evolutionary mechanisms underlying allopolyploidization in polyploid Prunus.

Polyploidy↗
PubMed2026

Phylogenomics reveals persistent gene-tree discordance in the Chenopodium album aggregate.

BACKGROUND AND AIMS: Complex genomic histories shaped by hybridisation and polyploidy can influence traits related to plant defence, stress tolerance and toxicity, particularly in Amaranthaceae, which includes crops such as quinoa and spinach. Within this family, white goosefoot (Chenopodium album), a widespread agricultural weed and traditional food resource, belongs to a diploid-polyploid aggregate with extensive phylogenetic discordance. Clarifying its evolutionary history provides context for interpreting ecologically and agronomically relevant trait variation across the aggregate. Building on the established genome-lineage framework, we tested whether discordance persists when constituent genome-lineage components are represented separately and whether the remaining signal is compatible with reticulate evolution. METHODS: We analysed 2,298 conserved nuclear BUSCO families across 27 assembly-level terminals using tree- and network-based approaches. Genome-lineage-aware analyses used 2,156 families after separating polyploid Chenopodium into A-H components, with Dysphania ambrosioides as outgroup. HyDe tested site-pattern asymmetry under global false-discovery-rate correction. KEY RESULTS: Assembly-level analyses grouped the Danish C. album aggregate accession Ca6-1 with hexaploid C. album sensu stricto dcCheAlbu1.1, whereas relationships among surrounding Chenopodium taxa were less stable. Genome-lineage-aware analyses recovered the expected B-, C- and D-affinity relationships, but substantial gene-family heterogeneity persisted. Reticulate network models fitted the assembly-level data better than bifurcating models, although inferred patterns differed between methods. HyDe detected significant site-pattern asymmetry in a small subset of loci, with most retained signal shared between the focal assemblies. CONCLUSIONS: Gene-tree discordance persists in the C. album aggregate after genome-lineage separation. The established genome-lineage framework captures the dominant phylogenomic structure, while residual heterogeneity is compatible with both tree-like and reticulate processes without identifying direct progenitors or a unique hybridisation history. This framework supports future analyses of lineage-specific and trait-associated loci related to plant defence, food quality and toxicity in C. album and related Amaranthaceae.

Chenopodium album↗
PubMed2026

Phylogenomic subsampling and upsampling for efficient evolutionary analyses of big data.

Long runtimes, high memory demands, and reliance on high-performance computing impede phylogenomic analyses. We review a scalable phylogenomic subsampling with upsampling (PSU) framework to address this challenge, which reduces runtime and memory requirements by orders of magnitude. In PSU, small subsamples of sites from a concatenated alignment are analyzed, which are expanded by upsampling before inference, and the resulting inferences are aggregated to obtain evolutionary estimates. PSU harnesses the fact that the computational cost of maximum likelihood analysis is strongly influenced by the number of distinct site patterns in the concatenated alignment, whereas statistical power depends primarily on the amount of evolutionary information represented by the total number of sites and substitutions. By reducing the former while restoring the latter through upsampling, PSU can approximate many full-alignment analyses at substantially lower computational cost. Analysis of simulated and empirical datasets shows that PSU can accurately estimate bootstrap support values, select the optimal substitution model, test evolutionary hypotheses, and infer branch lengths, divergence times, and associated uncertainty measures. PSU also provides distributions of inferred clade support across independent subsamples, enabling detection of conflicting phylogenetic signals that may remain hidden in conventional bootstrap analysis of concatenated alignments. Automated tuning of subsample size, the number of subsamples, and the number of upsampling replicates make PSU practical. We suggest that PSU is a general approach for scalable phylogenomic inference using a broad range of statistical methods. By enabling analyses of genome-scale alignments on commodity hardware, PSU broadens research access and reduces environmental and infrastructural costs of big-data phylogenomics.

Phylogeny↗
PubMed2026

Phylogenomic subsampling and upsampling for efficient evolutionary analyses of big data.

Long runtimes, high memory demands, and reliance on high-performance computing impede phylogenomic analyses. We review a scalable phylogenomic subsampling with upsampling (PSU) framework, in which small subsamples of sites from a concatenated alignment are expanded by upsampling before inference, and the resulting analyses are then aggregated to obtain evolutionary estimates. PSU harnesses the fact that the computational cost of maximum likelihood analysis is strongly influenced by the number of distinct site patterns in the concatenated alignment, whereas statistical power depends primarily on the amount of evolutionary information represented by the total number of sites and substitutions. By reducing the former while restoring the latter through upsampling, PSU can approximate many full-data analyses at substantially lower computational cost. Analysis of simulated and empirical datasets shows that PSU can accurately estimate bootstrap support values, select the optimal substitution model, test evolutionary hypotheses, and infer branch lengths, divergence times, and associated uncertainty measures, while reducing runtime and memory requirements by orders of magnitude. PSU also provides distributions of inferred clade support across independent subsamples, enabling detection of conflicting phylogenetic signals that may remain hidden in conventional bootstrap analysis. Automated tuning of subsample size, the number of subsamples, and the number of upsampling replicates make PSU practical across diverse datasets. We suggest that PSU is a general strategy for scalable phylogenomic inference using a broad range of statistical methods. By enabling analyses of genome-scale alignments on commodity hardware, PSU broadens research access and reduces environmental and infrastructural costs of big-data phylogenomics.

confidence limits↗