SEARCH · Search PubMed

Results for “Evolutionary investigation”

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.

At least 19 recordsLinked to original sources

PubMed2026

Untangling the Arisaema enigma: Investigating the complex evolutionary history and species relationships in North American Arisaema.

PREMISE: The evolutionary history of morphologically variable plant groups is often obscured by cryptic diversity, morphological convergence, and limited genetic data. Arisaema, a diverse genus within Araceae, exemplifies these challenges. Although some taxonomic treatments recognize only two species of North American Arisaema (A. dracontium and A. triphyllum), other studies have identified morphologically distinct groups within both taxa. Here, we reconstructed evolutionary relationships in North American Arisaema, assessed genetic structure and admixture, and tested the monophyly of proposed species. METHODS: We used 2b-RAD sequencing to generate genome-wide SNP data for 146 samples from 31 populations across the eastern United States. Phylogenetic relationships were inferred using maximum-likelihood and Bayesian approaches. Population structure and admixture were assessed using the program structure and principal component analysis (PCA). RESULTS: Both the Arisaema triphyllum and A. dracontium complexes formed well-supported monophyletic groups. Within the A. dracontium complex, we recovered three monophyletic lineages: A. dracontium, A. calciphilum, and A. macrospathum. In the A. triphyllum complex, A. quinatum, A. stewardsonii, and A. allegheniense consistently formed distinct groups. Relationships between A. pusillum and A. acuminatum, and among A. triphyllum s.s., A. purpurascens, and A. striatum were less clearly resolved, likely due to recent or incomplete divergence, gene flow, or polyploidy. CONCLUSIONS: The results support the monophyly of multiple newly proposed taxa within North American Arisaema, but additional sampling across the species' ranges is needed to fully resolve species boundaries. Our study provides the first evolutionary framework for this group, providing a foundation for future ecological, taxonomic, and conservation research in the genus.

Araceae↗

PubMed2026

Comparative dynamics of Japanese encephalitis virus adaptation in porcine macrophages and insect cells.

BACKGROUND: Japanese encephalitis virus (JEV) is a zoonotic mosquito-borne Orthoflavivirus that circulates primarily in birds and pigs. Previous observations of vector-free transmission between pigs indicates the possibility of single-host cycling in swine. Therefore, the aim of this work was to investigate the evolutionary pressure of single host cycling using a relevant primary cell culture model. METHODS: To investigate whether such single-host cycles affect viral infectivity, fitness and genomic adaptations, two strains and a reverse genetic cDNA-derived clone of JEV were serially passaged 12 times in primary porcine monocyte-derived macrophages (MDMs), in Aedes albopictus-derived C6/36 cells, and alternately between both cell types. Next-generation sequencing analysis was used to identify selected single nucleotide variants (SNVs) and haplotypes. Phenotype-to-genotype connections were confirmed using reverse genetics. RESULTS: For all viruses, serial passaging in MDMs - but not in C6/36 cells - led to a rapid increase in relative infectivity toward MDMs, accompanied by reduced plaque sizes in porcine endothelial cells. In contrast to C6/36 cells, MDM imposed a strong selective pressure, rapidly favoring selection of many SNVs and viral haplotypes. In addition, we identified a dominant selection of mutants with glutamic acid to lysine substitutions at positions 49 or 138 in the E protein, which explained the small plaque phenotype and caused viral sensitivity to heparin-mediated inhibition of attachment, indicating enhanced virus binding to glycosaminoglycans (GAG). The E138K mutant also explained the increased relative infectivity for MDM. CONCLUSION: This work demonstrates a high evolutionary pressure on JEV in MDM causing rapid selections of minor haplotypes. Furthermore, the efficient selection of E49K and E138K SNV, which were responsible for the phenotype, are likely caused by a selective pressure for GAG binding, observed in vitro with other mammalian cells.

Animals↗

PubMed2025

Ancient Introgression Explains Mitochondrial Genome Capture and Mitonuclear Discordance Among South American Collared Tropidurus Lizards.

Mitonuclear discordance-evolutionary discrepancies between mitochondrial and nuclear DNA phylogenies-can arise from various factors, including introgression, incomplete lineage sorting, recent or ancient demographic fluctuations, sex-biased dispersal asymmetries, among others. Understanding this phenomenon is crucial for accurately reconstructing evolutionary histories, as failing to account for discordance can lead to misinterpretations of species boundaries, phylogenetic relationships, and historical biogeographic patterns. We investigate the evolutionary drivers of mitonuclear discordance in the Tropidurus spinulosus species group, which contains nine species of lizards inhabiting open tropical and subtropical environments in South America. Using a combination of population genetic and phylogenomic approaches applied to mitochondrial and nuclear data, we identified different instances of gene flow that occurred in ancestral lineages of extant species. Our results point to a complex evolutionary history marked by prolonged isolation between species, demographic fluctuations, and potential episodes of secondary contact with genetic admixture. These conditions likely facilitated mitochondrial genome capture while diluting signals of nuclear introgression. Furthermore, we found no strong evidence supporting incomplete lineage sorting or natural selection as primary drivers of the observed mitonuclear discordance. Therefore, the unveiled patterns are most consistent with neutral demographic processes, coupled with ancient mitochondrial introgression, as the main factors underlying the mismatch between nuclear and mitochondrial phylogenies in this system. Future research could further explore the role of other demographic processes, such as asymmetric sex-biased dispersal, in shaping these complex evolutionary patterns.

Animals↗

PubMed2025

Unraveling evolutionary pathways: allopolyploidization and introgression in polyploid Prunus (Rosaceae).

Allopolyploidization, resulting from hybridization and subsequent whole-genome duplication (WGD), is a fundamental mechanism driving evolutionary diversification across various lineages within the Tree of Life. The polyploid Prunus (Rosaceae), significant for its economic and agricultural value, provides an ideal model for investigating the evolutionary dynamics associated with allopolyploidy. In this study, we utilized deep genome skimming (DGS) data to demonstrate a comprehensive analytical framework for elucidating the underlying allopolyploidy that includes a newly adapted tool (DGS-Tree2GD) tailored explicitly for accurately detecting WGD events. Additionally, we introduced two methods to evaluate the contribution of incomplete lineage sorting (ILS) to lineage diversification. Phylogenomic discordance analyses revealed that allopolyploidization, rather than ILS, played a dominant role in the origin and dynamics of polyploid Prunus. Moreover, we inferred that the uplift of the Himalayas from the Middle to Late Miocene was a key driver in the rapid diversification of the Maddenia clade, an endemic group in East Asia. This geological event facilitated extensive hybridization and allopolyploidization, particularly the introgression between the Himalayas-Hengduan and Central-Eastern China clades. This case study demonstrates the robustness and efficacy of our analytical approach in precisely identifying WGD events and elucidating the evolutionary mechanisms underlying allopolyploidization in polyploid Prunus.

Polyploidy↗

PubMed2026

Phylogenomics and evolution of the Lauraceae based on targeted capture data.

The family Lauraceae, a hyper-diverse magnoliid family comprising approximately 63 genera and over 3,000 species, plays a key ecological role in tropical and subtropical forests. Yet deep relationships among its nine tribes remain unresolved, likely due to limited sampling and complex evolutionary processes such as incomplete lineage sorting (ILS) and gene flow. To address these challenges, we generated datasets of 255 single-copy nuclear genes and chloroplast genomes using a newly designed Lauraceae-specific probe set, achieving the most comprehensive genus-level sampling (84%) to date. Phylogenomic analyses reconstructed a robust nuclear tree, which resolved the Neocinnamomeae as sister to the Caryodaphnopsideae and revealed pronounced gene tree conflict and pervasive cytonuclear discordance. To investigate the evolutionary processes underlying these patterns, comprehensive analyses were conducted. The results indicate that conflicting nuclear gene trees reflect the combined effects of ILS, gene tree estimation error, and gene flow, with ILS dominating across the core Lauraceae, whereas cytonuclear discordance is primarily driven by extensive gene flow. Diversification analyses further indicate that episodes of rapid lineage accumulation coincide with major gene flow events, suggesting a potential role of gene flow in the diversification of Lauraceae. Overall, this study provides a robust nuclear phylogenomic framework for Lauraceae and demonstrates that gene flow had profound effects on its evolutionary history, shedding light on the contribution of gene flow to the diversification of hyper-diverse tropical plant lineages.

Cytonuclear discordance↗

PubMed2026

Discovery and evolution of endogenous retroviruses in the genome of crab-eating macaque (Macaca fascicularis).

Endogenous retroviruses (ERVs) are a dynamic and biologically significant component of vertebrate genomes, with integration events spanning deep evolutionary time. The crab-eating macaque (Macaca fascicularis) is an important non-human primate model for biomedical research because of its close phylogenetic relationship to humans and its conservation status as an endangered species. However, the ERV complement of its genome has not been systematically characterized. Using the current highest-quality chromosome-level genome assembly for this species, we performed a genome-wide, homology-based survey of relatively intact ERV proviruses in M. fascicularis. We identified 106 proviral loci distributed across all chromosomes. Phylogenetic reconstruction based on conserved reverse transcriptase domains classified these elements into β-, γ-, and unclassified lineages, with β- and γ-retroviral lineages predominating. LTR divergence-based dating indicated that these proviruses represent multiple waves of historical retroviral activity and span a broad range of integration ages. This curated dataset provides a high-confidence reference set for investigating the evolutionary history and genomic impact of preserved ERV proviruses in an endangered primate model; however, it does not include degraded ERV fragments or solo LTRs.

Animals↗

PubMed2026

Genome-wide characterization of the tomato PERK gene family and its expression profiling under abiotic stresses.

UNLABELLED: This study presents the first systematic genome-wide characterization of the proline-rich extensin-like receptor kinases (PERK) gene family in tomato (Solanum lycopersicum) and their transcriptional responses under abiotic stresses. Using the latest SL4.0/ITAG4.0 genome assembly, we identified six SlPERK genes, all harboring the conserved Ser/Thr protein kinase domain. Evolutionary and structural analyses revealed strong purifying selection (Ka/Ks < 1), distinct exon-intron organizations, and the presence of stress- and hormone-responsive cis-regulatory elements in their promoters. Furthermore, post-transcriptional regulation by 57 miRNAs and complex protein-protein interaction networks were predicted. To validate their stress-responsive roles, two tomato cultivars (GMOTL-1 and Roma) were subjected to cold, heat, and salinity treatments. Quantitative RT-PCR analysis revealed cultivar-specific expression dynamics: SlPERK4 exhibited strong transient induction under cold and heat stress, while SlPERK6 was highly responsive to salinity. Notably, the GMOTL-1 cultivar displayed significantly higher and broader stress-responsive expression profiles compared to Roma, indicating a potential role of these SlPERK genes in cultivar-specific stress tolerance. These findings provide a comprehensive genomic resource and establish a critical foundation for the functional validation and molecular breeding for stress-resilience tomato cultivars. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at https://doi.org/10.1007/s13205-026-05044-y.

Abiotic stress↗

PubMed2026

Divergent PXR function in seals: Endocrine adaptation or functional loss?

Seals accumulate xenobiotics through dietary biomagnification and exposure to polluted marine environments, with contaminants concentrating in their blubber. Biotransformation mitigates xenobiotic toxicity by converting lipophilic compounds into excretable hydrophilic metabolites, a process coordinated by nuclear receptors including the Pregnane X Receptor (PXR), whose plastic ligand-binding domain enables broad xenobiotic sensing. By examining PXR in pinnipeds, we investigated the evolutionary conservation and functional characterization of PXR using genomic sequence analysis, protein structural prediction, and transactivation assays, revealing broadly conserved structural features alongside species-specific functional divergence in receptor responsiveness to environmental stressors. Specifically, the obtained results highlight divergent gene and functional landscapes with ORF-disrupting mutations identified in Monachus monachus and Neomonachus schauinslandi that abolish receptor activation toward known PXR ligands. In contrast, Leptonychotes weddelli retained an intact PXR ORF but showed reduced receptor activity, revealing functional divergence in PXR among pinnipeds.

Biotransformation↗

PubMed2026

Lineage-specific adaptation and resistance in Candida albicans.

Candida albicans exhibits substantial phenotypic and ecological diversity; however, the exact relationship between its population structure, adaptation to specific niches, and antifungal resistance remains incompletely understood. To investigate these evolutionary dynamics, we analyzed the whole-genome sequences from 591 publicly available isolates, integrating nuclear and mitochondrial phylogenomics with ecological and resistance-associated genomic analyses. Phylogenomic analyses resolved 18 core nuclear clusters together with multiple admixed lineages. Strong cytonuclear concordance was noted in the majority of the central lineages, contrasting with a higher discordance among the admixed groups, consistent with recurrent genetic exchange. The analysis revealed that geographic origin explains a larger fraction of genetic variance than anatomical niche, supporting a predominantly generalist population structure. A notable exception was Cluster N16 (Candida africana), which presented a strict genital origin in our dataset (n = 34). Additionally, although the mitochondrial genome exhibits strong purifying selection, candidate residues under diversifying selection correlated with specific niches (e.g., bloodstream) have been identified. Analysis of five resistance-associated genes (ERG11, UPC2, FKS1, TAC1 and FUR1) revealed that resistance-associated variants were generally rare but exhibited distinct gene-specific patterns. In case of ERG11 and FUR1 they were concentrated in a specific clade (N11, N17, and their admixed Group A) and exhibit gene-dependent zygosity patterns. In summary, the evolution of C. albicans appears to be driven by a predominantly clonal model punctuated by episodic genetic exchange, where both ecological adaptation and antifungal resistance mutations exhibit genomic signatures marked by lineage specificity.

Antifungal resistance↗

PubMed2026

Global lessons from antibiotic resistance: Metformin-hydrolysing genes in transposable elements, a new threat for type II diabetic patients?

OBJECTIVES: To investigate the evolutionary origin, genomic mobility, and potential dissemination of metformin-hydrolysing genes (mfmAB), and to assess whether environmental selection by metformin pollution may drive the emergence of transferable pharmaceutical-degrading traits analogous to antibiotic resistance. METHODS: Large-scale comparative genomics was performed using publicly available bacterial genomes carrying mfmAB homologs. Phylogenomic reconstruction, average nucleotide identity analysis, genomic context comparison, plasmid characterization, and insertion sequence mapping were used to infer evolutionary history and identify mechanisms of horizontal gene transfer. RESULTS: mfmAB homologs were identified in twelve Aminobacter and three Pseudomonas genomes within a conserved ∼8.2 kb gene cluster. Phylogenomic analyses showed that metformin-degrading capacity emerged independently in multiple Aminobacter lineages across distinct continents, consistent with convergent evolution under anthropogenic selective pressure. Genomic comparisons indicated a chromosomal origin of mfmAB, followed by mobilization onto conjugative plasmids through IS1182-mediated transposition. In Pseudomonas, additional IS3/IS6-mediated transposition events integrated mfmAB into diverse plasmid backbones, frequently within composite transposons also encoding guanylurea and biguanide degradation pathways (guuH, bguH). These findings reveal a dynamic modular assembly of metabolic functions facilitating adaptation to pharmaceutical pollutants. CONCLUSIONS: Metformin pollution appears to promote the emergence and mobilization of pharmaceutical-degrading genes through mechanisms paralleling antibiotic resistance evolution. Although no clinical impact has yet been demonstrated, the potential spread of such genes into human-associated microbiomes and their possible co-selection with antibiotic resistance determinants represent an emerging One Health concern. Environmental surveillance of pharmaceutical-degrading genes is warranted to anticipate future threats to drug efficacy.

Convergent evolution↗

PubMed2025

Genome-scale evolution and phylodynamics of swine influenza A viruses in China: a genomic epidemiology study.

BACKGROUND: Pigs are recognised as crucial intermediate hosts for the emergence of influenza viruses of pandemic potential. As the largest pork-producing nation, China hosts a complex ecosystem of swine influenza viruses (SIVs). We aimed to investigate the evolutionary processes, spatiotemporal dynamics, and biological characteristics of SIVs in China. METHODS: From Jan 15, 2016, to Dec 22, 2020, we collected nasal swabs from pigs at eight abattoirs and 16 swine farms in the Guangdong, Henan, and Shandong provinces of China, as part of SIV surveillance. SIVs were detected with RT-PCR. Positive samples underwent viral isolation and genome sequencing. We analysed evolution and spatiotemporal dynamics using the whole genomes of isolated SIVs, as well as genome sequences of SIV isolates from human infections worldwide retrieved from the Global Initiative on Sharing All Influenza Data and GenBank Flu databases up to April 28, 2024. Viral sequences without a sample collection area or date were excluded from the analysis. Viral receptor-binding properties and in-vitro replication of strains isolated in this study were evaluated with a solid-phase binding assay and various cell lines, including Madin-Darby canine kidney cells, porcine alveolar macrophages, primary porcine trachea epithelial cells, human bronchial epithelioid, and human lung adenocarcinoma epithelial (A549) cells. Viral replication and transmission studies were conducted in 33 guinea pigs and 13 pigs. Additionally, we collected serum samples from pig farm workers and members of the general public recruited by the Third Affiliated Hospital of Sun Yat-sen University between Feb 28 and May 11, 2023, to detect specific antibodies against Eurasian avian-like A(H1) and human-like A(H3N2) SIVs using the haemagglutination inhibition assay. FINDINGS: 23 (1·3%) of 1818 nasal swabs collected in abattoirs had SIVs; 22 (0·9%) of 2375 swabs from swine farms had SIVs. Further viral isolation yielded 39 strains of SIV. We identified 534 A(H1N1), 69 A(H1N2), and 92 A(H3N2) SIVs, representing 20 genotypes within the Eurasian avian-like lineage, 14 within the classical swine A(H1) lineage, and 16 within the human-like A(H3N2) lineage. The introduction of the A(H1N1)pdm/09 virus significantly influenced the internal gene pool of SIVs, enhancing genotypic diversity in China. Notably, the Eurasian avian-like A(H1), classical swine A(H1), and human-like A(H3N2) lineages showed human-mediated spread over long distances between provinces, with the Eurasian avian-like A(H1) lineage showing the most prevalent spread pathways. Eurasian avian-like A(H1) SIVs showed a preference for binding to sialic acid α-2,6 glycan receptors, predominantly found in humans, resulting in an increased production of progeny viruses in human airway epithelial cells, as well as effective transmission and infectivity among guinea pigs and pigs. Among 54 eligible serum samples collected from pig farm workers (24 from slaughterhouses and 30 from swine farms), 23 (43%) were seropositive for Eurasian avian-like A(H1) SIVs and 46 (85%) for human-like A(H3N2) SIVs. Among 100 eligible samples from members of the general public, 14 (14%) were seropositive for Eurasian avian-like A(H1) SIVs and 85 (85%) for human-like A(H3N2) SIVs. INTERPRETATION: This study elucidates the evolutionary processes and spatiotemporal patterns of SIVs, highlighting potential risks to public health. These findings are crucial for informing public health interventions that aim to prevent future SIV epidemics in China and other countries worldwide. FUNDING: Scientific Innovation Strategy-Construction of High-Level Academy of Agriculture Science-Distinguished Scholar (R2020PY-JC001).

Animals↗

PubMed2026

Spatiotemporal patterns of Rift Valley fever virus in Africa: a retrospective genomic epidemiology and phylodynamic modelling study.

BACKGROUND: Rift Valley fever virus (RVFV) is a mosquito-borne zoonotic pathogen causing outbreaks in humans and ruminants across Africa and the Arabian Peninsula. Originally restricted to the Great Rift Valley, RVFV has expanded geographically, prompting its classification by WHO as a pathogen of pandemic potential. We investigated the evolutionary and spatial dynamics of RVFV across Africa. METHODS: We used genomic data generated at the International Livestock Research Institute Nairobi genomic laboratory (BioProject PRJNA1106221) and combined with publicly available datasets retrieved from the National Center for Biotechnology (NCBI) GenBank nucleotide database. In retrieving RVFV genome sequences from the NCBI GenBank, we applied the search terms "Rift Valley fever virus segment L AND 6404[SLEN]", "Rift Valley fever virus segment M AND 3885[SLEN]", and "Rift Valley fever virus segment S AND 1520:1690[SLEN]" for L (Large), M (Medium), and S (Small) segments, respectively. For sequences without additional spatiotemporal information, we searched PubMed to extract the associated sequence metadata. We performed molecular clock analysis, phylogenetic inference, phylodynamic modelling (continuous phylogeographic reconstruction), and landscape phylogeography on the three RVFV genome segments (L, M, and S). We aimed to assess evolutionary rates, dispersal patterns, and environmental drivers. Focus was placed on lineage C, the most widely distributed variant. FINDINGS: The global dataset used in this study consisted of large (n=236), medium (n=237), and small (n=247), which were further filtered to exclude potential reassortants and vaccine strains. Genome sequences retrieved from NCBI GenBank database comprised large (n=180), medium (n=184), and small (n=202). The genome sequences from retrospective human and livestock isolates comprised large (n=56), medium (n=53), and small (n=45) collected in Burundi (2018), Kenya (2007, 2018, 2019, 2021, and 2022), and Rwanda (2018 and 2022). Our dataset revealed that RVFV exhibited low overall genetic diversity. Lineage C, however, showed evidence of active evolution, with substitution rates ranging from 3·58 × 10-4 to 9·76 × 10-4 substitutions per site per year. This lineage probably originated in Zimbabwe in the mid-1970s and has since expanded across eastern and southern Africa. Phylogeographic reconstructions revealed rapid spread, with diffusion coefficients exceeding 50 000 km2 per year. INTERPRETATION: Lineage C appears capable of establishing endemic transmission in new regions, with ongoing diversification observed during interepidemic periods. These observations reinforce the value of continuous genomic surveillance, particularly during cryptic transmission phases when adaptive mutations might emerge. Although further evidence is needed, observed trends in climate variability and land-use change point to the potential benefit of targeted surveillance in settings that could be at increased risk, including urban centres and wetlands. FUNDING: This work was supported by the German Federal Ministry for Economic Cooperation and Development, the Rockefeller Foundation, and the Africa Centres for Disease Control and Prevention.

Rift Valley fever virus↗

PubMed2025

A telomere-to-telomere reference genome assembly of the red silk cotton tree (Bombax ceiba).

Bombax ceiba, an important ornamental tree and potential fiber resource in the textile industry, is widely distributed in tropical and subtropical regions. In this study, we assembled a nearly gap-free telomere-to-telomere (T2T) genome of B. ceiba using Illumina, PacBio High-fidelity (HiFi), ONT ultra-long, and Hi-C sequencing technologies. The genome spanned approximately 807.89 Mb, with a scaffold N50 of 16.58 Mb, and 754.68 Mb (93.41%) of genomic sequences were anchored onto 48 pseudo-chromosomes. Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis revealed a completeness of 99.40%, identifying 1,378 single-copy and 213 duplicated genes out of 1,614. The genome contained 67.72% (547.11 Mb) repeat regions, with 39,708 predicted protein-coding genes. Collectively, our study provides valuable genomic data for investigating the evolutionary history of the Malvaceae family.

Genome, Plant↗

PubMed2026

Molting in Pancrustacea Is Characterized by Both Deeply Conserved and Recently Evolved Gene Modules.

Arthropods such as insects and crustaceans, which together form the monophyletic group Pancrustacea, possess a rigid chitinous exoskeleton that must be periodically shed through molting to allow growth and morphological change. Although molting is a deeply conserved developmental process across Arthropoda, our understanding of its molecular mechanisms is still largely derived from insect model species. Lineage-specific innovations and losses of molting-related genes raise fundamental questions about the extent of its conservation outside noninsect arthropods. Here, we investigate the evolutionary conservation of molting gene expression across five representative pancrustacean species using publicly available transcriptomic datasets. Changes in gene expression during molting are characterized by both deeply conserved and lineage-specific gene modules. Temporal gene expression analyses reveal that these lineage-specific signatures are not uniformly distributed across the molting process: the middle transitional phase is more lineage-specific, thereby exhibiting an inverse hourglass pattern. This is likely due to life-history-specific processes, development of the cuticle, and specialized structures of the exoskeleton. Overall, this study provides evidence for both the evolutionary conservation and divergence of this key postembryonic developmental process and highlights the modular architecture of the molting program.

Animals↗

PubMed2026

A UG5 reverse transcriptase-nitrilase antiviral module confers phage immunity in the plant symbiont Sinorhizobium meliloti.

Bacteriophages exert strong selective pressure on soil- and rhizosphere-associated bacteria, including plant-associated symbionts. Reverse transcriptase-associated defense systems of the UG family are widespread across bacterial lineages, yet their ecological roles remain largely undefined. Within this family, UG5 systems are distinguished by reverse transcriptases fused to or associated with a nitrilase domain. Here, we combine phylogenetic, metagenomic, and functional analyses to investigate the evolutionary context and antiviral activity of UG5-associated systems. Phylogenetic analysis of 728 nitrilase domains places UG5-associated nitrilases within a well-supported UG-related radiation encompassing the UG1, UG5, and UG6 families, with UG1 nested within a broader UG5 lineage. Metagenomic analysis further revealed UG5-associated reverse transcriptases in soil- and rhizosphere-derived metagenomes. Based on this observation, we characterized a UG5-large reverse transcriptase (RT)-associated system, here designated DRT11, encoded on the pSymA megaplasmid of Sinorhizobium meliloti RMO17, a nitrogen-fixing symbiont of Medicago sativa. Despite lacking the transmembrane protein typical of canonical UG5-large architectures, DRT11 confers protection against naturally occurring M. sativa rhizosphere phages with podovirus-like morphology. Phage infection assays reveal protection at low multiplicities of infection, consistent with an abortive-infection-like mechanism. Moreover, mutational analyses demonstrate that antiviral activity requires only the RT and its fused C-terminal nitrilase domain, establishing DRT11 as a minimal UG5-associated antiviral system.IMPORTANCEIn this study, we report the functional characterization of a UG5-large reverse transcriptase-associated defense system (DRT11) encoded on the pSymA megaplasmid of the nitrogen-fixing plant symbiont Sinorhizobium meliloti. Using a combination of phylogenetic, metagenomic, genomic, and experimental approaches, we demonstrate that DRT11 functions as a bona fide antiviral defense module, providing protection against naturally occurring rhizosphere phages through a minimal reverse transcriptase-nitrilase architecture. This work establishes direct functional evidence for antiviral activity within the UG5 family and clarifies the evolutionary placement of UG5-associated systems within the broader UG radiation.

Phylogeny↗

PubMed2026

Tandem duplication-driven expansion and UV-B stress adaptation of the LHC gene family in Artemisia annua L.

BACKGROUND: Artemisia annua L., is the primary natural source of the antimalarial drug artemisinin. In nature, fluctuating light is a major environmental stress that affects plant growth and artemisinin biosynthesis. Although the light-harvesting chlorophyll a/b-binding (LHC) superfamily plays a key role in mediating plant responses to fluctuating light, systematic research of this gene family in A. annua has not yet been conducted, limiting our understanding of light adaptation in this medicinally important species. RESULTS: This study investigated the evolutionary dynamics and functional adaptation of the light-harvesting chlorophyll a/b-binding (LHC) superfamily in A. annua, with a focus on the early light‑induced protein (ELIP) subfamily. Comparative genomics of 24 plant species showed that the LHC superfamily recently expanded in the examined Asteraceae lineages through duplication events. In A. annua, 229 LHC genes identified from four haplotype genomes comprised 205 allelic and 24 haplotype-specific loci, with the ELIP subfamily expanding significantly via tandem duplication. Notably, compared to non-Asteraceae plants, ELIPs exhibited a uniform single-exon architecture, indicating it is a genomic feature unique to Asteraceae plants. Population genomics of 41 individuals showed dynamic copy number variations ranging from 1 to 4 copies per locus. Interestingly, a structurally disrupted ELIP allele remained transcriptionally active and produced long aberrant transcripts, showing that this subfamily is still actively evolving. Under UV-B stress, AaELIP loci showed synchronized induction trend but differed in expression levels, suggesting a division into major and auxiliary roles within the expanded tandem cluster. Overall, while the response of ELIPs to light stress is evolutionarily conserved, this dramatic expansion and structural streamlining of AaELIPs may represent a key evolutionary adaptation that enhances the plant's ability to cope with intense light and radiation stress. CONCLUSIONS: Collectively, this study demonstrates a significant expansion of the LHC superfamily in A. annua, especially within the ELIP subfamily, as well as its robust response to UV-B treatment, underscoring the essential role of ELIPs in mediating light stress responses. These findings provide a valuable foundation for future research to uncover the molecular mechanisms underlying A. annua's adaptation to complex light environments.

Artemisia annua↗

PubMed2024

Comparative genomic analysis and functional investigations for MCs catabolism mechanisms and evolutionary dynamics of MCs-degrading bacteria in ecology.

Microcystins (MCs) significantly threaten the ecosystem and public health. Biodegradation has emerged as a promising technology for removing MCs. Many MCs-degrading bacteria have been identified, including an indigenous bacterium Sphingopyxis sp. YF1 that could degrade MC-LR and Adda completely. Herein, we gained insight into the MCs biodegradation mechanisms and evolutionary dynamics of MCs-degrading bacteria, and revealed the toxic risks of the MCs degradation products. The biochemical characteristics and genetic repertoires of strain YF1 were explored. A comparative genomic analysis was performed on strain YF1 and six other MCs-degrading bacteria to investigate their functions. The degradation products were investigated, and the toxicity of the intermediates was analyzed through rigorous theoretical calculation. Strain YF1 might be a novel species that exhibited versatile substrate utilization capabilities. Many common genes and metabolic pathways were identified, shedding light on shared functions and catabolism in the MCs-degrading bacteria. The crucial genes involved in MCs catabolism mechanisms, including mlr and paa gene clusters, were identified successfully. These functional genes might experience horizontal gene transfer events, suggesting the evolutionary dynamics of these MCs-degrading bacteria in ecology. Moreover, the degradation products for MCs and Adda were summarized, and we found most of the intermediates exhibited lower toxicity to different organisms than the parent compound. These findings systematically revealed the MCs catabolism mechanisms and evolutionary dynamics of MCs-degrading bacteria. Consequently, this research contributed to the advancement of green biodegradation technology in aquatic ecology, which might protect human health from MCs.

Humans↗

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