Search PubMedSearch

SEARCH · Search PubMed

Results for “chordate evolution”

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.

93 records · Page 3Linked to original sources

Evolutionary architecture and lineage-specific diversification of Forkhead box transcription factors in Perna viridis.

The Forkhead box (Fox) transcription factors are evolutionarily conserved regulators of development, cell cycle, and apoptosis across metazoans. This study provides the first comprehensive genome-wide analysis of the Fox gene family in the Asian green mussel (Perna viridis). We identified 28 Fox genes distributed across 10 chromosomes. Comparative analysis reveals the absence of the FoxI, FoxQ1, FoxR and FoxS subfamily, consistent with other bivalves and indicative of lineage-specific gene loss during molluscan evolution. Notably, gene duplications in the FoxAB, FoxD, FoxH, FoxN1-4, FoxQ2 and FoxQD subfamilies may reflect functional diversification associated with environmental adaptation. Exon-intron structural variability, including intron loss in several paralogues, suggests structural diversification and potential regulatory variation. Phylogenetic reconstruction confirmed the monophyly of core Fox classes while highlighting divergent expansion patterns in lophotrochozoans. Selection analyses showed strong purifying selection across duplicated Fox paralogs, supporting functional conservation after lineage-specific expansion. Gene Ontology enrichment linked Fox genes to stress response, apoptosis, and transcriptional regulation. By integrating phylogenetic, structural, and transcriptomic analyses, this study provides a genomic framework for understanding Fox gene organisation, evolution, and tissue-associated expression patterns in Perna viridis and establishes a comparative resource for future functional studies in bivalves.

Animals

Alternative genetic codes in bacteria and archaea identified with a fast k-mer-based algorithm.

The genetic code is conserved across all domains of life and is often described as universal. Nevertheless, many exceptions to the "universal" code have now been documented, most of these through manual or semiautomated inspection of highly conserved genes. Modern bioinformatics tools improved our ability to find alternative genetic codes but remain computationally expensive, preventing widespread use on thousands of new species identified by sequencing environmental samples. Here, I report a >100-fold accelerated method for inferring the genetic code directly from assembled genomes and apply it to thousands of previously uncharacterized assemblies from archaea and bacteria. I describe three candidate genetic code variations, one of which, an alternative genetic code used by a family of Asgard archaea, is a unique example of sense codon reassignments for this domain. Identifying genetic code variations is important for understanding evolution of the standard code and improving accuracy of protein databases and open reading frame identification.

Genetic Code

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

Building phenotypic character matrices for phylogenetic inference: exploration of 35 years of practice.

Recent methodological development in phylogenetic inference has focused predominantly on molecular data. However, renewed interest in other data types, particularly morphological data, has followed from the increased recognition of the power of total evidence and tip-dating approaches, including fossil data, for inference of time-scaled trees and rates of evolution. However, attention has largely focused on the improvement of models of morphological evolution and other analytical tools with much less discussion about data acquisition itself. Here we review past and current practice for describing and collecting morphological data for phylogenetic inference. We present a systematic review of 164 phylogenetic analyses conducted over the last 35 years and focused on a diverse group of extinct arthropods: trilobites. Trends in increasing matrix size, data type, and coding strategy are evident. Where present, polymorphic characters have been predominantly derived from discretized continuous characters, although increasingly practitioners are utilizing alternative approaches for the treatment of quantitative characters. Not surprisingly, traditional indices that describe character consistency are highly correlated with matrix size but show surprising variation at different taxonomic scales. More recent attempts to describe data quality using information theory imply that characters can have high information content even if data are missing for many tips, providing support against the exclusion of characters because of missing data. In consideration of this, as well as advances in the study of developmental biology and variational complexity, we identify several avenues for increasing the quality and quantity of morphological data going forward.

Phylogeny

First insights into the role of evolutionary history in shaping venom composition of Vipera ammodytes.

Understanding intraspecific venom variation requires distinguishing the contributions of neutral population history from natural selection. This study aims to determine whether venom variation in Vipera ammodytes species complex is structured across eight phylogenetic groups. Despite a complex evolutionary history, venom composition did not differ among phylogenetic groups within the analytical framework used, suggesting that shared ancestry alone does not explain venom variation. Whether local adaptation to environmental conditions explains the observed variation remains an open question for future studies.

Animals

Comparative genomics and full-length transcriptome profiling of wing morphs in Tetrix grossus (Orthoptera: Tetrigidae).

Wing polymorphism represents a paradigmatic dispersal-reproduction trade-off, yet its molecular basis remains uncharacterised in the phylogenetically distant pygmy grasshoppers (Tetrigidae). Here we integrate comparative genomics across ten orthopteran species with full-length transcriptomics of long-winged (FL) and short-winged (FS) Tetrix grossus. OrthoFinder recovered 118 orthogroups specific to T. grossus. Against a backdrop of pronounced gene-family contraction (36 expansions versus 222 contractions; net -186, mirrored at the ancestral Tetrix node, +37/-140), we identified an ancestral, Tetrix-specific expansion of hormone-regulation (12 genes; fold enrichment 7.93) and lipid/carbohydrate-metabolic families organised into syntenic clusters, alongside 513 positively selected genes enriched for integrin-mediated cell adhesion (6 genes), a process relevant to epithelial and appendage morphogenesis. Full-length transcriptomics of one long-winged (FL) and one short-winged (FS) adult female detected 7530 (FL) and 7515 (FS) expressed genes, with 794 FL- and 776 FS-restricted transcriptome-derived SNP-associated genes. The FL morph was enriched for an EGFR/Ras-Rho developmental-patterning axis and neuromuscular flight genes, whereas the FS morph was enriched for insulin/peptide-hormone response and growth-regulatory loci. Overall, we present genomic resources and testable hypotheses concerning the evolution and regulation of wing morphs in Tetrigidae rather than a validated genetic architecture of wing-morph determination.

Animals

Sexual selection purges mutation load, but not overall genetic diversity, decreasing vulnerability to extinction.

Theory suggests sexual selection will enhance population viability by purging deleterious alleles. However, direct genomic evidence for this fundamental idea is scarce and contradictory. We combined long-term experimental evolution with whole-genome resequencing to directly test how sexual selection affects mutation load, genomic divergence, and extinction risk in small populations (maximum Ne = 40) of Tribolium castaneum. After 156 generations, populations evolving under strong sexual selection carried substantially fewer deleterious alleles than populations under weak sexual selection, based on both individual-level estimates of missense and nonsense variants and population-level Rxy analyses, indicating more efficient purging of deleterious alleles. In contrast, nucleotide diversity and runs of homozygosity were similar across treatments, indicating that purging acted most strongly on deleterious variation, and that reduced mutation load in these small populations under strong sexual selection was not explained by demographic effects. Importantly, population-level mutation load estimates best explained extinction risk under inbreeding, directly linking sexual selection to purging and population viability. Genome scans of high and low sexual selection populations revealed peaks of divergence, which included genes involved in courtship, sex discrimination, and seminal fluid proteins. Our results provide direct genomic evidence that sexual selection can reduce mutation load without eroding standing genetic diversity and thus adaptive potential, while driving adaptive divergence in reproductive traits. This beneficial purging may help explain the widespread prevalence of sexual reproduction in nature despite inherent costs and have important ramifications as to how we manage populations of conservation concern.

Animals

The genomic origins and evolutionary path to a key innovation in the world's most venomous snakes.

Evolutionary innovation is a catalyst for the colonization of new environments and the adaptive radiations of major groups. Novel traits typically evolve through the modification of preexisting characters, but the genetic paths underlying their origin have been challenging to trace, and the general requirements for and relative order of different kinds of gene mutations have been difficult to assess. Here, we trace the genomic origins of four procoagulant venom toxins (factor X, factor V, group I phospholipase A2, and Kunitz-type toxins) that collectively underlie a novel, especially potent blood-clotting venom type in the recently evolved Australian brown snake and taipan clade. We find evidence for a previously unknown fifth toxin, coagulation factor VII, and show that the toxins evolved through two distinct genetic paths. The factor X and factor V toxins evolved through the sequential de novo co-option of ancestral clotting factor proteins that entailed their heterotopic expression in the venom gland, the fixation of segmental duplications containing each locus, and subsequent gain-of-function mutations that rendered factor X and factor V constitutively active. In contrast, the phospholipase A2 and Kunitz-type toxins evolved by modifying the functions of neurotoxins that were part of the venom arsenal. Our findings support models in which innovative mutations in single-copy genes precede gene duplication in the evolution of novel proteins and offer a rare view into the genesis of a complex trait that has played a central role in a major adaptive radiation.

Animals

Assembly and Characterization of the First Complete Mitochondrial Genome of Tussilago farfara L.: Insights into Biological Functions and Phylogenetic Relationships within the Asteraceae Family.

Tussilago farfara L., a member of the Asteraceae family, is an economically valuable species due to its edible and medicinal properties. To elucidate the structural characteristics, genetic mechanisms, and evolutionary pathways of the organelle genomes of T. farfara, we sequenced, assembled, and annotated its mitochondrial genome for the first time. The complete mitochondrial genome of T. farfara spans 306,024 bp and contains 33 mitochondrial protein-coding genes (PCGs), 3 rRNAs, and 22 tRNAs. Analysis of the nucleotide substitution rate and genetic diversity revealed that most mitochondrial genome genes may have undergone purifying selection, indicating a slow evolutionary rate and a relatively conserved genomic structure. We further identified 13 fragments of chloroplast-derived DNA integrated into the mitochondrial genome, evidencing intracellular gene transfer. Collinearity analysis showed that Arctium lappa shares the most extensive mitochondrial homologous sequences and the highest sequence similarity with T. farfara. Phylogenetic analysis based on the mitochondrial genome helped to clarify the evolutionary and taxonomic position of T. farfara within the Asteraceae family. The mitochondrial genome sequence of T. farfara provides a valuable genomic resource for species identification and for evolutionary studies within the Asteraceae family.

Genome, Mitochondrial

High levels of mitotic gene conversion are needed to effectively purge deleterious mutations in asexual organisms.

Self-fertilisation and asexual reproduction are both hypothesised to cause long-term extinction due to inefficient selection against deleterious mutations. Self-fertilisation can counter these effects through creating homozygous genotypes and purging deleterious mutations. Although complete asexuality lacks meiotic gene exchange, mitotic gene conversion creates homozygous regions that could limit deleterious mutation accumulation in an analogous manner. We compare mutation accumulation in self-fertilising and facultative sexual populations subject to mitotic gene conversion, and quantify the efficacy of purging in the latter. We first show analytically that purging is most effective with high levels of asexuality and gene conversion, and when deleterious mutations are recessive. We further show using simulations that, when mitotic gene conversion becomes sufficiently high in obligate asexuals, there is a reduction in the mutation count and a jump in homozygosity, reflecting purging. However, this mechanism is not necessarily as efficient at purging under high self-fertilisation, and elevated rates of mitotic gene conversion seem to be needed for widespread purging compared to empirical estimates. If gene conversion rates are allowed to evolve, then elevated rates that increase mean fitness can arise, but only if there is sufficient variance in the gene conversion rate. Conversely, if gene conversion rates are already high and rates are not constrained then they will slightly decrease, reducing mean fitness.

Self-fertilisation

Divergent evolutionary strategies in spider venoms: A comparative proteomic profiling of four sympatric species from Yunnan.

Spider venoms comprise complex cocktails of bioactive molecules evolved for predation and defense, representing a valuable resource for biological research and pharmaceutical discovery. In this study, we performed a systematic analysis of venom gland extracts from four common spider species indigenous to Yunnan, China: Agelena limbata, Hippasa lycosina, Lycosa grahami, and Sinopoda pengi. Using an integrated transcriptomic and proteomic targeted profiling approach, we successfully annotated 141 distinct toxins. Comparative analysis revealed significant interspecific heterogeneity, suggesting distinct evolutionary trajectories and "weapon system economics." Both A. limbata and L. grahami exhibited a "peptide-dominant" profile anchored by neurotoxic peptides and isomerases, optimized for rapid chemical paralysis. In contrast, S. pengi displayed a distinct "protein-dominant" signature enriched with high-molecular-weight enzymes and CAP superfamily proteins, likely functioning to facilitate tissue degradation and toxin diffusion. Occupying an intermediate position, H. lycosina demonstrated a hybrid composition. These findings suggest that although these species share the same geographical range, their venom systems have undergone divergent evolutionary adaptations driven by specific ecological niches and hunting strategies. This study represents the first systematic proteomic characterization of these venom components, providing a valuable reservoir of molecular candidates while highlighting the bioinformatic nuances of analyzing whole-gland homogenates.

Animals

Ambrosia beetle invasions are structured by inbreeding, intraspecific hybridisation, and bridgeheads.

When invasive populations establish in regions far from their origin, they may accumulate deleterious mutations that limit population viability and later expansion. Invasions stemming from such bridgehead populations may experience further sequential bottlenecks. However, deleterious mutations can be masked or eliminated when populations outbreed with other lineages. Here, we analyse global invasions of a species complex of persistently inbreeding ambrosia beetles, using genomic data (N=247) from invasive populations in Africa, North America and Australia, and from native populations in Asia. We mostly focus on one species of this complex (Euwallacea fornicatus) which poses a severe threat to tree species worldwide and is rapidly expanding its global range. We uncover a single lineage of this species across California, South Africa, and Western Australia, involving an invasive bridgehead and containing almost no nuclear genetic variation. In South Africa we identify a second lineage that has repeatedly hybridised with the first lineage. Genetic patterns in the native range indicate that such opportunistic outbreeding may be common. Despite lacking nuclear variation, the first lineage contained two CO1 haplotypes that were also observed in every hybrid lineage, pointing to heteroplasmy and possible hybrid origins of this lineage. Native populations had fewer missense mutations than invasive populations, indicating that opportunistic outbreeding may help purge fixed deleterious mutations when local lineage diversity is high. These findings highlight the importance of outbreeding even when inbreeding is common, and they demonstrate the biosecurity threat posed by subsequent gene flow into invasive populations.

Journal Article

Structural and tissue-specific organisation of endocrine Fgf19 and Fgf21 signalling in rainbow trout.

Endocrine fibroblast growth factors (FGF19 subfamily) play a key role in regulating metabolic homeostasis in vertebrates. However, their functional diversification in salmonids remains poorly understood. In this study, we conducted an integrative characterisation of Fgf19 and Fgf21 signalling in rainbow trout (Oncorhynchus mykiss) by combining phylogenetic, structural and expression analyses. Phylogenetic analyses revealed the conservation of single fgf19 and fgf21 genes, despite the extensive expansion of receptors post-Ss4R (salmonid-specific fourth-round whole genome duplication). Structural modelling and molecular dynamics simulations demonstrated the stable interactions of both ligands to multiple Fgfr isoforms, with receptor-specific energetic profiles and conserved core interaction residues. Tissue expression profiling revealed clear differences from mammalian models, such as predominant hepatic fgf19 expression and the absence of hepatic fgf21 under basal conditions. In addition, there were complex and tissue-dependent distributions of fgfr and klotho transcripts. These findings support a receptor-driven diversification model of endocrine Fgf signalling in salmonids, suggesting enhanced endocrine plasticity associated with the retention of receptors following post-genomic duplication. Taken together, our findings provide new insights into the structural and regulatory organisation of endocrine Fgf signalling, as well as its potential role in metabolic regulation in rainbow trout.

Animals

Epitranscriptomic erasers in bivalves: Evolutionary divergence and species-specific transcriptional plasticity of the ALKBH family under acute thermal stress.

The AlkB homolog (ALKBH) family of Fe(II)/α-ketoglutarate-dependent dioxygenases mediates nucleic acid demethylation, thereby governing RNA metabolism and genomic stability. Despite their pivotal roles in epitranscriptomic regulation across vertebrates, the evolutionary dynamics and functional significance of ALKBH proteins in bivalve mollusks remain largely unexplored. Here, we present a comprehensive phylogenomic analysis of 210 ALKBH genes identified across 35 bivalve species. Our analyses reveal a distinct evolutionary trajectory characterized by the lineage-specific loss of ALKBH4 and the restricted distribution of ALKBH5 to the Mytilidae family, contrasting sharply with vertebrate repertoires. Using the noble scallop (Chlamys nobilis) and Pacific oyster (Crassostrea gigas) as model systems, we demonstrate that ALKBH genes exhibit conserved spatiotemporal expression patterns, with pronounced enrichment in gonadal tissues and during metamorphic transitions, implicating these enzymes in gametogenesis and larval development. Furthermore, comparative thermal stress experiments reveal divergent transcriptional plasticity: the subtropical scallop C. nobilis mounts rapid, transient induction of ALKBH1/2/6 under heat shock, whereas the eurythermal oyster C. gigas maintains sustained ALKBH3 expression, potentially underpinning its superior thermal tolerance. Conversely, cold stress elicits bimodal regulation in C. nobilis, with ALKBH1/2 upregulation contrasting with ALKBH6/7/8 suppression. These findings illuminate the functional diversification of bivalve ALKBH genes and their potential utility as molecular biomarkers for assessing developmental competence and thermal resilience in shellfish aquaculture.

Animals

Robust error-minimization in the genetic code across physicochemical metrics and variant codes: A graph-theoretic analysis in GF(2)6.

The standard genetic code reduces the impact of point mutations, but the robustness of this property across physicochemical metrics, naturally occurring variant codes, and codon-reassignment mechanisms remains incompletely quantified. Embedding the 64 codons in GF(2)6 represents the hypercube Q6 as a coordinate-dependent subgraph of the encoding-independent single-nucleotide mutation graph H(3,4), and enables continuous &#x3c1;-interpolation between the two. Under a quartet-pattern shuffle null (n=10,000), the standard code is significantly low-cost across four established, code-independent physicochemical distance metrics with partially overlapping content (Grant ham p=0.0062; Miyata p<0.001; Woese polar requirement p=0.003; Kyte-Doolittle hydropathy p=0.001), and the signal strengthens monotonically as &#x3c1; moves Q6&#x2192;H(3,4). A structure-aware sensitivity analysis under the alignment-derived ProtSub matrix (Jia & Jernigan 2021) yields the most extreme percentile of any measure tested (p=0.0004; all five p-values pass Bonferroni at &#x3b1;=0.05). Across the 27 NCBI translation tables, near-optimality is preserved: 11 of 12 informative-distance variants retain top-5% placement after BH-FDR correction. Natural codon reassignments avoid disrupting codon-family connectivity: under the encoding-independent H(3,4) adjacency, observed events are topology-breaking at relative risk 0.32 versus the candidate landscape (permutation p&#x2264;10-4). The H(3,4) result is stable by construction; the Q6 decomposition is representation-specific and fails to show depletion under 8 of 24 base-to-bit encodings, so we report H(3,4) as the primary test and Q6 as a sensitivity. Event-level conditional-logit modelling shows that topology avoidance and local physicochemical cost provide complementary, only weakly correlated signal (rs=0.15), and that topology adds explanatory value beyond physicochemistry under both Q6 and encoding-independent H(3,4) adjacency. Retrospective reanalysis of nine genome-recoding datasets is consistent with codon-family topology operating as an evolutionary-trajectory constraint distinct from acute engineering fitness. The contribution is the second axis: code evolution is jointly constrained by physicochemical smoothness and codon-family topological integrity, and these two constraints are partly independent.

Codon reassignment

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

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

Animals

Evolutionary characterization and expression profiling of ACC and FASN genes in Chinese mitten crab Eriocheir sinensis.

Acetyl-CoA carboxylase (ACC) and fatty acid synthase (FASN) are rate-limiting enzymes in the fatty acid biosynthetic pathway, yet their evolutionary relationships, sequence features, and expression profiles remain poorly understood in crustaceans, particularly in the economically important Chinese mitten crab (Eriocheir sinensis). Here, we identified and systematically analyzed ACC and FASN genes in E. sinensis using comparative genomic analyses across 43 species. ACC was highly conserved as a single-copy gene in invertebrates, in contrast to the multiple paralogs observed in vertebrates. Similarly, FASN was generally maintained as a single-copy gene across most taxa but exhibited lineage-specific expansion in certain insect groups. Phylogenetic and structural analysis revealed strong conservation of both genes within crustaceans, supported by multiple conserved motifs and canonical functional domains. Expression profiling showed predominant expression in the hepatopancreas and midgut, suggesting their potential involvement in crustacean lipid metabolism. During the molting cycle, ACC and FASN exhibited higher expression levels during stages C and D, suggesting an increased capacity for fatty acid biosynthesis before molting. In addition, dietary lipid levels experiment revealed that ACC and FASN expression responded dynamically to dietary lipid availability, with increased expression at moderate lipid levels but reduced expression under excessive lipid supplementation, indicating a possible adaptive transcriptional response to lipid status. Collectively, this study provides insights into the evolutionary conservation and expression dynamics of ACC and FASN and improves our understanding of lipid metabolic adaptation in crustaceans.

Animals

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