Manipulating the EphB4-ephrinB2 axis to reduce metastasis in HNSCC.

The EphB4-ephrinB2 signaling axis has been heavily implicated in metastasis across numerous cancer types. Our emerging understanding of the dichotomous roles that EphB4 and ephrinB2 play in head and neck squamous cell carcinoma (HNSCC) poses a significant challenge to rational drug design. We find that EphB4 knockdown in cancer cells enhances metastasis in preclinical HNSCC models by augmenting immunosuppressive cells like T regulatory cells (Tregs) within the tumor microenvironment. EphB4 inhibition in cancer cells also amplifies their ability to metastasize through increased expression of genes associated with epithelial mesenchymal transition and hallmark pathways of metastasis. In contrast, vascular ephrinB2 knockout coupled with radiation therapy (RT) enhances anti-tumor immunity, reduces Treg accumulation into the tumor, and decreases metastasis. Notably, targeting the EphB4-ephrinB2 signaling axis with the engineered EphB4 ligands EFNB2-Fc-His and Fc-TNYL-RAW-GS reduces local tumor growth and distant metastasis in a preclinical model of HNSCC. Our data suggest that targeted inhibition of vascular ephrinB2 while avoiding inhibition of EphB4 in cancer cells could be a promising strategy to mitigate HNSCC metastasis.

Suppression of the host antiviral response by non-infectious varicella zoster virus extracellular vesicles.

Varicella zoster virus (VZV) reactivates from ganglionic sensory neurons to produce herpes zoster (shingles) in a unilateral dermatomal distribution, typically in the thoracic region. Reactivation not only heightens the risk of stroke and other neurological complications but also increases susceptibility to co-infections with various viral and bacterial pathogens at sites distant from the original infection. The mechanism by which VZV results in complications remote from the initial foci remains unclear. Small extracellular vesicles (sEVs) are membranous signaling structures that can deliver proteins and nucleic acids to modify the function of distal cells and tissues during normal physiological conditions. Although viruses have been documented to exploit the sEV machinery to propagate infection, the role of non-infectious sEVs released from VZV-infected neurons in viral spread and disease has not been studied. Using multi-omic approaches, we characterized the content of sEVs released from VZV-infected human sensory neurons (VZV sEVs). One viral protein was detected (immediate-early 62), as well as numerous immunosuppressive and vascular disease-associated host proteins and miRNAs that were absent in sEVs from uninfected neurons. Notably, VZV sEVs are non-infectious yet transcriptionally altered primary human cells, suppressing the antiviral type 1 interferon response and promoting neuroinvasion of a secondary pathogen in vivo. These results challenge our understanding of VZV infection, proposing that the virus may contribute to distant pathologies through non-infectious sEVs beyond the primary infection site. Furthermore, this study provides a previously undescribed immune-evasion mechanism induced by VZV that highlights the significance of non-infectious sEVs in early VZV pathogenesis. IMPORTANCE: Varicella zoster virus (VZV) is a ubiquitous human virus that predominantly spreads by direct cell-cell contact and requires efficient and immediate host immune evasion strategies to spread. The mechanisms of immune evasion prior to virion entry have not been fully elucidated and represent a critical gap in our complete understanding of VZV pathogenesis. This study describes a previously unreported antiviral evasion strategy employed by VZV through the exploitation of the infected host cell's small extracellular vesicle (sEV) machinery. These findings suggest that non-infectious VZV sEVs could travel throughout the body, affecting cells remote from the site of infection and challenging the current understanding of VZV clinical disease, which has focused on local effects and direct infection. The significance of these sEVs in early VZV pathogenesis highlights the importance of further investigating their role in viral spread and secondary disease development to reduce systemic complications following VZV infections.

Diverse Gene Regulatory Mechanisms Alter Rattlesnake Venom Gene Expression at Fine Evolutionary Scales.

Understanding and predicting the relationships between genotype and phenotype is often challenging, largely due to the complex nature of eukaryotic gene regulation. A step towards this goal is to map how phenotypic diversity evolves through genomic changes that modify gene regulatory interactions. Using the Prairie Rattlesnake (Crotalus viridis) and related species, we integrate mRNA-seq, proteomic, ATAC-seq and whole-genome resequencing data to understand how specific evolutionary modifications to gene regulatory network components produce differences in venom gene expression. Through comparisons within and between species, we find a remarkably high degree of gene expression and regulatory network variation across even a shallow level of evolutionary divergence. We use these data to test hypotheses about the roles of specific trans-factors and cis-regulatory elements, how these roles may vary across venom genes and gene families, and how variation in regulatory systems drive diversity in venom phenotypes. Our results illustrate that differences in chromatin and genotype at regulatory elements play major roles in modulating expression. However, we also find that enhancer deletions, differences in transcription factor expression, and variation in activity of the insulator protein CTCF also likely impact venom phenotypes. Our findings provide insight into the diversity and gene-specificity of gene regulatory features and highlight the value of comparative studies to link gene regulatory network variation to phenotypic variation.

Assessing Target Specificity of the Small Molecule Inhibitor MARIMASTAT to Snake Venom Toxins: A Novel Application of Thermal Proteome Profiling.

New treatments that circumvent the pitfalls of traditional antivenom therapies are critical to address the problem of snakebite globally. Numerous snake venom toxin inhibitors have shown promising cross-species neutralization of medically significant venom toxins in vivo and in vitro. The development of high-throughput approaches for the screening of such inhibitors could accelerate their identification, testing, and implementation and thus holds exciting potential for improving the treatments and outcomes of snakebite envenomation worldwide. Energetics-based proteomic approaches, including thermal proteome profiling and proteome integral solubility alteration (PISA) assays, represent "deep proteomics" methods for high throughput, proteome-wide identification of drug targets and ligands. In the following study, we apply thermal proteome profiling and PISA methods to characterize the interactions between venom toxin proteoforms in Crotalus atrox (Western Diamondback Rattlesnake) and the snake venom metalloprotease (SVMP) inhibitor marimastat. We investigate its venom proteome-wide effects and characterize its interactions with specific SVMP proteoforms, as well as its potential targeting of non-SVMP venom toxin families. We also compare the performance of PISA thermal window and soluble supernatant with insoluble precipitate using two inhibitor concentrations, providing the first demonstration of the utility of a sensitive high-throughput PISA-based approach to assess the direct targets of small molecule inhibitors for snake venom.

Selected humanization of yeast U1 snRNP leads to global suppression of pre-mRNA splicing and mitochondrial dysfunction in the budding yeast.

The recognition of the 5' splice site (5' ss) is one of the earliest steps of pre-mRNA splicing. To better understand, the mechanism and regulation of 5' ss recognition, we selectively humanized components of the yeast U1 (yU1) snRNP to reveal the function of these components in 5' ss recognition and splicing. We targeted U1C and Luc7, two proteins that interact with and stabilize the yU1 snRNA and the 5' ss RNA duplex. We replaced the zinc-finger (ZnF) domain of yeast U1C (yU1C) with its human counterpart, which resulted in a cold-sensitive growth phenotype and moderate splicing defects. We next added an auxin-inducible degron to yeast Luc7 (yLuc7) protein (to mimic the lack of Luc7Ls in human U1 snRNP). We found that Luc7-depleted yU1 snRNP resulted in the concomitant loss of Prp40 and Snu71 (two other essential yU1 snRNP proteins), and further biochemical analyses suggest a model of how these three proteins interact with each other in the U1 snRNP. The loss of these proteins resulted in a significant growth retardation accompanied by a global suppression of pre-mRNA splicing. The splicing suppression led to mitochondrial dysfunction as revealed by a release of Fe2+ into the growth medium and an induction of mitochondrial reactive oxygen species. Together, these observations indicate that the human U1C ZnF can substitute that of yeast, Luc7 is essential for the incorporation of the Luc7-Prp40-Snu71 trimer into yU1 snRNP, and splicing plays a major role in the regulation of mitochondrial function in yeast.

Serine protease inhibitor, SerpinA3n, regulates cardiac remodelling after myocardial infarction.

AIMS: Following myocardial infarction (MI), the heart repairs itself via a fibrotic repair response. The degree of fibrosis is determined by the balance between deposition of extracellular matrix (ECM) by activated fibroblasts and breakdown of nascent scar tissue by proteases that are secreted predominantly by inflammatory cells. Excessive proteolytic activity and matrix turnover has been observed in human heart failure, and protease inhibitors in the injured heart regulate matrix breakdown. Serine protease inhibitors (Serpins) represent the largest and the most functionally diverse family of evolutionary conserved protease inhibitors, and levels of the specific Serpin, SerpinA3, have been strongly associated with clinical outcomes in human MI as well as non-ischaemic cardiomyopathies. Yet, the role of Serpins in regulating cardiac remodelling is poorly understood. The aim of this study was to understand the role of Serpins in regulating scar formation after MI. METHODS AND RESULTS: Using a SerpinA3n conditional knockout mice model, we observed the robust expression of Serpins in the infarcted murine heart and demonstrate that genetic deletion of SerpinA3n (mouse homologue of SerpinA3) leads to increased activity of substrate proteases, poorly compacted matrix, and significantly worse post-infarct cardiac function. Single-cell transcriptomics complemented with histology in SerpinA3n-deficient animals demonstrated increased inflammation, adverse myocyte hypertrophy, and expression of pro-hypertrophic genes. Proteomic analysis of scar tissue demonstrated decreased cross-linking of ECM peptides consistent with increased proteolysis in SerpinA3n-deficient animals. CONCLUSION: Our study demonstrates a hitherto unappreciated causal role of Serpins in regulating matrix function and post-infarct cardiac remodelling.

Dermonecrosis caused by a spitting cobra snakebite results from toxin potentiation and is prevented by the repurposed drug varespladib.

Snakebite envenoming is a neglected tropical disease that causes substantial mortality and morbidity globally. The venom of African spitting cobras often causes permanent injury via tissue-destructive dermonecrosis at the bite site, which is ineffectively treated by current antivenoms. To address this therapeutic gap, we identified the etiological venom toxins in Naja nigricollis venom responsible for causing local dermonecrosis. While cytotoxic three-finger toxins were primarily responsible for causing spitting cobra cytotoxicity in cultured keratinocytes, their potentiation by phospholipases A2 toxins was essential to cause dermonecrosis in vivo. This evidence of probable toxin synergism suggests that a single toxin-family inhibiting drug could prevent local envenoming. We show that local injection with the repurposed phospholipase A2-inhibiting drug varespladib significantly prevents local tissue damage caused by several spitting cobra venoms in murine models of envenoming. Our findings therefore provide a therapeutic strategy that may effectively prevent life-changing morbidity caused by snakebite in rural Africa.

Assessing Heterogeneity in the N-Telopeptides of Type I Collagen by Mass Spectrometry.

Collagen cross-links created by the lysyl oxidase and lysyl hydroxylase families of enzymes are a significant contributing factor to the biomechanical strength and rigidity of tissues, which in turn influence cell signaling and ultimately cell phenotype. In the clinic, the proteolytically liberated N-terminal cross-linked peptide of collagen I (NTX) is used as a biomarker of bone and connective tissue turnover, which is altered in several disease processes. Despite the clinical utility of these collagen breakdown products, the majority of the cross-linked peptide species have not been identified in proteomic datasets. Here we evaluate several parameters for the preparation and identification of these peptides from the collagen I-rich Achilles tendon. Our refined approach involving chemical digestion for protein solubilization coupled with mass spectrometry allows for the identification of the NTX cross-links in a range of modification states. Based on the specificity of the enzymatic cross-linking reaction we utilized follow-up variable modification searches to facilitate identification with a wider range of analytical workflows. We then applied a spectral library approach to identify differences in collagen cross-links in bovine pulmonary hypertension. The presented method offers unique opportunities to understand extracellular matrix remodeling events in development, aging, wound healing, and fibrotic disease that modulate collagen architecture through lysyl-hydroxylase and lysyl-oxidase enzymes.

PPARα Agonism Enhances Immune Response to Radiotherapy While Dietary Oleic Acid Results in Counteraction.

PURPOSE: Head and neck cancer (HNC) improvements are stagnant, even with advances in immunotherapy. Our previous clinical trial data show that altered fatty acid (FA) metabolism correlates with outcome. We hypothesized that pharmacologic and dietary modulation of FA catabolism will affect therapeutic efficacy. EXPERIMENTAL DESIGN: We performed in vivo and in vitro experiments using PPARα agonism with fenofibrate (FF) or high oleic acid diets (OAD) with radiotherapy, generating metabolomic, proteomic, stable isotope tracing, extracellular flux analysis, and flow-cytometric data to investigate these alterations. RESULTS: FF improved antitumor efficacy of high dose per fraction radiotherapy in HNC murine models, whereas the OAD reversed this effect. FF-treated mice on the control diet had evidence of increased FA catabolism. Stable isotope tracing showed less glycolytic utilization by ex vivo CD8+ T cells. Improved efficacy correlated with intratumoral alterations in eicosanoid metabolism and downregulated mTOR and CD36. CONCLUSIONS: Metabolic intervention with increased FA catabolism improves the efficacy of HNC therapy and enhances antitumoral immune response.

Decoding the molecular interplay in the central dogma: An overview of mass spectrometry-based methods to investigate protein-metabolite interactions.

With the emergence of next-generation nucleotide sequencing and mass spectrometry-based proteomics and metabolomics tools, we have comprehensive and scalable methods to analyze the genes, transcripts, proteins, and metabolites of a multitude of biological systems. Despite the fascinating new molecular insights at the genome, transcriptome, proteome and metabolome scale, we are still far from fully understanding cellular organization, cell cycles and biology at the molecular level. Significant advances in sensitivity and depth for both sequencing as well as mass spectrometry-based methods allow the analysis at the single cell and single molecule level. At the same time, new tools are emerging that enable the investigation of molecular interactions throughout the central dogma of molecular biology. In this review, we provide an overview of established and recently developed mass spectrometry-based tools to probe metabolite-protein interactions-from individual interaction pairs to interactions at the proteome-metabolome scale.

ECM-Focused Proteomic Analysis of Ear Punch Regeneration in Acomys Cahirinus.

In mammals, significant injury is generally followed by the formation of a fibrotic scar which provides structural integrity but fails to functionally restore damaged tissue. Spiny mice of the genus Acomys represent the first example of full skin autotomy in mammals. Acomys cahirinus has evolved extremely weak skin as a strategy to avoid predation and is able to repeatedly regenerate healthy tissue without scar after severe skin injury or full-thickness ear punches. Extracellular matrix (ECM) composition is a critical regulator of wound repair and scar formation and previous studies have suggested that alterations in its expression may be responsible for the differences in regenerative capacity observed between Mus musculus and A. cahirinus , yet analysis of this critical tissue component has been limited in previous studies by its insolubility and resistance to extraction. Here, we utilize a 2-step ECM-optimized extraction to perform proteomic analysis of tissue composition during wound repair after full-thickness ear punches in A. cahirinus and M. musculus from weeks 1 to 4 post-injury. We observe changes in a wide range of ECM proteins which have been previously implicated in wound regeneration and scar formation, including collagens, coagulation and provisional matrix proteins, and matricryptic signaling peptides. We additionally report differences in crosslinking enzyme activity and ECM protein solubility between Mus and Acomys. Furthermore, we observed rapid and sustained increases in CD206, a marker of pro-regenerative M2 macrophages, in Acomys, whereas little or no increase in CD206 was detected in Mus. Together, these findings contribute to a comprehensive understanding of tissue cues which drive the regenerative capacity of Acomys and identify a number of potential targets for future pro-regenerative therapies.

Col1a2-Deleted Mice Have Defective Type I Collagen and Secondary Reactive Cardiac Fibrosis with Altered Hypertrophic Dynamics.

RATIONALE: The adult cardiac extracellular matrix (ECM) is largely comprised of type I collagen. In addition to serving as the primary structural support component of the cardiac ECM, type I collagen also provides an organizational platform for other ECM proteins, matricellular proteins, and signaling components that impact cellular stress sensing in vivo. OBJECTIVE: Here we investigated how the content and integrity of type I collagen affect cardiac structure function and response to injury. METHODS AND RESULTS: We generated and characterized Col1a2-/- mice using standard gene targeting. Col1a2-/- mice were viable, although by young adulthood their hearts showed alterations in ECM mechanical properties, as well as an unanticipated activation of cardiac fibroblasts and induction of a progressive fibrotic response. This included augmented TGFß activity, increases in fibroblast number, and progressive cardiac hypertrophy, with reduced functional performance by 9 months of age. Col1a2-loxP-targeted mice were also generated and crossed with the tamoxifen-inducible Postn-MerCreMer mice to delete the Col1a2 gene in myofibroblasts with pressure overload injury. Interestingly, while germline Col1a2-/- mice showed gradual pathologic hypertrophy and fibrosis with aging, the acute deletion of Col1a2 from activated adult myofibroblasts showed a loss of total collagen deposition with acute cardiac injury and an acute reduction in pressure overload-induce cardiac hypertrophy. However, this reduction in hypertrophy due to myofibroblast-specific Col1a2 deletion was lost after 2 and 6 weeks of pressure overload, as fibrotic deposition accumulated. CONCLUSIONS: Defective type I collagen in the heart alters the structural integrity of the ECM and leads to cardiomyopathy in adulthood, with fibroblast expansion, activation, and alternate fibrotic ECM deposition. However, acute inhibition of type I collagen production can have an anti-fibrotic and anti-hypertrophic effect.

From volcanoes to the bench: Advantages of novel hyperthermoacidic archaeal proteases for proteomics workflows.

Here we introduce hyperthermoacidic archaeal proteases (HTA-Proteases©) isolated from organisms that thrive in nearly boiling acidic volcanic springs and investigate their use for bottom-up proteomic experiments. We find that HTA-Proteases have novel cleavage specificities, show no autolysis, function in dilute formic acid, and store at ambient temperature for years. HTA-Proteases function optimally at 70-90 °C and pH of 2-4 with rapid digestion kinetics. The extreme HTA-Protease reaction conditions actively denature sample proteins, obviate the use of chaotropes, are largely independent of reduction and alkylation, and allow for a one-step/five-minute sample preparation protocol without sample manipulation, dilution, or additional cleanup. We find that brief one-step HTA-Protease protocols significantly increase proteome and protein sequence coverage with datasets orthogonal to trypsin. Importantly, HTA-Protease digests markedly increase coverage and identifications for ribonucleoproteins, histones, and mitochondrial membrane proteins as compared to tryptic digests alone. In addition to increased coverage in these classes, HTA-Proteases and simplified one-step protocols are expected to reduce technical variability and advance the fields of clinical and high-throughput proteomics. This work reveals significant utility of heretofore unavailable HTA-Proteases for proteomic workflows. We discuss some of the potential for these remarkable enzymes to empower new proteomics methods, approaches, and biological insights. SIGNIFICANCE: Here we introduce new capabilities for bottom-up proteomics applications with hyperthermoacidic archaeal proteases (HTA-Proteases©). HTA-Proteases have novel cleavage specificity, require no chaotropes, and allow simple one-step/five-minute sample preparations that promise to reduce variability between samples and laboratories. HTA-Proteases generate unique sets of observable peptides that are non-overlapping with tryptic peptides and significantly increase sequence coverage and available peptide targets relative to trypsin alone. HTA-Proteases show some bias for the detection and coverage of nucleic acid-binding proteins and membrane proteins relative to trypsin. These new ultra-stable enzymes function optimally in nearly boiling acidic conditions, show no autolysis, and do not require aliquoting as they are stable for years at ambient temperatures. Used independently or in conjunction with tryptic digests, HTA-Proteases offer increased proteome coverage, unique peptide targets, and brief one-step protocols amenable to automation, rapid turnaround, and high-throughput approaches.

Luteal Lipid Droplets: A Novel Platform for Steroid Synthesis.

Progesterone is an essential steroid hormone that is required to initiate and maintain pregnancy in mammals and serves as a metabolic intermediate in the synthesis of endogenously produced steroids, including sex hormones and corticosteroids. Steroidogenic luteal cells of the corpus luteum have the tremendous capacity to synthesize progesterone. These specialized cells are highly enriched with lipid droplets that store lipid substrate, which can be used for the synthesis of steroids. We recently reported that hormone-stimulated progesterone synthesis by luteal cells requires protein kinase A-dependent mobilization of cholesterol substrate from lipid droplets to mitochondria. We hypothesize that luteal lipid droplets are enriched with steroidogenic enzymes and facilitate the synthesis of steroids in the corpus luteum. In the present study, we analyzed the lipid droplet proteome, conducted the first proteomic analysis of lipid droplets under acute cyclic adenosine monophosphate (cAMP)-stimulated conditions, and determined how specific lipid droplet proteins affect steroidogenesis. Steroidogenic enzymes, cytochrome P450 family 11 subfamily A member 1 and 3 beta-hydroxysteroid dehydrogenase (HSD3B), were highly abundant on lipid droplets of the bovine corpus luteum. High-resolution confocal microscopy confirmed the presence of active HSD3B on the surface of luteal lipid droplets. We report that luteal lipid droplets have the capacity to synthesize progesterone from pregnenolone. Lastly, we analyzed the lipid droplet proteome following acute stimulation with cAMP analog, 8-Br-cAMP, and report increased association of HSD3B with luteal lipid droplets following stimulation. These findings provide novel insights into the role of luteal lipid droplets in steroid synthesis.

Selective targeting of IL2Rßγ combined with radiotherapy triggers CD8- and NK-mediated immunity, abrogating metastasis in HNSCC.

The implementation of cancer immunotherapies has seen limited clinical success in head and neck squamous cell carcinoma (HNSCC). Interleukin-2 (IL-2), which modulates the survival and functionality of lymphocytes, is an attractive target for new immunotherapies but one that is limited by presence of regulatory T cells (Tregs) expressing the high-affinity IL-2Rα. The bispecific immunocytokine PD1-IL2v preferentially delivers IL-2 signaling through IL-2Rßγ on PD-1-expressing cells. Selectively targeting the intermediate-affinity IL-2Rßγ can be leveraged to induce anti-tumor immune responses in effector T cells and natural killer (NK) cells while limiting the negative regulation of IL-2Rα activation on Tregs. Using radiation therapy (RT) in combination with PD1-IL2v improves local tumor control and survival, and controls metastatic spread in orthotopic HNSCC tumor models. PD1-IL2v drives systemic activation and expansion of circulating and tumor-infiltrating cytotoxic T cells and NK cells while limiting Treg-mediated immunosuppression. These data show that PD1-L2v induces durable systemic tumor control in HNSCC.

Proregenerative extracellular matrix hydrogel mitigates pathological alterations of pelvic skeletal muscles after birth injury.

Pelvic floor disorders, including pelvic organ prolapse and urinary and fecal incontinence, affect millions of women globally and represent a major public health concern. Pelvic floor muscle (PFM) dysfunction has been identified as one of the leading risk factors for the development of these morbid conditions. Childbirth, specifically vaginal delivery, has been recognized as the most important potentially modifiable risk factor for PFM injury; however, the precise mechanisms of PFM dysfunction after parturition remain elusive. In this study, we demonstrated that PFMs exhibit atrophy and fibrosis in parous women with symptomatic pelvic organ prolapse. These pathological alterations were recapitulated in a preclinical rat model of simulated birth injury (SBI). The transcriptional signature of PFMs after injury demonstrated an impairment in muscle anabolism, persistent expression of genes that promote extracellular matrix (ECM) deposition, and a sustained inflammatory response. We also evaluated the administration of acellular injectable skeletal muscle ECM hydrogel for the prevention of these pathological alterations. Treatment of PFMs with the ECM hydrogel either at the time of birth injury or 4 weeks after injury mitigated PFM atrophy and fibrosis. By evaluating gene expression, we demonstrated that these changes are mainly driven by the hydrogel-induced enhancement of endogenous myogenesis, ECM remodeling, and modulation of the immune response. This work furthers our understanding of PFM birth injury and demonstrates proof of concept for future investigations of proregenerative biomaterial approaches for the treatment of injured pelvic soft tissues.

The best of both worlds? Rattlesnake hybrid zones generate complex combinations of divergent venom phenotypes that retain high toxicity.

Studying the consequences of hybridization between closely related species with divergent traits can reveal patterns of evolution that shape and maintain extreme trophic adaptations. Snake venoms are an excellent model system for examining the evolutionary and ecological patterns that underlie highly selected polymorphic traits. Here we investigate hybrid venom phenotypes that result from natural introgression between two rattlesnake species that express highly divergent venom phenotypes: Crotalus o. concolor and C. v. viridis. Though not yet documented, interbreeding between these species may lead to novel venom phenotypes with unique activities that break the typical trends of venom composition in rattlesnakes. The characteristics of these unusual phenotypes could unveil the roles of introgression in maintaining patterns of venom composition and variation, including the near ubiquitous dichotomy between neurotoxic or degradative venoms observed across rattlesnakes. We use RADseq data to infer patterns of gene flow and hybrid ancestry between these diverged lineages and link these genetic data with analyses of venom composition, biological activity, and whole animal model toxicity tests to understand the impacts of introgression on venom composition. We find that introgressed populations express admixed venom phenotypes that do not sacrifice biological activity (lethal toxicity) or overall abundance of dominant toxins compared to parental venoms. These hybridized venoms therefore do not represent a trade-off in functionality between the typical phenotypic extremes but instead represent a unique combination of characters whose expression appears limited to the hybrid zone.

Single-Cell Heterogeneity in Snake Venom Expression Is Hardwired by Co-Option of Regulators from Progressively Activated Pathways.

The ubiquitous cellular heterogeneity underlying many organism-level phenotypes raises questions about what factors drive this heterogeneity and how these complex heterogeneous systems evolve. Here, we use single-cell expression data from a Prairie rattlesnake (Crotalus viridis) venom gland to evaluate hypotheses for signaling networks underlying snake venom regulation and the degree to which different venom gene families have evolutionarily recruited distinct regulatory architectures. Our findings suggest that snake venom regulatory systems have evolutionarily co-opted trans-regulatory factors from extracellular signal-regulated kinase and unfolded protein response pathways that specifically coordinate expression of distinct venom toxins in a phased sequence across a single population of secretory cells. This pattern of co-option results in extensive cell-to-cell variation in venom gene expression, even between tandemly duplicated paralogs, suggesting this regulatory architecture has evolved to circumvent cellular constraints. While the exact nature of such constraints remains an open question, we propose that such regulatory heterogeneity may circumvent steric constraints on chromatin, cellular physiological constraints (e.g., endoplasmic reticulum stress or negative protein-protein interactions), or a combination of these. Regardless of the precise nature of these constraints, this example suggests that, in some cases, dynamic cellular constraints may impose previously unappreciated secondary constraints on the evolution of gene regulatory networks that favors heterogeneous expression.

Snakes on a plain: biotic and abiotic factors determine venom compositional variation in a wide-ranging generalist rattlesnake.

BACKGROUND: Snake venoms are trophic adaptations that represent an ideal model to examine the evolutionary factors that shape polymorphic traits under strong natural selection. Venom compositional variation is substantial within and among venomous snake species. However, the forces shaping this phenotypic complexity, as well as the potential integrated roles of biotic and abiotic factors, have received little attention. Here, we investigate geographic variation in venom composition in a wide-ranging rattlesnake (Crotalus viridis viridis) and contextualize this variation by investigating dietary, phylogenetic, and environmental variables that covary with venom. RESULTS: Using shotgun proteomics, venom biochemical profiling, and lethality assays, we identify 2 distinct divergent phenotypes that characterize major axes of venom variation in this species: a myotoxin-rich phenotype and a snake venom metalloprotease (SVMP)-rich phenotype. We find that dietary availability and temperature-related abiotic factors are correlated with geographic trends in venom composition. CONCLUSIONS: Our findings highlight the potential for snake venoms to vary extensively within species, for this variation to be driven by biotic and abiotic factors, and for the importance of integrating biotic and abiotic variation for understanding complex trait evolution. Links between venom variation and variation in biotic and abiotic factors indicate that venom variation likely results from substantial geographic variation in selection regimes that determine the efficacy of venom phenotypes across populations and snake species. Our results highlight the cascading influence of abiotic factors on biotic factors that ultimately shape venom phenotype, providing evidence for a central role of local selection as a key driver of venom variation.

Quantification of snake venom proteomes by mass spectrometry-considerations and perspectives.

The advent of soft ionization mass spectrometry-based proteomics in the 1990s led to the development of a new dimension in biology that conceptually allows for the integral analysis of whole proteomes. This transition from a reductionist to a global-integrative approach is conditioned to the capability of proteomic platforms to generate and analyze complete qualitative and quantitative proteomics data. Paradoxically, the underlying analytical technique, molecular mass spectrometry, is inherently nonquantitative. The turn of the century witnessed the development of analytical strategies to endow proteomics with the ability to quantify proteomes of model organisms in the sense of "an organism for which comprehensive molecular (genomic and/or transcriptomic) resources are available." This essay presents an overview of the strategies and the lights and shadows of the most popular quantification methods highlighting the common misuse of label-free approaches developed for model species' when applied to quantify the individual components of proteomes of nonmodel species (In this essay we use the term "non-model" organisms for species lacking comprehensive molecular (genomic and/or transcriptomic) resources, a circumstance that, as we detail in this review-essay, conditions the quantification of their proteomes.). We also point out the opportunity of combining elemental and molecular mass spectrometry systems into a hybrid instrumental configuration for the parallel identification and absolute quantification of venom proteomes. The successful application of this novel mass spectrometry configuration in snake venomics represents a proof-of-concept for a broader and more routine application of hybrid elemental/molecular mass spectrometry setups in other areas of the proteomics field, such as phosphoproteomics, metallomics, and in general in any biological process where a heteroatom (i.e., any atom other than C, H, O, N) forms integral part of its mechanism.