Phylogeography of the cottonmouth, Agkistrodon piscivorus, using AFLP and venom protein profiles.

The objective of this study was to examine population structure in cottonmouths (Agkistrodon piscivorus) using Amplified Fragment Length Polymorphism (AFLP) and compare genetic and venom protein profiles in Texas. AFLP profiles using 622 fragments were generated for 105 individuals to understand the level of variation within Agkistrodon. In Texas, there was a significant lack of gene flow detected and support for the isolation of Concho Valley individuals. Cottonmouths showed the greatest genetic variation when compared to other Agkistrodon species but there was not complete support for two species of cottonmouths as currently proposed. RP-HPLC was used to examine venom protein profiles in 86 Texas cottonmouths. Relative peak heights were analyzed using PCA and the MANOVA demonstrated separation of populations based on profiles (p<0.001). Genetic and venom variation did not follow the same pattern indicating that there may be other selection pressures acting on the venom proteins

Trait differentiation and modular toxin expression in palm-pitvipers

Background Modularity is the tendency for systems to organize into semi-independent units and can be a key to the evolution and diversification of complex biological systems. Snake venoms are highly variable modular systems that exhibit extreme diversification even across very short time scales. One well-studied venom phenotype dichotomy is a trade-off between neurotoxicity versus hemotoxicity that occurs through the high expression of a heterodimeric neurotoxic phospholipase A2 (PLA2) or snake venom metalloproteinases (SVMPs). We tested whether the variation in these venom phenotypes could occur via variation in regulatory sub-modules through comparative venom gland transcriptomics of representative Black-Speckled Palm-Pitvipers (Bothriechis nigroviridis) and Talamancan Palm-Pitvipers (B. nubestris). Results We assembled 1517 coding sequences, including 43 toxins for B. nigroviridis and 1787 coding sequences including 42 toxins for B. nubestris. The venom gland transcriptomes were extremely divergent between these two species with one B. nigroviridis exhibiting a primarily neurotoxic pattern of expression, both B. nubestris expressing primarily hemorrhagic toxins, and a second B. nigroviridis exhibiting a mixed expression phenotype. Weighted gene coexpression analyses identified six submodules of transcript expression variation, one of which was highly associated with SVMPs and a second which contained both subunits of the neurotoxic PLA2 complex. The sub-module association of these toxins suggest common regulatory pathways underlie the variation in their expression and is consistent with known patterns of inheritance of similar haplotypes in other species. We also find evidence that module associated toxin families show fewer gene duplications and transcript losses between species, but module association did not appear to affect sequence diversification. Conclusion Sub-modular regulation of expression likely contributes to the diversification of venom phenotypes within and among species and underscores the role of modularity in facilitating rapid evolution of complex traits

Comparative venom-gland transcriptomics and venom proteomics of four Sidewinder Rattlesnake (\u3ci\u3eCrotalus cerastes\u3c/i\u3e) lineages reveal little differential expression despite individual variation

Changes in gene expression can rapidly influence adaptive traits in the early stages of lineage diversification. Venom is an adaptive trait comprised of numerous toxins used for prey capture and defense. Snake venoms can vary widely between conspecific populations, but the influence of lineage diversification on such compositional differences are unknown. To explore venom differentiation in the early stages of lineage diversification, we used RNA-seq and mass spectrometry to characterize Sidewinder Rattlesnake (Crotalus cerastes) venom. We generated the first venom-gland transcriptomes and complementary venom proteomes for eight individuals collected across the United States and tested for expression differences across life history traits and between subspecific, mitochondrial, and phylotranscriptomic hypotheses. Sidewinder venom was comprised primarily of hemorrhagic toxins, with few cases of differential expression attributable to life history or lineage hypotheses. However, phylotranscriptomic lineage comparisons more than doubled instances of significant expression differences compared to all other factors. Nevertheless, only 6.4% of toxins were differentially expressed overall, suggesting that shallow divergence has not led to major changes in Sidewinder venom composition. Our results demonstrate the need for consensus venom-gland transcriptomes based on multiple individuals and highlight the potential for discrepancies in differential expression between different phylogenetic hypotheses

Schwann Cell O-GlcNAc Glycosylation Is Required for Myelin Maintenance and Axon Integrity

Schwann cells (SCs), ensheathing glia of the peripheral nervous system, support axonal survival and function. Abnormalities in SC metabolism affect their ability to provide this support and maintain axon integrity. To further interrogate this metabolic influence on axon–glial interactions, we generated OGT-SCKO mice with SC-specific deletion of the metabolic/nutrient sensing protein O-GlcNAc transferase that mediates the O-linked addition of N-acetylglucosamine (GlcNAc) moieties to Ser and Thr residues. The OGT-SCKO mice develop tomaculous demyelinating neuropathy characterized by focal thickenings of the myelin sheath (tomacula), progressive demyelination, axonal loss, and motor and sensory nerve dysfunction. Proteomic analysis identified more than 100 O-GlcNAcylated proteins in rat sciatic nerve, including Periaxin (PRX), a myelin protein whose mutation causes inherited neuropathy in humans. PRX lacking O-GlcNAcylation is mislocalized within the myelin sheath of these mutant animals. Furthermore, phenotypes of OGT-SCKO and Prx-deficient mice are very similar, suggesting that metabolic control of PRX O-GlcNAcylation is crucial for myelin maintenance and axonal integrity. SIGNIFICANCE STATEMENT The nutrient sensing protein O-GlcNAc transferase (OGT) mediates post-translational O-linked N-acetylglucosamine (GlcNAc) modification. Here we find that OGT functions in Schwann cells (SCs) to maintain normal myelin and prevent axonal loss. SC-specific deletion of OGT (OGT-SCKO mice) causes a tomaculous demyelinating neuropathy accompanied with progressive axon degeneration and motor and sensory nerve dysfunction. We also found Periaxin (PRX), a myelin protein whose mutation causes inherited neuropathy in humans, is O-GlcNAcylated. Importantly, phenotypes of OGT-SCKO and Prx mutant mice are very similar, implying that compromised PRX function contributes to the neuropathy of OGT-SCKO mice. This study will be useful in understanding how SC metabolism contributes to PNS function and in developing new strategies for treating peripheral neuropathy by targeting SC function

Intraspecific venom variation of Mexican West Coast Rattlesnakes (Crotalus basiliscus) and its implications for antivenom production

14 páginas, 9 figuras, 3 tablasIntraspecific variation in snake venoms has been widely documented worldwide. However, there are few studies on this subject in Mexico. Venom characterization studies provide important data used to predict clinical syndromes, to evaluate the efficacy of antivenoms and, in some cases, to improve immunogenic mixtures in the production of antivenoms. In the present work, we evaluated the intraspecific venom variation of Crotalus basiliscus, a rattlesnake of medical importance and whose venom is used in the immunization of horses to produce one of the Mexican antivenoms. Our results demonstrate that there is variation in biological and biochemical activities among adult venoms and that there is an ontogenetic change from juvenile to adult venoms. Juvenile venoms were more lethal and had higher percentages of crotamine and crotoxin, while adult venoms had higher percentages of snake venom metalloproteases (SVMPs). Additionally, we documented crotoxin-like PLA2 variation in which specimens from Zacatecas, Sinaloa and Michoacán (except 1) lacked the neurotoxin, while the rest of the venoms had it. Finally, we evaluated the efficacy of three lots of Birmex antivenom and all three were able to neutralize the lethality of four representative venoms but were not able to neutralize crotamine. We also observed significant differences in the LD50 values neutralized per vial among the different lots. Based on these results, we recommend including venoms containing crotamine in the production of antivenom for a better immunogenic mixture and to improve the homogeneity of lots.This study was financially supported by DGAPA-PAPIIT (project IN211621), CONACYT (project264255); FORDECYT PRONACE (project 1715618/2020), FORDECYT (project 303045), Clemson University and the National Science Foundation (DEB 1822417) to CLP.Peer reviewe

Evidence for divergent patterns of local selection driving venom variation in Mojave Rattlesnakes (\u3ci\u3eCrotalus scutulatus\u3c/i\u3e)

Snake venoms represent an enriched system for investigating the evolutionary processes that lead to complex and dynamic trophic adaptations. It has long been hypothesized that natural selection may drive geographic variation in venom composition, yet previous studies have lacked the population genetic context to examine these patterns. We leverage range-wide sampling of Mojave Rattlesnakes (Crotalus scutulatus) and use a combination of venom, morphological, phylogenetic, population genetic, and environmental data to characterize the striking dichotomy of neurotoxic (Type A) and hemorrhagic (Type B) venoms throughout the range of this species. We find that three of the four previously identified major lineages within C. scutulatus possess a combination of Type A, Type B, and a ‘mixed’ Type A + B venom phenotypes, and that fixation of the two main venom phenotypes occurs on a more fine geographic scale than previously appreciated. We also find that Type A + B individuals occur in regions of inferred introgression, and that this mixed phenotype is comparatively rare. Our results support strong directional local selection leading to fixation of alternative venom phenotypes on a fine geographic scale, and are inconsistent with balancing selection to maintain both phenotypes within a single population. Our comparisons to biotic and abiotic factors further indicate that venom phenotype correlates with fang morphology and climatic variables. We hypothesize that links to fang morphology may be indicative of co-evolution of venom and other trophic adaptations, and that climatic variables may be linked to prey distributions and/or physiology, which in turn impose selection pressures on snake venoms

Genomic Adaptations to Salinity Resist Gene Flow in the Evolution of Floridian Watersnakes

The migration-selection balance often governs the evolution of lineages, and speciation with gene flow is now considered common across the tree of life. Ecological speciation is a process that can facilitate divergence despite gene flow due to strong selective pressures caused by ecological differences; however, the exact traits under selection are often unknown. The transition from freshwater to saltwater habitats provides strong selection targeting traits with osmoregulatory function. Several lineages of North American watersnakes (Nerodia spp.) are known to occur in saltwater habitat and represent a useful system for studying speciation by providing an opportunity to investigate gene flow and evaluate how species boundaries are maintained or degraded. We use double digest restriction-site associated DNA sequencing to characterize the migration-selection balance and test for evidence of ecological divergence within the Nerodia fasciata-clarkii complex in Florida. We find evidence of high intraspecific gene flow with a pattern of isolation-by-distance underlying subspecific lineages. However, we identify genetic structure indicative of reduced gene flow between inland and coastal lineages suggesting divergence due to isolation-by-environment. This pattern is consistent with observed environmental differences where the amount of admixture decreases with increased salinity. Furthermore, we identify significantly enriched terms related to osmoregulatory function among a set of candidate loci, including several genes that have been previously implicated in adaptation to salinity stress. Collectively, our results demonstrate that ecological differences, likely driven by salinity, cause strong divergent selection which promotes divergence in the N. fasciata-clarkii complex despite significant gene flow

Biological and proteolytic variation in the venom of Crotalus scutulatus scutulatus from Mexico

Rattlesnake venoms may be classified according to the presence/absence and relative abundance of the neurotoxic phospholipases A2s (PLA2s), such as Mojave toxin, and snake venom metalloproteinases (SVMPs). In Mexico, studies to determine venom variation in Mojave Rattlesnakes (Crotalus scutulatus scutulatus) are limited and little is known about the biological and proteolytic activities in this species. Tissue (34) and venom (29) samples were obtained from C. s. scutulatus from different locations within their distribution in Mexico. Mojave toxin detection was carried out at the genomic (by PCR) and protein (by ELISA) levels for all tissue and venom samples. Biological activity was tested on representative venoms by measuring LD50 and hemorrhagic activity. To determine the approximate amount of SVMPs, 15 venoms were separated by RP-HPLC and variation in protein profile and proteolytic activity was evaluated by SDS-PAGE (n = 28) and Hide Powder Azure proteolytic analysis (n = 27). Three types of venom were identified in Mexico which is comparable to the intraspecific venom diversity observed in the Sonoran Desert of Arizona, USA: Venom Type A ( Type II), with Mojave toxin, highly toxic, lacking hemorrhagic activity, and with scarce proteolytic activity; Type B ( Type I), without Mojave toxin, less toxic than Type A, highly hemorrhagic and proteolytic; and Type A + B, containing Mojave toxin, as toxic as venom Type A, variable in hemorrhagic activity and with intermediate proteolytic activity. We also detected a positive correlation between SVMP abundance and hemorrhagic and proteolytic activities. Although more sampling is necessary, our results suggest that venoms containing Mojave toxin and venom lacking this toxin are distributed in the northwest and southeast portions of the distribution in Mexico, respectively, while an intergradation in the middle of both zones is presentConsejo Nacional de Ciencia y Tecnologia/[221343]/CONACYT/MéxicoUCR::Vicerrectoría de Investigación::Unidades de Investigación::Ciencias de la Salud::Instituto Clodomiro Picado (ICP)UCR::Vicerrectoría de Docencia::Salud::Facultad de Microbiologí

Ubiquitous outflows in DEEP2 spectra of star-forming galaxies at z=1.4

Galactic winds are a prime suspect for the metal enrichment of the intergalactic medium and may have a strong influence on the chemical evolution of galaxies and the nature of QSO absorption line systems. We use a sample of 1406 galaxy spectra at z~1.4 from the DEEP2 redshift survey to show that blueshifted Mg II 2796, 2803 A absorption is ubiquitous in starforming galaxies at this epoch. This is the first detection of frequent outflowing galactic winds at z~1. The presence and depth of absorption are independent of AGN spectral signatures or galaxy morphology; major mergers are not a prerequisite for driving a galactic wind from massive galaxies. Outflows are found in coadded spectra of galaxies spanning a range of 30x in stellar mass and 10x in star formation rate (SFR), calibrated from K-band and from MIPS IR fluxes. The outflows have column densities of order N_H ~ 10^20 cm^-2 and characteristic velocities of ~ 300-500 km/sec, with absorption seen out to 1000 km/sec in the most massive, highest SFR galaxies. The velocities suggest that the outflowing gas can escape into the IGM and that massive galaxies can produce cosmologically and chemically significant outflows. Both the Mg II equivalent width and the outflow velocity are larger for galaxies of higher stellar mass and SFR, with V_wind ~ SFR^0.3, similar to the scaling in low redshift IR-luminous galaxies. The high frequency of outflows in the star-forming galaxy population at z~1 indicates that galactic winds occur in the progenitors of massive spirals as well as those of ellipticals. The increase of outflow velocity with mass and SFR constrains theoretical models of galaxy evolution that include feedback from galactic winds, and may favor momentum-driven models for the wind physics.Comment: Accepted by ApJ. 25 pages, 17 figures. Revised to add discussions of intervening absorbers and AGN-driven outflows; conclusions unchange