Mutual enlightenment: A toolbox of concepts and methods for integrating evolutionary and clinical toxinology via snake venomics and the contextual stance

Snakebite envenoming is a neglected tropical disease that may claim over 100,000 human lives annually worldwide. Snakebite occurs as the result of an interaction between a human and a snake that elicits either a defensive response from the snake or, more rarely, a feeding response as the result of mistaken identity. Snakebite envenoming is therefore a biological and, more specifically, an ecological problem. Snake venom itself is often described as a “cocktail”, as it is a heterogenous mixture of molecules including the toxins (which are typically proteinaceous) responsible for the pathophysiological consequences of envenoming. The primary function of venom in snake ecology is pre-subjugation, with defensive deployment of the secretion typically considered a secondary function. The particular composition of any given venom cocktail is shaped by evolutionary forces that include phylogenetic constraints associated with the snake’s lineage and adaptive responses to the snake’s ecological context, including the taxa it preys upon and by which it is predated upon. In the present article, we describe how conceptual frameworks from ecology and evolutionary biology can enter into a mutually enlightening relationship with clinical toxinology by enabling the consideration of snakebite envenoming from an “ecological stance”. We detail the insights that may emerge from such a perspective and highlight the ways in which the high-fidelity descriptive knowledge emerging from applications of -omics era technologies – “venomics” and “antivenomics” – can combine with evolutionary explanations to deliver a detailed understanding of this multifactorial health crisis.Ministerio de Ciencia e Innovacion/[BMC 2004-01432]//EspañaMinisterio de Ciencia e Innovacion/[BFU 2007-61563]//EspañaMinisterio de Ciencia e Innovacion/[BFU 2010-173730]//EspañaMinisterio de Ciencia e Innovacion/[BFU 2013-42833-P]//EspañaMinisterio de Ciencia e Innovacion/[BFU 2017-89103-P]//EspañaNorwegian Research Council/[No.287462.]/NFR/NoruegaNational Health and Medical Research Council/[Grant 13/093/002 AVRU]/AustraliaDBT/Wellcome Trust India Alliance/[IA/I/19/2/504647]//IndiaUCR::Vicerrectoría de Investigación::Unidades de Investigación::Ciencias de la Salud::Instituto Clodomiro Picado (ICP

Studying Venom Toxin Variation Using Accurate Masses from Liquid Chromatography–Mass Spectrometry Coupled with Bioinformatic Tools

This study provides a new methodology for the rapid analysis of numerous venom samples in an automated fashion. Here, we use LC-MS (Liquid Chromatography–Mass Spectrometry) for venom separation and toxin analysis at the accurate mass level combined with new in-house written bioinformatic scripts to obtain high-throughput results. This analytical methodology was validated using 31 venoms from all members of a monophyletic clade of Australian elapids: brown snakes (Pseudonaja spp.) and taipans (Oxyuranus spp.). In a previous study, we revealed extensive venom variation within this clade, but the data was manually processed and MS peaks were integrated into a time-consuming and labour-intensive approach. By comparing the manual approach to our new automated approach, we now present a faster and more efficient pipeline for analysing venom variation. Pooled venom separations with post-column toxin fractionations were performed for subsequent high-throughput venomics to obtain toxin IDs correlating to accurate masses for all fractionated toxins. This workflow adds another dimension to the field of venom analysis by providing opportunities to rapidly perform in-depth studies on venom variation. Our pipeline opens new possibilities for studying animal venoms as evolutionary model systems and investigating venom variation to aid in the development of better antivenoms

Association of Blood Biomarkers With Acute Sport-Related Concussion in Collegiate Athletes: Findings From the NCAA and Department of Defense CARE Consortium

Importance: There is potential scientific and clinical value in validation of objective biomarkers for sport-related concussion (SRC). Objective: To investigate the association of acute-phase blood biomarker levels with SRC in collegiate athletes. Design, Setting, and Participants: This multicenter, prospective, case-control study was conducted by the National Collegiate Athletic Association (NCAA) and the US Department of Defense Concussion Assessment, Research, and Education (CARE) Consortium from February 20, 2015, to May 31, 2018, at 6 CARE Advanced Research Core sites. A total of 504 collegiate athletes with concussion, contact sport control athletes, and non-contact sport control athletes completed clinical testing and blood collection at preseason baseline, the acute postinjury period, 24 to 48 hours after injury, the point of reporting being asymptomatic, and 7 days after return to play. Data analysis was conducted from March 1 to November 30, 2019. Main Outcomes and Measures: Glial fibrillary acidic protein (GFAP), ubiquitin C-terminal hydrolase-L1 (UCH-L1), neurofilament light chain, and tau were quantified using the Quanterix Simoa multiplex assay. Clinical outcome measures included the Sport Concussion Assessment Tool-Third Edition (SCAT-3) symptom evaluation, Standardized Assessment of Concussion, Balance Error Scoring System, and Brief Symptom Inventory 18. Results: A total of 264 athletes with concussion (mean [SD] age, 19.08 [1.24] years; 211 [79.9%] male), 138 contact sport controls (mean [SD] age, 19.03 [1.27] years; 107 [77.5%] male), and 102 non-contact sport controls (mean [SD] age, 19.39 [1.25] years; 82 [80.4%] male) were included in the study. Athletes with concussion had significant elevation in GFAP (mean difference, 0.430 pg/mL; 95% CI, 0.339-0.521 pg/mL; P < .001), UCH-L1 (mean difference, 0.449 pg/mL; 95% CI, 0.167-0.732 pg/mL; P < .001), and tau levels (mean difference, 0.221 pg/mL; 95% CI, 0.046-0.396 pg/mL; P = .004) at the acute postinjury time point compared with preseason baseline. Longitudinally, a significant interaction (group × visit) was found for GFAP (F7,1507.36 = 16.18, P < .001), UCH-L1 (F7,1153.09 = 5.71, P < .001), and tau (F7,1480.55 = 6.81, P < .001); the interaction for neurofilament light chain was not significant (F7,1506.90 = 1.33, P = .23). The area under the curve for the combination of GFAP and UCH-L1 in differentiating athletes with concussion from contact sport controls at the acute postinjury period was 0.71 (95% CI, 0.64-0.78; P < .001); the acute postinjury area under the curve for all 4 biomarkers combined was 0.72 (95% CI, 0.65-0.79; P < .001). Beyond SCAT-3 symptom score, GFAP at the acute postinjury time point was associated with the classification of athletes with concussion from contact controls (β = 12.298; 95% CI, 2.776-54.481; P = .001) and non-contact sport controls (β = 5.438; 95% CI, 1.676-17.645; P = .005). Athletes with concussion with loss of consciousness or posttraumatic amnesia had significantly higher levels of GFAP than athletes with concussion with neither loss of consciousness nor posttraumatic amnesia at the acute postinjury time point (mean difference, 0.583 pg/mL; 95% CI, 0.369-0.797 pg/mL; P < .001). Conclusions and Relevance: The results suggest that blood biomarkers can be used as research tools to inform the underlying pathophysiological mechanism of concussion and provide additional support for future studies to optimize and validate biomarkers for potential clinical use in SRC

Thrombosis Is Reduced by Inhibition of COX-1, but Unaffected by Inhibition of COX-2, in an Acute Model of Platelet Activation in the Mouse

This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

COMPLEMENT FACTOR B IS A DETERMINANT OF BOTH METABOLIC AND CARDIOVASCULAR FEATURES OF METABOLIC SYNDROME

CFB (complement factor B) is elevated in adipose tissue and serum from patients with type 2 diabetes mellitus and cardiovascular disease, but the causal relationship to disease pathogenesis is unclear. Cfb is also elevated in adipose tissue and serum of the spontaneously hypertensive rat, a well-characterized model of metabolic syndrome. To establish the role of CFB in metabolic syndrome, we knocked out the Cfb gene in the spontaneously hypertensive rat. Cfb−/− rats showed improved glucose tolerance and insulin sensitivity, redistribution of visceral to subcutaneous fat, increased adipocyte mitochondrial respiration, and marked changes in gene expression. Cfb−/− rats also had lower blood pressure, increased ejection fraction and fractional shortening, and reduced left ventricular mass. These changes in metabolism and gene expression, in adipose tissue and left ventricle, suggest new adipose tissue-intrinsic and blood pressure-independent mechanisms for insulin resistance and cardiac hypertrophy in the spontaneously hypertensive rat. In silico analysis of the human CFB locus revealed 2 cis-regulated expression quantitative trait loci for CFB expression significantly associated with visceral fat, circulating triglycerides and hypertension in genome-wide association studies. Together, these data demonstrate a key role for CFB in the development of spontaneously hypertensive rat metabolic syndrome phenotypes and of related traits in humans and indicate the potential for CFB as a novel target for treatment of cardiometabolic disease

Catch a tiger snake by its tail: Differential toxicity, co-factor dependence and antivenom efficacy in a procoagulant clade of Australian venomous snakes

A paradigm of venom research is adaptive evolution of toxins as part of a predator-prey chemical arms race. This study examined differential co-factor dependence, variations relative to dietary preference, and the impact upon relative neutralisation by antivenom of the procoagulant toxins in the venoms of a clade of Australian snakes. All genera were characterised by venoms rich in factor Xa which act upon endogenous prothrombin. Examination of toxin sequences revealed an extraordinary level of conservation, which indicates that adaptive evolution is not a feature of this toxin type. Consistent with this, the venoms did not display differences on the plasma of different taxa. Examination of the prothrombin target revealed endogenous blood proteins are under extreme negative selection pressure for diversification, this in turn puts a strong negative selection pressure upon the toxins as sequence diversification could result in a drift away from the target. Thus this study reveals that adaptive evolution is not a consistent feature in toxin evolution in cases where the target is under negative selection pressure for diversification. Consistent with this high level of toxin conservation, the antivenom showed extremely high-levels of cross-reactivity. There was however a strong statistical correlation between relative degree of phospholipid-dependence and clotting time, with the least dependent venoms producing faster clotting times than the other venoms even in the presence of phospholipid. The results of this study are not only of interest to evolutionary and ecological disciplines, but also have implications for clinical toxinology

Highly Evolvable: Investigating Interspecific and Intraspecific Venom Variation in Taipans (Oxyuranus spp.) and Brown Snakes (Pseudonaja spp.)

Snake venoms are complex mixtures of toxins that differ on interspecific (between species) and intraspecific (within species) levels. Whether venom variation within a group of closely related species is explained by the presence, absence and/or relative abundances of venom toxins remains largely unknown. Taipans (Oxyuranus spp.) and brown snakes (Pseudonaja spp.) represent medically relevant species of snakes across the Australasian region and provide an excellent model clade for studying interspecific and intraspecific venom variation. Using liquid chromatography with ultraviolet and mass spectrometry detection, we analyzed a total of 31 venoms covering all species of this monophyletic clade, including widespread localities. Our results reveal major interspecific and intraspecific venom variation in Oxyuranus and Pseudonaja species, partially corresponding with their geographical regions and phylogenetic relationships. This extensive venom variability is generated by a combination of the absence/presence and differential abundance of venom toxins. Our study highlights that venom systems can be highly dynamical on the interspecific and intraspecific levels and underscores that the rapid toxin evolvability potentially causes major impacts on neglected tropical snakebites

Distinct regulatory networks control toxin gene expression in elapid and viperid snakes

Background Venom systems are ideal models to study genetic regulatory mechanisms that underpin evolutionary novelty. Snake venom glands are thought to share a common origin, but there are major distinctions between venom toxins from the medically significant snake families Elapidae and Viperidae, and toxin gene regulatory investigations in elapid snakes have been limited. Here, we used high-throughput RNA-sequencing to profile gene expression and microRNAs between active (milked) and resting (unmilked) venom glands in an elapid (Eastern Brown Snake, Pseudonaja textilis), in addition to comparative genomics, to identify cis- and trans-acting regulation of venom production in an elapid in comparison to viperids (Crotalus viridis and C. tigris). Results Although there is conservation in high-level mechanistic pathways regulating venom production (unfolded protein response, Notch signaling and cholesterol homeostasis), there are differences in the regulation of histone methylation enzymes, transcription factors, and microRNAs in venom glands from these two snake families. Histone methyltransferases and transcription factor (TF) specificity protein 1 (Sp1) were highly upregulated in the milked elapid venom gland in comparison to the viperids, whereas nuclear factor I (NFI) TFs were upregulated after viperid venom milking. Sp1 and NFI cis-regulatory elements were common to toxin gene promoter regions, but many unique elements were also present between elapid and viperid toxins. The presence of Sp1 binding sites across multiple elapid toxin gene promoter regions that have been experimentally determined to regulate expression, in addition to upregulation of Sp1 after venom milking, suggests this transcription factor is involved in elapid toxin expression. microRNA profiles were distinctive between milked and unmilked venom glands for both snake families, and microRNAs were predicted to target a diversity of toxin transcripts in the elapid P. textilis venom gland, but only snake venom metalloproteinase transcripts in the viperid C. viridis venom gland. These results suggest differences in toxin gene posttranscriptional regulation between the elapid P. textilis and viperid C. viridis. Conclusions Our comparative transcriptomic and genomic analyses between toxin genes and isoforms in elapid and viperid snakes suggests independent toxin regulation between these two snake families, demonstrating multiple different regulatory mechanisms underpin a venomous phenotype

The Milky Way's circular velocity curve between 4 and 14 kpc from APOGEE data

We measure the Milky Way's rotation curve over the Galactocentric range 4 kpc <~ R <~ 14 kpc from the first year of data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE). We model the line-of-sight velocities of 3,365 stars in fourteen fields with b = 0 deg between 30 deg < l < 210 deg out to distances of 10 kpc using an axisymmetric kinematical model that includes a correction for the asymmetric drift of the warm tracer population (\sigma_R ~ 35 km/s). We determine the local value of the circular velocity to be V_c(R_0) = 218 +/- 6 km/s and find that the rotation curve is approximately flat with a local derivative between -3.0 km/s/kpc and 0.4 km/s/kpc. We also measure the Sun's position and velocity in the Galactocentric rest frame, finding the distance to the Galactic center to be 8 kpc < R_0 < 9 kpc, radial velocity V_{R,sun} = -10 +/- 1 km/s, and rotational velocity V_{\phi,sun} = 242^{+10}_{-3} km/s, in good agreement with local measurements of the Sun's radial velocity and with the observed proper motion of Sgr A*. We investigate various systematic uncertainties and find that these are limited to offsets at the percent level, ~2 km/s in V_c. Marginalizing over all the systematics that we consider, we find that V_c(R_0) 99% confidence. We find an offset between the Sun's rotational velocity and the local circular velocity of 26 +/- 3 km/s, which is larger than the locally-measured solar motion of 12 km/s. This larger offset reconciles our value for V_c with recent claims that V_c >~ 240 km/s. Combining our results with other data, we find that the Milky Way's dark-halo mass within the virial radius is ~8x10^{11} M_sun.Comment: submitted to Ap

Correlation between ontogenetic dietary shifts and venom variation in Australian brown snakes ( Pseudonaja )

Venom is a key evolutionary trait, as evidenced by its widespread convergent evolution across the animal kingdom. In an escalating prey-predator arms race, venoms evolve rapidly to guarantee predatory or defensive success. Variation in venom composition is ubiquitous among snakes. Here, we tested variation in venom activity on substrates relevant to blood coagulation among Pseudonaja (brown snake) species, Australian elapids responsible for the majority of medically important human envenomations in Australia. A functional approach was employed to elucidate interspecific variation in venom activity in all nine currently recognised species of Pseudonaja. Fluorometric enzymatic activity assays were performed to test variation in whole venom procoagulant activity among species. Analyses confirmed the previously documented ontogenetic shift from non-coagulopathic venom in juveniles to coagulopathic venom as adults, except for the case of P. modesta, which retains non-coagulopathic venom as an adult. These shifts in venom activity correlate with documented ontogenetic shifts in diet among brown snakes from specialisation on reptilian prey as juveniles (and throughout the life cycle of P. modesta), to a more generalised diet in adults that includes mammals. The results of this study bring to light findings relevant to both clinical and evolutionary toxinology