Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “phylogenetic diversity”

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.

259 records · Page 2Linked to original sources
PubMed2026

Viral surveillance beyond detection: JMTV and the need for ensemble approaches in emerging virus discovery.

The recent report by T. Murillo, L. E. Enrique Chaves-González, S. Temmam, S. Bermúdez, et al. (Microbiol Spectr 14:e04078-25, 2026, https://doi.org/10.1128/spectrum.04078-25) expands the known geographic and ecological range of Jingmen tick virus (JMTV) by detecting the virus in Amblyomma mixtum ticks collected from horses in Costa Rica. This is an important finding because A. mixtum can feed on wildlife, domestic animals, and humans, creating a possible interface for virus movement across various hosts. The study also places the Costa Rican virus in a wider phylogenetic context, linking it to JMTV diversity reported from other regions. However, the detection of viral RNA in ticks should not be interpreted as proof of local disease, human infection, or active transmission, especially in the absence of supporting results. Instead, it reflects an important signal for careful viral surveillance. Here, I discuss how JMTV illustrates the need for ensemble approaches that combine field sampling, phylogeny, segment-level genome analysis, serology, experimental validation, and data-driven virus discovery tools.

emerging viruses↗
PubMed2026

The Complete Mitochondrial Genome of a Newly Recorded Chinese Species of Diglyphus sabulosus (Hymenoptera: Eulophidae) and Insights into Its Phylogenetic Position.

Diglyphus Walker, 1844 is an economically important genus which many species acting as biocontrol agents against agromyzid leafminer pests, but there is a lack of mitogenomic data on the evolutionary relationships within this genus, hindering a comprehensive understanding of its evolutionary history. We used traditional morphological methods to identify species, and present the first complete mitochondrial genome sequence and characterization of features of Diglyphus sabulosus and further infer its phylogenetic position based on the amino acid sequences of 13 protein-coding genes (PCGs). The complete mitochondrial genome of D. sabulosus is 15,690 bp in length, including 13 PCGs, 22 transfer RNA genes, 2 ribosomal RNA genes and a control region. The AT content of the whole genome sequence was 81.0%, indicating a significant AT bias. All protein-coding genes have the typical ATN as the start codon and TAA as the stop codon. Phylogenetic analysis inferred from the amino acid sequences of 13 PCGs revealed that all species within the family Eulophidae constituted a monophyletic clade, supporting the monophyly of this family. D. sabulosus and D. poppoea form a well-supported sister group, representing the species with the closest phylogenetic relationship within the analyzed taxa. In this study, the mitogenome structure was analyzed and the taxonomic status of D. sabulosus was clarified, thus providing a theoretical basis for understanding the phylogenetic relationships of Diglyphus.

Animals↗
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

Target Capture of Ancient Shell DNA Enables Phylogenetic Reconstruction of Deep-Sea Molluscs.

Target capture is widely used to enrich endogenous DNA from calcium phosphate skeletal material in vertebrates, but its performance on calcium carbonate hard parts widely produced by invertebrates remains poorly understood. Here, we compared DNA recovery from four fresh and 12 ancient (eight radiocarbon-dated to 1671-1135&#x2009;years old before present) deep-sea vesicomyid clam shells, including species Archivesica marissinica, A. nanshaensis and A. okutanii, using whole-genome sequencing (WGS) or target capture of ultraconserved elements (UCEs). WGS achieved 16.65% on-target read recovery of UCEs from fresh soft tissue, but <&#x2009;1% from shell specimens. By contrast, UCE capture in the same specimen increased on-target reads by up to 155-fold, reaching 29.84% in fresh shells and up to 72-fold, reaching 19.89% in ancient shells. Target capture of UCEs recovered 142-1001 loci per sample compared to 0-230 with WGS alone. Ancient shells of A. marissinica and A. okutanii, based on reads mapped with bwa-mem2 and bbmap, exhibited characteristic post-mortem DNA damage signals, with average 5'-end C-to-T misincorporation rates of 3.46% and 15.97%, respectively, exceeding the levels observed in fresh A. marissinica shells (maximum 1.24%). UCE-based phylogenetic reconstructions incorporating shell ancient DNA recovered two major clades within Pliocardiinae, consistent with published phylogenomic trees. Together, these findings demonstrate that target-capture enrichment enables effective recovery of highly degraded DNA from ancient mollusc shells and supports robust phylogenetic inference at the intrageneric scale, expanding the utility of shells-one of the most abundant invertebrate remains-for evolutionary, biogeographic and conservation studies.

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

Genome-wide scans reveal candidate genes associated with wing morph differentiation in Tetrix japonica.

Wing dimorphism is an important dispersal-related trait in insects, but its genomic basis remains poorly understood in pygmy grasshoppers. Here, we integrated genome-wide single-nucleotide polymorphism (SNP) analyses, population structure inference, selection scans, and functional annotation to investigate genomic differentiation between long- and short-winged Tetrix japonica. Principal component analysis (PCA), ADMIXTURE, and phylogenetic analyses revealed weak genome-wide separation between morphs, indicating differentiation on a largely shared genetic background. Genome-wide scans based on the fixation index (FST), nucleotide diversity ratios, and Tajima's D, using 50-kb non-overlapping windows and empirical top-5% outlier thresholds, identified multiple candidate regions across seven chromosomes. The broader long- and short-winged candidate sets spanned 9.35&#xa0;Mb and 9.37&#xa0;Mb and directly overlapped 82 and 77 genes, respectively. Candidate genes were associated with signaling/hormone regulation, membrane transport, metabolism, cytoskeletal organization, extracellular matrix structure, and development. Short-winged candidate genes were significantly enriched for ABC-type transporter activity and ATP hydrolysis activity. Because all individuals originated from a single laboratory-maintained population with weak genome-wide structure, these regions should be regarded as candidate loci from a screening-stage analysis that require validation in independent populations and by functional assays, rather than as confirmed targets of selection.

Animals↗
PubMed2026

Fructophilic lactic acid bacteria as a window into multi-scale convergent evolution.

Fructophilic lactic acid bacteria (FLAB) are a group of lactic acid bacteria with unique growth characteristics, that is, poor growth on glucose. Their growth is enhanced in the presence of fructose or external electron acceptors. These organisms inhabit fructose-rich environments such as flowers, fruits, and pollinating insects, particularly honey bees. Apilactobacillus spp. and Fructobacillus spp. are representatives of FLAB, although they belong to phylogenetically distant clades. These organisms commonly possess markedly small genomes with a low number of coding DNA sequences. Furthermore, their genomes are characterized by a markedly reduced number of genes involved in carbohydrate transport and metabolism. Genome reduction in FLAB reflects convergent adaptation to fructose-rich environments rather than general genome streamlining. The two distinct FLAB genera, Fructobacillus and Apilactobacillus, independently lost more than 100 genes in statistically similar orders. In contrast, genes involved in carbohydrate and amino acid metabolism exhibited reversed orders of loss between the two genera. Furthermore, FLAB genomes lack an intact bifunctional alcohol/aldehyde dehydrogenase gene (adhE), which causes their poor growth on glucose. A comparative genomic study suggested the evolutionary process underlying adhE gene decay during adaptation to the fructose-rich environments, including pollinating insects. In conclusion, FLAB represent a unique example of habitat-driven convergent reductive evolution that can be investigated across multiple biological scales - from individual genes to whole genomes - in the diverse LAB group with a wide range of habitats, and partially share the fructophilic evolution with eukaryotic yeasts found in fructose-rich habitats.

Fructose↗
PubMed2026

Genome-wide identification and expression profiling of CSP and OBP genes in Stictocephala bisonia reveals candidate genes potentially associated with insecticide response.

Stictocephala bisonia is an important invasive agricultural pest. Due to the frequent application of insecticides in its habitat, this species is under intense selection pressure. Chemosensory proteins (CSPs) and odorant-binding proteins (OBPs) are known to play key roles in insecticide resistance, but their specific functions in S. bisonia remain unclear. In this study, we identified a total of 22 SbisCSPs and 16 SbisOBPs based on the S. bisonia genome. To screen for candidate genes potentially linked to insecticide resistance, we adopted a multi-criteria screening strategy that integrated phylogenetic analysis, molecular docking with three insecticides, and tissue-specific expression profiling. Phylogenetic analysis identified several SbisCSPs and SbisOBPs clustering with genes known to be involved in insecticide resistance, serving as an initial evolutionary filter. Molecular docking results indicated that &#x3bb;-Cyhalothrin exhibited the strong predicted binding affinity with most of SbisCSPs and SbisOBPs. Subsequent qPCR validation of seven prioritized candidates revealed distinct expression patterns: SbisCSP22 was highly expressed in adults and demonstrated strong binding affinity to all three insecticides tested, suggesting a potential role in mediating multi-insecticide response. Conversely, SbisCSP17 was significantly upregulated in larvae, clustered with genes known to mediate imidacloprid resistance, and exhibited strong binding affinity to imidacloprid. Given its larval-specific expression and the soil-dwelling behavior of larvae, we hypothesize that SbisCSP17 is a key candidate gene for larvae coping with soil-treated insecticides.

Animals↗
PubMed2026

Genome-wide SNP data support species boundaries in sympatric Polylepis Ruiz & Pav. (Rosaceae) species from Bolivia and Ecuador.

Species delimitation in the South American genus Polylepis is notoriously challenging due to high morphological similarity and phenotypic plasticity, likely driven by hybridization and gene flow. Previous phylogenetic studies suggested that genetic structure aligns more strongly with geography than with taxonomy, questioning existing species concepts and hampering conservation efforts. We used double-digest RAD sequencing (ddRADseq) to generate genome-wide SNP data for 11 Polylepis species sampled across multiple localities in Bolivia and Ecuador. Population genetic analyses, phylogenetic inference, and network approaches were combined to assess whether genetic structure aligns more closely with taxonomy or geography. Morphologically defined species formed largely cohesive genetic lineages across regions, with species identity explaining substantially more genetic variation than locality. While localized admixture and reticulation were detected among closely related taxa, widespread species showed strong genetic cohesion and clear separation from congeners. Our results indicate that the sampled Polylepis species from Bolivia and Ecuador maintain distinct genetic identities despite localized signals consistent with gene flow. This genome-wide support for current taxonomy highlights Polylepis as a valuable model for studying speciation under gene flow and indicates that multiple geographic sampling will be essential in reconstructing a robust phylogeny of the genus, with important implications for conservation planning in Andean montane forests.

Bolivia↗
PubMed2026

First identification and molecular subtyping of Blastocystis spp. in donkeys in Aksaray province, T&#xfc;rkiye.

Blastocystis is a common intestinal protist worldwide that can infect humans and animals. Although its molecular epidemiology in T&#xfc;rkiye is mostly focused primarily on humans and livestock, equids have received limited attention despite their traditional roles and frequent contact with humans and other animals in rural environments. This study aimed to determine the molecular prevalence and subtype (ST) distribution of Blastocystis spp. in donkeys in Aksaray Province, providing the first molecular data on donkeys in T&#xfc;rkiye. A total of 182 fresh fecal samples were collected from donkeys in nine villages within Aksaray province. Genomic DNA was extracted, and the small subunit ribosomal RNA (SSU rRNA) gene fragment of Blastocystis spp. was amplified via PCR analysis. Positive isolates were sequenced bidirectionally for identification and subsequent phylogenetic analysis of Blastocystis in donkeys. The overall molecular prevalence of Blastocystis spp. in donkeys was 4.4% (8/182). The infection rate was higher in young donkeys (under 3&#xa0;years old; 8.33%) than in adults (3&#xa0;years or older; 2.46%). However, this difference was not statistically significant. Sequence analysis of the positive PCR products revealed the presence of one known livestock-specific subtype, ST10. Phylogenetic analysis showed that the ST10 isolates characterized in this study clustered with isolates identified from different hosts. This study provides the first molecular data on Blastocystis presence in donkeys in T&#xfc;rkiye. The exclusive detection of ST10 suggests potential cross-species transmission, likely facilitated by the traditional practice of co-housing donkeys with other animals in confined barns. These findings indicate that donkeys may contribute to Blastocystis transmission, underscoring the importance of a "One Health" approach in future epidemiological surveillance.

Animals↗
PubMed2026

Bioprospecting microbial genomes to expand the biocatalytic toolbox of rubber oxygenases.

A set of rubber oxygenases was discovered through phylogenetic analysis and AI-based structural modeling of complexes of the putative enzymes with a substrate mimicking cis-1,4-polyisoprene. Sixteen candidate proteins were selected from thermophilic microorganisms, all sequence-related to the Latex clearing protein from Streptomyces sp. K30 (LcpK30). Sequence truncation and solubility tags were then evaluated to enhance protein expression, with the SUMO tag proving to be the most effective. Including LcpK30, nine heme-containing oxygenases were successfully expressed in E. coli NEB 10-beta cells, purified (35-157 mg L-1 yield) and characterized. Steady-state kinetics revealed significant rubber latex-degrading properties for six of them, with the truncated SUMO-fused LcpK30 (SUMO-LcpK30T) showing activity in agreement with literature. Notably, the catalytic efficiencies of all the expressed homologs lay within one order of magnitude and the oxygenase from Thermomonospora echinospora was found to be particularly promising in terms of activity, especially at high latex concentrations (more than 1% w/v). The analysis of reaction mixtures by both HPLC and HPLC-MS confirmed the oxidation of cis-1,4-polyisoprene to form the expected isoprenoid oligomers (n&#x202f;=&#x202f;2-12), whose distribution was consistent with the usual endo-type cleavage pattern in all but one case. This bioprospecting effort afforded a platform of new rubber-degrading enzymes with diverse efficiencies and product profiles, capable of adapting to targeted applications.

Oxygenases↗
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

Identification and characterization of G protein-coupled receptors in the nocturnal halictid bee Megalopta genalis.

G protein-coupled receptors (GPCRs) are one of the largest families of membrane proteins in insects, regulating vision, neural signal transduction, and various physiological behaviors. Megalopta genalis exhibits a unique facultatively eusocial lifestyle and possesses adaptations for nocturnal activity; however, its GPCR family has not yet been systematically characterized. In this study, we performed genome-wide identification, phylogenetic analysis, and expression profiling of GPCRs in M. genalis by integrating genomic annotation and transcriptomic analysis. The results showed that a total of 99 GPCRs were identified in the genome of M. genalis, which were classified into four major families. Here, we show that M. genalis has undergone lineage-specific GPCR repertoire remodeling, marked by the expansion of novel orphan receptors and the systematic loss of multiple receptor subtypes, such as the neuropeptide receptors MIP-R and NPFR. Moreover, opsins have formed a diverse array of combinations and non-GPCR odorant receptors have undergone significant expansion via tandem duplication. Together, these features may represent part of the molecular repertoire associated with the adaptation of M. genalis to a nocturnal lifestyle. Furthermore, transcriptomic analysis revealed distinct spatiotemporal expression divergence within each of the Mth/Mthl and Fz GPCR families, suggesting functional specialization across development and adult tissues. This study provides the first systematic identification and initial functional characterization of GPCRs in M. genalis, revealing an evolutionary pattern characterized by the coexistence of contraction and expansion within the GPCR family. These findings lay a foundation for further studies aimed at elucidating the roles of these GPCRs in regulating M. genalis physiology and behavior.

Animals↗
PubMed2026

Comparison of paralog identification methods and their impact on species tree topologies in target capture phylogenomics within the Sindora clade (Detarioideae: Leguminosae).

Target capture is a common method of generating high throughput DNA sequencing data for phylogenetic reconstruction of species relationships, for which single copy genes are usually most informative. However, a pervasive problem with target capture is that putatively single copy genes may in fact be paralogs resulting from gene duplication, which are problematic for phylogenetic inference because their evolutionary history may differ from the divergence history of species. Here, we use as a case study a target enrichment dataset of 88 species of Detarioideae (Leguminosae) with a focus on the Sindora clade to examine approaches for handling paralogs, including the built-in paralog handling functions in HybPiper and CAPTUS, plus subsequent steps using Putative Paralog Detection and the tree-based Yang & Smith orthology inference approach. We compare the paralogs flagged using these methods and verify their performance with BLAST mapping against a reference genome sequence of Sindora glabra, and then subsequently compare the species tree topologies produced across these methods. Our comparisons of paralogs flagged across the Sindora clade show that the Putative Paralog Detection pipeline was the most accurate in identifying paralogs in terms of its similarity to the BLAST mapping, followed by the built-in paralog identification function of CAPTUS. However, the results we recovered for the Detarioideae subfamily suggest that the largest differences in species tree topology resulted from the use of paralog-filtered alignments (such as with the Putative Paralog Detection pipeline and the Yang & Smith orthology inference approaches) rather than just by removing the sequences of identified paralogous genes. This was the true for HybPiper-assembled datasets but was not seen in CAPTUS-assembled datasets. In all comparisons, the topological differences caused by different paralog handling methods tended to be confined to clades where processes such as hybridisation and introgression are prevalent. Our study provides a roadmap to establish the best approach to identify, eliminate or separate paralogs in the absence of a chromosomally contiguous reference genome for a study group, and highlights the importance of careful data inspection and processing in addition to understanding the extent of paralogy and paralog characteristics (e.g. sequence divergence between copies) for their study group.

Phylogeny↗
PubMed2026

Evolutionary diversification of invertase paralogs couples carbon metabolism and sexual reproduction in fission yeasts.

Dynamic patterns of gene gain and loss play a major role in the diversification of eukaryotes, reflecting adaptation to a broad range of ecological contexts. Reconstructing the evolutionary history of genes provides a powerful framework for understanding how functional innovation shapes life-history traits. Here we report a comprehensive analysis of gene gain and loss across the fission yeast clade, whose evolutionary trajectory remains elusive. Reductive evolution of metabolic genes is a major contributor to species diversification, as observed in other fungal taxa. Notably, we uncovered an evolutionary scenario in which an ancestral gene duplication was followed by lineage-specific loss of one or the other paralog, except in S. pombe, which retained both. We demonstrate that these paralogs encode catalytically-active invertases, named Inv1 and Inv2, with distinct enzymatic properties, localization, regulation, and physiological roles. Inv1 is a secreted enzyme subject to glucose catabolite repression and is the sole invertase required for sucrose assimilation, resembling canonical yeast invertases. In contrast, Inv2 is intracellular, constitutively expressed, and required for inducing sexual differentiation in response to nutrient availability. Overall, these findings reveal an unexpected role for carbon metabolism in modulating the haploid-diploid cycle of fission yeasts, suggesting that diversification of core metabolic functions may contribute to adaptation to environments with distinct sugar compositions.

Evolution↗
PubMed2026

Uce-based phylogeny and classification of Megachilini.

The generic-level classification of the bee tribe Megachilini (Megachilidae) has remained controversial due to poor phylogenetic resolution at the base of the group, particularly among the brood parasitic genera and the numerous dauber ("Chalicodoma s. l.") lineages. We present a phylogenomic analysis of Megachilini based on ultraconserved elements (UCEs), sampling 52 ingroup taxa with emphasis on the dauber lineages. We also present a combined UCE&#xa0;+&#xa0;six-gene analysis to improve taxon coverage, resulting in a dataset with 127 ingroup taxa. Maximum likelihood, coalescent, and Bayesian analyses of multiple UCE matrices recover largely congruent topologies with substantially improved support relative to previous studies. Our results strongly support the monophyly of Megachilini, the early divergence of Noteriades and Gronoceras, and a single origin of brood parasitism. All remaining non-parasitic Megachilini form a moderately supported clade sister to the brood parasitic lineage. The leafcutter bees are monophyletic and nested within dauber lineages. Several major dauber clades are consistently recovered, including an exclusively Australian clade corresponding to the Hackeriapis group of subgenera, while several recognized subgenera are paraphyletic. The lineage known as Morphella, previously placed in synonymy with the subgenus Callomegachile, was not closely related to that subgenus and is here treated as a valid subgenus. Divergence-time analyses place the crown age of Megachilini in the late Eocene to early Oligocene, with major extant lineages diversifying during the Miocene. Limited morphological diagnosability of several clades indicates that splitting non-parasitic lineages into numerous genera would result in an impractical classification that would widen the gap between taxonomists and non-specialists and exacerbate the taxonomic impediment in bees. We therefore advocate retaining a single genus Megachile for non-parasitic Megachilini (excluding Noteriades and Gronoceras), as the classification best supported by phylogenomic evidence and most robust to future taxon sampling.

Animals↗
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↗