SEARCH · Search PubMed

Results for “experimental 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.

SourceResearch connectionFrom yearThrough year

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

272 records · Page 3Linked to original sources

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

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 + 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

Can't see the forest for the trees: The influence of marker type on inferred phylogenetic relationships in a cosmopolitan bat genus.

Fine-resolution information on species relationships and biological diversity is critically needed to guide conservation efforts amidst rapid environmental changes. Systematics, which forms the foundation of this knowledge, has been revolutionized by phylogenomics, utilizing genome-scale datasets. However, the use of diverse marker types, non-comparable taxon sampling, and outgroup selection can lead to conflicting phylogenetic hypotheses. These inconsistencies complicate study comparisons and hinder our ability to assess marker-specific impacts on phylogenetic resolution. The phylogenetic reconstruction of the bat genus Myotis, encompassing over 140 species and characterized by a rapid radiation in the last 20 million years, has been particularly influenced by these challenges. Achieving phylogenetic resolution in Myotis is particularly complex due to subtle interspecific differences in both morphological and molecular traits. Mitochondrial and nuclear markers often produce discordant trees, influenced by hybridization, introgression, and methodological variations. In this study, we employed a consistent taxonomic sample set of 44 Myotis taxa to evaluate the impact of five different genetic marker types on phylogenetic reconstruction. We observed significant discordance between topologies derived from conserved nuclear and mitochondrial markers and found that transposable elements were inadequate for resolving relationships across the entire genus. Our results also clarify the placement of previously problematic taxa within the genus. These findings emphasize the importance of aligning genetic marker choice with specific phylogenetic questions and highlight the influence of taxonomic and methodological variation on phylogenomic outcomes. This work provides a framework for improving phylogenetic inference in rapidly radiating groups and enhances our understanding of evolutionary history in Myotis.

Animals↗

PubMed2026

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↗

PubMed2026

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↗

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

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↗

PubMed2026

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↗

PubMed2026

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↗

PubMed2026

Spatial transcriptomics of Ciona adult brains reveals functional zonalization and insights into neural gland function.

The ascidian Ciona is a pivotal chordate model for illuminating the evolutionary origins of the vertebrate brain. Here, spatial transcriptomics of the adult Ciona neural complex, combined with image-based computational super-resolution mapping, resolved distinct tissue domains including the cerebral ganglion, neural gland, ciliated funnel, neural gland duct/dorsal strand, and body wall muscle. Within the cerebral ganglion, high-resolution mapping revealed clear molecular zonalization separating the cortex and medulla, alongside regional specialization within the cortex itself. The neural gland exhibited localized enrichment of genes associated with extracellular matrix and cell-cell interactions. These spatial features suggest that the neural gland functions as a homeostatic and signaling interface, reminiscent of primitive vertebrate meninges or choroid plexus. Overall, this spatially defined gene expression map provides a foundational framework for understanding functional regionalization in the tunicate brain and its evolutionary relationship to vertebrate nervous systems.

Ciona↗

PubMed2026

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↗

PubMed2026

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↗

PubMed2026

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↗

PubMed2026

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↗

PubMed2026

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↗

PubMed2026

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↗

PubMed2026

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↗

PubMed2026

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 ρ-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 ρ moves Q6→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 α=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≤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↗