Search PubMed⌕ Search

SEARCH · Search PubMed

Search Search PubMed

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 91 records · Page 5Linked to original sources
PubMed2026

Complete genome sequence of the Anaplasma phagocytophilum clinical isolate NCH-1.

Anaplasma phagocytophilum is an obligate intracellular gram-negative bacterium and etiologic agent of human granulocytic anaplasmosis. A. phagocytophilum genomic sequencing has historically been performed via short-read platforms. Our optimized bacterial isolation protocol combined with Nanopore sequencing produced a single, closed 1,481,805 bp circular A. phagocytophilum strain NCH-1 chromosome.

Anaplasma phagocytophilum↗
PubMed2026

Draft genome sequence of Enterococcus casseliflavus strain MBBL_MP4 isolated from healthy bovine milk.

We report the draft genome sequence of Enterococcus casseliflavus MBBL_MP4, recovered from healthy bovine milk. The 3.45-Mbp genome assembly comprises 27 contigs and indicates low pathogenic potential, with no acquired antimicrobial resistance or known virulence genes. This genome provides a valuable resource for the genomic characterization of bovine-associated E. casseliflavus.

Enterococcus casseliflavus↗
PubMed2026

Identification of a novel plant polerovirus in the whitefly Aleuroclava gordoniae.

We report the genome sequence of Aleuroclava gordoniae-associated polerovirus (AgAP), identified from the whitefly Aleuroclava gordoniae. The 5,650-nt AgAP genome contains 6 open reading frames. Phylogenetic analysis places AgAP within the genus Polerovirus, which comprises plant-infecting viruses. This study provides a genomic resource for further investigation of virus-insect associations.

plant virus↗
PubMed2026

Genetic mutations driving ciprofloxacin resistance in laboratory-evolved Salmonella Typhimurium.

Ciprofloxacin resistance in Salmonella Typhimurium is a significant public health concern, and the mechanisms by which the resistance evolves are poorly defined. Here, by serial passaging under antibiotic selection, we isolated ciprofloxacin-resistant S. Typhimurium mutants and subjected them to whole-genome sequencing to reveal the major mutations associated with resistance. The Low CipR mutant acquired four chromosomal mutations in ramR, icdA, lipB, and gyrA, and the High CipR mutant gained additional mutations in gyrB, yaiC, and corA. Functional characterization determined that mutations in ramR resulted in efflux pump upregulation, while disruptions in the TCA cycle caused by mutations in icdA and lipB led to metabolic alterations. These changes indirectly enhanced resistance by increasing the expression of the global regulator MarA and reducing OmpF-dependent membrane permeability. Despite the observation of the G105A substitution in GyrA, enzymatic assays confirmed the failure to support resistance to ciprofloxacin, possibly because the structural alteration remained minimal. GyrB488-489dup was associated with maintained supercoiling under ciprofloxacin and enhanced fluoroquinolone resistance, suggesting a major role in resistance evolution. Other mutations in yaiC impaired biofilm and, in corA, intracellular accumulation of magnesium, possibly stabilizing the bacterial cell envelope under antibiotic pressure. The findings provide novel explanations for the multifaceted mechanisms leading to ciprofloxacin resistance in Salmonella and suggest targets to combat antimicrobial resistance.IMPORTANCEAntibiotic resistance in Salmonella Typhimurium is an increasing public health concern, yet the genetic changes that allow bacteria to become resistant are not fully understood. In this study, we evolved ciprofloxacin-resistant Salmonella in the laboratory and identified the mutations that arise during resistance development. We found that resistance does not result from a single change but from multiple adaptations affecting drug efflux, metabolism, and the antibiotic target. Some mutations increased the activity of pumps that remove antibiotics from the cell, while others altered bacterial metabolism and reduced membrane permeability, making it harder for the drug to enter. A duplication in the DNA gyrase subunit GyrB played a particularly important role in maintaining DNA function under antibiotic stress. Together, these results reveal how diverse genetic changes cooperate to generate ciprofloxacin resistance and provide insights that may help guide strategies to combat drug-resistant Salmonella infections.

DNA gyrase↗
PubMed2026

Genome-resolved analysis reveals disruption of gut microbial vitamin B and K2 biosynthesis during Toxoplasma gondii infection in mice.

UNLABELLED: Toxoplasma gondii infection remodels the gut microbiome, yet its impact on microbial vitamin biosynthetic potential and host redox metabolism remains unclear. Here, we integrated mouse gut metagenomes with publicly available metagenome-assembled genomes (MAGs) to construct a genome-resolved atlas of B-vitamin and vitamin K2 biosynthesis. From 45,697 MAGs, we curated 4,771 representative genomes, of which 2,682 met high-quality criteria (completeness &#x2265;90%, contamination <5%). Functional annotation identified 229,717 vitamin-related genes corresponding to 177 Kyoto Encyclopedia of Genes and Genomes (KEGG) orthologs across de novo pathways for eight B vitamins, thiamine (B1), riboflavin (B2), niacin (B3), pantothenate (B5), pyridoxine (B6), biotin (B7), folate (B9), cobalamin (B12), and vitamin K2. Among the high-quality genomes, 1,665 encoded complete de novo pathways for at least one vitamin, highlighting functional specialization and community-level complementarity. Transcripts per million-normalized metagenomic read counts revealed significant differences in KEGG ortholog abundances across six of the nine vitamin pathways. Reanalysis of metagenomic data from infected mice (acute, chronic, and control; n = 10 per group) revealed a stage-dependent reduction in &#x3b1;-diversity of vitamin biosynthesis pathways during acute infection, and a clear &#x3b2;-diversity separation from chronic and control groups. Core niacin biosynthesis genes (nadB, nadA, nadC) displayed phylum-specific redistribution, indicating selective remodeling of microbial NAD+ precursor production under infection-induced metabolic stress. These results suggest that T. gondii infection disrupts cooperative vitamin biosynthetic networks while specifically modulating niacin pathways linked to host NAD+ metabolism. IMPORTANCE: Gut microbes can synthesize essential vitamins, but how infection alters this function is poorly understood. By integrating mouse gut metagenomes with genome-resolved microbial data, we show that Toxoplasma gondii infection reshapes the vitamin biosynthetic potential of the gut microbiome in a stage-dependent manner. Acute infection reduces the diversity of vitamin biosynthesis pathways and shifts the taxonomic distribution of key niacin biosynthesis genes involved in microbial NAD+ precursor production. These findings identify vitamin metabolism, especially niacin-related pathways, as a sensitive functional axis of microbiome remodeling during infection. Our work links microbial taxonomic changes to functional metabolic consequences and suggests that microbiome-mediated regulation of NAD+-related metabolism may contribute to host redox adaptation during T. gondii infection.

B vitamins↗
PubMed2026

Targeting Both Oncogenic Signaling and Dependence Receptor Function is Required to Fully Suppress MET Exon 14 Skipping-Driven tumorigenesis.

Receptor tyrosine kinases (RTKs) classically function as oncogenic drivers that promote survival and proliferation upon ligand binding. A subset of RTKs can also function as dependence receptors, inducing apoptosis in the absence of their ligands. Genetic alterations that enhance RTK signaling are well characterized in cancer and can be targeted with kinase inhibitors, which show limited efficacy in some clinical settings. Elucidation of whether oncogenic mutations can promote tumorigenesis by directly abolishing the pro-apoptotic activity of dependence receptors could help improve strategies to target RTKs. Here, we identified MET exon 14 skipping (METex14Del) as a paradigmatic example of an oncogenic alteration that drives tumorigenesis through genetic inactivation of the dependence receptor function of an RTK. METex14Del removed both the caspase cleavage site and adjacent CBL-binding motif, preventing generation of the pro-apoptotic p40MET fragment while sustaining oncogenic MET signaling. Uncoupling regulatory functions of MET using genome editing showed that loss of apoptosis capacity is a critical determinant of METex14Del-driven tumorigenesis. Combined-but not individual-mutation of the caspase and CBL sites was sufficient to recapitulate resistance to apoptosis and tumor growth induced by METex14Del in HGF-humanized mouse models. Importantly, inducible re-expression of p40MET in METex14Del-expressing cells restored apoptotic sensitivity, decreased tumor formation in vivo, and resensitized tumors to capmatinib. Together, these findings redefine RTKs as receptors with dual oncogenic and tumor-suppressive functions and show that disruption of dependence receptor-mediated apoptosis is an oncogenic mechanism. These results provide a conceptual framework explaining why therapies targeting only RTK signaling may fail and support strategies restoring dependence receptor function to achieve durable tumor suppression.

Journal Article↗
PubMed2026

Trio-based whole-exome sequencing identifies convergent epithelial junction-related pathways in syndromic hidradenitis suppurativa.

INTRODUCTION: Hidradenitis suppurativa (HS)-related autoinflammatory syndromes, simply termed as syndromic HS (sHS), represent a group of rare immune-mediated inflammatory disorders in which HS coexists with systemic or cutaneous autoinflammatory features like PASH (pyoderma gangrenosum-PG-, acne and HS), PAPASH (PASH, pyogenic arthritis), PASS (PG, acne, HS, and ankylosing spondylitis), and SAPHO syndrome (synovitis, acne, pustulosis, hyperostosis, and osteitis). In recent years, genetic studies identified several novel pathogenic variants underlying sHS; however, most investigations rely exclusively on affected individuals sequencing and the absence of parental genomic information limits the possibility to determine inheritance patterns. METHODS: To address these gaps, we performed trio-based whole-exome sequencing (WES) on five individuals diagnosed with sHS and their unaffected parents. RESULTS: The pathway related to epidermal adhesion and desmosome organization was the most represented across our cohort, encompassing seven genes: DSC3, DSG1, FAT1, LAMA3, MICALL2, PLEC and TJP2. Integrin-extracellular matrix (ECM) adhesion signaling pathway, represented by ten genes (CSPG4, FERMT3, ITGA3, LAMA3, LAMA5, LIMS2, LTBP3, PLEC, TGM2, TNC) was also retrieved. Also, variants affecting innate immune pathways, including cytokine signalling and antigen presentation, have been observed. CONCLUSION: Our exploratory findings suggest that genetically heterogeneous variants in syndromic HS converge on biological processes involving epithelial junction organisation, extracellular matrix interactions and innate immune regulation. Although not establishing a unique pathogenic mechanism, these observations identify epithelial barrier biology as a candidate pathway warranting validation in larger cohorts and functional studies.

Journal Article↗
PubMed2026

Comparative Transcriptomics Reveals Shared Downstream Pathways in Craniofacial Pathology.

Treacher Collins syndrome and Nager syndrome are craniofacial developmental disorders caused by defects in ribosome biogenesis and RNA splicing, respectively, yet they exhibit overlapping abnormalities affecting neural crest cell-derived craniofacial structures. To investigate shared downstream pathogenic mechanisms, we performed a comparative transcriptomic analysis of zebrafish polr1c and sf3b4 mutant models from our previous studies. Comparative analysis identified 17 shared differentially expressed genes (DEGs) between polr1c and sf3b4 mutants, with the majority of shared genes dysregulated in the same direction, indicating a coordinated rather than random transcriptional response. Gene ontology analysis identified ATP-dependent protein folding chaperone activity as the only shared molecular function, driven in part by upregulation of hsp90aa1.2, indicating a common proteostasis response. Because chaperone activity is linked to extracellular matrix (ECM) protein processing, we cross-referenced DEGs from both mutants against the curated zebrafish matrisome. Three of the 17 shared DEGs (serpinh1b, il11a, and lepa) were matrisome-associated and upregulated in both mutants. Serpinh1b, a collagen-specific chaperone, was strongly expressed in craniofacial cartilage and mesenchymal populations during pharyngeal arch development and exhibited nearly identical fold changes in both mutants. Il11a is of particular interest because its receptor, IL11RA, is known to be associated with human craniosynostosis, suggesting potential relevance to craniofacial development. Together, it is possible to hypothesize that shared chaperone-associated transcriptional changes, together with altered ECM-related gene expression, may contribute to polr1c- and sf3b4-associated craniofacial disorders, warranting further functional validation.

Extracellular Matrix↗
PubMed2026

PaNDA: Efficient Optimization of Phylogenetic Diversity in Networks.

Phylogenetic diversity (PD) plays an important role in biodiversity, conservation, and evolutionary studies by measuring the diversity of a set of taxa based on their phylogenetic relationships. In phylogenetic trees, a subset of k taxa with maximum PD can be found by a simple and efficient greedy algorithm. However, this algorithmic tractability is lost when considering phylogenetic networks, which incorporate reticulate evolutionary events such as hybridization and horizontal gene transfer. To address this challenge, we introduce PaNDA (Phylogenetic Network Diversity Algorithms), the first software package and interactive graphical user-interface for exploring, visualizing, and maximizing diversity in phylogenetic networks. PaNDA includes a novel algorithm to find a subset of k taxa with maximum diversity, running in polynomial time for networks of bounded scanwidth, a measure of tree-likeness of a network that grows slower than the well-known level measure. This algorithm considers the variant of PD on networks in which the branch lengths of all paths from the root to the selected taxa contribute towards their diversity. We demonstrate the scalability of this algorithm on simulated networks, successfully analyzing level-15 networks with up to 200 taxa in seconds. We also provide a proof-of-concept analysis using a phylogenetic network on Xiphophorus species, illustrating how the tool can support diversity studies based on real genomic data. The software is easily installable and freely available at https://github.com/nholtgrefe/panda. Additionally, we extend the definition of PD to semi-directed phylogenetic networks, which are mixed graphs increasingly used in phylogenetic analysis to model uncertainty of the root location. We prove that finding a subset of k taxa with maximum diversity remains NP-hard on semi-directed networks, but do present a polynomial-time algorithm for networks with bounded level.

network↗
PubMed2026

Age- and sex-adjusted genomic differences between Korean and Beat AML cohorts.

Genomic profiling plays a central role in risk stratification and therapeutic decision-making in acute myeloid leukemia (AML), yet the clinical implications of population-specific genomic architectures remain incompletely defined. We conducted a prospective, multicenter study of 603 adults with newly diagnosed AML in Korea, integrating targeted sequencing of 83 recurrently mutated genes with comprehensive clinical annotation across treatment intensities, including allogeneic hematopoietic stem cell transplantation (allo-HSCT). For contextual comparison, genomic profiles were evaluated against the Beat AML cohort. The overall genomic landscape was broadly conserved, supporting shared core disease biology across populations. However, RUNX1::RUNX1T1, CEBPA, GATA2, KIT, and DDX41 mutations were more frequent in the Korean cohort, whereas FLT3 and NPM1 mutations were less common. These differences translated into a distinct distribution of European LeukemiaNet (ELN) 2022 risk categories, with implications for therapeutic stratification. Notably, most DDX41 alterations were germline (3.2%), highlighting the need for systematic germline evaluation with implications for genetic counseling and donor selection. Although unadjusted overall survival appeared longer in the Korean cohort, this difference was not significant after adjustment for key clinical variables. These findings indicate that population-specific genomic distributions reshape the clinical application of risk stratification and support population-aware precision medicine strategies in AML.

Journal Article↗
PubMed2026

Immune dysregulation in depression and psychosis: summary of current evidence and future perspectives.

Despite compelling epidemiological, genetic and cellular evidence linking immune dysregulation to depression and schizophrenia (and other psychotic disorders), causality remains contested and no immune biomarker has yet demonstrated robust clinical utility. Emerging methodological approaches - from target trial emulation on observational data to functional genomics - offer a potential path towards precision immunopsychiatry and stratified immunomodulatory treatment.

Neuroimmunology↗
PubMed2026

Context-dependent functional diversity of dorsomedial posterior parietal neurons revealed by single-unit fMRI mapping during naturalistic viewing.

The dorsomedial posterior parietal cortex (dmPPC) plays an important role in episodic processing by integrating sensory, cognitive, and motor information across distributed brain systems. However, how individual dmPPC neurons participate in large-scale functional organization during naturalistic experience remains poorly understood. To address this question, we combined single-unit electrophysiology and awake fMRI in five rhesus macaques of both sexes viewing identical naturalistic video stimuli. Using single-unit fMRI mapping, we generated whole-brain neuron-BOLD functional maps by correlating individual neuronal activity with voxel-wise fMRI signals across the brain. We found that neuron-BOLD functional maps exhibited strong context-dependent organization, with neurons recorded during the same video context showing substantially greater similarity than neurons recorded during different video conditions. Compared with neuronal spiking activity or critical fMRI frames alone, neuron-BOLD functional maps more robustly captured contextual structure. Despite this shared large-scale organization, a substantial subset of neighboring neurons recorded simultaneously from the same electrode displayed markedly distinct whole-brain association patterns, revealing substantial local functional heterogeneity within the dmPPC. This local heterogeneity was not readily explained by waveform-based putative cell class or by opposing neuronal firing dynamics. In addition, distributed cortical and medial temporal regions exhibited highly context-dependent neuron-BOLD association patterns during naturalistic viewing. Together, these findings demonstrate that dmPPC neurons participate in dynamic and heterogeneous large-scale functional organization during naturalistic episodic processing. More broadly, this study establishes single-unit fMRI mapping as a framework for linking single-neuron activity to distributed whole-brain dynamics across contextual conditions.Significance Statement Using single-unit fMRI mapping, this study examined how individual dorsomedial posterior parietal cortex (dmPPC) neurons relate to large-scale brain activity during naturalistic video viewing in macaque monkeys. We found that neuron-BOLD functional maps exhibit strong context-dependent organization and capture contextual structure more robustly than neuronal spiking activity or fMRI frames alone. Despite this shared organization, a substantial subset of neighboring dmPPC neurons displayed markedly distinct whole-brain association patterns, revealing local functional heterogeneity that was not readily explained by waveform-based putative cell class or opposing firing dynamics. These findings provide insight into how local neuronal populations participate in distributed brain-wide functional organization during naturalistic episodic processing.

Journal Article↗
PubMed2026

Longitudinal Repeated Protein Measurements in a Multiethnic Cohort Identify Novel Diabetes Biomarkers That Reveal Unique Disease Pathways.

There is up to a fourfold increase in diabetes biomarkers identified with longitudinal repeated versus single time point proteomic measurements. The increase in biomarkers identified with longitudinal repeated measurements is supported by a similar proportion being nominated as causal for type 2 diabetes with Mendelian randomization. Proteins unique to the longitudinal repeated analyses highlighted biological pathways (e.g., posttranslational protein modification and cellular structure and cycle regulation) that were distinct from pathways enriched among the shared proteins (e.g., small-molecule metabolic and catabolic processes). Longitudinal protein measurements identify additional novel disease biomarkers and disparate biological pathways compared with single measurement analyses.

Journal Article↗
PubMed2026

Azithromycin-resistant Salmonella enterica Typhi with AcrB R717L/Q mutations in the United States.

BACKGROUND AND OBJECTIVES: Azithromycin is a critical oral treatment for typhoid fever caused by Salmonella enterica serovar Typhi (Salmonella Typhi), since XDR has rendered other first-line treatment options ineffective. Azithromycin resistance conferred by amino acid changes in AcrB, an AcrAB-TolC efflux pump component, represents an emerging public health concern. Leveraging phenotypic and genotypic data from U.S. Salmonella Typhi surveillance systems, this study describes the prevalence, phenotype and genomic epidemiology of Salmonella Typhi with AcrB mutations in U.S. patients since the first detection in 2015. METHODS: AST and WGS data of >3000 Salmonella Typhi isolates were used to identify all cases with an AcrB mutation in the United States (2015-2025). We calculated annual prevalence and MIC ranges. Phylogenetic analysis was used to contextualize U.S. cases of Salmonella Typhi with an AcrB mutation within all globally reported cases. RESULTS: While the prevalence of AcrB mutations in the United States is low (1.5%), it has risen significantly in recent years, from 0.2% in 2016-2022 to 2.2% in 2023-2025. This increase is predominantly driven by clonal expansion of existing strains circulating in South Asia. AcrB mutations do not reliably confer resistance to azithromycin (MIC&#x200a;&#x2265;&#x200a;32&#x2005;mg/L), complicating clinical interpretation. CONCLUSIONS: The prevalence of AcrB mutations in Salmonella Typhi is increasing in the United States, and likely globally, given that U.S. data function as an informal proxy for regions without routine surveillance infrastructure. Clinical outcomes data are needed to inform Salmonella Typhi treatment guidelines and potentially amend clinical breakpoints for azithromycin.

Journal Article↗
PubMed2026

ER proteostasis failure in HYOU1 deficiency alters B cells, neutrophils, and interferon signalling.

Hypoxia upregulated 1 (HYOU1) is a stress-inducible ER chaperone. We investigated 2 unrelated patients carrying biallelic HYOU1 variants and presenting with primary immunodeficiency. Patient 1, homozygous for p.Pro444His, displayed failure to thrive, hypoglycemia, B cell lymphopenia, and neutropenia. Patient 2, compound heterozygous for p.Arg262Gln and p.Pro757_Glu758insAla, exhibited recurrent infections, enteropathy, and hypogammaglobulinemia. In Patient 1, while HYOU1 transcription was preserved, the protein was severely reduced. Tunicamycin treatment of dermal fibroblasts showed a blunted unfolded protein response and defective induction of ER stress-responsive genes. Immunophenotyping showed near-absence of circulating B cells, and single-cell RNA sequencing of bone marrow identified an arrest at the pro-B cell stage. Neutrophils displayed hypogranulation and dysregulated IFN- and apoptosis-associated transcriptional signatures, unresponsive to G-CSF. HYOU1 deficiency hence results in ER stress-induced proteostasis failure that simultaneously impairs adaptive immunity through B cell developmental arrest and innate immunity through neutrophil dysfunction and IFN pathway imbalance. This work expands the spectrum of HYOU1 deficiency and further identifies ER proteostasis as a central determinant of immune homeostasis.

Journal Article↗
PubMed2026

Stage-specific ROMO1 in rheumatoid arthritis: predictive immune insights into the MIF pathway and HLA-DR/IL2RA axis via integrated GWAS, transcriptomic, single-cell, and spatial profiling.

Emerging evidence links reactive oxygen species modulator 1 (ROMO1), a key mitochondrial ROS regulator, to rheumatoid arthritis (RA) pathogenesis. However, its exact mechanism remains elusive given the conflicting evidence about its specific function. We used a four-level integrative framework combining multi-omics data and literature&#x2011;supported mechanistic inference. At the genetic level, Mendelian randomization (MR) was performed to explore potential causal relationships between ROMO1, IL2RA, HLA-DR, MIF, and RA risk, followed by differential expression analysis and machine learning-based feature selection to identify key mROS genes. The temporal expression dynamics of ROMO1 were assessed in RA progression. At the cellular and tissue levels, we integrated single-cell RNA sequencing and spatial transcriptomics to map cell-type-specific expression and synovial localization of ROMO1-related immune cells and pathways. Finally, our multi-omics findings were contextualized with literature-supported mechanistic inference. (1) MR results were consistent with a potential protective effect of ROMO1 on RA (OR&#x2009;=&#x2009;0.52) and its potential regulation of risk factors IL2RA (OR&#x2009;=&#x2009;0.46) and HLA-DR (OR&#x2009;=&#x2009;0.40). Conversely, IL2RA (OR&#x2009;=&#x2009;1.42), HLA-DR (OR&#x2009;=&#x2009;1.88), and MIF (OR&#x2009;=&#x2009;1.17) were positively associated with RA risk. Additionally, ROMO1 was identified as a top candidate diagnostic predictor with stage-specific dynamics: downregulated in the early but upregulated in the late/remission stages. (2) Single-cell RNA sequencing showed ROMO1's cell-specific expression in CD14+&#x2009;HLA-DR+&#x2009;CD74+&#x2009;monocytes and CD4+&#x2009;IL2RA+&#x2009;T cells. Cell communication analysis further suggested that these cells may participate in MIF pathway regulation. Spatial transcriptomics subsequently identified that ROMO1-related cells localized to synovial pathological regions, with MIF pathway changes correlated with RA progression. (3) Finally, literature-supported mechanistic inference suggests that ROMO1 may modulate mROS levels to promote anti-inflammatory M2 macrophage polarization, which could theoretically contribute to reduced systemic inflammation and the alleviation of multi-organ decline in RA. This integrated multi-omics investigation, supported by literature-based mechanistic inference, suggests ROMO1 as a stage-dependent biomarker candidate and potential immune regulator in RA.

Humans↗
PubMed2026

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

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

Animals↗