Search PubMedSearch

SEARCH · Search PubMed

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

163 records · Page 6Linked to original sources

A multi-model genome-wide association study identifies genetic variants underlying resistance to Largemouth Bass Ranavirus (LMBV) in Micropterus salmoides.

Largemouth bass (Micropterus salmoides) is an economically important freshwater aquaculture species, yet recurrent outbreaks of Largemouth Bass Ranavirus (LMBV) continue to impair production and cause substantial losses. The genetic basis of host variation in LMBV resistance remains insufficiently characterized. Here, we applied a multi-model genome-wide association study (GWAS) to identify loci associated with resistance following a controlled challenge with the LMBV-23PY strain. Whole-genome resequencing was performed for 146 phenotyped fish, including 72 susceptible and 74 resistant individuals. After stringent quality control, 877,262 high-quality variants were retained and tested using six GWAS models. Across binary survival status and survival time phenotypes, 32 shared suggestive variants were consistently detected across models, representing suggestive loci for LMBV-23PY resistance. Genes within ±50 kb of these loci were annotated, and functional enrichment highlighted immune- and redox-related biological processes. Three prioritized candidates-GSTT3L (glutathione S-transferase theta-3-like), CGRP2 (calcitonin gene-related peptide 2), and NPPC (natriuretic peptide C)-were associated with pathways involved in oxidative stress responses and immune regulation. Collectively, these results provide insight into the genetic architecture of LMBV-23PY resistance in largemouth bass and identify suggestive variants and associated candidate genes for downstream validation, functional interrogation, and the development of marker-assisted and genome-enabled breeding strategies.

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

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

A conserved distal-tail helical extension defines a tailspike attachment architecture in Gram-negative siphophages.

Rapid growth of bacteriophage genome collections has outpaced functional annotation of tail-tip proteins, limiting comparative analysis of host-recognition structures. Starting from a shared distal-tail gene organization in the Salmonella phages 9NA and Jersey, I developed a morphogenetic bioinformatic framework integrating gene synteny, sequence comparison, profile hidden Markov model (HMM) screening, structural evidence, structure-aware searching, and AlphaFold modeling. Comparison with the experimentally characterized lambda and Sf11 tail assemblies identified a predominantly alpha-helical C-terminal extension of the distal-tail (DT) protein associated with tailspike attachment, termed the distal-tail helical extension (DT-helix). Screening 541,986 proteins from 5167 complete NCBI RefSeq tailed-phage genomes, followed by evidence-based evaluation of sequence, genomic context, and structural architecture, identified 165 curated DT-helical-extension-associated phages. Their DT proteins segregated into six sequence groups. In the four principal multi-member groups, cognate tailspikes showed group-specific conservation in proximal N-terminal regions but substantially greater downstream diversity, consistent with sequence constraint at the DT-tailspike attachment boundary. A complementary ProstT5/Foldseek search supported the established groups but revealed no convincing additional highly divergent family. Together with the experimentally characterized Sf11 attachment interface, these findings define a recurrent morphogenetic architecture linking conserved distal-tail scaffolds to more variable receptor-binding proteins across siphophages infecting Gram-negative bacteria. Although universal exchangeability is not established, the identified scaffold-receptor-binding boundaries provide a framework for molecular characterization and rational phage engineering. Accession-level information for the 165 curated phages is available through PhageTailDB.

Viral Tail Proteins

A 4-year longitudinal wastewater surveillance of five gastroenteritis viruses and the correlation with clinical cases in Alberta, Canada.

Viruses are common causes of acute gastroenteritis worldwide. They are detected in large quantities in raw sewage making them amenable to wastewater-based surveillance (WBS). To monitor the prevalence of gastroenteritis viruses in wastewater and assess their correlation with clinical cases, wastewater samples collected between July 2020 and June 2024 from 12 wastewater treatment plants across Alberta, Canada were analyzed for norovirus (NoV) GI & GII, rotavirus (RoV), adenovirus (AdV), sapovirus (SaV) and astrovirus (AsV). Among the 5726 wastewater samples tested, AdV (80.5%) had the highest detection rate followed by NoV GII (75.6%), SaV (63.2%), NoV GI (59.4%), RoV (42.2%) and AsV (20.2%). Winter and spring seasonality was found for NoV and RoV in both wastewater and clinical disease. Public health interventions especially in the 1st year of the COVID-19 pandemic had a significant impact on their burden with marked reduction in wastewater detected viruses and clinical cases. NoV showed a strong correlation between its level in wastewater and the number of clinical cases, while moderate correlation was observed for the other four viruses. Cross-correlation analysis showed that changes of viral RNA concentration in wastewater lagged behind reported gastroenteritis cases by approximately 6 days to 3 weeks. To our knowledge, this is the longest multi-region WBS study monitoring multiple gastroenteritis viruses spanning both COVID-19 pandemic and post-pandemic periods. The data obtained from this study supported WBS as a complementary tool to track population-based circulation of gastroenteritis viruses, providing actionable public health data.

Clinical cases

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

Studies on functional differentiation of xpr1a and xpr1b genes in zebrafish.

Xenotropic and polytropic retrovirus receptor 1 (XPR1) is known to be involved in various biological processes, including phosphate homeostasis, cellular signaling, brain and vascular mineralization, whereas its specific contribution to bone development remains incompletely characterized. Due to genome duplication in teleosts, zebrafish Danio rerio possess two paralogous genes of XPR1 namely xpr1a and xpr1b, whose functional divergence remains unclear. The amino acid sequence similarity between zebrafish xpr1a and xpr1b was 83.26%. In situ hybridization demonstrated overlapping localization in the head and spinal cord at 24-48 hpf, while diverged by 72 hpf, with xpr1a becoming restricted to the head while xpr1b persisted in both regions. CRISPR/Cas9 was used to generate xpr1a and xpr1b mutants. The xpr1a mutants are comparatively healthy, viable but with mild growth reduction, whereas the xpr1b mutants display high mortality, reduced body length and severe vertebral deformities. Interestingly, all the double mutants died at the embryonic stage. Moreover, to further investigate the molecular and regulatory mechanisms, we conducted comparative transcriptome analysis on bone and brain tissues from xpr1b+/+ and xpr1b-/- zebrafish. In bone tissue, 6749 DEGs were identified, comprising 3846 upregulated and 2903 downregulated genes. These DEGs were mainly enriched in the MAPK signaling pathway, Wnt signaling pathway, cysteine and methionine metabolism, and ECM-receptor interaction. RT-qPCR validated results showed that seven osteogenesis-related genes (col1a1a, sp7, runx2b, col1a2, col1a1b, alp1 and entpd5), and two phosphate homeostasis related genes (slc20a2 and pdgfba), which are essential for skeletal mineralization and phosphate homeostasis, exhibited significantly downregulated expression in bone tissue of xpr1b mutant zebrafish. These results highlight the pivotal role of xpr1b in regulating skeletal mineralization and phosphate metabolism, thereby elucidating the functional specialization of XPR1 paralogs while providing a theoretical basis for understanding bone developmental mechanism in teleost vertebrates.

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

Pre-clinical immunogenicity and safety evaluation of H2 strain Hepatitis A Inactivated Vaccine in rhesus macaques.

BACKGROUND: Hepatitis A is a viral infection of the liver that can cause mild to severe illness. Currently, two types of HAV vaccines are used worldwide, inactivated hepatitis A vaccines, which are used in most countries, and live attenuated vaccines (H2 and L-A-1 strain), which are mainly used in China. The major disadvantage of live attenuated virus to cause secondary infections among contacts and mutation shifts of the live vaccine strain. The H2 strain was selected for the development of an inactivated hepatitis A vaccine to further reduce biosafety risks. Rhesus macaques high genomic homology with humans and the incubation period after hepatitis A vaccination and human natural infections are similar. We use rhesus macaques to assess immunogenicity and safety of the H2 strain Hepatitis A Inactivated Vaccine. METHODS: The vaccine was assessed in rhesus macaques, divided into four groups (n = 10 per group): the control group (adjuvant buffer; aluminum content 0.35 mg/mL; 2 mL per dose), the low-dose group (320EU, 0.5 mL of 640EU/mL with aluminum content 0.35 mg/mL), the medium-dose group (640EU, 1 mL of 640EU/mL with aluminum content 0.35 mg/mL), and the high-dose group (1280EU, 2 mL of 640EU/mL with aluminum content 0.35 mg/mL). Animals were injected intramuscularly at multiple sites in the hind limbs and received four inoculations at 4-week intervals. Test items including Clinical indicators, immunogenicity indicators and Histopathological examination. RESULTS: No abnormalities were observed in any group in terms of general clinical condition throughout the study period except for slight decreases in body temperature after immunization. Hematological parameters, serum biochemistry indices, and histopathological findings showed fluctuated to different degrees of fluctuation after immunization across all groups. Immunogenicity assessments showed that the inactivated hepatitis A vaccine (H2) induced both humoral and cellular immune responses effectively, and the levels of antibodies increased with certain dose- and time-response trends. CONCLUSION: The inactivated hepatitis A vaccine (H2 strain, human diploid cell) was safe and immunogenic in non-human primates. The results provide strong preclinical support for the further clinical development of this vaccine candidate.

Animals

Efficacy of NAs in combination with Peg-IFN for functional cure in patients with CHB: a meta-analysis of RCTs.

BACKGROUND: The goal of treating chronic hepatitis B (CHB) is to achieve functional cure. We aimed to evaluate the efficacy of the combination of nucleoside (acid) analogues (NAs) and pegylated interferon (Peg-IFN) on the functional cure of CHB patients at the EOT and at the EOF. METHOD: PubMed, Embase, and Web of Science were searched systematically up to March 15, 2025. Sixteen RCTs on CHB patients receiving combination or monotherapy were included. RESULT: Compared with NAs treatment, the NAs combined with Peg-IFN treatment group significantly improved HBsAg clearance rate (RR: 14.05, 95% CI 6.13-32.20) and HBsAg seroconversion rate (RR: 12.82, 95% CI 5.08-32.33) at the EOT. Moreover, compared with NAs treatment, the NAs combined with Peg-IFN treatment group significantly improved HBsAg clearance rate (RR: 7.70, 95% CI 4.24-13.98), HBsAg seroconversion rate (RR: 11.93, 95% CI 5.14-27.67) at the EOF. However, compared with Peg-IFN monotherapy, the combination therapy group did not show any improvement in HBsAg clearance rate, HBsAg seroconversion rate, and the rate of qHBsAg decrease > 1 log10 IU/mL at the EOT and the EOF. Moreover, in terms of safety, the Peg-IFN monotherapy group showed more ALT flares (ALT > 5 × ULN) during the treatment process. CONCLUSION: Compared with NAs therapy, Peg-IFN combined with NAs therapy significantly improves functional cure rate. However, combination therapy shows no additional advantage over Peg-IFN monotherapy in achieving functional cure. The Peg-IFN monotherapy group has a higher incidence of ALT flares than the combination therapy group.

Humans

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