Search PubMed⌕ Search

SEARCH · Search PubMed

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

At least 19 recordsLinked to original sources
PubMed2025

Molecular evolution of the members of the Snq2/Pdr18 subfamily of Pdr transporters in the Hemiascomycete yeasts.

The transporters of the ATP-Binding Cassette (ABC) Superfamily involved in the Multidrug Resistance (MDR) phenomena are also known as ABC-Pleiotropic Drug Resistance (PDR) proteins. The homologs of the Saccharomyces cerevisiae SNQ2 and PDR18 genes were identified in 171 yeast genomes, representing 68 different hemiascomycetous species. All early-divergent yeast species analyzed in this work lack Snq2/Pdr18 homologs, suggesting that the origin of these ABC-PDR genes in hemiascomycete yeasts resulted from a horizontal transfer event. The evolutionary pathway of the Snq2/Pdr18 protein subfamily in pathogenic Candida species was also reconstructed, revealing a main gene lineage leading to the Candida albicans SNQ2 gene. The results indicate that, after the gene duplication event at the origin of the SNQ2/PDR18 paralogs, the PDR18 ortholog has been under strong diversifying selection and suggest that a small portion of the sequence of the SNQ2 ancestral ortholog might have been under mild positive selection. The results also showed that strong positive selection was exerted over one of the two paralogs generated by the Whole Genome Duplication (WGD) event, corresponding to the duplicate at the origin of a "short-lived" WGD sublineage.

Evolution, Molecular↗
PubMed2026

Mapping the Molecular Evolution and Role of Wild Rice GLYIII Protein-Encoding Genes in Abiotic Stress Response.

To address the need for sustainable food production amid rapid global climate change, developing rice varieties that grow optimally even under harsh conditions is essential. An effective approach in this direction would be to harness the stress resilience traits of the crop wild relatives (CWRs) of rice. Among the various crucial stress-responsive genes, the Glyoxalase III (GLYIII) gene family is of utmost importance for its ability to detoxify the toxic glycolytic byproduct, methylglyoxal (MG), in a less energy-intensive, single-step process, as well as for its multifaceted cytoprotective role. In our study, a comprehensive genome-wide search across the Oryza genus revealed that GLYIII genes are conserved across wild rice genotypes. Their number has expanded during domestication, driven by gene duplications. Interestingly, only a few orthologous pairs showed positive selection, suggesting that the functions of most others need to be constrained and or conserved.We found that higher GLYIII activity, Total Antioxidant Capacity, endogenous glutathione (GSH) levels, and free radical scavenging activity contributes to the stress resilience of wild rices O. punctata, O. meridionalis, and O. nivara, in addition to other factors. , , . , . Our qRT-PCR analysis revealed differential expression of the OpGLYIII, OmGLYIII, and OnGLYIII genes across different developmental stages and in response to various abiotic stresses. Furthermore, we report that wild rice GLYIII proteins, specifically OpGLYIII-3, OmGLYIII-3, and OnGLYIII-5, exhibit high catalytic efficiency over a broad pH range and at higher temperatures under in vitro assay conditions. Overexpression of these proteins was found to impart substantial stress resilience to the transformed E. coli cells. These findings collectively suggest that GLYIII proteins constitute a key component of the abiotic stress response machinery in wild rice.

Oryza↗
PubMed2026

Advances in CRISPR Base Editing: From Molecular Evolution to Therapeutic Applications in Genomic Medicine.

CRISPR-Cas9 systems revolutionized gene editing, but inherent drawbacks, namely DNA double-strand breaks (DSBs) and the difficulty of achieving precise repairs (due to low HDR efficiency), led researchers to invent new, more accurate gene editing tools. Base editing represents a significant leap forward, enabling targeted single-nucleotide conversions directly on the DNA without DSBs or donor templates. The core technology involves fusing catalytically dead or nickase Cas proteins to DNA deaminase enzymes. Cytosine base editors (CBEs) convert C•G to T•A pairs, while adenine base editors (ABEs) change A•T to G•C. These editors exploit the deaminase function within the R-loop structure formed by Cas binding and co-opt endogenous DNA repair mechanisms for precision. While offering improved efficiency and editing precision, base editing faces persistent challenges, such as off-target effects, bystander edits, delivery and ethical concerns. Continuous engineering efforts have refined these tools, enhancing accuracy, expanding targetability and reducing unwanted edits. The base editing arsenal has also broadened to include C-to-G base editors (CGBEs), dual A&C editors and versions targeting organelles. Successful preclinical studies demonstrating the correction of mutations responsible for the disease have paved the way for clinical trials, which are now testing therapies for conditions like sickle cell disease, β-thalassaemia and hypercholesterolemia using various delivery systems. This review explores CRISPR base editing's origins, mechanisms of action, potential therapies and current restrictions, pointing to its broadening impact on medical genetics.

Humans↗
PubMed2026

Comparative Analysis of Chloroplast Genomes Reveals Molecular Evolution and Phylogenetic Relationships in Fraxinus (Fraxinus mandshurica).

Fraxinus mandshurica (Manchurian ash) is an ecologically and economically valuable hardwood tree native to Northeast Asia, yet its genomic resources remain limited. We assembled its complete chloroplast (cp) genome (155,559 bp) using hybrid PacBio and Illumina sequencing and performed comparative, phylogenetic, and evolutionary analyses. The cp genome exhibits a typical quadripartite structure encoding 132 gene copies, comprising 114 unique genes (80 protein-coding, 30 tRNA, and 4 rRNA genes), with 18 genes duplicated in the inverted repeat (IR) regions. Simple sequence repeat analysis revealed dominance of mononucleotide A/T repeats. Phylogenetic analysis of 53 complete cp genomes strongly supported the monophyly of Oleaceae and resolved F. mandshurica as sister to the North American F. nigra, consistent with previously proposed Miocene intercontinental dispersal scenarios between East Asia and North America. Most protein-coding genes were under strong purifying selection (Ka/Ks << 1), whereas petB, rpl2, and several ndh genes showed elevated Ka/Ks values that are suggestive of altered selective constraint but are based on very few substitutions and are therefore not, on their own, evidence of positive selection. Nucleotide diversity (Pi) analysis identified 15 hypervariable intergenic spacers (mean Pi = 0.067), among which trnM-CAU-rps14, ndhJ-ndhK, and petL-petG represent promising candidate barcode regions requiring further validation. This study provides a high-quality, fully annotated cp genome of F. mandshurica and a valuable genomic resource for future phylogenetic, population genetic, and conservation studies of this important genus.

Fraxinus↗
PubMed2026

Spatial proteomics reveals four-stage molecular evolution in cancer immunotherapy-related gastritis.

BACKGROUND: Immune checkpoint inhibitors (ICIs) have transformed cancer treatment, yet immune-related adverse events (irAEs) including immunotherapy-related gastritis (IRAEG) pose significant clinical challenges-often necessitating treatment interruption that may compromise antitumor efficacy. IRAEG presents with atypical symptoms, lacks specific biomarkers, and shows histopathological overlap with other forms of gastritis, complicating diagnosis and management. Despite increasing clinical recognition, a systematic understanding of spatial molecular alterations across the full disease course remains limited. Here, we used spatial proteomics to map the molecular landscape of IRAEG during disease progression and to define stage-specific patterns of molecular evolution relevant to cancer immunotherapy management. METHODS: We analyzed tissue samples from seven patients, including four non-immunotherapy-related gastritis controls and three cancer patients who developed IRAEG following ICI therapy for solid tumors, sampled longitudinally across four disease stages: baseline (G1), acute severe inflammation (G2), early recovery (G3), and complete recovery (G4). Using laser capture microdissection coupled with data-independent acquisition mass spectrometry, we profiled 177 spatially resolved gastric tissue regions. Multiplex immunohistochemistry and immunofluorescence characterized features of the immune microenvironment, while Gene Ontology, KEGG pathway analysis, Gene Set Variation Analysis, and xCell inference enabled functional, metabolic, and immune profiling. Key immune and NET-related findings were further validated by multiplex immunofluorescence in an independent, expanded cohort of IRAEG and non-immunotherapy-related gastritis samples. RESULTS: IRAEG was characterized by widespread HLA molecule activation and enhanced antigen processing, resembling the immune phenotype observed in organ transplant rejection. The acute G2 stage exhibited excessive neutrophil extracellular trap formation, profound metabolic suppression, and collapse of immune homeostasis-features that may inform early intervention strategies to preserve ICI treatment continuity. During early recovery (G3), inflammatory injury transitioned toward repair, marked by activation of fatty acid metabolism and PPAR signaling. Notably, even at complete clinical recovery (G4), more than 1,000 proteins remained differentially expressed, reflecting sustained enhancement of metabolic and immune functions and establishing a distinct molecular "memory" state with implications for ICI rechallenge decisions. CONCLUSIONS: These findings define four molecularly distinct stages of IRAEG progression and recovery. The stage-specific signatures identified here serve as candidate biomarkers for diagnosis, disease staging, and therapeutic response assessment, and may guide clinical decisions regarding irAE management, treatment modification, and safe ICI rechallenge to support continued antitumor therapy.

Humans↗
PubMed2025

Molecular Evolution and Zoonotic Potential of Muju Virus (Orthohantavirus puumalaense) in Craseomys regulus, Republic of Korea.

Orthohantavirus puumalaense causes hemorrhagic fever with renal syndrome in Europe, with Puumala virus (PUUV) as its primary representative. Muju virus (MUJV), harbored by Craseomys regulus, an Arvicolinae rodent species endemic to the Republic of Korea (ROK), is also a genotype of O. puumalaense. However, their genomic diversity and zoonotic potential remain largely unknown. To investigate their prevalence, 185 voles were collected from 23 regions of the ROK between 2012 and 2023. Serological assays detected anti-PUUV immunoglobulin G antibodies in five samples (3.1%), whereas reverse-transcription polymerase chain reaction confirmed MUJV RNA in identical specimens (2.7%). Amplicon-based nanopore sequencing facilitates near-complete genome recovery, enabling high-resolution comparative analysis. Phylogenetic analysis revealed distinct genetic lineages in Gangwon and Jeollabuk Provinces. Evolutionary rate estimates indicated greater sequence divergence in the S and L segments than in the M segment. A zoonotic risk assessment revealed that most MUJV variants exhibited moderate-to-high spillover potential. The molecular detection of MUJV in Cheorwon, Gangwon Province, expands its known geographic range and provides the first molecular evidence of MUJV circulation in this region. These findings highlight the need for continued surveillance and seroprevalence studies of MUJV to assess its potential for human exposure and public health relevance in the ROK.

Animals↗
PubMed2026

Molecular evolution and immune expression analysis of ELF transcription factors in Lethenteron reissneri.

As important members of the ETS superfamily, the E74-like factor (ELF) transcription factor family regulates gene transcription through a conserved ETS domain and plays critical roles in immune regulation. However, the evolutionary characteristics and functions of this family in lampreys (Lethenteron reissneri) remain unclear. In this study, the ELF gene family of lampreys (Lr-ELF1, Lr-ELF2, Lr-ELF3, and Lr-ELF5) was systematically identified, and their molecular evolutionary features and immune response functions were investigated. Phylogenetic analysis revealed evolutionary characteristics reflecting the transition from jawless to jawed vertebrates. Domain architecture, gene structure, and three-dimensional structural analyses indicated that these genes appear to be conserved among vertebrates, with their three-dimensional structures showing high similarity to the core secondary structural elements of human homologous proteins. Synteny analysis demonstrated significant differences in the genomic neighborhoods of ELF genes between lampreys and jawed vertebrates. Quantitative real-time PCR (qRT-PCR) was performed with three biological and three technical replicates; relative expression levels were calculated using the &#x394;Ct method, and statistical analysis was carried out with GraphPad Prism 9. Quantitative real-time PCR (qRT-PCR) results suggested that the ELF gene family may be involved in immune defense. This study not only enriches our understanding of the evolution of ELF genes but also provides new insights into the roles of lamprey ELFs in immune defense.

Animals↗
PubMed2026

Nipah virus in the era of global connectivity: molecular evolution, transmission risk, and preparedness strategies.

Nipah virus (NiV) is a highly pathogenic zoonotic RNA virus belonging to the genus Henipavirus within the family Paramyxoviridae, representing a continuing global health concern due to its high case fatality rate and potential for epidemic expansion in the era of increasing international connectivity. The virus demonstrates strong evolutionary adaptability driven by the absence of proofreading mechanisms during RNA replication, enabling genetic diversification that may influence host range, virulence, and transmission dynamics. Molecular pathogenesis of NiV is primarily mediated through interaction of viral glycoproteins with ephrin-B2 and ephrin-B3 receptors, facilitating host cell entry, endothelial damage, and neuroinvasion. Immune evasion facilitated by the action of accessory proteins encoded by the P gene (P, V, W, and C) acts to suppress innate antiviral immunity through the inhibition of interferon induction and JAK/STAT signaling. Human-to-human transmission of Nipah virus remains limited, with epidemiological evidence indicating basic reproduction numbers generally below unity; however, respiratory involvement and healthcare-associated exposure may enhance cluster outbreaks. Global travel, ecological disruption, and fragmented surveillance systems contribute to spillover risk, particularly in South and Southeast Asia where fruit bats of the genus Pteropus serve as natural reservoirs. Despite advances in vaccine technology, including subunit, viral vector, mRNA-based platforms, and monoclonal antibody therapies, no licensed prophylactic or therapeutic agent is currently available for human use. Global preparedness remains challenged by the scarcity of high-containment biosafety facilities, limited research funding, and absence of integrated One Health surveillance networks. Ethical considerations surrounding wildlife population control further complicate disease mitigation strategies. Emerging genomic surveillance, artificial intelligence-assisted predictive modeling, and regional data-sharing frameworks are essential for early detection and response. Strengthening molecular research on viral-host interactions and transmission determinants will be critical for preventing future Nipah virus outbreaks in an increasingly interconnected world.

Genomic surveillance↗
PubMed2025

Phylogenetic diversity and molecular evolution of Hantaan virus harbored by Apodemus chejuensis on Jeju Island, Republic of Korea, 2022-2023.

BACKGROUND: Hantaan virus (HTNV), hosted by Apodemus spp., is a well-recognized causative agent of hemorrhagic fever with renal syndrome (HFRS) and poses a crucial global public health concern. Based on the current evidence, HTNV carried by A. chejuensis is proposed as the likely etiological agent of HFRS on Jeju Island, Republic of Korea (ROK). METHODOLOGY/PRINCIPAL FINDINGS: In this study, 50 small mammals were collected from five locations in Seogwipo-si and Jeju-si on Jeju Island, ROK, during 2022-2023. Serological and molecular analyses revealed HTNV prevalence rates of 34% (16/47) and 27.7% (13/47), respectively. Using a multiplex polymerase chain reaction-based nanopore sequencing approach, nine complete HTNV genomes were sequenced from the lung tissues of A. chejuensis, representing the first comprehensive genomic characterization of HTNV from Seogwipo-si (Hogeun-dong) and Jeju-si (Sangdae-ri). Phylodynamic analyses suggest evolutionary divergence and phylogeographic diversity, with four unique amino acid substitutions identified in HTNV genomes from Seogwipo-si. CONCLUSION/SIGNIFICANCE: This study provides important insights into the genomic surveillance, genetic diversity, and evolutionary dynamics of orthohantaviruses, which are essential for guiding effective public health strategies to control and prevent future HFRS outbreaks in the ROK.

Animals↗
PubMed2026

Molecular Evolution and Expression Analysis of the ADH Gene Family in Apple Bud Mutants.

Alcohol dehydrogenase (ADH) catalyzes the reduction of aldehydes to alcohols, key precursor substrates for volatile ester biosynthesis, which determines the characteristic aroma of apple fruit. However, a comprehensive genome-wide investigation of the ADH gene family in apple has been lacking. In this study, we systematically identified ADH genes in the apple genome using integrated bioinformatics approaches, including phylogenetic analysis, synteny evaluation, promoter cis-element prediction, codon usage bias assessment, and protein interaction network modeling. Expression patterns were examined through transcriptomic data and validated by RT-qPCR analysis across different organs and among 'Red Delicious' and its four bud mutant lines. We identified 44 ADH genes, with 12 forming a prominent cluster on chromosome 1. RT-qPCR analysis revealed that MdADH20 was dramatically upregulated in the 'Red Chief' mutant (relative expression of 59.38), suggesting its pivotal role. Phylogenetic analysis revealed a close evolutionary relationship with wild strawberry. The encoded proteins were generally stable and predominantly localized to the cytoplasm. Promoter analysis showed enrichment of growth/development-related and ARE elements, while codon usage analysis identified AGA, GCU, GUU, and CUU as preferred codons. Protein interaction prediction suggested MdADH19 and MdADH20 as hub proteins. Expression profiling and RT-qPCR further identified MdADH20 as a core candidate gene, characterized by its stable and high expression, particularly in the 'Red Delicious' mutant. Its central position in the predicted protein-protein interaction network suggests a potential regulatory role in the aroma biosynthesis pathway of apple fruit. This study provides the first systematic genome-wide characterization of the apple ADH gene family, establishing a theoretical groundwork for deciphering aroma biosynthesis mechanisms and offering potential target genes for flavor improvement through bud mutation breeding strategies.

ADH gene family↗
PubMed2026

Phylogeography and molecular evolution of Newcastle disease virus across a century of global surveillance.

Newcastle disease virus (NDV) remains one of the most economically important avian pathogens worldwide, causing recurrent outbreaks in poultry despite decades of vaccination and disease control efforts. Since the first reported outbreak of NDV a hundred years ago, numerous molecular epidemiological studies have been conducted globally across diverse geographic and production settings. Following a century of NDV circulation and evolution, the present study aimed to compile all publicly available NDV sequence data and perform a comprehensive global analysis of the genetic diversity, phylogenetic relationship, and global spatiotemporal distribution of NDV over a 100-year timescale. All publicly available NDV complete genome and full-length fusion (F) gene sequences were retrieved from GenBank up to February 2026. Following rigorous quality control, phylogenetic analyses were performed based on complete genomes and F gene datasets. Phylogenetic analysis identified two genotypes within Class I and 20 genotypes within Class II NDVs, with extensive diversification at the sub-genotype level. Genotype XIII exhibited the greatest sub-genotypic diversity, while genotype VII represented the most globally disseminated genotype, reported across 36 countries. Chronological assessment based on the earliest available reports indicated an increasing number of recognized genotypes from the 1930s to recently described sub-genotypes such as XIII.2.3 and XXII.2.2. Regional diversity analysis revealed the highest genotype diversity in Western Africa, Eastern Asia, and Southern Asia. Comparative residue analysis demonstrated substantial genotype-specific variation within critical functional domains of the fusion protein, including cleavage sites, neutralizing epitopes, and heptad repeat regions. Overall, this study provides the first comprehensive 100-year global overview of NDV evolution and phylogeography. The findings highlight continuous viral diversification, broad geographic dissemination of multiple genotypes, and ongoing molecular variation, emphasizing the need for sustained genomic surveillance and periodic evaluation of vaccine compatibility with emerging NDV genotypes.

100-years of data↗
PubMed2026

Gene-level complexity explains genome-wide variation in the distribution of fitness effects.

The distribution of fitness effects (DFE)-describing how harmful, neutral, or beneficial new mutations are-is central to understanding how populations evolve. Although the DFE varies across genomes and species, it remains unclear which aspects of genomic organization drive this variation. Here, we inferred gene-level selective constraints across the genomes of Mus musculus castaneus, Drosophila melanogaster and Saccharomyces cerevisiae using a combination of population genetics and machine learning trained on diverse gene features. Many gene features were predictive of selective constraint, with conservation, gene structure, and expression being the most informative. These selective constraints delineated gene classes with distinct DFEs. Genes with higher connectivity and expression-features reflecting how many traits a gene influences-experienced stronger and less dispersed deleterious effects with increasing selective constraint. Between species, the rate of adaptation decreased with increasing organismal complexity, whereas across the genome it did not decrease monotonically with selective constraint, but tended to be higher at intermediate levels. While between-species comparisons of DFE parameters were less consistent with predictions of Fisher's geometric model (FGM) based on organismal complexity, variation in DFE parameters across the genome aligned more closely with FGM when complexity was considered at the gene level. Our results suggest that gene-level complexity, captured by genomic feature proxies, provides a more informative definition of complexity for DFE variation than organism-level labels, and highlight the value of using gene features collectively to link genomic architecture, fitness landscapes, and patterns of molecular evolution.

Animals↗
PubMed2025

Comparative Phylogenetics Reveal Clade-specific Drivers of Recombination Rate Evolution Across Vertebrates.

Meiotic recombination is an integral cellular process, required for the production of viable gametes. Recombination rate is a fundamental genomic parameter, modulating genomic responses to selection. Our increasingly detailed understanding of its molecular underpinnings raises the prospect that we can gain insight into trait divergence by examining the molecular evolution of recombination genes from a pathway perspective, as in mammals, where protein-coding changes in later stages of the recombination pathway are connected to divergence in intra-clade recombination rate. Here, we leverage increased availability of avian and teleost genomes to reconstruct the evolution of the recombination pathway across two additional vertebrate clades: birds, which have higher and more variable rates of recombination and similar divergence times to mammals, and teleost fish, which have much deeper divergence times. Rates of molecular evolution of recombination genes are highly correlated between vertebrate clades and significantly elevated compared to control panels, suggesting that they experience similar selective pressures. Avian recombination genes are significantly more likely to exhibit signatures of positive selection than other clades, unrestricted to later stages of the pathway. Signatures of positive selection in genes linked to recombination rate variation in mammalian populations and those with signatures of positive selection across the avian phylogeny are highly correlated. In contrast, teleost fish recombination genes have significantly less evidence of positive selection despite high intra-clade recombination rate variability. Gaining clade-specific understanding of patterns of variation in recombination genes can elucidate drivers of recombination rate and thus, factors influencing genetic diversity, selection efficacy, and species divergence.

Animals↗
PubMed2026

Genome evolution and long-term demographic history in true crocodiles.

Reference-quality genomes remain scarce for true crocodiles (Crocodylus), limiting comparative analyses of genome evolution and demographic history. Here, we generated and analyzed 2 long-read genomes, 1 for Crocodylus intermedius and 1 for C. niloticus, to investigate genome architecture, coalescent effective population size (Ne), and patterns of molecular evolution across crocodilians. Comparative analyses revealed broadly similar repeat landscapes in both species and extensive macro-synteny with Alligator sinensis, indicating strong structural conservation across crocodilian genomes. Using phased diploid assemblies and MSMC2, we reconstructed historical Ne trajectories and found marked differences between species. Crocodylus intermedius exhibited persistently low Ne throughout most of the late Quaternary, with a pronounced decline during the Late Pleistocene-early Holocene transition. In contrast, C. niloticus showed substantially larger Ne over comparable time intervals. Genome-wide codon-based analyses identified significant heterogeneity in dN/dS (&#x3c9;) among crocodilian lineages. Crocodylus niloticus showed the lowest genome-wide &#x3c9;, whereas elevated values in C. intermedius and other lineages were consistent with reduced long-term efficacy of purifying selection under smaller historical population sizes. Branch-site tests identified candidate genes under positive selection in both focal species, with functional categories related to ion transport, endocrine regulation, and cellular signaling. Together, these results provide genomic resources for Crocodylus and support an association between long-term demographic history and genome-wide patterns of molecular evolution across crocodilians.

Animals↗
PubMed2025

Unequally Abundant Chromosomes and Unusual Collections of Transferred Sequences Characterize Mitochondrial Genomes of Gastrodia (Orchidaceae), One of the Largest Mycoheterotrophic Plant Genera.

The mystery of genomic alternations in heterotrophic plants is among the most intriguing in evolutionary biology. Compared to plastid genomes (plastomes) with parallel size reduction and gene loss, mitochondrial genome (mitogenome) variation in heterotrophic plants remains underexplored in many aspects. To further unravel the evolutionary outcomes of heterotrophy, we present a comparative mitogenomic study with 13 de novo assemblies of Gastrodia (Orchidaceae), one of the largest fully mycoheterotrophic plant genera, and its relatives. Analyzed Gastrodia mitogenomes range from 0.56 to 2.1 Mb, each consisting of numerous, unequally abundant chromosomes or contigs. Size variation might have evolved through chromosome rearrangements followed by stochastic loss of "dispensable" chromosomes, with deletion-biased mutations. The discovery of a hyper-abundant (&#x223c;15 times intragenomic average) chromosome in two assemblies represents the hitherto most extreme copy number variation in any mitogenomes, with similar architectures discovered in two metazoan lineages. Transferred sequence contents highlight asymmetric evolutionary consequences of heterotrophy: despite drastically reduced intracellular plastome transfers convergent across heterotrophic plants, their rarity of horizontally acquired sequences sharply contrasts parasitic plants, where massive transfers from their hosts prevail. Rates of sequence evolution are markedly elevated but not explained by copy number variation, extending prior findings of accelerated molecular evolution from parasitic to heterotrophic plants. Putative evolutionary scenarios for these mitogenomic convergence and divergence fit well with the common (e.g. plastome contraction) and specific (e.g. host identity) aspects of the two heterotrophic types. These idiosyncratic mycoheterotrophs expand known architectural variability of plant mitogenomes and provide mechanistic insights into their content and size variation.

Genome, Mitochondrial↗
PubMed2026

The genetic code at the balance point of error and demand.

The origin and organizing principles of the genetic code remain central problems in molecular evolution. The low probability of the natural codon-to-amino acid mapping arising by chance has spurred the hypothesis that its structure is optimized for robustness to mutations and translational errors. For the construction of effective molecular machines, the repertoire of encoded amino acids must also be diverse enough in physicochemical features. Here, we examine whether the standard genetic code can be understood as a near-optimal solution balancing these two objectives: minimizing error load and aligning codon assignments with the naturally occurring amino acid composition. Using simulated annealing, we explore this trade-off across a broad range of parameters. We find that the standard genetic code resides near an optimum in the fitness landscape of possible genetic codes. The degeneracy of the code plays a dual role, minimizing mistranslation errors while matching codon multiplicity to amino acid usage frequencies. As a result, uniform codon usage alone is sufficient to recover the empirical amino acid composition, without any additional bias. It is a highly effective solution that balances fidelity against resource availability constraints. A comparative analysis of natural variants also reveals a functional decoupling: error robustness acts as a rigid global constraint determined by code topology, whereas compositional alignment serves as a more flexible variable that adapts to lineage-specific demands. These results support a multi-objective optimization framework in which the genetic code reflects a balance between translational fidelity and proteomic demand.

Genetic Code↗
PubMed2025

Compensatory Evolution Following Deleterious Episodes of GC-biased Gene Conversion in Rodents.

GC-biased gene conversion (gBGC) is a widespread evolutionary force associated with meiotic recombination that favors the accumulation of deleterious AT to GC substitutions in proteins, moving them away from their fitness optimum. In many mammals, recombination hotspots have a rapid turnover, leading to episodic gBGC, with the accumulation of deleterious mutations stopping when the recombination hotspot dies. Selection is therefore expected to act to repair the damage caused by gBGC episodes through compensatory evolution. However, this process has never been studied or quantified so far. Here, we analyzed the nucleotide substitution pattern in coding sequences of a highly diversified group of Murinae rodents. Using phylogenetic analyses of about 70,000 coding exons, we identified numerous exon-specific, lineage-specific gBGC episodes, characterized by a clustering of synonymous AT to GC substitutions and by an increasing rate of nonsynonymous AT to GC substitutions, many of which are potentially deleterious. Analyzing the molecular evolution of the affected exons in downstream lineages, we found evidence for pervasive compensatory evolution after deleterious gBGC episodes. Compensation appears to occur rapidly after the end of the episode and to be driven by the standing genetic variation rather than new mutations. Our results demonstrate the impact of gBGC on the evolution of amino-acid sequences and underline the key role of epistasis in protein adaptation. This study contributes to a growing body of literature emphasizing that adaptive mutations, which arise in response to environmental changes, are just 1 subset of beneficial mutations, alongside mutations resulting from oscillations around the fitness optimum.

Gene Conversion↗
PubMed2026

Comparative genomic and proteomic analysis reveals orthogroup structured evolution of tick protease inhibitors.

Protease inhibitors (PIs) play central roles in regulating endogenous proteolysis and host-parasite interactions in ticks. However, the evolutionary architecture underlying their diversification across tick lineages remains insufficiently resolved. Here, we performed a genome-wide comparative analysis of predicted proteomes from 14 tick species to systematically characterize PI repertoires. In total, 4931 putative PIs were identified and grouped into 20 families using the MEROPS classification system. Further, PI families such as Antistasin, WAP-type, and Pacifastin, which have not previously been systematically reported in tick genomes, were classified. Orthogroup inference demonstrated that PI expansion is structured at the level of evolutionary lineages rather than uniformly across families. By stratifying orthogroups according to duplication burden and taxonomic conservation, we identified a broadly conserved single-copy core under strong purifying selection. Motif level analysis of serpin reactive center loops further revealed conservation of inhibitory specificity within single copy orthogroups and diversification of key functional residues in duplication-associated lineages. Integration of secretion prediction and tissue-resolved proteomics from Hyalomma anatolicum and Rhipicephalus microplus demonstrated that evolutionary stratification is reflected at the protein level. Together, these findings provide an orthogroup-resolved evolutionary framework linking duplication dynamics, molecular evolution, and tissue-level protein deployment. This integrative approach offers a systematic basis for prioritizing conserved and diversified PI lineages for future functional and anti-tick intervention studies.

Animals↗